Solar origin of heliospheric magnetic field inversions: evidence for coronal loop opening within pseudostreamers

The orientation of the heliospheric magnetic field (HMF) in near‒Earth space is generally a good indicator of the polarity of HMF foot points at the photosphere. There are times, however, when the HMF folds back on itself (is inverted), as indicated by suprathermal electrons locally moving sunward, even though they must ultimately be carrying the heat flux away from the Sun. Analysis of the near‒Earth solar wind during the period 1998–2011 reveals that inverted HMF is present approximately 5.5% of the time and is generally associated with slow, dense solar wind and relatively weak HMF intensity. Inverted HMF is mapped to the coronal source surface, where a new method is used to estimate coronal structure from the potential‒field source‒surface model. We find a strong association with bipolar streamers containing the heliospheric current sheet, as expected, but also with unipolar or pseudostreamers, which contain no current sheet. Because large‒scale inverted HMF is a widely accepted signature of interchange reconnection at the Sun, this finding provides strong evidence for models of the slow solar wind which involve coronal loop opening by reconnection within pseudostreamer belts as well as the bipolar streamer belt. Occurrence rates of bipolar‒ and pseudostreamers suggest that they are equally likely to result in inverted HMF and, therefore, presumably undergo interchange reconnection at approximately the same rate. Given the different magnetic topologies involved, this suggests the rate of reconnection is set externally, possibly by the differential rotation rate which governs the circulation of open solar flux.

Some results from an intercomparison of the climates simulated by 14 atmospheric general circulation models

Some climatological information from 14 atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind, and meridional stream function together with latitudinal distributions of mean sea level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated; the general improvement in the ability of models to simulate certain aspects of the climate is shown; consideration is given to the effect of increasing resolution on simulated climate; and approaches to understanding and reducing model deficiencies are discussed. The information presented here is a subset of a more voluminous compilation which is available in report form (Boer et al., 1991). This report contains essentially the same text, but results from all 14 models are presented together with additional results in the form of geographical distributions of surface variables and certain difference statistics.

An intercomparison of the climates simulated by 14 atmospheric general circulation models

Climatological information from fourteen atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind and meridional stream function together with latitudinal distributions of mean sea-level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated, general improvement in the ability of models to simulate certain aspects of the climate is shown, consideration is given to the effect of increasing resolution on simulated climate and approaches to the understanding and reduction of model deficiencies are discussed.

Method of determining a geophysical-scale sea ice rheology from laboratory measurements

We present a methodology that allows a sea ice rheology, suitable for use in a General Circulation Model (GCM), to be determined from laboratory and tank experiments on sea ice when combined with a kinematic model of deformation. The laboratory experiments determine a material rheology for sea ice, and would investigate a nonlinear friction law of the form τ ∝ σ n⅔, instead of the more familiar Amonton's law, τ = μσn (τ is the shear stress, μ is the coefficient of friction and σ n is the normal stress). The modelling approach considers a representative region R containing ice floes (or floe aggregates), separated by flaws. The deformation of R is imposed and the motion of the floes determined using a kinematic model, which will be motivated from SAR observations. Deformation of the flaws is inferred from the floe motion and stress determined from the material rheology. The stress over R is then determined from the area-weighted contribution from flaws and floes

Representation of model error in a convective-scale ensemble prediction system

In this paper ensembles of forecasts (of up to six hours) are studied from a convection-permitting model with a representation of model error due to unresolved processes. The ensemble prediction system (EPS) used is an experimental convection-permitting version of the UK Met Office’s 24- member Global and Regional Ensemble Prediction System (MOGREPS). The method of representing model error variability, which perturbs parameters within the model’s parameterisation schemes, has been modified and we investigate the impact of applying this scheme in different ways. These are: a control ensemble where all ensemble members have the same parameter values; an ensemble where the parameters are different between members, but fixed in time; and ensembles where the parameters are updated randomly every 30 or 60 min. The choice of parameters and their ranges of variability have been determined from expert opinion and parameter sensitivity tests. A case of frontal rain over the southern UK has been chosen, which has a multi-banded rainfall structure. The consequences of including model error variability in the case studied are mixed and are summarised as follows. The multiple banding, evident in the radar, is not captured for any single member. However, the single band is positioned in some members where a secondary band is present in the radar. This is found for all ensembles studied. Adding model error variability with fixed parameters in time does increase the ensemble spread for near-surface variables like wind and temperature, but can actually decrease the spread of the rainfall. Perturbing the parameters periodically throughout the forecast does not further increase the spread and exhibits “jumpiness” in the spread at times when the parameters are perturbed. Adding model error variability gives an improvement in forecast skill after the first 2–3 h of the forecast for near-surface temperature and relative humidity. For precipitation skill scores, adding model error variability has the effect of improving the skill in the first 1–2 h of the forecast, but then of reducing the skill after that. Complementary experiments were performed where the only difference between members was the set of parameter values (i.e. no initial condition variability). The resulting spread was found to be significantly less than the spread from initial condition variability alone.

The physics of orographic gravity wave drag

The drag and momentum fluxes produced by gravity waves generated in flow over orography are reviewed, focusing on adiabatic conditions without phase transitions or radiation effects, and steady mean incoming flow. The orographic gravity wave drag is first introduced in its simplest possible form, for inviscid, linearized, non-rotating flow with the Boussinesq and hydrostatic approximations, and constant wind and static stability. Subsequently, the contributions made by previous authors (primarily using theory and numerical simulations) to elucidate how the drag is affected by additional physical processes are surveyed. These include the effect of orography anisotropy, vertical wind shear, total and partial critical levels, vertical wave reflection and resonance, non-hydrostatic effects and trapped lee waves, rotation and nonlinearity. Frictional and boundary layer effects are also briefly mentioned. A better understanding of all of these aspects is important for guiding the improvement of drag parametrization schemes.

Can dynamically downscaled windstorm footprints be improved by observations through a probabilistic approach?

Windstorms are a main feature of the European climate and exert strong socioeconomic impacts. Large effort has been made in developing and enhancing models to simulate the intensification of windstorms, resulting footprints, and associated impacts. Simulated wind or gust speeds usually differ from observations, as regional climate models have biases and cannot capture all local effects. An approach to adjust regional climate model (RCM) simulations of wind and wind gust toward observations is introduced. For this purpose, 100 windstorms are selected and observations of 173 (111) test sites of the German Weather Service are considered for wind (gust) speed. Theoretical Weibull distributions are fitted to observed and simulated wind and gust speeds, and the distribution parameters of the observations are interpolated onto the RCM computational grid. A probability mapping approach is applied to relate the distributions and to correct the modeled footprints. The results are not only achieved for single test sites but for an area-wide regular grid. The approach is validated using root-mean-square errors on event and site basis, documenting that the method is generally able to adjust the RCM output toward observations. For gust speeds, an improvement on 88 of 100 events and at about 64% of the test sites is reached. For wind, 99 of 100 improved events and ~84% improved sites can be obtained. This gives confidence on the potential of the introduced approach for many applications, in particular those considering wind data.

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.

The low-latitude boundary layer: Application of ISTP advances to past data

The destruction of the four Cluster craft was a major loss to the planned ISTP effort, of which studies of the magnetopause and low-latitude boundary layer (LLBL) were an important part. While awaiting the re-flight mission, Cluster-II, we have been applying advances in our understanding made using other ISTP craft (like Polar and Wind) and using ground-based facilities (in particular the EISCAT incoherent scatter radars and the SuperDARN HF coherent radars) to measurements of the LLBL made in 1984 and 1985 by the AMPTE-UKS and -IRM spacecraft pair. In particular, one unexplained result of the AMPTE mission was that the electron characteristics could, in nearly all cases, order independent measurements near the magnetopause, such as the magnetic field, ion temperatures and the plasma flow. Studies of the cusp have shown that the precipitation is ordered by the time-elapsed since the field line was opened by reconnection. This insight has allowed us to reanalyse the AMPTE data and show that the ordering by the transition parameter is also due to the variation of time elapsed since reconnection, with the important implication that reconnection usually coats most of the dayside magnetopause with at least some newly-opened field lines. In addition, we can use the electron characteristics to isolate features like RDs, slow-mode shocks and slow-mode expansion fans. The ion characteristics can be used to compute the reconnection rate. We here retrospectively apply these new techniques, developed in the ISTP era, to a much-studied flux transfer event observed by the AMPTE satellites. As a result, we gain new understanding of its cause and structure.

Ionospheric convection response to slow, strong variations in a northward interplanetary magnetic field: A case study for January 14, 1988

We analyze ionospheric convection patterns over the polar regions during the passage of an interplanetary magnetic cloud on January 14, 1988, when the interplanetary magnetic field (IMF) rotated slowly in direction and had a large amplitude. Using the assimilative mapping of ionospheric electrodynamics (AMIE) procedure, we combine simultaneous observations of ionospheric drifts and magnetic perturbations from many different instruments into consistent patterns of high-latitude electrodynamics, focusing on the period of northward IMF. By combining satellite data with ground-based observations, we have generated one of the most comprehensive data sets yet assembled and used it to produce convection maps for both hemispheres. We present evidence that a lobe convection cell was embedded within normal merging convection during a period when the IMF By and Bz components were large and positive. As the IMF became predominantly northward, a strong reversed convection pattern (afternoon-to-morning potential drop of around 100 kV) appeared in the southern (summer) polar cap, while convection in the northern (winter) hemisphere became weak and disordered with a dawn-to-dusk potential drop of the order of 30 kV. These patterns persisted for about 3 hours, until the IMF rotated significantly toward the west. We interpret this behavior in terms of a recently proposed merging model for northward IMF under solstice conditions, for which lobe field lines from the hemisphere tilted toward the Sun (summer hemisphere) drape over the dayside magnetosphere, producing reverse convection in the summer hemisphere and impeding direct contact between the solar wind and field lines connected to the winter polar cap. The positive IMF Bx component present at this time could have contributed to the observed hemispheric asymmetry. Reverse convection in the summer hemisphere broke down rapidly after the ratio |By/Bz| exceeded unity, while convection in the winter hemisphere strengthened. A dominant dawn-to-dusk potential drop was established in both hemispheres when the magnitude of By exceeded that of Bz, with potential drops of the order of 100 kV, even while Bz remained northward. The later transition to southward Bz produced a gradual intensification of the convection, but a greater qualitative change occurred at the transition through |By/Bz| = 1 than at the transition through Bz = 0. The various convection patterns we derive under northward IMF conditions illustrate all possibilities previously discussed in the literature: nearly single-cell and multicell, distorted and symmetric, ordered and unordered, and sunward and antisunward.

Flux transfer events at the dayside magnetopause: transient reconnection or magnetosheath dynamic pressure pulses?

The suggestion is discussed that characteristic particle and field signatures at the dayside magnetopause, termed “flux transfer events” (FTEs), are, in at least some cases, due to transient solar wind and/or magnetosheath dynamic pressure increases, rather than time-dependent magnetic reconnection. It is found that most individual cases of FTEs observed by a single spacecraft can, at least qualitatively, be explained by the pressure pulse model, provided a few rather unsatisfactory features of the predictions are explained in terms of measurement uncertainties. The most notable exceptions to this are some “two-regime” observations made by two satellites simultaneously, one on either side of the magnetopause. However, this configuration has not been frequently achieved for sufficient time, such observations are rare, and the relevant tests are still not conclusive. The strongest evidence that FTEs are produced by magnetic reconnection is the dependence of their occurrence on the north-south component of the interplanetary magnetic field (IMF) or of the magnetosheath field. The pressure pulse model provides an explanation for this dependence (albeit qualitative) in the case of magnetosheath FTEs, but this does not apply to magnetosphere FTEs. The only surveys of magnetosphere FTEs have not employed the simultaneous IMF, but have shown that their occurrence is strongly dependent on the north-south component of the magnetosheath field, as observed earlier/later on the same magnetopause crossing (for inbound/outbound passes, respectively). This paper employs statistics on the variability of the IMF orientation to investigate the effects of IMF changes between the times of the magnetosheath and FTE observations. It is shown that the previously published results are consistent with magnetospheric FTEs being entirely absent when the magnetosheath field is northward: all crossings with magnetosphere FTEs and a northward field can be attributed to the field changing sense while the satellite was within the magnetosphere (but close enough to the magnetopause to detect an FTE). Allowance for the IMF variability also makes the occurrence frequency of magnetosphere FTEs during southward magnetosheath fields very similar to that observed for magnetosheath FTEs. Conversely, the probability of attaining the observed occurrence frequencies for the pressure pulse model is 10−14. In addition, it is argued that some magnetosheath FTEs should, for the pressure pulse model, have been observed for northward IMF: the probability that the number is as low as actually observed is estimated to be 10−10. It is concluded that although the pressure model can be invoked to qualitatively explain a large number of individual FTE observations, the observed occurrence statistics are in gross disagreement with this model.

Vertical ion flow in the polar ionosphere

Outflowing ions from the polar ionosphere fall into two categories: the classical polar wind and the suprathermal ion flows. The flows in both these categories vary a great deal with altitude. The classical polar wind is supersonic at high altitude: at ∼3 RE geocentric, the observed polar wind is H+ dominated and has a Mach number of 2.5–5.1. At 400–600 km, thermal and suprathermal upward O+ ion fluxes frequently occur at the poleward edge of the nightside auroral oval during magnetically active times. Above 500 km, ions are accelerated transverse to the local geomagnetic field. At 1400 km, transversely accelerated ions are frequently observed in winter nights but rarely appear in the summer. In the dayside cleft above ∼2000 km, ions of all species are transversely heated and upwell with significant number and heat fluxes, forming a cleft ion fountain as they convect across the polar cap. Upwelling ions are observed most (least) frequently in the summer (winter). At yet higher altitudes, energetic (>10 eV to several kiloelectron volts) upflowing H+ and O+ ions are frequently observed, their active time occurrence frequency being as high as 0.7 at auroral latitudes and 0.3 in the polar cap. Their composition, intensity, and angular characteristics vary quantitatively with solar activity, being O+ dominant and more intense near solar maximum. Their resulting ion outflow is dominated by ions below 1 keV and reaches 3.5×10^26 O+ and 7×10^25 H+ ions s^{−1} at magnetically active times (Kp≥5) near solar maximum. In comparison, the estimated polar wind ion outflow at times of moderate solar activity is 7×10^25H+ and 4×10^24 He+ ions s^{−1}. The estimated <10-eV cleft ion fountain flow is 3.8×10^25 O+ and 8.6×10^23 H+ ions s^{−1} near solar maximum.

Low-energy ion flows into the magnetosphere

Ion flows from the ionosphere into the magnetosphere fall into two main categories: cold (<1eV), “classical” polar wind and heated (>1eV), suprathermal ion outflows. A wealth of new understanding of these outflows has resulted from the Dynamics Explorer Mission. This review describes both the confirmation of the predicted classical polar wind as well as the revelation of a great variety of low-energy suprathermal outflows: the cleft ion fountain, the nightside auroral fountaion (X-events, toroids and field-aligned flows) and polar cap outflows. The main emphasis is placed on flows at energies below about 50eV, observed by the Retarding Ion Mass Spectrometer (RIMS) on board the Dynamics Explorer 1 satellite; limited comparisons are made with results from other instruments which sample different energy ranges.

A new source of suprathermal O+ions near the dayside polar cap boundary

A new dayside source of O+ ions for the polar magnetosphere is described, and a statistical survey presented of upward flows of O+ ions using 2 years of data from the retarding ion mass spectrometer (RIMS) experiment on board DE 1, at geocentric distances below 3 RE and invariant latitudes above 40°. The flows are classified according to their spin angle distributions. It is believed that the spacecraft potential near perigee is generally less than +2 V, in which case the entire O+ population at energies below about 60 eV is sampled. Examples are given of field-aligned flow and of transversely accelerated “core” O+ ions; in the latter events a large fraction of the total O+ ion population has been transversely accelerated, and in some extreme cases all the observed ions (of all ion species) have been accelerated, and no residual cold population is observed (“toroidal” distributions). However, by far the most common type of O+ upflow seen by DE RIMS lies near the dayside polar cap boundary (particularly in the prenoon sector) and displays an asymmetric spin angle distribution. In such events the ions carry an upward heat flux, and strong upflow of all species is present (H+, He+, O+, O++, and N+ have all been observed with energies up to about 30 eV, but with the majority of ions below about 2 eV); hence, these have been termed upwelling ion events. The upwelling ions are embedded in larger regions of classical light ion polar wind and are persistently found under the following conditions: at geocentric distances greater than 1.4 RE; at all Kp in summer, but only at high Kp in winter. Low-energy conical ions (<30 eV) are only found near the equatorial edge of the events, the latitude of which moves equatorward with increasing Kp and is highly correlated with the location of field-aligned currents. The RIMS data are fully consistent with a “mass spectrometer effect,” whereby light ions and the more energetic O+ ions flow into the lobes and mantle and hence the far-tail plasma sheet, but lower-energy O+ is swept across the polar cap by the convection electric field, potentially acting as a source for the nightside auroral acceleration regions. The occurrence probability of upwelling ion events, as compared to those of low-altitude transversely accelerated core ions and of field-aligned flow, suggests this could be the dominant mechanism for supplying the nightside auroral acceleration region, and subsequently the ring current and near-earth plasma sheet, with ionospheric O+ ions. It is shown that the total rate of O+ outflow in upwelling ion events (greater than 10^25 s^{−1}) is sufficient for the region near the dayside polar cap boundary to be an important ionospheric heavy ion source.

Departures from diffusive equilibrium in the topside F-layer from satellite soundings

A method for the detection of O+ ion fluxes from topside soundings is described. The shape of the plasma scale-height profile is altered by such flows only at heights near the F2-peak, where ion-neutral drag is large. Model profiles are used to relate changes in scale height to the ratio (φ/φL) where φ is the field-aligned O+ flux (relative to the neutral air) and φL is the limiting value set by frictional drag. Values of (φ/φL) can then be determined to within a few per cent from experimental soundings, using the plasma temperature and its gradient (as deduced from the observed profile) and the MSIS model neutral temperature. It was found that 3700 topside profiles show departures from diffusive equilibrium, out of 10,000 used to obtain the global morphology of (φ/φL) near the sunspot minimum. Results reveal dynamic ion-flow effects such as the transequatorial breeze and the effects of the polar wind and protonospheric replenishment light-ion flows can be inferred.