Predicting space climate change

The recent decline in the open magnetic flux of the Sun heralds the end of the Grand Solar Maximum (GSM) that has persisted throughout the space age, during which the largest‐fluence Solar Energetic Particle (SEP) events have been rare and Galactic Cosmic Ray (GCR) fluxes have been relatively low. In the absence of a predictive model of the solar dynamo, we here make analogue forecasts by studying past variations of solar activity in order to evaluate how long‐term change in space climate may influence the hazardous energetic particle environment of the Earth in the future. We predict the probable future variations in GCR flux, near‐Earth interplanetary magnetic field (IMF), sunspot number, and the probability of large SEP events, all deduced from cosmogenic isotope abundance changes following 24 GSMs in a 9300‐year record.

Solar cycle 24: implications for energetic particles and long-term space climate change

The recent solar minimum was the longest and deepest of the space age, with the lowest average sunspot numbers for nearly a century. The Sun appears to be exiting a grand solar maximum (GSM) of activity which has persisted throughout the space age, and is headed into a significantly quieter period. Indeed, initial observations of solar cycle 24 (SC24) continue to show a relatively low heliospheric magnetic field strength and sunspot number (R), despite the average latitude of sunspots and the inclination of the heliospheric current sheet showing the rise to solar maximum is well underway. We extrapolate the available SC24 observations forward in time by assuming R will continue to follow a similar form to previous cycles, despite the end of the GSM, and predict a very weak cycle 24, with R peaking at ∼65–75 around the middle/end of 2012. Similarly, we estimate the heliospheric magnetic field strength will peak around 6nT. We estimate that average galactic cosmic ray fluxes above 1GV rigidity will be ∼10% higher in SC24 than SC23 and that the probability of a large SEP event during this cycle is 0.8, compared to 0.5 for SC23. Comparison of the SC24 R estimates with previous ends of GSMs inferred from 9300 years of cosmogenic isotope data places the current evolution of the Sun and heliosphere in the lowest 5% of cases, suggesting Maunder Minimum conditions are likely within the next 40 years.

Memory effects in liquid crystal elastomers

Free-standing monodomain liquid crystal elastomer samples are shown to have a complete memory of the orientational configuration at the time of cross-linking. This memory is demonstrated through samples in which the parent polymer system is first aligned in a magnetic field prior to cross-linking. These films show reversible nematic-isotropic phase transitions and x-ray scattering patterns characteristic of nematic phases. The liquid crystal elastomer films exhibit a remarkable memory effect, in that the sample may be held at temperatures well above the nematic-isotropic transition for extended periods ( > 2 weeks), but on cooling into the liquid crystal phase region, both the original director alignment and the degree of preferred orientation are recovered. It is demonstrated that these novel memory effects are equilibrium in nature. The origins of this phenomena in terms of coupling between the mesogenic side-chains and the polymer network are discussed.

Inner plasma structure of the low-latitude reconnection layer

We report a clear transition through a reconnection layer at the low-latitude magnetopause which shows a complete traversal across all reconnected field lines during northwestward interplanetary magnetic field (IMF) conditions. The associated plasma populations confirm details of the electron and ion mixing and the time history and acceleration through the current layer. This case has low magnetic shear with a strong guide field and the reconnection layer contains a single density depletion layer on the magnetosheath side which we suggest results from nearly field-aligned magnetosheath flows. Within the reconnection boundary layer, there are two plasma boundaries, close to the inferred separatrices on the magnetosphere and magnetosheath sides (Ssp and Ssh) and two boundaries associated with the Alfvén waves (or Rotational Discontinuities, RDsp and RDsh). The data are consistent with these being launched from the reconnection site and the plasma distributions are well ordered and suggestive of the time elapsed since reconnection of the field lines observed. In each sub-layer between the boundaries the plasma distribution is different and is centered around the current sheet, responsible for magnetosheath acceleration. We show evidence for a velocity dispersion effect in the electron anisotropy that is consistent with the time elapsed since reconnection. In addition, new evidence is presented for the occurrence of partial reflection of magnetosheath electrons at the magnetopause current layer.

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

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

Classification and fault detection methods for fuel cell monitoring and quality control

In this paper, various types of fault detection methods for fuel cells are compared. For example, those that use a model based approach or a data driven approach or a combination of the two. The potential advantages and drawbacks of each method are discussed and comparisons between methods are made. In particular, classification algorithms are investigated, which separate a data set into classes or clusters based on some prior knowledge or measure of similarity. In particular, the application of classification methods to vectors of reconstructed currents by magnetic tomography or to vectors of magnetic field measurements directly is explored. Bases are simulated using the finite integration technique (FIT) and regularization techniques are employed to overcome ill-posedness. Fisher's linear discriminant is used to illustrate these concepts. Numerical experiments show that the ill-posedness of the magnetic tomography problem is a part of the classification problem on magnetic field measurements as well. This is independent of the particular working mode of the cell but influenced by the type of faulty behavior that is studied. The numerical results demonstrate the ill-posedness by the exponential decay behavior of the singular values for three examples of fault classes.

Ultralow-frequency modulation of whistler-mode wave growth

Measurements from ground-based magnetometers and riometers at auroral latitudes have demonstrated that energetic (~30-300keV) electron precipitation can be modulated in the presence of magnetic field oscillations at ultra-low frequencies. It has previously been proposed that an ultra-low frequency (ULF) wave would modulate field and plasma properties near the equatorial plane, thus modifying the growth rates of whistler-mode waves. In turn, the resulting whistler-mode waves would mediate the pitch-angle scattering of electrons resulting in ionospheric precipitation. In this paper, we investigate this hypothesis by quantifying the changes to the linear growth rate expected due to a slow change in the local magnetic field strength for parameters typical of the equatorial region around 6.6RE radial distance. To constrain our study, we determine the largest possible ULF wave amplitudes from measurements of the magnetic field at geosynchronous orbit. Using nearly ten years of observations from two satellites, we demonstrate that the variation in magnetic field strength due to oscillations at 2mHz does not exceed ±10% of the background field. Modifications to the plasma density and temperature anisotropy are estimated using idealised models. For low temperature anisotropy, there is little change in the whistler-mode growth rates even for the largest ULF wave amplitude. Only for large temperature anisotropies can whistler-mode growth rates be modulated sufficiently to account for the changes in electron precipitation measured by riometers at auroral latitudes.

Comparison of the open-closed separatrix in a global magnetospheric simulation with observations: the role of the ring current

The development of global magnetospheric models, such as Space Weather Modeling Framework (SWMF), which can accurately reproduce and track space weather processes has high practical utility. We present an interval on 5 June 1998, where the location of the polar cap boundary, or open-closed field line boundary (OCB), can be determined in the ionosphere using a combination of instruments during a period encompassing a sharp northward to southward interplanetary field turning. We present both point- and time-varying comparisons of the observed and simulated boundaries in the ionosphere and find that when using solely the coupled ideal magnetohydrodynamic magnetosphere-ionosphere model, the rate of change of the OCB to a southward turning of the interplanetary field is significantly faster than that computed from the observational data. However, when the inner magnetospheric module is incorporated, the modeling framework both qualitatively, and often quantitatively, reproduces many elements of the studied interval prior to an observed substorm onset. This result demonstrates that the physics of the inner magnetosphere is critical in shaping the boundary between open and closed field lines during periods of southward interplanetary magnetic field (IMF) and provides significant insight into the 3-D time-dependent behavior of the Earth's magnetosphere in response to a northward-southward IMF turning. We assert that during periods that do not include the tens of minutes surrounding substorm expansion phase onset, the coupled SWMF model may provide a valuable and reliable tool for estimating both the OCB and magnetic field topology over a wide range of latitudes and local times.

Electron trapping in shear Alfvén waves that power the aurora

Results from 1D Vlasov drift-kinetic plasma simulations reveal how and where auroral electrons are accelerated along Earth’s geomagnetic field. In the warm plasma sheet, electrons become trapped in shear Alfven waves, preventing immediate wave damping. As waves move to regions with larger vTe=vA, their parallel electric field decreases, and the trapped electrons escape their influence. The resulting electron distribution functions compare favorably with in situ observations, demonstrating for the first time a self-consistent link between Alfven waves and electrons that form aurora.

Near-Earth initiation of a terrestrial substorm

Despite the characterization of the auroral substorm more than 40 years ago, controversy still surrounds the processes triggering substorm onset initiation. That stretching of the Earth's magnetotail following the addition of new nightside magnetic flux from dayside reconnection powers the substorm is well understood; the trigger for explosive energy release at substorm expansion phase onset is not. Using ground-based data sets with unprecedented combined spatial and temporal coverage, we report the discovery of new localized and contemporaneous magnetic wave and small azimuthal scale auroral signature of substorm onset. These local auroral arc undulations and magnetic field signatures rapidly evolve on second time scales for several minutes in advance of the release of the auroral surge. We also present evidence from a conjugate geosynchronous satellite of the concurrent magnetic onset in space as the onset of magnetic pulsations in the ionosphere, to within technique error. Throughout this time period, the more poleward arcs that correspond to the auroral oval which maps to the central plasma sheet remain undisturbed. There is good evidence that flows from the midtail crossing the plasma sheet can generate north-south auroral structures, yet no such auroral forms are seen in this event. Our observations present a severe challenge to the standard hypothesis that magnetic reconnection in stretched magnetotail fields triggers onset, indicating substorm expansion phase initiation occurs on field lines that are close to the Earth, as bounded by observations at geosynchronous orbit and in the conjugate ionosphere.

The 22-Year Hale Cycle in cosmic ray flux: evidence for direct heliospheric modulation

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 as 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, which in turn causes the peak to be more dome-shaped 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, we also demonstrate that these polarity-dependent heliospheric differences are evident during the space-age but are much less clear in earlier data: using geomagnetic reconstructions, we show 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.

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

We investigate a coronal mass ejection (CME) propagating toward 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 is 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 three-dimensional (3-D) magnetospheric space weather simulation from the Community Coordinated Modeling Center. 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 values were then compared with the observations from specific equatorially located ground stations and showed 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 a 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.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr – Part 1: A new geomagnetic data composite

We present a new composite of geomagnetic activity which is designed to be as homogeneous in its construction as possible. This is done by only combining data that, by virtue of the locations of the source observatories used, have similar responses to solar wind and IMF (interplanetary magnetic field) variations. This will enable us (in Part 2, Lockwood et al., 2013a) to use the new index to reconstruct the interplanetary magnetic field, B, back to 1846 with a full analysis of errors. Allowance is made for the effects of secular change in the geomagnetic field. The composite uses interdiurnal variation data from Helsinki for 1845–1890 (inclusive) and 1893–1896 and from Eskdalemuir from 1911 to the present. The gaps are filled using data from the Potsdam (1891–1892 and 1897–1907) and the nearby Seddin observatories (1908–1910) and intercalibration achieved using the Potsdam–Seddin sequence. The new index is termed IDV(1d) because it employs many of the principles of the IDV index derived by Svalgaard and Cliver (2010), inspired by the u index of Bartels (1932); however, we revert to using one-day (1d) means, as employed by Bartels, because the use of near-midnight values in IDV introduces contamination by the substorm current wedge auroral electrojet, giving noise and a dependence on solar wind speed that varies with latitude. The composite is compared with independent, early data from European-sector stations, Greenwich, St Petersburg, Parc St Maur, and Ekaterinburg, as well as the composite u index, compiled from 2–6 stations by Bartels, and the IDV index of Svalgaard and Cliver. Agreement is found to be extremely good in all cases, except two. Firstly, the Greenwich data are shown to have gradually degraded in quality until new instrumentation was installed in 1915. Secondly, we infer that the Bartels u index is increasingly unreliable before about 1886 and overestimates the solar cycle amplitude between 1872 and 1883 and this is amplified in the proxy data used before 1872. This is therefore also true of the IDV index which makes direct use of the u index values.

Implications of the recent low solar minimum for the solar wind during the Maunder minimum

The behavior of the Sun and near-Earth space during grand solar minima is not understood; however, the recent long and low minimum of the decadal-scale solar cycle gives some important clues, with implications for understanding the solar dynamo and predicting space weather conditions. The speed of the near-Earth solar wind and the strength of the interplanetary magnetic field (IMF) embedded within it can be reliably reconstructed for before the advent of spacecraft monitoring using observations of geomagnetic activity that extend back to the mid-19th century. We show that during the solar cycle minima around 1879 and 1901 the average solar wind speed was exceptionally low, implying the Earth remained within the streamer belt of slow solar wind flow for extended periods. This is consistent with a broader streamer belt, which was also a feature of the recent low minimum (2009), and yields a prediction that the low near-Earth IMF during the Maunder minimum (1640-1700), as derived from models and deduced from cosmogenic isotopes, was accompanied by a persistent and relatively constant solar wind of speed roughly half the average for the modern era.

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

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