Anomalous variations in the structure of the ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region at geomagnetic midlatitudes of the Southern and Northern hemispheres in going from summer to winter conditions at high solar activity

The structure and dynamics of the ionosphere and plasmasphere at high solar activity under quiet geomagnetic conditions of June 2–3, 1979, and January 5–6, 1980, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The plasma drift velocity, determined by comparing the calculated and measured heights of the F ₂ layer maximum (hmF2), and the correction of [N₂] and [O₂], found in the NRLMSISE-00 model, make it possible to coordinate the electron densities (NmF2) calculated at the hmF2 height and the measured anomalous variations in NmF2 over the Argentine Islands ionosonde as well as the calculated and measured NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that, if the interference of the diffusion velocities of O⁺(⁴S) and H⁺ ions is taken into account, the additional heating of plasmaspheric electrons leads to an increase in the flux of O⁺(⁴S) ions from the topside ionosphere to lower F ₂ layer altitudes, as a result of which an anomalous nighttime increase in NmF2 6, observed on January 6, 1980, over Millstone Hill station, is mainly produced. The second component of the formation of anomalous night-time NmF2 is the plasma drift along the magnetic field caused by the neutral wind, which shifts O⁺(⁴S) ions to higher altitudes where the recombination rate of O⁺(⁴S) with N₂ and O₂ is lower and slows down a decrease in NmF2 in the course of time. It has been shown that the influence of electronically excited O⁺ ions and vibrationally excited N₂ and O₂ molecules on electron density (N _e ) considerably differs under winter and summer conditions. This difference forms significant part of the winter anomaly in N _e at heights of the F ₂ region and topside ionosphere over Millstone Hill station.     

High-speed correction factor to the O+-O resonance charge exchange collision frequency

The high-speed correction factor to the O⁺-O collision frequency, resulting from drift velocities between ions and neutrals, is calculated by solving the integral expression in this factor both numerically and analytically. Although the analytic solution is valid for either small or large drift velocities between ions and neutrals, for temperatures of interest and all drift velocities considered, agreement is found between analytic and detailed numerical integration results within less than 1% error. Let T _r designate the average of the ion and neutral temperatures in K, and u=v _d /<alpha>, where v _d is the relative drift velocity in cm s^?1, and <alpha>=4.56×10³\sqrtT _r cm s^?1 is the thermal velocity of the O⁺-O system. Then, as u ranges from 0 to 2, the correction factor multiplying the collision frequency increases monotonically from 1 to about 1.5. An interesting result emerging from this calculation is that the correction factor for temperatures of aeronomical interest is to a good approximation independent of the temperature, depending only on the scaled velocity u.     

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.     

Data modeling and assimilation studies with the MU radar

We report initial results of data modeling and assimilation studies for several MU radar experiments. Various inputs to a one-dimensional ionospheric model are adjusted to provide agreement with observation and also to learn the sensitivity of the model to their variations. Certain observations are also used directly in the model to anchor or constrain its behavior. In particular, studies of the electron density from 100 to 500 km altitude in the ionosphere are carried out with the help of a theoretical model of O⁺, NO⁺, O⁺₂ and N⁺₂ densities and MU radar observations of the power, ion-drift and plasma-temperature profiles. Four typical cases are selected to study quantitatively the effects of the (A) perpendicular-north component of the plasma drift (15 December 1986), (B) atmospheric composition (7 October 1986), (C) solar EUV flux (2 August 1989) and (D) upper-boundary O⁺ density (5 October 1989) on the model N_mF2, h_mF2 and N_e profile, as well as on the neutral wind calculation from h_mF2 and drift data. It is found that the measured vertical ion drift explains quantitatively well the measured h_mF2 (particularly at low solar activity) while the model gives a better match with the measured N_e when it uses the h_mF2-based wind rather than the measured plasma drift. Different model values of the atmospheric O/N₂ ratio and EUV flux and different values of the upper-bound O⁺ density may modify not only N_mF2 markedly but also h_mF2: a lower O/N₂ ratio results in higher h_mF2; the EUVAC model gives higher h_mF2 at high solar activity than does the EUV91 model; with a smaller upper-bound O⁺ density, h_mF2 is lower by day but little changed by night. We specifically note that the meridional wind needed by the model to reproduce the observed h_mF2 differed according to how well the model reproduced the observed N_mF2. The uncertainties in the MSIS86 and EUV model predictions are also discussed. It is found that if the MSIS and EUV91 models are used together, the model gives an N_mF2 higher than that measured at high solar activity. Thus the O/N₂ ratio needs to be reduced from the MSIS value if EUV91 is used. If EUVAC is used, no large modification is required. At equinox for low solar activity, modeling with either EUV model produces N_mF2 values lower than those measured, and so the true O/N₂ ratio may be higher than that given by MSIS model.     

Photochemistry of Ions at D-region Altitudes of the Ionosphere: A Review

The current state of knowledge of the D-region ion photochemistry is reviewed. Equations determining production rates of electrons and positive ions by photoionization of atmospheric neutral species are presented and briefly discussed. Considerable attention is given to the progress in the chemistry of O⁺(⁴S), O⁺(²D), O⁺(²P), N⁺, N₂ ⁺, O₂ ⁺, NO⁺, N₄ ⁺, O₄ ⁺, NO⁺(N₂), NO⁺(CO₂), NO⁺(CO₂)₂, NO⁺(H₂O) _n for n = 1–3, NO⁺(H₂O)(N₂), NO⁺(H₂O)₂(N₂), NO⁺(H₂O)(CO₂), NO⁺(H₂O)₂(CO₂), O₂ ⁺(H₂O), H₃O⁺(OH), H⁺(H₂O) _n for n = 1–8, O^?, O₂ ^?, O₃ ^?, O₄ ^?, OH^?, CO₃ ^?, CO₄ ^?, NO₂ ^?, NO₃ ^?, ONOO^?, Cl^?, Cl^?(H₂O), Cl^?(CO₂), HCO₃ ^?, CO₃ ^?(H₂O), CO₃ ^?(H₂O)₂, NO₃ ^?(H₂O), NO₃ ^?(H₂O)₂, OH^?(H₂O), and OH^?(H₂O)₂ ions. The analysis of the D-region rocket ion mass spectrometer measurements shows that, among these ions, O₂ ⁺, NO⁺, NO⁺(H₂O), and H⁺(H₂O) _n for n = 1–7 can make the main contribution to the total positive ion number density, and O^?, O₂ ^?, Cl^?, OH^?(H₂O), CO₃ ^?, HCO₃ ^?, NO₃ ^?, ONOO^?, CO₄ ^?, NO₃ ^?(H₂O), NO₃ ^?(H₂O)₂, and ³⁵Cl^?(CO₂) ions can be responsible for the main contribution to the total negative ion number density. Photodetachment of electrons from O^?, Cl^?, O₂ ^?, O₃ ^?, OH^?, NO₂ ^?, and NO₃ ^?, dissociative electron photodetachment of O₄ ^? and OH^?(H₂O), and photodissociation of O₃ ^?, O₄ ^?, CO₃ ^?, CO₄ ^?, ONOO^?, HCO₃ ^?, CO₃ ^?(H₂O), NO₃ ^?(H₂O), O₂ ⁺(H₂O), O₄ ⁺, N₄ ⁺, NO⁺(H₂O), NO⁺(H₂O)₂, H⁺(H₂O) _n for n = 2–4, NO⁺(N₂), and NO⁺(CO₂) are studied, and the photodetachment and photodissociation rate coefficients are calculated using the current state of knowledge on the cross sections of these processes and fluxes of solar radiation.     

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.     

Anomalous variations in the ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript>-layer structure at geomagnetic midlatitudes of the Southern and Northern hemispheres at the transition from summer to winter conditions under low solar activity

The structure and dynamics of the ionosphere and plasmasphere at low solar activity under quiet geomagnetic conditions on January 15–17, 1985, and July 10–13, 1986, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The detected correction of the model input parameters makes it possible to coordinate the measured and calculated anomalous variations in the electron density NmF2 at the height hmF2 of the ionospheric F2 layer over Argentine Islands ionosonde as well as the calculated and measured values of NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that vibrationally excited N₂ and O₂ molecules almost do not influence the formation of the winter anomaly under the conditions of low solar activity. A difference between the influence of electronically excited O⁺ on N _e ions under winter and summer conditions forms not more than 11% of the N _e winter anomaly event in the F ₂ layer and topside ionosphere. The model without electronically excited O⁺ ions reduces the duration of the N _e winter anomaly event. It has been shown that the seasonal variations in the composition of the neutral atmosphere form mainly the NmF2 winter anomaly event over the Millstone Hill radar at low solar activity.     

Ion Chemistry of the Ionosphere at E- and F-Region Altitudes: A Review

The current state of knowledge of E- and F-region ion chemistry is reviewed. Considerable attention is given to the progress in the chemistry of unexcited N₂ ⁺, O₂ ⁺, NO⁺, O⁺(⁴S), N⁺, H⁺, He⁺, Fe⁺, Mg⁺, Na⁺, Ca⁺, and K⁺ ions and electronically excited O⁺(²D), O⁺(²P), O⁺(⁴P), and $ {\text{O}}^{ + } (^{2} {\text{P}}^{*} ) $ ions. Achievements in our understanding of the role of vibrationally excited N₂ ⁺, O₂ ⁺, and NO⁺ ions in the ionosphere are discussed.     

Using tomography of GPS TEC to routinely determine ionospheric average electron density profiles

This paper introduces a technique that calculates average electron density (N_e) profiles over a wide geographic area of coverage, using tomographic ionospheric N_e profiles. These N_e profiles, which can provide information of the N_e distribution up to global positioning system (GPS) orbiting altitude (with the coordination of space-based GPS tomographic profiles), can be incorporated into the next generation of the international reference ionosphere (IRI) model. An additional advantage of tomography is that it enables accurate modeling of the topside ionosphere. By applying the tomographic reconstruction approach to ground-based GPS slant total electron content (STEC), we calculate 3-h average N_e profiles over a wide region. Since it uses real measurement data, tomographic average N_e profiles describe the ionosphere during quiet and disturbed periods. The computed average N_e profiles are compared with IRI model profiles and average N_e profiles obtained from ground-based ionosondes.     

Comparison of models and data at Millstone Hillduring the 5–11 June 1991 storm

We compare measurements of the ionospheric F region at Millstone Hillduring the severe geomagnetic disturbances of 5–11 June 1991 with results from the IZMIRANand FLIP time-dependent mathematical models of the Earths ionosphere and plasmasphere. Somecomparisons are also made with the Millstone Hill semi-empirical model which was previouslyused to model this storm. New rate coefficients from recent laboratory measurements of the O⁺+N₂ and O⁺+O₂ loss rates are included in theIZMIRAN and Millstone Hill models. The laboratory measurements show that vibrationallyexcited N₂ and O₂ (N₂(v) and O₂(v)) are both important at high temperatures such as found in the thermosphere during disturbedconditions at summer solar maximum. Increases in the O⁺+N₂ loss ratedue to N₂(v) result in a factor ∼2 reduction in the daytime F₂ peak electron density. On some days inclusion of N₂(v) improves theagreement between the models and the data, and on other days it worsens it. In the present workwe show for the first time the significant effect that the increase in the O⁺recombination rate due to O₂(v) may have on the calculated NmF₂. There are considerable uncertainties in the model calculations during the unusual,extremely disturbed conditions found during the daytime on 6 June. The results illustratedifficulties involved and the current state of the art in modelling severe disturbances, and thusprovide a benchmark against which future progress can be gauged.     

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.     

Variations of daytime and nighttime electron temperature and heat flux in the upper ionosphere,topside ionosphere and lower plasmasphere for low and high solar activity

A database of the electron temperature (T_e) comprising of most of the available LEO satellite measurements is used for studying the solar activity variations of T_e. The T_e data are grouped for two levels of solar activity (low LSA and high HSA), five altitude ranges between 350 and 2000 km, and day and night. By fitting a theoretical expression to the T_e values we obtain variation of T_e along magnetic field lines and heat flux for LSA and HSA. We have found that T_e increases with increase in solar activity at low and mid-latitudes during nighttime at all altitudes studied. During daytime the T_e response to solar activity depends on latitude, altitude, and season. This analysis shows existence of anti-correlation between T_e and solar activity at mid-latitudes below 700 km during the equinox and winter day hours. Heat fluxes show small latitudinal dependence for daytime but substantial for nighttime.     

The role of vibrationally excited nitrogen in the ionosphere

Theoretical and experimental aspects of the production, transformation, diffusion and loss of N₂ in the upper atmosphere are considered. The N₂-CO₂ near-resonant system in theD andE regions is taken into account. We describe our understanding of the methods necessary to find the vibrational populations of N₂ and CO₂ (asymmetric mode of CO₂). The calculations of the vibrational temperatures in theD, E, andF regions for the mid-latitude ionosphere and an aurora are presented. The connection between the excited species and the 4.26-m radiation intensities is considered. The models for the rate coefficient of the reaction of O⁺ with N₂ and the electron density decrease resulting from N₂ in the F region are discussed.     

Ionospheric disturbances during the severe magnetic storm of November 7–10, 2004

Results of the study of the behavior of the F ₂ region and topside ionosphere during the magnetic storm on November 7–10, 2004, which was a superposition of two sequent Severe magnetic disturbances (Kp = 9–) are presented. The observations were conducted by the incoherent scatter radar at Kharkov. Considerable effects of a negative ionospheric disturbance are registered, including a decrease in the electron density in the F ₂-layer maximum by a factor of 6–7 and of the total electron content up to a height of 1000 km by a factor of 2, a lifting up of the ionospheric F ₂ layer by 300 km at night and by 150–180 km in the daytime, unusual nighttime heating of the plasma with an increase of the ion and electron temperatures up to 2000 and 3000 K, respectively, and a decrease in the relative density of hydrogen ions N(H⁺)/N _e by a factor of up to 3.5 because of the emptying of the magnetic flux tube passing over Kharkov. The effects usually observed in the high-latitude ionosphere, including the coherent echoes, are detected during the main phase of the storm. The results obtained manifest a shift of the large-scale structures of the high-latitude ionosphere (the auroral oval, main ionospheric trough, hot zone, etc.) down to latitudes close to the latitude of the Kharkov radar.     

Electron loss and the determination of electron concentrations in theD-region

The most meaningful way to compare observations of the daytimeD-region under all ionization conditions, for the purpose of improving our understanding of this region, would appear to be through use of the effective electron recombination coefficient, =q/[e]², whereq is the ionization production rate, and where [e] is the electron concentration. This ratio apparently increases monotonically with decreasing altitude, and is much less variable than its componentsq and [e]. This parameter appears eminently suitable as a best first order solution to [e] after determiningq. For nighttime conditions, is more sensitive toq. However, an accurate but simple formula exists for the calculation of electron concentrations.     

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.     

The role of vibrationally excited oxygen and nitrogen in the ionosphere during the undisturbed and geomagnetic storm period of 6–12 April 1990

We present a comparison of the observed behavior of the F-region ionosphere over Millstone Hill during the geomagnetically quiet and storm periods of 6–12 April 1990 with numerical model calculations from the IZMIRAN time-dependent mathematical model of the Earths ionosphere and plasmasphere. The major enhancement to the IZMIRAN model developed in this study is the use of a new loss rate of O⁺(⁴S) ions as a result of new high-temperature flowing afterglow measurements of the rate coefficients K₁ and K₂ for the reactions of O⁺(⁴S) with N₂ and O₂. The deviations from the Boltzmann distribution for the first five vibrational levels of O₂(v) were calculated, and the present study suggests that these deviations are not significant. It was found that the difference between the non-Boltzmann and Boltzmann distribution assumptions of O₂(v) and the difference between ion and neutral temperature can lead to an increase of up to about 3% or a decrease of up to about 4% of the calculated NmF2 as a result of a respective increase or a decrease in K₂. The IZMIRAN model reproduces major features of the data. We found that the inclusion of vibrationally excited N₂(v > 0) and O₂(v > 0) in the calculations improves the agreement between the calculated NmF2 and the data on 6, 9, and 10 April. However, both the daytime and nighttime densities are reproduced by the IZMIRAN model without the vibrationally excited nitrogen and oxygen on 8 and 11 April better than the IZMIRAN model with N₂(v > 0) and O₂(v > 0). This could be due to possible uncertainties in model neutral temperature and densities, EUV fluxes, rate coefficients, and the flow of ionization between the ionosphere and plasmasphere, and possible horizontal divergence of the flux of ionization above the station. Our calculations show that the increase in the O⁺ + N₂ rate factor due to N₂(v > 0) produces a 5–36% decrease in the calculated daytime peak density. The increase in the O⁺ + O₂ loss rate due to vibrationally excited O₂ produces 8–46% reductions in NmF2. The effects of vibrationally excited O₂ and N₂ on N_e and T_e are most pronounced during the daytime.     

The role of vibrationally excited nitrogen and oxygen in the ionosphere over Millstone Hill during 16-23 March, 1990

We present a comparison of the observed behavior of the F region ionosphere over Millstone Hill during the geomagnetically quiet and storm period on 16/23 March, 1990, with numerical model calculations from the time-dependent mathematical model of the Earths ionosphere and plasmasphere. The effects of vibrationally excited N₂(v) and O₂(v) on the electron density and temperature are studied using the N₂(v) and O₂(v) Boltzmann and non-Boltzmann distribution assumptions. The deviations from the Boltzmann distribution for the first five vibrational levels of N₂(v) and O₂(v) were calculated. The present study suggests that these deviations are not significant at vibrational levels v = 1 and 2, and the calculated distributions of N₂(v) and O₂(v) are highly non-Boltzmann at vibrational levels v > 2. The N₂(v) and O₂(v) non-Boltzmann distribution assumption leads to the decrease of the calculated daytime NmF2 up to a factor of 1.44 (maximum value) in comparison with the N₂(v) and O₂(v) Boltzmann distribution assumption. The resulting effects of N₂(v > 0) and O₂(v) > 0) on the NmF2 is the decrease of the calculated daytime NmF2 up to a factor of 2.8 (maximum value) for Boltzmann populations of N₂(v) and O₂(v) and up to a factor of 3.5 (maximum value) for non-Boltzmann populations of N₂(v) and O₂(v). This decrease in electron density results in the increase of the calculated daytime electron temperature up to about 1040/1410 K (maximum value) at the F2 peak altitude giving closer agreement between the measured and modeled electron temperatures. Both the daytime and nighttime densities are not reproduced by the model without N₂(v > 0) and O₂(v > 0), and inclusion of vibrationally excited N2 and O2 brings the model and data into better agreement. The effects of vibrationally excited O2 and N2 on the electron density and temperature are most pronounced during daytime.