Simulations of the seasonal variations of the thermosphere and ionosphere using a coupled,three-dimensional,global model,including variations of the interplanetary magnetic field

The University College London Thermospheric Model and the Sheffield University Ionospheric Convection Model have been integrated and improved to produce a self-consistent coupled global thermospheric/high latitude ionospheric model. The neutral thermospheric equations for wind velocity, composition, density and energy are solved, including their full interactions with the evolution of high latitude ion drift and plasma density, as these respond to convection, precipitation, solar photoionisation and changes of the thermosphere, particularly composition and wind velocity. Four 24 h Universal Time (UT) simulations have been performed. These correspond to positive and negative values of the IMF BY component at high solar activity, for a level of moderate geomagnetic activity, for each of the June and December solstices. In this paper we will describe the seasonal and IMF reponses of the coupled ionosphere/thermosphere system, as depicted by these simulations. In the winter polar region the diurnal migration of the polar convection pattern into and out of sunlight, together with ion transport, plays a major role in the plasma density structure at F-region altitudes. In the summer polar region an increase in the proportion of molecular to atomic species, created by the global seasonal thermospheric circulation and augmented by the geomagnetic forcing, controls the plasma densities at all Universal Times. The increased destruction of F-region ions in the summer polar region reduces the mean level of ionization to similar mean levels seen in winter, despite the increased level of solar insolation. In the upper thermosphere in winter for BY negative, a tongue of plasma is transported anti-sunward over the dusk side of the polar cap. To effect this transport, co-rotation and plasma convection work in the same sense. For IMF BY positive, plasma convection and co-rotation tend to oppose so that, despite similar cross-polar cap electric fields, a smaller polar cap plasma tongue is produced, distributed more centrally across the polar cap. In the summer polar cap, the enhanced plasma destruction due to enhancement of neutral molecular species and thus a changed ionospheric composition, causes F-region plasma minima at the same locations where the polar cap plasma maxima are produced in winter.     

Studies of ionospheric F-region irregularities from geomagnetic mid-latitude conjugate regions

Scintillation data from near Boston, U.S.A., and spread-F data from Argentine Islands, Antarctica are used to investigate the diurnal and seasonal variations of the simultaneous occurrence of medium-scale (~ 1–10 km) irregularities in the electron concentration in the F-region of the ionosphere at conjugate magnetic mid-latitude regions. It is found that these two stations near 52° CGL observe similar irregularity occurrence on ~75% of occasions at night when the data are considered on an hour by hour basis. During solstices, the relationship is dominated by occasions when irregularities are absent from both ends of the geomagnetic field lines; however, at equinoxes, periods of the simultaneous occurrence and non-occurrence of irregularities are approximately equally frequent. During periods of high geomagnetic activity, processes associated with the convection electric field and particle precipitation are likely to be important for the formation and transport of irregularities over these higher mid-latitude observatories. These processes are likely to occur simultaneously in conjugate regions. On days following geomagnetic activity, two processes may be operating that enhance the probability of the temperature-gradient instability, and hence lead to the formation of irregularities. These are the presence of stable auroral red arcs which occur simultaneously in conjugate locations, and the negative F-region storm effects whereby latitudinal plasma concentration gradients are increased; these effects are only similar in conjugate regions. During very quiet geomagnetic periods, F-region irregularities are occasionally observed, but seldom simultaneously at the two ends of the field lines. There is also an anomalous peak in the occurrence of irregularities over Argentine Islands associated with local sunrise in winter. No explanation is offered for these observations. Photo-electrons from the conjugate hemisphere appear to have no effect on irregularity occurrence.     

Vertical plasma drifts in the post-sunset F-region at the magnetic equator

HF doppler observations of vertical plasma drifts in the post-sunset equatorial F-region at Trivandrum (dip 0.9°S), conducted over a range of solar and geomagnetic conditions, are presented. The observations show that under magnetically quiet conditions, the characteristic post-sunset enhancement in the vertical plasma drift is quite sensitive to solar activity; the peak velocity drops by about a factor of 3 as the solar flux index (S_10.7) changes from about 125 to 70. It is found that the drift velocity enhancement has strong magnetic activity dependence only during high solar activity; the drift velocity drops by more than a factor of 2 from quiet to moderate activity, but builds back to the quiet day level for high magnetic activity. The occurrence of equatorial spread-F (ESF) is seen to be closely linked to the post-sunset enhancement in the vertical drift velocity, both showing essentially the same dependence on solar and magnetic activities. A comparison with Jicamarca observations shows that while the gross characteristics of the drift velocity pattern are about the same for the two stations, there are significant differences in the detailed variations, particularly for magnetically disturbed conditions.     

Antarctic polar cap total electron content observations from Casey Station

In early 1990 a modified JMR-1 satellite receiver system was installed at Casey Station, Antarctica (g.g. 66.28°S, 110.54° E, -80.4°A, magnetic midnight 1816UT, L = 37.8), in order to monitor the differential phase between the 150 and 400 MHz signals from polar orbiting NNSS satellites. Total electron content (TEC) was calculated using the differential phase and Casey ionosonde foF2 data, and is presented here for near sunspot maximum in August 1990 and exactly one year later. The data are used to investigate long-lived ionization enhancements at invariant latitudes polewards of − 80° A, and the ‘polar hole’, a region from −70 to − 80° A on the nightside of the polar cap where reduced electron densitiy exists because of the long transport time of plasma from the dayside across the polar cap. A comparison is made between the Casey TEC data and the Utah State University Time Dependent Ionospheric Model (TDIM) which uses as variables the solar index (F 10.7), season (summer, winter or equinox), global magnetic index (K_p), IMF B_y direction, and universal time (UT) [sojkaet al. (1991) Adv. Space Res.11(10), 39].     

Climatologies of semi-diurnal and diurnal tides in the middle atmosphere (70–110 km) at middle latitudes (40–55°)

The radars utilized are meteor (2), medium-frequency (2) and the new low-frequency (1) systems: analysis techniques have been exhaustively studied internally and comparatively and are not thought to affect the results. Emphasis is placed upon the new height-time contours of 24, 12 h tidal amplitudes and phases which best display height and seasonal structures; where possible high resolution (10 d) is used (Saskatoon) but all stations provide monthly mean resolution. At these latitudes the semi-diurnal tide is generally larger than the diurnal (10–30 m s⁻¹ vs. < 10 ms⁻¹), and displays less month to month variability. The semi-diurnal tide does show significant regular seasonal structure; wavelengths are generally small (⩽50 km) in winter, large in summer (≲ 100 km), and these states are separated by rapid equinoctial transitions. There is some evidence for less regularity toward 40°C. Coupling with mean winds is apparent. The diurnal tide has weaker seasonal variations; however there is a tendency for vertical wavelengths and amplitudes to be larger during summer months. On occasions in winter and fall wavelengths may be less than 50 km. Again the seasonal transitions are in phase with reversals of the zonal wind. Agreement with new numerical models is to be shown encouraging.     

Low latitude thermospheric meridional winds between 250 and 450 km altitude: AE-E satellite data

The daily variations of the meridional wind at ±18° latitude have been obtained for summer and winter between 1977 and 1979 using the in situ measurements from the Atmosphere Explorer-E (AE-E) satellite. The AE-E altitude increased from about 250 to about 450 km during this period, with solar activity increasing simultaneously. Data are presented at three altitudes, around 270, 350 and 440 km. It was possible to average the data to obtain the 24 h variations of the meridional wind simultaneously at northern and southern latitudes and thereby study the seasonal variation of the meridional wind in the altitude range covered. Two features are found showing significant seasonal variation: (a) a late afternoon maximum of the poleward wind occurring only in winter at 1800 LT at all three altitudes; (b) a night-time maximum in the equatorward wind—the summer equatorward wind abating earlier (near 2130 LT) and more rapidly than the winter wind (after 2300 LT). Furthermore, in summer the night-time wind reaches higher amplitudes than in winter. The night-time feature is consistent with the observed seasonal variation of the equatorial midnight temperature maximum, which occurs at or before midnight in summer and after midnight in winter, showing a stronger maximum in summer. The observed night-time abatement and seasonal variations in the night-time winds are in harmony with ground based observations at 18° latitude (Arecibo). The time difference found between summer and winter abatements of the night-time equatorward wind are in large part due to a difference between the phases of the summer and winter diurnal (fundamental) components, and diurnal amplitudes are larger in summer than in winter at all threee altitudes. However, the higher harmonics play an important role, their amplitudes being roughly 50% of the diurnal and in some instances larger. The 24 h variation is mainly diurnal at all altitudes in both summer and winter, except in winter around 2700 km altitude where the semi- and ter-diurnal components are approximately equal to or larger than the diurnal.     

foF2 hysteresis variations and the semi-annual geomagnetic wave

Seasonal and solar cycle variations of the foF2 hysteresis magnitude are investigated. Data for the noon foF2 monthly medians for Slough (51.48°N, 0.57°W), the monthly means for the sunspot numbers, and for the geomagnetic activity index aa(N) for the northern hemisphere for the period 1933–1986, covering solar cycle from 17 to 21, are used. It is found that: (1) the greatest negative amplitudes of the foF2 hysteresis variation are near the equinoxes, and (2) the solar cycle average noon foF2 hysteresis magnitude is linearly correlated with the solar cycle average semi-annual geomagnetic amplitude of the aa-index. These results support the hypothesis that the foF2 hysteresis is due to the geomagnetic activity variation during the sunspot cycle.     

Seasonal variations of day-time ionisation flows inferred from a comparison of calculated and observed NmF2

This paper reports on a comparison of calculated and observed monthly mean day-time ionospheric F2-peak density (NmF2) at a chain of stations from Japan to Australia for both solar minimum (1976) and solar maximum (1980). Nm values are calculated using the MSIS model for the observed peak heights (hmF2) and a simplified version of the continuity equation for day-time equilibrium conditions. The observed NmF2 values are always higher than the calculated ones in winter. This implies that a substantial downward flow of ionisation from above into the winter ionosphere is induced by the strongly poleward winter neutral wind which drives the ionisation down the field lines, lowering the peak height hmF2. In summer, winds are smaller, and the fluxes are more upward in comparison to winter. The seasonal variation of the ionisation fluxes and neutral winds are estimated for solar minimum, and compared with results of detailed calculations.     

Modelling studies of the conjugate-hemisphere differences in ionospheric ionization at equatorial anomaly latitudes

The relative importance of the equatorial plasma fountain (caused by vertical E x B drift at the equator) and neutral winds in leading to the ionospheric variations at equatorial-anomaly latitudes, with particular emphasis on conjugate-hemisphere differences, is investigated using a plasmasphere model. Values of ionospherec electron content (IEC) and peak electron density (Nmax) computed at conjugate points in the magnetic latitude range 10–30° at longitude 158°W reproduce the observed seasonal, solar activity, and latitudinal variations of IEC and Nmax, including the conjugate-hemisphere differences. The model results show that the plasma fountain, in the absence of neutral winds, produces almost identical effects at conjugate points in all seasons; neutral winds cause conjugate-hemisphere differences by modulating the fountain and moving the ionospheres at the conjugate hemispheres to different altitudes.At equinox., the neutral winds, mainly the zonal wind, modulate the fountain to supply more ionization to the northern hemisphere during evening and night-time hours and, at the same time, cause smaller chemical loss in the southern hemisphere by raising the ionosphere. The gain of ionization through the reduction in chemical loss is greater than that supplied by the fountain and causes stronger premidnight enhancements. in IEC and Nmax (with delayed peaks) in the southern hemisphere at all latitudes (10–30°). The same mechanism, but with the hemispheres of more flux and less chemical loss interchanged, causes stronger daytime IEC in the northern hemisphere at all latitudes. At solstice, the neutral winds, mainly the meridional wind, modulate the fountain differently at different altitudes and latitudes with a general interhemispheric flow from the summer to the winter hemisphere at altitudes above the F-region peaks. The interhemispheric flow causes stronger premidnight enhancements in IEC and Nmax and stronger daytime Nmax in the winter hemisphere, especially at latitudes equatorward of the anomaly crest. The altitude and latitude distributions of the daytime plasma flows combined with the longer daytime period can cause stronger daytime IEC in the summer hemisphere at all latitudes.     

Improved world-wide maps of monthly median foF2

New global maps of monthly median values of foF2 have been prepared using over 45,000 station months of foF2 observations, semi-empirical model values of foF2 in the mid-latitude ocean areas and empirical model values for the equatorial and high latitude regions. These observations have been carefully screened and mapped, using the Jones-Galley technique, to produce monthly median maps of foF2 for each hour, each month and for high and low levels of solar activity.     

Faraday polarization fluctuations associated with amplitude scintillations at Lunping

The diurnal, seasonal and solar cycle variations of Faraday polarization fluctuations (FPF) associated with amplitude scintillations observed at Lunping, Taiwan (25.0°N, 121.2°E : geographic) during the period 1978–1981 are presented. The occurrence of polarization fluctuations is maximum in the premidnight hours. FPFs occur either simultaneously or with a time lag after the onset of amplitude scintillations. There is an increase in FPF activity with an increase in sunspot activity. Occurrence of FPF peaks in the equinoxes. There had been a moderate activity in summer while the winter occurrence is a minimum. The seasonal occurrence pattern compared with reports from other locations indicates a longitudinal control on FPF activity. The maximum probable duration of FPF ranges from 15 to 30 min. It was found that the association of FPF with range spread-F is much better than that with frequency spread-F. Large ambient ionization densities corresponding to plasma frequencies greater than 15 MHz appear to create a favourable environment for the occurrence of FPF.     

Low latitude electrodynamic plasma drifts : a review

Radar and radio measurements have provided detailed information on the dependence of F-region electrodynamic drifts on height, season, solar cycle and magnetic activity. Recently, satellite ion drift and electric field probes have determined the variation of low latitude ionospheric drifts over a large range of altitudes and latitudes. The general characteristics of the quiet time plasma can be explained as resulting from E- and F-region dynamo and interhemispheric coupling processes. The low latitude and equatorial zonal and upward/poleward components of the plasma drift respond differently to geomagnetic activity. Disturbance dynamo effects are responsible for the drift perturbations following periods of enhanced magnetic activity. The prompt penetration of high latitude electric fields to lower latitudes produces large perturbations on the upward/poleward drifts, but has no significant effect on the low latitude and equatorial zonal drifts. A number of processes such as ‘overshielding’, ‘fossil wind’ and magnetic reconfiguration were suggested as being responsible for the direct penetration of high latitude electric fields to lower latitudes. Detailed low latitude and global numerical models were used to study the characteristics of low latitude and equatorial plasma drifts and their response to changes in the polar cap potential drop or in the high latitude field-aligned currents. These models can reproduce the latitudinal variation of the perturbation electric fields and their diurnal variations, but are still unable to account for several aspects of the experimental data as a result of the complexity of the high latitude and magnetospheric processes involved.     

The seasonal variations in the nighttime 6300 Å emission at midlatitudes from ISIS-II spacecraft

Optical limb observations at F-region heights from the ISIS-II satellite have been used to study the seasonal variations in the 6300 Å limb emission for nighttime conditions and the aeronomic implications. The observations were carried out over the American zone at northern midlatitudes, and refer mainly to the period 1973–1975 of low solar activity.The observed seasonal variations in the emission seem to be mainly controlled by the electron density at F-region heights for nighttime and quiet geomagnetic conditions. The winter minimum is found to be deeper than the summer minimum. The obervations give clear evidence of semiannual variation in the emission. The phase variations agree closely with that of the semiannual variations in electron density and neutral atmospheric density at F-region heights. However, the amplitude variations of the semiannual variations are found to be larger than suggested by the observed F-region electron density. The observations during highly disturbed conditions possibly show the presence of gravity waves with wavelengths around 500 km, which could transport auroral energy to lower latitudes. The midlatitude enhancements observed during disturbed conditions seem to be related to the inward movement of the plasmapause.     

Geomagnetic quiet daily variations in the Australian region information from a new station at Charters Towers (20.1°s)

The establishment of the new magnetic observatory at Charters Towers is described. Hourly values of the magnetic elements are analysed from two data sets (1984–1985 and 1986–1987) to provide solar diurnal variations during different months. These are compared with other stations in Australia and Papua New Guinea to study the behaviour of the focus of the Sq current systems. The unar geomagnetic variations are consistent with other Australian stations. An oceanic lunar tide is detected in the vertical element Z.     

Asymmetries in mesospheric tidal structure

Radar wind measurements made at Adelaide (35°S, 138°E) and Kyoto (35°N, 136°E) are used to construct climatologies of solar tidal wind motions in the 80–185 km region. The climatologies, in the form of contour plots of amplitude and phase of the diurnal (24 h) and semidiurnal (12 h) tides, show that there are significant asymmetries between Adelaide and Kyoto. The amplitude of the diurnal tide is significantly larger at Adelaide than at Kyoto. At both stations the phase changes in a systematic way with lime such that the phases of the zonal wind components tend to be in anti-phase at the solstices. At Adelaide, there is more evidence of the propagating (1,1) diurnal mode. At both stations, the semidiurnal tide is strongest and has the longest vertical wavelengths (>100 km) in late summer; short vertical wavelength (~ 50–80 km) oscillations are most in evidence in winter. In order to place the Adelaide and Kyoto observations in context they are compared with observations made at other latitudes and with recent numerical simulations. There is encouraging agreement between the observations and models, especially for the semidiurnal tide.     

Electron density near the plasmapause measured over one year by GEOS-2: a statistical analysis

This paper presents the results of a statistical analysis made under the 1 h average regime from the electron density data obtained over one year by the relaxation sounder on board the satellite GEOS-2 of ESA. The electron density diurnal variations, monthly and annually averaged, are sorted out. A comparison between monthly averaged, daily electron density profiles obtained over all the year has revealed the existence of a seasonal variation of plasma density in the equatorial region, the electron density being larger, on average, in summer than in winter. This seasonal variation is superimposed on the variation related to geomagnetic activity. The asymmetry of the Earth's internal magnetic field is mentioned as being a potential candidate for explaining this seasonal modulation. The annually averaged, daily profile is given. It is found to reach a maximum at 1600 LT, which indicates a plasmaspheric bulge in the predusk, rather than dusk, sector. This is found to be significantly earlier than the average LT location of the stagnation point (≳ 1800 LT) inferred from previous empirical studies for comparable geomagnetic activity levels. This latter feature is interpreted as resulting from the asymmetry of the individual daily density profiles with respect to their maximum, which has been previously reported, and whose effect, after averaging individual daily profiles having their maxima distributed around dusk, is to shift the maximum of the average profile toward the predusk sector.     

Structure in the seasonal variations of Es at southern latitudes

The diurnal variations of the seasonal characteristics of sporadic-E occurrence have been studied by analyzing a large data set of ionosonde parameters for two southern hemisphere stations. The seasonal patterns are found to display anomalous short-term variations apparently not associated with solar control or the effects of dynamic meteorology.     

The trouble with thermospheric vertical winds: geomagnetic,seasonal and solar cycle dependence at high latitudes

The vertical wind component is frequently used to determine the zero-velocity baseline for measurements of thermospheric winds by Fabry-Perot and other interferometers. For many of the upper atmospheric emission lines from which Doppler shifts are determined, for example for the OI 630 nm emission, available laboratory sources are not convenient for long-term use at remote automatic observatories. Therefore, the assumption that the long-term average vertical wind is zero is frequently used to create a baseline from which the Doppler shifts corresponding with the line-of-sight wind from other observing directions can then be calculated. A data base consisting of 1242 nights of thermospheric wind measurements from Kiruna (68°N, 20°E), a high-latitude site, has been analysed. There are many interesting short-term fluctuations of the vertical wind which will be discussed in future papers. However, the mean vertical wind at Kiruna also has a systematic variation dependent on geomagnetic activity, season and solar cycle. This means that the assumption that the average value of the vertical wind is zero over the observing period cannot be used in isolation to determine the instrument reference or baseline. Despite this note of caution, even within the auroral oval, the assumption of a zero mean vertical wind can be used to derive a baseline which is probably valid within 5 ms⁻¹ during periods of quiet geomagnetic activity (K_p < 2), near winter solstice. During other seasons, and during periods of elevated geomagnetic activity, a systematic error in excess of 10 ms⁻¹ may occur.     

On the IRI model predictions for the low-latitude ionosphere

Measurements of ionospheric electron density vertical profiles, carried out at a magnetic equatorial station located at Fortaleza (4°S, 38°W; dip latitude 2°S) in Brazil, are analyzed and compared with low-latitude electron density profiles predicted by the International Reference Ionosphere (IRI) model. The analysis performed here covers periods of high (1979/1980) and low (1986) solar activities, considering data obtained under magnetically quiet conditions representative of the summer, winter and equinox seasons. Some discrepancies are found to exist between the observed and the IRI model-predicted ionospheric electron density profiles. For high solar activity conditions the most remarkable one is the observed fast upward motion of the F-layer just after sunset, not considered in the IRI model and which precedes the occurrence of nighttime ionospheric plasma irregularities. These discrepancies are attributed mainly to dynamical effects associated with the low latitude E × B electromagnetic plasma drifts and the thermospheric neutral winds, which are not satisfactorily reproduced either in the CCIR numerical maps or in the IRI profile shapes. In particular, the pre-reversal enhancement in the vertical E × B plasma drifts around sunset hours has a great influence on the nighttime spatial distribution of the low-latitude ionospheric plasma. Also, the dynamical control exerted by the electromagnetic plasma drifts and by the thermospheric neutral winds on the low-latitude ionospheric plasma is strongly dependent on the magnetic declination angle at a given longitude. These important longitudinal and latitudinal dependences must be considered for improvement of IRI model predictions at low latitudes.     

A comparison of northern and southern hemisphere TEC storm behaviour

Changes in total electron content during magnetic storms are compared at stations with similar geographic and geomagnetic latitudes and eastward declinations in the northern and southern hemispheres.Mean patterns are obtained from 58 storms at ±35° and 28 storms at ± 20° latitude. The positive storm phase is generally larger (and earlier) in the southern hemisphere, while negative storm effects are larger in the north. These changes reduce the normal asymmetry in TEC between the two hemispheres. Composition changes calculated from the MSIS86 atmospheric model agree well with the maximum decreases in TEC in both seasons (when changes in the F-layer height are ignored). Recovery occurs with a time constant of about 35 h; this is 50% longer than in the MSIS86 model. There is a marked diurnal variation at 35°S, with a rapid overnight decay and enhanced values of TEC in the afternoon. This pattern is inverted (and weaker) at 35°N, where night-time decay is consistently slower than on undisturbed nights. These results require a diurnal change in composition of opposite sign in the two hemispheres, or enhanced westward winds at night changing to eastward near sunrise. There is some evidence for both these mechanisms. Following a night-time sudden commencement there is a large annual effect with daytime TEC increasing for storms near the June solstice and decreasing near December. Storms occurring between November and April tend to give large, irregular increases in TEC for several days, particularly at low latitudes. In summer and winter at both stations, the mean size of the negative phase does not increase for storms with K_p> 6. The size of the positive phase is proportional to the size of the change in a_p in winter, while in summer a positive phase is seen only for the larger storms.