Field-aligned plasma flow in the quiet, mid-latitude ionosphere deduced from topside soundings

A method for quantifying diffusive flows of O+ ions in the topside ionosphere from satellite soundings is described. A departure from diffusive equilibrium alters the shape of the plasma scale-height profile near the F2-peak where ion-neutral frictional drag is large. The effect enables the evaluation of , the field-aligned flux of O+ ions relative to the neutral oxygen atom gas, using MSIS model values for the neutral thermospheric densities and temperature. Upward flow values are accurate to within about 10%, the largest sources of error being the MSIS prediction for the concentration of oxygen atoms and the plasma temperature gradient deduced from the sounding. Downward flux values are only determined to within 20%. From 60,000 topside soundings, taken at the minimum and rising phase of the solar cycle, a total of 1098 mean scale-height profiles are identified for which no storm sudden commencement had occurred in the previous 12 days and for which Kp was less than 2o, each mean profile being an average of about six soundings. A statistical study ofdeduced from these profiles shows the diurnal cycle of O+ flow in the quiet, topside ionosphere at mid-latitudes and its seasonal variations. The differences betweenand ion flux observations from incoherent scatter radars are considered using the meridional thermospheric winds predicted by a global, three-dimensional model. The mean interhemispheric flow from summer to winter is compared with predictions by a numerical model of the protonospheric coupling of conjugate ionospheres for up to 6 days following a geomagnetic storm. The observed mean (of order 3 × 1016 ions day−1 along a flux tube of area 1 m2 at 1000 km) is larger than predicted for day 6 and the suggested explanation is a decrease in upward flows from the winter, daytime ionosphere between the sixth and twelfth days.

Extracting planetary waves from geomagnetic time series using Empirical Mode Decomposition

Empirical Mode Decomposition is presented as an alternative to traditional analysis methods to decompose geomagnetic time series into spectral components. Important comments on the algorithm and its variations will be given. Using this technique, planetary wave modes of 5-, 10-, and 16-day mean periods can be extracted from magnetic field components of three different stations in Germany. In a second step, the amplitude modulation functions of these wave modes can be shown to contain significant contribution from solar cycle variation through correlation with smoothed sunspot numbers. Additionally, the data indicate connections with geomagnetic jerk occurrences, supported by a second set of data providing reconstructed near-Earth magnetic field for 150 years. Usually attributed to internal dynamo processes within the Earth's outer core, the question of who is impacting whom will be briefly discussed here.

The Maunder minimum (1645-1715) was indeed a grand minimum: a reassessment of multiple datasets

Aims. Although the time of the Maunder minimum (1645–1715) is widely known as a period of extremely low solar activity, it is still being debated whether solar activity during that period might have been moderate or even higher than the current solar cycle (number 24). We have revisited all existing evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods. We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of the sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results. The level of solar activity during the Maunder minimum is reassessed on the basis of all available datasets. Conclusions. We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely lower than during the Dalton minimum of around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected with a high confidence level.

Tests of sunspot number sequences: 3. Effects of regression procedures on the calibration of historic sunspot data

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 RB 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 RB are then re-scaled to the full observed RGO group number RA using a variety of regression techniques. It is found that a very high correlation between RA and RB (rAB > 0.98) does not prevent large errors in the intercalibration (for example sunspot maximum values can be over 30 % too large even for such levels of rAB). In generating the backbone sunspot number (RBB), Svalgaard and Schatten (2015, this issue) 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.

Observations of the diurnal cycle of outgoing longwave radiation from the Geostationary Earth Radiation Budget instrument

The Geostationary Earth Radiation Budget instrument on Meteosat-8, located over Africa, provides unprecedented temporal sampling (~17 minutes) of the broadband emitted thermal and reflected solar radiances. We analyse the diurnal cycle of the outgoing longwave radiation (OLR) fluxes derived from the thermal radiances for July 2006. Principal component (PC) analysis separates the signals of the surface temperature response to solar heating and of the development of convective clouds. The first two PCs explain most of the OLR variations: PC1 (surface heating) explains 82.3% of the total variance and PC2 (cloud development) explains 12.8% of the variance. Convection is initiated preferentially over mountainous regions and the cloud then advects downstream in the ambient flow. Diurnal variations are much weaker over the oceans, but a coherent signal over the Gulf of Guinea suggests that the cloudiness is modulated by the diurnally varying contrast between the Gulf and the adjacent land mass.

Influence of the seasonal cycle on the termination of El Nino events in a coupled general circulation model

In this study, the mechanisms leading to the El Nino peak and demise are explored through a coupled general circulation model ensemble approach evaluated against observations. The results here suggest that the timing of the peak and demise for intense El Nino events is highly predictable as the evolution of the coupled system is strongly driven by a southward shift of the intense equatorial Pacific westerly anomalies during boreal winter. In fact, this systematic late-year shift drives an intense eastern Pacific thermocline shallowing, constraining a rapid El Nino demise in the following months. This wind shift results from a southward displacement in winter of the central Pacific warmest SSTs in response to the seasonal evolution of solar insolation. In contrast, the intensity of this seasonal feedback mechanism and its impact on the coupled system are significantly weaker in moderate El Nino events, resulting in a less pronounced thermocline shallowing. This shallowing transfers the coupled system into an unstable state in spring but is not sufficient to systematically constrain the equatorial Pacific evolution toward a rapid El Nino termination. However, for some moderate events, the occurrence of intense easterly wind anomalies in the eastern Pacific during that period initiate a rapid surge of cold SSTs leading to La Nina conditions. In other cases, weaker trade winds combined with a slightly deeper thermocline allow the coupled system to maintain a broad warm phase evolving through the entire spring and summer and a delayed El Nino demise, an evolution that is similar to the prolonged 1986/87 El Nino event. La Nina events also show a similar tendency to peak in boreal winter, with characteristics and mechanisms mainly symmetric to those described for moderate El Nino cases.

The importance of the diurnal and annual cycle of air traffic for contrail radiative forcing

Air traffic condensation trails, or contrails, are believed to have a net atmospheric warming effect(1), although one that is currently small compared to that induced by other sources of human emissions. However, the comparably large growth rate of air traffic requires an improved understanding of the resulting impact of aircraft radiative forcing on climate(2). Contrails have an effect on the Earth's energy balance similar to that of high thin ice clouds(3). Their trapping of outgoing longwave radiation emitted by the Earth and atmosphere (positive radiative forcing) is partly compensated by their reflection of incoming solar radiation (negative radiative forcing). On average, the longwave effect dominates and the net contrail radiative forcing is believed to be positive(1,2,4). Over daily and annual timescales, varying levels of air traffic, meteorological conditions, and solar insolation influence the net forcing effect of contrails. Here we determine the factors most important for contrail climate forcing using a sophisticated radiative transfer model(5,6) for a site in southeast England, located in the entrance to the North Atlantic flight corridor. We find that night-time flights during winter (December to February) are responsible for most of the contrail radiative forcing. Night flights account for only 25 per cent of daily air traffic, but contribute 60 to 80 per cent of the contrail forcing. Further, winter flights account for only 22 per cent of annual air traffic, but contribute half of the annual mean forcing. These results suggest that flight rescheduling could help to minimize the climate impact of aviation.

The accuracy of using the Ulysses result of the spatial invariance of the radial heliospheric field to compute the open solar flux

We survey observations of the radial magnetic field in the heliosphere as a function of position, sunspot number, and sunspot cycle phase. We show that most of the differences between pairs of simultaneous observations, normalized using the square of the heliocentric distance and averaged over solar rotations, are consistent with the kinematic "flux excess" effect whereby the radial component of the frozen-in heliospheric field is increased by longitudinal solar wind speed structure. In particular, the survey shows that, as expected, the flux excess effect at high latitudes is almost completely absent during sunspot minimum but is almost the same as within the streamer belt at sunspot maximum. We study the uncertainty inherent in the use of the Ulysses result that the radial field is independent of heliographic latitude in the computation of the total open solar flux: we show that after the kinematic correction for the excess flux effect has been made it causes errors that are smaller than 4.5%, with a most likely value of 2.5%. The importance of this result for understanding temporal evolution of the open solar flux is reviewed.

Solar induced climate effects

Solar electromagnetic radiation powers Earth’s climate system and, consequently, it is often naively assumed that changes in this solar output must be responsible for changes in Earth’s climate. However, the Sun is close to a blackbody radiator and so emits according to its surface temperature and the huge thermal time constant of the outer part of the Sun limits the variability in surface temperature and hence output. As a result, on all timescales of interest, changes in total power output are limited to small changes in effective surface temperature (associated with magnetic fields) and potential, although as yet undetected, solar radius variations. Larger variations are seen in the UV part of the spectrum which is emitted from the lower solar atmosphere (the chromosphere) and which influences Earth’s stratosphere. There is interest in“top-down” mechanisms whereby solar UV irradiance modulates stratospheric temperatures and winds which, in turn, may influence the underlying troposphere where Earth’s climate and weather reside. This contrasts with “bottom-up” effects in which the small total solar irradiance (dominated by the visible and near-IR) variations cause surface temperature changes which drive atmospheric circulations. In addition to these electromagnetic outputs, the Sun modulates energetic particle fluxes incident on the Earth. Solar Energetic Particles (SEP) are emitted by solar flares and from the shock fronts ahead of supersonic (and super-Alfvenic) ejections of material from the solar atmosphere. These SEPs enhance the destruction of polar stratospheric ozone which could be an additional form of top-down climate forcing. Even more energetic are Galactic Cosmic Rays (GCRs). These particles are not generated by the Sun, rather they originate at the shock fronts emanating from violent galactic events such as supernovae explosions; however, the expansion of the solar magnetic field into interplanetary space means that the Sun modulates the number of GCRs reaching Earth. These play a key role in enabling Earth’s global electric (thunderstorm) circuit and it has been proposed that they also modulate the formation of clouds. Both electromagnetic and corpuscular solar effects are known to vary over the solar magnetic cycle which is typically between 10 and 14 yrs in length (with an average close to 11 yrs). The solar magnetic field polarity at any one phase of one of these activity cycles is opposite to that at the same phase of the next cycle and this influences some phenomena, for example GCRs, which therefore show a 22 yr (“Hale”) cycle on average. Other phenomena, such as irradiance modulation, do not depend on the polarity of the magnetic field and so show only the basic 11-yr activity cycle. However, any effects on climate are much more significant for solar drifts over centennial timescales. This chapter discusses and evaluates potential effects on Earth’s climate system of variations in these solar inputs. Because of the great variety of proposed mechanisms, the wide range of timescales studied (from days to millennia) and the many debates (often triggered by the application of inadequate statistical methods), the literature on this subject is vast, complex, divergent and rapidly changing: consequently the number of references cited in this review is very large (yet still only a small fraction of the total).

The last glacial cycle: transient simulations with an AOGCM

A number of transient climate runs simulating the last 120kyr have been carried out using FAMOUS, a fast atmosphere-ocean general circulation model (AOGCM). This is the first time such experiments have been done with a full AOGCM, providing a three-dimensional simulation of both atmosphere and ocean over this period. Our simulation thus includes internally generated temporal variability over periods from days to millennia, and physical, detailed representations of important processes such as clouds and precipitation. Although the model is fast, computational restrictions mean that the rate of change of the forcings has been increased by a factor of 10, making each experiment 12kyr long. Atmospheric greenhouse gases (GHGs), northern hemisphere ice sheets and variations in solar radiation arising from changes in the Earth's orbit are treated as forcing factors, and are applied either separately or combined in different experiments. The long-term temperature changes on Antarctica match well with reconstructions derived from ice-core data, as does variability on timescales longer than 10 kyr. Last Glacial Maximum (LGM) cooling on Greenland is reasonably well simulated, although our simulations, which lack ice-sheet meltwater forcing, do not reproduce the abrupt, millennial scale climate shifts seen in northern hemisphere climate proxies or their slower southern hemisphere counterparts. The spatial pattern of sea surface cooling at the LGM matches proxy reconstructions reasonably well. There is significant anti-correlated variability in the strengths of the Atlantic Meridional Overturning Circulation (AMOC) and the Antarctic Circumpolar Current (ACC) on timescales greater than 10kyr in our experiments. We find that GHG forcing weakens the AMOC and strengthens the ACC, whilst the presence of northern hemisphere ice-sheets strengthens the AMOC and weakens the ACC. The structure of the AMOC at the LGM is found to be sensitive to the details of the ice-sheet reconstruction used. The precessional component of the orbital forcing induces ~20kyr oscillations in the AMOC and ACC, whose amplitude is mediated by changes in the eccentricity of the Earth's orbit. These forcing influences combine, to first order, in a linear fashion to produce the mean climate and ocean variability seen in the run with all forcings.

The global atmospheric water cycle

Water vapour plays a key role in the Earth's energy balance. Almost 50% of the absorbed solar radiation at the surface is used to cool the surface, through evaporation, and warm the atmosphere, through release of latent heat. Latent heat is the single largest factor in warming the atmosphere and in transporting heat from low to high latitudes. Water vapour is also the dominant greenhouse gas and contributes to a warming of the climate system by some 24°C (Kondratev 1972). However, water vapour is a passive component in the troposphere as it is uniquely determined by temperature and should therefore be seen as a part of the climate feedback system. In this short overview, we will first describe the water on planet Earth and the role of the hydrological cycle: the way water vapour is transported between oceans and continents and the return of water via rivers to the oceans. Generally water vapour is well observed and analysed; however, there are considerable obstacles to observing precipitation, in particular over the oceans. The response of the hydrological cycle to global warming is far reaching. Because different physical processes control the change in water vapour and evaporation/precipitation, this leads to a more extreme distribution of precipitation making, in general, wet areas wetter and dry areas dryer. Another consequence is a transition towards more intense precipitation. It is to be expected that the changes in the hydrological cycle as a consequence of climate warming may be more severe that the temperature changes.

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].

Ionospheric measurements of relative coronal brightness during the total solar eclipses of 11 August, 1999 and 9 July, 1945

Swept-frequency (1-10 MHz) ionosonde measurements were made at Helston, Cornwall (50 degrees 06'N, 5 degrees 18'W) during the total solar eclipse on August 11, 1999. Soundings were made every three minutes. We present a method for estimating the percentage of the ionising solar radiation which remains unobscured at any time during the eclipse by comparing the variation of the ionospheric E-layer with the behaviour of the layer during a control day. Application to the ionosonde date for II August, 1999, shows that the flux of solar ionising radiation fell to a minimum of 25 +/- 2% of the value before and after the eclipse. For comparison, the same technique was also applied to measurements made during the total solar eclipse of 9 July, 1945, at Sormjole (63 degrees 68'N, 20 degrees 20'E) and yielded a corresponding minimum of 16 +/- 2%. Therefore the method can detect variations in the fraction of solar emissions that originate from the unobscured corona and chromosphere. We discuss the differences between these two eclipses in terms of the nature of the eclipse, short-term fluctuations, the sunspot cycle and the recently-discovered long-term change in the coronal magnetic field.

Long-term drift of the coronal source magnetic flux and the total solar irradiance

We test the method of Lockwood et al. [1999] for deriving the coronal source flux from the geomagnetic aa index and show it to be accurate to within 12% for annual means and 4.5% for averages over a sunspot cycle. Using data from four solar constant monitors during 1981-1995, we find a linear relationship between this magnetic flux and the total solar irradiance. From this correlation, we show that the 131% rise in the mean coronal source field over the interval 1901-1995 corresponds to a rise in the average total solar irradiance of {\Delta}I = 1.65 +/- 0.23 Wm^{-2}.