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.     

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.     

Physical mechanism of strong negative storm effects in the daytime ionospheric F2 region observed with EISCAT

A self-consistent method for daytime F-region modelling was applied to EISCAT observations during two periods comprising the very disturbed days 3 April 1992 and 10 April 1990. The observed strong N_e decrease at F2-layer heights originated from different physical mechanisms in the two cases. The negative F2-layer storm effect with an N_mF2 decrease by a factor of 6.4 on 3 April 1992 was produced by enhanced electric fields (E 85 mV/m) and strong downward plasma drifts, but without any noticeable changes in thermos-pheric parameters. The increase of the O⁺ + N₂ reaction rate resulted in a strong enrichment of the ionosphere with molecular ions even at F2-layer heights. The enhanced electric field produced a wide mid-latitude daytime trough on 03 April 1992 not usually observed during similar polarization jet events. The other strong negative storm effect on 10 April 1990 with a complete disappearance of the F2-layer maximum at the usual heights was attributed mainly to changes in neutral composition and temperature. A small value for the shape parameter S in the neutral temperature profile and a low neutral temperature at 120 km indicate strong cooling of the lower thermosphere. We propose that this cooling is due to increased nitric oxide concentration usually observed at these heights during geomagnetic storms.     

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.     

Localized structure in the cusp and high-latitude ionosphere: a modelling study

The ionospheric signature of a flux transfer event (FTE) seen in EISCAT radar data has been used as the basis for a modelling study using a new numerical model of the high-latitude ionosphere developed at the University of Sheffield, UK. The evolution of structure in the high-latitude ionosphere is investigated and examined with respect to the current views of polar patch formation and development. A localized velocity enhancement, of the type associated with FTEs, is added to the plasma as it passes through the cusp. This is found to produce a region of greatly enhanced ion temperature. The new model can provide greater detail during this event as it includes anisotropic temperature calculations for the O⁺ ions. This illustrates the uneven partitioning of the energy during an event of this type. O⁺ ion temperatures are found to become increasingly anisotropic, with the perpendicular temperature being substantially larger than the parallel component during the velocity enhancement. The enhanced temperatures lead to an increase in the recombination rate, which results in an alteration of the ion concentrations. A region of decreased O⁺ and increased molecular ion concentration develops in the cusp. The electron temperature is less enhanced than the ions. As the new model has an upper boundary of 10 000 km the topside can also be studied in great detail. Large upward fluxes are seen to transport plasma to higher altitudes, contributing to the alteration of the ion densities. Plasma is stored in the topside ionosphere and released several hours after the FTE has finished as the flux tube convects across the polar cap. This mechanism illustrates how concentration patches can be created on the dayside and be maintained into the nightside polar cap.     

Current state-of-the-art for the measurement of non-Maxwellian plasma parameters with the EISCAT UHF Facility

New results on the information that can be extracted from simulated non-Maxwellian incoherent radar spectra are presented. The cases of a pure ionosphere and of a composite ionosphere typical of a given altitude of the auroral F region are considered. In the case of a pure ionosphere of NO⁺ or O⁺ ions it has been shown that the electron temperature and the electron density can be derived from a Maxweilian analysis of radar spectra measured at aspect angles of 0° or 21° respectively; the ion temperature and ion temperature anisotropy can be derived from a non- constraining model such as the ID Raman fitting of a complementary measurement made at an aspect angle larger than 0° for the NO⁺ ions, or at an aspect angle larger than 21° for the O⁺ ions. Moreover with such measurements at large aspect angles, the shape of the velocity ion distribution functions can simultaneously be inferred. The case of a composite ionosphere of atomic O⁺ and molecular NO⁺ions is a difficult challenge which requires simultaneously a complementary measurement of the electron temperature to provide the ion composition and the electron density from the incoherent radar spectra at a specific aspect angle of 21°; hence, a model dependent routine is necessary to derive the ion temperatures and ion temperature anisotropies. In the case where the electron temperature is not given, a routine which depends on ion distribution models is required first: the better the ion distribution models are, the more accurately derived the plasma parameters will be. In both cases of a composite ionosphere, the 1D Raman fitting can be used to keep a check on the validity of the results provided by the ion distribution model dependent routine.     

Equatorial and low latitude ionosphere during intense geomagnetic storms

An investigation is made in order to analyse the role of neutral gas composition in the equatorial and low latitude ionosphere during intense geomagnetic storms. To this end data taken by the Dynamic Explorer 2 satellite at 280–300 km (molecular nitrogen N₂ and atomic oxygen O concentrations, electron density and vertical plasma drifts) are used. The sudden commencements of the events considered occurred at 11:38 UT on March 1, 1982, 18:41 UT on November 20, 1982 and 16:14 UT on February 4, 1983. Vertical plasma drifts are the most important contributor to the initial storm time response of the equatorial F region. Neutral composition changes (expressed as an increase in the molecular species, mainly N₂) possibly play a predominant role in the equatorial and low latitude (10–20°) decreases of electron density at heights near F2-region maximum during the main and recovery phases of intense geomagnetic storms. Delayed increases of electron density observed at daytime during the recovery phase may be also attributed to increases in atomic oxygen. At low latitudes possibly a combined effect of O increase and upward plasma drift due to enhanced equatorward winds is the responsible mechanism for the maintenance of enhanced electron density values.     

Vertical circulation and thermospheric composition: a modelling study

The coupled thermosphere-ionosphere-plasmasphere model CTIP is used to study the global three-dimensional circulation and its effect on neutral composition in the midlatitude F-layer. At equinox, the vertical air motion is basically up by day, down by night, and the atomic oxygen/molecular nitrogen [O/N₂] concentration ratio is symmetrical about the equator. At solstice there is a summer-to-winter flow of air, with downwelling at subauroral latitudes in winter that produces regions of large [O/N₂] ratio. Because the thermospheric circulation is influenced by the high-latitude energy inputs, which are related to the geometry of the Earth’s magnetic field, the latitude of the downwelling regions varies with longitude. The downwelling regions give rise to large F2-layer electron densities when they are sunlit, but not when they are in darkness, with implications for the distribution of seasonal and semiannual variations of the F2-layer. It is also found that the vertical distributions of O and N₂ may depart appreciably from diffusive equilibrium at heights up to about 160 km, especially in the summer hemisphere where there is strong upwelling.     

Comparison of the measured and modelled electron densities and temperatures in the ionosphere and plasmasphere during 20/30 January, 1993

We present a comparison of the electron density and temperature behaviour in the ionosphere and plasmasphere measured by the Millstone Hill incoherent-scatter radar and the instruments on board of the EXOS-D satellite with numerical model calculations from a time-dependent mathematical model of the Earths ionosphere and plasmasphere during the geomagnetically quiet and storm period on 20/30 January, 1993. We have evaluated the value of the additional heating rate that should be added to the normal photoelectron heating in the electron energy equation in the daytime plasmasphere region above 5000 km along the magnetic field line to explain the high electron temperature measured by the instruments on board of the EXOS-D satellite within the Millstone Hill magnetic field flux tube in the Northern Hemisphere. The additional heating brings the measured and modelled electron temperatures into agreement in the plasmasphere and into very large disagreement in the ionosphere if the classical electron heat flux along magnetic field line is used in the model. A new approach, based on a new effective electron thermal conductivity coefficient along the magnetic field line, is presented to model the electron temperature in the ionosphere and plasmasphere. This new approach leads to a heat flux which is less than that given by the classical Spitzer-Harm theory. The evaluated additional heating of electrons in the plasmasphere and the decrease of the thermal conductivity in the topside ionosphere and the greater part of the plasmasphere found for the first time here allow the model to accurately reproduce the electron temperatures observed by the instruments on board the EXOS-D satellite in the plasmasphere and the Millstone Hill incoherent-scatter radar in the ionosphere. The effects of the daytime additional plasmaspheric heating of electrons on the electron temperature and density are small at the F-region altitudes if the modified electron heat flux is used. 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 the first vibrational levels of N₂ and O₂ and the second level of O₂, and the calculated distributions of N₂(v) and O₂(v) are highly non-Boltzmann at vibrational levels v > 2. The resulting effect of N₂(v > 0) and O₂(v > 0) on NmF2 is the decrease of the calculated daytime NmF2 up to a factor of 1.5. The modelled electron temperature is very sensitive to the electron density, and this decrease in electron density results in the increase of the calculated daytime electron temperature up to about 580 K at the F2 peak altitude giving closer agreement between the measured and modelled electron temperatures. Both the daytime and night-time densities are not reproduced by the model without N₂(v > 0) and O₂(v > 0), and inclusion of vibrationally excited N₂ and O₂ brings the model and data into better agreement.     

Local enhancements in the horizontal distribution of the electron density at heights of the topside ionosphere near the dayside polar cap boundary

The structure of the electron density horizontal distribution at heights of the topside ionosphere in the Northern Hemisphere under quiet magnetic conditions in the polar peak region has been analyzed based on the measurements conducted on board the STSAT-1, DMSP F13, and DMSP F15 satellites during the period of moderate solar activity. The sector near 1000 MLT (magnetic local time) for the spring equinox season has been studied most thoroughly. It has been obtained that the polar peak consists of localized irregularities in the electron density of a specific form, inside which the electron density exceeds its background values by several tens of percent. The characteristic dimensions of the irregularities vary from a few degrees to several tens of degrees in the meridional and zonal directions, respectively. The irregularities are centered along the dayside boundary of the polar cap and coincide with the region of “soft” electron precipitation (with the average energy lower than 500 eV and with the flux density above 10⁷ electrons cm^?2 s^?1 sr^?1).     

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.     

Annual and semiannual variations in the ionospheric F2-layer: II. Physical discussion

The companion paper by Zou et al. shows that the annual and semiannual variations in the peak F2-layer electron density (NmF2) at midlatitudes can be reproduced by a coupled thermosphere-ionosphere computational model (CTIP), without recourse to external influences such as the solar wind, or waves and tides originating in the lower atmosphere. The present work discusses the physics in greater detail. It shows that noon NmF2 is closely related to the ambient atomic/molecular concentration ratio, and suggests that the variations of NmF2 with geographic and magnetic longitude are largely due to the geometry of the auroral ovals. It also concludes that electric fields play no important part in the dynamics of the midlatitude thermosphere. Our modelling leads to the following picture of the global three-dimensional thermospheric circulation which, as envisaged by Duncan, is the key to explaining the F2-layer variations. At solstice, the almost continuous solar input at high summer latitudes drives a prevailing summer-to-winter wind, with upwelling at low latitudes and throughout most of the summer hemisphere, and a zone of downwelling in the winter hemisphere, just equatorward of the auroral oval. These motions affect thermospheric composition more than do the alternating day/night (up-and-down) motions at equinox. As a result, the thermosphere as a whole is more molecular at solstice than at equinox. Taken in conjunction with the well-known relation of F2-layer electron density to the atomic/molecular ratio in the neutral air, this explains the F2-layer semiannual effect in NmF2 that prevails at low and middle latitudes. At higher midlatitudes, the seasonal behaviour depends on the geographic latitude of the winter downwelling zone, though the effect of the composition changes is modified by the large solar zenith angle at midwinter. The zenith angle effect is especially important in longitudes far from the magnetic poles. Here, the downwelling occurs at high geographic latitudes, where the zenith angle effect becomes overwhelming and causes a midwinter depression of electron density, despite the enhanced atomic/molecular ratio. This leads to a semiannual variation of NmF2. A different situation exists in winter at longitudes near the magnetic poles, where the downwelling occurs at relatively low geographic latitudes so that solar radiation is strong enough to produce large values of NmF2. This circulation-driven mechanism provides a reasonably complete explanation of the observed pattern of F2 layer annual and semiannual quiet-day variations.     

The formation of electron-temperature hot spots in the main ionospheric trough by the internal processes

A mathematical model of the convecting high-latitude ionosphere is described which produces three-dimensional distributions of electron density, positive-ion velocity and electron and ion temperatures at the F-layer altitudes. The results of simulation of the behaviour of the high-latitude ionosphere, in particular, the heat regime of the F-layer, are presented and analysed. From our study, it was found that electron-temperature hot spots in the main ionospheric trough can arise owing to internal ionospheric processes, and not due to effects of any external causes. Three conditions, to be satisfied simultaneously, are necessary for the formation of the considered electron-temperature hot spots: first, low values of electron density; second, solar illumination of the upper F region and darkness of the lower F region; third, low values of neutral-component densities. These conditions are valid in the main ionospheric trough near the terminator on the nightside when the density of the neutral atmosphere is not high. The physical processes which lead to the formation of the electron-temperature hot spots are the heat transfer from the upper into the lower F region, the reduced heat capacity of electron gas and the weakened cooling of electron gas due to inelastic collisions with neutral atoms and molecules. Also investigated is the influence of seasonal and solar-activity variations on the efficiency of the identified mechanism responsible for the formation of the electron temperature peaks in the main ionospheric trough by the internal processes.     

Variations in ionospheric parameters during solar flares

Results of studying the ionospheric response to solar flares, obtained based on the incoherent scatter radar observations of the GPS signals and as a result of the model simulations, are presented. The method, based on the effect of partial “shadowing” of the atmosphere by the globe, has been used to analyze the GPS data. This method made it possible to estimate the value of a change in the electron content in the upper ionosphere during the solar flare of July 14, 2000. It has been shown that a flare can cause a decrease in the electron content at heights of the upper ionosphere (h > 300 km) according to the GPS data. Similar effects in the formation of a negative disturbance in the ionospheric F region were also observed during the solar flares of May 21 and 23, 1967, at the Arecibo incoherent scatter radar. The mechanism by which negative disturbances are formed in the upper ionosphere during solar flares has been studied based on the theoretical model of the ionosphere-plasmasphere coupling. It has been shown that an intense ejection of O⁺ ions into the above located plasmasphere under the action of a sharp increase in the ion production rate and the thermal expansion of the ionospheric plasma cause the formation of a negative disturbance in the electron concentration in the upper ionosphere.     

Configuration of the upper boundary of the ionosphere

Variations of the upper boundary of the ionosphere (UBI) are investigated based on three sources of information: (i) ionosonde-derived parameters: critical frequency foF2, propagation factor M3000F2, and sub-peak thickness of the bottomside electron density profile; (ii) total electron content (TEC) observations from signals of the Global Positioning System (GPS) satellites; (iii) model electron densities of the International Reference Ionosphere (IRI^*) extended towards the plasmasphere. The ionospheric slab thickness is calculated as ratio of TEC to the F2 layer peak electron density, NmF2, representing a measure of thickness of electron density profile in the bottomside and topside ionosphere eliminating the plasmaspheric slab thickness of GPS-TEC with the IRI^* code. The ratio of slab thickness to the real thickness in the topside ionosphere is deduced making use of a similar ratio in the bottomside ionosphere with a weight Rw. Model weight Rw is represented as a superposition of the base-functions of local time, geomagnetic latitude, solar and magnetic activity. The time-space variations of domain of convergence of the ionosphere and plasmasphere differ from an average value of UBI at ∼1000 km over the earth. Analysis for quiet monthly average conditions and during the storms (September 2002, October–November 2003, November 2004) has shown shrinking UBI altitude at daytime to 400 km. The upper ionosphere height is increased by night with an ‘ionospheric tail’ which expands from 1000 km to more than 2000 km over the earth under quiet and disturbed space weather. These effects are interposed on a trend of increasing UBI height with solar activity when both the critical frequency foF2 and the peak height hmF2 are growing during the solar cycle.     

ELF wave generation in the ionosphere using pulse modulated HF heating: initial tests of a technique for increasing ELF wave generation efficiency

This paper describes the results of a preliminary study to determine the effective heating and cooling time constants of ionospheric currents in a simulated modulated HF heating, ‘beam painting’ configuration. It has been found that even and odd harmonics of the fundamental ELF wave used to amplitude modulate the HF heater are sourced from different regions of the ionosphere which support significantly different heating and cooling time constants. The fundamental frequency and its odd harmonics are sourced in a region of the ionosphere where the heating and cooling time constants are about equal. The even harmonics on the other hand are sourced from regions of the ionosphere characterised by ratios of cooling to heating time constant greater than ten. It is thought that the even harmonics are sourced in the lower ionosphere (around 65 km) where the currents are much smaller than at the higher altitudes around 78 km where the currents at the fundamental frequency and odd harmonics maximise.     

Formation Mechanisms of the Spring–Autumn Asymmetry of the Midlatitudinal <Emphasis Type="Italic">NmF</Emphasis>2 under Daytime Quiet Geomagnetic Conditions at Low Solar Activity

Formation mechanism of the spring–autumn asymmetry of the F2-layer peak electron number density of the midlatitudinal ionosphere, NmF2, under daytime quiet geomagnetic conditions at low solar activity are studied. We used the ionospheric parameters measured by the ionosonde and incoherent scatter radar at Millstone Hill on March 3, 2007, March 29, 2007, September 12, 2007, and September 18, 1984. The altitudinal profiles of the electron density and temperature were calculated for the studied conditions using a one-dimensional, nonstationary, ionosphere–plasmasphere theoretical model for middle geomagnetic latitudes. The study has shown that there are two main factors contributing to the formation of the observed spring–autumn asymmetry of NmF2: first, the spring–autumn variations of the plasma drift along the geomagnetic field due to the corresponding variations in the components of the neutral wind velocity, and, second, the difference between the composition of the neutral atmosphere under the spring and autumn conditions at the same values of the universal time and the ionospheric F2-layer peak altitude. The seasonal variations of the rate of O⁺(⁴S) ion production, which are associated with chemical reactions with the participation of the electronically excited ions of atomic oxygen, does not significantly affect the studied NmF2 asymmetry. The difference in the degree of influence of O⁺(⁴S) ion reactions with vibrationally excited N₂ and O₂ on NmF2 under spring and autumn conditions does not significantly change the spring–autumn asymmetry of NmF2.     

化学物质释放激发中低纬扩展F的数值模拟

利用离子和电子动量方程、连续性方程以及电流连续性方程建立了适合描述中低纬spread-F发展的物理模型,并对模型进行了数值求解,讨论了利用H₂O释放来激发电离层Rayleigh-Taylor不稳定性的可能性. 结果表明,电离层处于不稳定状态时,H₂O在电离层底部释放后,造成电子的大量消耗,增强了峰值高度以下电子的密度梯度,有利于spread-F的发展,在spread-F发展的过程中,释放中心附近会形成电子密度的消耗区,两侧出现密度的增强区;而电离层比较稳定时,初始扰动会逐渐稳定下来,但化学物质释放仍能造成电子较长时间、较大范围的扰动.     

A self-consistent estimate of O+ + N2– rate coefficient and total EUV solar flux with λ < 1050 Å using EISCAT observations

There are differences between existing models of solar EUV with < 1050 Å and between laboratory measurements of the O⁺ + N₂ – reaction rate coefficient, both parameters being crucial for the F2-region modeling. Therefore, indirect aeronomic estimates of these parameters may be useful for qualifying the existing EUV models and the laboratory measured O⁺ + N₂ – rate coefficient. A modified self-consistent method for daytime F2-region modeling developed by Mikhailov and Schlegel was applied to EISCAT observations (32 quiet summer and equinoctial days) to estimate the set of main aeronomic parameters. Three laboratory measured temperature dependencies for the O⁺ + N₂ – rate coefficient were used in our calculations to find self-consistent factors both for this rate coefficient and for the solar EUV flux model from Nusinov. Independent of the rate coefficient used, the calculated values group around the temperature dependence recently measured by Hierl et al. in the 850–1400 K temperature range. Therefore, this rate coefficient may be considered as the most preferable and is recommended for aeronomic calculations. The calculated EUV flux shows a somewhat steeper dependence on solar activity than both, the Nusinov and the EUVAC models predict. In practice both EUV models may be recommended for the F2-region electron density calculations with the total EUV flux shifted by ±25% for the EUVAC and Nusinov models, correspondingly.     

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.