Long-term variations in the open solar flux and possible links to Earth’s climate

Recent paleoclimate studies provide strong evidence for an association between cosmogenic isotope production and Earth’s climate throughout the holecene. These isotopes are generated by the bombardment of Earth’s atmosphere by galactic cosmic rays, the fluxes of which vary in approximately inverse proportion to the total open magnetic flux of the Sun. This paper discusses how results from the Ulysses spacecraft allow us to quantify the open solar flux from observations of near-Earth interplanetary space and to study its long-term variations using the homogeneous record of geomagnetic activity. A study of the results and of their accuracy is presented. The two proposed mechanisms that could lead to the open solar flux being a good proxy for solar-induced climate change are discussed: the first is the modulation of the production of some types of cloud by the air ions produced by cosmic rays; the second is a variation in the total or spectral solar irradiance, in association with changes in the open flux. Some implications for our understanding of anthropogenic climate change are discussed.

An evaluation of the correlation between open solar flux and total solar irradiance

The correlation between the coronal source flux F_{S} and the total solar irradiance I_{TS} is re-evaluated in the light of an additional 5 years' data from the rising phase of solar cycle 23 and also by using cosmic ray fluxes detected at Earth. Tests on monthly averages show that the correlation with F_{S} deduced from the interplanetary magnetic field (correlation coefficient, r = 0.62) is highly significant (99.999%), but that there is insufficient data for the higher correlation with annual means (r = 0.80) to be considered significant. Anti-correlations between I_{TS} and cosmic ray fluxes are found in monthly data for all stations and geomagnetic rigidity cut-offs (r ranging from −0.63 to −0.74) and these have significance levels between 85% and 98%. In all cases, the t is poorest for the earliest data (i.e., prior to 1982). Excluding these data improves the anticorrelation with cosmic rays to r = −0:93 for one-year running means. Both the interplanetary magnetic field data and the cosmic ray fluxes indicate that the total solar irradiance lags behind the open solar flux with a delay that is estimated to have an optimum value of 2.8 months (and is within the uncertainty range 0.8-8.0 months at the 90% level).

Long-term changes in the solar photosphere associated with changes in the coronal source flux

Using sunspot observations from Greenwich and Mount Wilson, we show that the latitudinal spread of sunspot groups has increased since 1874, in a manner that closely mirrors the long-term (similar to 100 year) changes in the coronal source flux, F-s, as inferred from geomagnetic activity. This latitude spread is shown to be well correlated with the flux emergence rate required by the model of the coronal source flux variation by Solanki er al. [2000]. The time constant for the decay of this open flux is found to be 3.6 +/-0.8 years. Using this value, and quantifying the photospheric flux emergence rate using the latitudinal spread of sunspot groups, the model reproduces the observed coronal source flux variation. The ratio of the 100-year drift to the solar cycle amplitude for the flux emergence rate is found to be half of the same ratio for F-s.

Long term changes in EUV and X-ray emissions from the solar corona and chromosphere as measured by the response of the Earth’s ionosphere during total solar eclipses from 1932 to 1999

Measurements of the ionospheric E region during total solar eclipses in the period 1932-1999 have been used to investigate the fraction of Extreme Ultra Violet and soft X-ray radiation, phi, that is emitted from the limb corona and chromosphere. The relative apparent sizes of the Moon and the Sun are different for each eclipse, and techniques are presented which correct the measurements and, therefore, allow direct comparisons between different eclipses. The results show that the fraction of ionising radiation emitted by the limb corona has a clear solar cycle variation and that the underlying trend shows this fraction has been increasing since 1932. Data from the SOHO spacecraft are used to study the effects of short-term variability and it is shown that the observed long-term rise in phi has a negligible probability of being a chance occurrence.

Solar causes of the long-term increase in geomagnetic activity

We analyze the causes of the century-long increase in geomagnetic activity, quantified by annual means of the aa index, using observations of interplanetary space, galactic cosmic rays, the ionosphere, and the auroral electrojet, made during the last three solar cycles. The effects of changes in ionospheric conductivity, the Earth's dipole tilt, and magnetic moment are shown to be small; only changes in near-Earth interplanetary space make a significant contribution to the long-term increase in activity. We study the effects of the interplanetary medium by applying dimensional analysis to generate the optimum solar wind-magnetosphere energy coupling function, having an unprecedentedly high correlation coefficient of 0.97. Analysis of the terms of the coupling function shows that the largest contributions to the drift in activity over solar cycles 20-22 originate from rises in the average interplanetary magnetic field (IMF) strength, solar wind concentration, and speed; average IMF orientation has grown somewhat less propitious for causing geomagnetic activity. The combination of these factors explains almost all of the 39% rise in aa observed over the last three solar cycles. Whereas the IMF strength varies approximately in phase with sunspot numbers, neither its orientation nor the solar wind density shows any coherent solar cycle variation. The solar wind speed peaks strongly in the declining phase of even-numbered cycles and can be identified as the chief cause of the phase shift between the sunspot numbers and the aa index. The rise in the IMF magnitude, the largest single contributor to the drift in geomagnetic activity, is shown to be caused by a rise in the solar coronal magnetic field, consistent with a rise in the coronal source field, modeled from photospheric observations, and an observed decay in cosmic ray fluxes.

Predicting solar disturbance effects on navigation systems

A variety of operational systems are vulnerable to disruption by solar disturbances brought to the Earth by the solar wind. Of particular importance to navigation systems are energetic charged particles which can generate temporary malfunctions and permanent damage in satellites. Modern spacecraft technology may prove to be particularly at risk during the next maximum of the solar cycle. In addition, the associated ionospheric disturbances cause phase shifts of transionospheric and ionosphere-reflected signals, giving positioning errors and loss of signal for GPS and Loran-C positioning systems and for over-the-horizon radars. We now have sufficient understanding of the solar wind, and how it interacts with the Earth's magnetic field, to predict statistically the likely effects on operational systems over the next solar cycle. We also have a number of advanced ways of detecting and tracking these disturbances through space but we cannot, as yet, provide accurate forecasts of individual disturbances that could be used to protect satellites and to correct errors. In addition, we have recently discovered long-term changes in the Sun, which mean that the number and severity of the disturbances to operational systems are increasing.

Fathoming solar effects in Neptune’s atmosphere

Long-duration observations of Neptune’s brightness in two visible wavelengths provide a disk-averaged estimate of its atmospheric aerosol. Brightness variations were previously associated with the 11-year solar cycle, through solar-modulated mechanisms linked with either ultra-violet (UV) or galactic cosmic ray (GCR) effects on atmospheric particles. Here we use a recently extended brightness dataset (1972-2014), with physically realistic modelling to show that rather than alternatives, UV and GCR are likely to be modulating Neptune’s atmosphere in combination. The importance of GCR is further supported by the response of Neptune's atmosphere to an intermittent 1.5 to 1.9 year periodicity, which occurred preferentially in GCR (not UV) during the mid-1980s. This periodicity was detected both at Earth, and in GCR measured by Voyager 2, then near Neptune. A similar coincident variability in Neptune’s brightness suggests nucleation onto GCR ions. Both GCR and UV mechanisms may occur more rapidly than the subsequent atmospheric particle transport.

The role of eddies in driving the tropospheric response to stratospheric heating perturbations

A simplified general circulation model has been used to investigate the chain of causality whereby changes in tropospheric circulation and temperature are produced in response to stratospheric heating perturbations. Spinup ensemble experiments have been performed to examine the evolution of the tropospheric circulation in response to such perturbations. The primary aim of these experiments is to investigate the possible mechanisms whereby a tropospheric response to changing solar activity over the 11-yr solar cycle could be produced in response to heating of the equatorial lower stratosphere. This study therefore focuses on a stratospheric heating perturbation in which the heating is largest in the tropics. For comparison, experiments are also performed in which the stratosphere is heated uniformly at all latitudes and in which it is heated preferentially in the polar region. Thus, the mechanisms discussed have a wider relevance for the impact of stratospheric perturbations on the troposphere. The results demonstrate the importance of changing eddy momentum fluxes in driving the tropospheric response. This is confirmed by the lack of a similar response in a zonally symmetric model with fixed eddy forcing. Furthermore, it is apparent that feedback between the tropospheric eddy fluxes and tropospheric circulation changes is required to produce the full model response. The quasigeostrophic index of refraction is used to diagnose the cause of the changes in eddy behavior. It is demonstrated that the latitudinal extent of stratospheric heating is important in determining the direction of displacement of the tropospheric jet and storm track.

Stratospheric temperature trends: impact of ozone variability and the QBO

In most climate simulations used by the Intergovernmental Panel on Climate Change 2007 fourth assessment report, stratospheric processes are only poorly represented. For example, climatological or simple specifications of time-varying ozone concentrations are imposed and the quasi-biennial oscillation (QBO) of equatorial stratospheric zonal wind is absent. Here we investigate the impact of an improved stratospheric representation using two sets of perturbed simulations with the Hadley Centre coupled ocean atmosphere model HadGEM1 with natural and anthropogenic forcings for the 1979–2003 period. In the first set of simulations, the usual zonal mean ozone climatology with superimposed trends is replaced with a time series of observed zonal mean ozone distributions that includes interannual variability associated with the solar cycle, QBO and volcanic eruptions. In addition to this, the second set of perturbed simulations includes a scheme in which the stratospheric zonal wind in the tropics is relaxed to appropriate zonal mean values obtained from the ERA-40 re-analysis, thus forcing a QBO. Both of these changes are applied strictly to the stratosphere only. The improved ozone field results in an improved simulation of the stepwise temperature transitions observed in the lower stratosphere in the aftermath of the two major recent volcanic eruptions. The contribution of the solar cycle signal in the ozone field to this improved representation of the stepwise cooling is discussed. The improved ozone field and also the QBO result in an improved simulation of observed trends, both globally and at tropical latitudes. The Eulerian upwelling in the lower stratosphere in the equatorial region is enhanced by the improved ozone field and is affected by the QBO relaxation, yet neither induces a significant change in the upwelling trend.

Reconciling the electron counterstreaming and dropout occurrence rates with the heliospheric flux budget

Counterstreaming electrons (CSEs) are treated as signatures of closed magnetic flux, i.e., loops connected to the Sun at both ends. However, CSEs at 1 AU likely fade as the apex of a closed loop passes beyond some distance R, owing to scattering of the sunward beam along its continually increasing path length. The remaining antisunward beam at 1 AU would then give a false signature of open flux. Subsequent opening of a loop at the Sun by interchange reconnection with an open field line would produce an electron dropout (ED) at 1 AU, as if two open field lines were reconnecting to completely disconnect from the Sun. Thus EDs can be signatures of interchange reconnection as well as the commonly attributed disconnection. We incorporate CSE fadeout into a model that matches time-varying closed flux from interplanetary coronal mass ejections (ICMEs) to the solar cycle variation in heliospheric flux. Using the observed occurrence rate of CSEs at solar maximum, the model estimates R ∼ 8–10 AU. Hence we demonstrate that EDs should be much rarer than CSEs at 1 AU, as EDs can only be detected when the juncture points of reconnected field lines lie sunward of the detector, whereas CSEs continue to be detected in the legs of all loops that have expanded beyond the detector, out to R. We also demonstrate that if closed flux added to the heliosphere by ICMEs is instead balanced by disconnection elsewhere, then ED occurrence at 1 AU would still be rare, contrary to earlier expectations.

Coronal mass ejections and magnetic flux buildup in the heliosphere

To test for magnetic flux buildup in the heliosphere from coronal mass ejections (CMEs), we simulate heliospheric flux as a constant background open flux with a time-varying interplanetary CME (ICME) contribution. As flux carried by ejecta can only contribute to the heliospheric flux budget while it remains closed, the ICME flux opening rate is an important factor. Two separate forms for the ICME flux opening rate are considered: (1) constant and (2) exponentially decaying with time. Coronagraph observations are used to determine the CME occurrence rates, while in situ observations are used to estimate the magnetic flux content of a typical ICME. Both static equilibrium and dynamic simulations, using the constant and exponential ICME flux opening models, require flux opening timescales of ∼50 days in order to match the observed doubling in the magnetic field intensity at 1 AU over the solar cycle. Such timescales are equivalent to a change in the ICME closed flux of only ∼7–12% between 1 and 5 AU, consistent with CSE signatures; no flux buildup results. The dynamic simulation yields a solar cycle flux variation with high variability that matches the overall variability of the observed magnetic field intensity remarkably well, including the double peak forming the Gnevyshev gap.

Cosmic rays, clouds, and climate

It has been proposed that Earth's climate could be affected by changes in cloudiness caused by variations in the intensity of galactic cosmic rays in the atmosphere. This proposal stems from an observed correlation between cosmic ray intensity and Earth's average cloud cover over the course of one solar cycle. Some scientists question the reliability of the observations, whereas others, who accept them as reliable, suggest that the correlation may be caused by other physical phenomena with decadal periods or by a response to volcanic activity or El Niño. Nevertheless, the observation has raised the intriguing possibility that a cosmic ray–cloud interaction may help explain how a relatively small change in solar output can produce much larger changes in Earth's climate. Physical mechanisms have been proposed to explain how cosmic rays could affect clouds, but they need to be investigated further if the observation is to become more than just another correlation among geophysical variables.

Mid-nineteenth century smoke concentrations near London

It has been proposed that Earth's climate could be affected by changes in cloudiness caused by variations in the intensity of galactic cosmic rays in the atmosphere. This proposal stems from an observed correlation between cosmic ray intensity and Earth's average cloud cover over the course of one solar cycle. Some scientists question the reliability of the observations, whereas others, who accept them as reliable, suggest that the correlation may be caused by other physical phenomena with decadal periods or by a response to volcanic activity or El Niño. Nevertheless, the observation has raised the intriguing possibility that a cosmic ray–cloud interaction may help explain how a relatively small change in solar output can produce much larger changes in Earth's climate. Physical mechanisms have been proposed to explain how cosmic rays could affect clouds, but they need to be investigated further if the observation is to become more than just another correlation among geophysical variables.

Coupled chemistry climate model simulations of stratospheric temperatures and their trends for the recent past

Temperature results from multi-decadal simulations of coupled chemistry climate models for the recent past are analyzed using multi-linear regression including a trend, solar cycle, lower stratospheric tropical wind, and volcanic aerosol terms. The climatology of the models for recent years is in good agreement with observations for the troposphere but the model results diverge from each other and from observations in the stratosphere. Overall, the models agree better with observations than in previous assessments, primarily because of corrections in the observed temperatures. The annually averaged global and polar temperature trends simulated by the models are generally in agreement with revised satellite observations and radiosonde data over much of their altitude range. In the global average, the model trends underpredict the radiosonde data slightly at the top of the observed range. Over the Antarctic some models underpredict the temperature trend in the lower stratosphere, while others overpredict the trends

An update of observed stratospheric temperature trends

An updated analysis of observed stratospheric temperature variability and trends is presented on the basis of satellite, radiosonde, and lidar observations. Satellite data include measurements from the series of NOAA operational instruments, including the Microwave Sounding Unit covering 1979–2007 and the Stratospheric Sounding Unit (SSU) covering 1979–2005. Radiosonde results are compared for six different data sets, incorporating a variety of homogeneity adjustments to account for changes in instrumentation and observational practices. Temperature changes in the lower stratosphere show cooling of 0.5 K/decade over much of the globe for 1979–2007, with some differences in detail among the different radiosonde and satellite data sets. Substantially larger cooling trends are observed in the Antarctic lower stratosphere during spring and summer, in association with development of the Antarctic ozone hole. Trends in the lower stratosphere derived from radiosonde data are also analyzed for a longer record (back to 1958); trends for the presatellite era (1958–1978) have a large range among the different homogenized data sets, implying large trend uncertainties. Trends in the middle and upper stratosphere have been derived from updated SSU data, taking into account changes in the SSU weighting functions due to observed atmospheric CO2 increases. The results show mean cooling of 0.5–1.5 K/decade during 1979–2005, with the greatest cooling in the upper stratosphere near 40–50 km. Temperature anomalies throughout the stratosphere were relatively constant during the decade 1995–2005. Long records of lidar temperature measurements at a few locations show reasonable agreement with SSU trends, although sampling uncertainties are large in the localized lidar measurements. Updated estimates of the solar cycle influence on stratospheric temperatures show a statistically significant signal in the tropics (30N–S), with an amplitude (solar maximum minus solar minimum) of 0.5 K (lower stratosphere) to 1.0 K (upper stratosphere).