Magnetic discontinuities in the near-earth solar wind: evidence of in-transit turbulence or remnants of coronal structure?

Fluctuations in the solar wind plasma and magnetic field are well described by the sum of two power law distributions. It has been postulated that these distributions are the result of two independent processes: turbulence, which contributes mainly to the smaller fluctuations, and crossing the boundaries of flux tubes of coronal origin, which dominates the larger variations. In this study we explore the correspondence between changes in the magnetic field with changes in other solar wind properties. Changes in density and temperature may result from either turbulence or coronal structures, whereas changes in composition, such as the alpha-to-proton ratio are unlikely to arise from in-transit effects. Observations spanning the entire ACE dataset are compared with a null hypothesis of no correlation between magnetic field discontinuities and changes in other solar wind parameters. Evidence for coronal structuring is weaker than for in-transit turbulence, with only ∼ 25% of large magnetic field discontinuities associated with a significant change in the alpha-to-proton ratio, compared to ∼ 40% for significant density and temperature changes. However, note that a lack of detectable alpha-to-proton signature is not sufficient to discount a structure as having a solar origin.

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.

A 27 day persistence model of near-Earth solar wind conditions: a long lead-time forecast and a benchmark for dynamical models

Geomagnetic activity has long been known to exhibit approximately 27 day periodicity, resulting from solar wind structures repeating each solar rotation. Thus a very simple near-Earth solar wind forecast is 27 day persistence, wherein the near-Earth solar wind conditions today are assumed to be identical to those 27 days previously. Effective use of such a persistence model as a forecast tool, however, requires the performance and uncertainty to be fully characterized. The first half of this study determines which solar wind parameters can be reliably forecast by persistence and how the forecast skill varies with the solar cycle. The second half of the study shows how persistence can provide a useful benchmark for more sophisticated forecast schemes, namely physics-based numerical models. Point-by-point assessment methods, such as correlation and mean-square error, find persistence skill comparable to numerical models during solar minimum, despite the 27 day lead time of persistence forecasts, versus 2–5 days for numerical schemes. At solar maximum, however, the dynamic nature of the corona means 27 day persistence is no longer a good approximation and skill scores suggest persistence is out-performed by numerical models for almost all solar wind parameters. But point-by-point assessment techniques are not always a reliable indicator of usefulness as a forecast tool. An event-based assessment method, which focusses key solar wind structures, finds persistence to be the most valuable forecast throughout the solar cycle. This reiterates the fact that the means of assessing the “best” forecast model must be specifically tailored to its intended use.

Ionospheric signatures of pulsed reconnection at the Earth's magnetopause

THE plasma precipitating into the Earth's dayside auroral atmosphere has characteristics which show that it originates from the shocked solar-wind plasma of the magnetosheath1'2. The particles of the magnetosheath plasma precipitate down a funnel-shaped region (cusp) of open field lines resulting from reconnection of the geomagnetic field with the interplanetary magnetic field3. Although the cusp has long been considered a well defined spatial structure maintained by continuous reconnection, it has recently been suggested4–6 that reconnection instead may take place in a series of discontinuous events; this is the ‘pulsating cusp model’. Here we present coordinated radar and satellite observations of a series of discrete, poleward-moving plasma structures that are consistent with the pulsating-cusp model.

Near-earth cosmic ray decreases associated with remote coronal mass ejections

Galactic cosmic ray (GCR) flux is modulated by both particle drift patterns and solar wind structures on a range of timescales. Over solar cycles, GCR flux varies as a function of the total open solar magnetic flux and the latitudinal extent of the heliospheric current sheet. Over hours, drops of a few percent in near-Earth GCR flux (Forbush decreases, FDs) are well known to be associated with the near-Earth passage of solar wind structures resulting from corotating interaction regions (CIRs) and transient coronal mass ejections (CMEs). We report on four FDs seen at ground-based neutron monitors which cannot be immediately associated with significant structures in the local solar wind. Similarly, there are significant near-Earth structures which do not produce any corresponding GCR variation. Three of the FDs are during the STEREO era, enabling in situ and remote observations from three well-separated heliospheric locations. Extremely large CMEs passed the STEREO-A spacecraft, which was behind the West limb of the Sun, approximately 2–3 days before each near- Earth FD. Solar wind simulations suggest that the CMEs combined with pre-existing CIRs, enhancing the pre-existing barriers to GCR propagation. Thus these observations provide strong evidence for the modulation of GCR flux by remote solar wind structures.

Convectively coupled equatorial waves: A new methodology for identifying wave structures in observational data

Convectively coupled equatorial waves are fundamental components of the interaction between the physics and dynamics of the tropical atmosphere. A new methodology, which isolates individual equatorial wave modes, has been developed and applied to observational data. The methodology assumes that the horizontal structures given by equatorial wave theory can be used to project upper- and lower-tropospheric data onto equatorial wave modes. The dynamical fields are first separated into eastward- and westward-moving components with a specified domain of frequency–zonal wavenumber. Each of the components for each field is then projected onto the different equatorial modes using the y structures of these modes given by the theory. The latitudinal scale yo of the modes is predetermined by data to fit the equatorial trapping in a suitable latitude belt y = ±Y. The extent to which the different dynamical fields are consistent with one another in their depiction of each equatorial wave structure determines the confidence in the reality of that structure. Comparison of the analyzed modes with the eastward- and westward-moving components in the convection field enables the identification of the dynamical structure and nature of convectively coupled equatorial waves. In a case study, the methodology is applied to two independent data sources, ECMWF Reanalysis and satellite-observed window brightness temperature (Tb) data for the summer of 1992. Various convectively coupled equatorial Kelvin, mixed Rossby–gravity, and Rossby waves have been detected. The results indicate a robust consistency between the two independent data sources. Different vertical structures for different wave modes and a significant Doppler shifting effect of the background zonal winds on wave structures are found and discussed. It is found that in addition to low-level convergence, anomalous fluxes induced by strong equatorial zonal winds associated with equatorial waves are important for inducing equatorial convection. There is evidence that equatorial convection associated with Rossby waves leads to a change in structure involving a horizontal structure similar to that of a Kelvin wave moving westward with it. The vertical structure may also be radically changed. The analysis method should make a very powerful diagnostic tool for investigating convectively coupled equatorial waves and the interaction of equatorial dynamics and physics in the real atmosphere. The results from application of the analysis method for a reanalysis dataset should provide a benchmark against which model studies can be compared.

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.

Evaluation of the Met Office global forecast model using Geostationary Earth Radiation Budget (GERB) data

Simulations of the top-of-atmosphere radiative-energy budget from the Met Office global numerical weather-prediction model are evaluated using new data from the Geostationary Earth Radiation Budget (GERB) instrument on board the Meteosat-8 satellite. Systematic discrepancies between the model simulations and GERB measurements greater than 20 Wm-2 in outgoing long-wave radiation (OLR) and greater than 60 Wm-2 in reflected short-wave radiation (RSR) are identified over the period April-September 2006 using 12 UTC data. Convective cloud over equatorial Africa is spatially less organized and less reflective than in the GERB data. This bias depends strongly on convective-cloud cover, which is highly sensitive to changes in the model convective parametrization. Underestimates in model OLR over the Gulf of Guinea coincide with unrealistic southerly cloud outflow from convective centres to the north. Large overestimates in model RSR over the subtropical ocean, greater than 50 Wm-2 at 12 UTC, are explained by unrealistic radiative properties of low-level cloud relating to overestimation of cloud liquid water compared with independent satellite measurements. The results of this analysis contribute to the development and improvement of parametrizations in the global forecast model.

Water wave scattering by finite arrays of circular structures

The scattering of small amplitude water waves by a finite array of locally axisymmetric structures is considered. Regions of varying quiescent depth are included and their axisymmetric nature, together with a mild-slope approximation, permits an adaptation of well-known interaction theory which ultimately reduces the problem to a simple numerical calculation. Numerical results are given and effects due to regions of varying depth on wave loading and free-surface elevation are presented.

Data assimilation: From photon counts to Earth System forecasts

Data assimilation – the set of techniques whereby information from observing systems and models is combined optimally – is rapidly becoming prominent in endeavours to exploit Earth Observation for Earth sciences, including climate prediction. This paper explains the broad principles of data assimilation, outlining different approaches (optimal interpolation, three-dimensional and four-dimensional variational methods, the Kalman Filter), together with the approximations that are often necessary to make them practicable. After pointing out a variety of benefits of data assimilation, the paper then outlines some practical applications of the exploitation of Earth Observation by data assimilation in the areas of operational oceanography, chemical weather forecasting and carbon cycle modelling. Finally, some challenges for the future are noted.

Preferred structures in large-scale circulation and the effect of doubling greenhouse gas concentration in HadCM3

Preferred structures in the surface pressure variability are investigated in and compared between two 100-year simulations of the Hadley Centre climate model HadCM3. In the first (control) simulation, the model is forced with pre-industrial carbon dioxide concentration (1×CO2) and in the second simulation the model is forced with doubled CO2 concentration (2×CO2). Daily winter (December-January-February) surface pressures over the Northern Hemisphere are analysed. The identification of preferred patterns is addressed using multivariate mixture models. For the control simulation, two significant flow regimes are obtained at 5% and 2.5% significance levels within the state space spanned by the leading two principal components. They show a high pressure centre over the North Pacific/Aleutian Islands associated with a low pressure centre over the North Atlantic, and its reverse. For the 2×CO2 simulation, no such behaviour is obtained. At higher-dimensional state space, flow patterns are obtained from both simulations. They are found to be significant at the 1% level for the control simulation and at the 2.5% level for the 2×CO2 simulation. Hence under CO2 doubling, regime behaviour in the large-scale wave dynamics weakens. Doubling greenhouse gas concentration affects both the frequency of occurrence of regimes and also the pattern structures. The less frequent regime becomes amplified and the more frequent regime weakens. The largest change is observed over the Pacific where a significant deepening of the Aleutian low is obtained under CO2 doubling.