Showing 1–33 of 33 results for author: Owens, M J

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.

Comments: 31 pages, 10 figures

On the Origin of the sudden Heliospheric Open Magnetic Flux Enhancement during the 2014 Pole Reversal

Authors: Stephan G. Heinemann, Mathew J. Owens, Manuela Temmer, James A. Turtle, Charles N. Arge, Carl J. Henney, Jens Pomoell, Eleanna Asvestari, Jon A. Linker, Cooper Downs, Ronald M. Caplan, Stefan J. Hofmeister, Camilla Scolini, Rui F. Pinto, Maria S. Madjarska

Abstract: Coronal holes are recognized as the primary sources of heliospheric open magnetic flux (OMF). However, a noticeable gap exists between in-situ measured OMF and that derived from remote sensing observations of the Sun. In this study, we investigate the OMF evolution and its connection to solar structures throughout 2014, with special emphasis on the period from September to October, where a sudden… ▽ More Coronal holes are recognized as the primary sources of heliospheric open magnetic flux (OMF). However, a noticeable gap exists between in-situ measured OMF and that derived from remote sensing observations of the Sun. In this study, we investigate the OMF evolution and its connection to solar structures throughout 2014, with special emphasis on the period from September to October, where a sudden and significant OMF increase was reported. By deriving the OMF evolution at 1au, modeling it at the source surface, and analyzing solar photospheric data, we provide a comprehensive analysis of the observed phenomenon. First, we establish a strong correlation between the OMF increase and the solar magnetic field derived from a Potential Field Source Surface (PFSS) model ($cc_{\mathrm{Pearson}}=0.94$). Moreover, we find a good correlation between the OMF and the open flux derived from solar coronal holes ($cc_{\mathrm{Pearson}}=0.88$), although the coronal holes only contain $14-32\%$ of the Sun's total open flux. However, we note that while the OMF evolution correlates with coronal hole open flux, there is no correlation with the coronal hole area evolution ($cc_{\mathrm{Pearson}}=0.0$). The temporal increase in OMF correlates with the vanishing remnant magnetic field at the southern pole, caused by poleward flux circulations from the decay of numerous active regions months earlier. Additionally, our analysis suggests a potential link between the OMF enhancement and the concurrent emergence of the largest active region in solar cycle 24. In conclusion, our study provides insights into the strong increase in OMF observed during September to October 2014. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: accepted in ApJ

Tracking solar radio bursts using Bayesian multilateration

Authors: L. A. Cañizares, S. T. Badman, S. A. Maloney, M. J. Owens, D. M. Weigt, E. P. Carley, P. T. Gallagher

Abstract: Solar radio bursts (SRBs), are emitted by electrons propagating through the corona and interplanetary space. Tracking such bursts is key to understanding the properties of accelerated electrons and radio wave propagation as well as the local plasma environment that they propagate through. Here, we present a novel multilateration algorithm called BayEsian LocaLisation Algorithm (BELLA). In addition… ▽ More Solar radio bursts (SRBs), are emitted by electrons propagating through the corona and interplanetary space. Tracking such bursts is key to understanding the properties of accelerated electrons and radio wave propagation as well as the local plasma environment that they propagate through. Here, we present a novel multilateration algorithm called BayEsian LocaLisation Algorithm (BELLA). In addition, apparent SRB positions from BELLA are compared with comparable localisation methods and the predictions of solar wind models. BELLA uses Bayesian inference to create probabilistic distributions of source positions and their uncertainties. This facilitates the estimation of algorithmic, instrumental, and physical uncertainties in a quantitative manner. We validated BELLA using simulations and a Type III SRB observed by STEREO A/B and Wind. BELLA tracked the Type III source from $\sim$ 10--150 $R_{sun}$ (2-0.15 MHz) along a spiral trajectory. This allowed for an estimate of an apparent solar wind speed of $v_{sw} \sim$ 400 km s$^{-1}$ and a source longitude of $φ_0 \sim$ 30deg. We compared these results with well-established methods of positioning: Goniopolarimetric (GP), analytical time-difference-of-arrival (TDOA), and Solar radio burst Electron Motion Tracker (SEMP). We found them to be in agreement with the results obtained by BELLA. Additionally, the results aligned with solar wind properties assimilated by the Heliospheric Upwind Extrapolation with time dependence (HUXt) model. We have validated BELLA and used it to identify apparent source positions as well as velocities and densities of the solar wind. Furthermore, we identified higher than expected electron densities, suggesting that the true emission sources were at lower altitudes than those identified by BELLA, an effect that may be due to appreciable scattering of electromagnetic waves by electrons in interplanetary space. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 13 pages, 12 figures, in press, accepted by Astronomy and Astrophysics for publication

Coronal Models and Detection of Open Magnetic Field

Authors: Eleanna Asvestari, Manuela Temmer, Ronald M. Caplan, Jon A. Linker, Stephan G. Heinemann, Rui F. Pinto, Carl J. Henney, Charles N. Arge, Mathew J. Owens, Maria S. Madjarska, Jens Pomoell, Stefan J. Hofmeister, Camilla Scolini, Evangelia Samara

Abstract: A plethora of coronal models, from empirical to more complex magnetohydrodynamic (MHD) ones, are being used for reconstructing the coronal magnetic field topology and estimating the open magnetic flux. However, no individual solution fully agrees with coronal hole observations and in situ measurements of open flux at 1~AU, as there is a strong deficit between model and observations contributing to… ▽ More A plethora of coronal models, from empirical to more complex magnetohydrodynamic (MHD) ones, are being used for reconstructing the coronal magnetic field topology and estimating the open magnetic flux. However, no individual solution fully agrees with coronal hole observations and in situ measurements of open flux at 1~AU, as there is a strong deficit between model and observations contributing to the known problem of the missing open flux. In this paper we investigate the possible origin of the discrepancy between modeled and observed magnetic field topology by assessing the effect on the simulation output by the choice of the input boundary conditions and the simulation set up, including the choice of numerical schemes and the parameter initialization. In the frame of this work, we considered four potential field source surface based models and one fully MHD model, different types of global magnetic field maps and model initiation parameters. After assessing the model outputs using a variety of metrics, we conclude that they are highly comparable regardless of the differences set at initiation. When comparing all models to coronal hole boundaries extracted by extreme ultraviolet (EUV) filtergrams we find that they do not compare well. This miss-match between observed and modeled regions of open field is a candidate contributing to the open flux problem. △ Less

Submitted 7 November, 2023; originally announced November 2023.

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.

Comments: Accepted for publication in Advances in Space Research

Modelling Cosmic Radiation Events in the Tree-ring Radiocarbon Record

Authors: Qingyuan Zhang, Utkarsh Sharma, Jordan A. Dennis, Andrea Scifo, Margot Kuitems, Ulf Buentgen, Mathew J. Owens, Michael W. Dee, Benjamin J. S. Pope

Abstract: Annually-resolved measurements of the radiocarbon content in tree-rings have revealed rare sharp rises in carbon-14 production. These 'Miyake events' are likely produced by rare increases in cosmic radiation from the Sun or other energetic astrophysical sources. The radiocarbon produced is not only circulated through the Earth's atmosphere and oceans, but also absorbed by the biosphere and locked… ▽ More Annually-resolved measurements of the radiocarbon content in tree-rings have revealed rare sharp rises in carbon-14 production. These 'Miyake events' are likely produced by rare increases in cosmic radiation from the Sun or other energetic astrophysical sources. The radiocarbon produced is not only circulated through the Earth's atmosphere and oceans, but also absorbed by the biosphere and locked in the annual growth rings of trees. To interpret high-resolution tree-ring radiocarbon measurements therefore necessitates modelling the entire global carbon cycle. Here, we introduce 'ticktack', the first open-source Python package that connects box models of the carbon cycle with modern Bayesian inference tools. We use this to analyse all public annual 14C tree data, and infer posterior parameters for all six known Miyake events. They do not show a consistent relationship to the solar cycle, and several display extended durations that challenge either astrophysical or geophysical models. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Comments: Accepted Proceedings of the Royal Society A. 19 pages 6 figures body, 12 pages appendices which are supplementary material in the published version

Coronal Hole Detection and Open Magnetic Flux

Authors: J. A. Linker, S. G. Heinemann, M. Temmer, M. J. Owens, R. M. Caplan, C. N. Arge, E. Asvestari, V. Delouille, C. Downs, S. J. Hofmeister, I. C. Jebaraj, M. Madjarska, R. Pinto, J. Pomoell, E. Samara, C. Scolini, B. Vrsnak

Abstract: Many scientists use coronal hole (CH) detections to infer open magnetic flux. Detection techniques differ in the areas that they assign as open, and may obtain different values for the open magnetic flux. We characterize the uncertainties of these methods, by applying six different detection methods to deduce the area and open flux of a near-disk center CH observed on 9/19/2010, and applying a sin… ▽ More Many scientists use coronal hole (CH) detections to infer open magnetic flux. Detection techniques differ in the areas that they assign as open, and may obtain different values for the open magnetic flux. We characterize the uncertainties of these methods, by applying six different detection methods to deduce the area and open flux of a near-disk center CH observed on 9/19/2010, and applying a single method to five different EUV filtergrams for this CH. Open flux was calculated using five different magnetic maps. The standard deviation (interpreted as the uncertainty) in the open flux estimate for this CH was about 26%. However, including the variability of different magnetic data sources, this uncertainty almost doubles to 45%. We use two of the methods to characterize the area and open flux for all CHs in this time period. We find that the open flux is greatly underestimated compared to values inferred from in-situ measurements (by 2.2-4 times). We also test our detection techniques on simulated emission images from a thermodynamic MHD model of the solar corona. We find that the methods overestimate the area and open flux in the simulated CH, but the average error in the flux is only about 7%. The full-Sun detections on the simulated corona underestimate the model open flux, but by factors well below what is needed to account for the missing flux in the observations. Under-detection of open flux in coronal holes likely contributes to the recognized deficit in solar open flux, but is unlikely to resolve it. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 28 pages, 10 figures, submitted to ApJ

Multi-spacecraft Study of the Solar Wind at Solar Minimum: Dependence on Latitude and Transient Outflows

Authors: R. Laker, T. S. Horbury, S. D. Bale, L. Matteini, T. Woolley, L. D. Woodham, J. E. Stawarz, E. E. Davies, J. P. Eastwood, M. J. Owens, H. O'Brien, V. Evans, V. Angelini, I. Richter, D. Heyner, C. J. Owen, P. Louarn, A. Federov

Abstract: The recent launches of Parker Solar Probe (PSP), Solar Orbiter (SO) and BepiColombo, along with several older spacecraft, have provided the opportunity to study the solar wind at multiple latitudes and distances from the Sun simultaneously. We take advantage of this unique spacecraft constellation, along with low solar activity across two solar rotations between May and July 2020, to investigate h… ▽ More The recent launches of Parker Solar Probe (PSP), Solar Orbiter (SO) and BepiColombo, along with several older spacecraft, have provided the opportunity to study the solar wind at multiple latitudes and distances from the Sun simultaneously. We take advantage of this unique spacecraft constellation, along with low solar activity across two solar rotations between May and July 2020, to investigate how the solar wind structure, including the Heliospheric Current Sheet (HCS), varies with latitude. We visualise the sector structure of the inner heliosphere by ballistically mapping the polarity and solar wind speed from several spacecraft onto the Sun's source surface. We then assess the HCS morphology and orientation with the in situ data and compare with a predicted HCS shape. We resolve ripples in the HCS on scales of a few degrees in longitude and latitude, finding that the local orientation of sector boundaries were broadly consistent with the shape of the HCS but were steepened with respect to a modelled HCS at the Sun. We investigate how several CIRs varied with latitude, finding evidence for the compression region affecting slow solar wind outside the latitude extent of the faster stream. We also identified several transient structures associated with HCS crossings, and speculate that one such transient may have disrupted the local HCS orientation up to five days after its passage. We have shown that the solar wind structure varies significantly with latitude, with this constellation providing context for solar wind measurements that would not be possible with a single spacecraft. These measurements provide an accurate representation of the solar wind within $\pm 10^{\circ}$ latitude, which could be used as a more rigorous constraint on solar wind models and space weather predictions. In the future, this range of latitudes will increase as SO's orbit becomes more inclined. △ Less

Submitted 22 June, 2021; v1 submitted 27 February, 2021; originally announced March 2021.

Comments: Accepted version

Journal ref: A&A 652, A105 (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.

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.

Comments: Published in Monthly Notices of the Royal Astronomical Society

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.

Comments: Accepted for publication in J. Space Weather Space Clim., contains 23 Figures plus has an Appendix containing 4 Figures

In situ multi-spacecraft and remote imaging observations of the first CME detected by Solar Orbiter and BepiColombo

Authors: E. E. Davies, C. Möstl, M. J. Owens, A. J. Weiss, T. Amerstorfer, J. Hinterreiter, M. Bauer, R. L. Bailey, M. A. Reiss, R. J. Forsyth, T. S. Horbury, H. O'Brien, V. Evans, V. Angelini, D. Heyner, I. Richter, H-U. Auster, W. Magnes, W. Baumjohann, D. Fischer, D. Barnes, J. A. Davies, R. A. Harrison

Abstract: On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ i… ▽ More On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ in longitude. The source of the CME, which was launched on April 15th, was an almost entirely isolated streamer blowout. STEREO-A observed the event remotely from -75.1$^\circ$ longitude, which is an exceptionally well suited viewpoint for heliospheric imaging of an Earth directed CME. The configuration of the four spacecraft has provided an exceptionally clean link between remote imaging and in situ observations of the CME. We have used the in situ observations of the CME at Solar Orbiter, Wind, and BepiColombo, and the remote observations of the CME at STEREO-A in combination with flux rope models to determine the global shape of the CME and its evolution as it propagated through the inner heliosphere. A clear flattening of the CME cross-section has been observed by STEREO-A, and further confirmed by comparing profiles of the flux rope models to the in situ data, where the distorted flux rope cross-section qualitatively agrees most with in situ observations of the magnetic field at Solar Orbiter. Comparing in situ observations of the magnetic field between spacecraft, we find that the dependence of the maximum (mean) magnetic field strength decreases with heliocentric distance as $r^{-1.24 \pm 0.50}$ ($r^{-1.12 \pm 0.14}$), in disagreement with previous studies. Further assessment of the axial and poloidal magnetic field strength dependencies suggests that the expansion of the CME is likely neither self-similar nor cylindrically symmetric. △ Less

Submitted 23 February, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Journal ref: A&A 656, A2 (2021)

Parker Solar Probe Observations of Suprathermal Electron Flux Enhancements Originating from Coronal Hole Boundaries

Authors: Allan R Macneil, Mathew J Owens, Laura Berčič, Adam J Finley

Abstract: Reconnection between pairs of solar magnetic flux elements, one open and the other a closed loop, is theorised to be a crucial process for both maintaining the structure of the corona and producing the solar wind. This 'interchange reconnection' is expected to be particularly active at the open-closed boundaries of coronal holes (CHs). Previous analysis of solar wind data at 1AU indicated that pea… ▽ More Reconnection between pairs of solar magnetic flux elements, one open and the other a closed loop, is theorised to be a crucial process for both maintaining the structure of the corona and producing the solar wind. This 'interchange reconnection' is expected to be particularly active at the open-closed boundaries of coronal holes (CHs). Previous analysis of solar wind data at 1AU indicated that peaks in the flux of suprathermal electrons at slow-fast stream interfaces may arise from magnetic connection to the CH boundary, rather than dynamic effects such as compression. Further, offsets between the peak and stream interface locations are suggested to be the result of interchange reconnection at the source. As a preliminary test of these suggestions, we analyse two solar wind streams observed during the first Parker Solar Probe (PSP) perihelion encounter, each associated with equatorial CH boundaries (one leading and one trailing with respect to rotation). Each stream features a peak in suprathermal electron flux, the locations and associated plasma properties of which are indicative of a solar origin, in agreement with previous suggestions from 1AU observations. Discrepancies between locations of the flux peaks and other features suggest these peaks may too be shifted by source region interchange reconnection. Our interpretation of each event is compatible with a global pattern of open flux transport, although random footpoint motions or other explanations remain feasible. These exploratory results highlight future opportunities for statistical studies regarding interchange reconnection and flux transport at CH boundaries with modern near-Sun missions. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Using gradient boosting regression to improve ambient solar wind model predictions

Authors: R. L. Bailey, M. A. Reiss, C. N. Arge, C. Möstl, M. J. Owens, U. V. Amerstorfer, C. J. Henney, T. Amerstorfer, A. J. Weiss, J. Hinterreiter

Abstract: Studying the ambient solar wind, a continuous pressure-driven plasma flow emanating from our Sun, is an important component of space weather research. The ambient solar wind flows in interplanetary space determine how solar storms evolve through the heliosphere before reaching Earth, and especially during solar minimum are themselves a driver of activity in the Earth's magnetic field. Accurately f… ▽ More Studying the ambient solar wind, a continuous pressure-driven plasma flow emanating from our Sun, is an important component of space weather research. The ambient solar wind flows in interplanetary space determine how solar storms evolve through the heliosphere before reaching Earth, and especially during solar minimum are themselves a driver of activity in the Earth's magnetic field. Accurately forecasting the ambient solar wind flow is therefore imperative to space weather awareness. Here we present a machine learning approach in which solutions from magnetic models of the solar corona are used to output the solar wind conditions near the Earth. The results are compared to observations and existing models in a comprehensive validation analysis, and the new model outperforms existing models in almost all measures. In addition, this approach offers a new perspective to discuss the role of different input data to ambient solar wind modeling, and what this tells us about the underlying physical processes. The final model discussed here represents an extremely fast, well-validated and open-source approach to the forecasting of ambient solar wind at Earth. △ Less

Submitted 23 March, 2021; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: 15 pages, 6 figures, 3 tables

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.

Comments: main text 14 pages, references 6 pages, 3 figures, and accepted for publication in the ApJ

Journal ref: The Astrophysical Journal, 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.

Journal ref: J. Space Weather Space Climate, 8, A58, 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.

Journal ref: J. Space Weather & Space Climate, 8, A53, 2018

The State of the Solar Wind, Magnetosphere, and Ionosphere During the Maunder Minimum

Authors: Pete Riley, Roberto Lionello, Jon A. Linker, Mathew J. Owens

Abstract: Both direct observations and reconstructions from various datasets, suggest that conditions were radically different during the Maunder Minimum (MM) than during the space era. Using an MHD model, we develop a set of feasible solutions to infer the properties of the solar wind during this interval. Additionally, we use these results to drive a global magnetospheric model. Finally, using the 2008/20… ▽ More Both direct observations and reconstructions from various datasets, suggest that conditions were radically different during the Maunder Minimum (MM) than during the space era. Using an MHD model, we develop a set of feasible solutions to infer the properties of the solar wind during this interval. Additionally, we use these results to drive a global magnetospheric model. Finally, using the 2008/2009 solar minimum as an upper limit for MM conditions, we use results from the International Reference Ionosphere (ILI) model to speculate on the state of the ionosphere. The results describe interplanetary, magnetospheric, and ionospheric conditions that were substantially different than today. For example: (1) the solar wind density and magnetic field strength were an order of magnitude lower; (2) the Earth's magnetopause and shock standoff distances were a factor of two larger; and (3) the maximum electron density in the ionosphere was substantially lower. △ Less

Submitted 19 March, 2018; originally announced March 2018.

Comments: 4 pages, 1 figure, IAU340 Conference

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.

Journal ref: J. Space Weather Space Clim., doi 10.1051/swsc/2017048, 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.

Comments: In press, accepted for publication by J. Space Weather Space Clim. Ms# swsc170035 (accepted version, 28 pages, 2 Tables, 12 Figures)

The Open Flux Problem

Authors: J. A. Linker, R. M. Caplan, C. Downs, P Riley, Z Mikic, R. Lionello, C. J. Henney, C. N. Arge, Y. Liu, M. L. Derosa, A. Yeates, M. J. Owens

Abstract: The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field reg… ▽ More The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field regions in the model should approximately correspond to coronal holes observed in emission, and (2) the magnitude of the open magnetic flux in the model should match that inferred from in situ spacecraft measurements. In this study, we calculate both MHD and PFSS solutions using fourteen different magnetic maps produced from five different types of observatory magnetograms, for the time period surrounding July, 2010. We have found that for all of the model/map combinations, models that have coronal hole areas close to observations underestimate the interplanetary magnetic flux, or, conversely, for models to match the interplanetary flux, the modeled open field regions are larger than coronal holes observed in EUV emission. In an alternative approach, we estimate the open magnetic flux entirely from solar observations by combining automatically detected coronal holes for Carrington rotation 2098 with observatory synoptic magnetic maps. This approach also underestimates the interplanetary magnetic flux. Our results imply that either typical observatory maps underestimate the Sun's magnetic flux, or a significant portion of the open magnetic flux is not rooted in regions that are obviously dark in EUV and X-ray emission. △ Less

Submitted 5 September, 2017; v1 submitted 7 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.

Comments: In press with Solar Physics

Report number: Ms. # SOLA-D-15-00228R1

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.

Comments: Solar Physics, in Press. 26 pages, 8 figures

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.

Comments: In press, Astrophys. J., 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.

Comments: 35 pages, 11 figures

Report number: Ms #: SOLA-D-15-00189R3

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.

Comments: 14 pages, 1 Table, 11 figures

Using Coordinated Observations in Polarised White Light and Faraday Rotation to Probe the Spatial Position and Magnetic Field of an Interplanetary Sheath

Authors: Ming Xiong, Jackie A. Davies, Xueshang Feng, Mathew J. Owens, Richard A. Harrison, Chris J. Davis, Ying D. Liu

Abstract: Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White-light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons, and therefore depend on both the viewing geometry and the electron density. The Faraday rotation (FR) of radio wa… ▽ More Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White-light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons, and therefore depend on both the viewing geometry and the electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component $B_\parallel$, and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modelling of an Earth-directed shock and synthesise the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance $I$ decreases with heliocentric distance $r$ roughly according to the expression $I \propto r^{-3}$. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site $\mathbf{r}_{\rm sheath}$ and the mass of plasma at that position $M_{\rm sheath}$ can be inferred from the polarisation of the shock-associated enhancement in WL radiance. From the FR measurements, the local $B_{\parallel {\rm sheath}}$ at $\mathbf{r}_{\rm sheath}$ can then be estimated. Simultaneous observations in polarised WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties. △ Less

Submitted 15 August, 2013; originally announced August 2013.

Comments: 20 pages, 11 figures, 1 table. To be published by Astrophys. J., 2013

Journal ref: Astrophysical Journal, 777:32 (14pp), 2013

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.

Comments: 16 pages, 7 figures

Tracking the momentum flux of a CME and quantifying its influence on geomagnetically induced currents at Earth

Authors: Neel P. Savani, A. Vourlidas, A. Pulkkinen, T. Nieves-Chinchilla, B. Lavraud, M. J. Owens

Abstract: We investigate a CME propagating towards Earth on 29 March 2011. This event is specifically chosen for its predominately northward directed magnetic field, so that the influence from the momentum flux onto Earth can be isolated. We focus our study on understanding how a small Earth-directed segment propagates. Mass images are created from the white-light cameras onboard STEREO which are also conve… ▽ More We investigate a CME propagating towards Earth on 29 March 2011. This event is specifically chosen for its predominately northward directed magnetic field, so that the influence from the momentum flux onto Earth can be isolated. We focus our study on understanding how a small Earth-directed segment propagates. Mass images are created from the white-light cameras onboard STEREO which are also converted into mass height-time maps (mass J-maps). The mass tracks on these J-maps correspond to the sheath region between the CME and its associated shock front as detected by in situ measurements at L1. A time-series of mass measurements from the STEREO COR-2A instrument are made along the Earth propagation direction. Qualitatively, this mass time-series shows a remarkable resemblance to the L1 in situ density series. The in situ measurements are used as inputs into a 3D magnetospheric space weather simulation from CCMC. These simulations display a sudden compression of the magnetosphere from the large momentum flux at the leading edge of the CME and predictions are made for the time-derivative of the magnetic field (dB/dt) on the ground. The predicted dB/dt were then compared with observations from specific equatorially-located ground stations and show notable similarity. This study of the momentum of a CME from the Sun down to its influence on magnetic ground stations on Earth is presented as preliminary proof of concept, such that future attempts may try to use remote sensing to create density and velocity time-series as inputs to magnetospheric simulations. △ Less

Submitted 8 March, 2013; originally announced March 2013.

Comments: Accepted for publication 8th March 2013. Submitted 18th Dec 2012. 19 Pages, 10 figures, 2 Appendices

Journal ref: Space Weather 2013

Effects of Thomson-Scattering Geometry on White-Light Imaging of an Interplanetary Shock: Synthetic Observations from Forward Magnetohydrodynamic Modelling

Authors: Ming Xiong, J. A. Davies, M. M. Bisi, M. J. Owens, R. A. Fallows, G. D. Dorrian

Abstract: Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using for… ▽ More Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time-elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east-west asymmetry in its brightness and degree of polarization. Therefore, measurement of the interplanetary polarized brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging. △ Less

Submitted 2 May, 2012; originally announced May 2012.

Comments: 27 pages, 8 figures, 4 tables, Solar Physics

Journal ref: Solar Physics, 285: 369-389, 2013

Forward modelling to determine the observational signatures of white-light imaging and interplanetary scintillation for the propagation of an interplanetary shock in the ecliptic plane

Authors: Ming Xiong, A. R. Breen, M. M. Bisi, M. J. Owens, R. A. Fallows, G. D. Dorrian, J. A. Davies, P. Thomasson

Abstract: Recent coordinated observations of interplanetary scintillation (IPS) and stereoscopic heliospheric imagers (HIs) are significant to continuously track the propagation and evolution of solar eruptions throughout interplanetary space. In order to obtain a better understanding of the observational signatures in these two remote-sensing techniques, the magnetohydrodynamics of the macro-scale interpla… ▽ More Recent coordinated observations of interplanetary scintillation (IPS) and stereoscopic heliospheric imagers (HIs) are significant to continuously track the propagation and evolution of solar eruptions throughout interplanetary space. In order to obtain a better understanding of the observational signatures in these two remote-sensing techniques, the magnetohydrodynamics of the macro-scale interplanetary disturbance and the radio-wave scattering of the micro-scale electron-density fluctuation are coupled and investigated using a newly-constructed multi-scale numerical model. This model is then applied to a case of an interplanetary shock propagation within the ecliptic plane. The shock could be nearly invisible to an HI, once entering the Thomson-scattering sphere of the HI. The asymmetry in the optical images between the western and eastern HIs suggests the shock propagation off the Sun-Earth line. Meanwhile, an IPS signal, strongly dependent on the local electron density, is insensitive to the density cavity far downstream of the shock front. When this cavity (or the shock nose) is cut through by an IPS ray-path, a single speed component at the flank (or the nose) of the shock can be recorded; when an IPS ray-path penetrates the sheath between the shock nose and this cavity, two speed components at the sheath and flank can be detected. Moreover, once a shock front touches an IPS ray-path, the derived position and speed at the irregularity source of this IPS signal, together with an assumption of a radial and constant propagation of the shock, can be used to estimate the later appearance of the shock front in the elongation of the HI field of view. The results of synthetic measurements from forward modelling are helpful in inferring the in-situ properties of coronal mass ejection from real observational data via an inverse approach. △ Less

Submitted 13 August, 2010; originally announced August 2010.

Comments: 31 pages, 8 figures, Journal of Atmospheric and Solar-Terrestrial Physics, accepted

Journal ref: Journal of Atmospheric and Solar-Terrestrial Physics, 73(2011), 1270-1280

Numerical Simulation of the May 12, 1997 CME Event - the Role of Magnetic Reconnection

Authors: O. Cohen, G. D. R. Attrill, N. A. Schwadron, N. U. Crooker, M. J. Owens, C. Downs, T. I. Gombosi

Abstract: We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the May 12, 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux… ▽ More We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the May 12, 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance ($\approx 1R_\odot$) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure. △ Less

Submitted 22 June, 2010; originally announced June 2010.

Comments: Accepted to JGR Space Physics, 15 pages, 8 figures

Impact of coronal mass ejections, interchange reconnection, and disconnection on heliospheric magnetic field strength

Authors: N. U. Crooker, M. J. Owens

Abstract: An update of Owens et al. (2008) shows that the relationship between the coronal mass ejection (CME) rate and the heliospheric magnetic field strength predicts a field floor of less than 4 nT at 1 AU. This implies that the record low values measured during this solar minimum do not necessarily contradict the idea that open flux is conserved. The results are consistent with the hypothesis that CMEs… ▽ More An update of Owens et al. (2008) shows that the relationship between the coronal mass ejection (CME) rate and the heliospheric magnetic field strength predicts a field floor of less than 4 nT at 1 AU. This implies that the record low values measured during this solar minimum do not necessarily contradict the idea that open flux is conserved. The results are consistent with the hypothesis that CMEs add flux to the heliosphere and interchange reconnection between open flux and closed CME loops subtracts flux. An existing model embracing this hypothesis, however, overestimates flux during the current minimum, even though the CME rate has been low. The discrepancy calls for reasonable changes in model assumptions. △ Less

Submitted 13 April, 2010; originally announced April 2010.

Comments: in press, ASP Conference Series, SOHO-23: Understanding a Peculiar Solar Minimum, edited by S. Cranmer, T. Hoeksema, and J. Kohl