Ensemble downscaling in coupled solar wind-magnetosphere modeling for space weather forecasting

Advanced forecasting of space weather requires simulation of the whole Sun-to-Earth system, which necessitates driving magnetospheric models with the outputs from solar wind models. This presents a fundamental difficulty, as the magnetosphere is sensitive to both large-scale solar wind structures, which can be captured by solar wind models, and small-scale solar wind “noise,” which is far below typical solar wind model resolution and results primarily from stochastic processes. Following similar approaches in terrestrial climate modeling, we propose statistical “downscaling” of solar wind model results prior to their use as input to a magnetospheric model. As magnetospheric response can be highly nonlinear, this is preferable to downscaling the results of magnetospheric modeling. To demonstrate the benefit of this approach, we first approximate solar wind model output by smoothing solar wind observations with an 8 h filter, then add small-scale structure back in through the addition of random noise with the observed spectral characteristics. Here we use a very simple parameterization of noise based upon the observed probability distribution functions of solar wind parameters, but more sophisticated methods will be developed in the future. An ensemble of results from the simple downscaling scheme are tested using a model-independent method and shown to add value to the magnetospheric forecast, both improving the best estimate and quantifying the uncertainty. We suggest a number of features desirable in an operational solar wind downscaling scheme.

The first real-time worldwide ionospheric predictions network: an advance in support of spaceborne experimentation, on-line model validation, and space weather

We report on the first realtime ionospheric predictions network and its capabilities to ingest a global database and forecast F-layer characteristics and "in situ" electron densities along the track of an orbiting spacecraft. A global network of ionosonde stations reported around-the-clock observations of F-region heights and densities, and an on-line library of models provided forecasting capabilities. Each model was tested against the incoming data; relative accuracies were intercompared to determine the best overall fit to the prevailing conditions; and the best-fit model was used to predict ionospheric conditions on an orbit-to-orbit basis for the 12-hour period following a twice-daily model test and validation procedure. It was found that the best-fit model often provided averaged (i.e., climatologically-based) accuracies better than 5% in predicting the heights and critical frequencies of the F-region peaks in the latitudinal domain of the TSS-1R flight path. There was a sharp contrast however, in model-measurement comparisons involving predictions of actual, unaveraged, along-track densities at the 295 km orbital altitude of TSS-1R In this case, extrema in the first-principle models varied by as much as an order of magnitude in density predictions, and the best-fit models were found to disagree with the "in situ" observations of Ne by as much as 140%. The discrepancies are interpreted as a manifestation of difficulties in accurately and self-consistently modeling the external controls of solar and magnetospheric inputs and the spatial and temporal variabilities in electric fields, thermospheric winds, plasmaspheric fluxes, and chemistry.

Incoherent scatter radar observations related to magnetospheric dynamics

Over the past decade incoherent scatter radars have provided fundamental observations of velocities and plasma parameters in the high-latitude ionosphere which relate to the dynamical processes responsible for the excitation of flow in the coupled solar wind-magnetosphere-ionosphere system. These observations have played a central role in inspiring a change of paradigm from a picture of quasi-steady flows parameterised by the direction of the interplanetary magnetic field to a picture of inherently time-dependent flows driven by coupling processes at the magnetopause and in the tail. Flows and particle precipitation in the dayside ionosphere are reasonably well understood in principle in terms of the effects of time-dependent reconnection at the magnetopause, though coordinated high- and low-altitude observations are lacking. Related phenomena also appear to occur in the tail, forming the “equatorward-drifting arcs” which are present during quiet times, as well as during the growth and early expansion phases of substorms. At expansion onset, the substorm bulge forms well equatorward of the arc formation region, and may take ∼ 10 min or more to reach it in its poleward expansion. Nightside ionospheric flows are then considerably perturbed by the effects of strong precipitation-induced conductivity gradients.

Testing substorm theories: The need for multipoint observations

The requirement for multipoint observations to test theories of magnetospheric substorms is reviewed. A wide variety of such theories have been proposed, but these cannot be properly evaluated because we do not understand how the various features of a substorm are causally linked. In terms of explaining certain substorm features, some theories may be mutually-exclusive rivals. But this is not always the case, making it possible that theories may be either combined into a synthesis model or loosely connected in a more modular view of substorms. Some key questions are defined which require multipoint in-situ measurements, combined with remote sensing observations, of the development and relationship of the major substorm features.

Opportunities for magnetospheric research with coordinated cluster and ground-based observations

ESA’s first multi-satellite mission Cluster is unique in its concept of 4 satellites orbiting in controlled formations. This will give an unprecedented opportunity to study structure and dynamics of the magnetosphere. In this paper we discuss ways in which ground-based remote-sensing observations of the ionosphere can be used to support the multipoint in-situ satellite measurements. There are a very large number of potentially useful configurations between the satellites and any one ground-based observatory; however, the number of ideal occurrences for any one configuration is low. Many of the ground-based instruments cannot operate continuously and Cluster will take data only for a part of each orbit, depending on how much high-resolution (‘burst-mode’) data are acquired. In addition, there are a great many instrument modes and the formation, size and shape of the cluster of the four satellites to consider. These circumstances create a clear and pressing need for careful planning to ensure that the scientific return from Cluster is maximised by additional coordinated ground-based observations. For this reason, ESA established a working group to coordinate the observations on the ground with Cluster. We will give a number of examples how the combined spacecraft and ground-based observations can address outstanding questions in magnetospheric physics. An online computer tool has been prepared to allow for the planning of conjunctions and advantageous constellations between the Cluster spacecraft and individual or combined ground-based systems. During the mission a ground-based database containing index and summary data will help to identify interesting datasets and allow to select intervals for coordinated studies. We illustrate the philosophy of our approach, using a few important examples of the many possible configurations between the satellite and the ground-based instruments.

Radar observations of non-thermal plasmas at different aspect angles

Data are presented from the EISCAT CP-3-E experiment which show the presence of non-thermal plasma over a range of latitudes. The O+ ion-velocity distribution function is almost toroidal when the electric field reaches values of 125 mV m−1. The ion temperature derived from such data assuming a Maxwellian distribution function will overestimate the true ion temperature when the observing angle is large with respect to the magnetic field, and underestimate the temperature when the aspect angle is small. When the expressions for the distribution function are extended to include mixed ion composition, an improvement is sometimes found in fitting the observed data, and estimates of the composition can be made. Such an analysis suggests that N2+ can occasionally form a significant part of the total ion density in a narrow height region centred at 275 km.

Incoherent scatter radar observations of non-Maxwellian ion velocity distributions in the auroral F-region

Observations by the EISCAT experiments “POLAR” and Common Programme CP-3 reveal non-Maxwellian ion velocity distributions in the auroral F-region ionosphere. Analysis of data from three periods is presented. During the first period, convection velocities are large (≈2 km s-1) and constant over part of a CP-3 latitude scan; the second period is one of POLAR data containing a short-lived (<1 min.) burst of rapid (>1.5 km s-1) flow. We concentrate on these two periods as they allow the study of a great many features of the ion-neutral interactions which drive the plasma non-thermal and provide the best available experimental test for models of the 3-dimensional ion velocity distribution function. The third period is included to illustrate the fact that non-thermal plasma frequently exists in the auroral ionosphere: the data, also from the POLAR experiment, cover a three-hour period of typical auroral zone flow and analysis reveals that the ion distribution varies from Maxwellian to the threshold of a toroidal form.

Ion flows and heating at a contracting polar-cap boundary: GISMOS observations indicating viscous-like interaction on the flanks of the magnetotail

This paper complements that in this issue by Clauer et al. concerning the international GISMOS campaign of 3–5 June 1987. From a detailed study of the EISCAT data, the polar-cap boundary, as defined by an almost shear east-west convection reversal, is found to contract across the EISCAT field of view between 04 and 07 MLT. An annulus of enhanced ion temperature and non-thermal plasma is observed immediately equatorward of the contracting boundary due to the lag in the response of the neutral-wind pattern to the change in ion flows. The ion flow inside the polar cap and at the boundary is shown to be relatively smooth, compared with that in the auroral oval, at 15-second resolution. The flow at the boundary is directed poleward, with velocities which exceed that of the boundary itself. The effect of velocity shears on the beamswinging technique used to derive the ion flows has been analysed in detail and it is found that spurious flows across a moving boundary can be generated. However, these are much smaller than the observed flows into the polar cap and cannot explain the 7 kV potential difference across the observed segment of the cap boundary between 04:30–06:30 UT. The ion temperature enhancements at the two observing azimuths is used to define the boundary orientation. The results are consistent with recent observations of slow anti-sunward flow of closed field lines on the flanks of the geomagnetic tail, which appears to be generated by some form of “viscous” coupling to the magnetosheath plasma.

Observations of large field-aligned flows of thermal plasma in the auroral ionosphere

Large upward field-aligned ion flows have previously been observed in the high latitude ionosphere in response to frictional heating of the local ion population. Results from a recent experiment using the EISCAT radar show similar features but allow, for the first time, determination of the field-aligned profiles of plasma parameters during these events. The upflows occur during frictional heating. The flows are shown to be transient plasma upwellings, from regions where the ion temperature has been elevated by the motion of a convection shear over the observed field line.

Observations at the magnetopause and in the auroral ionosphere of momentum transfer from the solar wind

Recent radar studies of field-perpendicular flows in the auroral ionosphere, in conjunction with observations of the interplanetary medium immediately upstream of the Earth's bow shock, have revealed direct control of dayside convection by the Bz component of the interplanetary magnetic field (IMF). The ionospheric flows begin to respond to both northward and southward turnings of the IMF impinging upon the magnetopause after a delay of only a few minutes in the early afternoon sector, rising to about 15 minutes nearer dawn and dusk. In both the polar cap and the auroral oval, the subsequent rise and decay times are of order 5–10 minutes. We conclude there is very little convection “flywheel” effect in the dayside polar ionosphere and that only newly-opened flux tubes impart significant momentum to the ionosphere, in a relatively narrow region immediately poleward of the cusp. These findings concerning the effects of quasi-steady reconnection have important implications for any ionospheric signatures of transient reconnection which should be considerably shorter-lived than thought hitherto. In order to demonstrate the difficulty of uniquely identifying a Flux Transfer Event (FTE) in ground-based magnetometer data, we present observations of an impulsive signature, identical with that expected for an FTE if data from only one station is studied, following an observed magnetopause compression when the IMF was purely northward. We also report new radar observations of a viscous-like interaction, consistent with an origin on the flanks of the magnetotail and contributing an estimated 15–30kV to the total cross-cap potential during quiet periods.

Direct observations of the full Dungey convection cycle in the polar ionosphere for southward interplanetary magnetic field conditions

Tracking the formation and full evolution of polar cap ionization patches in the polar ionosphere, we directly observe the full Dungey convection cycle for southward interplanetary magnetic field (IMF) conditions. This enables us to study how the Dungey cycle influences the patches’ evolution. The patches were initially segmented from the dayside storm enhanced density plume at the equatorward edge of the cusp, by the expansion and contraction of the polar cap boundary due to pulsed dayside magnetopause reconnection, as indicated by in situ Time History of Events and Macroscale Interactions during Substorms(THEMIS) observations. Convection led to the patches entering the polar cap and being transported antisunward, while being continuously monitored by the globally distributed arrays of GPS receivers and Super Dual Auroral Radar Network radars. Changes in convection over time resulted in the patches following a range of trajectories, each of which differed somewhat from the classical twin-cell convection streamlines. Pulsed nightside reconnection, occurring as part of the magnetospheric substorm cycle, modulated the exit of the patches from the polar cap, as confirmed by coordinated observations of the magnetometer at Tromsø and European Incoherent Scatter Tromsø UHF radar. After exiting the polar cap, the patches broke up into a number of plasma blobs and returned sunward in the auroral return flow of the dawn and/or dusk convection cell. The full circulation time was about 3 h.

Validation of a priori CME arrival predictions made using real-time heliospheric imager observations

Between December 2010 and March 2013, volunteers for the Solar Stormwatch (SSW) Citizen Science project have identified and analyzed coronal mass ejections (CMEs) in the near real-time Solar Terrestrial Relations Observatory Heliospheric Imager observations, in order to make “Fearless Forecasts” of CME arrival times and speeds at Earth. Of the 60 predictions of Earth-directed CMEs, 20 resulted in an identifiable Interplanetary CME (ICME) at Earth within 1.5–6 days, with an average error in predicted transit time of 22 h, and average transit time of 82.3 h. The average error in predicting arrival speed is 151 km s−1, with an average arrival speed of 425km s−1. In the same time period, there were 44 CMEs for which there are no corresponding SSW predictions, and there were 600 days on which there was neither a CME predicted nor observed. A number of metrics show that the SSW predictions do have useful forecast skill; however, there is still much room for improvement. We investigate potential improvements by using SSW inputs in three models of ICME propagation: two of constant acceleration and one of aerodynamic drag. We find that taking account of interplanetary acceleration can improve the average errors of transit time to 19 h and arrival speed to 77 km s−1.

In situ observations in cyclones over Fram Strait

During the international FRAMZY expedition in March 2002 in-situ observations of Fram Strait cyclones were made by aircraft, ship and automatic buoys in order to study the interaction between cyclones and sea ice. The atmospheric characteristics of the observed cyclones are presented in this paper. The cyclones were generated in the baroclinic zone at the ice edge and moved NNE-ward along the ice edge. This was supported by warm air advection from WSW by an upper-level wave. The cyclones were rather small (diameter 200– 700 km) and shallow (1–1.5 km e-folding height for the horizontal pressure and temperature difference) with life times between 12 and 36 hours. In spite of the small space and time scales, remarkable extremes were observed within the cyclones. Winds reached maxima above 20 ms−1 lasting for only a few hours. The transition from the cold to the advancing warm air over sea ice occurred within narrow (5–30 km) frontal zones in which vorticity and convergence reached maxima on the order of 10−3 s−1. It is discussed whether the sea ice in spite of its inertia is able to react on these strong sub cyclone-scale processes and, thus, these processes have to be taken into account in models in order to simulate the cyclone-sea ice interaction properly.

Understanding three-dimensional effects in polarized observations with the ground-based ADMIRARI radiometer during the CHUVA campaign

Measurements of down-welling microwave radiation from raining clouds performed with the Advanced Microwave Radiometer for Rain Identification (ADMIRARI) radiometer at 10.7-21-36.5 GHz during the Global Precipitation Measurement Ground Validation ""Cloud processes of the main precipitation systems in Brazil: A contribution to cloud resolving modeling and to the Global Precipitation Measurement"" (CHUVA) campaign held in Brazil in March 2010 represent a unique test bed for understanding three-dimensional (3D) effects in microwave radiative transfer processes. While the necessity of accounting for geometric effects is trivial given the slant observation geometry (ADMIRARI was pointing at a fixed 30 elevation angle), the polarization signal (i.e., the difference between the vertical and horizontal brightness temperatures) shows ubiquitousness of positive values both at 21.0 and 36.5 GHz in coincidence with high brightness temperatures. This signature is a genuine and unique microwave signature of radiation side leakage which cannot be explained in a 1D radiative transfer frame but necessitates the inclusion of three-dimensional scattering effects. We demonstrate these effects and interdependencies by analyzing two campaign case studies and by exploiting a sophisticated 3D radiative transfer suited for dichroic media like precipitating clouds.