The consequences of high latitude particle precipitation on global thermospheric dynamics

A previous comparison of experimental measurements of thermospheric winds with simulations using a global self-consistent three-dimensional time-dependent model confirmed a necessity for a high latitude source of energy and momentum acting in addition to solar u.v. and e.u.v. heating. During quiet geomagnetic conditions, the convective electric field over the polar cap and auroral oval seemed able to provide adequate momentum input to explain the thermospheric wind distribution observed in these locations. However, it seems unable to provide adequate heating, by the Joule mechanism, to complete the energy budget of the thermosphere and, more importantly, unable to provide the high latitude input required to explain mean meridional winds at mid-latitudes. In this paper we examine the effects of low energy particle precipitation on thermospheric dynamics and energy budget. Modest fluxes over the polar cap and auroral oval, of the order of 0.4 erg cm ⁻²/s, are consistent with satellite observations of the particles themselves and with photometer observations of the OI and OII airglow emissions. Such particle fluxes, originating in the dayside magnetosheath cusp region and in the nightside central plasma sheet, heat the thermosphere and modify mean meridional winds at mid-latitudes without enhancing the OI 557.7 line, or the ionization of the lower thermosphere (and thus enhancing the auroral electrojets), neither of which would be consistent with observations during quiet geomagnetic conditions.     

The generation of vertical thermospheric winds and gravity waves at auroral latitudes—II. Theory and numerical modelling of vertical winds

Recent observations of strong vertical thermospheric winds and the associated horizontal wind structures, using the 01(3P-1D)nm emission line, by ground-based Fabry-Perot interferometers in Northern Scandinavia have been described in an accompanying paper (Paper I). The high latitude thermosphere at a height of 200–300 km displays strong vertical winds (30–50m ms^?1)of a persistent nature in the vicinity of the auroral oval even during relatively quiet geomagnetic conditions. During an auroral substorm, the vertical (upward) wind in the active region, including that invaded by a Westward Travelling Surge, may briefly(10–30 min)exceed 150 m s^?1. Very large and rapid changes of horizontal wind structure (up to 500 m^?1 in 30 min) usually accompany such large impulsive vertical winds. Magnetospheric energy and momentum sources generate large vertical winds of both a quasi-steady nature and of a strongly time-dependent nature. The thermospheric effects of these sources can be evaluated using the UCL three-dimensional, time-dependent thermospheric model. The auroral oval is, under average geomagnetic conditions, a stationary source of significant vertical winds (10–40 m s^?1). In large convective events (directly driven by a strong momentum coupling from the solar wind) the magnitude may increase considerably. Auroral substorms and Westward Travelling Surges appear to be associated with total energy disposition rates of several tens to more than 100 erg cm^?2s^?1, over regions of a few hours local time, and typically 2–5° of geomagnetic latitude (approximately centred on magnetic midnight). Such deposition rates are needed to drive observed time-dependent vertical (upward) winds of the order of 100–200m s^?1.The response of the vertical winds to significant energy inputs is very rapid, and initially the vertical lifting of the atmosphere absorbs a large fraction (30% or more) of the total substorm input. Regions of strong upward winds tend to be accompanied in space (and time) by regions of rather lower downward winds, and the equatorward propagation of thermospheric waves launched by auroral substorms is extremely complex.     

Electrical coupling of the E- and F-regions and its effect on F-region drifts and winds

Electric currents, generated by thermospheric winds, flow along the geomagnetic field lines linking the E-and F-regions. Their effects on the electric field distribution are investigated by solving the electrical and dynamical equations. The input data include appropriate models of the F-region tidal winds, the thermospheric pressure distribution and the E-and F-layer concentrations. At the magnetic equator, the calculated neutral air wind at 240 km height has a prevailling eastward component of 55 m sec^-1 and the west-east and vertical ion drifts agree in their general form with incoherent scatter data from Jicamarca     

The westward thermospheric jet-stream of the evening auroral oval

One of the most consistent and often dramatic interactions between the high latitude ionosphere and the thermosphere occurs in the vicinity of the auroral oval in the afternoon and evening period. Ionospheric ions, convected sunward by the influence of the magnetospheric electric field, create a sunward jet-stream in the thermosphere, where wind speeds of up to 1 km s^?1 can occur. This jet-stream is nearly always present in the middle and upper thermosphere (above 200 km altitude), even during periods of very low geomagnetic activity. However, the magnitude of the winds in the jet-stream, as well as its location and range in latitude, each depend on geomagnetic activity. On two occasions, jet-streams of extreme magnitude have been studied using simultaneous ground-based and satellite observations, probing both the latitudinal structure and the local time dependence. The observations have then been evaluated with the aid of simulations using a global, three-dimensional, time-dependent model of thermospheric dynamics including the effects of magnetospheric convection and particle precipitation. The extreme events, where sunward winds of above 800 ms^?1 are generated at relatively low geomagnetic latitudes (60–70°) require a greatly expanded auroral oval and large cross-polar cap electric field ( ~ 150 kV). These in turn are generated by a persistent strong Interplanetary Magnetic Field, with a large southward component. Global indices such as K_p are a relatively poor indicator of the magnitude and extent of the jet-stream winds.     

The calculated and observed diurnal variation of the ionosphere over Millstone Hill on 23–24 March 1970

A numerical model of current F-region theory is use to calculate the diurnal variation of the mid-latitude ionospheric F-region over Millstone Hill on 23–24 March 1970, during quiet geomagnetic conditions. From the solar EUV flux, the model calculates at each altitude and time step primary photoelectron spectra and ionization rates of various ion species. The photoelectron transport equation is solved for the secondary ionization rates, photoelectron spectra, and various airglow excitation rates. Five ion continuity equations that include the effects of transport by diffusion, magnetospheric-ionospheric plasma transport, electric fields, and neutral winds are solved for the ion composition and electron density. The electron and ion temperatures are also calculated using the heating rates determined from chemical reactions, photoelectron collisions, and magnetospheric-ionospheric energy transport. The calculations are performed for a diurnal cycle considering a stationary field tube co-rotating with the Earth; only the vertical plasma drift caused by electric fields perpendicular to the geomagnetic field line is allowed but not the horizontal drift. The boundary conditions used in the model are determined from the incoherent scatter radar measurements of T_e, T_i and O⁺ flux at 800km over Millstone Hill (Evans, 1971a). The component of the neutral thermospheric winds along the geomagnetic field has an important influence on the overall ionospheric structure. It is determined from a separate dynamic model of the neutral thermosphere, using incoherent scatter radar measurements.The calculated diurnal variation of the ionospheric structure agrees well with the values measured by the incoherent scatter radar when certain restrictions are placed on the solar EUV flux and model neutral atmospheric compositions. Namely, the solar EUV fluxes of Hinteregger (1970) are doubled and an atomic oxygen concentration of at least $\text{10}{}^{11}\text{cm}{}^3$ at 120 km is required for the neutral model atmosphere. Calculations also show that the topside thermal structure of the ionosphere is primarily maintained by a flow of heat from the magnetosphere and the night-time F2-region is maintained in part by neutral winds, diffusion, electric fields, and plasma flow from the magnetosphere. The problem of maintaining the calculated night-time ionosphere at the observed values is also discussed.     

The generation of vertical thermospheric winds and gravity waves at auroral latitudes—I. Observations of vertical winds

Observations of vertical and horizontal thermospheric winds, using the OI (3P-1D) 630 nm emission line, by ground-based Fabry-Perot interferometers in Northern Scandinavia and in Svalbard (Spitzbergen) have identified sources of strong vertical winds in the high latitude thermosphere. Observations from Svalbard (78.2N 15.6E) indicate a systematic diurnal pattern of strong downward winds in the period 06.00 U.T. to about 18.00 U.T., with strong upward winds between 20.00 U.T. and 05.00 U.T. Typical velocities of 30 m s^?1 downward and 50 m s^?1 upward occur, and there is day to day variability in the magnitude (30–80 m s^?1) and phase (+/- 3 h) in the basically diurnal variation. Strong and persistent downward winds may also occur for periods of several hours in the afternoon and evening parts of the auroral oval, associated with the eastward auroral electrojet (northward electric fields and westward ion drifts and winds), during periods of strong geomagnetic disturbances. Average downward values of 30–50 m s^?1 have been observed for periods of 4–6 h at times of large and long-lasting positive bay disturbances in this region. It would appear that the strong vertical winds of the polar cap and disturbed dusk auroral oval are not in the main associated with propagating wave-like features of the wind field. A further identified source is strongly time-dependent and generates very rapid upward vertical motions for periods of 15–30 min as a result of intense local heating in the magnetic midnight region of the auroral oval during the expansion phase of geomagnetic disturbances, and accompanying intense magnetic and auroral disturbances. In the last events, the height-integrated vertical wind (associated with a mean altitude of about 240 km) may exceed 100–150 m s^?1. These disturbances also invariably cause major time-dependent changes of the horizontal wind field with, for example, horizontal wind changes exceeding 500 m s^?1 within 30 min. The changes of vertical winds and the horizontal wind field are highly correlated, and respond directly to the local geomagnetic energy input. In contrast to the behaviour observed in the polar cap or in the disturbed afternoon auroral oval, the ‘expansion phase’ source, which corresponds to the classical ‘auroral substorm’, generates strong time-dependent wind features which may propagate globally. This source thus directly generates one class of thermospheric gravity waves. In this first paper we will consider the experimental evidence for vertical winds. In a second paper we will use a three-dimensional time-dependent model to identify the respective roles of geomagnetic energy and momentum in the creation of both classes of vertical wind sources, and consider their propagation and effects on global thermospheric dynamics.     

Horizontal and vertical winds and temperatures in the equatorial thermosphere: Measurements from Natal,Brazil during August–September 1982

Fabry-Perot interferometer measurements of Doppler shifts and widths of the 630.0 nm nightglow line have been used to determine the neutral winds and temperatures in the equatorial thermosphere over Natal, Brazil during August–September 1982. During this period, in the early night (2130 U.T.) the average value of the horizontal wind vector was 95 m s^?1 at 100° azimuth, and the temperature varied from a low of 950 K during geomagnetically quiet conditions to a high of ~ 1400 K during a storm (6 September). The meridional winds were small, ?, 50 m s^?1, and the eastward zonal winds reached a maximum value 1–3 h after sunset, in qualitative agreement with TGCM predictions. On 26 August, an observed persistent convergence in the horizontal meridional flow was accompanied by a downward vertical velocity and an increase in the thermospheric temperature measured overhead. Oscillations with periods of 40–45 min in both the zonal and vertical wind velocities were observed during the geomagnetic storm of 6 September, suggesting gravity wave modulation of the equatorial thermospheric flow.     

Solar transients disturbing the terrestrial magnetic environment at higher latitudes

Geomagnetic field variations during five major Solar Energetic Particle (SEP) events of solar cycle 23 have been investigated in the present study. The SEP events of 1 October 2001, 4 November 2001, 22 November 2001, 21 April 2002 and 14 May 2005 have been selected to study the geomagnetic field variations at two high-latitude stations, Thule (77.5^° N, 69.2^° W) and Resolute Bay (74.4^° E, 94.5^° W) of the northern polar cap. We have used the GOES proton flux in seven different energy channels (0.8–4 MeV, 4–9 MeV, 9–15 MeV, 15–40 MeV, 40–80 MeV, 80–165 MeV, 165–500 MeV). All the proton events were associated with geoeffective or Earth directed CMEs that caused intense geomagnetic storms in response to geospace. We have taken high-latitude indices, AE and PC, under consideration and found fairly good correlation of these with the ground magnetic field records during the five proton events. The departures of the H component during the events were calculated from the quietest day of the month for each event and have been represented as ΔH _THL and ΔH _RES for Thule and Resolute Bay, respectively. The correspondence of spectral index, inferred from event integrated spectra, with ground magnetic signatures ΔH _THL and ΔH _RES along with Dst and PC indices have been brought out. From the correlation analysis we found a very strong correlation to exist between the geomagnetic field variation (ΔHs) and high-latitude indices AE and PC. To find the association of geomagnetic storm intensity with proton flux characteristics we derived the correspondence between the spectral indices and geomagnetic field variations (ΔHs) along with the Dst and AE index. We found a strong correlation (0.88) to exist between the spectral indices and ΔHs and also between spectral indices and AE and PC.     

The composition,structure, temperature and dynamics of the upper thermosphere in the polar regions during October to December 1981

During the period October to December 1981, the Dynamics Explorer-2 (DE-2) spacecraft successively observed the South polar and the North polar regions, and recorded the temperature, composition and dynamical structure of the upper thermosphere. In October 1981, perigee was about 310 km altitude, in the vicinity of the South Pole, with the satellite orbit in the 09.00–21.00 L.T. plane. During late November and December, the perigee had precessed to the region of the North Pole, with the spacecraft sampling the upper thermosphere in the 06.00 18.00 L.T. plane. DE-2 observed the meridional wind with a Fabry-Perot interferometer (FPI), the zonal wind with the wind and temperature spectrometer (WATS), the neutral temperature with the FPI, and the neutral atmosphere composition and density with the neutral atmosphere composition spectrometer (NACS). A comparison between the South (summer) Pole and the North (winter) Pole data shows considerable seasonal differences in all neutral atmosphere parameters. The region of the summer pole, under similar geomagnetic and solar activity conditions, and at a level of about 300 km, is about 300 K warmer than that of the winter pole, and the density of atomic oxygen is strongly depleted (and nitrogen enhanced) around the summer pole (compared with the winter pole). Only part of the differences in temperature and composition structure can be related to the seasonal variation of solar insolation, however, and both polar regions display structural variations (with latitude and Universal Time) which are unmistakeable characteristics of strong magnetospheric forcing. The magnitude of the neutral atmosphere perturbations in winds, temperature, density and composition within both summer and winter polar regions all increase with increasing levels of geomagnetic activity.The UCL 3-dimensional time dependent global model has been used to simulate the diurnal, seasonal and geomagnetic response of the neutral thermosphere, attempting to follow the major features of the solar and geomagnetic inputs to the thermosphere which were present during the late 1981 period.In the UCL model, geomagnetic forcing is characterized by semi-empirical models of the polar electric field which show a dependence on the Y component of the Interplanetary Magnetic Field, due to Heppner and Maynard (1983), It is possible to obtain an overall agreement, in both summer and winter hemispheres, with the thermospheric wind structure at high latitudes, and to explain the geomagnetic control of the combined thermal and compositional structure both qualitatively and quantitatively. To obtain such agreement, however, it is essential to enhance the polar ionosphere as a consequence of magnetospheric particle precipitation, reflecting both widespread auroral (kilovolt) electrons, and “soft” cusp and polar cap sources. Geomagnetic forcing of the high latitude thermosphere cannot be explained purely by a polar convective electric field, and the thermal as well as ionising properties of these polar and auroral electron sources are crucial components of the total geomagnetic input.     

A model of thermospheric temperature and composition

An empirical model of thermospheric temperature (T_∞T₁₂₀, and s) and composition (H, He, N, O, N₂, O₂, and Ar) was derived from measurements of 8 satellites (AE-C, AE-E, AEROS-A, AEROS-B, ARIEL-3, ESRO-4, OGO-6, and SAN MARCO-3) and 4 incoherent scatter stations (Arecibo, Jicamarca, Millstone Hill, and St Santin). The altitude covered extends from 120 km up to about 600 km over the time period 1967 to 1976. The analytical framework used in the model resembles closely the MSIS setup: time independent terms, solar flux terms, geomagnetic activity (K_p) effect, annual (semiannual) and diurnal (semidiurnal, terdiurnal) variations, longitudinal terms, the U.T. effect, and corrections compensating for deviations from diffusive equilibrium at altitudes below 200 km. The model describes quiet to medium disturbed geomagnetic conditions (K_p ? 4) at solar fluxes (10.7cm) ranging from 60 to 180 × 10^?22 Wm^?2Hz^?1. To get an impression of the accuracy presently obtained, the model is compared with MSIS, Jacchia (1977), and the models of Thuillier (T_∞ and Engebretson (N). The best agreement is found for the temperature and the constituents He, O, and N₂ with increasing deviations in the order of H, N, Ar, and O₂.     

The Doppler Imaging System: Initial observations of the auroral thermosphere

During December 1982, a novel Fabry-Perot interferometer—a Doppler Imaging System (DIS)— was used at Kiruna Geophysical Institute (KGI), Sweden (67.8°N, 21.2°E) to complement a series of coordinated observations of global thermospheric dynamics utilizing a number of conventional ground-based Fabry-Perot interferometers and the NASA Dynamics Explorer satellite. The DIS is an interferometer with two unique attributes : it has a luminosity or étendue more than one hundred times that of the conventional Fabry-Perot interferometer, and it is also capable of deducing a two-dimensional velocity field of a suitable line-emitting areal source by independently measuring the Doppler shift at a large number of points within the field of view. On 17 December 1982, a very large geomagnetic Storm Sudden Commencement (08.05 U.T.) preceded a major geomagnetic disturbance. During this disturbance, Northern Scandinavia was influenced by a strong eastward auroral electrojet for an extended period (10–19 U.T.). The DIS was able to observe the dynamical response of the upper thermosphere to this event in conjunction with a second Fabry-Perot interferometer (FPI) at KGI. Westward thermospheric winds of about 900 m s^?1 occurred during the disturbance and, at the peak of the disturbance, the combined DIS and FPI observations indicate that the thermospheric flow was quite chaotic. Fluctuations of the order of ± 150 ms^?1, associated with spatial scales of the order of 100 or 200 km occurred within the mean westward flow inside the 800 km diameter region observed from Kiruna.     

The behaviour of the midlatitude F-region at the time of four magnetospheric substorms during the night of 29/30 October 1973

The data from observations of the geomagnetic field, ionospheric parameters and atmospheric emissions, carried out at four midlatitude station in Bulgaria are analysed. The observations refer to the geomagnetic disturbance on 28/30 October 1973 (K_p^max = 7) and also to a very quiet period before it. It is shown that all four geomagnetic substorms during the night of 29/30 October influenced the midlatitude F-region. This is indicated by a lowering of the height of the F-region by ca. 50–70 km. Owing to this downward drift of ionisation the dissociative recombination and the intensity of the red line is accordingly increased. As an explanation of this phenomenon we suggest the action of the electric fields, which can at the same time be transported from the magnetosphere to the ionosphere.     

The solar wind-magnetosphere energy coupling and magnetospheric disturbances

The recent finding of the solar wind-magnetosphere energy coupling function ε has advanced significantly our understanding of magnetospheric disturbances. It is shown that the magnetosphere-ionosphere coupling system responds somewhat differently to three different input energy flux levels of ε. As ε increases from < 10¹⁷ erg s^?1 to > 0¹⁹ erg s^?1, typical responses of the magnetosphere-ionosphere coupling system are:ε < 10¹⁷ erg s^?1: an enhancement of the S_q^p, etc,ε ≈ 10¹⁸ erg s^?1: substorm onset,10¹⁸ erg s^?1 < ε < 10¹⁹ erg s^?1f: a typical substorm,ε >10¹⁹ erg s^?1: an abnormal growth of the ring current belt, resulting in a magnetospheric storm.It is stressed that the magnetospheric substorm results as a direct response of the magnetosphere to a rise and fall of ε above ≈ 10¹⁸ erg s^?1, so that it is not caused by a sudden conversion of magnetic energy accumulated prior to substorm onset. The variety of the development of the main phase of geomagnetic storms is also primarily controlled by ε.     

Theoretical studies of storm effects in ionospheric total electron content

Investigations have been made of the effects produced by thermospheric winds, composition changes and magnetospheric electric fields on the ionospheric F-layer during disturbed conditions. The results of the computations suggest that a combination of realistic temperature and electric field changes would explain fairly satisfactorily the observed changes in total electron content.     

An estimation of the electric potential field generated by the diurnal atmospheric tide with the first negative mode excited in the lower ionosphere

The electric potential field in the ionosphere, which is generated by the atmospheric tide with the diurnal first negative mode excited in the lower ionosphere, is estimated. An equivalent mode is used instead of the exact mode in the following approximations : (i) the horizontal structure of the energy input is composed of the first two terms of the series expansion of the Hough function with the associated Legendre functions, and (ii) the Coriolis force is assumed to be constant in the equation of motion and the value appropriate to colatitude 45 being used. It is found that the equivalent mode gives a good approximation of the exact mode at midlatitudes in the neutral atmosphere. The result shows that the oscillation with its temperature variation of about 40 K in the excitation region generates the horizontal electric potential field of 10^?3 V/m in magnitude. The horizontal pattern of the field distribution is similar to that estimated by many authors. The maximum vertical electric current is about 3 × 10^?9 A/m² in the dynamo layer. It is found that the “virtual current,” introduced by Volland (1970) as a vertical current in and above the dynamo layer, in order to obtain a better agreement between theory and observations of the location of the diurnal symmetric Sq current vortex, is actually to be expected. The vertical current is very small, but this small steady current is necessary for the electric potential field of the dynamo layer to be mapped into the magnetosphere along the lines of the geomagnetic field.     

Effects of convection electric field on the distribution of ring current type protons

The topology of the boundaries of penetration (or inversely the boundaries of the forbidden regions) of 90° pitch angle equatorial protons with energies less than 100 keV are explored for an equatorial convection E-field which is directed in general from dawn to dusk. Due to the dependence of drift path on energy (or magnetic moment) complex structural features are expected in the proton energy spectra detected on satellites since the penetration distance of a proton is not a monotonically increasing or decreasing function of energy. During a storm when the convection E is enhanced, model calculations predict elongations of the forbidden regions analogous to plasmatail extensions of the plasmasphere. Following a reduction in the convection field, spiral-structured forbidden regions can occur. Structural features inherent to large scale convection field changes may be seen in the noselike proton spectrograms observed near dusk by instrumentation on the satellite Explorer 45 (S³) (Smith and Hoffman, 1974). These nose events are modelled by using an electric field model developed originally by Volland (1973). The strength of the field is related to K_p through night-time equatorial plasmapause measurements.     

On the dispersal of artificially-injected gases in the night-time atmosphere

We have studied the extent to which various transport processes affect the dispersal of a gas artificially injected into the night-time atmosphere at F-region altitudes. In addition to diffusion, we have found that nonlinear acceleration, viscous stress, and thermospheric winds affect the dispersal of the injected gas. The magnitude of the effect depends on the atmospheric density, which is a function of solar activity. For an injected H₂ gas, non-linear acceleration and viscous stress rapidly become more important than diffusion above about 300 km for low solar activity (T_∞ = 750K), 340 km for moderate solar activity (T_∞ = 1000K), and 400 km for high solar activity (T_∞ = 1500K). For an injected H₂O gas, the corresponding altitudes are 350, 400, and 470 km for low, moderate and high solar activity, respectively. The effect of nonlinear acceleration and viscous stress is to retard the expansion of the injected gas. Thermospheric winds of 150–400 m s^?1 are important at altitudes near and below the F-region peak electron density. These winds act to transport the injected gas in the wind direction and this affects the shape and temporal development of the subsequent ionospheric hole. Because the H₂O diffusion coefficient is smaller than the H₂ diffusion coefficient, winds are more important for H₂O than for H₂.     

Jupiter's radiation belts

A model for the production and loss of energetic electrons in Jupiter's radiation belt is presented. It is postulated that the electrons originate in the solar wind and are diffused in toward the planet by perturbations which violate the particles' third adiabatic invariant. At large distances, magnetic perturbations, electric fields associated with magnotospheric convection, or interchange instabilities driven by thermal plasma gradients may drive the diffusion. Inside about 10 R_J the diffusion is probably driven by electric fields associated with the upper atmosphere dynamo which is driven by neutral winds in the ionosphere. The diurnal component of the dynamo wind fields produces a dawn-dusk asymmetry in the decimetric radiation from the electrons in the belts, and the lack of obvious measured asymmetries in the decimetric radiation measurements provides estimates of upper limits for these Jovian ionospheric neutral winds. The average diurnal winds are less than or comparable to those on earth, but only modest fluctuating winds are required to drive the energetic electron diffusion referred to above.The winds required to diffuse the energetic particles across the orbit of the satellite lo in a time equal to their drift period are also estimated. If Io is non-conducting, modest winds are required, but if Io is conducting, only small winds are needed. It is concluded that both protons and electrons are diffused in from the solar wind to small distances without serious losses occurring due to the particles being swept up by the satellites.Consideration of proton and electron diffusion in energy shows that once the electrons become relativistic, the ratio of proton to electron energy increases. Thus, if protons and electrons have the same energy in the solar wind, when the electrons reach nMeV, the protons will be nMeV if n ? 1 or n² MeV if n ? 1. If the proton-to-electron energy ratio is initially, e.g., 5, then these figures are 5n and 5n², respectively.