-
Slow Solar Wind Connection Science during Solar Orbiter's First Close Perihelion Passage
Authors:
Stephanie L. Yardley,
Christopher J. Owen,
David M. Long,
Deborah Baker,
David H. Brooks,
Vanessa Polito,
Lucie M. Green,
Sarah Matthews,
Mathew Owens,
Mike Lockwood,
David Stansby,
Alexander W. James,
Gherado Valori,
Alessandra Giunta,
Miho Janvier,
Nawin Ngampoopun,
Teodora Mihailescu,
Andy S. H. To,
Lidia van Driel-Gesztelyi,
Pascal Demoulin,
Raffaella D'Amicis,
Ryan J. French,
Gabriel H. H. Suen,
Alexis P. Roulliard,
Rui F. Pinto
, et al. (54 additional authors not shown)
Abstract:
The Slow Solar Wind Connection Solar Orbiter Observing Plan (Slow Wind SOOP) was developed to utilise the extensive suite of remote sensing and in situ instruments on board the ESA/NASA Solar Orbiter mission to answer significant outstanding questions regarding the origin and formation of the slow solar wind. The Slow Wind SOOP was designed to link remote sensing and in situ measurements of slow w…
▽ More
The Slow Solar Wind Connection Solar Orbiter Observing Plan (Slow Wind SOOP) was developed to utilise the extensive suite of remote sensing and in situ instruments on board the ESA/NASA Solar Orbiter mission to answer significant outstanding questions regarding the origin and formation of the slow solar wind. The Slow Wind SOOP was designed to link remote sensing and in situ measurements of slow wind originating at open-closed field boundaries. The SOOP ran just prior to Solar Orbiter's first close perihelion passage during two remote sensing windows (RSW1 and RSW2) between 2022 March 3-6 and 2022 March 17-22, while Solar Orbiter was at a heliocentric distance of 0.55-0.51 and 0.38-0.34 au from the Sun, respectively. Coordinated observation campaigns were also conducted by Hinode and IRIS. The magnetic connectivity tool was used, along with low latency in situ data, and full-disk remote sensing observations, to guide the target pointing of Solar Orbiter. Solar Orbiter targeted an active region complex during RSW1, the boundary of a coronal hole, and the periphery of a decayed active region during RSW2. Post-observation analysis using the magnetic connectivity tool along with in situ measurements from MAG and SWA/PAS, show that slow solar wind, with velocities between 210 and 600 km/s, arrived at the spacecraft originating from two out of the three of the target regions. The Slow Wind SOOP, despite presenting many challenges, was very successful, providing a blueprint for planning future observation campaigns that rely on the magnetic connectivity of Solar Orbiter.
△ Less
Submitted 20 April, 2023; v1 submitted 19 April, 2023;
originally announced April 2023.
-
Re-calibration of the Sunspot Number: Status Report
Authors:
F. Clette,
L. Lefèvre,
T. Chatzistergos,
H. Hayakawa,
V. M. Carrasco,
R. Arlt,
E. W. Cliver,
T. Dudok de Wit,
T. Friedli,
N. Karachik,
G. Kopp,
M. Lockwood,
S. Mathieu,
A. Muñoz-Jaramillo,
M. Owens,
D. Pesnell,
A. Pevtsov,
L. Svalgaard,
I. G. Usoskin,
L. van Driel-Gesztelyi,
J. M. Vaquero
Abstract:
We report progress on the ongoing recalibration of the Wolf sunspot number (SN) and Group sunspot number (GN) following the release of version 2.0 of SN in 2015. This report constitutes both an update of the efforts reported in the 2016 Topical Issue of Solar Physics and a summary of work by the International Space Science Institute (ISSI) International Team formed in 2017 to develop optimal SN an…
▽ More
We report progress on the ongoing recalibration of the Wolf sunspot number (SN) and Group sunspot number (GN) following the release of version 2.0 of SN in 2015. This report constitutes both an update of the efforts reported in the 2016 Topical Issue of Solar Physics and a summary of work by the International Space Science Institute (ISSI) International Team formed in 2017 to develop optimal SN and GN re-construction methods while continuing to expand the historical sunspot number database. Significant progress has been made on the database side while more work is needed to bring the various proposed SN and (primarily) GN reconstruction methods closer to maturity, after which the new reconstructions (or combinations thereof) can be compared with (a) ``benchmark'' expectations for any normalization scheme (e.g., a general increase in observer normalization factors going back in time), and (b) independent proxy data series such as F10.7 and the daily range of variations of Earth's undisturbed magnetic field. New versions of the underlying databases for SN and GN will shortly become available for years through 2022 and we anticipate the release of next versions of these two time series in 2024.
△ Less
Submitted 5 January, 2023;
originally announced January 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.
-
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.
-
Cosmic Meteorology
Authors:
Mike Lockwood,
Mat Owens
Abstract:
Mike Lockwood and Mathew Owens discuss how eclipse observations are aiding the development of a climatology of near-Earth space
Mike Lockwood and Mathew Owens discuss how eclipse observations are aiding the development of a climatology of near-Earth space
△ Less
Submitted 25 May, 2021;
originally announced May 2021.
-
Evolving Solar Wind Flow Properties of Magnetic Inversions Observed by Helios
Authors:
Allan R Macneil,
Mathew J Owens,
Robert T Wicks,
Mike Lockwood
Abstract:
In its first encounter at solar distances as close as r = 0.16AU, Parker Solar Probe (PSP) observed numerous local reversals, or inversions, in the heliospheric magnetic field (HMF), which were accompanied by large spikes in solar wind speed. Both solar and in situ mechanisms have been suggested to explain the existence of HMF inversions in general. Previous work using Helios 1, covering 0.3-1AU,…
▽ More
In its first encounter at solar distances as close as r = 0.16AU, Parker Solar Probe (PSP) observed numerous local reversals, or inversions, in the heliospheric magnetic field (HMF), which were accompanied by large spikes in solar wind speed. Both solar and in situ mechanisms have been suggested to explain the existence of HMF inversions in general. Previous work using Helios 1, covering 0.3-1AU, observed inverted HMF to become more common with increasing r, suggesting that some heliospheric driving process creates or amplifies inversions. This study expands upon these findings, by analysing inversion-associated changes in plasma properties for the same large data set, facilitated by observations of 'strahl' electrons to identify the unperturbed magnetic polarity. We find that many inversions exhibit anti-correlated field and velocity perturbations, and are thus characteristically Alfvénic, but many also depart strongly from this relationship over an apparent continuum of properties. Inversions depart further from the 'ideal' Alfvénic case with increasing r, as more energy is partitioned in the field, rather than the plasma, component of the perturbation. This departure is greatest for inversions with larger density and magnetic field strength changes, and characteristic slow solar wind properties. We find no evidence that inversions which stray further from 'ideal' Alfvénicity have different generation processes from those which are more Alfvénic. Instead, different inversion properties could be imprinted based on transport or formation within different solar wind streams.
△ Less
Submitted 5 January, 2021;
originally announced January 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.
-
The Solar Corona during the Total Eclipse on 16 June 1806: Graphical Evidence of the Coronal Structure during the Dalton Minimum
Authors:
Hisashi Hayakawa,
Mathew J. Owens,
Michael Lockwood,
Mitsuru Sôma
Abstract:
Visible coronal structure, in particular the spatial evolution of coronal streamers, provides indirect information about solar magnetic activity and the underlying solar dynamo. Their apparent absence of structure observed during the total eclipses of throughout the Maunder Minimum has been interpreted as evidence of a significant change in the solar magnetic field from that during modern cycles.…
▽ More
Visible coronal structure, in particular the spatial evolution of coronal streamers, provides indirect information about solar magnetic activity and the underlying solar dynamo. Their apparent absence of structure observed during the total eclipses of throughout the Maunder Minimum has been interpreted as evidence of a significant change in the solar magnetic field from that during modern cycles. Eclipse observations available from the more recent Dalton Minimum may be able to provide further information, sunspot activity being between the levels seen during recent cycles and in the Maunder Minimum. Here, we show and examine two graphical records of the total solar eclipse on 1806 June 16, during the Dalton Minimum. These records show significant rays and streamers around an inner ring. The ring is estimated to be ~ 0.44 R_S in width and the streamers in excess of 11.88 R_S in length. In combination with records of spicules or prominences, these eclipse records visually contrast the Dalton Minimum with the Maunder Minimum in terms of their coronal structure and support the existing discussions based on the sunspot observations. These eclipse records are broadly consistent with the solar cycle phase in the modelled open solar flux and the reconstructed slow solar wind at most latitudes.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
The Evolution of Inverted Magnetic Fields Through the Inner Heliosphere
Authors:
Allan Macneil,
Mathew Owens,
Robert Wicks,
Mike Lockwood,
Sarah Bentley,
Mathew Lang
Abstract:
Local inversions are often observed in the heliospheric magnetic field (HMF), but their origins and evolution are not yet fully understood.Parker Solar Probe has recently observed rapid, Alfvenic, HMF inversions in the inner heliosphere, known as 'switchbacks', which have been interpreted as the possible remnants of coronal jets. It has also been suggested that inverted HMF may be produced by near…
▽ More
Local inversions are often observed in the heliospheric magnetic field (HMF), but their origins and evolution are not yet fully understood.Parker Solar Probe has recently observed rapid, Alfvenic, HMF inversions in the inner heliosphere, known as 'switchbacks', which have been interpreted as the possible remnants of coronal jets. It has also been suggested that inverted HMF may be produced by near-Sun interchange reconnection; a key process in mechanisms proposed for slow solar wind release. These cases suggest that the source of inverted HMF is near the Sun, and it follows that these inversions would gradually decay and straighten as they propagate out through the heliosphere. Alternatively, HMF inversions could form during solar wind transit, through phenomena such velocity shears, draping over ejecta, or waves and turbulence. Such processes are expected to lead to a qualitatively radial evolution of inverted HMF structures. Using Helios measurements spanning 0.3-1 AU, we examine the occurrence rate of inverted HMF, as well as other magnetic field morphologies, as a function of radial distance A, and find that it continually increases. This trend may be explained by inverted HMF observed between 0.3-1 AU being primarily driven by one or more of the above in-transit processes, rather than created at the Sun. We make suggestions as to the relative importance of these different processes based on the evolution of the magnetic field properties associated with inverted HMF. We also explore alternative explanations outside of our suggested driving processes which may lead to the observed trend.
△ Less
Submitted 11 April, 2020;
originally announced April 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.
-
A New Calibrated Sunspot Group Series Since 1749: Statistics of Active Day Fractions
Authors:
I. G. Usoskin,
G. A. Kovaltsov,
M. Lockwood,
K. Mursula,
M. Owens,
S. K. Solanki
Abstract:
Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a nov…
▽ More
Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a novel method to calibrate the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900-1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds [S_S] defined such that the observer is assumed to miss all of the groups with an area smaller than S_S and report all the groups larger than S_S. Next, using a Monte-Carlo method we construct, from the reference data set, a correction matrix for each observer. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot group records for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series displays secular minima around 1800 (Dalton minimum) and 1900 (Gleissberg minimum), as well as the Modern grand maximum of activity in the second half of the 20th century. The uniqueness of the grand maximum is confirmed for the last 250 years. It is shown that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the 18th and 19th centuries in some reconstructions.
△ Less
Submitted 20 December, 2015;
originally announced December 2015.
-
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.
-
The Maunder minimum (1645--1715) was indeed a Grand minimum: A reassessment of multiple datasets
Authors:
Ilya G. Usoskin,
Rainer Arlt,
Eleanna Asvestari,
Ed Hawkins,
Maarit Käpylä,
Gennady A. Kovaltsov,
Natalie Krivova,
Michael Lockwood,
Kalevi Mursula,
Jezebel O'Reilly,
Matthew Owens,
Chris J. Scott,
Dmitry D. Sokoloff,
Sami K. Solanki,
Willie Soon,
José M. Vaquero
Abstract:
Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity durin…
▽ More
Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum.
Methods: We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum.
Results: The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets.
Conclusions: We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level.
△ Less
Submitted 18 July, 2015;
originally announced July 2015.
-
Further considerations of cosmic ray modulation of infra-red radiation in the atmosphere
Authors:
Karen Aplin,
Mike Lockwood
Abstract:
Understanding effects of ionisation in the lower atmosphere is a new interdisciplinary area, crossing traditionally distinct scientific boundaries. Following the paper of Erlykin et al. (Astropart. Phys. 57--58 (2014) 26--29) we develop the interpretation of observed changes in long-wave (LW) radiation (Aplin and Lockwood, Env. Res. Letts. 8, 015026 (2013)), by taking account of cosmic ray ionisat…
▽ More
Understanding effects of ionisation in the lower atmosphere is a new interdisciplinary area, crossing traditionally distinct scientific boundaries. Following the paper of Erlykin et al. (Astropart. Phys. 57--58 (2014) 26--29) we develop the interpretation of observed changes in long-wave (LW) radiation (Aplin and Lockwood, Env. Res. Letts. 8, 015026 (2013)), by taking account of cosmic ray ionisation yields and atmospheric radiative transfer. To demonstrate this, we show that the thermal structure of the whole atmosphere needs to be considered along with the vertical profile of ionisation. Allowing for ionisation by all components of a cosmic ray shower and not just by the muons, reveals that the effect we have detected is certainly not inconsistent with laboratory observations of the LW absorption cross section. The analysis presented here, although very different from that of Erlykin et al., does come to the same conclusion that the events detected were not caused by individual cosmic ray primaries -- not because it is impossible on energetic grounds, but because events of the required energy are too infrequent for the 12/hr rate at which they were seen by the AL experiment. The present paper numerically models the effect of three different scenario changes to the primary GCR spectrum which all reproduce the required magnitude of the observed effect. However, they cannot solely explain the observed delay in the peak effect which, if confirmed, would appear to open up a whole new and interesting area in the study of water oligomers and their effects on LW radiation. We argue that a technical artefact in the earlier experiment is highly unlikely and that our initial observations merit both follow-up experiments and more rigorous, self-consistent, three-dimensional radiative transfer modelling.
△ Less
Submitted 5 March, 2015;
originally announced March 2015.
-
The 22-Year Hale Cycle in Cosmic Ray Flux - Evidence for Direct Heliospheric Modulation
Authors:
Simon R Thomas,
Mathew J Owens,
Mike Lockwood
Abstract:
The ability to predict times of greater galactic cosmic ray (GCR) fluxes is important for reducing the hazards caused by these particles to satellite communications, aviation, or astronauts. The 11-year solar cycle variation in cosmic rays is highly correlated with the strength of the heliospheric magnetic field. Differences in GCR flux during alternate solar cycles yield a 22-year cycle, known as…
▽ More
The ability to predict times of greater galactic cosmic ray (GCR) fluxes is important for reducing the hazards caused by these particles to satellite communications, aviation, or astronauts. The 11-year solar cycle variation in cosmic rays is highly correlated with the strength of the heliospheric magnetic field. Differences in GCR flux during alternate solar cycles yield a 22-year cycle, known as the Hale Cycle, which is thought to be due to different particle drift patterns when the northern solar pole has predominantly positive (denoted a qA>0 cycle) or negative (qA<0) polarities. This results in the onset of the peak cosmic ray flux at Earth occurring earlier during qA>0 cycles than for qA<0 cycles and hence the peak being more domed for qA>0 and more sharply peaked for qA<0. In this study, we demonstrate that properties of the large-scale heliospheric magnetic field are different during the declining phase of the qA<0 and qA>0 solar cycles, when the difference in GCR flux is most apparent. This suggests that particle drifts may not be the sole mechanism responsible for the Hale Cycle in GCR flux at Earth. However, it is also demonstrated that these polarity-dependent heliospheric differences are evident during the space-age but much less clear in earlier data: using geomagnetic reconstructions, it is shown that for the period of 1905 - 1965, alternate polarities do not give as significant a difference during the declining phase of the solar cycle. Thus we suggest that the 22-year cycle in cosmic ray flux is at least partly the result of direct modulation by the heliospheric magnetic field and that this effect may be primarily limited to the grand solar maximum of the space-age.
△ Less
Submitted 27 July, 2013; v1 submitted 18 April, 2013;
originally announced April 2013.