The origin of grand minima in the sunspot cycle

One of the most striking aspects of the 11-year sunspot cycle is that there have been times in the past when some cycles went missing, a most well-known example of this being the Maunder minimum during 1645-1715. Analyses of cosmogenic isotopes (C14 and Be10) indicated that there were about 27 grand minima in the last 11,000 yr, implying that about 2.7% of the solar cycles had conditions appropriate for forcing the Sun into grand minima. We address the question how grand minima are produced and specifically calculate the frequency of occurrence of grand minima from a theoretical dynamo model. We assume that fluctuations in the poloidal field generation mechanism and the meridional circulation produce irregularities of sunspot cycles. Taking these fluctuations to be Gaussian and estimating the values of important parameters from the data of last 28 solar cycles, we show from our flux transport dynamo model that about 1-4% of the sunspot cycles may have conditions suitable for inducing grand minima.

💡 Research Summary

The paper tackles one of the most intriguing puzzles in solar physics: why the regular 11‑year sunspot cycle sometimes disappears for several decades, producing so‑called grand minima such as the Maunder Minimum (1645‑1715). Cosmogenic isotope records (¹⁴C and ¹⁰Be) indicate that roughly 27 grand minima have occurred over the past 11 000 years, corresponding to about 2.7 % of all solar cycles. The authors ask whether a standard flux‑transport dynamo, perturbed by realistic stochastic fluctuations, can reproduce this frequency without invoking exotic external triggers.

To answer this, they adopt a kinematic flux‑transport dynamo model in which the toroidal field (the source of sunspots) is generated by differential rotation, while the poloidal field is regenerated through the Babcock‑Leighton (BL) mechanism. Two ingredients are allowed to fluctuate: (i) the efficiency of the BL poloidal source, denoted α, and (ii) the speed of the meridional circulation, v. Both are assumed to follow independent Gaussian distributions. The mean values and standard deviations of α and v are estimated from the most recent 28 solar cycles (approximately 300 years) using observed sunspot numbers and helioseismic flow measurements.

The governing equations are the standard axisymmetric mean‑field dynamo equations for the toroidal field Bφ(r,θ,t) and the vector potential A(r,θ,t) of the poloidal field. The α‑effect enters as a source term in the A‑equation, while v appears in the advection terms that transport magnetic flux from the surface to the tachocline. By randomly drawing α and v for each synthetic cycle, the authors generate thousands of realizations of the solar cycle using a Monte‑Carlo approach. A grand minimum is defined operationally as a period of three or more consecutive cycles in which the sunspot number falls below 5 % of the long‑term mean.

The statistical outcome is strikingly simple: when both α and v dip below their respective means simultaneously, the dynamo’s feedback loop weakens enough to drive the toroidal field to near‑zero values, producing a grand‑minimum episode. In the ensemble of simulations, such simultaneous low‑α/low‑v events occur in roughly 1 %–4 % of all cycles, depending on the exact choice of standard deviations. This range comfortably brackets the 2.7 % figure derived from isotope data, suggesting that the observed grand‑minimum frequency can be explained purely by internal stochastic variability. Moreover, the model reproduces the observed spread in grand‑minimum durations: when the low‑α/low‑v condition persists, the minimum can last for several decades up to a few centuries, consistent with the 70‑year duration of the Maunder Minimum.

The authors discuss several implications. First, the result supports the view that grand minima are not the consequence of rare external events (e.g., planetary alignments) but are an intrinsic property of a nonlinear dynamo subject to ordinary fluctuations. Second, the agreement between model and data validates the flux‑transport paradigm as a realistic description of the solar magnetic cycle. Third, the study highlights the importance of the meridional flow speed as a regulator of cycle amplitude; a slower flow lengthens the residence time of magnetic flux in the tachocline, allowing diffusion to erode the field more effectively.

Limitations are acknowledged. The Gaussian assumption for α and v may underestimate the probability of extreme excursions, as turbulent convection can produce heavy‑tailed statistics. The BL mechanism itself may be modulated by other processes (e.g., turbulent pumping, magnetic helicity constraints) that are not included. The model is axisymmetric and kinematic, ignoring back‑reaction of the magnetic field on flows, which could introduce additional nonlinearities. Future work is suggested to explore non‑Gaussian noise, fully dynamic MHD simulations, and longer observational constraints (e.g., extending the isotope record with higher temporal resolution).

In conclusion, by coupling realistic stochastic fluctuations to a well‑established flux‑transport dynamo, the authors demonstrate that the Sun naturally spends about 1 %–4 % of its cycles in a grand‑minimum state. This quantitative match with the paleoclimate record provides a compelling, self‑consistent explanation for the intermittent disappearance of sunspots and underscores the role of internal dynamo variability in shaping long‑term solar activity.