This paper presents a new solar Extreme Ultraviolet (EUV) flux model for aeronomic calulations (EUVAC), which is based on the measured F74113 solar EUV reference spectrum. The model provides fluxes in the 37 wavelength bins that are in widespread use. This paper also presents cross sections to be used with the EUVAC flux model to calculate photoionization rates. The flux scaling for solar activity is accomplished using a proxy-based on the F10.7 index and its 81-day average together with the measured solar flux variation from the EUVS instrument on the Atmosphere Explorer E satellite. This new model produces 50-575 A integrated EUV fluxes in good agreement with rocket observations. The solar cycle variation of the chromospheric fluxes agrees well with the measured variation of the Lyman alpha flux between 1982 and 1988. In addition, the theoretical photoelectron fluxes, calculated using the new EUV flux model, are in good agreement with the solar minimum photoelectron fluxes from the Atmosphere Explorer E satellite and also with the solar maximum photoelectron fluxes from the Dynamics Explorer satellite. Its relative simplicity coupled with its ability to reproduce the 50-575 A solar EUV flux as well as the measured photoelectron spectrum makes the model well suited for aeronomic applications. However, EUVAC is not designed to accurately predict the solar flux variability for numerous individual lines.

EUVAC: A solar EUV Flux Model for aeronomic calculations

In this paper, our current understanding and recent advances in the study of ionospheric storms is reviewed, with emphasis on the F2-region. Ionospheric storms represent an extreme form of space weather with important effects on ground- and space-based technological systems. These phenomena are driven by highly variable solar and magnetospheric energy inputs to the Earth's upper atmosphere, which continue to provide a major difficulty for attempts now being made to simulate the detailed storm response of the coupled neutral and ionized upper atmospheric constituents using increasingly sophisticated global first principle physical models. Several major programs for coordinated theoretical and experimental study of these storms are now underway. These are beginning to bear fruit in the form of improved physical understanding and prediction of ionospheric storm effects at high, middle, and low latitude.

Ionospheric Storms — A Review

Four numerical simulations have been performed, at equinox, using a coupled thermosphere-ionosphere model, to illustrate the response of the upper atmosphere to geomagnetic storms. The storms are characterized by an increase in magnetospheric energy input at high latitude for a 12-hour period; each storm commences at a different universal time (UT). The initial response at high latitude is that Joule heating raises the temperature of the upper thermosphere and ion drag drives high-velocity neutral winds. The heat source drives a global wind surge, from both polar regions, which propagates to low latitudes and into the opposite hemisphere. The surge has the character of a large-scale gravity wave with a phase speed of about 600 m/s. Behind the surge a global circulation of magnitude 100 m/s is established at middle latitudes, indicating that the wave and the onset of global circulation are manifestations of the same phenomena. A dominant feature of the response is the penetration of the surge into the opposite hemisphere where it drives poleward winds for a few hours. The global wind surge has a preference for the night sector and for the longitude of the magnetic pole and therefore depends on the UT start time of the storm. A second phase of the meridional circulation develops after the wave interaction but is also restricted, in this case by the buildup of zonal winds via the Coriolis interaction. Conservation of angular momentum may limit the buildup of zonal wind in extreme cases. The divergent wind field drives upwelling and composition change on both height and pressure surfaces. The composition bulge responds to both the background and the storm-induced horizontal winds; it does not simply rotate with Earth. During the storm the disturbance wind modulates the location of the bulge; during the recovery the background winds induce a diurnal variation in its position. Equatorward winds in sunlight produce positive ionospheric changes during the main driving phase of the storm. Negative ionospheric phases are caused by increases of molecular nitrogen in regions of sunlight, the strength of which depends on longitude and the local time of the sector during the storm input. Regions of positive phase in the ionosphere persist in the recovery period due to decreases in mean molecular mass in regions of previous downwelling. Ion density changes, expressed as a ratio of disturbed to quiet values, exhibit a diurnal variation that is driven by the location of the composition bulge; this variation explains the ac component of the local time variation of the observed negative storm phase.

Response of the thermosphere and ionosphere to geomagnetic storms

The terrestrial ring current is an electric current flowing toroidally around the Earth, centered at the equatorial plane and at altitudes of ;10,000 - 60,000 km. Changes in this current are responsible for global decreases in the Earth's surface magnetic field, which are known as geomagnetic storms. Intense geomagnetic storms have severe effects on technological systems, such as disturbances or even permanent damage to tele- communication and navigation satellites, telecommuni- cation cables, and power grids. The main carriers of the storm ring current are positive ions, with energies from ;1 keV to a few hundred keV, which are trapped by the geomagnetic field and undergo an azimuthal drift. The ring current is formed by the injection of ions originating in the solar wind and the terrestrial ionosphere. The injection process involves electric fields, associated with enhanced magnetospheric convection and/or magneto- spheric substorms. The quiescent ring current is carried mainly by protons of predominantly solar wind origin, while geospace activity tends to increase the abundance (both absolute and relative) of O 1 ions, which are of ionospheric origin. During intense magnetic storms, the O 1 abundance increases dramatically, resulting in a rapid intensification of the ring current and an O 1 dominance around storm maximum. This compositional change affects, among other processes, the decay of the ring current through the species- and energy-dependent charge exchange and wave-particle scattering loss. En- ergetic neutral atoms, products of charge exchange, en- able global imaging of the ring current and are the most promising diagnostic tool of ring current evolution. This review will cover the origin of ring current particles, their transport and acceleration, the effects of composi- tional variations in the ring current, the effects of sub- storms on ring current growth, and the dynamics of ring current decay with an emphasis on the process of charge exchange and the potential for wave scattering loss.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/1999RG900009

The terrestrial ring current: Origin, formation, and decay

The Earth's ultraviolet airglow contains fundamental diagnostic information about the state of its upper atmosphere and ionosphere. Our understanding of the excitation and emission processes which are responsible for the airglow has undergone dramatic evolution from the earliest days of space research through the past several years during which a wealth of new information has been published from high-resolution spectroscopy and imaging experiments. This review of the field begins with an overview of the phenomenology: how the Earth looks in the ultraviolet. Next the basic processes leading to excitation of atomic and molecular energy states are discussed. These concepts are developed from first principles and applied to selected examples of day and night airglow; a detailed review of radiation transport theory is included. This is followed by a comprehensive examination of the current status of knowledge of individual emission features seen in the airglow, in which atomic physics issues as well as relevant atmospheric observations of major and minor neutral and ionic constituents are addressed. The use of airglow features as remote sensing observables is then examined for the purpose of selecting those species most useful as diagnostics of the state of the thermosphere and ionosphere. Imaging of the plasmasphere and magnetosphere is also briefly considered. A summary of upcoming UV remote sensing missions is provided.

Ultraviolet spectroscopy and remote sensing of the upper atmosphere

This paper presents calculations performed on the photoelectron impact to EUV photoionization rate ratios for atmospheric O, O2, and N2, and also for the production of N(+) by photodissociative ionization of N2. It was found that the ratios vary greatly with altitude. At high altitudes in the absence of photoelectron transport, the O(+) and N2(+) ionization rate ratios are about 0.35, but they increase with increasing optical depth; in the vicinity of the ionization peak, the photoelectron rate can exceed the EUV ionization rate for O(+) and N2(+). At high altitudes, the O2(+) ionization rate ratio is about half that of O(+); this ratio increases with increasing optical depth to reach a peak of about 0.4. There are also seasonal variations in the ratios at high altitudes depending on the magnitude of the conjugate photoelectron flux.

Ratios of photoelectron to EUV ionization rates for aeronomic studies

The degree to which the O(+) + N2 reaction rate is increased as a result of enhanced vibrational excitation of the N2 molecule in the thermosphere was investigated. It is found that the reaction rate may be sufficiently elevated in summer at solar maximum to decrease the peak O(+) density by a factor of 2, but there is only a small reduction in the winter peak density. Therefore the vibrational excitation of N2 acts to increase the magnitude of the seasonal anomaly. This work emphasizes the need for more laboratory work to clear up uncertainties in some of the key parameters.

A factor of 2 reduction in theoretical F2 peak electron density due to enhanced vibrational excitation of N2 in summer at solar maximum

An efficient two-stream auroral electron model is used to study the deposition of auroral energy and the dependence of auroral emission rates on characteristic energy. This model incorporates the concept of average energy loss to reduce the computation time. This simple two-stream model produces integrated emission rates that are in excellent agreement with the much more complex multistream model of Strickland et al. (1983) but disagrees with a recent study by Rees and Lummerzheim (1989) that indicates that the N2 second positive emission rate is a strongly decreasing function of the characteristic energy of the precipitating flux. These calculations reveal that a 10 keV electron will undergo approximately 160 ionizing collisions, with an average energy loss per collision of 62 eV before thermalizing. The secondary electrons are created with an average energy of 42 eV. When all processes including the backscattered escape fluxes are taken into account, the average energy loss per electron-ion pair is 35 eV in good agreement with laboratory results.

Auroral modeling of the 3371 Å emission rate: Dependence on characteristic electron energy

It is shown that fluxes of precipitating energetic O(+) that have been observed by satellites in the topside ionosphere can explain the magnitude of the N2+(1 N) (first negative) 3914-angstroms and N2(2 P) (second positive) 3371-angstroms emission rate observed during a mid-latitude aurora over Logan, Utah (41 degrees N, 111 degrees W), on 21-22 September 1982. Heavy particle precipitation has previously been invoked to explain the anomalously high populations in the upper vibrational levels of the N2+(1 N) system that are evident in the observed emissions. An improved model is used to investigate the impact of precipitating heavy ions on the atmosphere. The nocturnal thermospheric heating rate and ionization rate during one of these events can be equivalent to the daytime EUV heating and ionization rates. The observed spectrum can be explained by energetic O(+) precipitation to within the uncertainties of the inputs to the model.

The role of energetic O(+) precipitation in a mid-latitude aurora

The sensitivity of selected auroral emissions to anticipated changes in the neutral atmosphere was investigated from the results of a series of sensitivity studies conducted using an auroral emission code developed by Richards and Torr (1990). In particular, the behavior of OI 1356 and two Lyman Birge Hopfield (LBH) bands and their ratios to each other with changing atmospheric composition was examined. It was found that, for anticipated average uncertainties in the neutral atmosphere (factor 2 at auroral altitudes), the resultant change in the modeled intensities is comparable to or less than the uncertainty in the neutral atmosphere. The variation in the I 1356/I 1838 ratio over the equivalent of a solar cycle is less than 50 percent, and the summer-to-winter changes are approximately a factor of 2.

D. G. Torr

Papers

Ratios of photoelectron to EUV ionization rates for aeronomic studies

A factor of 2 reduction in theoretical F2 peak electron density due to enhanced vibrational excitation of N2 in summer at solar maximum

Auroral modeling of the 3371 Å emission rate: Dependence on characteristic electron energy

The role of energetic O(+) precipitation in a mid-latitude aurora

The dependence of modeled OI 1356 and N2 Lyman Birge Hopfield auroral emissions on the neutral atmosphere