Modelling the peak of the ionospheric E-layer

Calculations with a full time-varying model are used to study changes in the height and density of the E-layer peak, caused by known changes in the neutral atmosphere. Agreement with mean observed values of N_mE requires an increase of 10% in calculated ion densities, and an increase of 33% in the solar-maximum EUV model at λ<150 Å. At a fixed site, changes with the solar zenith angle χ agree well with the simple Chapman theory during most of the daylight hours. Simple modifications to the Chapman equations give improved accuracy near sunrise and sunset. When corrected for changes in χ, model results for summer and equinox show a decrease in the peak density N_mE at increasing latitudes. The overall change agrees well with experimental data, as summarised in the IRI model. Known changes in the neutral atmosphere also reproduce the increase in N_mE in winter, at latitudes up to 30°. The continuing increase at higher winter latitudes, in the IRI model, requires a major reduction in NO densities in winter. A suitable compromise is suggested. Equations fitted to the model results then provide a simpler and better behaved replacement for the IRI equations. Calculations at night show that known sources of ionisation, largely from starlight, can produce observed peak densities using current chemistry. There is an appreciable change with latitude, as starlight production increases in the southern hemisphere. The improbably large solar cycle change built into the IRI model, at night, cannot be reproduced and is not found in recent data. A new, simpler model is suggested. Changes in zenith angle and atmospheric composition cause the peak height (h_mE) to vary between 105 and 120 km, as a function of time, latitude, season and solar flux. These changes are approximated by simple equations that should be definitely preferable over the single, fixed height used in the IRI models.     

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.     

Model simulation of the equatorial electrojet in the Peruvian and Philippine sectors

Between 100 and 120 km height at the Earth's magnetic equator, the equatorial electrojet (EEJ) flows as an enhanced eastward current in the daytime E region ionosphere, which can induce a magnetic perturbation on the ground. Calculating the difference between the horizontal components of magnetic perturbation (H) at magnetometers near the equator and about 6–9° away from the equator, ΔH, provides us with information about the strength of the EEJ. The NCAR Thermosphere–Ionosphere–Electrodynamics General Circulation Model (TIE-GCM) is capable of simulating the EEJ current and its magnetic perturbation on the ground. The simulated diurnal, seasonal (March equinox, June solstice, December solstice), and solar activity (F_10.7=80, 140 and 200 units) variations of ΔH in the Peruvian (76°W) and Philippine (121°E) sectors, and the relation of ΔH to the ionospheric vertical drift velocity, are presented in this paper. Results show the diurnal, seasonal and solar activity variations are captured well by the model. Agreements between simulated and observed magnitudes of ΔH and its linear relationship to vertical drift are improved by modifying the standard daytime E region photoionization in the TIE-GCM in order to better simulate observed E region electron densities.     

Dynamics of the midlatitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at sunrise

The behavior of the F2 layer at sunrise has been studied based on vertical-incidence ionospheric sounding data in Almaty (76°55′E, 43°15′N). Records with small amplitudes of electron density background fluctuations were selected in order to exactly estimate the onsets of a pronounced increase in the electron density at different altitudes. It has been indicated that the electron density growth rate is a function of altitude; in this case, the growth rate at the F2 layer maximum is much lower than such values at fixed altitudes of ～30–55 km below the layer maximum. The solar zenith angle (χ) and the blanketing layer thickness (h ₀) at the beginning of a pronounced increase in the electron density at altitude h are linearly related to the h value, and these quantities vary within ～90° < χ < 100° and 180 km < h ₀ < 260 km, respectively.     

Comparison of satellite measurements of the low-latitude nighttime upper ionosphere with IRI

In this paper, we report the results of our comparison study between satellite measurements and the International Reference Ionosphere (IRI) model on the seasonal and longitudinal changes of the low-latitude nighttime topside ionosphere during the period of solar maximum from June 2000 to July 2001. Satellite measurements were made by KOMPSAT-1 and DMSP F15 at 685 km altitude and 840 km altitude, respectively. The results show that the IRI2001 model gives reasonable density estimations for the summer hemisphere and the March equinox at both altitudes. However, the observed wintertime densities are smaller than the predictions of the IRI2001 model, especially at a higher (840 km) altitude, manifesting strong hemispheric asymmetries. The observed electron temperatures generally reside between the two estimations of IRI2001, one based on the Aeros–ISIS data and the other based on Intercosmos, and the latter estimation better represents the observations. With more or less monotonic increase with latitude, the temperature profiles of the IRI2001 model do not predict the enhancement seen around 15° magnetic latitude of the winter hemisphere. Longitudinal variation, probably caused by the zonal winds, is seen in all seasons at both altitudes, while the IRI2001 model does not show a large variation. The observed density and temperature show significant changes according to the F10.7 values in the whole low-latitude region from 40°S to 40°N geomagnetic latitude. The effect is manifested as increases in the density and temperature, but not in the hemispheric asymmetry or in the longitudinal variation.     

Annual and seasonal variations in the low-latitude topside ionosphere

Annual and seasonal variations in the low-latitude topside ionosphere are investigated using observations made by the Hinotori satellite and the Sheffield University Plasmasphere Ionosphere Model (SUPIM). The observed electron densities at 600 km altitude show a strong annual anomaly at all longitudes. The average electron densities of conjugate latitudes within the latitude range ±25° are higher at the December solstice than at the June solstice by about 100% during daytime and 30% during night-time. Model calculations show that the annual variations in the neutral gas densities play important roles. The model values obtained from calculations with inputs for the neutral densities obtained from MSIS86 reproduce the general behaviour of the observed annual anomaly. However, the differences in the modelled electron densities at the two solstices are only about 30% of that seen in the observed values. The model calculations suggest that while the differences between the solstice values of neutral wind, resulting from the coupling of the neutral gas and plasma, may also make a significant contribution to the daytime annual anomaly, the E × B drift velocity may slightly weaken the annual anomaly during daytime and strengthen the anomaly during the post-sunset period. It is suggested that energy sources, other than those arising from the 6% difference in the solar EUV fluxes at the two solstices due to the change in the Sun-Earth distance, may contribute to the annual anomaly. Observations show strong seasonal variations at the solstices, with the electron density at 600 km altitude being higher in the summer hemisphere than in the winter hemisphere, contrary to the behaviour in NmF2. Model calculations confirm that the seasonal behaviour results from effects caused by transequatorial component of the neutral wind in the direction summer hemisphere to winter hemisphere.     

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.     

Results of foF2 and Ne-h profiles at low latitude using recent digital ionosonde observations and their comparison with IRI-2000

Using digital ionosonde observations at low-latitude station, Delhi (28.6 N, 77.2 E, mag. dip 42.4 N), the diurnal and seasonal variations of the critical frequency of F2 layer (foF2) are analyzed from August 2000 to July 2001 during a high solar activity period. Also, noontime bottomside electron density (Ne-h) profiles, below the F2-peak, are derived from ionogram, using the POLAN (Report UAG-93, WDC-A, for Solar Terrestrial Physics, Boulder, Co.) program during the same period, and these profiles are then normalized to the peak height and density (hmF2, NmF2) of the F2-region. These observations are used to assess the predictability of the International Reference Ionosphere, IRI-2000 model (Radio Sc. 36(2) (2001) 261). Results show in general, a large variability, (1σ, σ is standard deviation), in foF2 during nighttime than daytime during winter and equinox, the variability of foF2 about the mean is about ±25% by night and ±15% by day. The IRI model shows a fairly good agreement with foF2 observations during daytime, however during nighttime, the discrepancies between the two exist. Comparative studies of the normalized observed profiles with those obtained with the IRI model (Bilitza, 2001) using both the options namely: Gulyaeva's (Adv. Space Res. 7 (1987) 39) model and B0-Table (Adv. Space Res. 25(1) (2000) 89), show that during all the seasons, in general, the B0-Tab option, reveals a better agreement with the observations, while the IRI model using Gulyaeva's option, overestimates the electron density distribution during summer and equinox, however, during winter, the model is close to the observations. The comparisons of average profile shape parameters (B0,B1) derived from noontime observed profiles, with those obtained, using B0-Tab option, in the IRI model, show a good agreement during all the seasons. However, B0, B1 obtained, using Gulyaeva's option in the IRI model, show a disagreement with the derived B0, B1 values during all the seasons, except during winter, for B0 parameter.     

Anisotropy and temporal variations in the fast azimuth of the fast shear wave beneath southern Kamchatka during 1993–2002

Parameters of split shear waves from local earthquakes in the area of the PET IRIS station (town of Petropavlovsk-Kamchatski) were measured over the period 1993–2002 for the study of anisotropic properties of rocks in the subduction zone and variations in the fast azimuth of the fast shear wave (?). The dominating fast shear wave polarization directions were oriented in 1993–2002 along N90°E ± 20° in agreement with the direction of the Pacific plate motion. The normalized shear wave delay times δt _SS increase to a depth of 150 km. The values of δt _SS are largest (up to 20 ms/km) for earthquakes at depths of 50–60 and 90–150 km and smallest (up to 6 ms/km) for earthquakes at depths greater than 200 km. The fast azimuths for events with H < 80 km are described in terms of a horizontal transversely isotropic (HTI) model of the medium, with the axis oriented northward. Temporal variations in the fast azimuths with an amplitude of up to 90° and a predominant period of about 400–600 days are observed for events at depths of 80–120 km. The anisotropy of rocks is described by effective models of the orthorhombic and HTI symmetries. The predominant fast shear wave fast azimuths from events at depths of 120–310 km vary with time: the polarization axis was oriented to the north in 1993–1995, to the north and east in 1996–1998, to the east in 1999–2000, and to the northeast and southeast in 2001–2002. The anisotropy of rocks can be described in terms of the HTI model with the symmetry axis subparallel to the focal zone dip.     

Atmospheric temperature profiling by joint Raman,Rayleigh and Fe Boltzmann lidar measurements

Simultaneous and complete temperature profiles from near ground to about 100 km are essential for studying the dynamical coupling between different atmospheric layers. They are acquired by combining three different lidar techniques at Wuhan, China (30.5°N, 114.4°E). The atmospheric temperatures from about 3 to 25 km are calculated from the nitrogen molecule density profiles obtained from the N₂ vibrational Raman backscatter, while the atmospheric temperatures between 30 and ∼75 km are calculated by the standard Rayleigh scattering method. The temperatures in the 80–100 km altitude region are derived from the Fe Boltzmann technique. The temperature profiles measured by our lidar systems exhibit good agreement when compared with the radiosonde and satellite data, as well as the model. A Lomb–Scargle spectral analysis of the normalized temperature perturbations in the altitude range from 4 to 60 km shows that the spectral slopes of the vertical wave number spectra tended to −3 for large vertical wave numbers. This is consistent with the model predictions of saturated gravity wave spectra.     

Response of the midlatitude ionosphere to extreme geomagnetic storms of the 23rd solar cycle

The ionospheric response in the Irkutsk region (52.3° N, 104.3° E) to the extreme geomagnetic storms of solar cycle 23 was studied based on the data of the Irkutsk incoherent scatter radar (ISR) and DPS-4 vertical sounding digital ionosonde. The deviations of parameters from the undisturbed level, i.e., from the monthly medians or the values obtained on a quiet day, were considered as an ionospheric response. Values of the electron concentration maximum (N _mF2) and electron temperature (T _e) at a height of 350 km were chosen as parameters. The ionospheric response is interpreted in the scope of the concept of a thermospheric storm and penetration of the magnetospheric electric field.     

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.     

Solar activity dependence of electron and ion temperatures using SROSS C2 RPA data and comparison with IRI model

Electron and ion temperature (T_e and T_i) data observed using RPA on board SROSS C2 satellite are investigated for the variation with local time, season, latitude (0–30°N geographic) over a half of a solar cycle (1995–2000). The nighttime T_e (∼1000 K) is independent of the season and the solar flux whereas T_i exhibits positive correlation with the solar activity during all three seasons. In the early morning hours during summer, T_e is higher by ∼500 K than other seasons in all three levels of solar activity. During winter and equinox in the early morning hours, T_e and T_i are higher during low solar activity, showing a negative correlation with solar flux. During daytime, the T_i increases with the solar flux in winter and summer solstice, but is independent in equinox. IRI underestimates T_e and T_i during the morning period by 50–75% in the equatorial and near-equatorial stations during all levels of solar activities.     

Comparison between IRI predictions and digital ionosonde measurements of hmF2 at New Delhi during low and moderate solar activity

This paper deals with the diurnal and seasonal variations of height of the peak electron density of the F2-layer (hmF2) derived from digital ionosonde measurements at a low–middle-latitude station, New Delhi (28.6°N, 77.2°E, dip 42.4°N). Diurnal and seasonal variations of hmF2 are examined and comparisons of the observations are made with the predictions of the International Reference Ionosphere (IRI-2001) model. Our study shows that during both the moderate and low solar activity periods, the diurnal pattern of median hmF2 reveals a more or less similar trend during all the seasons with pre-sunrise and daytime peaks during winter and equinox except during summer, where the pre-sunrise peak is absent. Comparison of observed median hmF2 values with the IRI during moderate and low solar activity periods, in general, reveals an IRI overestimation in hmF2 during all the seasons for local times from about 06 LT till midnight hours except during summer for low solar activity, while outside this time period, the observed hmF2 values are close to the IRI predictions. The hmF2 representation in the IRI model does not reproduce pre-sunrise peaks occurring at about 05 LT during winter and equinox as seen in the observations during both the solar activity periods. The noontime observed median hmF2 values increase by about 10–25% from low (2004–2005) to high solar activity (2001–2002) during winter and equinox, while the IRI in the same time period and seasons shows an increase of about 10–20%. During summer, however, the observed noontime median hmF2 values show a little increase with the solar activity, as compared to the IRI with an increase of about 12%.     

Height profiles of the amplitudes of large-scale traveling ionospheric disturbances

Nighttime height profiles of the amplitudes of large-scale traveling ionospheric disturbances (LSTIDs) obtained from the data of vertical sounding in Almaty (76°55′ E, 43°15′ N) for the period 2000–2007 are analyzed. The height profiles are plotted using the time variations in electron density N _h (t) at a series of heights for the F region in the ionosphere with a height step of 10 km. In total, observations were conducted during 1166 nights, among which 581 nights are characterized by wave activity. Nights with the maximum amplitude of variations in N _h (t) exceeding 25% are selected for analysis. The total number of such nights is 63; LSTIDs have been recorded in both magnetically quiet and active periods. The regressive ratios between the height of the F-region maximum and the height that corresponds to the maximum absolute amplitude of a wave, as well as between the values of the maximum amplitude at a height profile and the value of the amplitude of variations in N _m F(t) at the layer maximum, are obtained.     

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.     

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.     

Trimpi perturbations from large ionisation enhancement patches

A number of increasingly sophisticated and realistic models have been developed in order to investigate the interaction between sub-ionospherically propagating VLF waves and regions of ionisation enhancement (LIE¹) in the D-region caused by lightning-induced electron precipitation enhancements (LEP). This LEP-produced LIE can result in phase and amplitude perturbations on received VLF radio signals that are referred to as Trimpis or more precisely, classic Trimpis, to distinguish them from “early/fast Trimpis” or “VLF sprites” which are not caused by LEP and are not considered here. It is important, for comparison with experimentally observed Trimpi effects, that the spatial extent of the D-region electron density (N_e) perturbation is modeled accurately. Here, it is argued that most previous modeling has used patch (LIE) sizes that are typically up to 100 km in both latitudinal and longitudinal extent, which are generally smaller than those that actually occur for real lightning induced electron precipitation events. It would also appear that maximum ΔN_e values assumed have often been too large, and the patches (LIEs) have been incorrectly modelled as circular rather than elliptical in horizontal extent. Consequently, in the present work, Trimpi perturbations are determined for LIEs with smaller maximum ΔN_e, larger spatial extent and elliptical shape. Calculations of VLF Trimpis have been made as a function of the horizontal coordinates of the LIE centre, over the whole rectangular corridor linking transmitter and receiver. The Trimpi modelling program is fully 3D, and takes account of modal mixing at the LIE. The underlying theory assumes weak Born scattering, but the code calculates a non-Born skin depth attenuation function for the LIE in question. The LIE is modelled as an electron density enhancement with a Gaussian profile in all coordinates. Results for a large elliptical LIE ∼ 200 × 600 km show that significant Trimpis, ∼−0.4 dB in amplitude and ∼+4° in phase are predicted, using modest maximum ΔN_e values ∼ 1.5 el/cc. Such an electron density enhancement is well within the range that would be expected to result from experimentally observed fluxes of electron precipitation following wave particle interactions with whistler-mode waves.     

Statistical properties of nonlinear wave signatures in OH and O2 airglow brightness data observed at lower midlatitudes

Among 2187 nights of airglow observations of the OH(6-2) and O₂b(0-1) bands from Argentina (mainly from El Leoncito, 32°S 69°W), 132 show airglow brightness jumps (ABJs) of short duration (16 min median). ABJs are supposed to be related to mesospheric bores or similar nonlinear waves. Several occurrence patterns were identified, which a successful explanation must take into account. ABJs occur preferably in the OH layer at 87 km, and are less likely in the O₂ layer at 95 km, maybe because ducts prefer lower altitudes. The seasonal distribution of nights when ABJs are observed only in the OH layer clearly shows a winter maximum centered around solstice, and equinox minima. In contrast, the seasonal distribution of ABJ nights in O₂ is flat. Most ABJs simultaneously present in OH and O₂ show anticorrelated variation between both layers. ABJ nights tend to occur in clusters lasting several days, which probably reflects duct lifetime.