Solar-Wind Density Fluctuations Mapped Worldwide via Scintillation Data, 2008–2022

Billions of tons of plasma leave the Sun every hour, streaming outward as the solar wind and inflating a magnetic cocoon—the heliosphere—that encloses every planet in the solar system. That invisible ocean is neither calm nor constant. Embedded within it are ripples of electron density that swell and ebb with the Sun’s eleven-year activity cycle, and a new fifteen-year study has now charted those ripples with unprecedented continuity. By watching the radio emission from dozens of distant quasars flicker as their signals skim past the Sun, Munetoshi Tokumaru of the Institute for Space–Earth Environmental Research at Nagoya University has mapped the global distribution of solar-wind density fluctuations from the deep minimum between Solar Cycles 23 and 24, through the maximum of Cycle 24, and into the rising phase of Cycle 25. The results, published in Solar Physics, capture the heliosphere breathing with the solar cycle—and reveal that the recent, unusually weak cycles thinned it in a historically distinctive way.

The technique behind the measurement, known as interplanetary scintillation, or IPS, is elegant in its simplicity. When radio waves from a compact, point-like source such as a quasar sweep through the solar wind, they encounter countless small-scale irregularities in electron density. Each irregularity scatters the waves slightly, and the overlapping scattered wavelets interfere with one another, producing a rapidly shifting diffraction pattern that sweeps across Earth. Radio telescopes positioned in the pattern’s path record a characteristic twinkling—a stochastic flutter in signal strength on timescales of a second or so—that carries coded information about the plasma it traversed. The stronger the density fluctuations along the line of sight, and the faster the solar wind carries them past the observer, the more violently the signal flickers. First exploited in the 1960s by Antony Hewish and colleagues, IPS remains one of the only ground-based methods capable of sampling the solar wind close to the Sun and across the full range of heliographic latitudes, something no operating spacecraft has been able to do since the Ulysses mission ended.

The pivotal quantity in the new analysis is the turnover distance. Far from the Sun, density fluctuations are modest, and the phase variations they imprint on passing radio waves are small—well below a radian of phase variance—so scattered waves interfere constructively with the direct wave. In this weak-scattering regime the scintillation index, the fractional strength of the flicker, grows in near-proportion to the fluctuation amplitude as the line of sight approaches the Sun, where fluctuations fall off roughly as the inverse square of radial distance; the observed scintillation slope of about −1.5 in power-law terms matches precisely that uniformly expanding wind once the Fresnel propagation filter is accounted for. But close to the Sun the fluctuations grow so large that the weak-scattering condition breaks: scattered waves can no longer lock constructively onto the direct beam, destructive interference takes over, and the scintillation index collapses instead of climbing. Plotted against the Sun-source offset, the index therefore rises to a peak and falls away. The location of that peak—the turnover distance—marks the weak-to-strong scattering transition and, at a fixed observing frequency, is fixed uniquely by the fluctuation level along the line of view, independent of the source being observed.

Toyokawa Observatory’s 327-megahertz Solar Wind Imaging Facility—SWIFT—is the workhorse behind the dataset. Its phased-array antenna, spanning 106 meters north–south and 38 meters east–west, is the most sensitive element of a three- to four-station network that Nagoya University has operated since the early 1980s; sister telescopes roughly 100 kilometers away record the same sources simultaneously, and cross-correlating their signals yields solar-wind speeds. Each clear day, every program source was observed around meridian transit, at solar elongation angles below 90 degrees, so that the scintillation index could be sampled as a function of how closely the line of sight passed to the Sun. From the 2008–2022 archive, turnover distances were determined for 36 compact sources lying at a spread of heliographic latitudes. Plotting those distances together traces a contour of constant fluctuation level across the plane of the sky—an annual silhouette of the Sun’s turbulent wind, assembled one quasar at a time.

What makes the turnover distance so valuable is its stubbornness. The scintillation index itself is contaminated by confounding factors—slow degradation of receiver sensitivity over decades, bandwidth effects, and crucially the apparent angular size of each source, since extended sources scintillate more weakly than point-like ones. The turnover position is essentially immune to such effects: the weak-to-strong scattering transition does not depend on source size and survives even when instrumental sensitivity drifts, so measurements from quasars of different brightness and angular extent can be safely combined into one record. That robustness let the study build a coherent fifteen-year dataset and compare it directly with the landmark 327-megahertz IPS surveys made with the Ooty Radio Telescope in India during Cycles 21 and 22, which established that the fluctuation distribution at solar maximum was nearly spherically symmetric—a circle of radius about 46 solar radii—while at minimum it collapsed into an ellipse, 44 solar radii along the equator but only 28 over the poles.

Compiled year by year, the contours show the heliosphere swelling and contracting nearly in lockstep with the sunspot record. Around the maximum of Cycle 24, the equatorial radius of the fluctuation contour peaked at 46 solar radii in 2015 and 2016—precisely the value Ooty had measured at the maxima of Cycles 21 and 22—while the polar radius climbed to about 42 solar radii, a pole-to-equator ratio of roughly 0.9. In 2010 the contour was almost perfectly circular, and a nearly uniform geometry recurred around 2015, mirroring the visible corona at solar maximum, when bright helmet streamers blanket every latitude and the dark polar coronal holes largely vanish. Yearly-mean turnover distances reached 43 to 45 solar radii in 2015 and again in 2022, and both the contour radii and the means lagged the sunspot-number peak by one to two years—evidence that the heliosphere’s response to solar activity is smoothed and delayed by the wind’s own dynamics.

At the minima the picture contracts sharply. In 2008, near the deep Cycle 23/24 minimum, the polar radius of the contour shrank to 28 solar radii, the smallest value in the dataset; the same figure reappeared in 2018 at the Cycle 24/25 minimum, when the equatorial radius had fallen to 38–39 solar radii. Only with the rise of Cycle 25 did the geometry re-inflate, the equatorial radius returning to 46 solar radii in 2022. The persistent polar asymmetry is the fingerprint of the coronal holes: those dark, open-field caps emit fast, comparatively smooth wind carrying fewer density irregularities, depressing the fluctuation level over the poles. Translated into fluctuation amplitudes at roughly 0.2 astronomical units from the Sun, the analysis implies values at solar minimum of about 70 percent of the solar-maximum level at the equator but only about 56 percent at the poles—a latitude-dependent thinning never before tracked so continuously through two consecutive minima.

The historical comparison is where the study delivers its surprise. During the minima of Cycles 20/21 and 21/22, Ooty’s contours formed a pronounced ellipse—44 solar radii along the equator, 28 over the poles—a pole-to-equator ratio of about 0.6, implying that polar fluctuations were suppressed by a factor of roughly 2.5 relative to the equator by the prevailing coronal holes. The recent minima paint a different portrait. At both the Cycle 23/24 and Cycle 24/25 minima, Toyokawa recorded equatorial turnover distances of 38–39 solar radii, well below Ooty’s 44, and polar values of 31–32 solar radii, slightly above Ooty’s 28; the pole-to-equator ratio consequently rose to about 0.8, a distribution still latitude-dependent but markedly more uniform than in earlier cycles. The depressed equatorial values indicate an equatorial solar wind that was significantly rarefied—a hallmark of the anomalously weak Cycle 24. Coherence-bandwidth measurements at Nancy decades earlier had already shown that fluctuations are latitude-invariant at maximum and latitude-dependent at minimum, and the new data refine that picture: the rarefaction was largely confined to the slow equatorial wind, since pulsar soundings through the polar coronal hole in 2018–2019 found polar densities essentially unchanged from the declining phases of Cycles 20 and 22, echoing in situ evidence that only the low-speed wind thinned.

The dataset also yielded an unexpected statistical relationship. Tokumaru characterized the sharpness of each yearly scintillation peak with a curvature coefficient—values between −3 and −5—that reflects the radial gradient of the fluctuation distribution at the turnover distance, and found a significant negative correlation between sharpness and turnover distance: a correlation coefficient of −0.59 with a p-value statistically indistinguishable from zero. In physical terms, when fluctuations are strong and the turnover lies far from the Sun, the peak is sharp and the radial gradient steep; when fluctuations weaken and the turnover moves inward, the peak broadens and the gradient flattens. The relation held for high- and low-latitude sources alike, although at solar minima the high-latitude peaks were systematically broader, signaling that the fluctuation profile depends on latitude precisely when polar coronal holes dominate. A weaker but still significant correlation between peak sharpness and latitude itself (coefficient 0.33, p = 0.003) reinforced that conclusion and carries direct implications for how large-scale density structure forms across the heliosphere.

The practical stakes are considerable. The density of the solar wind and the turbulence it carries govern how coronal mass ejections expand on their way to Earth, how energetic particles are scattered through interplanetary space, and ultimately how geomagnetic storms unfold, so a cycle-resolved picture of the global wind feeds directly into space-weather forecasting. Since Ulysses ended, no spacecraft has sampled the wind over the poles, leaving ground-based scintillation as the only eye on the heliosphere’s full geometry; computer-assisted tomography built on these kinds of measurements already reconstructs three-dimensional density maps and has revealed rarefied corridors within the slow wind. The new record also sets a benchmark for the current cycle: by 2022, turnover distances had rebounded to their solar-maximum extent of 46 solar radii, hinting at a vigorous recovery. Whether the approaching maximum matches or exceeds Cycle 24’s subdued peak, the flicker of distant quasars will keep recording it—one tiny flutter at a time.