-
A Multi-Model Ensemble System for the outer Heliosphere (MMESH): Solar Wind Conditions near Jupiter
Authors:
M. J. Rutala,
C. M. Jackman,
M. J. Owens,
C. Tao,
A. R. Fogg,
S. A. Murray,
L. Barnard
Abstract:
How the solar wind influences the magnetospheres of the outer planets is a fundamentally important question, but is difficult to answer in the absence of consistent, simultaneous monitoring of the upstream solar wind and the large-scale dynamics internal to the magnetosphere. To compensate for the relative lack of in-situ data, propagation models are often used to estimate the ambient solar wind c…
▽ More
How the solar wind influences the magnetospheres of the outer planets is a fundamentally important question, but is difficult to answer in the absence of consistent, simultaneous monitoring of the upstream solar wind and the large-scale dynamics internal to the magnetosphere. To compensate for the relative lack of in-situ data, propagation models are often used to estimate the ambient solar wind conditions at the outer planets for comparison to remote observations or in-situ measurements. This introduces another complication: the propagation of near-Earth solar wind measurements introduces difficult-to-assess uncertainties. Here, we present the Multi-Model Ensemble System for the outer Heliosphere (MMESH) to begin to address these issues, along with the resultant multi-model ensemble (MME) of the solar wind conditions near Jupiter. MMESH accepts as input any number of solar wind models together with contemporaneous in-situ spacecraft data. From these, the system characterizes typical uncertainties in model timing, quantifies how these uncertainties vary under different conditions, attempts to correct for systematic biases in the input model timing, and composes a MME with uncertainties from the results. For the case of the Jupiter-MME presented here, three solar wind propagation models were compared to in-situ measurements from the near-Jupiter spacecraft Ulysses and Juno which span diverse geometries and phases of the solar cycle, amounting to more than 14,000 hours of data over 2.5 decades. The MME gives the most-probable near-Jupiter solar wind conditions for times within the tested epoch, outperforming the input models and returning quantified estimates of uncertainty.
△ Less
Submitted 29 February, 2024;
originally announced February 2024.
-
CME Propagation Through the Heliosphere: Status and Future of Observations and Model Development
Authors:
M. Temmer,
C. Scolini,
I. G. Richardson,
S. G. Heinemann,
E. Paouris,
A. Vourlidas,
M. M. Bisi,
writing teams,
:,
N. Al-Haddad,
T. Amerstorfer,
L. Barnard,
D. Buresova,
S. J. Hofmeister,
K. Iwai,
B. V. Jackson,
R. Jarolim,
L. K. Jian,
J. A. Linker,
N. Lugaz,
P. K. Manoharan,
M. L. Mays,
W. Mishra,
M. J. Owens,
E. Palmerio
, et al. (9 additional authors not shown)
Abstract:
The ISWAT clusters H1+H2 have a focus on interplanetary space and its characteristics, especially on the large-scale co-rotating and transient structures impacting Earth. SIRs, generated by the interaction between high-speed solar wind originating in large-scale open coronal magnetic fields and slower solar wind from closed magnetic fields, are regions of compressed plasma and magnetic field follo…
▽ More
The ISWAT clusters H1+H2 have a focus on interplanetary space and its characteristics, especially on the large-scale co-rotating and transient structures impacting Earth. SIRs, generated by the interaction between high-speed solar wind originating in large-scale open coronal magnetic fields and slower solar wind from closed magnetic fields, are regions of compressed plasma and magnetic field followed by high-speed streams that recur at the ca. 27 day solar rotation period. Short-term reconfigurations of the lower coronal magnetic field generate flare emissions and provide the energy to accelerate enormous amounts of magnetised plasma and particles in the form of CMEs into interplanetary space. The dynamic interplay between these phenomena changes the configuration of interplanetary space on various temporal and spatial scales which in turn influences the propagation of individual structures. While considerable efforts have been made to model the solar wind, we outline the limitations arising from the rather large uncertainties in parameters inferred from observations that make reliable predictions of the structures impacting Earth difficult. Moreover, the increased complexity of interplanetary space as solar activity rises in cycle 25 is likely to pose a challenge to these models. Combining observational and modeling expertise will extend our knowledge of the relationship between these different phenomena and the underlying physical processes, leading to improved models and scientific understanding and more-reliable space-weather forecasting. The current paper summarizes the efforts and progress achieved in recent years, identifies open questions, and gives an outlook for the next 5-10 years. It acts as basis for updating the existing COSPAR roadmap by Schrijver+ (2015), as well as providing a useful and practical guide for peer-users and the next generation of space weather scientists.
△ Less
Submitted 9 August, 2023;
originally announced August 2023.
-
SIR-HUXt -- a particle filter data assimilation scheme for assimilating CME time-elongation profiles
Authors:
Luke Barnard,
Mathew Owens,
Chris Scott,
Matthew Lang,
Mike Lockwood
Abstract:
We present the development of SIR-HUXt, the integration of a sequential importance resampling (SIR) data assimilation scheme with the HUXt solar wind model. SIR-HUXt is designed to assimilate the time-elongation profiles of CME fronts in the low heliosphere, such as those typically extracted from heliospheric imager data returned by the STEREO, Parker Solar Probe, and Solar Orbiter missions. We us…
▽ More
We present the development of SIR-HUXt, the integration of a sequential importance resampling (SIR) data assimilation scheme with the HUXt solar wind model. SIR-HUXt is designed to assimilate the time-elongation profiles of CME fronts in the low heliosphere, such as those typically extracted from heliospheric imager data returned by the STEREO, Parker Solar Probe, and Solar Orbiter missions. We use Observing System Simulation Experiments to explore the performance of SIR-HUXt for a simple synthetic CME scenario of a fully Earth directed CME flowing through a uniform ambient solar wind, where the CME is initialised with the average observed CME speed and width. These experiments are performed for a range of observer locations, from 20 deg to 90 deg behind Earth, spanning the L5 point where ESA's future Vigil space weather monitor will return heliospheric imager data for operational space weather forecasting.
We show that SIR-HUXt performs well at constraining the CME speed, and has some success at constraining the CME longitude. The CME width is largely unconstrained by the SIR-HUXt assimilations, and more experiments are required to determine if this is due to this specific CME scenario, or is a general feature of assimilating time-elongation profiles. Rank-histograms suggest that the SIR-HUXt ensembles are well calibrated, with no clear indications of bias or under/over dispersion. Improved constraints on the initial CME speed lead directly to improvements in the CME transit time to Earth and arrival speed. For an observer in the L5 region, SIR-HUXt returned a 69% reduction in the CME transit time uncertainty, and a 63% reduction in the arrival speed uncertainty. This suggests SIR-HUXt has potential to improve the real-world representivity of HUXt simulations, and therefore has potential to reduce the uncertainty of CME arrival time hindcasts and forecasts.
△ Less
Submitted 5 October, 2022;
originally announced October 2022.
-
HUXt -- An open source, computationally efficient reduced-physics solar wind model, written in Python
Authors:
Luke Barnard,
Mathew Owens
Abstract:
HUXt is an open source numerical model of the solar wind written in Python. It is based on the solution of the 1D inviscid Burger's equation. This reduced-physics approach produces solar wind flow simulations that closely emulate the flow produced by 3-D magnetohydrodynamic solar wind models at a small fraction of the computational expense. While not intended as a replacement for 3-D MHD, the simp…
▽ More
HUXt is an open source numerical model of the solar wind written in Python. It is based on the solution of the 1D inviscid Burger's equation. This reduced-physics approach produces solar wind flow simulations that closely emulate the flow produced by 3-D magnetohydrodynamic solar wind models at a small fraction of the computational expense. While not intended as a replacement for 3-D MHD, the simplicity and computational efficiency of HUXt offers several key advantages that enable experiments and the use of techniques that would otherwise be cost prohibitive. For example, large ensembles can easily be run with modest computing resources, which are useful for exploring and quantifying the uncertainty in space weather predictions, as well as for the application of some data assimilation methods.
We present the developments in the latest version of HUXt, v4.0, and discuss our plans for future developments and applications of the model. The three key developments in v4.0 are: a restructuring of the models solver to enable fully time-dependent boundary conditions, such that HUXt can in principle be initialised with in-situ observations from any of the fleet of heliospheric monitors; new functionality to trace streaklines through the HUXt flow solutions, which can be used to track features such as the Heliospheric Current Sheet; introduction of a small test-suite so that we can better ensure the reliability and reproducibility of HUXt simulations for all users across future versions. Other more minor developments are discussed in the article.
Future applications of HUXt are discussed, including the development data assimilation schemes for assimilation of both remote sensing and in-situ plasma measures. We discuss the progress of transitioning HUXt into an operational model at the UK's Met Office Space Weather Operations Center as part of the UK governments SWIMMR programme.
△ Less
Submitted 2 October, 2022;
originally announced October 2022.
-
Solar Energetic Particle Ground-Level Enhancements and the Solar Cycle
Authors:
Mathew Owens,
Luke Barnard,
Benjamin Pope,
Mike Lockwood,
Ilya Usoskin,
Eleanna Asvestari
Abstract:
Severe geomagnetic storms appear to be ordered by the solar cycle in a number of ways. They occur more frequently close to solar maximum and declining phase, are more common in larger solar cycles and show different patterns of occurrence in odd- and even-numbered solar cycles. Our knowledge of the most extreme space weather events, however, comes from the spikes in cosmogenic-isotope ($^{14}$C,…
▽ More
Severe geomagnetic storms appear to be ordered by the solar cycle in a number of ways. They occur more frequently close to solar maximum and declining phase, are more common in larger solar cycles and show different patterns of occurrence in odd- and even-numbered solar cycles. Our knowledge of the most extreme space weather events, however, comes from the spikes in cosmogenic-isotope ($^{14}$C, $^{10}$Be and $^{36}$Cl) records that are attributed to significantly larger solar energetic particle (SEP) events than have been observed during the space age. Despite both storms and SEPs being driven by solar eruptive phenomena, the event-by-event correspondence between extreme storms and extreme SEPs is low. Thus it should not be assumed a priori that the solar cycle patterns found for storms also hold for SEPs and the cosmogenic-isotope events. In this study we investigate the solar cycle trends in the timing and magnitude of the 67 SEP ground-level enhancements (GLEs) recorded by neutron monitors since the mid 1950s. Using a number of models of GLE occurrence probability, we show that GLEs are around a factor four more likely around solar maximum than around solar minimum, and that they preferentially occur earlier in even-numbered solar cycles than in odd-numbered cycles. There are insufficient data to conclusively determine whether larger solar cycles produce more GLEs. Implications for putative space-weather events in the cosmogenic-isotope records are discussed. We find that GLEs tend to cluster within a few tens of days, likely due to particularly productive individual active regions, and with approximately 11-year separations, owing to the solar cycle ordering. But these timescales do not explain cosmogenic-isotope spikes which require multiple extreme SEP events over consecutive years.
△ Less
Submitted 27 July, 2022; v1 submitted 26 July, 2022;
originally announced July 2022.
-
Quantifying the uncertainty in CME kinematics derived from geometric modelling of Heliospheric Imager data
Authors:
Luke Barnard,
Mathew Owens,
Christopher J. Scott,
Mike Lockwood,
Curt A. de Koning,
Tanja Amerstorfer,
Jürgen Hinterreiter,
Christian Möstl,
Jackie Davies,
Pete Riley
Abstract:
Geometric modelling of Coronal Mass Ejections (CMEs) is a widely used tool for assessing their kinematic evolution. Furthermore, techniques based on geometric modelling, such as ELEvoHI, are being developed into forecast tools for space weather prediction. These models assume that solar wind structure does not affect the evolution of the CME, which is an unquantified source of uncertainty. We use…
▽ More
Geometric modelling of Coronal Mass Ejections (CMEs) is a widely used tool for assessing their kinematic evolution. Furthermore, techniques based on geometric modelling, such as ELEvoHI, are being developed into forecast tools for space weather prediction. These models assume that solar wind structure does not affect the evolution of the CME, which is an unquantified source of uncertainty. We use a large number of Cone CME simulations with the HUXt solar wind model to quantify the scale of uncertainty introduced into geometric modelling and the ELEvoHI CME arrival times by solar wind structure. We produce a database of simulations, representing an average, a fast, and an extreme CME scenario, each independently propagating through 100 different ambient solar wind environments. Synthetic heliospheric imager observations of these simulations are then used with a range of geometric models to estimate the CME kinematics. The errors of geometric modelling depend on the location of the observer, but do not seem to depend on the CME scenario. In general, geometric models are biased towards predicting CME apex distances that are larger than the true value. For these CME scenarios, geometric modelling errors are minimised for an observer in the L5 region. Furthermore, geometric modelling errors increase with the level of solar wind structure in the path of the CME. The ELEvoHI arrival time errors are minimised for an observer in the L5 region, with mean absolute arrival time errors of $8.2\pm1.2$~h, $8.3\pm1.0$~h, and $5.8\pm0.9$~h for the average, fast, and extreme CME scenarios.
△ Less
Submitted 26 November, 2021;
originally announced November 2021.
-
Drag-based CME modeling with heliospheric images incorporating frontal deformation: ELEvoHI 2.0
Authors:
J. Hinterreiter,
T. Amerstorfer,
M. Temmer,
M. A. Reiss,
A. J. Weiss,
C. Möstl,
L. A. Barnard,
J. Pomoell,
M. Bauer,
U. V. Amerstorfer
Abstract:
The evolution and propagation of coronal mass ejections (CMEs) in interplanetary space is still not well understood. As a consequence, accurate arrival time and arrival speed forecasts are an unsolved problem in space weather research. In this study, we present the ELlipse Evolution model based on HI observations (ELEvoHI) and introduce a deformable front to this model. ELEvoHI relies on heliosphe…
▽ More
The evolution and propagation of coronal mass ejections (CMEs) in interplanetary space is still not well understood. As a consequence, accurate arrival time and arrival speed forecasts are an unsolved problem in space weather research. In this study, we present the ELlipse Evolution model based on HI observations (ELEvoHI) and introduce a deformable front to this model. ELEvoHI relies on heliospheric imagers (HI) observations to obtain the kinematics of a CME. With the newly developed deformable front, the model is able to react to the ambient solar wind conditions during the entire propagation and along the whole front of the CME. To get an estimate of the ambient solar wind conditions, we make use of three different models: Heliospheric Upwind eXtrapolation model (HUX), Heliospheric Upwind eXtrapolation with time dependence model (HUXt), and EUropean Heliospheric FORecasting Information Asset (EUHFORIA). We test the deformable front on a CME first observed in STEREO-A/HI on February 3, 2010 14:49 UT. For this case study, the deformable front provides better estimates of the arrival time and arrival speed than the original version of ELEvoHI using an elliptical front. The new implementation enables us to study the parameters influencing the propagation of the CME not only for the apex, but for the entire front. The evolution of the CME front, especially at the flanks, is highly dependent on the ambient solar wind model used. An additional advantage of the new implementation is given by the possibility to provide estimates of the CME mass.
△ Less
Submitted 18 August, 2021;
originally announced August 2021.
-
Why are ELEvoHI CME arrival predictions different if based on STEREO-A or STEREO-B heliospheric imager observations?
Authors:
Jürgen Hinterreiter,
Tanja Amerstorfer,
Martin A. Reiss,
Christian Möstl,
Manuela Temmer,
Maike Bauer,
Ute V. Amerstorfer,
Rachel L. Bailey,
Andreas J. Weiss,
Jackie A. Davies,
Luke A. Barnard,
Mathew J. Owens
Abstract:
Accurate forecasting of the arrival time and arrival speed of coronal mass ejections (CMEs) is a unsolved problem in space weather research. In this study, a comparison of the predicted arrival times and speeds for each CME based, independently, on the inputs from the two STEREO vantage points is carried out. We perform hindcasts using ELlipse Evolution model based on Heliospheric Imager observati…
▽ More
Accurate forecasting of the arrival time and arrival speed of coronal mass ejections (CMEs) is a unsolved problem in space weather research. In this study, a comparison of the predicted arrival times and speeds for each CME based, independently, on the inputs from the two STEREO vantage points is carried out. We perform hindcasts using ELlipse Evolution model based on Heliospheric Imager observations (ELEvoHI) ensemble modelling. An estimate of the ambient solar wind conditions is obtained by the Wang-Sheeley-Arge/Heliospheric Upwind eXtrapolation (WSA/HUX) model combination that serves as input to ELEvoHI. We carefully select 12 CMEs between February 2010 and July 2012 that show clear signatures in both STEREO-A and STEREO-B HI time-elongation maps, that propagate close to the ecliptic plane, and that have corresponding in situ signatures at Earth. We find a mean arrival time difference of 6.5 hrs between predictions from the two different viewpoints, which can reach up to 9.5 hrs for individual CMEs, while the mean arrival speed difference is 63 km s$^{-1}$. An ambient solar wind with a large speed variance leads to larger differences in the STEREO-A and STEREO-B CME arrival time predictions ($cc~=~0.92$). Additionally, we compare the predicted arrivals, from both spacecraft, to the actual in situ arrivals at Earth and find a mean absolute error of 7.5 $\pm$ 9.5 hrs for the arrival time and 87 $\pm$ 111 km s$^{-1}$ for the arrival speed. There is no tendency for one spacecraft to provide more accurate arrival predictions than the other.
△ Less
Submitted 15 February, 2021;
originally announced February 2021.
-
Semi-annual, annual and Universal Time variations in the magnetosphere and in geomagnetic activity: 4. Polar Cap motions and origins of the Universal Time effect
Authors:
Mike Lockwood,
Carl Haines,
Luke A. Barnard,
Mathew J. Owens,
Chris J. Scott,
Aude Chambodut,
Kathryn A. McWilliams
Abstract:
We use the am, an, as and the a-sigma geomagnetic indices to the explore a previously overlooked factor in magnetospheric electrodynamics, namely the inductive effect of diurnal motions of the Earth's magnetic poles toward and away from the Sun caused by Earth's rotation. Because the offset of the (eccentric dipole) geomagnetic pole from the rotational axis is roughly twice as large in the souther…
▽ More
We use the am, an, as and the a-sigma geomagnetic indices to the explore a previously overlooked factor in magnetospheric electrodynamics, namely the inductive effect of diurnal motions of the Earth's magnetic poles toward and away from the Sun caused by Earth's rotation. Because the offset of the (eccentric dipole) geomagnetic pole from the rotational axis is roughly twice as large in the southern hemisphere compared to the northern, the effects there are predicted to be roughly twice the amplitude. Hemispheric differences have previously been discussed in terms of polar ionospheric conductivities, effects which we allow for by studying the dipole tilt effect on time-of-year variations of the indices. The electric field induced in a geocentric frame is shown to also be a significant factor and gives a modulation of the voltage applied by the solar wind flow in the southern hemisphere of typically a 30% diurnal modulation for disturbed intervals rising to 76% in quiet times. Motion towards/away from the Sun reduces/enhances the directly-driven ionospheric voltages and reduces/enhances the magnetic energy stored in the near-Earth tail: 10% of the effect being directly-driven and 90% being in tail energy storage/release. Combined with the effect of solar wind dynamic pressure and dipole tilt on the pressure balance in the near-Earth tail, the effect provides an excellent explanation of how the observed Russell-McPherron pattern in the driving power input into the magnetosphere is converted into the equinoctial pattern in average geomagnetic activity (after correction is made for dipole tilt effects on ionospheric conductivity), added to a pronounced UT variation with minimum at 02-10UT. In addition, we show that the predicted and observed UT variations in average geomagnetic activity has implications for the occurrence of the largest events that also show the nett UT variation.
△ Less
Submitted 24 December, 2020;
originally announced December 2020.
-
A homogeneous aa index: 2. hemispheric asymmetries and the equinoctial variation
Authors:
Mike Lockwood,
Ivan D. Finch,
Aude Chambodut,
Luke A. Barnard,
Mathew J. Owens,
Ellen Clarke
Abstract:
Paper 1 [Lockwood et al., 2018] generated annual means of a new version of the $aa$ geomagnetic activity index which includes corrections for secular drift in the geographic coordinates of the auroral oval, thereby resolving the difference between the centennial-scale change in the northern and southern hemisphere indices, $aa_N$ and $aa_S$. However, other hemispheric asymmetries in the $aa$ index…
▽ More
Paper 1 [Lockwood et al., 2018] generated annual means of a new version of the $aa$ geomagnetic activity index which includes corrections for secular drift in the geographic coordinates of the auroral oval, thereby resolving the difference between the centennial-scale change in the northern and southern hemisphere indices, $aa_N$ and $aa_S$. However, other hemispheric asymmetries in the $aa$ index remain: in particular, the distributions of 3-hourly $aa_N$ and $aa_S$ values are different and the correlation between them is not high on this timescale ($r = 0.66$). In the present paper, a location-dependant station sensitivity model is developed using the $am$ index (derived from a much more extensive network of stations in both hemispheres) and used to reduce the difference between the hemispheric $aa$ indices and improve their correlation (to $r = 0.79$) by generating corrected 3-hourly hemispheric indices, $aa_{HN}$ and $aa_{HS}$, which also include the secular drift corrections detailed in Paper 1. These are combined into a new, 'homogeneous' $aa$ index, $aa_H$. It is shown that $aa_H$, unlike $aa$, reveals the 'equinoctial'-like time-of-day/time-of-year pattern that is found for the $am$ index.
△ Less
Submitted 17 December, 2018; v1 submitted 24 November, 2018;
originally announced November 2018.
-
A homogeneous aa index: 1. Secular variation
Authors:
Mike Lockwood,
Aude Chambodut,
Luke A. Barnard,
Mathew J. Owens,
Ellen Clarke,
Véronique Mendel
Abstract:
Originally complied for 1868-1967 and subsequently continued so that it now covers 150 years, the $aa$ index has become a vital resource for studying space climate change. However, there have been debates about the inter-calibration of data from the different stations. In addition, the effects of secular change in the geomagnetic field have not previously been allowed for. As a result, the compone…
▽ More
Originally complied for 1868-1967 and subsequently continued so that it now covers 150 years, the $aa$ index has become a vital resource for studying space climate change. However, there have been debates about the inter-calibration of data from the different stations. In addition, the effects of secular change in the geomagnetic field have not previously been allowed for. As a result, the components of the 'classical' $aa$ index for the southern and northern hemispheres ($aa_S$ and $aa_N$) have drifted apart. We here separately correct both $aa_S$ and $aa_N$ for both these effects using the same method as used to generate the classic $aa$ values but allowing $δ$, the minimum angular separation of each station from a nominal auroral oval, to vary as calculated using the IGRF-12 and gufm1 models of the intrinsic geomagnetic field. Our approach is to correct the quantized aK-values for each station, originally scaled on the assumption that $δ$ values are constant, with time-dependent scale factors that allow for the drift in $δ$. This requires revisiting the intercalibration of successive stations used in making the $aa_S$ and $aa_N$ composites. These intercalibrations are defined using independent data and daily averages from 11 years before and after each station change and it is shown that they depend on the time of year. This procedure produces new homogenized hemispheric aa indices, $aa_{HS}$ and $aa_{HN}$, which show centennial-scale changes that are in very close agreement. Calibration problems with the classic $aa$ index are shown to have arisen from drifts in $δ$ combined with simpler corrections which gave an incorrect temporal variation and underestimate the rise in $aa$ during the 20th century by about 15%.
△ Less
Submitted 17 December, 2018; v1 submitted 24 November, 2018;
originally announced November 2018.
-
Space climate and space weather over the past 400 years: 2. Proxy indicators of geomagnetic storm and substorm
Authors:
Mike Lockwood,
Mathew J. Owens,
Luke A. Barnard,
Chris J. Scott,
Clare E. Watt,
Sarah Bentley
Abstract:
Using the reconstruction of power input to the magnetosphere given in Paper 1 (arXiv:1708.04904), we reconstruct annual means of geomagnetic indices over the past 400 years to within a 1-sigma error of +/-20 pc. In addition, we study the behaviour of the lognormal distribution of daily and hourly values about these annual means and show that we can also reconstruct the fraction of geomagnetically-…
▽ More
Using the reconstruction of power input to the magnetosphere given in Paper 1 (arXiv:1708.04904), we reconstruct annual means of geomagnetic indices over the past 400 years to within a 1-sigma error of +/-20 pc. In addition, we study the behaviour of the lognormal distribution of daily and hourly values about these annual means and show that we can also reconstruct the fraction of geomagnetically-active (storm-like) days and (substorm-like) hours in each year to accuracies of 50-60 pc. The results are the first physics-based quantification of the space weather conditions in both the Dalton and Maunder minima. We predict terrestrial disturbance levels in future repeats of these minima, allowing for the weakening of Earth's dipole moment.
△ Less
Submitted 8 January, 2018;
originally announced January 2018.
-
Space Climate and Space Weather over the past 400 years: 1. The Power input to the Magnetosphere
Authors:
Mike Lockwood,
Mathew J. Owens,
Luke A. Barnard,
Chris J. Scott,
Clare E. Watt
Abstract:
Using information on geomagnetic activity, sunspot numbers and cosmogenic isotopes, supported by historic eclipse images and in conjunction with models, it has been possible to reconstruct annual means of solar wind speed and number density and heliospheric magnetic field (HMF) intensity since 1611, when telescopic observations of sunspots began. These models are developed and tuned using data rec…
▽ More
Using information on geomagnetic activity, sunspot numbers and cosmogenic isotopes, supported by historic eclipse images and in conjunction with models, it has been possible to reconstruct annual means of solar wind speed and number density and heliospheric magnetic field (HMF) intensity since 1611, when telescopic observations of sunspots began. These models are developed and tuned using data recorded by near-Earth interplanetary spacecraft and by solar magnetograms over the past 53 years. In this paper, we use these reconstructions to quantify power input into the magnetosphere over the past 400 years. For each year, both the annual mean power input is computed and its distribution in daily means. This is possible because the distribution of daily values divided by the annual mean is shown to maintain the same lognormal form with a constant variance. This study is another important step towards the development of a physics-based, long-term climatology of space weather conditions.
△ Less
Submitted 17 August, 2017; v1 submitted 16 August, 2017;
originally announced August 2017.
-
Tests of sunspot number sequences: 4. Discontinuities around 1946 in various sunspot number and sunspot group number reconstructions
Authors:
Mike Lockwood,
Mathew J. Owens,
Luke A. Barnard
Abstract:
We use 5 test data series to quantify putative discontinuities around 1946 in 5 annual-mean sunspot number or group number sequences. The series tested are: the original and new versions of the Wolf/Zurich/International sunspot number composite [$R_{ISNv1}$ and $R_{ISNv2}$] ; the corrected version of $R_{ISNv1}$ [$R_C$]; the backbone group number [$R_{BB}$]; and the group number composite […
▽ More
We use 5 test data series to quantify putative discontinuities around 1946 in 5 annual-mean sunspot number or group number sequences. The series tested are: the original and new versions of the Wolf/Zurich/International sunspot number composite [$R_{ISNv1}$ and $R_{ISNv2}$] ; the corrected version of $R_{ISNv1}$ [$R_C$]; the backbone group number [$R_{BB}$]; and the group number composite [$R_{UEA}$]. The test data are: the group number $N_G$ and total sunspot area $A_G$ from the RGO photoheliographic data; the CaK index from re-analysis of MWO CaII K spectroheliograms; the group number from the MWO sunspot drawings, $N_{MWO}$; and ionospheric critical frequencies at Slough [$foF2$]. The test data all vary with sunspot numbers, in some cases non-linearly. Tests use both before-and-after fit-residual comparison and correlation methods, applied to intervals iterated to minimise errors and eliminate the effect of the discontinuity date. It is not assumed that the correction required is by a constant factor, nor even linear in sunspot number. A non-linear correction is required by $R_C$, $R_{BB}$ and $R_{ISNv1}$, but not by $R_{ISNv2}$ or $R_{UEA}$. The test datasets give very similar results in all cases. By multiplying the probability distribution functions together we obtain the optimum correction for each data series that must be applied to pre-discontinuity data to make them consistent with the post-discontinuity data. It is shown that, on average, values for 1932-1943 are too small (relative to later values) by about 12.3% for $R_{ISNv1}$ but are too large for $R_{ISNv2}$ and $R_{BB}$ by 3.8% and 5.2%. The correction applied to generate $R_C$ from $R_{ISNv1}$ reduces this average factor to 0.5% but does not remove the non-linear variation, and other errors remain uncorrected. A test is provided by $R_{UEA}$, which is identical to the RGO $N_G$ values over the interval used.
△ Less
Submitted 17 May, 2016;
originally announced May 2016.
-
Tests of Sunspot Number Sequences: 2. Using Geomagnetic and Auroral Data
Authors:
Mike Lockwood,
Mathew J. Owens,
Luke A Barnard,
Chris J. Scott,
Ilya G. Usoskin,
Heikki Nevanlinna
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 6 May, 2016;
originally announced May 2016.
-
An Assessment of Sunspot Number Data Composites over 1845-2014
Authors:
Mike Lockwood,
Mathew J. Owens,
Luke A. Barnard,
Ilya G. Usoskin
Abstract:
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 in…
▽ More
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.
△ Less
Submitted 15 April, 2016;
originally announced April 2016.
-
Tests of Sunspot Number Sequences: 1. Using Ionosonde Data
Authors:
M. Lockwood,
C. J. Scott,
M. J. Owens,
L. Barnard,
D. M. Willis
Abstract:
More than 70 years ago it was recognised that ionospheric F2-layer critical frequencies $foF2$ had a strong relationship to sunspot number. Using historic datasets from the Slough and Washington ionosondes, we evaluate the best statistical fits of $foF2$ to sunspot numbers (at each Universal Time [UT] separately) in order to search for drifts and abrupt changes in the fit residuals over Solar Cycl…
▽ More
More than 70 years ago it was recognised that ionospheric F2-layer critical frequencies $foF2$ had a strong relationship to sunspot number. Using historic datasets from the Slough and Washington ionosondes, we evaluate the best statistical fits of $foF2$ to sunspot numbers (at each Universal Time [UT] separately) in order to search for drifts and abrupt changes in the fit residuals over Solar Cycles 17 - 21. Polynomial fits are made both with and without allowance for the white-light facular area, which has been reported as being associated with cycle-to-cycle changes in the sunspot number - $foF2$ relationship. Over the interval studied here, the ISN, $R$, the backbone group number $Rbb$, and the corrected number $Rc$ largely differ in their allowance for the 'Waldmeier discontinuity' around 1945 (the correction factor for which for $R$, $Rbb$ and $Rc$ is, respectively, zero, effectively over 20%, and explicitly 11.6%). It is shown that for Solar Cycles 18 - 21, all three sunspot data sequences perform well, but that the fit residuals are lowest and most uniform for $Rbb$. We here use $foF2$ for those UTs for which $R$, $Rbb$, and $Rc$ all give correlations exceeding 0.99 for intervals both before and after the Waldmeier discontinuity. The error introduced by the Waldmeier discontinuity causes $R$ to underestimate the fitted values based on the $foF2$ data for 1932 - 1945 but $Rbb$ overestimates them by almost the same factor, implying that the correction for the Waldmeier discontinuity inherent in $Rbb$ is too large by a factor of two. Fit residuals are smallest and most uniform for $Rc$ and the ionospheric data support the optimum discontinuity multiplicative correction factor derived from the independent Royal Greenwich Observatory (RGO) sunspot group data for the same interval.
△ Less
Submitted 24 January, 2016;
originally announced January 2016.
-
Tests of sunspot number sequences: 3. Effects of regression procedures on the calibration of historic sunspot data
Authors:
M. Lockwood,
M. J Owens,
L. Barnard,
I. G. Usoskin
Abstract:
We use sunspot group observations from the Royal Greenwich Observatory (RGO) to investigate the effects of intercalibrating data from observers with different visual acuities. The tests are made by counting the number of groups $R_B$ above a variable cut-off threshold of observed total whole-spot area (uncorrected for foreshortening) to simulate what a lower acuity observer would have seen. The sy…
▽ More
We use sunspot group observations from the Royal Greenwich Observatory (RGO) to investigate the effects of intercalibrating data from observers with different visual acuities. The tests are made by counting the number of groups $R_B$ above a variable cut-off threshold of observed total whole-spot area (uncorrected for foreshortening) to simulate what a lower acuity observer would have seen. The synthesised annual means of $R_B$ are then re-scaled to the observed RGO group number $R_A$ using a variety of regression techniques. It is found that a very high correlation between $R_A$ and $R_B$ ($r_{AB}$ > 0.98) does not prevent large errors in the intercalibration (e.g. sunspot maximum values can be over 30% too large even for such levels of $r_{AB}$). In generating the backbone sunspot number, Svalgaard and Schatten [2015] force regression fits to pass through the scatter plot origin which generates unreliable fits (the residuals do not form a normal distribution) and causes sunspot cycle amplitudes to be exaggerated in the intercalibrated data. It is demonstrated that the use of Quantile-Quantile (Q-Q) plots to test for a normal distribution is a useful indicator of erroneous and misleading regression fits. Ordinary least squares linear fits, not forced to pass through the origin, are sometimes reliable (although the optimum method used is shown to be different when matching peak and average sunspot group numbers). However other fits are only reliable if non-linear regression is used. From these results it is entirely possible that the inflation of solar cycle amplitudes in the backbone group sunspot number as one goes back in time, relative to related solar-terrestrial parameters, is entirely caused by the use of inappropriate and non-robust regression techniques to calibrate the sunspot data.
△ Less
Submitted 24 January, 2016; v1 submitted 27 October, 2015;
originally announced October 2015.