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.

The cloud chamber and CTR Wilson’s legacy to atmospheric science

2011 is the centenary year of the short paper (Wilson,1911) first describing the cloud chamber, the device for visualising high-energy charged particles which earned the Scottish physicist Charles Thomas Rees (‘CTR’) Wilson the 1927 Nobel Prize for physics. His many achievements in atmospheric science, some of which have current relevance, are briefly reviewed here. CTR Wilson’s lifetime of scientific research work was principally in atmospheric electricity at the Cavendish Laboratory, Cambridge; he was Reader in Electrical Meteorology from 1918 and Jacksonian Professor from 1925 to 1935. However, he is immortalised in physics for his invention of the cloud chamber, because of its great significance as an early visualisation tool for particles such as cosmic rays1 (Galison, 1997). Sir Lawrence Bragg summarised its importance:

Electromagnetic atmosphere-plasma coupling: the global atmospheric electric circuit

A description is given of the global atmospheric electric circuit operating between the Earth’s surface and the ionosphere. Attention is drawn to the huge range of horizontal and vertical spatial scales, ranging from 10−9 m to 1012 m, concerned with the many important processes at work. A similarly enormous range of time scales is involved from 10−6 s to 109 s, in the physical effects and different phenomena that need to be considered. The current flowing in the global circuit is generated by disturbed weather such as thunderstorms and electrified rain/shower clouds, mostly occurring over the Earth’s land surface. The profile of electrical conductivity up through the atmosphere, determined mainly by galactic cosmic ray ionization, is a crucial parameter of the circuit. Model simulation results on the variation of the ionospheric potential, ∼250 kV positive with respect to the Earth’s potential, following lightning discharges and sprites are summarized. Experimental results comparing global circuit variations with the neutron rate recorded at Climax, Colorado, are then discussed. Within the return (load) part of the circuit in the fair weather regions remote from the generators, charge layers exist on the upper and lower edges of extensive layer clouds; new experimental evidence for these charge layers is also reviewed. Finally, some directions for future research in the subject are suggested.

Whistler mode wave growth and propagation in the prenoon magnetosphere

Pitch-angle scattering of electrons can limit the stably trapped particle flux in the magnetosphere and precipitate energetic electrons into the ionosphere. Whistler-mode waves generated by a temperature anisotropy can mediate this pitch-angle scattering over a wide range of radial distances and latitudes, but in order to correctly predict the phase-space diffusion, it is important to characterise the whistler-mode wave distributions that result from the instability. We use previously-published observations of number density, pitch-angle anisotropy and phase space density to model the plasma in the quiet pre-noon magnetosphere (defined as periods when AE<100nT). We investigate the global propagation and growth of whistler-mode waves by studying millions of growing ray paths and demonstrate that the wave distribution at any one location is a superposition of many waves at different points along their trajectories and with different histories. We show that for observed electron plasma properties, very few raypaths undergo magnetospheric reflection, most rays grow and decay within 30 degrees of the magnetic equator. The frequency range of the wave distribution at large L can be adequately described by the solutions of the local dispersion relation, but the range of wavenormal angle is different. The wave distribution is asymmetric with respect to the wavenormal angle. The numerical results suggest that it is important to determine the variation of magnetospheric parameters as a function of latitude, as well as local time and L-shell.

A Riemann-Hilbert problem with a vanishing coefficient and applications to Toeplitz operators

We study the homogeneous Riemann-Hilbert problem with a vanishing scalar-valued continuous coefficient. We characterize non-existence of nontrivial solutions in the case where the coefficient has its values along several rays starting from the origin. As a consequence, some results on injectivity and existence of eigenvalues of Toeplitz operators in Hardy spaces are obtained.

Further evidence for variable synchrotron emission in XTE J1118+480 in outburst

We present simultaneous multicolor infrared and optical photometry of the black hole X-ray transient XTE J1118+480 during its short 2005 January outburst, supported by simultaneous X-ray observations. The variability is dominated by short timescales, ~10 s, although a weak superhump also appears to be present in the optical. The optical rapid variations, at least, are well correlated with those in X-rays. Infrared JHKs photometry, as in the previous outburst, exhibits especially large-amplitude variability. The spectral energy distribution (SED) of the variable infrared component can be fitted with a power law of slope α=-0.78+/-0.07, where F_ν~ν^α. There is no compelling evidence for evolution in the slope over five nights, during which time the source brightness decayed along almost the same track as seen in variations within the nights. We conclude that both short-term variability and longer timescale fading are dominated by a single component of constant spectral shape. We cannot fit the SED of the IR variability with a credible thermal component, either optically thick or thin. This IR SED is, however, approximately consistent with optically thin synchrotron emission from a jet. These observations therefore provide indirect evidence to support jet-dominated models for XTE J1118+480 and also provide a direct measurement of the slope of the optically thin emission, which is impossible, based on the average spectral energy distribution alone.

Multiwavelength Observations of EXO 0748-676. I. Reprocessing of X-Ray Bursts

We present high time-resolution multiwavelength observations of X-ray bursts in the low-mass X-ray binary UY Vol. Strong reprocessed signals are present in the ultraviolet and optical, lagged and smeared with respect to the X-rays. The addition of far-ultraviolet coverage for one burst allows much tighter constraints on the temperature and geometry of the reprocessing region than previously possible. A blackbody reprocessing model for this burst suggests a rise in temperatures during the burst from 18,000 to 35,000 K and an emitting area comparable to that expected for the disk and/or irradiated companion star. The lags are consistent with those expected. The single-zone blackbody model cannot reproduce the ratio of optical to ultraviolet flux during the burst, however. The discrepancy seems too large to explain with deviations from a local blackbody spectrum and more likely indicates that a range of reprocessing temperatures are required. Comparable results are derived from other bursts, and in particular the lag and smearing both appear shorter when the companion star is on the near side of the disk as predicted. The burst observed by HST also yielded a spectrum of the reprocessed light. It is dominated by continuum, with a spectral shape consistent with the temperatures derived from lightcurve modeling. Taken as a whole, our observations confirm the standard paradigm of prompt reprocessing distributed across the disk and companion star, with the response dominated by a thermalized continuum rather than by emission lines.

The heliospheric magnetic field

The heliospheric magnetic field (HMF) is the extension of the coronal magnetic field carried out into the solar system by the solar wind. It is the means by which the Sun interacts with planetary magnetospheres and channels charged particles propagating through the heliosphere. As the HMF remains rooted at the solar photosphere as the Sun rotates, the large-scale HMF traces out an Archimedean spiral. This pattern is distorted by the interaction of fast and slow solar wind streams, as well as the interplanetary manifestations of transient solar eruptions called coronal mass ejections. On the smaller scale, the HMF exhibits an array of waves, discontinuities, and turbulence, which give hints to the solar wind formation process. This review aims to summarise observations and theory of the small- and large-scale structure of the HMF. Solar-cycle and cycle-to-cycle evolution of the HMF is discussed in terms of recent spacecraft observations and pre-spaceage proxies for the HMF in geomagnetic and galactic cosmic ray records.

Evidence for solar wind modulation of lightning

The response of lightning rates over Europe to arrival of high speed solar wind streams at Earth is investigated using a superposed epoch analysis. Fast solar wind stream arrival is determined from modulation of the solar wind V y component, measured by the Advanced Composition Explorer spacecraft. Lightning rate changes around these event times are determined from the very low frequency arrival time difference (ATD) system of the UK Met Office. Arrival of high speed streams at Earth is found to be preceded by a decrease in total solar irradiance and an increase in sunspot number and Mg II emissions. These are consistent with the high speed stream's source being co-located with an active region appearing on the Eastern solar limb and rotating at the 27 d period of the Sun. Arrival of the high speed stream at Earth also coincides with a small (~1%) but rapid decrease in galactic cosmic ray flux, a moderate (~6%) increase in lower energy solar energetic protons (SEPs), and a substantial, statistically significant increase in lightning rates. These changes persist for around 40 d in all three quantities. The lightning rate increase is corroborated by an increase in the total number of thunder days observed by UK Met stations, again persisting for around 40 d after the arrival of a high speed solar wind stream. This result appears to contradict earlier studies that found an anti-correlation between sunspot number and thunder days over solar cycle timescales. The increase in lightning rates and thunder days that we observe coincides with an increased flux of SEPs which, while not being detected at ground level, nevertheless penetrate the atmosphere to tropospheric altitudes. This effect could be further amplified by an increase in mean lightning stroke intensity that brings more strokes above the detection threshold of the ATD system. In order to remove any potential seasonal bias the analysis was repeated for daily solar wind triggers occurring during the summer months (June to August). Though this reduced the number of solar wind triggers to 32, the response in both lightning and thunder day data remained statistically significant. This modulation of lightning by regular and predictable solar wind events may be beneficial to medium range forecasting of hazardous weather.

Detection of lower tropospheric responses to solar energetic particles at midlatitudes

Solar energetic particles (SEPs) occasionally contribute additional atmospheric ionization beyond that arising from the usual galactic cosmic ray background. During an SEP event associated with a solar flare on April 11, 2013, the vertical ionization rate profile obtained using a balloon-borne detector showed enhanced ionization with a 26% increase at 20 km, over Reading, United Kingdom. Fluctuations in atmospheric electrical parameters were also detected at the surface, beneath the balloon’s trajectory. As no coincident changes in geomagnetism occurred, the electrical fluctuations are very likely to be associated with increased ionization, as observed by the balloon measurements. The lack of response of surface neutron monitors during this event indicates that energetic particles that are not detected at the surface by neutron monitors can nevertheless enter and influence the atmosphere’s weather-generating regions.

Vertical profile measurements of lower troposphere ionisation

Vertical soundings of the atmospheric ion production rate have been obtained from Geiger counters integrated with conventional meteorological radiosondes. In launches made from Reading (UK) during 2013-2014, the Regener-Pfotzer ionisation maximum was at an altitude equivalent to a pressure of (63.1±2.4) hPa, or, expressed in terms of the local air density, (0.101±0.005) kgm−3. The measured ionisation profiles have been evaluated against the Usoskin-Kovaltsov model and, separately, surface neutron monitor data from Oulu. Model ionisation rates agree well with the observed cosmic ray ionisation below 20 km altitude. Above 10 km, the measured ionisation rates also correlate well with simultaneous neutron monitor data, although, consistently with previous work, measured variability at the ionisation maximum is greater than that found by the neutron monitor. However, in the lower atmosphere (below 5 km altitude), agreement between the measurements and simultaneous neutron monitor data is poor. For studies of transient lower atmosphere phenomena associated with cosmic ray ionisation, this indicates the need for in situ ionisation measurements and improved lower atmosphere parameterisations.

Modulation of UK lightning by heliospheric magnetic field polarity

Observational studies have reported solar magnetic modulation of terrestrial lightning on a range of time scales, from days to decades. The proposed mechanism is two-step: lightning rates vary with galactic cosmic ray (GCR) flux incident on Earth, either via changes in atmospheric conductivity and/or direct triggering of lightning. GCR flux is, in turn, primarily controlled by the heliospheric magnetic field (HMF) intensity. Consequently, global changes in lightning rates are expected. This study instead considers HMF polarity, which doesnʼt greatly affect total GCR flux. Opposing HMF polarities are, however, associated with a 40–60% difference in observed UK lightning and thunder rates. As HMF polarity skews the terrestrial magnetosphere from its nominal position, this perturbs local ionospheric potential at high latitudes and local exposure to energetic charged particles from the magnetosphere. We speculate as to the mechanism(s) by which this may, in turn, redistribute the global location and/or intensity of thunderstorm activity.

Phase calibration of the EISCAT Svalbard Radar interferometer using optical satellite signatures

The link between natural ion-line enhancements in radar spectra and auroral activity has been the subject of recent studies but conclusions have been limited by the spatial and temporal resolution previously available. The next challenge is to use shorter sub-second integration times in combination with interferometric programmes to resolve spatial structure within the main radar beam, and so relate enhanced filaments to individual auroral rays. This paper presents initial studies of a technique, using optical and spectral satellite signatures, to calibrate the received phase of a signal with the position of the scattering source along the interferometric baseline of the EISCAT Svalbard Radar. It is shown that a consistent relationship can be found only if the satellite passage through the phase fringes is adjusted from the passage predicted by optical tracking. This required adjustment is interpreted as being due to the vector between the theoretical focusing points of the two antennae, i.e. the true radar baseline, differing from the baseline obtained by survey between the antenna foot points. A method to obtain a measurement of the true interferometric baseline using multiple satellite passes is outlined.

Twenty-three cycles of changing open solar magnetic flux

This paper presents a comparison of various estimates of the open solar flux, deduced from measurements of the interplanetary magnetic field, from the aa geomagnetic index and from photospheric magnetic field observations. The first two of these estimates are made using the Ulysses discovery that the radial heliospheric field is approximately independent of heliographic latitude, the third makes use of the potential-field source surface method to map the total flux through the photosphere to the open flux at the top of the corona. The uncertainties associated with using the Ulysses result are 5%, but the effects of the assumptions of the potential field source surface method are harder to evaluate. Nevertheless, the three methods give similar results for the last three solar cycles when the data sets overlap. In 11-year running means, all three methods reveal that 1987 marked a significant peak in the long-term variation of the open solar flux. This peak is close to the solar minimum between sunspot cycles 21 and 22, and consequently the mean open flux (averaged from minimum to minimum) is similar for these two cycles. However, this similarity between cycles 21 and 22 in no way implies that the open flux is constant. The long-term variation shows that these cycles are fundamentally different in that the average open flux was rising during cycle 21 (from consistently lower values in cycle 20 and toward the peak in 1987) but was falling during cycle 22 (toward consistently lower values in cycle 23). The estimates from the geomagnetic aa index are unique as they extend from 1842 onwards (using the Helsinki extension). This variation gives strong anticorrelations, with very high statistical significance levels, with cosmic ray fluxes and with the abundances of the cosmogenic isotopes that they produce. Thus observations of photospheric magnetic fields, of cosmic ray fluxes, and of cosmogenic isotope abundances all support the long-term drifts in open solar flux reported by Lockwood et al. [1999a, 1999b].