Christopher A. Jones
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Christopher A. Jones has a total of 162 co-authors. Here are the top co-authors by number of shared publications:
K. Zhang
Department of Mathematics, University of Exeter, Exeter EX4 4QE, UK
13 shared publications- Steven M. Tobias
School of Physics and Astronomy University of Edinburgh Edinburgh UK
8 shared publications - Rainer Hollerbach
School of Mathematics, University of Leeds , Leeds LS2 9JT
7 shared publications C. F. Barenghi
Department of Mathematics , University of Newcastle Upon Tyne , NE1 7RU, UK
4 shared publications- Christopher J. Davies
University of Leeds
3 shared publications - Céline Guervilly
School of Mathematics Statistics and Physics Newcastle University Newcastle upon Tyne NE17RU UK
3 shared publications
Latest Publications
- Open Access
- 7 October 2025
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 10 (10)
Model of deep zonal flows in giant planets
A mechanism by which the surface zonal flows of giant planets can be gradually attenuated with depth is explored. The zonal flow is driven by an imposed forcing in a thin layer near the surface. A meridional circulation is set up, analogous to the Ferrel-like cells observed in Jupiter's atmosphere. Acting on a stably stratified thin surface layer, the meridional flow induces a horizontal temperature anomaly which leads to a gradual reduction of the zonal winds with depth, governed by the thermal wind equation. Our model is a Boussinesq plane layer, with gravity acting parallel to the rotation axis. A suite of fully three-dimensional time-dependent numerical simulations has been performed to investigate the model behavior. Below the forced stable layer, convection is occurring, typically in the form of tall thin cells. The fluctuating components of the three-dimensional flow can be driven by either the convection or the Reynolds stresses associated with the jet shear flow. When fluctuations are mainly driven by convection in the form of tall thin columns and the forcing is relatively weak, the horizontal temperature anomaly persists much deeper into the interior than when it is driven by shear flow. The zonal jets can therefore extend deep into the interior, consistent with the Juno gravity data.
- 7 October 2025
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 10 (10)
Open AccessModel of deep zonal flows in giant planets
A mechanism by which the surface zonal flows of giant planets can be gradually attenuated with depth is explored. The zonal flow is driven by an imposed forcing in a thin layer near the surface. A meridional circulation is set up, analogous to the Ferrel-like cells observed in Jupiter's atmosphere. Acting on a stably stratified thin surface layer, the meridional flow induces a horizontal temperature anomaly which leads to a gradual reduction of the zonal winds with depth, governed by the thermal wind equation. Our model is a Boussinesq plane layer, with gravity acting parallel to the rotation axis. A suite of fully three-dimensional time-dependent numerical simulations has been performed to investigate the model behavior. Below the forced stable layer, convection is occurring, typically in the form of tall thin cells. The fluctuating components of the three-dimensional flow can be driven by either the convection or the Reynolds stresses associated with the jet shear flow. When fluctuations are mainly driven by convection in the form of tall thin columns and the forcing is relatively weak, the horizontal temperature anomaly persists much deeper into the interior than when it is driven by shear flow. The zonal jets can therefore extend deep into the interior, consistent with the Juno gravity data.
- 4 July 2025
- Published by Taylor & Francis in Geophysical & Astrophysical Fluid Dynamics
- Vol. 119 (3-4) , 247-252
Introduction
- 4 July 2025
- Published by Taylor & Francis in Geophysical & Astrophysical Fluid Dynamics
- Vol. 119 (3-4) , 247-252
Introduction
- Open Access
Low inertia reversing geodynamos
Convection driven geodynamo models in rotating spherical geometry have regimes in which reversals occur. However, reversing dynamo models are usually found in regimes where the kinetic and magnetic energy is comparable, so that inertia is playing a significant role in the dynamics. In the Earth's core, the Rossby number is very small, and the magnetic energy is much larger than the kinetic energy. Here we investigate dynamo models in the strong field regime, where magnetic forces have a significant effect on convection. In the core, the strong field is achieved by having the magnetic Prandtl number Pm small, but the Ekman number E extremely small. In simulations, very small E is not possible, but the strong field regime can be reached by increasing Pm. However, if Pm is raised while the fluid Prandtl number is fixed at unity, the most common choice, the Peclet number number becomes small, so that the linear terms in the heat (or composition) equation dominate, which is also far from Earth-like behaviour. Here we increase Pr and Pm together, so that nonlinearity is important in the heat equation and the dynamo is strong field. We find that Earth-like reversals are possible at numerically achievable parameter values, and the simulations have Earth-like magnetic fields away from the times at which it reverses. The magnetic energy is much greater than the kinetic energy except close to reversal times.
Open AccessLow inertia reversing geodynamos
Convection driven geodynamo models in rotating spherical geometry have regimes in which reversals occur. However, reversing dynamo models are usually found in regimes where the kinetic and magnetic energy is comparable, so that inertia is playing a significant role in the dynamics. In the Earth's core, the Rossby number is very small, and the magnetic energy is much larger than the kinetic energy. Here we investigate dynamo models in the strong field regime, where magnetic forces have a significant effect on convection. In the core, the strong field is achieved by having the magnetic Prandtl number Pm small, but the Ekman number E extremely small. In simulations, very small E is not possible, but the strong field regime can be reached by increasing Pm. However, if Pm is raised while the fluid Prandtl number is fixed at unity, the most common choice, the Peclet number number becomes small, so that the linear terms in the heat (or composition) equation dominate, which is also far from Earth-like behaviour. Here we increase Pr and Pm together, so that nonlinearity is important in the heat equation and the dynamo is strong field. We find that Earth-like reversals are possible at numerically achievable parameter values, and the simulations have Earth-like magnetic fields away from the times at which it reverses. The magnetic energy is much greater than the kinetic energy except close to reversal times.
- 14 December 2023
- Published by Springer Nature in Nature Astronomy
- Vol. 7 (12) , 1417-1418
The deep winds of Jupiter
- 14 December 2023
- Published by Springer Nature in Nature Astronomy
- Vol. 7 (12) , 1417-1418
The deep winds of Jupiter
- Open Access
- 24 November 2023
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 976
A two-phase pure slurry model for planetary cores: one-dimensional solutions and implications for Earth's F-layer
We develop and analyse a continuum model of two-phase slurry dynamics for planetary cores. Mixed solid–liquid slurry regions may be ubiquitous in the upper cores of small terrestrial bodies and have also been invoked to explain anomalous seismic structure in the F-layer at the base of Earth's liquid iron core. These layers are expected to influence the dynamics and evolution of planetary cores, including their capacity to generate global magnetic fields; however, to date, models of two-phase regions in planetary cores have largely ignored the complex fluid dynamics that arises from interactions between phases. As an initial application of our model, and to focus on fundamental fluid dynamical processes, we consider a non-rotating and non-magnetic slurry comprised of a single chemical component with a temperature that is tied to the liquidus. We study one-dimensional solutions in a configuration set up to mimic Earth's F-layer, varying gravitational strength $R$ , the solid/liquid viscosity ratio $\lambda _{\mu }$ and the interaction parameter $K$ , which measures friction between the phases. We develop scalings describing behaviour in the limit $R \gg 1$ and $\lambda _{\mu } \gg 1$ , which are in excellent agreement with our numerical results. Application to Earth's core, where $R \sim 10^{28}$ and $\lambda _{\mu } \sim 10^{22}$ , suggests that a pure iron slurry F-layer would contain a mean solid fraction of at most 5 %.- 24 November 2023
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 976
Open AccessA two-phase pure slurry model for planetary cores: one-dimensional solutions and implications for Earth's F-layer
We develop and analyse a continuum model of two-phase slurry dynamics for planetary cores. Mixed solid–liquid slurry regions may be ubiquitous in the upper cores of small terrestrial bodies and have also been invoked to explain anomalous seismic structure in the F-layer at the base of Earth's liquid iron core. These layers are expected to influence the dynamics and evolution of planetary cores, including their capacity to generate global magnetic fields; however, to date, models of two-phase regions in planetary cores have largely ignored the complex fluid dynamics that arises from interactions between phases. As an initial application of our model, and to focus on fundamental fluid dynamical processes, we consider a non-rotating and non-magnetic slurry comprised of a single chemical component with a temperature that is tied to the liquidus. We study one-dimensional solutions in a configuration set up to mimic Earth's F-layer, varying gravitational strength $R$ , the solid/liquid viscosity ratio $\lambda _{\mu }$ and the interaction parameter $K$ , which measures friction between the phases. We develop scalings describing behaviour in the limit $R \gg 1$ and $\lambda _{\mu } \gg 1$ , which are in excellent agreement with our numerical results. Application to Earth's core, where $R \sim 10^{28}$ and $\lambda _{\mu } \sim 10^{22}$ , suggests that a pure iron slurry F-layer would contain a mean solid fraction of at most 5 %. - 18 May 2023
- Published by Springer Nature in Nature Astronomy
- Vol. 7 (7) , 825-835
Jupiter’s cloud-level variability triggered by torsional oscillations in the interior
- 18 May 2023
- Published by Springer Nature in Nature Astronomy
- Vol. 7 (7) , 825-835
Jupiter’s cloud-level variability triggered by torsional oscillations in the interior
- Open Access
- 8 November 2021
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 930, A13
Fully developed anelastic convection with no-slip boundaries
Anelastic convection at high Rayleigh number in a plane parallel layer with no slip boundaries is considered. Energy and entropy balance equations are derived, and they are used to develop scaling laws for the heat transport and the Reynolds number. The appearance of an entropy structure consisting of a well-mixed uniform interior, bounded by thin layers with entropy jumps across them, makes it possible to derive explicit forms for these scaling laws. These are given in terms of the Rayleigh number, the Prandtl number and the bottom to top temperature ratio, which also measures how much the density varies across the layer. The top and bottom boundary layers are examined and they are found to be very different, unlike in the Boussinesq case. Elucidating the structure of these boundary layers plays a crucial part in determining the scaling laws. Physical arguments governing these boundary layers are presented, concentrating on the case in which the boundary layers are so thin that temperature and density vary little across them, even though there may be substantial temperature and density variations across the whole layer. Different scaling laws are found, depending on whether the viscous dissipation is primarily in the boundary layers or in the bulk. The cases of both high and low Prandtl number are considered. Numerical simulations of no-slip anelastic convection up to a Rayleigh number of $10^7$ have been performed and our theoretical predictions are compared with the numerical results.- 8 November 2021
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 930, A13
Open AccessFully developed anelastic convection with no-slip boundaries
Anelastic convection at high Rayleigh number in a plane parallel layer with no slip boundaries is considered. Energy and entropy balance equations are derived, and they are used to develop scaling laws for the heat transport and the Reynolds number. The appearance of an entropy structure consisting of a well-mixed uniform interior, bounded by thin layers with entropy jumps across them, makes it possible to derive explicit forms for these scaling laws. These are given in terms of the Rayleigh number, the Prandtl number and the bottom to top temperature ratio, which also measures how much the density varies across the layer. The top and bottom boundary layers are examined and they are found to be very different, unlike in the Boussinesq case. Elucidating the structure of these boundary layers plays a crucial part in determining the scaling laws. Physical arguments governing these boundary layers are presented, concentrating on the case in which the boundary layers are so thin that temperature and density vary little across them, even though there may be substantial temperature and density variations across the whole layer. Different scaling laws are found, depending on whether the viscous dissipation is primarily in the boundary layers or in the bulk. The cases of both high and low Prandtl number are considered. Numerical simulations of no-slip anelastic convection up to a Rayleigh number of $10^7$ have been performed and our theoretical predictions are compared with the numerical results. - Open Access
- 12 October 2020
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 904
Solitary magnetostrophic Rossby waves in spherical shells
Finite-amplitude hydromagnetic Rossby waves in the magnetostrophic regime are studied. We consider the slow mode, which travels in the opposite direction to the hydrodynamic or fast mode, in the presence of a toroidal magnetic field and zonal flow by means of quasi-geostrophic models for thick spherical shells. The weakly nonlinear long waves are derived asymptotically using a reductive perturbation method. The problem at the first order is found to obey a second-order ordinary differential equation, leading to a hypergeometric equation for a Malkus field and a confluent Heun equation for an electrical wire field, and is non-singular when the wave speed approaches the mean flow. Investigating its neutral non-singular eigensolutions for different basic states, we find the evolution is described by the Korteweg–de Vries equation. This implies that the nonlinear slow wave forms solitons and solitary waves. These may take the form of a coherent eddy, such as a single anticyclone. We speculate on the relation of the anticyclone to the asymmetric gyre seen in the Earth's fluid core, and in state-of-the-art dynamo direct numerical simulations.
- 12 October 2020
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 904
Open AccessSolitary magnetostrophic Rossby waves in spherical shells
Finite-amplitude hydromagnetic Rossby waves in the magnetostrophic regime are studied. We consider the slow mode, which travels in the opposite direction to the hydrodynamic or fast mode, in the presence of a toroidal magnetic field and zonal flow by means of quasi-geostrophic models for thick spherical shells. The weakly nonlinear long waves are derived asymptotically using a reductive perturbation method. The problem at the first order is found to obey a second-order ordinary differential equation, leading to a hypergeometric equation for a Malkus field and a confluent Heun equation for an electrical wire field, and is non-singular when the wave speed approaches the mean flow. Investigating its neutral non-singular eigensolutions for different basic states, we find the evolution is described by the Korteweg–de Vries equation. This implies that the nonlinear slow wave forms solitons and solitary waves. These may take the form of a coherent eddy, such as a single anticyclone. We speculate on the relation of the anticyclone to the asymmetric gyre seen in the Earth's fluid core, and in state-of-the-art dynamo direct numerical simulations.
- Open Access
- 4 May 2020
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 894
Viscous and inviscid strato-rotational instability
We examine the critical viscous mode of the Taylor–Couette strato-rotational instability, concentrating on cases where the buoyancy frequency $N$ and the inner cylinder rotation rate $\unicode[STIX]{x1D6FA}_{in}$ are comparable, giving a detailed account for $N=\unicode[STIX]{x1D6FA}_{in}$ . The ratio of the outer to the inner cylinder rotation rates $\unicode[STIX]{x1D707}=\unicode[STIX]{x1D6FA}_{out}/\unicode[STIX]{x1D6FA}_{in}$ and the ratio of the inner to the outer cylinder radius $\unicode[STIX]{x1D702}=r_{in}/r_{out}$ satisfy $0<\unicode[STIX]{x1D707}<1$ and $0<\unicode[STIX]{x1D702}<1$ . We find considerable variation in the structure of the mode, and the critical Reynolds number $Re_{c}$ at which the flow becomes unstable. For $N=\unicode[STIX]{x1D6FA}_{in}$ , we classify different regions of the $\unicode[STIX]{x1D702}\unicode[STIX]{x1D707}$ -plane by the critical viscous mode of each region. We find that there is a triple point in the $\unicode[STIX]{x1D702}\unicode[STIX]{x1D707}$ -plane where three different viscous modes all onset at the same Reynolds number. We also find a discontinuous change in $Re_{c}$ along a curve in the $\unicode[STIX]{x1D702}\unicode[STIX]{x1D707}$ -plane, on one side of which exist closed unstable domains where the flow can restabilise when the Reynolds number is increased. A new form of viscous instability occurring for wide gaps has been detected. We show for the first time that there is a region of the parameter space for which the critical viscous mode at the onset of instability corresponds to the inviscid radiative instability of Le Dizès & Riedinger (J. Fluid Mech., vol. 660, 2010, pp. 147–161). Focusing on small-to-moderate wavenumbers, we demonstrate that the viscous and inviscid systems are not always correlated. We explore which viscous modes relate to inviscid modes and which do not. For asymptotically large vertical wavenumbers, we have extended the inviscid analysis of Park & Billant (J. Fluid Mech., vol. 725, 2013, pp. 262–280) to cover the cases where $N$ and $\unicode[STIX]{x1D6FA}_{in}$ are comparable.- 4 May 2020
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 894
Open AccessViscous and inviscid strato-rotational instability
We examine the critical viscous mode of the Taylor–Couette strato-rotational instability, concentrating on cases where the buoyancy frequency $N$ and the inner cylinder rotation rate $\unicode[STIX]{x1D6FA}_{in}$ are comparable, giving a detailed account for $N=\unicode[STIX]{x1D6FA}_{in}$ . The ratio of the outer to the inner cylinder rotation rates $\unicode[STIX]{x1D707}=\unicode[STIX]{x1D6FA}_{out}/\unicode[STIX]{x1D6FA}_{in}$ and the ratio of the inner to the outer cylinder radius $\unicode[STIX]{x1D702}=r_{in}/r_{out}$ satisfy $0<\unicode[STIX]{x1D707}<1$ and $0<\unicode[STIX]{x1D702}<1$ . We find considerable variation in the structure of the mode, and the critical Reynolds number $Re_{c}$ at which the flow becomes unstable. For $N=\unicode[STIX]{x1D6FA}_{in}$ , we classify different regions of the $\unicode[STIX]{x1D702}\unicode[STIX]{x1D707}$ -plane by the critical viscous mode of each region. We find that there is a triple point in the $\unicode[STIX]{x1D702}\unicode[STIX]{x1D707}$ -plane where three different viscous modes all onset at the same Reynolds number. We also find a discontinuous change in $Re_{c}$ along a curve in the $\unicode[STIX]{x1D702}\unicode[STIX]{x1D707}$ -plane, on one side of which exist closed unstable domains where the flow can restabilise when the Reynolds number is increased. A new form of viscous instability occurring for wide gaps has been detected. We show for the first time that there is a region of the parameter space for which the critical viscous mode at the onset of instability corresponds to the inviscid radiative instability of Le Dizès & Riedinger (J. Fluid Mech., vol. 660, 2010, pp. 147–161). Focusing on small-to-moderate wavenumbers, we demonstrate that the viscous and inviscid systems are not always correlated. We explore which viscous modes relate to inviscid modes and which do not. For asymptotically large vertical wavenumbers, we have extended the inviscid analysis of Park & Billant (J. Fluid Mech., vol. 725, 2013, pp. 262–280) to cover the cases where $N$ and $\unicode[STIX]{x1D6FA}_{in}$ are comparable. - Open Access
- 1 January 2020
- Published by Elsevier in Earth and Planetary Science Letters
- Vol. 530
Characterising Jupiter's dynamo radius using its magnetic energy spectrum
- 1 January 2020
- Published by Elsevier in Earth and Planetary Science Letters
- Vol. 530
Open AccessCharacterising Jupiter's dynamo radius using its magnetic energy spectrum