-
The solar dynamo begins near the surface
Authors:
Geoffrey M Vasil,
Daniel Lecoanet,
Kyle Augustson,
Keaton J Burns,
Jeffrey S Oishi,
Benjamin P Brown,
Nicholas Brummell,
Keith Julien
Abstract:
The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pin…
▽ More
The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pinpoints low-latitude torsional oscillations to the Sun's outer 5-10%, the "Near-Surface Shear Layer". Within this zone, inwardly increasing differential rotation coupled with a poloidal magnetic field strongly implicates the Magneto-Rotational Instability prominent in accretion-disk theory and observed in laboratory experiments. Together, these two facts prompt the general question: Is it possible that the solar dynamo is a near-surface instability? Here, we report strong affirmative evidence in stark contrast to traditional paradigms focusing on the deeper tachocline. Simple analytic estimates show that the near-surface magneto-rotational instability better explains the spatiotemporal scales of the torsional oscillations and inferred subsurface magnetic field amplitudes. State-of-the-art numerical simulations corroborate these estimates and, strikingly, reproduce hemispherical magnetic current helicity laws. The dynamo resulting from a well-understood near-surface phenomenon improves prospects for accurate predictions of full magnetic cycles and space weather, impacting Earth's electromagnetic infrastructure.
△ Less
Submitted 13 April, 2024; v1 submitted 11 April, 2024;
originally announced April 2024.
-
Low-Frequency Internal Gravity Waves are Pseudo-incompressible
Authors:
Bradley W. Hindman,
Keith Julien
Abstract:
Starting from the fully compressible fluid equations in a plane-parallel atmosphere, we demonstrate that linear internal gravity waves are naturally pseudo-incompressible in the limit that the wave frequency $ω$ is much less than that of surface gravity waves, i.e., $ω\ll \sqrt{g k_h}$ where $g$ is the gravitational acceleration and $k_h$ is the horizontal wavenumber. We accomplish this by perform…
▽ More
Starting from the fully compressible fluid equations in a plane-parallel atmosphere, we demonstrate that linear internal gravity waves are naturally pseudo-incompressible in the limit that the wave frequency $ω$ is much less than that of surface gravity waves, i.e., $ω\ll \sqrt{g k_h}$ where $g$ is the gravitational acceleration and $k_h$ is the horizontal wavenumber. We accomplish this by performing a formal expansion of the wave functions and the local dispersion relation in terms of a dimensionless frequency $\varepsilon = ω/ \sqrt{g k_h}$. Further, we show that in this same low-frequency limit, several forms of the anelastic approximation, including the Lantz-Braginsky-Roberts (LBR) formulation, poorly reproduce the correct behavior of internal gravity waves. The pseudo-incompressible approximation is achieved by assuming that Eulerian fluctuations of the pressure are small in the continuity equation. Whereas, in the anelastic approximation Eulerian density fluctuations are ignored. In an adiabatic stratification, such as occurs in a convection zone, the two approximations become identical. But, in a stable stratification, the differences between the two approximations are stark and only the pseudo-incompressible approximation remains valid.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
The Puzzling Structure of Solar Convection: Window into the Dynamo
Authors:
Nicholas A. Featherstone,
Kyle C. Augustson,
Jonathan M. Aurnou,
Catherine Blume,
Benjamin P. Brown,
Nicholas Brummell,
Keaton J. Burns,
Michael A. Calkins,
Maria Camisassa,
Mausumi Dikpati,
Yuhong Fan,
J. R. Fuentes,
Gustavo Guerrero,
Bradley W. Hindman,
Keith Julien,
Irina N. Kitiashvili,
Lydia Korre,
Daniel Lecoanet,
Bhishek Manek,
Loren Matilsky,
Mark Miesch,
Nicholas J. Nelson,
Jeffrey S. Oishi,
Whitney T. Powers,
Matthias Rempel
, et al. (3 additional authors not shown)
Abstract:
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective…
▽ More
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements.
△ Less
Submitted 15 May, 2023;
originally announced May 2023.
-
The rotational influence on solar convection
Authors:
Geoffrey M. Vasil,
Keith Julien,
Nicholas A. Featherstone
Abstract:
This paper considers the dominant dynamical, thermal and rotational balances within the solar convection zone. The reasoning is such that: Coriolis forces balance pressure gradients. Background vortex stretching, baroclinic torques and nonlinear advection balance jointly. Turbulent fluxes convey what part of the solar luminosity that radiative diffusion cannot. These four relations determine estim…
▽ More
This paper considers the dominant dynamical, thermal and rotational balances within the solar convection zone. The reasoning is such that: Coriolis forces balance pressure gradients. Background vortex stretching, baroclinic torques and nonlinear advection balance jointly. Turbulent fluxes convey what part of the solar luminosity that radiative diffusion cannot. These four relations determine estimates for the dominant length scales and dynamical amplitudes strictly in terms of known physical quantities. We predict that the dynamical Rossby number for convection is less than unity below the near-surface shear layer, indicating strong rotational constraint. We also predict a characteristic convection length scale of roughly 30 Mm throughout much of the convection zone. These inferences help explain recent observations that reveal weak flow amplitudes at 100-200 Mm scales.
△ Less
Submitted 29 October, 2020;
originally announced October 2020.
-
Connections between Non-Rotating, Slowly Rotating, and Rapidly Rotating Turbulent Convection Transport Scalings
Authors:
Jonathan M. Aurnou,
Susanne Horn,
Keith Julien
Abstract:
Buoyancy-driven convection is likely the dominant driver of turbulent motions in the universe, and thus, is widely studied by physicists, engineers, geophysicists and astrophysicists. Maybe unsurprisingly, these different communities discuss the gross convective behaviors in different ways, often without significant cross-talk existing between them. Here, we seek to draw connections between these…
▽ More
Buoyancy-driven convection is likely the dominant driver of turbulent motions in the universe, and thus, is widely studied by physicists, engineers, geophysicists and astrophysicists. Maybe unsurprisingly, these different communities discuss the gross convective behaviors in different ways, often without significant cross-talk existing between them. Here, we seek to draw connections between these communities. We do so by carrying out a set of basic scale estimations for how heat and fluid momentum transport should behave in non-rotating, slowly rotating and rapidly rotating buoyancy-driven convective environments. We find that slowly and rapidly rotating scalings can be inter-related via one parameter, the so-called convective Rossby number $\RoC$, a dissipation-free parameter measuring the importance of buoyancy driving relative to rotation. Further, we map between non-flux-based and the flux-based, buoyancy-driven scalings used by different groups. In doing so, these scalings show that there are clean connections between the different communities' approaches and that a number of the seemingly different scalings are actually synonymous with one another.
△ Less
Submitted 7 September, 2020;
originally announced September 2020.
-
Inertia-less convectively-driven dynamo models in the limit of low Rossby number and large Prandtl number
Authors:
Michael A. Calkins,
Keith Julien,
Steven M. Tobias
Abstract:
Compositional convection is thought to be an important energy source for magnetic field generation within planetary interiors. The Prandtl number, $Pr$, characterizing compositional convection is significantly larger than unity, suggesting that the inertial force may not be important on the small scales of convection as long as the buoyancy force is not too strong. We develop asymptotic dynamo mod…
▽ More
Compositional convection is thought to be an important energy source for magnetic field generation within planetary interiors. The Prandtl number, $Pr$, characterizing compositional convection is significantly larger than unity, suggesting that the inertial force may not be important on the small scales of convection as long as the buoyancy force is not too strong. We develop asymptotic dynamo models for the case of small Rossby number and large Prandtl number in which inertia is absent on the convective scale. The relevant diffusivity parameter for this limit is the compositional Roberts number, $q = D/η$, which is the ratio of compositional and magnetic diffusivities. Dynamo models are developed for both order one $q$ and the more geophysically relevant low $q$ limit. For both cases the ratio of magnetic to kinetic energy densities, $M$, is asymptotically large and reflects the fact that Alfvén waves have been filtered from the dynamics. Along with previous investigations of asymptotic dynamo models for $Pr=O(1)$, our results show that the ratio $M$ is not a useful indicator of dominant force balances in the momentum equation since many different asymptotic limits of $M$ can be obtained without changing the leading order geostrophic balance. Furthermore, the present models show that inertia is not a requirement for driving low $q$, large-scale dynamos.
△ Less
Submitted 23 February, 2017; v1 submitted 13 September, 2016;
originally announced September 2016.
-
Convection-driven kinematic dynamos at low Rossby and magnetic Prandtl numbers: single mode solutions
Authors:
Michael A. Calkins,
Keith Julien,
Steven M. Tobias,
Jonathan M. Aurnou,
Philippe Marti
Abstract:
The onset of dynamo action is investigated within the context of a newly developed low Rossby, low magnetic Prandtl number, convection-driven dynamo model. This multiscale model represents an asymptotically exact form of an $α^2$ mean field dynamo model in which the small-scale convection is represented explicitly by finite amplitude, single mode solutions. Both steady and oscillatory convection a…
▽ More
The onset of dynamo action is investigated within the context of a newly developed low Rossby, low magnetic Prandtl number, convection-driven dynamo model. This multiscale model represents an asymptotically exact form of an $α^2$ mean field dynamo model in which the small-scale convection is represented explicitly by finite amplitude, single mode solutions. Both steady and oscillatory convection are considered for a variety of horizontal planforms. The kinetic helicity is observed to be a monotonically increasing function of the Rayleigh number. As a result, very small magnetic Prandtl number dynamos can be found for sufficiently large Rayleigh numbers. All dynamos are found to be oscillatory with an oscillation frequency that increases as the strength of the convection is increased and the magnetic Prandtl number is reduced. Kinematic dynamo action is strongly controlled by the profile of the helicity; single mode solutions which exhibit boundary layer behavior in the helicity show a decrease in the efficiency of dynamo action due to the enhancement of magnetic diffusion in the boundary layer regions. For a given value of the Rayleigh number, lower magnetic Prandtl number dynamos are excited for the case of oscillatory convection in comparison to steady convection. With regards to planetary dynamos, these results suggest that the low magnetic Prandtl number dynamos typical of liquid metals are more easily driven by thermal convection than by compositional convection.
△ Less
Submitted 31 January, 2016; v1 submitted 19 November, 2015;
originally announced November 2015.
-
Baroclinic Vorticity Production in Protoplanetary Disks; Part I: Vortex Formation
Authors:
Mark R. Petersen,
Keith Julien,
Glen R. Stewart
Abstract:
The formation of vortices in protoplanetary disks is explored via pseudo-spectral numerical simulations of an anelastic-gas model. This model is a coupled set of equations for vorticity and temperature in two dimensions which includes baroclinic vorticity production and radiative cooling. Vortex formation is unambiguously shown to be caused by baroclinicity because (1) these simulations have zer…
▽ More
The formation of vortices in protoplanetary disks is explored via pseudo-spectral numerical simulations of an anelastic-gas model. This model is a coupled set of equations for vorticity and temperature in two dimensions which includes baroclinic vorticity production and radiative cooling. Vortex formation is unambiguously shown to be caused by baroclinicity because (1) these simulations have zero initial perturbation vorticity and a nonzero initial temperature distribution; and (2) turning off the baroclinic term halts vortex formation, as shown by an immediate drop in kinetic energy and vorticity. Vortex strength increases with: larger background temperature gradients; warmer background temperatures; larger initial temperature perturbations; higher Reynolds number; and higher resolution. In the simulations presented here vortices form when the background temperatures are $\sim 200K$ and vary radially as $r^{-0.25}$, the initial vorticity perturbations are zero, the initial temperature perturbations are 5% of the background, and the Reynolds number is $10^9$. A sensitivity study consisting of 74 simulations showed that as resolution and Reynolds number increase, vortices can form with smaller initial temperature perturbations, lower background temperatures, and smaller background temperature gradients. For the parameter ranges of these simulations, the disk is shown to be convectively stable by the Solberg-Høiland criteria.
△ Less
Submitted 16 November, 2006;
originally announced November 2006.
-
Baroclinic Vorticity Production in Protoplanetary Disks; Part II: Vortex Growth and Longevity
Authors:
Mark R. Petersen,
Glen R. Stewart,
Keith Julien
Abstract:
The factors affecting vortex growth in convectively stable protoplanetary disks are explored using numerical simulations of a two-dimensional anelastic-gas model which includes baroclinic vorticity production and radiative cooling. The baroclinic feedback, where anomalous temperature gradients produce vorticity through the baroclinic term and vortices then reinforce these temperature gradients,…
▽ More
The factors affecting vortex growth in convectively stable protoplanetary disks are explored using numerical simulations of a two-dimensional anelastic-gas model which includes baroclinic vorticity production and radiative cooling. The baroclinic feedback, where anomalous temperature gradients produce vorticity through the baroclinic term and vortices then reinforce these temperature gradients, is found to be an important process in the rate of growth of vortices in the disk. Factors which strengthen the baroclinic feedback include fast radiative cooling, high thermal diffusion, and large radial temperature gradients in the background temperature. When the baroclinic feedback is sufficiently strong, anticyclonic vortices form from initial random perturbations and maintain their strength for the duration of the simulation, for over 600 orbital periods.
Based on both simulations and a simple vortex model, we find that the local angular momentum transport due to a single vortex may be inward or outward, depending its orientation. The global angular momentum transport is highly variable in time, and is sometimes negative and sometimes positive. This result is for an anelastic gas model, and does not include shocks that could affect angular momentum transport in a compressible-gas disk.
△ Less
Submitted 5 December, 2006; v1 submitted 16 November, 2006;
originally announced November 2006.
