On I(5577 Å) and I (7620 Å) auroral emissions and atomic oxygen densities

A model of auroral electron deposition processes has been developed using Monte Carlo techniques to simulate electron transport and energy loss. The computed differential electron flux and pitch angle were compared with in situ auroral observations to provide a check on the accuracy of the model. As part of the energy loss process, a tally was kept of electronic excitation and ionization of the important atomic and molecular states. The optical emission rates from these excited states were computed and compared with auroral observations of (3914 Å), (5577 Å), (7620 Å) and (N₂VK). In particular, the roles played by energy transfer from N₂(A³⁺_u) and by other processes in the excitation of O(¹S) and O₂(b¹⁺_g) were investigated in detail. It is concluded that the N₂(A³⁺_u) mechanism is dominant for the production of OI(5577 Å) in the peak emission region of normal aurora, although the production efficiency is much smaller than the measured laboratory value; above 150 km electron impact on atomic oxygen is dominant. Atomic oxygen densities in the range of 0.75±0.25 MSIS-86 [O] were derived from the optical comparisons for auroral latitudes in mid-winter for various levels of solar and magnetic activity.     

Subauroral red arcs as a conjugate phenomenon: comparison of OV1-10 satellite data with numerical calculations

This study compares the OV1-10 satellite measurements of the integral airglow intensities at 630 nm in the SAR arc regions observed in the northern and southern hemisphere as a conjugate phenomenon, with the model results obtained using the time-dependent one-dimensional mathematical model of the Earth ionosphere and plasmasphere (the IZMIRAN model) during the geomagnetic storm of the period 15–17 February 1967. The major enhancements to the IZMIRAN model developed in this study are the inclusion of He+ ions (three major ions: O⁺ H⁺ and He⁺ and three ion temperatures), the updated photochemistry and energy balance equations for ions and electrons, the diffusion of NO⁺ and O⁺₂ ions and O(¹D) and the revised electron cooling rates arising from their collisions with unexcited N₂, O₂ molecules and N₂ molecules at the first vibrational level. The updated model includes the option to use the models of the Boltzmann or non-Boltzmann distributions of vibrationally excited molecular nitrogen. Deviations from the Boltzmann distribution for the first five vibrational levels of N₂ were calculated. The calculated distribution is highly non-Boltzmann at vibrational levels v > 2 and leads to a decrease in the calculated electron density and integral intensity at 630 nm in the northern and southern hemispheres in comparison with the electron density and integral intensity calculated using the Boltzmann vibrational distribution of N₂. It is found that the intensity at 630 nm is very sensitive to the oxygen number densities. Good agreement between the modeled and measured intensities is obtained provided that at all altitudes of the southern hemisphere a reduction of about factor 1.35 in MSIS-86 atomic oxygen densities is included in the IZMIRAN model with the non-Boltzm-ann vibrational distribution of N². The effect of using of the O(¹D) diffusion results in the decrease of 4–6% in the calculated integral intensity of the northern hemisphere and 7–13% in the calculated integral intensity of the southern hemisphere. It is found that the modeled intensities of the southern hemisphere are more sensitive to the assumed values of the rate coefficients of O⁺(⁴S) ions with vibrationally excited nitrogen molecules and quenching of O⁺(²D) by atomic oxygen than the modeled intensities of the northern hemisphere.     

Interference filter spectral imaging of twilight O+(2P-2D) emission

A spectral imager specifically designed to measure the O⁺(²P-²D) emission in the thermosphere during twilight has been constructed and tested in Toronto (43.8°N, 79.3°W), and found to show promise for long-term and campaign-mode operations. A modification of the mesopause oxygen rotational temperature imager (MORTI), it consists basically of a narrow-band interference filter (0.14 nm bandwidth) to separate wavelengths as a function of off-axis angle, a lens to focus the spectrum into a series of concentric rings, and a focal plane array (CCD) to record the spectral images in digital form. The instrument was built with two fields of view, one for the zenith and one for 20° above the horizon, movable to track the azimuth of the Sun, in order to provide appropriate data for inversion. Data gathered during June 1991 provided measurements of the column-integrated emission rate with a precision of about 3%. An atomic oxygen profile was deduced that showed good agreement with that predicted by the MSIS-90 model atmosphere. Geomagnetically induced variations of the O⁺ lines, calcium spectra resulting from meteor showers, and OH nightglow were also observed.     

Rocket observations of atomic oxygen density and airglow emission rate in the WAVE2000 campaign

Atomic oxygen density and airglow volume emission rate profiles measured in the rocket experiment S-310-29 carried out as a part of the Waves in Airglow Structures Experiment over Kagoshima in 2000 (WAVE2000) campaign are presented, and the excitation processes of the atomic oxygen 557.7 nm line and the molecular oxygen atmospheric band airglow emissions are discussed. The volume emission rate profiles calculated from the measured atomic oxygen densities using the methods and parameters proposed by the ETON campaign (EATON model) are found to well represent the shapes of measured twin-peak emission rate profiles seen during this campaign suggesting that the EATON model is valid for a perturbed atmosphere. There is some discrepancy in the modeled absolute values for the emissions. Applying the model to the current O density measurements results in predicted emission rates that are a factor of 3.2 and 1.5 too high for the 557.7 nm line and Atmospheric band, respectively. This suggests either that the atomic oxygen densities of the present campaign are too large by a factor of 1.2 (=1.5^1/2) to 1.5 (=3.2^1/3), or that those of the ETON campaign were too small by the same factor or that the combined errors in both campaigns can account for the discrepancy (the modeled volume emission rates of the 557.7 nm line and atmospheric band are roughly in proportion to [O]³ and [O]², respectively). Our present data findings do not favor the 2-step more than the 1-step excitation process for the atmospheric band because a calculation of the quantum efficiency based on the observed O density does not show a steep gradient around 100 km.     

An anomalous subauroral red arc on 4 August, 1972: comparison of ISIS-2 satellite data with numerical calculations

This study compares the Isis II satellite measurements of the electron density and temperature, the integral airglow intensity and volume emission rate at 630 nm in the SAR arc region, observed at dusk on 4 August, 1972, in the Southern Hemisphere, during the main phase of the geomagnetic storm. The model results were obtained using the time dependent one-dimensional mathematical model of the Earth’s ionosphere and plasmasphere (the IZMIRAN model). The major enhancement to the IZMIRAN model developed in this study to explain the two component 630 nm emission observed is the analytical yield spectrum approach to calculate the fluxes of precipitating electrons and the additional production rates of N⁺₂, O⁺₂, O⁺(⁴S), O⁺(²D), O⁻(²P), and O⁺(²P) ions, and O(¹D) in the SAR arc regions in the Northern and Southern Hemispheres. In order to bring the measured and modelled electron temperatures into agreement, the additional heating electron rate of 1.05 eV cm⁻³ s⁻¹ was added in the energy balance equation of electrons at altitudes above 5000 km during the main phase of the geomagnetic storm. This additional heating electron rate determines the thermally excited 630 nm emission observed. The IZMIRAN model calculates a 630 nm integral intensity above 350 km of 4.1 kR and a total 630 nm integral intensity of 8.1 kR, values which are slightly lower compared to the observed 4.7 kR and 10.6 kR. We conclude that the 630 nm emission observed can be explained considering both the soft energy electron excited component and the thermally excited component. It is found that the inclusion of N₂(v > 0) and O₂(v > 0) in the calculations of the O⁺(⁴S) loss rate improves the agreement between the calculated N_e and the data on 4 August, 1972. The N₂(v > 0) and O₂(v > 0) effects are enough to explain the electron density depression in the SAR arc F-region and above F2 peak altitude. Our calculations show that the increase in the O⁺ + N₂ rate factor due to the vibrationally excited nitrogen produces the 5–19% reductions in the calculated quiet daytime peak density and the 16–24% decrease in NmF2 in the SAR arc region. The increase in the O⁺ + N₂ loss rate due to vibrationally excited O₂ produces the 7–26% decrease in the calculated quiet daytime peak density and the 12–26% decrease in NmF2 in the SAR arc region. We evaluated the role of the electron cooling rates by low-lying electronic excitation of O₂(a¹δ_g) and O₂(b¹σ_g⁺), and rotational excitation of O₂, and found that the effect of these cooling rates on T_e can be considered negligible during the quiet and geomagnetic storm period 3–4 August, 1972. The energy exchange between electron and ion gases, the cooling rate in collisions of O(³P) with thermal electrons with excitation of O(¹D), and the electron cooling rates by vibrational excitation of O₂ and N₂ are the largest cooling rates above 200 km in the SAR arc region on 4 August, 1972. The enhanced IZMIRAN model calculates also number densities of N₂(B³π_g⁺), N₂(C³π_u), and N₂(A³σ_u⁺) at several vibrational levels, O(¹S), and the volume emission rate and integral intensity at 557.7 nm in the region between 120 and 1000 km. We found from the model that the integral integral intensity at 557.7 nm is much less than the integral intensity at 630 nm.     

Planetary distribution of the intensity of Auroral luminosity obtained using a model of Aurora precipitation

A model of auroral precipitation (AP) developed on the basis of statistical processing of DMSP F6 and F7 satellite data (Vorobjev and Yagodkina, 2005, 2007) was used for the calculation of the global distribution of the auroral luminosity in different spectral ranges. The algorithm for the calculation of the integral intensity in bands N₂ LBH (170.0 nm), ING N ₂ ⁺ (391.4 nm), 1PG N₂ (669.0 nm), and (OI) 557.7-nm emission is shown in detail. The processes of formation of electronically excited atoms O(¹S) as a result of the transport of excitation energy from metastable state N₂(A³Σ _u ⁺ ), excitation of O(³P) by primary and secondary electrons, and dissociative recombination were taken into account to calculate the intensity of emission at 557.7 nm. A high correlation between the model distribution of the auroral luminosity in the UV spectral range and the observations of the Polar satellite is demonstrated.     

Determination of the lifetime of molecules O 2 + and estimation of the flux of penetrating electrons generating a cloud of excess concentration of O 2 + in the ionosphere at altitudes above 220 km

The occurrence of anomalous (nonthermal) profiles of green emission of oxygen atoms detected with a Fabry-Perot spectrometer in auroras with the effect of a rapid decrease in the intensity of the wings of their dissociative component has been investigated. Based on an analysis of these measured profiles, it has been found that the characteristic time of recombination of a molecular oxygen ion at altitudes of 200–400 km is about 5–7 s. It appears that these molecular ions occur in a horizontally limited region of the auroral ionosphere as a result of ionization by a space localized flux of soft electrons with energies of 0.2–0.4 keV penetrating up to altitudes of 200 km. The estimation of the electron flux produces a value of 10¹⁰–10¹³ electrons cm^?2 s^?1. They generate the excess concentration n(O ₂ ⁺ ) ~ 5.6 × 10⁵ cm^?3.     

The role of vibrationally excited nitrogen in the formation of the mid-latitude negative ionospheric storms

Millstone Hill ionospheric storm time measurements of the electron density and temperature during the ionospheric storms (15-16 June 1965; 29–30 September 1969 and 17–18 August 1970) are compared with model results. The model of the Earth’s ionosphere and plasmasphere includes interhemispheric coupling, the H⁺, O⁺(⁴S), O⁺(²D), O⁺(²P), NO⁺, O⁺₂ and N⁺₂ ions, electrons, photoelectrons, the electron and ion temperature, vibrationally excited N₂ and the components of thermospheric wind.In order to model the electron temperature at the time of the 16 June 1965 negative storm, the heating rate of the electron gas by photoelectrons in the energy balance equation was multiplied by the factors 5–30 at he altitude above 700 km for the period 4.50-12.00 LT, 16 June 1965. The [O]/[N₂] MSIS-86 decrease and vibrationally excited N₂ effects are enough to account for the electron density depressions at Millstone Hill during the three storms. The factor of 2 (for 27–30 September 1969 magnetic storm) and the & actor 2.7 (for 16–18 August 1970 magnetic storm) reduction in the daytime peak density due to enhanced vibrationally excited N₂ is brought about by the increase in the O⁺+N₂ rate factor.     

Potassium Doppler-resonance lidar for the study of the mesosphere and lower thermosphere at the Arecibo Observatory

We have developed a lidar to study the temperature structure of the nighttime mesopause region over the Arecibo Observatory (18.35°N, 66.75°W) by measuring the lineshape of the fluorescence spectrum of atomic potassium that is deposited in the mesosphere and lower thermosphere (MLT) by meteors. To demonstrate how the potassium lidar can enhance MLT studies at Arecibo, we show recent results for: (1) comparisons with airglow temperature measurements; (2) simultaneous operations with stratospheric and mesospheric temperature profiling by Rayleigh lidar; (3) simultaneous observations of K, Ca⁺, and E-region electron density profiles; and (4) occurrences of sporadic K layers, and relationships to sporadic E layers.     

Model of vibrational level populations of Herzberg states of oxygen molecules at heights of the lower thermosphere and mesosphere

The paper presents a model of the kinetics of electronically excited O₂(c¹Σ_u⁻,v), O₂(A′³Δ_u,v), O₂(A³Σ_u⁺,v) molecules at heights of the lower thermosphere and mesosphere with allowance for electronic excitation transfer processes during molecular collisions. The model is used to calculate the relative O₂(A³Σ_u⁺,v) and O₂(A′³Δ_u,v) populations at heights of 80–110 km. The calculated populations are compared with the available literature results on experimental estimates, and good agreement is obtained. It is shown how the increase in the quenching rates of the considered states by oxygen atoms affects the calculation results.     

IR band of O<Subscript>2</Subscript> at 1.27 μm as the tracer of O<Subscript>3</Subscript> in the mesosphere and lower thermosphere: Correction of the method

The problem of systematic overestimation (20–50%) of the retrieved ozone concentrations in the altitude range of 60–80 km in the TIMED–SABER satellite experiment in the daytime has been solved. The reason for overestimation is the neglect of the electronic vibrational kinetics of photolysis products of ozone and molecular oxygen O₂(b¹Σ_g ⁺, ν) and O₂(a¹Δ_g, ν). The IR emission band of O₂(a¹Δ_g, ν = 0) at 1.27 μm can be correctly used in remote sensing in order to obtain the ozone altitude profile in the altitude range of 50–88 km only with the use of a complete model of electronic vibrational kinetics of O₂ and O₃ photolysis products (YM2011) in the Earth’s mesosphere and lower thermosphere. Alternative ozone tracers have been considered, and an optimum tracer in the altitude range of 50–100 km such as O₂(b¹Σ_g ⁺, ν = 1) molecule emissions has been proposed.     

The green oxygen line 5577 Å in the mediumE-F andF regions during the night

Summary In this paper the contribution of theF region to the green oxygen line 5577 Å glow is investigated. It is shown that at calm nocturnal conditions the ions' O ₂ ⁺ dissociative recombination leads to a simultaneous excitation of the red and green lines in this region. It has been found that the intensity of the green line in theF region is about one-third that of the red line. By this non-negligible component we can explain some unexplained phenomena in the behaviour of the green line during the night.     

Global model of the thermosphere-ionosphere-protonosphere system

In this paper the formulation of the problem and preliminary numerical computation results of the thermosphere-ionosphere-protonosphere system parameters are discussed.The model constructed describes time-dependent distributions of the multicomponent near-earth space plasma parameters by means of numerical integration of the appropriate three-dimensional plasma hydrodynamic equations. In the thermospheric block of the model, global distribution of neutral gas temperature and N₂, O₂, O concentrations, as well as three-dimensional circulation of the neutral gas are calculated in the range of height from 80 km to 520 km. In the ionospheric section of the model, global time-dependent distribution of ion and electron temperatures, as well as molecular and atomic O⁺, H⁺ ion concentrations are calculated. Global two-dimensional distribution of electric potential is calculated taking into account computed thermosphere and ionosphere parameters.The inputs needed for our global model are the solar EUV spectrum; the auroral precipitation pattern; the distribution of the field-aligned currents and the model of the geomagnetic field.Preliminary results are obtained without regard to electromagnetic plasma drift for the solar minimum, low geomagnetic activity and spring equinox conditions. Global distributions of the calculated parameters in the magnetic dipole latitude-longitude frame are presented for 1200 UT. In the summary ignored processes and future direction are discussed.     

Lower thermospheric nitric oxide concentrations derived from WINDII observations of the green nightglow continuum at 553.1 nm

Vertical profiles of nitric oxide in the altitude range 90 to 105 km are derived from 553 nm nightglow continuum measurements made with the Wind Imaging Interferometer (WINDII) on the Upper Atmosphere Research Satellite (UARS). The profiles are derived under the assumption that the continuum emission is due entirely to the NO+O air afterglow reaction. Vertical profiles of the atomic oxygen density, which are required to determine the nitric oxide concentrations, are derived from coordinated WINDII measurements of the atomic oxygen OI 557.7 nm nightglow emission. Data coverage for local solar times ranging from 20 h to 04 h, and latitudes ranging from 42°S to 42°N, is achieved by zonally averaging and binning data obtained on 18 nights during a two-month period extending from mid-November 1992 until mid-January 1993. The derived nitric oxide concentrations are significantly smaller than those obtained from rocket measurements of the airglow continuum but they do compare well with model expectations and nitric oxide densities measured using the resonance fluorescence technique on the Solar Mesosphere Explorer satellite. The near-global coverage of the WINDII observations and the similarities to the nitric oxide global morphology established from other satellite measurements strongly suggests that the NO+O reaction is the major source of the continuum near 553 nm and that there is no compelling reason to invoke additional sources of continuum emission in this immediate spectral region.     

GCITEM-IGGCAS: A new global coupled ionosphere–thermosphere-electrodynamics model

The Global Coupled Ionosphere–Thermosphere-Electrodynamics Model developed at Institute of Geology and Geophysics, Chinese Academy of Sciences (GCITEM-IGGCAS), is introduced in this paper. This new model self-consistently calculates the time-dependent three-dimensional (3-D) structures of the main thermospheric and ionospheric parameters in the height range from 90 to 600 km, including neutral number density of major species O₂, N₂, and O and minor species N(²D), N(⁴S), NO, He and Ar; ion number densities of O⁺ ,O₂⁺, N₂⁺, NO⁺, N⁺ and electron; neutral, electron and ion temperature; and neutral wind vectors. The mid- and low-latitude electric fields can also be self-consistently calculated. GCITEM-IGGCAS is a full 3-D code with 5° latitude by 7.5° longitude cells in a spherical geographical coordinate system, which bases on an altitude grid. We show two simulations in this paper: a March Equinox one and a June Solstice one, and compare their simulation results to MSIS00 and IRI2000 empirical model. GCITEM-IGGCAS can reproduce the main features of the thermosphere and ionosphere in both cases.     

Anomalous electron density events in the quiet summer ionosphere at solar minimum over Millstone Hill

This study compares the observed behavior of the F region ionosphere over Millstone Hill with calculations from the IZMIRAN model for solar minimum for the geomagnetically quiet period 23–25 June 1986, when anomalously low values of hmF2(<200 km) were observed. We found that these low values of hmF2 (seen as a G condition on ionograms) exist in the ionosphere due to a decrease of production rates of oxygen ions resulting from low values of atomic oxygen density. Results show that determination of a G condition using incoherent scatter radar data is sensitive both to the true concentration of O⁺ relative to the molecular ions, and to the ion composition model assumed in the data reduction process. The increase in the O⁺ + N₂ loss rate due to vibrationally excited N₂ produces a reduction in NmF2 of typically 5–10%, but as large as 15%, bringing the model and data into better agreement. The effect of vibrationally excited NO⁺ ions on electron densities is negligible.     

Estimation of atomic oxygen concentrations from measured intensities of infrared nitric oxide radiation

The vibrational distribution of nitric oxide in the polar ionosphere computed according to the one-dimensional non-steady model of chemical and vibrational kinetics of the upper atmosphere has been compared with experimental data from rocket measurement. Some input parameters of the model have been varied to obtain the least-averaged deviation of the calculated population from experimental one. It is shown that the least deviation of our calculations from experimental measurements depends sufficiently on both the surprisal parameter of the production reaction of metastable atomic nitrogen with molecular oxygen and the profile of atomic oxygen concentration. The best agreement with the MSIS-83 profile was obtained for the value of surprisal parameter corresponding to recent laboratory estimations. The measured depression of level v = 2 is obtained in the calculation that uses sufficiently increased concentrations of atomic oxygen. It is pointed out that similar measurements of infrared radiation intensities could be used to estimate the atomic oxygen concentrations during auroral disturbances of the upper atmosphere.     

Non-thermal profile of the 557.7 nm emission with Doppler shift of dissociative component

A rare case, when non-thermal profiles of the [OI] 557.7 nm line with the dissociative components shifted relative to an ordinary Doppler kernel appear in auroras, is considered. Based on an analysis of these profiles, it has been indicated that the dissociative component is shifted because the electric field is present during the recombination of O ₂ ⁺ ion with background electrons in the ionospheric F region. The electric field component along the line of sight of the Fabry-Pérot interferometer (24 mV m^?1) has been estimated using the Doppler shift of the dissociative component of the 557.7 nm profile emission as an example.     

The morphology of oxygen greenline dayglow emission

In this study, the morphology of the oxygen greenline dayglow emission is presented. The volume emission rate profiles are obtained by using Solomons glow model. The glow model is updated in terms of recent cross sections, reaction rate coefficients and quantum yield of greenline emission. Throughout most of the thermosphere the modelled and observed emission rates are in reasonably good agreement. In the region between 98 and 120 km, the modelled emission rates are substantially higher (about a factor of 1.7) than the observed emission rates. This discrepancy is discussed in terms of scaling of solar fluxes which accounts the variation of solar activity for the day on which calculations are made. The modelled morphology of greenline emission is compared with those cases where WINDII data is available. The modelled and observed morphology is in reasonably good agreement at most of the latitudes above 120 km. In the mesosphere the qualitative nature of morphology is very similar to those of WINDII observation except the modelled emission rates are about a factor of 1.7 higher than the observed emission rates.