Abstract
The open solar flux (OSF) is the integrated unsigned magnetic flux leaving the top of the solar atmosphere to form the heliospheric magnetic field. As the OSF modulates the intensity of galactic cosmic rays at Earth, the production rate of cosmogenic isotopes – such as 14C and 10Be stored in tree rings and ice sheets – is closely related to the OSF. Thus on the basis of cosmogenic isotope data, OSF can be reconstructed over millennia. As sunspots are related to the production of OSF, this provides the possibility of reconstructing sunspot number (SSN) and hence properties of the solar cycles prior to the first sunspot telescopic observations in 1610. However, while models exist for estimating OSF on the basis of SSN, the hysteresis present in OSF and the lack of a priori knowledge of the start/end dates of individual solar cycles means that directly inverting these models is not possible. We here describe a new method that uses a forward model of OSF to estimate SSN and solar cycle start/end dates through a Monte Carlo approach. The method is tested by application to geomagnetic reconstructions of OSF over the period 1845-present, and compared to the known SSN record for this period. There is a substantial improvement in reconstruction of both the SSN time series and the solar cycle start/end dates compared with existing OSF-SSN regression methods. This suggests that more accurate solar-cycle information can be extracted from cosmogenic isotope records by forward modelling, and also provides a means to assess the level of agreement between independent SSN and OSF reconstructions. We find the geomagnetic OSF and observed SSN agree very well after 1875, but do differ during the early part of the geomagnetic record, though still agree within the larger observational uncertainties.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Knowledge of space climate is a necessary first step to understanding past and future space weather (Temmer, Veronig, and Hanslmeier, 2003; Barnard et al., 2011; Lockwood et al., 2017; Riley and Love, 2017; Lockwood et al., 2018; Hayakawa et al., 2020a; Owens et al., 2021; Pevtsov et al., 2023). Historical space climate can be reconstructed from a range of disparate observational sources (e.g. Usoskin, 2023). In order of increasing time span these are:
-
Spacecraft data. High-resolution in-situ measurements of the near-Earth solar wind and heliospheric magnetic field have been made since 1964 (King and Papitashvili, 2005), though with patchy coverage prior to 1995 (Lockwood et al., 2019).
-
Neutron monitors. High-resolution measurements of the atmospheric neutron intensity have been made back to 1951 (Simpson, 2000). This provides information about the energy-integrated flux (above a threshold set by the instrument’s geomagnetic latitude) of galactic cosmic rays (GCRs) at the top of the atmosphere, which is modulated by the solar and geomagnetic magnetic fields (e.g. Potgieter, 2013).
-
Geomagnetic data. Ground-based magnetometer data measure the disturbance to the terrestrial system resulting from near-Earth solar-wind conditions (Feynman and Crooker, 1978). This can be used to reconstruct near-Earth magnetic-field intensity (\(B\)) and solar-wind speed (\(V\)) at annual time scale, back to around 1845 (Lockwood, Stamper, and Wild, 1999; Svalgaard and Cliver, 2010; Lockwood and Owens, 2011; Lockwood et al., 2022).
-
Sunspot observations. Sunspot number (SSN) is a proxy for photospheric magnetic flux intensity. Semi-homogeneous monthly SSN records can be reconstructed back to 1749 (e.g. Clette et al., 2023, and references therein) with a reasonable degree of certainty. Prior to that, data coverage is poorer (Hayakawa et al., 2022) and producing a consistent SSN record is difficult. However, incomplete records exist back to 1610 (e.g. Vokhmyanin, Arlt, and Zolotova, 2020; Arlt and Vaquero, 2020; Usoskin, 2023).
-
Cosmogenic isotopes. In addition to energetic neutrons, GCR interaction with the atmosphere produces radioactive isotopes that are not otherwise naturally occurring, such as 14C, which gets stored in tree rings that can be annually dated, and 10Be, which is deposited in ice sheets (e.g. Beer et al., 2013; Muscheler et al., 2016). Ice sheet dating for more recent times is possible by counting annual striations in ice colour, but for earlier times requires ice sheet modelling, constrained by layers of ash from identified volcanic eruptions. Recently, annual 14C has been used to derive sunspot number on a millennial timescale (Brehm et al., 2021; Usoskin et al., 2021).
In addition to these quasi-continuous time series, there also exist episodic measurements of the corona, in the form of solar eclipse observations. Total solar eclipses enable the solar coronal streamers to be observed with the naked eye, serving as independent spot references for estimates of the large-scale solar magnetic field. Solar coronal streamers typically concentrate near the solar equator around solar minimum and radiate all the directions near solar maximum (Loucif and Koutchmy, 1989; Riley et al., 2015). Dateable graphical records of eclipses are confirmed up to 1706, including during the Dalton Minimum and Maunder Minimum (Vaquero, 2003; Hayakawa et al., 2020b, 2021; Lockwood and Owens, 2021; Hayakawa, Sôma, and Daigo, 2022).
While information about space climate can be extracted from each of these data sources, they reflect very different facets of solar activity. Spacecraft and geomagnetic reconstructions pertain to local, near-Earth solar-wind conditions. Conversely, GCRs traverse the whole heliosphere to reach Earth and thus solar modulation of GCRs is a measure of the global heliosphere (Potgieter, 2013). Similarly, sunspot number provides information about the synoptic magnetic field in the solar photosphere, rather than the near-Earth conditions specifically. Eclipse observations provide mainly qualitative evidence about the state of the coronal magnetic field, although the width of the streamer belt can be estimated (Lockwood et al., 2022). In order to assess and compare these different data sources, it is necessary to find some common ground.
The open solar flux (OSF) is the total magnetic flux that is dragged out from the corona by the solar wind to form the heliospheric magnetic flux (HMF, Lockwood, 2013; Owens and Forsyth, 2013). Thus it is useful to define a theoretical ‘source surface’, typically around 2 to 3 solar radii, where the field is everywhere radial (Arge and Pizzo, 2000). This means the OSF can be simply defined as the integrated unsigned magnetic flux (with both toward and away polarities) threading the source surface. OSF is closely related to the solar modulation of GCRs (Herbst et al., 2010; Usoskin et al., 2002, 2021). It can also be computed from local (e.g. near-Earth) HMF observations (Owens et al., 2008), owing to the observed latitudinal invariance in the radial component of the HMF (Smith and Balogh, 1995; Lockwood et al., 2004). Finally, OSF can be modelled by a continuity equation, using sunspot number as the source of new OSF (Solanki, Schüssler, and Fligge, 2000; Krivova, Balmaceda, and Solanki, 2007; Owens and Lockwood, 2012; Krivova et al., 2021).
While OSF serves as a useful common factor relevant to the different space climate proxies, it is important to relate this variable to other commonly used measures of solar activity. Of particular interest is SSN, as it most readily enables solar-cycle information (e.g. cycle amplitude, cycle length, cycle start/end dates) to be estimated. Furthermore, it is SSN, not OSF, that is used as a general measure of solar activity levels, as well as for input to models of total solar irradiance (Lean, 2000; Krivova, Balmaceda, and Solanki, 2007; Solanki, Krivova, and Haigh, 2013), which in turn is used by climate models to incorporate the effect of solar forcing (e.g. Solanki, Krivova, and Haigh, 2013; Owens et al., 2017a).
Here we demonstrate a new method for reconstructing SSN purely from OSF, without the need to prescribe solar cycle start/end dates, which have been required by previous approaches (Owens and Lockwood, 2012). This can be particularly important during times of very low solar activity, when the simple correlation between OSF and SSN breaks down (Owens, Usoskin, and Lockwood, 2012). In Section 2 we describe the OSF data used in this study, while Section 3 outlines the sunspot data and the method for determining the solar cycle start/end dates. An existing model for reconstructing OSF from a known SSN record and known solar cycle start dates is demonstrated in Section 4. A new approach to the inverse problem – estimating SSN from OSF – is outlined in Section 5. To assess the performance of this new method, we outline benchmark methods, based on regression, in Section 6. The results are evaluated in Section 7. Implications of these results and future use of this methodology are discussed in Section 8.
2 OSF Data
In principle, OSF can be estimated from in-situ observations of the radial component of the HMF, \(B_{\mathrm{R}}\), made at a heliocentric distance, \(r\), by computing \(\mathrm{OSF} = 4 \pi r^{2} |B_{\mathrm{R}}|\). This assumes that local measurements of \(B_{\mathrm{R}}\) are representative of the whole sphere of radius \(r\). Longitudinal structure is accommodated by integrating over a solar rotation (≈ 27 days, for spacecraft in near-Earth space). Integration in latitude is achieved by assuming no variation in \(B_{R}\) with latitude, an assumption found to be valid to considerable accuracy by measurements from the Ulysses spacecraft (Smith and Balogh, 1995; Lockwood et al., 2004). However, calculating OSF in this way is further complicated by the existence of inverted HMF (Kahler, Crooker, and Gosling, 1998; Crooker et al., 2004); flux tubes which thread a sphere at the point of observation, typically \(r = 1\) AU, multiple times, but only thread the coronal source surface once. This produces an ‘excess flux’ in the heliosphere compared to the corona, which may partially explain the difference between in-situ and photospheric OSF estimates (Linker et al., 2017; Wallace et al., 2019). Excess flux is caused by any processes generating inverted HMF (Lockwood, Owens, and Macneil, 2019), including coronal interchange reconnection (Owens, Crooker, and Lockwood, 2013), large-scale solar-wind stream shears (Lockwood, Owens, and Rouillard, 2009), turbulence (Bruno and Carbone, 2005) and kinetic-scale ‘switchbacks’ (e.g. Mozer et al., 2021).
The most accurate OSF estimate is obtained by combining in-situ HMF and suprathermal electron observations of the field-aligned beam, or ‘strahl’ (Pilipp et al., 1987; Gosling et al., 1987), which allows inverted HMF to be directly detected and accounted for. The methodology for calculating OSF in this way is described in Owens et al. (2017b). We here use the 27-day averages of OSF provided as supplementary material to Frost et al. (2022). These span 1995 to 2021 and are here further averaged to annual resolution. This strahl-based OSF time series is shown in red in Figure 1 and referred to as OSFS.
To estimate OSF from in-situ observations prior to 1994, when strahl observations are not routinely available, observations of \(B_{\mathrm{R}}\) must be used in isolation. Frost et al. (2022) showed that the removal of inverted HMF can be approximated by taking 20-hour averages of \(B_{\mathrm{R}}\) before integrating over the source surface (i.e. before computing OSF = \(4 \pi r^{2} |B_{\mathrm{R}}|\)). This estimate is referred to as OSFBR. Using in-situ measurements of the near-Earth HMF \(B_{\mathrm{R}}\), as compiled by the OMNI database of near-Earth spacecraft observations (King and Papitashvili, 2005), OSFBR is shown as the black line in Figure 1. At the annual time scale, OSFS and OSFBR agree extremely well, with a linear correlation coefficient of \(r_{L} = 0.99\) and a mean absolute error (MAE) of \(0.175 \times 10 ^{14}\) Wb.
Geomagnetic reconstructions of near-Earth solar-wind conditions, discussed further below, do not directly reconstruct \(|B_{\mathrm{R}}|\). Instead, they provide the near-Earth \(B\) and \(V\) at annual resolution. It is important to test how well OSF can be reconstructed from such parameters. To do so, we follow a similar approach here as in Lockwood et al. (2022). The first step is to estimate annual \(|B_{R}|\) from the annual \(B\) and \(V\) values. If an ideal Parker spiral (Parker, 1958) is assumed, then:
where \(\gamma \) is the angle of the Parker spiral HMF to the radial direction, given by:
where \(\Omega \) is the mean angular rotation speed of the Sun, taken to be \(2.865 \times 10^{-6}\text{ rad}\mbox{ s}^{-1}\).
The blue line in Figure 1a shows OSFPS, the OSF estimate from annual OMNI \(B\) and \(V\) values, assuming an ideal Parker spiral to compute the equivalent \(|B_{\mathrm{R}}|\). While the relative variations are in good agreement, there is a systematic overestimate compared with both OSFS and OSFBR. This is expected, as inverted HMF means that for a given \(B\), the effective \(B_{\mathrm{R}}\) is lower than that for an ideal Parker spiral. As discussed above, where vector HMF data are available, this can be accounted for by taking 20-hour average of \(B_{\mathrm{R}}\) before computing the modulus, which results in about a 75% reduction from 1-minute \(B_{\mathrm{R}}\) (Frost et al., 2022). To account for this, we compute the total (or orthogonal) least-squares linear regression (Markovsky and Van Huffel, 2007) between OMNI OSFBR and OMNI OSFPS, which is given by:
where \({\mathrm{OSF_{PS*}}}\) is the best estimate of OSF from Parker spiral analysis of \(B\) and \(V\), and uncertainty ranges are the 95% confidence level determined by a bootstrap resampling process (Press et al., 1989). This is shown by the blue-shaded region in Figure 1b. As this is a simple regression, it is important not to apply it outside the valid parameter range. Fortunately, the deep 2009/2010 solar minimum and the high 1991 solar maximum mean that OSF conditions have been sampled over a wide range of conditions during the OMNI era, meaning valid parameter space is not a major concern for subsequent use with geomagnetic data.
Applying this linear correction to the OMNI \(B\) and \(V\) data produces OMNI \({\mathrm{OSF_{PS}*}}\), shown as the blue line in Figure 1c. Again, the 95% confidence interval is computed from a bootstrap method and is shown as a blue shaded area. Comparing \({\mathrm{OSF_{PS}*}}\) with \({\mathrm{OSF_{BR}}}\) gives \(r_{L} = 0.98\) and MAE = \(0.31 \times 10^{14}\) Wb.
Having established the procedure by processing the OMNI data, we can now apply these methods to the geomagnetic reconstructions of annual \(B\) and \(V\). Geomagnetic (GEO) reconstructions, along with their uncertainties, are provided by Lockwood et al. (2022). Figures 2a and b show the geomagnetic \(B\) and \(V\) reconstructions, respectively. The pink-shaded regions show the prediction intervals at the 95% confidence level computed from both the uncertainty in the regression between geomagnetic indices and solar-wind parameters, and the intrinsic uncertainty in the geomagnetic indices themselves. The latter results in a wider uncertainty range in the 19th century. Geomagnetic \(B\) and \(V\) are converted to \({\mathrm{OSF_{PS}}}\) using the Parker spiral equations 1 and 2, then corrected to \({\mathrm{OSF_{PS}*}}\) using the regression in Equation 3. Note that in Lockwood et al. (2022), geomagnetic OSF estimates were scaled to OMNI estimates of OSF made using a kinematic correction (Lockwood, Owens, and Rouillard, 2009), rather than to \({\mathrm{OSF_{S}}}\) and \({\mathrm{OSF_{BR}}}\). This produces slightly higher values than those shown here. This GEO \({\mathrm{OSF_{PS}*}}\) is shown in red in Figure 2c. There is very good agreement with OMNI \({\mathrm{OSF_{BR}}}\), with a linear correlation coefficient of 0.91. The GEO \({\mathrm{OSF_{PS}*}}\) series is used in the following sections to test the SSN reconstruction procedure.
3 Solar Cycle Start/End Dates
To both demonstrate the existing OSF forward model and validate the new SSN reconstruction method, we require accurate SSN and accurate solar cycle start/end dates. SSN is provided by the Sunspot Index and Long-term Solar Observations (SILSO) Version 2.0 (SILSO) sunspot number record, which provides monthly values back to 1749 (Clette et al., 2016). There are alternative sunspot series (e.g. Chatzistergos et al., 2017), which may differ in the cycle amplitude, particularly prior to 1830 (e.g. Hayakawa et al., 2022). However, the cycles themselves are defined robustly throughout.
Precise determination of the solar cycle start/end dates from SSN alone can be ambiguous, as it requires pinpointing the minimum in a long and shallow variation at solar minimum. However, the position of sunspots – specifically their heliographic latitude, \(\theta \) – can provide an additional constraint (Owens et al., 2011). We here update the solar cycle start/end dates using the average latitude of sunspots. Figure 3 shows sunspot area and latitude data from a range of observatories, compiled and inter-calibrated by Mandal et al. (2020). The familiar ‘butterfly diagram’ (Maunder, 1904) of sunspot occurrence with latitude and time is shown in Figure 3a. Integrating the total sunspot area over all latitudes gives the time series shown in Figure 3b, which follows a very similar pattern to the SILSO SSN record, shown in Figure 4b. Figure 3c shows \(\langle |\theta |\rangle \), the mean absolute latitude of all sunspots, weighted by area. At the start of each cycle, there is a sharp rise in \(\langle |\theta |\rangle \), as spots from the new cycle emerge at high latitude. We use the time when \(\langle |\theta |\rangle \) exceeds 20∘ to define the cycle start/end dates. These times are shown as vertical red dashed lines and listed in Table 1. It can be seen that these dates occur during the times of sunspot number minimum, as expected. By using the rise in \(\langle |\theta |\rangle \), rather than the appearance of the first spot group of the new cycle at high latitudes, the cycle start/end times are not as variable because they are not influenced by the formation of just one group.
Prior to 1874, inter-calibrated routine observations of sunspot latitude are not available (Mandal et al., 2020) and cycle start/end dates must be estimated from sunspot number alone. For these cycles, we define the solar cycle start/end dates as the time of the lowest 13-month smoothed SILSO SSN value. During extended minima, when this value remains at zero for multiple months (e.g. around 1810, in the middle of the Dalton minimum) we take the mid-point of the zero SSN interval. These dates are also shown in Table 1. The mean absolute error between the two methods is 0.3 years. Solar-cycle phase is then defined as 0 at the start of a cycle, increasing linearly to 1 at the end of a cycle. In order to compute the solar-cycle phase during the current solar cycle (25), we use a Solar Cycle 26 start year of 2031.0. This is 11 years from the Solar Cycle 25 start year, which agrees well with the observed \(\langle |\theta |\rangle \) progression to date.
4 Reconstructing OSF from SSN
Before outlining the new method for reconstructing SSN and solar-cycle phase from OSF, we first demonstrate the inverse problem: reconstructing the OSF using the known SSN record and associated solar-cycle phase as input. We follow the approach first established by Solanki, Schüssler, and Fligge (2000), wherein the rate of change of the OSF with time \(t\) can be modelled as a continuity equation:
where \(S\) is the OSF source term (i.e. rate of production of new OSF) and \(L\) is the OSF loss rate. \(S\) is expected to be closely related to the sunspot number (Solanki, Schüssler, and Fligge, 2000; Owens and Lockwood, 2012; Goelzer et al., 2013; Krivova et al., 2021). OSF increases when new magnetic flux emerges through the source surface. Thus a significant contribution to \(S\) is expected to be the magnetic flux carried out from the corona by coronal mass ejections over the solar cycle (Owens and Crooker, 2006).
A number of functional forms have been used to relate SSN and \(S\). Here we use the empirically optimised form and coefficients from Owens et al. (2016):
The exact details of this equation and coefficients do not greatly affect the resulting OSF reconstruction, as \(L\) is subsequently computed to provide the best agreement between OSF and \(S\). Specifically, Owens and Lockwood (2012) used Equation 4 with known OSF and SSN (and hence \(S\)) time series to compute the required \(L\) as a function of time. We here repeat this method with the newly updated data: The GEO \({\mathrm{OSF_{PS}*}}\) series and \(S\) computed using Equation 5 and the SILSO sunspot record. They are shown in Figures 4a and b, respectively. The resulting OSF loss term is shown in Figure 4c, where \(L\), the absolute OSF loss, has been converted to \(\chi \), the percentage loss. I.e. \(L = \chi \) OSF.
While peak SSN varies by around a factor 2.5 between cycles, and peak \(S\) by around a factor 2, \(\chi \) is remarkably consistent between cycles, both in amplitude and waveform. For example, Cycle 12 peaked around 1885 with a maximum SSN of about 100 and hence a peak OSF source rate of around \(1 \times 10 ^{15}\) Wb/yr. Conversely, Cycle 19 peaked around 1960 with SSN around 270, and hence a peak OSF source rate of around \(1.75 \times 10 ^{15}\) Wb/yr. These are approximately a factor two variations. In contrast, \(\chi \) peaks at around 150 %/yr for both Cycles 12 and 19, despite the large difference in the SSN peaks for these two cycles.
Owens and Lockwood (2012) exploited the cyclic nature of \(\chi \) to reconstruct OSF from knowledge of SSN and solar-cycle phase alone (i.e. without needing any further information to estimate \(\chi \)). Using the GEO \({\mathrm{OSF_{PS}*}}\) and SILSO SSN time series, the cycle-averaged profiles of \(\chi \) and SSN are shown in Figure 5. Panels a and b show the SSN variation with solar-cycle phase. There is large variability in SSN amplitude between individual cycles. Panels c and d show the same SSN variation normalised by the maximum 13-month smoothed SSN value in that cycle. This is equivalent to the Solar Activity Index used by Owens et al. (2022) and is tabulated in Table 2. From these panels it can be seen that the quasi-sinusoidal waveform of SSN is fairly consistent between cycles, despite the variation in amplitude (e.g. Pesnell, 2012). However, the rising phase is generally shorter for higher cycles (the so-called Waldmeier rule; Hathaway, 2010), which is not captured by an average variation over all cycles.. Panels e and f show the cycle-averaged profile of \(\chi \), the fractional OSF loss rate. It is also tabulated in Table 2. The waveform is different to that of SSN, being more saw-tooth in nature and closely mimicking the heliospheric current sheet tilt angle (Alanko-Huotari et al., 2007; Owens, Crooker, and Lockwood, 2011; Asvestari and Usoskin, 2016), as expected from OSF disconnection observations (Sheeley, Knudson, and Wang, 2001).
We now use SSN and solar-cycle phase (computed from the dates in Table 1) as input to the OSF model. \(S\) is determined from SSN using Equation 5 and \(\chi \) is computed from the solar-cycle phase using the cycle-averaged average profile, shown in Figure 5f and Table 2. These time series are used with a numerical method, with a 1-year time step and an initial value for the OSF. We use \(8\times 10^{14}\) Wb, though the choice of this value is lost within a few time steps and so does not greatly affect the resulting time series. At each time step, \(L\) is computed using \(\chi \) and the OSF at the previous time step. We then apply Equation 4 to compute the OSF. See Owens and Lockwood (2012) for more details.
The resulting OSF variation is shown in Figure 6. There is good agreement between the SSN-based model and observed OSF (GEO \({\mathrm{OSF_{PS}*}}\)), with a linear correlation coefficient of 0.84 (\(p < 0.00001\)), and a mean absolute error of \(0.8 \times 10^{14}\) Wb (i.e. around 10%). Note that the largest disagreement is prior to 1880: this could be from an overestimate of the sunspot numbers at this time, but more likely is an underestimate in the early geomagnetic observations. In fact, the reconstructions in Figure 12 of Lockwood et al. (2022) show that using using different combinations of geomagnetic indicies to reconstruct \(B\) and \(V\) can reduce this disagreement.
5 Reconstructing SSN from OSF
We now consider the inverse problem: reconstructing SSN using only the observed OSF. This is again achieved by assuming a purely cyclic variation in the fractional OSF loss rate, \(\chi \). The time dependence of Equation 4 and the lack of knowledge of solar-cycle phase – as the cycle start/end dates are not known a priori – means the model cannot be directly inverted. We will later demonstrate that assuming the time of minimum OSF is the cycle start date is only accurate to within around 0.5 years. Furthermore, this error is expected to increase significantly during times of very low OSF (Owens, Usoskin, and Lockwood, 2012). Instead, we take a Monte Carlo approach of producing random realisations of SSN times series and forward modelling the associated OSF for comparison with observed OSF. It is presumed that the best match to the observed OSF results from the SSN time series closest to the true value.
An example is shown in Figure 7, considering a 22-year window spanning 1920 to 1942. Figure 7a shows the result of selecting a random start date for the first solar cycle of 1912.8, and a random length and amplitude for this cycle of 11 years and 203, respectively. This is followed by a cycle of length 11.6 and amplitude 68, again selected randomly. Finally, this is followed by a third cycle of 9.4 years and amplitude 177. This reaches the end of the window of interest. (In some instances, a fourth cycle is required to span the 22-year window.) The solar-cycle length and amplitudes are used with the average SSN waveform (Figure 5d) to produce the synthetic SSN time series shown in red. This SSN time series is then used as input to the OSF model (Equation 4). The resulting OSF is shown in red in Figure 7b. There is little correspondence with the observed OSF (GEO OSF\(_{\rm PS}*\)), shown in black, as reflected by the high mean absolute error (MAE) of \(1.8\times 10^{14}\) Wb. Thus we would assume the synthetic SSN time series is not a good reproduction of the true SSN variation.
A second attempt at a random SSN time series is shown in Figure 7c. This produces an even worse match to the observed OSF, with \(\text{MAE} = 2.2\times 10^{14}\) Wb. It can be seen that this is largely due to setting the wrong solar cycle start date and hence a synthetic SSN time series out of phase with the true SSN (though the true SSN is obviously not available to the reconstruction process). Finally, the bottom panel shows a synthetic SSN time series which produces good match to the observed OSF, with \(\text{MAE} = 0.6\times 10^{14}\) Wb. Thus out of these three random SSN realisations, we would select the one in the bottom panels. As the SSN is known in this example, we can confirm that it is indeed the closest to the true value.
More specifically, the implementation of the process can be described as follows: We first split the observed OSF time series (in this case GEO OSF\(_{\rm PS}*\)) into overlapping 22-year windows, spaced at 1-year intervals. For each window we then:
-
a.
Draw a random solar-cycle phase for the start of the window from a uniform distribution between 0 and 1.
-
b.
Select a random solar-cycle length from a Gaussian distribution with mean of 10.5 years and a standard deviation of 2 years.
-
c.
Combine the starting solar-cycle phase with the cycle length to determine the start and end years for this cycle.
-
d.
Select additional random cycle lengths and append them to the previous cycles until the the whole 22-year window is filled.
-
e.
Compute solar-cycle phase over the whole 22-year window from the cycle start/end times.
-
f.
For each cycle, select random peak SSN amplitudes from a Gaussian distribution with a mean of 140 and a standard deviation of 70.
-
g.
Use the average sunspot-cycle shape (Table 2) with the amplitudes and solar-cycle phase to construct synthetic SSN time series spanning the 22-year window. Three examples are shown in red in the left-hand panels of Figure 7.
-
h.
Use the synthetic SSN time series to forward model the OSF using Equation 4 and \(\chi \) determined from Table 2. Three examples are shown in red in the right-hand panel of Figure 7.
-
i.
Compute the mean absolute error (MAE) between the model and observed OSF.
-
j.
Repeat steps (a) to (i) 10k times to create 10k OSF realisations from 10k random synthetic SSN time series.
-
k.
On the basis of the lowest MAE between synthetic and observed OSF, identify the ‘best’ OSF realisation from the 10k available, and store the corresponding synthetic SSN time series for this 22-year window.
Once the entire OSF dataset has been modelled in this way, the over-lapping nature of the windows mean that each individual year will have 22 independent ‘best’ estimates of OSF and SSN. There are also 22 sets of solar cycle start dates. We produce consensus values of OSF and SSN for each year by taking the average values across all 22 windows, weighted according to their MAE (specifically, by \(\exp [{-\mathrm{MAE}}]\)). Consensus solar cycle start dates are produced using a method explained in Section 7.
In Section 7 we apply this method to the GEO \({\mathrm{OSF_{PS}*}}\) series, as well as independently to the upper and lower 95% confidence bounds of the GEO \({\mathrm{OSF_{PS}*}}\) series. This gives the uncertainty in the reconstructed SSN record that is a result of the uncertainty in the input OSF data.
6 Benchmark Reconstruction Methods
In order to assess the performance of the forward-modelling approach outlined in Section 5, it is useful to compare it with existing benchmark models for estimating SSN from OSF. Svalgaard and Cliver (2005) and Wang, Lean, and Sheeley (2005) proposed a linear relation between OSF and the square root of SSN. As we require the inverse relation, we consider the linear regression between SSN and OSF2. An orthogonal least-squares fit using SILSO and GEO OSF\(_{\rm PS}*\) gives:
where uncertainties are the 95% confidence interval from a bootstrapping analysis. The best fit is shown as the blue line in Figure 8a.
Figure 9a shows the SSN time series which results from this OSF2 regression applied to GEO \({\mathrm{OSF_{PS}*}}\). The general trends in cycle-to-cycle amplitude are present, but there are obvious issues, in particular the solar minima often produce unphysical negative SSN values (though most of these are consistent with SSN = 0 to within uncertainty). Figure 8d shows the same data as a scatter plot. The linear correlation coefficient is 0.78, with a mean absolute error of 38.4. More advanced forms of regression, such as Gaussian processes, could be used to further optimise the agreement, but a fundamental issue is the underlying assumption that a single value of OSF maps to a single value of SSN. This assumption means that there is no allowance for time history, which at least partially explains the large range of OSF for low SSN values; zero SSN after a period of previously high SSN does not result in as low OSF as an extended period of zero SSN. And this degeneracy is expected to increase further at even lower solar activity considered here: During grand solar minima, such as the Maunder minimum, it has been proposed that OSF peaks at solar minimum, rather than solar maximum (Owens, Usoskin, and Lockwood, 2012). Hence there are considerable limitations to any simple regression-based method.
A more nuanced approach to the regression between OSF and SSN is to independently account for the short (year-to-year) and long (cycle-to-cycle) timescale variations in OSF (Usoskin and Kovaltsov, 2004). Considering first the cycle-to-cycle variations, Figure 8b shows the scatter plot of 11-year running means of SSN and OSF, denoted \(\langle \mathrm{SSN}\rangle _{11}\) and \(\langle \mathrm{OSF}\rangle _{11}\), respectively. The best linear regression from an orthogonal least-squares fit is given by:
where uncertainties are again the 95% confidence interval from a bootstrapping analysis.
We then fit the shorter-term variations by considering the anomalies from the 11-year mean values. The scatter plot of \(\mathrm{SSN} - \langle \mathrm{SSN}\rangle _{11}\) versus \(\mathrm{OSF} - \langle \mathrm{OSF}\rangle _{11}\) is shown in Figure 8c. The best fit is given by:
Combining Equations 7 and 8 with the observed GEO \({\mathrm{OSF_{PS}*}}\) results in the SSN reconstruction shown in Figure 8e. This is referred to as the ‘1 + 11y regression’. The linear correlation coefficient is 0.80, with a mean absolute error of 30.3. Using a Meng’s z-test (Meng, Rosenthal, and Rubin, 1992), this correlation is not found to be significantly different from the OSF2 regression in Figure 8d. Figure 9b shows that the OSF2 and 1 + 11y regression time series are similar. The 1 + 11y regression still produces unphysical negative sunspot numbers, though these are again mostly consistent with zero to within regression and OSF uncertainties, as shown by Figure 8e.
To obtain benchmark estimates of the solar cycle start/end dates, we assume that the solar cycle start coincides with minima in annual OSF. These dates are given in Table 1. From these we compute \(\Delta t\), the absolute timing error compared with the observed cycle start/end dates derived from sunspot number latitude. The cumulative distribution function of \(\Delta t\) is shown in Figure 10. The median \(\Delta t\) is 0.5 years, with an inter-quartile range (IQR) spanning 0.2 to 0.60 years.
7 Results
We now apply the forward-modelling methodology to the observed OSF (GEO OSF\(_{\rm PS}*\)) series. The model OSF is shown in Figure 9c and the associated model SSN in Figure 9b. The fact that the model OSF agrees well with the observed OSF is to be expected, as the model specifically selects the random SSN realisations that produce the best match to the observed OSF. Indeed, given enough free parameters in the construction of the underlying SSN series, a perfect match with the observed OSF could likely be produced. Times when the model is unable to reproduce the observed OSF record could result from one or more of three possibilities: 1. the assumption of the cycle-average SSN waveform is insufficient, such as the single year spikes in OSF around cycle maxima; 2. the assumed model loss term does not permit the required rate of change of OSF; or 3. the geomagnetic OSF record is inaccurate because of magnetometer calibration errors.
The test of the reconstruction methodology is the obtained SSN time series, shown in 9b. By eye, the model reconstructs the timing of the SSN maxima and minima very well. As shown in Figure 8f, we find a linear correlation coefficient of 0.95 and a MAE of 15.9, greatly improved over both regression methods. A Meng’s z-test finds the increase in correlation compared with regression to be significant at the 99.99% level.
For the forward modelling SSN, the relative amplitude of the cycle maxima is also well matched, and the minima are far improved over the regression methods. There is a slight tendency for the reconstruction to underestimate the cycle peaks. This is likely due to the average SSN waveform used by the model, which produces a smooth, quasi-sinusoidal SSN variation. The largest differences between the model and observed SSN are the early cycle peaks between 1845 and 1880, and – to a lesser extent – the later peaks in 1970 and 1980. Here the model values are between a factor 1.2 to 1.7 lower than observed. For the early cycles, the uncertainty in the geomagnetic OSF record is expected to be a significant contribution to the disagreement. In other words, our results suggest the GEO OSF is too low at this time, assuming the SSN record to be correct, which may be also not exactly true (see, e.g. discussion in Usoskin, 2023). However, for the 1970 peak, OSF uncertainty cannot completely explain the SSN underestimate. This is discussed further in Section 8. However, in general, the reconstructed and observed SSN agree within observational uncertainties.
As the SSN reconstruction works on a 22-year window which slides forward 1 year at a time, we obtain 22 estimates of each solar cycle start/end date. In order to produce consensus estimates of the solar cycle start dates, a Gaussian kernel density with a width of 0.05 years is applied to the solar-cycle start date in each 22-year window and then summed across all windows. A range of kernel widths were tried and 0.05 years was found to be broad enough to combine information from multiple windows, but narrow enough to resolve individual cycles. The resulting kernel-density time series is shown in red in Figure 9. It can be considered to be the model (unnormalised) probability of a cycle starting on a given date. Note that the units are arbitrary. Taking maxima in the kernel density to be the model cycle start/end dates results in the values listed in Table 1. The cumulative distribution of the solar cycle timing error, \(\Delta t\), is shown as the red line in Figure 10. The \(\Delta t\) median value is 0.29 years, with an IQR of spanning 0.14 to 0.34 years (cf. the regression approach, with a median \(\Delta t\) of 0.5 years, and an IQR spanning 0.2 to 0.60 years). Note that the forward-model approach additionally gives estimates of the Cycle 9 and 25 start dates, even though the input OSF data did not contain these minima. This is because the solar-cycle phase is fitted to the available data, enabling extrapolation beyond the period of observations. In principle, the width of each kernel-density peak also provides information about uncertainty in the reconstructed cycle start date, but that information is not exploited here.
8 Discussion
This study has described a method for reconstructing sunspot number (SSN) from open solar flux (OSF) estimates by the use of a forward model. The method was tested over the period 1850 to present, when the OSF is known from a newly calibrated geomagnetic OSF estimate and the reconstructed SSN can be verified against the observed record. The primary motivation for developing this method is to reconstruct SSN prior to 1610 using cosmogenic isotope estimates of OSF.
Compared with simple regressions between SSN and OSF, the forward modelling approach significantly improves the reconstruction of SSN. The linear correlation coefficient for forward modelling of SSN is 0.95, compared to 0.78 or 0.80 for regression methods. And the mean absolute error in SSN is more than halved (15.9 versus 38.4 or 40.7). The estimates of the solar cycle start/end dates are also substantially improved, with the median timing error being reduced from 0.5 to 0.29 years. As the forward modelling approach has a theoretical underpinning, it is also expected to be more reliable during periods of extremely low SSN and OSF – grand solar minima – when a simple regression is expected to further break down (Owens, Usoskin, and Lockwood, 2012). Though this cannot be directly tested with the available data.
Future work will apply this new method to OSF inferred from cosmogenic isotope records, such as annual 14C and 10Be (McCracken, 2007; Muscheler et al., 2007; McCracken and Beer, 2015; Muscheler et al., 2016; Usoskin et al., 2021; Nguyen et al., 2022). It may also be possible to adapt the forward-modelling approach to directly link SSN to the cosmogenic isotope production rate, without the intermediate step of computing OSF. By inferring the solar-cycle phase as part of the reconstruction process, this would enable known heliospheric current sheet inclination effects (Alanko-Huotari et al., 2007) to be accounted for. Additionally, by working backwards from the present era, the effect of odd/even cycle drift effects can also be incorporated, which is presently difficult to achieve, for example the solar-cycle parity is unknown prior to the Dalton minimum (Usoskin, Mursula, and Kovaltsov, 2001; Owens et al., 2015). These advances in interpretation of cosmogenic isotope abundance data may enable reconstruction of SSN from proxy data that is sufficiently accurate to provide a link between observed SSN records across gaps in observational coverage, such as the early 18th century (Hayakawa et al., 2022).
However, as shown in this study, even in the modern era, when the SSN is known from direct observations, the method can still provide valuable insight. It provides an assessment of when the independent OSF and SSN records agree within the model and observational uncertainties. For example, Solar Cycle 20 – spanning approximately 1965 to 1975 – has a higher SSN than the low-and-flat OSF variation would suggest. While violation of the model assumptions at this time cannot be ruled out, they do appear to hold well for other cycles. The observational uncertainties in geomagnetic solar-wind reconstructions and sunspot-number records are also very low at this time, and thus unlikely to be the entire source of the disagreement. We note that the geomagnetic (and in-situ) estimates of OSF are based purely on local, near-Earth solar wind conditions. Thus while the uncertainty in measuring annual \(B\), \(V\) or \(B_{\mathrm{R}}\) is low, the extrapolation to global OSF may be less certain. Unfortunately, there were no 1 AU out-of-ecliptic observations available at the time to test this hypothesis (Owens et al., 2008). However, careful analysis of the Pioneer 6, 7, and 11 data in the more distant heliosphere may provide some information. There is also a disagreement between model and observed SSN in the very early geomagnetic OSF record. During Cycle 9 – spanning approximately 1844 to 1856 – the observed SSN is higher than the geomagnetic OSF suggests. However, the uncertainty in the geomagnetic \(V\) and \(B\) records, from which OSF is estimated, is much larger than the later record. Thus observational error is the most likely explanation for this discrepancy at this time.
The method used here can also bring the cosmogenic isotope data to bear on debates about calibration differences in the various data series. For example, Figure 9 shows that the forward modelling is giving SSN peaks that are lower than the SILSO values in Solar Cycles 9 – 11, with the largest difference in Cycle 9. This is most likely due to the calibration of the magnetometer data which underpins the geomagnetic OSF data, but could also reflect in the increasing and accumulating uncertainties in the calibration of the SSN record as we go back in time (e.g. Bhattacharya et al., 2023). Version 3 of the SILSO SSN record is currently under construction and may shed some light on this journal. The forward modelling approach using cosmogenic isotope data will also allow us to test for the origin of the early differences.
Data Availability
OMNI data are available from https://omniweb.gsfc81.nasa.gov/. Sunspot data are provided by the Royal Observatory of Belgium SILSO and available from www.sidc.be/silso/DATA/SN_m_tot_V2.0.csv. Sunspot area and latitude data are from http://www2.mps.mpg.de/projects/sun-climate/data/indivi_group_area_1874_2023_Mandal.txt. All analysis and visualisation code is packaged with all required data here: https://github.com/University-of-Reading-Space-Science/SSNfromOSF. The Github repository contains the raw and corrected geomagnetic OSF estimate.
References
Alanko-Huotari, K., Usoskin, I.G., Mursula, K., Kovaltsov, G.A.: 2007, Cyclic variations of the heliospheric tilt angle and cosmic ray modulation. Adv. Space Res. 40, 1064. DOI.
Arge, C.N., Pizzo, V.J.: 2000, Improvement in the prediction of solar wind conditions using near-real time solar magnetic field updates. J. Geophys. Res. 105, 10465.
Arlt, R., Vaquero, J.M.: 2020, Historical sunspot records. Living Rev. Solar Phys. 17(1), 1. DOI.
Asvestari, E., Usoskin, I.G.: 2016, An empirical model of heliospheric cosmic ray modulation on long-term time scale. J. Space Weather Space Clim. 6, A15. DOI.
Barnard, L., Lockwood, M., Hapgood, M.A., Owens, M.J., Davis, C.J., Steinhilber, F.: 2011, Predicting space climate change. Geophys. Res. Lett. 381. DOI.
Beer, J., McCracken, K.G., Abreu, J., Heikkilä, U., Steinhilber, F.: 2013, Cosmogenic radionuclides as an extension of the neutron monitor era into the past: potential and limitations. Space Sci. Rev. 176, 89. DOI.
Bhattacharya, S., Lefèvre, L., Hayakawa, H., Jansen, M., Clette, F.: 2023, Scale transfer in 1849: Heinrich Schwabe to Rudolf Wolf. Solar Phys. 298(1), 12. DOI.
Brehm, N., Bayliss, A., Christl, M., Synal, H.-A., Adolphi, F., Beer, J., Kromer, B., Muscheler, R., Solanki, S.K., Usoskin, I., Bleicher, N., Bollhalder, S., Tyers, C., Wacker, L.: 2021, Eleven-year solar cycles over the last millennium revealed by radiocarbon in tree rings. Nat. Geosci. 14(1), 10. DOI.
Bruno, R., Carbone, V.: 2005, The solar wind as a turbulence laboratory. Living Rev. Solar Phys. 2, 4. DOI.
Chatzistergos, T., Usoskin, I.G., Kovaltsov, G.A., Krivova, N.A., Solanki, S.K.: 2017, New reconstruction of the sunspot group numbers since 1739 using direct calibration and “backbone” methods. Astron. Astrophys. 602, A69. DOI.
Clette, F., Lefèvre, L., Cagnotti, M., Cortesi, S., Bulling, A.: 2016, The revised Brussels–Locarno sunspot number (1981 – 2015). Solar Phys. 291(9), 2733. DOI.
Clette, F., Lefèvre, L., Chatzistergos, T., Hayakawa, H., Carrasco, V.M.S., Arlt, R., Cliver, E.W., Dudok de Wit, T., Friedli, T.K., Karachik, N., Kopp, G., Lockwood, M., Mathieu, S., Muñoz-Jaramillo, A., Owens, M., Pesnell, D., Pevtsov, A., Svalgaard, L., Usoskin, I.G., van Driel-Gesztelyi, L., Vaquero, J.M.: 2023, Recalibration of the sunspot-number: status report. Solar Phys. 298(3), 44. DOI.
Crooker, N.U., Kahler, S.W., Larson, D.E., Lin, R.P.: 2004, Large-scale magnetic field inversions at sector boundaries. J. Geophys. Res. 109. DOI.
Feynman, J., Crooker, N.U.: 1978, The solar wind at the turn of the century. Nature 275, 626. DOI.
Frost, A.M., Owens, M., Macneil, A., Lockwood, M.: 2022, Estimating the open solar flux from in-situ measurements. Solar Phys. 297(7), 82. DOI.
Goelzer, M.L., Smith, C.W., Schwadron, N.A., McCracken, K.G.: 2013, An analysis of heliospheric magnetic field flux based on sunspot number from 1749 to today and prediction for the coming solar minimum. J. Geophys. Res. 118(12), 7525. DOI.
Gosling, J.T., Baker, D.N., Bame, S.J., Feldman, W.C., Zwickl, R.D.: 1987, Bidirectional solar wind electron heat flux events. J. Geophys. Res. 92, 8519. DOI.
Hathaway, D.H.: 2010, The solar cycle. Living Rev. Solar Phys. 7, 1. DOI.
Hayakawa, H., Sôma, M., Daigo, R.: 2022, Analyses of historical solar eclipse records in Hokkaido Island in the 18 – 19th centuries. Publ. Astron. Soc. Japan 74(6), 1275. DOI.
Hayakawa, H., Owens, M.J., Lockwood, M., Sôma, M.: 2020b, The solar corona during the total eclipse on 1806 June 16: graphical evidence of the coronal structure during the dalton minimum. Astrophys. J. 900(2), 114. DOI.
Hayakawa, H., Ribeiro, P., Vaquero, J.M., Gallego, M.C., Knipp, D.J., Mekhaldi, F., Bhaskar, A., Oliveira, D.M., Notsu, Y., Carrasco, V.M.S., Caccavari, A., Veenadhari, B., Mukherjee, S., Ebihara, Y.: 2020a, The extreme space weather event in 1903 October/November: an outburst from the quiet sun. Astrophys. J. Lett. 897(1), L10. DOI.
Hayakawa, H., Lockwood, M., Owens, M.J., Sôma, M., Besser, B.P., van Driel-Gesztelyi, L.: 2021, Graphical evidence for the solar coronal structure during the Maunder minimum: comparative study of the total eclipse drawings in 1706 and 1715. J. Space Weather Space Clim. 11, 1. DOI.
Hayakawa, H., Hattori, K., Sôma, M., Iju, T., Besser, B.P., Kosaka, S.: 2022, An overview of sunspot observations in 1727 – 1748. Astrophys. J. 941(2), 151. DOI.
Herbst, K., Kopp, A., Heber, B., Steinhilber, F., Fichtner, H., Scherer, K., Matthiä, D.: 2010, On the importance of the local interstellar spectrum for the solar modulation parameter. J. Geophys. Res. 115(D1). DOI.
Kahler, S., Crooker, N.U., Gosling, J.T.: 1998, Properties of interplanetary magnetic sector boundaries based on electron heat-flux flow directions. J. Geophys. Res. 103, 20603. DOI.
King, J.H., Papitashvili, N.E.: 2005, Solar wind spatial scales in and comparisons of hourly wind and ACE plasma and magnetic field data. J. Geophys. Res. 110. DOI.
Krivova, N.A., Balmaceda, L., Solanki, S.K.: 2007, Reconstruction of solar total irradiance since 1700 from the surface magnetic flux. Astron. Astrophys. 467, 335. DOI.
Krivova, N.A., Solanki, S.K., Hofer, B., Wu, C.-J., Usoskin, I.G., Cameron, R.: 2021, Modelling the evolution of the Sun’s open and total magnetic flux. Astron. Astrophys. 650, A70. DOI.
Lean, J.: 2000, Evolution of the Sun’s spectral irradiance since the Maunder minimum. Geophys. Res. Lett. 27(16), 2425. DOI.
Linker, J.A., Caplan, R.M., Downs, C., Riley, P., Mikic, Z., Lionello, R., Henney, C.J., Arge, C.N., Liu, Y., Derosa, M.L., Yeates, A., Owens, M.J.: 2017, The open flux problem. Astrophys. J. 848(1). DOI.
Lockwood, M.: 2013, Reconstruction and prediction of variations in the open solar magnetic flux and interplanetary conditions. Living Rev. Solar Phys. 10(1), 4. DOI.
Lockwood, M., Owens, M.J.: 2011, Centennial changes in the heliospheric magnetic field and open solar flux: the consensus view from geomagnetic data and cosmogenic isotopes and its implications. J. Geophys. Res. 116(A4), A04109. DOI.
Lockwood, M., Owens, M.: 2021, Cosmic meteorology. Astron. Geophys. 62(3), 3.12. DOI.
Lockwood, M., Owens, M.J., Macneil, A.: 2019, On the origin of ortho-gardenhose heliospheric flux. Solar Phys. 294(6), 85. DOI.
Lockwood, M., Owens, M.J., Rouillard, A.P.: 2009, Excess open solar magnetic flux from satellite data: 2. A survey of kinematic effects. J. Geophys. Res. 114(A11), 1. DOI.
Lockwood, M., Stamper, R., Wild, M.N.: 1999, A doubling of the Sun’s coronal magnetic field during the past 100 years. Nature 399, 437. DOI.
Lockwood, M., Forsyth, R.J., Balogh, A., McComas, D.J.: 2004, Open solar flux estimates from near-Earth measurements of the interplanetary magnetic field: comparison of the first two perihelion passes of the Ulysses spacecraft. Ann. Geophys. 22, 1395.
Lockwood, M., Owens, M.J., Barnard, L.A., Scott, C.J., Watt, C.E.: 2017, Space climate and space weather over the past 400 years: 1. The power input to the magnetosphere. J. Space Weather Space Clim. 7, A25. DOI.
Lockwood, M., Owens, M.J., Barnard, L.A., Scott, C.J., Watt, C.E., Bentley, S.: 2018, Space climate and space weather over the past 400 years: 2. Proxy indicators of geomagnetic storm and substorm occurrence. J. Space Weather Space Clim. 8, A12. DOI.
Lockwood, M., Bentley, S.N., Owens, M.J., Barnard, L.A., Scott, C.J., Watt, C.E., Allanson, O., Freeman, M.P.: 2019, The development of a space climatology: 3. Models of the evolution of distributions of space weather variables with timescale. Space Weather 17(1), 180. DOI.
Lockwood, M., Owens, M.J., Barnard, L.A., Scott, C.J., Frost, A.M., Yu, B., Chi, Y.: 2022, Application of historic datasets to understanding open solar flux and the 20th-century grand solar maximum. 1. Geomagnetic, ionospheric, and sunspot observations. Front. Astron. Space Sci. 9.
Loucif, M.L., Koutchmy, S.: 1989, Solar cycle variations of coronal structures. Astron. Astrophys. Suppl. Ser. 77, 45.
Mandal, S., Krivova, N.A., Solanki, S.K., Sinha, N., Banerjee, D.: 2020, Sunspot area catalog revisited: daily cross-calibrated areas since 1874. Astron. Astrophys. 640, A78. DOI.
Markovsky, I., Van Huffel, S.: 2007, Overview of total least-squares methods. Signal Process. 87(10), 2283. DOI.
Maunder, E.W.: 1904, Note on the distribution of sun-spots in heliographic latitude, 1874 to 1902. Mon. Not. Roy. Astron. Soc. 64(8), 747. DOI.
McCracken, K.G.: 2007, Heliomagnetic field near Earth, 1428 – 2005. J. Geophys. Res. 112, 9106. DOI.
McCracken, K.G., Beer, J.: 2015, The annual cosmic-radiation intensities 1391 – 2014; the annual heliospheric magnetic field strengths 1391 – 1983; and identification of solar cosmic ray events in the cosmogenic record 1800 – 1983. Solar Phys. 290(10), 3051. DOI.
Meng, X., Rosenthal, R., Rubin, D.B.: 1992, Comparing correlated correlation coefficients. Psychol. Bull. 111(1), 172. DOI.
Mozer, F.S., Bale, S.D., Bonnell, J.W., Drake, J.F., Hanson, E.L.M., Mozer, M.C.: 2021, On the origin of switchbacks observed in the solar wind. Astrophys. J. 919(1), 60. DOI.
Muscheler, R., Joos, F., Beer, J., Müller, S.A., Vonmoos, M., Snowball, I.: 2007, Solar activity during the last 1000yr inferred from radionuclide records. Quat. Sci. Rev. 26(1), 82. DOI.
Muscheler, R., Adolphi, F., Herbst, K., Nilsson, A.: 2016, The revised sunspot record in comparison to cosmogenic radionuclide-based solar activity reconstructions. Solar Phys. 291(9), 3025. DOI.
Nguyen, L., Suttie, N., Nilsson, A., Muscheler, R.: 2022, A novel Bayesian approach for disentangling solar and geomagnetic field influences on the radionuclide production rates. Earth Planets Space 74(1), 130. DOI.
Owens, M.J., Crooker, N.U.: 2006, Coronal mass ejections and magnetic flux buildup in the heliosphere. J. Geophys. Res. 111(A10), A10104. DOI.
Owens, M.J., Crooker, N.U., Lockwood, M.: 2011, How is open solar magnetic flux lost over the solar cycle? J. Geophys. Res. 116(A4), A04111. DOI.
Owens, M.J., Crooker, N.U., Lockwood, M.: 2013, Solar origin of heliospheric magnetic field inversions: evidence for coronal loop opening within pseudostreamers. J. Geophys. Res. 118, 1868. DOI.
Owens, M.J., Forsyth, R.J.: 2013, The heliospheric magnetic field. Living Rev. Solar Phys. 10, 5. DOI.
Owens, M.J., Lockwood, M.: 2012, Cyclic loss of open solar flux since 1868: the link to heliospheric current sheet tilt and implications for the Maunder minimum. J. Geophys. Res. 117(A4), A04102. DOI.
Owens, M.J., Usoskin, I., Lockwood, M.: 2012, Heliospheric modulation of galactic cosmic rays during grand solar minima: past and future variations. Geophys. Res. Lett. 39(19), n/a. DOI.
Owens, M.J., Arge, C.N., Crooker, N.U., Schwadron, N.A., Horbury, T.S.: 2008, Estimating total heliospheric magnetic flux from single-point in situ measurements. J. Geophys. Res. 113(A12), A12103. DOI.
Owens, M.J., Lockwood, M., Barnard, L., Davis, C.J.: 2011, Solar cycle 24: implications for energetic particles and long-term space climate change. Geophys. Res. Lett. 38(19), 1. DOI.
Owens, M.J., McCracken, K.G., Lockwood, M., Barnard, L.: 2015, The heliospheric Hale cycle over the last 300 years and its implications for a “lost” late 18th century solar cycle. J. Space Weather Space Clim. 5, A30. DOI.
Owens, M.J., Cliver, E.W., McCracken, K.G., Beer, J., Balogh, A., Barnard, L., Lockwood, M., Rouillard, A.P., Wang, Y.-M., Passos, S., Riley, P., Usoskin, I.: 2016, Near-Earth heliospheric magnetic field intensity since 1750. Part 1: sunspot and geomagnetic reconstructions. J. Geophys. Res. 121, 6048. DOI.
Owens, M.J., Lockwood, M., Hawkins, E., Usoskin, I., Jones, G.S., Barnard, L., Schurer, A., Fasullo, J.: 2017a, The Maunder minimum and the little ice age: an update from recent reconstructions and climate simulations. J. Space Weather Space Clim. 7, A33. DOI.
Owens, M.J., Lockwood, M., Riley, P., Linker, J.: 2017b, Sunward strahl: a method to unambiguously determine open solar flux from in situ spacecraft measurements using suprathermal electron data. J. Geophys. Res. 122(11), 910. DOI.
Owens, M.J., Lockwood, M., Barnard, L.A., Scott, C.J., Haines, C., Macneil, A.: 2021, Extreme space-weather events and the solar cycle. Solar Phys. 296(5), 82. DOI.
Owens, M.J., Chakraborty, N., Turner, H., Lang, M., Riley, P., Lockwood, M., Barnard, L.A., Chi, Y.: 2022, Rate of change of large-scale solar-wind structure. Solar Phys. 297(7), 83. DOI.
Parker, E.N.: 1958, Dynamics of the interplanetary gas and magnetic fields. Astrophys. J. 128, 664. DOI.
Pesnell, W.D.: 2012, Solar cycle predictions (invited review). Solar Phys. 281(1), 507. DOI.
Pevtsov, A., Nandy, D., Usoskin, I.G., Pevtsov, A.A., Corti, C., Lefevre, L., Owens, M.J., Li, G., Krivova, N.A., Saha, C., Perri, B., Brun, S.A., Strugarek, A., Dayeh, M., Nagovitsyn, Y., Erdelyi, R.: 2023, Long-term Solar Variability: ISWAT S1 Cluster Review for COSPAR Space Weather Roadmap. Adv. Space Res. in press.
Pilipp, W.G., Muehlhaeuser, K.-H., Miggenrieder, H., Rosenbauer, H., Schwenn, R.: 1987, Variations of electron distribution functions in the solar wind. J. Geophys. Res. 92, 1103. DOI.
Potgieter, M.: 2013, Solar modulation of cosmic rays. Living Rev. Solar Phys. 10, 3.
Press, W.H., Flannery, B.P., Teukolsky, S.A., Vetterling, W.T.: 1989, Numerical Recipes, Cambridge University Press, New York.
Riley, P., Love, J.J.: 2017, Extreme geomagnetic storms: probabilistic forecasts and their uncertainties. Space Weather 15(1), 53. DOI.
Riley, P., Lionello, R., Linker, J.A., Cliver, E., Balogh, A., Beer, J., Charbonneau, P., Crooker, N., DeRosa, M., Lockwood, M., Owens, M., McCracken, K., Usoskin, I., Koutchmy, S.: 2015, Inferring the structure of the solar corona and inner heliosphere during the Maunder minimum using global thermodynamic magnetohydrodynamic simulations. Astrophys. J. 802(2), 105. DOI.
Sheeley, N.R. Jr, Knudson, T.N., Wang, Y.-M.: 2001, Coronal inflows and the Sun’s nonaxisymmetric open flux. Astrophys. J. 546(2), L131.
Simpson, J.A.: 2000, The cosmic ray nucleonic component: the invention and scientific uses of the neutron monitor-(keynote lecture). Space Sci. Rev. 93, 11. DOI.
Smith, E.J., Balogh, A.: 1995, Ulysses observations of the radial magnetic field. Geophys. Res. Lett. 22, 3317. DOI.
Solanki, S.K., Krivova, N.A., Haigh, J.D.: 2013, Solar irradiance variability and climate. Annu. Rev. Astron. Astrophys. 51(1), 311.
Solanki, S.K., Schüssler, M., Fligge, M.: 2000, Evolution of the Sun’s large-scale magnetic field since the Maunder minimum. Nature 408, 445. DOI.
Svalgaard, L., Cliver, E.W.: 2005, The IDV index: its derivation and use in inferring long-term variations of the interplanetary magnetic field strength. J. Geophys. Res. 110, 12103. DOI.
Svalgaard, L., Cliver, E.W.: 2010, Heliospheric magnetic field 1835 – 2009. J. Geophys. Res. 115, 9111. DOI.
Temmer, M., Veronig, A., Hanslmeier, A.: 2003, Does solar flare activity lag behind sunspot activity? Solar Phys. 215(1), 111. DOI.
Usoskin, I.G.: 2023, A history of solar activity over millennia. Living Rev. Solar Phys. 20(1), 2. DOI.
Usoskin, I.G., Kovaltsov, G.A.: 2004, Long-term solar activity: direct and indirect study. Solar Phys. 224(1), 37. DOI.
Usoskin, I.G., Mursula, K., Kovaltsov, G.A.: 2001, Was one sunspot cycle lost in late XVIII century? Astron. Astrophys. 370(2), L31.
Usoskin, I.G., Mursula, K., Solanki, S.K., Schüssler, M., Kovaltsov, G.A.: 2002, A physical reconstruction of cosmic ray intensity since 1610. J. Geophys. Res. 107(A11), 13. DOI.
Usoskin, I.G., Solanki, S.K., Krivova, N.A., Hofer, B., Kovaltsov, G.A., Wacker, L., Brehm, N., Kromer, B.: 2021, Solar cyclic activity over the last millennium reconstructed from annual 14C data. Astron. Astrophys. 649, A141. DOI.
Vaquero, J.M.: 2003, The solar corona in the eclipse of 24 June 1778. Solar Phys. 216(1), 41. DOI.
Vokhmyanin, M., Arlt, R., Zolotova, N.: 2020, Sunspot positions and areas from observations by Thomas Harriot. Solar Phys. 295(3), 39. DOI.
Wallace, S., Arge, C.N., Pattichis, M., Hock-Mysliwiec, R.A., Henney, C.J.: 2019, Estimating total open heliospheric magnetic flux. Solar Phys. 294(2), 19. DOI.
Wang, Y.-M., Lean, J.L., Sheeley, N.R. Jr: 2005, Modeling the sun’s magnetic field and irradiance since 1713. Astrophys. J. 625(1), 522.
Acknowledgments
We have benefited from sunspot data provided by Royal Observatory of Belgium SILSO, sunspot latitude and area data collated by Mandal et al. and OMNI data provided by NASA/SPDF. We benefited from useful discussions as part of the International Space Science Institute (ISSI, Bern) team “Magnetic open flux and solar wind structuring in interplanetary space” (2019 – 2021), led by Manuela Temmer. IU acknowledges the ISSI visiting fellowship.
Funding
This work was part-funded by Science and Technology Facilities Council (STFC) grant number ST/V000497/1. I.U. acknowledges the Academy of Finland (projects ESPERA No. 321882). H.H. thanks JSPS Grant-in-Aids JP20K22367, JP20H05643, JP21K13957, and JP22K02956, JSPS Overseas Challenge Program for Young Researchers, JSPS Overseas Challenge Program for Young Researchers, the ISEE director’s leadership fund for FYs 2021 – 2023, Young Leader Cultivation (YLC) programme of Nagoya University, Tokai Pathways to Global Excellence (Nagoya University) of the Strategic Professional Development Program for Young Researchers (MEXT), and the young researcher units for the advancement of new and undeveloped fields by Institute for Advanced Research of Nagoya University (the Program for Promoting the Enhancement of Research Universities).
Author information
Authors and Affiliations
Contributions
MO wrote the manuscript with input and contributions from all co-authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Ken McCracken is retired.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Owens, M.J., Lockwood, M., Barnard, L.A. et al. Reconstructing Sunspot Number by Forward-Modelling Open Solar Flux. Sol Phys 299, 3 (2024). https://doi.org/10.1007/s11207-023-02241-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11207-023-02241-3