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.     

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.     

Modelling F2-layer seasonal trends and day-to-day variability driven by coupling with the lower atmosphere

This paper presents results from the TIME-GCM-CCM3 thermosphere–ionosphere–lower atmosphere flux-coupled model, and investigates how well the model simulates known F2-layer day/night and seasonal behaviour and patterns of day-to-day variability at seven ionosonde stations. Of the many possible contributors to F2-layer variability, the present work includes only the influence of ‘meteorological’ disturbances transmitted from lower levels in the atmosphere, solar and geomagnetic conditions being held at constant levels throughout a model year.In comparison to ionosonde data, TIME-GCM-CCM3 models the peak electron density (NmF2) quite well, except for overemphasizing the daytime summer/winter anomaly in both hemispheres and seriously underestimating night NmF2 in summer. The peak height hmF2 is satisfactorily modelled by day, except that the model does not reproduce its observed semiannual variation. Nighttime values of hmF2 are much too low, thus causing low model values of night NmF2. Comparison of the variations of NmF2 and the neutral [O/N₂] ratio supports the idea that both annual and semiannual variations of F2-layer electron density are largely caused by changes of neutral composition, which in turn are driven by the global thermospheric circulation.Finally, the paper describes and discusses the characteristics of the F2-layer response to the imposed ‘meteorological’ disturbances. The ionospheric response is evaluated as the standard deviations of five ionospheric parameters for each station within 11-day blocks of data. At any one station, the patterns of variability show some coherence between different parameters, such as peak electron density and the neutral atomic/molecular ratio. Coherence between stations is found only between the closest pairs, some 2500 km apart, which is presumably related to the scale size of the ‘meteorological’ disturbances. The F2-layer day-to-day variability appears to be related more to variations in winds than to variations of thermospheric composition.     

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.     

The features and a possible mechanism of semiannual variation in the peak electron density of the low latitude F2 layer

Ionospheric data observed in 30 stations located in 3 longitude sectors (East Asia/Australia Sector, Europe/Africa Sector and America/East Pacific Ocean Sector) during 1974–1986 are used to analyse the characteristics of semiannual variation in the peak electron density of F2 layer (NmF2). The results indicate that the semiannual variation of NmF2 mainly presents in daytime. In nighttime, except in the region of geomagnetic equator between the two crests of ionospheric equatorial anomaly, NmF2 has no obvious semiannual variation. In the high latitude region, only in solar maxima years and in daytime, there are obvious semiannual variations of NmF2. The amplitude distribution of the semiannual variation of daytime NmF2 with latitude has a “double-humped structure”, which is very similar to the ionospheric equatorial anomaly. There is asymmetry between the Southern and the Northern Hemispheres of the profile of the amplitude of semiannual variation of NmF2 and longitudinal difference. A new possible mechanism of semiannual variation of NmF2 is put forward in this paper. The semiannual variation of the diurnal tide in the lower thermosphere induces the semiannual variation of the amplitude of the equatorial electrojet. This causes the semiannual variation of the amplitude of ionospheric equatorial anomaly through fountain effect. This process induces the semiannual variation of the low latitude NmF2.     

Daytime F2-layer positive storm effect at middle and lower latitudes

Daytime F2-layer positive storm effects at middle and lower latitudes in the winter thermosphere are analyzed using AE-C, ESRO-4 neutral gas composition data, ground-based ionosonde observations and model calculations. Different longitudinal sectors marked by the storm onset as ‘night-time’ and ‘daytime’ demonstrate different F2-layer positive storm mechanisms. Neutral composition changes in the ‘night-time’ sector with increased [O] and [N₂] absolute concentrations, while (N₂/O)_storm/(N₂/O)_quiet\approx1 at F2-layer heights, are shown to contribute largely to the background N_mF2 increase at lower latitudes lasting during daytime hours. Storm-induced surges of the equatorward wind give rise to an additional N_mF2 increase above this background level. The mid-latitude F2-layer positive storm effect in the ‘daytime’ sector is due to the vertical plasma drift increase, resulting from the interaction of background (poleward) and storm-induced (equatorward) thermospheric winds, but not to changes of [O] and [N₂] concentrations.     

Regional morphology features of electron concentration disturbances at the midlatitude <Emphasis Type="Italic">F</Emphasis>2-layer maximum during magnetic superstorm of July 15, 2000

Using data from ground-based ionospheric sounding stations, we studied the morphologic features of the disturbance pattern of the electron concentration at the midlatitude F2-layer maximum (NmF2) in the period of a magnetic superstorm, which began on July 15, 2000. In the Southern (winter) Hemisphere in the latitudinal sector, where the main storm phase began after sunrise, negative NmF disturbances were observed at quite high midlatitudes both day and night; whereas large positive NmF disturbances took place at lower midlatitudes in nighttime hours. In the Northern (summer) Hemisphere at latitudes where the main storm phase occurred in the local evening, only long-term negative disturbances were observed in daytime and nighttime hours; whereas at latitudes where the main storm phase began in the afternoon, NmF2 experienced both negative and positive disturbances. Based on analysis of data of KOMPSAT-l, ROCSAT-1, DMSP F13, F14, and F15 satellites, we present clear arguments for the viewpoint of many authors that it is just the enhancement of the eastward electric field in the evening sector that led to formation of the large-scale trough in the nighttime low-latitude upper ionosphere. This field enhancement was due to penetration of the magnetospheric electric field to low latitudes, not to the dynamo action of the disturbed neutral wind. It is also shown that, due to equatorward expansion of the magnetospheric convection system during the main storm phase, the plasmapause and the main ionospheric trough were shifted to a magnetic latitude of 40° (L ∼ 1.7).     

Spectral energy contributions of quasi-periodic oscillations (2–35 days) to the variability of the f oF2

The relative contributions of quasi-periodic oscillations from 2 to 35 days to the variability of f_oF2 at middle northern latitudes between 42°N and 60°N are investigated. The f_oF2 hourly data for the whole solar cycle 21 (1976–1986) for four European ionospheric stations Rome (41.9°N, 12.5°E), Poitiers (46.5°N, 0.3°E), Kaliningrad (54.7°N, 20.6°E) and Uppsala (59.8°N, 17.6°E) are used for analysis. The relative contributions of different periodic bands due to planetary wave activity and solar flux variations are evaluated by integrated percent contributions of spectral energy for these bands. The observations suggest that a clearly expressed seasonal variation of percent contributions exists with maximum at summer solstice and minimum at winter solstice for all periodic bands. The contributions for summer increase when the latitude increases. The contributions are modulated by the solar cycle and simultaneously influenced by the long-term geomagnetic activity variations. The greater percentage of spectral energy between 2 to 35 days is contributed by the periodic bands related to the middle atmosphere planetary wave activity.     

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.     

Seasonal changes in the amplitude and phase of diurnal and semidiurnal variations in parameters of the midlatitude F2 layer under minimum solar activity conditions according to the data of an ionospheric station in Irkutsk (52.5°N, 104.0°E)

We performed a statistical and spectral analysis of variations in two main parameters of the ionospheric F2 layer: critical frequency (f ₀F2) and peak height (h _m F2), recorded at an ionospheric station in Irkutsk (52.5°N, 104.0°E) in the period from December 1, 2006, to January 31, 2008, under low solar activity conditions. It was found that the f ₀F2 and h _m F2 variations contained quasi-harmonic oscillations with periods T _n = 24/n h (n = 1−7). We studied the seasonal changes in the mean and median values of monthly f ₀F2 and h _m F2 time series, their spectra, as well as the amplitudes and phases of the diurnal (n = 1) and semidiurnal (n = 2) variations. It is shown that the amplitude of the diurnal f ₀F2 variations was maximal in October–March 2007 and minimal in May–August 2007. The diurnal f ₀F2 variations were maximal at noon in the winter months and at 1600 LT in the summer months. The semidiurnal f ₀F2 variations had two maxima: a primary maximum in December and January and a secondary maximum in May–July. The maxima of semidiurnal f ₀F2 variations were shifted from 0000 and 1200 LT in winter to 0900 and 2100 LT in summer.     

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.     

Ionospheric slab thickness and its seasonal variations observed by GPS

The ionospheric slab thickness, the ratio of the total electron content (TEC) to the F2-layer peak electron density (NmF2), is closely related to the shape of the ionospheric electron density profile Ne (h) and the TEC. Therefore, the ionospheric slab thickness is a significant parameter representative of the ionosphere. In this paper, the continuous GPS observations in South Korea are firstly used to study the equivalent slab thickness (EST) and its seasonal variability. The averaged diurnal medians of December–January–February (DJF), March–April–May (MAM), June–July–August (JJA) and September–October–November (SON) in 2003 have been considered to represent the winter, spring, summer and autumn seasons, respectively. The results show that the systematic diurnal changes of TEC, NmF2 and EST significantly appeared in each season and the higher values of TEC and NmF2 are observed during the equinoxes (semiannual anomaly) as well as in the mid-daytime of each season. The EST is significantly smaller in winter than in summer, but with a consistent variation pattern. During 14–16 LT in daytime, the larger EST values are observed in spring and autumn, while the smaller ones are in summer and winter. The peaks of EST diurnal variation are around 10–18 LT which are probably caused by the action of the thermospheric wind and the plasmapheric flow into the F2-region.     

Letter to the editor

Using a method suggested by the authors earlier, the long-term trends of the F2-layer critical frequency, foF2 are derived for a set of ionospheric stations with a wide latitudinal and longitudinal coverage. All the trends are found to be negative. A pronounced dependence on geomagnetic latitude is found, the trend magnitude increasing with the latter. No globe scale longitudinal effect in trends is detected. For the majority of the stations there is also a pronounced seasonal effect, the trend magnitude being higher in summer than in winter.     

The F-layer at sunrise

It was noticed 60 years ago that at sunrise (i) the ionospheric critical frequency foF2 increases faster in winter than in summer and (ii) the increase begins at a greater solar zenith angle in winter. It was later suggested that this “seasonal sunrise anomaly” is due to a seasonal change of atomic/molecular ratio in the neutral air in the F2-layer. This paper uses the Sheffield University plasmasphere-ionosphere computational model (SUPIM) with the MSIS thermosphere model to examine the relative importance at sunrise of production, loss and diffusion processes, and the effect of neutral air winds. The results show that both (i) and (ii) can be explained in terms of neutral composition changes.     

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.     

Studie über Höhenschwankungen der F2-Schicht: I

Zusammenfassung Die mittags gemessenen scheinbaren Höhen (hF2) der F2-Schicht unterliegen jährlichen Schwankungen, die örtlich verschieden sind. Die Jahresmittelwerte und die Amplitude der ersten Harmonischen (einfache jährliche Periode) zeigen eine deutliche Abhängigkeit von der geographischen Breite des Beobachtungsortes. Während die Phase der einfachen Periode Höhenmaxima zum höchsten Sonnenstand ergibt, ist die Phase der doppelten Periode um rund 180° verschoben, so dass in äquatorialen Gebieten die Höhenmaxima zu den Solstitien auftreten.

---

The virtual height of the ionospheric F2-layer varies with the location of the ionospheric observatory and with the season. The annual mean of these heights and likewise the annual component of their variation varies with the geographical latitude. The greatest variation appears in high latitude, on the other hand the annual variation disappears on the equator. There remains only a half-yearly component. The maxima of F₂-heights appear in high latitude at summer-solstice and on the equator at both solstices.

     

Seasonal variation of the total electron content, maximum electron density and equivalent slab thickness at a South-American station

Total electron content (TEC) and foF2 ionosonde data obtained at Tucumán (26.9°S; 65.4°W) from April 1982 to March 1983 (high solar activity period) are analyzed to show the seasonal variation of TEC, NmF2 (proportional to square of foF2) and the equivalent slab thickness EST. Bimonthly averages of the monthly median for January–February, April–May, July–August and October–November have been considered to represent summer, autumn, winter and spring seasons, respectively. The results show that the higher values of TEC and maximum electron density of F2-layer NmF2 are observed during the equinoxes (semiannual anomaly). During daytime, both in TEC and in NmF2 the seasonal or winter anomaly can be seen. At nighttime, this effect is not observed. Also, the observed NmF2 values are used to check the validity of International Reference Ionosphere (IRI) to predict the seasonal variability of this parameter. In general, it is found that averaged monthly medians (obtained with the IRI model) overestimate averaged monthly median data for some hours of the day and underestimate for the other hours.     

Gravity wave intensity and momentum fluxes in the mesosphere over Shigaraki, Japan (35°N, 136°E) during 1986–1997

Averaged seasonal variations of wind perturbation intensities and vertical flux of horizontal momentum produced by internal gravity waves (IGWs) with periods 0.2/1 h and 1/6 h are studied at the altitudes 65/80 km using the MU radar measurement data from the middle and upper atmosphere during 1986/1997 at Shigaraki, Japan (35°N, 136°E). IGW intensity has maxima in winter and summer, winter values having substantial interannual variations. Mean wave momentum flux is directed to the west in winter and to the east in summer, opposite to the mean wind in the middle atmosphere. Major IGW momentum fluxes come to the mesosphere over Shigaraki from the Pacific direction in winter and continental Asia in summer.     

A study of the noont critical frequencies of theE andF1 layers of the ionosphere

Summary The mean monthly noon critical frequencies of theE andF1 layers of the ionosphere at a number of stations in different latitudes and their variation with sunspot number have been studied in this paper. It is found that while theE layer approximates to a chapman region, theF1 layer is markedly affected by other agenices, somewhat similar toF2. In high sunspot years,foF1 shows two maxima at middle latitudes with a minimum at the equator.     

Thermospheric Meridional Wind Control on Equatorial Scintillations and the Role of the Evening F-Region Height Rise, E?×?B Drift Velocities and F2-Peak Density Gradients

The possible role, on L-band scintillation activity, played by the nighttime magnetic meridional component of the thermospheric horizontal neutral winds, the post-sunset F-layer base height, the electrical field pre-reversal enhancement (PRE) and the latitudinal gradients of the F2-layer peak density is analyzed, considering different cases of scintillation occurrence (and their latitudinal extent) during August and September 2002. The meridional winds were derived over low-latitudes from a modified form of the nonlinear time-dependent servo-model. A chain of two scintillation monitors and three digital ionosondes was operational in Brazil and used to collect, respectively, global positioning system signal amplitude scintillation and ionospheric height (h′F; hpF2) and frequency (foF2) parameters. From the overall behavior in the 2 months analyzed, the results suggest that high near sunset upward vertical plasma drifts are conducive for the generation of spread-F irregularities, whereas large poleward meridional winds tend to suppress the development of plasma bubble irregularities and the occurrence of their associated scintillations. Even when generated, a reduced fountain effect, due to weak electric field PRE, acts for the bubbles to be expanded less effectively to higher latitudes. The results also reveal that high F-layer base and peak heights (at equatorial and off-equatorial latitudes), and intense gradients in the F2-peak density between the dip equator and the equatorial anomaly crests, are favorable conditions for the generation of F-region irregularities and increased scintillation activity. Other distinct features of the controlling factors in the cases of occurrence and non-occurrence of equatorial scintillations are presented and discussed.