Open Solar Flux Research Articles

Assessment of different sunspot number series using the cosmogenic isotope 44 Ti in meteorites

Many sunspot number series exist suggesting different levels of solar activity during the past centuries. Their reliability can be assessed only by comparing them with alternative indirect proxies. We test different sunspot number series against the updated record of cosmogenic radionuclide 44Ti measured in meteorites. Two bounding scenarios of solar activity changes have been considered: the HH-scenario (based on the series by Svalgaard and Schatten), in particular, predicting moderate activity during the Maunder minimum, and the LL-scenario (based on the RG series by Lockwood et al.) predicting moderate activity for the 18th–19th centuries and the very low activity level for the Maunder minimum. For each scenario, the magnetic open solar flux, the heliospheric modulation potential and the expected production of 44Ti were computed. The calculated production rates were compared with the corresponding measurements of 44Ti activity in stony meteorites fallen since 1766. The analysis reveals that the LL-scenario is fully consistent with the measured 44Ti data, in particular, recovering the observed secular trend between the 17th century and the Modern grand maximum. On the contrary, the HH-scenario appears significantly inconsistent with the data, mostly due to the moderate level of activity during the Maunder minimum. It is concluded that the HH-scenario sunspot number reconstruction significantly overestimates solar activity prior to the mid-18th century, especially during the Maunder minimum. The exact level of solar activity after 1750 cannot be distinguished with this method, since both H- and L- scenarios appear statistically consistent with the data.

Tests of Sunspot Number Sequences: 2. Using Geomagnetic and Auroral Data

We compare four sunspot-number data sequences against geomagnetic and terrestrial auroral observations. The comparisons are made for the original SIDC composite of Wolf-Zurich-International sunspot number [$R_{ISNv1}$], the group sunspot number [$R_{G}$] by Hoyt and Schatten (Solar Phys., 1998), the new "backbone" group sunspot number [$R_{BB}$] by Svalgaard and Schatten (Solar Phys., 2016), and the "corrected" sunspot number [$R_{C}$] by Lockwood at al. (J.G.R., 2014). Each sunspot number is fitted with terrestrial observations, or parameters derived from terrestrial observations to be linearly proportional to sunspot number, over a 30-year calibration interval of 1982-2012. The fits are then used to compute test sequences, which extend further back in time and which are compared to $R_{ISNv1}$, $R_{G}$, $R_{BB}$, and $R_{C}$. To study the long-term trends, comparisons are made using averages over whole solar cycles (minimum-to-minimum). The test variations are generated in four ways: i) using the IDV(1d) and IDV geomagnetic indices (for 1845-2013) fitted over the calibration interval using the various sunspot numbers and the phase of the solar cycle; ii) from the open solar flux (OSF) generated for 1845 - 2013 from four pairings of geomagnetic indices by Lockwood et al. (Ann. Geophys., 2014) and analysed using the OSF continuity model of Solanki at al. (Nature, 2000) which employs a constant fractional OSF loss rate; iii) the same OSF data analysed using the OSF continuity model of Owens and Lockwood (J.G.R., 2012) in which the fractional loss rate varies with the tilt of the heliospheric current sheet and hence with the phase of the solar cycle; iv) the occurrence frequency of low-latitude aurora for 1780-1980 from the survey of Legrand and Simon (Ann. Geophys., 1987). For all cases, $R_{BB}$ exceeds the test terrestrial series by an amount that increases as one goes back in time.

AN ASSESSMENT OF SUNSPOT NUMBER DATA COMPOSITES OVER 1845–2014

New sunspot data composites, some of which are radically different in the character of their long-term variation, are evaluated over the interval 1845-2014. The method commonly used to calibrate historic sunspot data, relative to modern-day data, is "daisy-chaining", whereby calibration is passed from one data subset to the neighbouring one, usually using regressions of the data subsets for the intervals of their overlap. Recent studies have illustrated serious pitfalls in these regressions and the resulting errors can be compounded by their repeated use as the data sequence is extended back in time. Hence the recent composite data series by Usoskin et al. (2016), $R_{UEA}$, is a very important advance because it avoids regressions, daisy-chaining and other common, but invalid, assumptions: this is achieved by comparing the statistics of "active day" fractions to those for a single reference dataset. We study six sunspot data series including $R_{UEA}$ and the new "backbone" data series $R_{BB}$, recently generated by Svalgaard and Schatten (2016) by employing both regression and daisy-chaining. We show that all six can be used with a continuity model to reproduce the main features of the open solar flux variation for 1845-2014, as reconstructed from geomagnetic activity data. However, some differences can be identified that are consistent with tests using a basket of other proxies for solar magnetic fields. Using data from a variety of sunspot observers, we illustrate problems with the method employed in $R_{BB}$ which cause it to increasingly overestimate sunspot numbers going back in time and we recommend using $R_{UEA}$ because it employs more robust procedures that avoid such problems.

COMPARISON OF CORONAL EXTRAPOLATION METHODS FOR CYCLE 24 USING HMI DATA

ABSTRACT Two extrapolation models of the solar coronal magnetic field are compared using magnetogram data from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager instrument. The two models, a horizontal current–current sheet–source surface (HCCSSS) model and a potential field–source surface (PFSS) model, differ in their treatment of coronal currents. Each model has its own critical variable, respectively, the radius of a cusp surface and a source surface, and it is found that adjusting these heights over the period studied allows for a better fit between the models and the solar open flux at 1 au as calculated from the Interplanetary Magnetic Field (IMF). The HCCSSS model provides the better fit for the overall period from 2010 November to 2015 May as well as for two subsets of the period: the minimum/rising part of the solar cycle and the recently identified peak in the IMF from mid-2014 to mid-2015 just after solar maximum. It is found that an HCCSSS cusp surface height of 1.7 R ⊙ provides the best fit to the IMF for the overall period, while 1.7 and 1.9 R ⊙ give the best fits for the two subsets. The corresponding values for the PFSS source surface height are 2.1, 2.2, and 2.0 R ⊙ respectively. This means that the HCCSSS cusp surface rises as the solar cycle progresses while the PFSS source surface falls.

An empirical model of heliospheric cosmic ray modulation on long-term time scale

Galactic Cosmic Rays (GCRs) entering the heliosphere are subject to modulation processes due to variable solar magnetic activity. Finding a relationship between cosmic-ray variations and the heliospheric parameters is important for reconstruction of solar activity in the past. Here, we develop a semi-empirical model describing the heliospheric modulation of GCRs in terms of heliospheric parameters such as the open solar magnetic flux, the tilt angle of the heliospheric current sheet and the polarity of the large scale solar magnetic field. Our model is fitted using annual data obtained for the period 1976–2013, which includes the very weak solar minimum during 2008–2010. The model shows a good agreement with the data, and therefore, can be used for reconstructions of the modulation potential at different levels of solar activity. The model’s validity is also tested using the cosmogenic radionuclides 14 C and 10 Be stored in terrestrial archives. The tilt angle used to fit the parameters in our semi-empirical modulation model is reconstructed by a mathematical model described here.

NEAR-EARTH COSMIC RAY DECREASES ASSOCIATED WITH REMOTE CORONAL MASS EJECTIONS

Galactic cosmic ray (GCR) flux is modulated by both particle drift patterns and solar wind structures on a range of timescales. Over solar cycles, GCR flux varies as a function of the total open solar magnetic flux and the latitudinal extent of the heliospheric current sheet. Over hours, drops of a few percent in near-Earth GCR flux (Forbush decreases, FDs) are well known to be associated with the near-Earth passage of solar wind structures resulting from corotating interaction regions (CIRs) and transient coronal mass ejections (CMEs). We report on four FDs seen at ground-based neutron monitors which cannot be immediately associated with significant structures in the local solar wind. Similarly, there are significant near-Earth structures which do not produce any corresponding GCR variation. Three of the FDs are during the STEREO era, enabling in situ and remote observations from three well-separated heliospheric locations. Extremely large CMEs passed the STEREO-A spacecraft, which was behind the West limb of the Sun, approximately 2-3 days before each near-Earth FD. Solar wind simulations suggest that the CMEs combined with pre-existing CIRs, enhancing the pre-existing barriers to GCR propagation. Thus these observations provide strong evidence for the modulation of GCR flux by remote solar wind structures.

The heliospheric Hale cycle over the last 300 years and its implications for a “lost” late 18th century solar cycle

A Hale cycle, one complete magnetic cycle of the Sun, spans two complete Schwabe cycles (also referred to as sunspot and, more generally, solar cycles). The approximately 22-year Hale cycle is seen in magnetic polarities of both sunspots and polar fields, as well as in the intensity of galactic cosmic rays reaching Earth, with odd- and even-numbered solar cycles displaying qualitatively different waveforms. Correct numbering of solar cycles also underpins empirical cycle-to-cycle relations which are used as first-order tests of stellar dynamo models. There has been much debate about whether the unusually long solar cycle 4 (SC4), spanning 1784–1799, was actually two shorter solar cycles combined as a result of poor data coverage in the original Wolf sunspot number record. Indeed, the group sunspot number does show a small increase around 1794–1799 and there is evidence of an increase in the mean latitude of sunspots at this time, suggesting the existence of a cycle “4b”. In this study, we use cosmogenic radionuclide data and associated reconstructions of the heliospheric magnetic field (HMF) to show that the Hale cycle has persisted over the last 300 years and that data prior to 1800 are more consistent with cycle 4 being a single long cycle (the “no SC4b” scenario). We also investigate the effect of cycle 4b on the HMF using an open solar flux (OSF) continuity model, in which the OSF source term is related to sunspot number and the OSF loss term is determined by the heliospheric current sheet tilt, assumed to be a simple function of solar cycle phase. The results are surprising; Without SC4b, the HMF shows two distinct peaks in the 1784–1799 interval, while the addition of SC4b removes the secondary peak, as the OSF loss term acts in opposition to the later rise in sunspot number. The timing and magnitude of the main SC4 HMF peak is also significantly changed by the addition of SC4b. These results are compared with the cosmogenic isotope reconstructions of HMF and historical aurora records. These data marginally favour the existence of SC4b (the “SC4b” scenario), though the result is less certain than that based on the persistence of the Hale cycle. Thus while the current uncertainties in the observations preclude any definitive conclusions, the data favour the “no SC4b” scenario. Future improvements to cosmogenic isotope reconstructions of the HMF, through either improved modelling or additional ice cores from well-separated geographic locations, may enable questions of the existence of SC4b and the phase of Hale cycle prior to the Maunder minimum to be settled conclusively.

Interannual to multidecadal Euro‐Atlantic blocking variability during winter and its relationship with extreme low temperatures in Europe

AbstractThe dominant modes of blocking frequency variability in the Atlantic‐European region are evaluated for the 1871–2010 period. An Empirical Orthogonal Function (EOF) analysis of a two‐dimensional blocking indicator field reveals three dominant EOFs, describing about 35% of interannual to multidecadal blocking variability. The first EOF captures an out‐of‐phase blocking frequency anomaly over Greenland and Western Europe regions. The corresponding principal component time series is strongly correlated with the North Atlantic Oscillation index but shows also significant correlations with indices of the East Atlantic, Scandinavian, and East Atlantic‐Western Russia patterns. The second EOF shows a dominant center over the North Sea region as well as a less pronounced center with anomalies of the same sign over southeastern Greenland. The multidecadal variations of this mode of blocking variability are related with a basin wide North Atlantic sea surface temperature anomaly which projects partly on the Atlantic Multidecadal Oscillation (AMO). The third mode is an east‐west dipole of blocking frequency anomalies from Scandinavian and southern Greenland regions and shows enhanced variability at ~20 year time scales. The coherent variations of the time coefficients of this pattern with open solar flux suggest a possible solar influence on blocking variability at these time scales. Furthermore, the dominant patterns of blocking variability are related with distinct anomaly patterns in the occurrence of extreme low temperature events over Europe at interannual to multidecadal time scales. AMO as well as the solar signals was detected also in the corresponding extreme low temperature blocking patterns. We argue that multivariate analysis of blocking indicators gives additional information about blocking and related extreme climate phenomena variability and predictability comparative with classical sectorial approach.

Variations of solar, interplanetary, and geomagnetic parameters with solar magnetic multipole fields during Solar Cycles 21–24

In this study we compare the temporal variations of the solar, interplanetary, and geomagnetic (SIG) parameters with that of open solar magnetic flux from 1976 to 2012 (from Solar Cycle 21 to the early phase of Cycle 24) for a purpose of identifying their possible relationships. By the open flux, we mean the average magnetic field over the source surface (2.5 solar radii) times the source area as defined by the potential field source surface (PFSS) model of the Wilcox Solar Observatory (WSO). In our result, most SIG parameters except the solar wind dynamic pressure show rather poor correlations with the open solar magnetic field. Good correlations are recovered when the contributions from individual multipole components are counted separately. As expected, solar activity indices such as sunspot number, total solar irradiance, 10.7cm radio flux, and solar flare occurrence are highly correlated with the flux of magnetic quadrupole component. The dynamic pressure of solar wind is strongly correlated with the dipole flux, which is in anti-phase with Solar Cycle (SC). The geomagnetic activity represented by the Ap index is correlated with higher order multipole components, which show relatively a slow time variation with SC. We also found that the unusually low geomagnetic activity during SC 23 is accompanied by the weak open solar fields compared with those in other SCs. It is argued that such dependences of the SIG parameters on the individual multipole components of the open solar magnetic flux may clarify why some SIG parameters vary in phase with SC and others show seemingly delayed responses to SC variation.

Centennial variations in sunspot number, open solar flux and streamer belt width: 3. Modeling

AbstractFrom the variation of near‐Earth interplanetary conditions, reconstructed for the mid‐19th century to the present day using historic geomagnetic activity observations, Lockwood and Owens (2014) have suggested that Earth remains within a broadened streamer belt during solar cycles when the Open Solar Flux (OSF) is low. From this they propose that the Earth was immersed in almost constant slow solar wind during the Maunder minimum (c. 1650–1710). In this paper, we extend continuity modeling of the OSF to predict the streamer belt width using both group sunspot numbers and corrected international sunspot numbers to quantify the emergence rate of new OSF. The results support the idea that the solar wind at Earth was persistently slow during the Maunder minimum because the streamer belt was broad.

Centennial variations in sunspot number, open solar flux, and streamer belt width: 2. Comparison with the geomagnetic data

AbstractWe investigate the relationship between interdiurnal variation geomagnetic activity indices, IDV and IDV(1d), corrected sunspot number, RC , and the group sunspot number RG . RC uses corrections for both the “Waldmeier discontinuity,” as derived in Paper 1, and the “Wolf discontinuity” revealed by Leussu et al. (2013). We show that the simple correlation of the geomagnetic indices with RCn or RGn masks a considerable solar cycle variation. Using IDV(1d) or IDV to predict or evaluate the sunspot numbers, the errors are almost halved by allowing for the fact that the relationship varies over the solar cycle. The results indicate that differences between RC and RG have a variety of causes and are highly unlikely to be attributable to errors in either RG alone, as has recently been assumed. Because it is not known if RC or RG is a better predictor of open flux emergence before 1874, a simple sunspot number composite is suggested which, like RG , enables modeling of the open solar flux for 1610 onward in Paper 3 but maintains the characteristics of RC .

Centennial variations in sunspot number, open solar flux, and streamer belt width: 1. Correction of the sunspot number record since 1874

AbstractWe analyze the widely used international/Zürich sunspot number record, R, with a view to quantifying a suspected calibration discontinuity around 1945 (which has been termed the “Waldmeier discontinuity”). We compare R against the composite sunspot group data from the Royal Greenwich Observatory network and the Solar Optical Observing Network, using both the number of sunspot groups, NG, and the total area of the sunspots, AG. In addition, we compare R with the recently developed interdiurnal variability geomagnetic indices IDV and IDV(1d). In all four cases, linearity of the relationship with R is not assumed and care is taken to ensure that the relationship of each with R is the same before and after the putative calibration change. It is shown the probability that a correction is not needed is of order 10−8 and that R is indeed too low before 1945. The optimum correction to R for values before 1945 is found to be 11.6%, 11.7%, 10.3%, and 7.9% using AG, NG, IDV, and IDV(1d), respectively. The optimum value obtained by combining the sunspot group data is 11.6% with an uncertainty range 8.1–14.8% at the 2σ level. The geomagnetic indices provide an independent yet less stringent test but do give values that fall within the 2σ uncertainty band with optimum values are slightly lower than from the sunspot group data. The probability of the correction needed being as large as 20%, as advocated by Svalgaard (2011), is shown to be 1.6 × 10−5.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr – Part 4: Near-Earth solar wind speed, IMF, and open solar flux

Abstract. In the concluding paper of this tetralogy, we here use the different geomagnetic activity indices to reconstruct the near-Earth interplanetary magnetic field (IMF) and solar wind flow speed, as well as the open solar flux (OSF) from 1845 to the present day. The differences in how the various indices vary with near-Earth interplanetary parameters, which are here exploited to separate the effects of the IMF and solar wind speed, are shown to be statistically significant at the 93% level or above. Reconstructions are made using four combinations of different indices, compiled using different data and different algorithms, and the results are almost identical for all parameters. The correction to the aa index required is discussed by comparison with the Ap index from a more extensive network of mid-latitude stations. Data from the Helsinki magnetometer station is used to extend the aa index back to 1845 and the results confirmed by comparison with the nearby St Petersburg observatory. The optimum variations, using all available long-term geomagnetic indices, of the near-Earth IMF and solar wind speed, and of the open solar flux, are presented; all with ±2σ uncertainties computed using the Monte Carlo technique outlined in the earlier papers. The open solar flux variation derived is shown to be very similar indeed to that obtained using the method of Lockwood et al. (1999).

MAGNETIC FLUX DENSITY IN THE HELIOSPHERE THROUGH SEVERAL SOLAR CYCLES

We studied the magnetic flux density carried by solar wind to various locations in the heliosphere, covering a heliospheric distance range of 0.3-5.4 AU and a heliolatitudinal range from 80° south to 80° north. Distributions of the radial component of the magnetic field, B{sub R} , were determined over long intervals from the Helios, ACE, STEREO, and Ulysses missions, as well as from using the 1 AU OMNI data set. We show that at larger distances from the Sun, the fluctuations of the magnetic field around the average Parker field line distort the distribution of B{sub R} to such an extent that the determination of the unsigned, open solar magnetic flux density from the average (|B{sub R} |) is no longer justified. We analyze in detail two methods for reducing the effect of fluctuations. The two methods are tested using magnetic field and plasma velocity measurements in the OMNI database and in the Ulysses observations, normalized to 1 AU. It is shown that without such corrections for the fluctuations, the magnetic flux density measured by Ulysses around the aphelion phase of the orbit is significantly overestimated. However, the matching between the in-ecliptic magnetic flux density at 1 AU (OMNI data)more » and the off-ecliptic, more distant, normalized flux density by Ulysses is remarkably good if corrections are made for the fluctuations using either method. The main finding of the analysis is that the magnetic flux density in the heliosphere is fairly uniform, with no significant variations having been observed either in heliocentric distance or heliographic latitude.« less

Solar cycle evolution of dipolar and pseudostreamer belts and their relation to the slow solar wind

Dipolar streamers are coronal structures formed by open solar flux converging from coronal holes of opposite polarity. Thus, the dipolar streamer belt traces the coronal foot point of the heliospheric current sheet, and it is strongly associated with the origin of slow solar wind. Pseudostreamers, on the other hand, separate converging regions of open solar flux from coronal holes of the same polarity and do not contain current sheets. They have recently received a great deal of interest as a possible additional source of slow solar wind. Here we add to that growing body of work by using the potential‐field source‐surface model to determine the occurrence and location of dipolar and pseudostreamers over the last three solar cycles. In addition to providing new information about pseudostreamer morphology, the results help explain why the observations taken during the first Ulysses perihelion pass in 1995 showed noncoincidence between dipolar streamer belt and the locus of slowest flow. We find that Carrington rotation averages of the heliographic latitudes of dipolar and pseudostreamer belts are systematically shifted away from the equator, alternately in opposite directions, with a weak solar cycle periodicity, thus keeping slow wind from the web of combined streamer belts approximately symmetric about the equator. The largest separation of dipolar and pseudostreamer belts occurred close to the Ulysses pass, allowing a unique opportunity to see that slow wind from pseudostreamer belts north of the southward displaced dipolar belt was responsible for the noncoincident pattern.

M-dwarf stellar winds: the effects of realistic magnetic geometry on rotational evolution and planets

We perform three-dimensional numerical simulations of stellar winds of early-M dwarf stars. Our simulations incorporate observationally reconstructed large-scale surface magnetic maps, suggesting that the complexity of the magnetic field can play an important role in the angular momentum evolution of the star, possibly explaining the large distribution of periods in field dM stars, as reported in recent works. In spite of the diversity of the magnetic field topologies among the stars in our sample, we find that stellar wind flowing near the (rotational) equatorial plane carries most of the stellar angular momentum, but there is no preferred colatitude contributing to mass loss, as the mass flux is maximum at different colatitudes for different stars. We find that more non-axisymmetric magnetic fields result in more asymmetric mass fluxes and wind total pressures $p_{\rm tot}$ (defined as the sum of thermal, magnetic and ram pressures). Because planetary magnetospheric sizes are set by pressure equilibrium between the planet's magnetic field and $p_{\rm tot}$, variations of up to a factor of $3$ in $p_{\rm tot}$ (as found in the case of a planet orbiting at several stellar radii away from the star) lead to variations in magnetospheric radii of about 20 percent along the planetary orbital path. In analogy to the flux of cosmic rays that impact the Earth, which is inversely modulated with the non-axisymmetric component of the total open solar magnetic flux, we conclude that planets orbiting M dwarf stars like DT~Vir, DS~Leo and GJ~182, which have significant non-axisymmetric field components, should be the more efficiently shielded from galactic cosmic rays, even if the planets lack a protective thick atmosphere/large magnetosphere of their own.

Solar origin of heliospheric magnetic field inversions: Evidence for coronal loop opening within pseudostreamers

The orientation of the heliospheric magnetic field (HMF) in near‐Earth space is generally a good indicator of the polarity of HMF foot points at the photosphere. There are times, however, when the HMF folds back on itself (is inverted), as indicated by suprathermal electrons locally moving sunward, even though they must ultimately be carrying the heat flux away from the Sun. Analysis of the near‐Earth solar wind during the period 1998–2011 reveals that inverted HMF is present approximately 5.5% of the time and is generally associated with slow, dense solar wind and relatively weak HMF intensity. Inverted HMF is mapped to the coronal source surface, where a new method is used to estimate coronal structure from the potential‐field source‐surface model. We find a strong association with bipolar streamers containing the heliospheric current sheet, as expected, but also with unipolar or pseudostreamers, which contain no current sheet. Because large‐scale inverted HMF is a widely accepted signature of interchange reconnection at the Sun, this finding provides strong evidence for models of the slow solar wind which involve coronal loop opening by reconnection within pseudostreamer belts as well as the bipolar streamer belt. Occurrence rates of bipolar‐ and pseudostreamers suggest that they are equally likely to result in inverted HMF and, therefore, presumably undergo interchange reconnection at approximately the same rate. Given the different magnetic topologies involved, this suggests the rate of reconnection is set externally, possibly by the differential rotation rate which governs the circulation of open solar flux.

Reconstruction and Prediction of Variations in the Open Solar Magnetic Flux and Interplanetary Conditions

Historic geomagnetic activity observations have been used to reveal centennial variations in the open solar flux and the near-Earth heliospheric conditions (the interplanetary magnetic field and the solar wind speed). The various methods are in very good agreement for the past 135 years when there were sufficient reliable magnetic observatories in operation to eliminate problems due to site-specific errors and calibration drifts. This review underlines the physical principles that allow these reconstructions to be made, as well as the details of the various algorithms employed and the results obtained. Discussion is included of: the importance of the averaging timescale; the key differences between “range” and “interdiurnal variability” geomagnetic data; the need to distinguish source field sector structure from heliosphericallyimposed field structure; the importance of ensuring that regressions used are statistically robust; and uncertainty analysis. The reconstructions are exceedingly useful as they provide calibration between the in-situ spacecraft measurements from the past five decades and the millennial records of heliospheric behaviour deduced from measured abundances of cosmogenic radionuclides found in terrestrial reservoirs. Continuity of open solar flux, using sunspot number to quantify the emergence rate, is the basis of a number of models that have been very successful in reproducing the variation derived from geomagnetic activity. These models allow us to extend the reconstructions back to before the development of the magnetometer and to cover the Maunder minimum. Allied to the radionuclide data, the models are revealing much about how the Sun and heliosphere behaved outside of grand solar maxima and are providing a means of predicting how solar activity is likely to evolve now that the recent grand maximum (that had prevailed throughout the space age) has come to an end.

Heliospheric modulation of galactic cosmic rays during grand solar minima: Past and future variations

Galactic cosmic ray flux at Earth is modulated by the heliospheric magnetic field. Heliospheric modulation potential, Φ, during grand solar minima is investigated using an open solar flux (OSF) model with OSF source based on sunspot number, R, and OSF loss on heliospheric current sheet inclination. Changing dominance between source and loss means Φ varies in‐ (anti‐) phase with R during strong (weak) cycles, in agreement with Φ estimates from ice core records of10Be concentration, which are in‐phase during most of the last 300 years, but anti‐phase during the Maunder Minimum. Model results suggest “flat” OSF cycles, such as solar cycle 20 result from OSF source and loss terms temporarily balancing throughout the cycle. Thus even if solar activity continues to decline steadily, the long‐term drop in OSF through SC21 to SC23 may plateau during SC24, though reemerge in SC25 with the inverted phase relation.

Solar-cycle precursors and predictions

AbstractThe sunspot number data during the past 400 years indicates that both the profile and the amplitude of the solar cycle have large variations. Some precursors of the solar cycle were identified aiming to predict the solar cycle. The polar field and the geomagnetic index are two precursors which are received the most attention. The geomagnetic variations during the solar minima are potentially caused by the solar polar field by the connection of the solar open flux. The robust prediction skill of the polar field indicates that the memory of the dynamo process is less than 11 yrs based on the frame of the Babcock-Leighton flux transport dynamo. One possible reason to get the short magnetic memory is the high magnetic diffusivity in the convective zone. Our recent studies show that the radial downward pumping is another possible reason. Based upon the mechanism, we well simulate the cycle irregularities during RGO time period. This opens the possibility to set up a standard dynamo based model to predict the solar cycle. In the end, the no correlation between the polar field and the preceding cycle strength due to two nonlinearities involved in the sunspot emergence will be stressed.