The chain of causality that makes such reconstructions possible begins far beyond Earth. The Sun’s magnetic field, dragged outward by the solar wind, fills the heliosphere, the vast protective bubble surrounding the solar system. Galactic cosmic rays, high-energy particles arriving from outside, must fight their way through this magnetic shielding to reach Earth. When the Sun is magnetically active, the shielding strengthens and fewer cosmic rays penetrate; when the Sun quiets, the floodgates loosen. Upon entering the atmosphere, cosmic rays collide with nitrogen, oxygen and argon, spawning showers of secondary particles that forge rare radionuclides such as carbon-14 and beryllium-10. Carbon-14, once formed, is incorporated into carbon dioxide, absorbed by trees during photosynthesis and permanently archived in annual growth rings. Because a magnetically active Sun suppresses cosmic ray influx, the abundance of these isotopes in wood and ice is inversely correlated with solar activity, giving scientists a decipherable record stretching back thousands of years.

Deciphering it, however, is harder than it sounds. Previous reconstruction efforts relied on statistical regressions that mapped isotope-derived quantities onto sunspot numbers, and those regressions carried hidden dangers. As the new study’s authors, led by Chitradeep Saha of the University of Reading, point out, even a regression with an impressively high correlation coefficient can misfire if the data contain non-linearities, zero-level offsets, or uneven variance across amplitudes. The most notorious failure mode appears during grand minima, the extended intervals when solar activity collapses, such as the Maunder minimum of the seventeenth century. When regression equations calibrated on modern data are extrapolated to these unusually quiet conditions, they can yield sunspot numbers below zero. Values are often simply clipped to zero, but that crude fix distorts cycle averages and inflates estimates of the Sun’s total irradiance, which in turn muddies attempts to understand the Sun’s influence on past climate.

The Reading-led team, which also includes Mathew Owens, Mike Lockwood and Luke Barnard, together with colleagues at ETH Zurich, the University of Lancashire, the University of Oslo, Lund University and the University of Groningen, took a different path. Rather than inverting the physics with a regression, they ran it forward, over and over, in a Monte Carlo framework. The method, implemented in publicly released code under the name PRISM, begins by generating an ensemble of thousands of hypothetical sunspot cycles drawn from statistical priors: cycle amplitudes sampled from a log-normal distribution, cycle lengths from a Gaussian distribution centred on 10.5 years, and a random start offset for each window of time. Each trial cycle is then passed through a sequence of two semi-empirical forward models that translate sunspot number into open solar flux, the total magnetic flux threading the outer boundary of the corona, and then into the heliospheric modulation potential, a quantity describing how much energy cosmic rays lose as they traverse the heliosphere.

The forward models rest on decades of established solar physics. The first couples sunspot number to the emergence of new magnetic flux through a empirically optimised source function, balanced against a phase-dependent loss rate derived from solar cycles 13 through 24. The second model computes the modulation potential from the open flux together with the tilt and polarity of the heliospheric current sheet, following a formulation calibrated by Owens and colleagues in 2024. Crucially, the open solar flux evolves with memory: it accumulates from past sunspot activity and decays through magnetic reconnection, introducing a hysteresis that makes the inverse problem fundamentally non-unique. Multiple distinct sunspot histories can produce statistically indistinguishable modulation records, which is precisely why simple deterministic inversions break down.

To handle that non-uniqueness, the team employed Approximate Bayesian Computation, a statistical technique that sidesteps the need for an explicit likelihood function. In each sliding ten-to-fifteen-year window, ten thousand Monte Carlo realisations of sunspot cycles are propagated through the forward models and compared directly against the observed modulation potential using a weighted Euclidean distance. The best two percent of candidates, some two hundred realisations, are retained as samples from the approximate posterior distribution, and their spread provides rigorous, quantified uncertainty bounds reported as 68 percent highest-density intervals. Because the sunspot cycle amplitudes are constrained to be non-negative by construction, the resulting reconstruction can never produce the negative values that haunted earlier regression-based approaches, and it requires no post-hoc correction.

To test the method, the researchers applied it to two annual-resolution records of the modulation potential. The first, spanning 1845 to 2020, was derived from geomagnetic observations of open solar flux by Owens and colleagues. When used as the inversion target, the method recovered sunspot numbers and open flux in close agreement with the direct instrumental record maintained by SILSO, with a mean absolute error of just 19.25 megavolts, about three percent of the mean modulation potential. That success validated the technique and justified applying it to the second, far longer dataset: a radiocarbon-based modulation potential record covering 971 to 1932, reconstructed by Nicolas Brehm of ETH Zurich and colleagues from tree-ring carbon-14 measurements. Before feeding the tree-ring record into the inversion, the team cross-calibrated it against the geomagnetic record over their overlapping decades, applying an additive correction of 65.76 megavolts that statistical tests showed produced a near-perfectly symmetric, homoskedastic residual distribution.

The millennial-scale reconstruction that emerged is rich with detail. It captures the familiar grand minima and maxima of the past thousand years, including the Maunder minimum, the Spörer minimum, the Dalton minimum, and the double-peaked grand maximum the Sun passed through between 1900 and 2020. Crucially, the reconstructed open solar flux never falls to zero. Even during the deepest quiet of the Spörer minimum, the unsigned open flux dropped only to about 1.21 times ten to the fourteen webers in 1443, the lowest value in the entire record and well below anything observed in the telescopic era, yet still decisively nonzero. During the Maunder minimum the flux averaged around 2.63 times ten to the fourteen webers. This confirms that the solar dynamo never fully shuts down during grand minima but instead idles in a reduced, finite activity state, a conclusion consistent with flux transport dynamo models in which meridional plasma circulation sustains weak cycles and eventually drives recovery.

The new record also documents some striking extremes at the other end of the scale. The peak annual open solar flux of the twentieth century, reached in 1958, was about 10.96 times ten to the fourteen webers, the highest since the year 1200. Only three years in the entire millennium exceeded it, the largest peaking at roughly 14.21 times ten to the fourteen webers in 981 CE, shortly after the record begins. In other words, the era of telescopic observation has sampled a large fraction, but not all, of the Sun’s true dynamic range; deeper minima and higher maxima both occurred before instruments existed to see them. The reconstruction also flags three intervals around 993, 1052 and 1279 CE where proposed Miyake events, extreme solar particle storms recorded as abrupt radiocarbon spikes, contaminate the cosmic-ray-based record, and the authors conservatively mask these windows from their results.

Beyond its intrinsic appeal as a thousand-year diary of solar magnetism, the work has immediate practical value. Annually resolved, physically consistent sunspot numbers feed directly into reconstructions of total and spectral solar irradiance, which in turn constrain climate models exploring the Sun’s role in terrestrial temperature variability over past centuries. The uncertainty-quantified cycle amplitudes during grand minima provide empirical targets for solar dynamo theorists probing the minimum operating point of the solar cycle engine. The authors note that their forward model templates were built from modern, regular solar cycles, so reconstructions of grand minimum cycles should be treated as indicative rather than definitive, and that the additive cross-calibration between the two modulation potential datasets is itself a simplification. Even so, the framework offers what the team describes as a probabilistic ensemble of physically admissible solar histories rather than a single deterministic answer, and it opens the door to pushing the same technique further back in time as longer and older cosmogenic isotope records become available. The Sun, it turns out, kept meticulous records all along; the trick was learning to read them without breaking the laws of physics.

Subject of Research: Physics-constrained reconstruction of annually resolved sunspot numbers from millennial-scale heliospheric modulation potential records.

Article Title: Physics-Constrained Reconstructions of Sunspot Number from Millennial-Scale Annual Heliospheric Modulation Potential

Article References: Saha, C., Owens, M., Lockwood, M., Barnard, L., Brehm, N., Dalla, S., Herbst, K., Muscheler, R., & Wang, J. (2026). Physics-Constrained Reconstructions of Sunspot Number from Millennial-Scale Annual Heliospheric Modulation Potential. Solar Physics, 301(9), Article 144. https://doi.org/10.1007/s11207-026-02731-0

Image Credits: AI Generated

DOI: 10.1007/s11207-026-02731-0

Keywords: sunspot number, solar cycle, heliospheric modulation potential, cosmogenic isotopes, radiocarbon, open solar flux, Maunder minimum, Spörer minimum, approximate Bayesian computation, space climate, solar dynamo, galactic cosmic rays

Cross-helicity is a measure of the correlation between the velocity and magnetic fluctuations in a turbulent conducting fluid. Unlike the more familiar magnetic helicity, which quantifies the knottedness and linkage of magnetic field lines and has been studied since the foundational work of H.K. Moffatt in the late 1960s, cross-helicity captures how intimately the swirling plasma motion is coupled to the magnetic structures it carries. When velocity and magnetic fluctuations are strongly correlated, cross-helicity becomes large, and its effects on transport, diffusion, and dynamo action can be profound. Recent reviews by Nobumitsu Yokoi have highlighted what he calls unappreciated cross-helicity effects in plasma physics, including anti-diffusion effects in both dynamo action and momentum transport, hinting that this quantity could reshape how physicists model magnetized turbulence.

The new analysis is grounded in a Kolmogorov-like phenomenology developed within the framework of what Bershadskii terms distributed chaos. In classical turbulence theory, the energy spectrum of a turbulent flow follows well-known power-law scaling relations, first articulated by Andrei Kolmogorov in 1941. But chaotic systems, including turbulent magnetohydrodynamic flows, can produce frequency power spectra whose decaying parts follow exponential rather than power-law forms. Exponential power spectra have been observed in a striking variety of chaotic settings, from magnetically confined laboratory plasmas studied by J.E. Maggs and G.J. Morales, to time series of sunspot numbers examined as far back as R.N. Bracewell’s pioneering analyses in the 1950s. In the distributed chaos picture, the exponential decay of the spectrum encodes the underlying deterministic but unpredictable dynamics of the system.

Bershadskii demonstrates, using results from laboratory dynamo experiments and direct numerical simulations of magnetohydrodynamic turbulence, that when the magnetohydrodynamic turbulence is sufficiently strong in chaotic, swirling flows, cross-helicity comes to dominate the decaying portion of the frequency power spectra of the generated magnetic field. The simulations draw on sparse-mode spectral methods of the kind developed by M. Meneguzzi and colleagues, while the laboratory context connects to experiments such as the Taylor-Green dynamo, in which bistability and chaos have already been documented by R.K. Yadav and collaborators. The essential message is that the interaction between swirl and field, quantified by cross-helicity, is not a minor correction to solar dynamo theory but a leading-order feature when turbulence is vigorous.

To test whether this picture applies to the real Sun, the study turns to observational records of the full-disc solar magnetic field and to the international sunspot number series maintained by the Solar Influences Data Analysis Center in Belgium, along with mean-field measurements from Stanford’s Wilcox Solar Observatory and vector observations from the SOLIS facility. The past two solar cycles have been notably weak in magnetic activity, a period that researchers such as K. Mursula and colleagues have described as a transition away from the Modern Maximum toward a weaker Sun. Bershadskii finds that the full-disc magnetic field during these cycles exhibits deterministic chaotic behavior, with the chaotic dynamics concentrated around the solar equator, precisely the region where sunspot activity emerges and migrates in the familiar butterfly diagram.

The finding that the field’s dynamics are concentrated near the equator resonates with independent evidence about the character of solar interior turbulence. Helioseismic analysis by S. Hanasoge, H. Hotta, and K.R. Sreenivasan has shown that turbulence in the Sun is suppressed on large scales and confined to equatorial regions, a result that surprised much of the community. A dynamo operating in a chaotic, swirling, equatorially concentrated flow environment is exactly the kind of system in which cross-helicity dominance of the spectral decay would be expected, according to the theoretical framework of the new paper. The consistency between theory, simulation, experiment, and observation strengthens the case that the distributed chaos approach captures something real about the solar dynamo.

Deterministic chaos is a concept with a long scientific lineage, tracing back to Edward Lorenz’s celebrated 1963 discovery of deterministic nonperiodic flow in atmospheric convection. A chaotic system obeys definite physical laws, yet its future evolution becomes practically unpredictable beyond a certain horizon because tiny uncertainties grow exponentially. The solar cycle, with its roughly eleven-year rhythm, its irregular amplitudes, and its occasional deep minima, has long tempted scientists to search for such behavior. Studies by C. Letellier and colleagues found evidence for low-dimensional chaos in sunspot cycles, while N. Platt, E.A. Spiegel, and C. Tresser modeled the intermittent character of the cycle, and P. Mininni and collaborators developed stochastic relaxation oscillator models. The new work adds a spectral fingerprint, the exponential decay associated with distributed chaos and cross-helicity, to this body of evidence.

The technical core of the argument lies in how the decaying parts of frequency power spectra distinguish between different underlying regimes. Purely stochastic processes tend to produce power-law spectra, whereas deterministic chaos in turbulent flows produces exponential spectral decay whose slope can be predicted by the Kolmogorov-like phenomenology of distributed chaos. When Bershadskii applies this diagnostic to the full-disc solar magnetic field data from the last two weak cycles, the observed spectra match the theoretical expectations for a system in which cross-helicity dominates. In other words, the correlation between the Sun’s plasma motions and its magnetic fluctuations is strong enough that it governs how magnetic energy is distributed across frequencies, a conclusion with direct implications for how dynamo models should be formulated and how the limits of solar activity prediction should be understood.

The implications extend beyond academic curiosity. Solar magnetic activity drives space weather that can damage satellites, disrupt radio communications, and stress terrestrial power grids. Forecasting the amplitude and timing of future cycles is a major practical goal, and approaches based on mean-field dynamo theory continue to be refined by groups such as V. Obridko and colleagues. If the full-disc magnetic field is genuinely governed by deterministic chaos with cross-helicity-dominated spectral dynamics, then there are hard, quantifiable limits to how far ahead reliable prediction can reach, however good the models become. At the same time, recognizing the system as chaotic rather than purely stochastic opens the possibility of short-term predictability techniques borrowed from nonlinear dynamics, exploiting the deterministic structure that chaos still preserves.

There is also a broader astrophysical payoff. Magnetic cycles and dynamos are observed on stars throughout the galaxy, and the physics uncovered in the solar context applies to magnetized plasmas everywhere, from planetary interiors to accretion disks and the interstellar medium. The study’s synthesis, linking swirling MHD turbulence, laboratory dynamo experiments, numerical simulation, and a half-century of solar monitoring, illustrates how a single theoretical lens, the distributed chaos framework enriched by cross-helicity, can bring coherence to phenomena that have often been treated separately. As the Sun continues its journey away from the Modern Maximum into a period of weaker activity, the chaotic, equator-concentrated dynamics identified in this work offer both a warning about the limits of prediction and an invitation to look more closely at the turbulent, helical plasma churning beneath the solar surface.

Subject of Research: Cross-helicity and deterministic chaotic dynamics of the full-disc solar magnetic field generated by a turbulent magnetohydrodynamic dynamo.

Article Title: Cross-Helicity and Chaotic Dynamics of Full-Disc Solar Magnetic Field

Article References: Bershadskii, A. (2026). Cross-Helicity and Chaotic Dynamics of Full-Disc Solar Magnetic Field. Solar Physics, 301(9), Article 141. https://doi.org/10.1007/s11207-026-02725-y

Image Credits: AI Generated

DOI: 10.1007/s11207-026-02725-y

Keywords: solar dynamo, cross-helicity, magnetohydrodynamic turbulence, distributed chaos, full-disc magnetic field, sunspot number, solar cycle, deterministic chaos, power spectra, Kolmogorov phenomenology, swirling flows, solar physics