Solar activity effects in the ionospheric D region

Variations in the D-region electron concentration within the solar activity cycle are considered. It is demonstrated that conclusions of various authors, who have analyzed various sets of experimental data on [e], differ significantly. The most reliable seem to be the conclusions based on analysis of the [e] measurements carried out by the Faraday rotation method and on the theoretical concepts on the D-region photochemistry. Possible QBO effects in the relation of [e] to solar activity are considered and an assumption is made that such effects may be the reason for the aforementioned disagreement in conclusions on the [e] relation to solar indices.Paper Presented at the Second IAGA/ICMA (IAMAS) Workshop on Solar Activity Forcing of the Middle Atmosphere, Prague, August 1997     

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.     

Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model

Current theories of F-layer storms are discussed using numerical simulations with the Upper Atmosphere Model, a global self-consistent, time dependent numerical model of the thermosphere-ionosphere-plasmasphere-magnetosphere system including electrodynamical coupling effects. A case study of a moderate geomagnetic storm at low solar activity during the northern winter solstice exemplifies the complex storm phenomena. The study focuses on positive ionospheric storm effects in relation to thermospheric disturbances in general and thermospheric composition changes in particular. It investigates the dynamical effects of both neutral meridional winds and electric fields caused by the disturbance dynamo effect. The penetration of short-time electric fields of magnetospheric origin during storm intensification phases is shown for the first time in this model study. Comparisons of the calculated thermospheric composition changes with satellite observations of AE-C and ESRO-4 during storm time show a good agreement. The empirical MSISE90 model, however, is less consistent with the simulations. It does not show the equatorward propagation of the disturbances and predicts that they have a gentler latitudinal gradient. Both theoretical and experimental data reveal that although the ratio of [O]/[N₂] at high latitudes decreases significantly during the magnetic storm compared with the quiet time level, at mid to low latitudes it does not increase (at fixed altitudes) above the quiet reference level. Meanwhile, the ionospheric storm is positive there. We conclude that the positive phase of the ionospheric storm is mainly due to uplifting of ionospheric F2-region plasma at mid latitudes and its equatorward movement at low latitudes along geomagnetic field lines caused by large-scale neutral wind circulation and the passage of travelling atmospheric disturbances (TADs). The calculated zonal electric field disturbances also help to create the positive ionospheric disturbances both at middle and low latitudes. Minor contributions arise from the general density enhancement of all constituents during geomagnetic storms, which favours ion production processes above ion losses at fixed height under day-light conditions.     

Semiannual and annual variations in the height of the ionospheric F2-peak

Ionosonde data from sixteen stations are used to study the semiannual and annual variations in the height of the ionospheric F2-peak, hmF2. The semiannual variation, which peaks shortly after equinox, has an amplitude of about 8 km at an average level of solar activity (10.7 cm flux = 140 units), both at noon and midnight. The annual variation has an amplitude of about 11 km at northern midlatitudes, peaking in early summer; and is larger at southern stations, where it peaks in late summer. Both annual and semiannual amplitudes increase with increasing solar activity by day, but not at night. The semiannual variation in hmF2 is unrelated to the semiannual variation of the peak electron density NmF2, and is not reproduced by the CTIP and TIME-GCM computational models of the quiet-day thermosphere and ionosphere. The semiannual variation in hmF2 is approximately isobaric, in that its amplitude corresponds quite well to the semiannual variation in the height of fixed pressure-levels in the thermosphere, as represented by the MSIS empirical model. The annual variation is not isobaric. The annual mean of hmF2 increases with solar 10.7 cm flux, both by night and by day, on average by about 0.45 km/flux unit, rather smaller than the corresponding increase of height of constant pressure-levels in the MSIS model. The discrepancy may be due to solar-cycle variations of thermospheric winds. Although geomagnetic activity, which affects thermospheric density and temperature and therefore hmF2 also, is greatest at the equinoxes, this seems to account for less than half the semiannual variation of hmF2. The rest may be due to a semiannual variation of tidal and wave energy transmitted to the thermosphere from lower levels in the atmosphere.     

Annual and semiannual variations in the ionospheric F2-layer. I. Modelling

Annual, seasonal and semiannual variations of F2-layer electron density (NmF2) and height (hmF2) have been compared with the coupled thermosphere-ionosphere-plasmasphere computational model (CTIP), for geomagnetically quiet conditions. Compared with results from ionosonde data from midlatitudes, CTIP reproduces quite well many observed features of NmF2, such as the dominant winter maxima at high midlatitudes in longitude sectors near the magnetic poles, the equinox maxima in sectors remote from the magnetic poles and at lower latitudes generally, and the form of the month-to-month variations at latitudes between about 60°N and 50°S. CTIP also reproduces the seasonal behaviour of NmF2 at midnight and the summer-winter changes of hmF2. Some features of the F2-layer, not reproduced by the present version of CTIP, are attributed to processes not included in the modelling. Examples are the increased prevalence of the winter maxima of noon NmF2 at higher solar activity, which may be a consequence of the increase of F2-layer loss rate in summer by vibrationally excited molecular nitrogen, and the semiannual variation in hmF2, which may be due to tidal effects. An unexpected feature of the computed distributions of NmF2 is an east-west hemisphere difference, which seems to be linked to the geomagnetic field configuration. Physical discussion is reserved to the companion paper by Rishbeth et al.     

New electron energy transfer rates for vibrational excitation of N2

In this paper we present the results of a study of the electron cooling rate, the production rates of vibrationally excited N₂(v), and the production frequency of the N₂ vibrational quanta arising from the collisions of electrons with unexcited N₂(0) and vibrationally excited N₂(1) molecules as functions of the electron temperature. The electron energy transfer rates for vibrational excitation of N₂ have been calculated and fit to analytical expressions by use of the revised vibrationally excited N₂ cross sections. These new analytical expressions are available to the researcher for quick reference and accurate computer modeling with a minimum of calculations.     

Trends in the ionospheric E and F regions over Europe

Continuous observations in the ionospheric E and F regions have been regularly carried out since the fifties of this century at many ionosonde stations. Using these data from 31 European stations long-term trends have been derived for different parameters of the ionospheric E layer (h E, foE), F1 layer (foF1) and F2 layer (hmF2,foF2). The detected trends in the E and F1 layers (lowering of the E region height hE; increase of the peak electron densities of the E and F1 layers, foE and foF1) are in qualitative agreement with model predictions of an increasing atmospheric greenhouse effect. In the F2 region, however, the results are more complex. Whereas in the European region west of 30° E negative trends in hmF2 (peak height of the F2 layer) and in the peak electron density (foF2) have been found, in the eastern part of Europe (east of 30° E) positive trends dominate in both parameters. These marked longitudinal differences cannot be explained by an increasing greenhouse effect only, here probably dynamical effects in the F2 layer seem to play an essential role.     

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.     

On the onset and meridional propagation of the ionospheric F2-region response to geomagnetic storms

The meridional propagation velocities of the ionospheric F2-region response to 268 geomagnetic storms are calculated. Ionospheric vertical sounding data of 1 h time resolution from several stations located in a longitude sector approximately centred along the great circle that contains both the geomagnetic poles and the geographic poles are used.Most meridional propagation velocities from high to low latitudes are less than 600 m/s. The smaller velocities are typical of global neutral meridional wind circulation and the larger are representative of traveling atmospheric disturbances.Simultaneous disturbances at several locations are more frequent during positive phases than during negative phases. Negative phase meridional propagation velocities associated with meridional neutral winds are less frequent in the southern hemisphere when compared with corresponding velocities observed in the northern hemisphere. This may be related to the fact that the distance between the geomagnetic pole and the equator is smaller in the northern hemisphere.Most negative phase onsets are within the 06–10 LT interval. For middle geomagnetic latitudes a “forbidden time interval” between 11 and 14 LT is present. The positive phase onsets show the “dusk effect”.     

Global transport and localized layering of metallic ions in the upper atmosphere

A numerical model has been developed which is capable of simulating all phases of the life cycle of metallic ions, and results are described and interpreted herein for the typical case of Fe⁺ ions. This cycle begins with the initial deposition of metallics through meteor ablation and sputtering, followed by conversion of neutral Fe atoms to ions through photoionization and charge exchange with ambient ions. Global transport arising from daytime electric fields and poleward/downward diffusion along geomagnetic field lines, localized transport and layer formation through descending convergent nulls in the thermospheric wind field, and finally annihilation by chemical neutralization and compound formation are treated. The model thus sheds new light on the interdependencies of the physical and chemical processes affecting atmospheric metallics. Model output analysis confirms the dominant role of both global and local transport to the ions life cycle, showing that upward forcing from the equatorial electric field is critical to global movement, and that diurnal and semidiurnal tidal winds are responsible for the formation of dense ion layers in the 90–250 km height region. It is demonstrated that the assumed combination of sources, chemical sinks, and transport mechanisms actually produces F-region densities and E-region layer densities similar to those observed. The model also shows that zonal and meridional winds and electric fields each play distinct roles in local transport, whereas the ion distribution is relatively insensitive to reasonable variations in meteoric deposition and chemical reaction rates.     

On the fast movements of the thermosphere: an interpretation

We examine the oscillations of the meridional neutral wind in the F region as seen by the EISCAT radar. We propose an interpretation in term of eddies (tourbillons) of typical size of a few tens to a few hundreds of kilometers. The observed rotation velocity is a few hundreds of meters per second. We suggest that the tourbillons are a common feature of thermospheric movements. We propose an optical experiment to check the validity of this assumption.     

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.     

First observation of quasi-2-day oscillations in ionospheric plasma frequency at fixed heights

The existence and development of the quasi-2-day oscillations in the plasma frequency variations of the F region at northern middle latitudes are investigated. A new approach to study the quasi-2-day oscillations is presented, using a methodology that allows us to do such a study at fixed heights. The hourly values of plasma frequency at fixed heights, from 170 km to 220 km at 10 km step, obtained at the Observatori de lEbre station (40.8°N, 0.5°E) during 1995 are used for analysis. It is found that quasi-2-day oscillations exist and persisted in the ionospheric plasma frequency variations over the entire year 1995 for all altitudes investigated. The dominant period of oscillation ranges from 42 to 56 h. The amplitude of oscillation is from 0.1 MHz to 1 MHz. The activity of the quasi-2-day oscillation is better expressed during the summer half year when several enhancements, about 15–30 days in duration, were observed. The largest enhancements of the oscillation occurred during early June, July and early August; i. e., near and after the summer solstice when the 2-day wave in the middle neutral atmosphere typically displays its largest activity in the Northern Hemisphere. The results obtained may help us understand better the possible influencing mechanisms between the 2-day wave in the middle neutral atmosphere and the ionospheric quasi-2-day oscillations.     

Combined CUTLASS, EISCAT and ESR observations of ionospheric plasma flows at the onset of an isolated substorm

On August 21st 1998, a sharp southward turning of the IMF, following on from a 20 h period of northward directed magnetic field, resulted in an isolated substorm over northern Scandinavia and Svalbard. A combination of high time resolution and large spatial scale measurements from an array of coherent scatter and incoherent scatter ionospheric radars, ground magnetometers and the Polar UVI imager has allowed the electrodynamics of the impulsive substorm electrojet region during its first few minutes of evolution at the expansion phase onset to be studied in great detail. At the expansion phase onset the substorm onset region is characterised by a strong enhancement of the electron temperature and UV aurora. This poleward expanding auroral structure moves initially at 0.9 km s^-1 poleward, finally reaching a latitude of 72.5°. The optical signature expands rapidly westwards at ~6 km s^-1, whilst the eastward edge also expands eastward at ~0.6 km s^-1. Typical flows of 600 m s^-1 and conductances of 2 S were measured before the auroral activation, which rapidly changed to ~100 m s^-1 and 10–20 S respectively at activation. The initial flow response to the substorm expansion phase onset is a flow suppression, observed up to some 300 km poleward of the initial region of auroral luminosity, imposed over a time scale of less than 10 s. The high conductivity region of the electrojet acts as an obstacle to the flow, resulting in a region of low-electric field, but also low conductivity poleward of the high-conductivity region. Rapid flows are observed at the edge of the high-conductivity region, and subsequently the high flow region develops, flowing around the expanding auroral feature in a direction determined by the flow pattern prevailing before the substorm intensification. The enhanced electron temperatures associated with the substorm-disturbed region extended some 2° further poleward than the UV auroral signature associated with it.     

The role of water-vapour photodissociation on the formation of a deep minimum in mesopause ozone

A one-dimensional atmospheric photochemical model with an altitude grid of about 1.5 km was used to examine the structure of the global mean vertical ozone profile and its night-time-to-daytime variation in the upper atmosphere. Two distinct ozone layers are predicted, separated by a sharp drop in the ozone concentration near the mesopause. This naturally occurring mesopause ozone deep minimum is primarily produced by the rapid increase in the destruction of water vapour, and hence increase in HO_x, at altitudes between 80 and 85 km, a region where water-vapour photodissociation by ultraviolet radiation of the solar Lyman-alpha line is significant, and where the supply of water vapour is maintained by methane oxidation even for very dry conditions at the tropospheric-stratospheric exchange region. The model indicates that the depth of the mesopause ozone minimum is limited by the efficiency with which inactive molecular hydrogen is produced, either by the conversion of atomic hydrogen to molecular hydrogen via one of the reaction channels of H with HO₂, or by Lyman-alpha photodissociation of water vapour via the channel that leads to the production of molecular hydrogen. The ozone concentration rapidly recovers above 85 km due to the rapid increase in O produced by the photodissociation of O₂ by absorption of ultraviolet solar radiation in the Schumann-Runge bands and continuum. Above 90 km, there is a decrease in ozone due to photolysis as the production of ozone through the three-body recombination of O₂ and O becomes slower with decreasing pressure. The model also predicts two peaks in the night-time/daytime ozone ratio, one near 75 km and the other near 110 km, plus a strong peak in the night-time/daytime ratio of OH near 110 km. Recent observational evidence supports the predictions of the model.     

A distribution law for relative humidity in the upper troposphere and lower stratosphere derived from three years of MOZAIC measurements

Data from three years of MOZAIC measurements made it possible to determine a distribution law for the relative humidity in the upper troposphere and lower stratosphere. Data amounting to 13.5% of the total were obtained in regions with ice supersaturation. Troposphere and stratosphere are distinguished by an ozone concentration of 130 ppbv as threshold. The probability of measuring a certain amount of ice supersaturation in the troposphere decreases exponentially with the degree of ice supersaturation. The probability of measuring a certain relative humidity in the stratosphere (both with respect to water and ice) decreases exponentially with the relative humidity. A stochastic model that naturally leads to the exponential distribution is provided. Mean supersaturation in the troposphere is about 15%, whereas ice nucleation requires 30% supersaturation on the average. This explains the frequency of regions in which aircraft induce persistent contrails but which are otherwise free of clouds. Ice supersaturated regions are 3-4 K colder and contain more than 50% more vapour than other regions in the upper troposphere. The stratospheric air masses sampled are dry, as expected, having mean relative humidity over water of 12% and over ice of 23%, respectively. However, 2% of the stratospheric data indicate ice supersaturation. As the MOZAIC measurements have been obtained on commercial flights mainly between Europe and North America, the data do not provide a complete global picture, but the exponential character of the distribution laws found is probably valid globally. Since water vapour is the most important greenhouse gas and since it might enhance the anthropogenic greenhouse effects via positive feedback mechanisms, it is important to represent its distribution correctly in climate models. The discovery of the distribution law of the relative humidity makes possible simple tests to show whether the hydrological cycle in climate models is represented in an adequate way or not.     

Case studies of the propagation characteristics of auroral TIDS with EISCAT CP2 data using maximum entropy cross-spectral analysis

In this paper case studies of propagation characteristics of two TIDs are presented which are induced by atmospheric gravity waves in the auroral F-region on a magnetic quiet day. By means of maximum entropy cross-spectral analysis of EISCAT CP2 data, apparent full wave-number vectors of the TIDs are obtained as a function of height. The analysis results show that the two events considered can be classified as moderately large-scale TID and medium-scale TID, respectively. One exhibits a dominant period of about 72 min, a mean horizontal phase speed of about 180 m/s (corresponding to a horizontal wavelength of about 780 km) directed south-eastwards and a vertical phase speed of 55 m/s for a height of about 300 km. The other example shows a dominant period of 44 min, a mean horizontal phase velocity of about 160 m/s (corresponding to a horizontal wavelength of about 420 km) directed southwestwards, and a vertical phase velocity of about 50 m/s at 250 km altitude.Max-Planck-Institut für Aeronomie, Germany. On leave from Department of Space Physics, Wuhan University, China     

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.     

曲靖非相干散射雷达电离层F区日间电子温度变化特征初步分析

利用曲靖非相干散射雷达2017—2018年春夏季观测数据首次分析了电离层日间150～450km电子温度的地方时与高度变化特征及其与电子密度的相关性.发现hmF2及以上的电子温度在日出日落时具有两个峰值,在11∶00—16∶00LT之间变化较小,高度越高午后上升的时间越早;从150km开始迅速增加,在约220km达到最大值,然后开始降低,在约300～350km达到最小值,最后单调上升;200km以下电子温度与电子密度成正相关(主要由热传导控制),200～450km之间存在明显的反相关(光电离过程占主导),电子-离子温度差与电子密度对数之间存在近似线性关系,电子温度逐日变化与光电离因子的变化趋势相似,这种相关性在中午与午后更明显;以上结果与其他非相干散射雷达观测和电离层模型计算结果基本一致,但也存在一些差别,需要结合更多数据深入分析.