[Paper Review] Cycling Dynamo in a Young Sun: Grand Minima and Equatorward Propagation

Kyle Augustson, S. Brun+2 more|arXiv (Cornell University)|Oct 31, 2013

Advanced Thermodynamics and Statistical Mechanics1 references5 citations

TL;DR

This study presents a 3D anelastic MHD simulation of a young Sun rotating three times faster than the present-day Sun, using the ASH code to model convective dynamos. It reproduces key solar magnetic features: 6.2-year magnetic polarity cycles with equatorward migration of toroidal fields, Lorentz-force-driven magnetic energy cycles, and a 20-year grand minimum, demonstrating self-consistent dynamo action with realistic solar-like behavior despite rapid rotation and scale parameterization.

ABSTRACT

We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles.

Motivation & Objective

• To investigate how rapid rotation affects large-scale dynamo action in a young Sun-like star.
• To determine whether global-scale MHD simulations can reproduce key solar magnetic phenomena such as polarity cycles, equatorward migration, and grand minima.
• To examine the role of Lorentz-force feedback on differential rotation and its impact on magnetic energy and polarity cycles.
• To explore the influence of near-surface shear layers on dynamo wave propagation and field organization.
• To assess whether parameterized convection and rotation can still yield realistic solar-like magnetic behavior in global-scale simulations.

Proposed method

• The simulation uses the ASH (Anelastic Spherical Harmonic) code to solve the anelastic MHD equations in a spherical shell spanning 3.8 density scale heights of the solar convection zone.
• Spherical harmonic decomposition is used for horizontal derivatives, while fourth-order non-uniform finite differences resolve radial derivatives.
• A stream function formalism maintains solenoidality of mass flux and magnetic fields, and slope-limited diffusion is applied to stabilize the system.
• The model includes a weak negative radial gradient in angular velocity in the upper 10% of the domain, simulating a near-surface shear layer.
• The simulation runs for 100 years, with time-averaged and time-resolved analysis of magnetic energy, toroidal field structure, and polarity evolution.
• Key diagnostics include time-latitude diagrams of azimuthally-averaged toroidal field, energy cycle analysis, and meridional field structure visualization.

Experimental results

Research questions

• RQ1 Can a rapidly rotating, young Sun-like star sustain a self-consistent dynamo that reproduces the 11-year sunspot cycle and associated magnetic phenomena?
• RQ2 What mechanisms drive the equatorward migration of toroidal magnetic field structures in a global MHD simulation?
• RQ3 How does Lorentz-force feedback on differential rotation influence the periodicity and amplitude of magnetic energy cycles?
• RQ4 Can the simulation reproduce a grand minimum with reduced magnetic activity and polarity cycle disruption, similar to the Maunder Minimum?
• RQ5 To what extent do near-surface shear layers and spatial separation of equatorial and polar dynamos contribute to polarity reversals and cycle organization?

Key findings

• The simulation exhibits a 6.2-year magnetic polarity cycle, with the dipole moment returning to its initial orientation after each cycle, indicating stable dynamo action.
• Magnetic energy cycles with a period of 12.4 years (two cycles per polarity reversal) arise from nonlinear feedback between differential rotation and Lorentz forces.
• Equatorward migration of toroidal field structures is driven by the location of maximum latitudinal shear and a weak negative radial shear in the upper convection zone.
• A 20-year grand minimum occurs, during which polarity reversals cease and magnetic energy drops by 25% globally (60% in lower latitudes), followed by recovery of regular cycles.
• The Lorentz force significantly quenches differential rotation, reducing its role in toroidal field generation during polarity reversals.
• The simulation captures the spatial and temporal coherence of solar-like magnetic cycles, including poleward advection of reversed field after equatorial cancellation, consistent with observed dynamo wave behavior.