The 16-day planetary wave in the mesosphere and lower thermosphere

A meteor radar located at Sheffield in the UK has been used to measure wind oscillations with periods in the range 10–28 days in the mesosphere/lower-thermosphere region at 53.5°N, 3.9°W from January 1990 to August 1994. The data reveal a motion field in which wave activity occurs over a range of frequencies and in episodes generally lasting for less than two months. A seasonal cycle is apparent in which the largest observed amplitudes are as high as 14 ms⁻¹ and are observed from January to mid-April. A minimum in activity occurs in late June to early July. A second, smaller, maximum follows in late summer/autumn where amplitudes reach up to 7–10 ms⁻¹. Considerable interannual variability is apparent but wave activity is observed in the summers of all the years examined, albeit at very small amplitudes near mid summer. This behaviour suggests that the equatorial winds in the mesopause region do not completely prevent inter-hemispheric ducting of the wave from the winter hemisphere, or that it is generated in situ.     

Transient eastward-propagating long-period waves observed over the South Pole

Observations of the horizontal wind field over the South Pole were made during 1995 using a meteor radar. These data have revealed the presence of a rich spectrum of waves over the South Pole with a distinct annual occurrence. Included in this spectrum are long-period waves, whose periods are greater than one solar day, which are propagating eastward. These waves exhibit a distinct seasonal occurrence where the envelope of wave periods decreases from a period of 10 days near the fall equinox to a minimum of 2 days near the winter solstice and then progresses towards a period near 10 days at the spring equinox. Computation of the meridional gradient of quasi-geostrophic potential vorticity has revealed a region in the high-latitude upper mesosphere which could support an instability and serve as a source for these waves. Estimation of the wave periods which would be generated from an instability in this region closely resembles the observed seasonal variation in wave periods over the South Pole. These results are consistent with the hypothesis that the observed eastward propagating long-period waves over the South Pole are generated by an instability in the polar upper mesosphere. However, given our limited data set we cannot rule out a stratospheric source. Embedded in this spectrum of eastward propagating waves during the austral winter are a number of distinct wave events. Eight such wave events have been identified and localized using a constant-Q filter bank. The periods of these wave events ranges from 1.7 to 9.8 days and all exist for at least 3 wave periods. Least squares analysis has revealed that a number of these events are inconsistent with a wave propagating zonally around the geographic pole and could be related to waves propagating around a dynamical pole which is offset from the geographic pole. Additionally, one event which was observed appears to be a standing oscillation.     

Vertical velocities associated with gravity waves measured in the mesosphere and lower thermosphere with the EISCAT VHF radar

The EISCAT VHF radar (69.4°N, 19.1°E) has been used to record vertical winds at mesopause heights on a total of 31 days between June 1990 and January 1993. The data reveal a motion field dominated by quasi-monochromatic gravity waves with representative apparent periods of 30–40 min, amplitudes of up to 2.5 m s^–1 and large vertical wavelength. In some instances waves appear to be ducted. Vertical profiles of the vertical-velocity variance display a variety of forms, with little indication of systematic wave growth with height. Daily mean variance profiles evaluated for consecutive days of recording show that the general shape of the variance profiles persists over several days. The mean variance evaluated over a 10 km height range has values from 1.2 m²s^–2 to 6.5 m²s^–2 and suggests a semi-annual seasonal cycle with equinoctial minima and solsticial maxima. The mean vertical wavenumber spectrum evaluated at heights up to 86 km has a slope (spectral index) of -1.36 ± 0.2, consistent with observations at lower heights but disagreeing with the predictions of a number of saturation theories advanced to explain gravity-wave spectra. The spectral slopes evaluated for individual days have a range of values, and steeper slopes are observed in summer than in winter. The spectra also appear to be generally steeper on days with lower mean vertical-velocity variance.     

The summertime 12-h wind oscillation with zonal wavenumber s = 1 in the lower thermosphere over the South Pole

Meteor radar measurements of winds near 95 km in four azimuth directions from the geographic South Pole are analyzed to reveal characteristics of the 12-h oscillation with zonal wavenumber one (s = 1). The wind measurements are confined to the periods from 19 January 1995 through 26 January 1996 and from 21 November 1996 through 27 January 1997. The 12-h s = 1 oscillation is found to be a predominantly summertime phenomenon, and is replaced in winter by a spectrum of oscillations with periods between 6 and 11.5 h. Both summers are characterized by minimum amplitudes (5–10 ms^–1) during early January and maxima (15–20 ms^–1) in November and late January. For 10-day means of the 12-h oscillation, smooth evolutions of phase of order 4–6 h occur during the course of the summer. In addition, there is considerable day-to-day variability (±5–10 ms^–1 in amplitude) with distinct periods (i.e., 5 days and 8 days) which suggests modulation by planetary-scale disturbances. A comparison of climatological data from Scott Base, Molodezhnaya, and Mawson stations suggests that the 12-h oscillation near 78°S is s = 1, but that at 68°S there is probably a mixture between s = 1 and other zonal wavenumber oscillations (most probably s = 2). The mechanism responsible for the existence of the 12-h s = 1 oscillation has not yet been identified. Possible origins discussed herein include in situ excitation, nonlinear interaction between the migrating semidiurnal tide and a stationary s = 1 feature, and thermal excitation in the troposphere.     

Wavelet analysis on transient behaviour of tidal amplitude fluctuations observed by meteor radar in the lower thermosphere above Bulgaria

On the basis of bispectral analysis applied to the hourly data set of neutral wind measured by meteor radar in the MLT region above Bulgaria it was demonstrated that nonlinear processes are frequently and regularly acting in the mesopause region. They contribute significantly to the short-term tidal variability and are apparently responsible for the observed complicated behavior of the tidal characteristics. A Morlet wavelet transform is proposed as a technique for studying nonstationary signals. By simulated data it was revealed that the Morlet wavelet transform is especially convenient for analyzing signals with: (1) a wide range of dominant frequencies which are localized in different time intervals; (2) amplitude and frequency modulated spectral components, and (3) singular, wave-like events, observed in the neutral wind of the MLT region and connected mainly with large-scale disturbances propagated from below. By applying a Morlet wavelet transform to the hourly values of the amplitudes of diurnal and semidiurnal tides the basic oscillations with periods of planetary waves (1.5/20 days), as well as their development in time, are obtained. A cross-wavelet analysis is used to clarify the relation between the tidal and mean neutral wind variability. The results of bispectral analysis indicate which planetary waves participated in the nonlinear coupling with the atmospheric tides, while the results of cross-wavelet analysis outline their time intervals if these interactions are local.     

Global empirical wind model for the upper mesosphere/lower thermosphere. I. Prevailing wind

An updated empirical climatic zonally averaged prevailing wind model for the upper mesosphere/lower thermosphere (70/110 km), extending from 80°N to 80°S is presented. The model is constructed from the fitting of monthly mean winds from meteor radar and MF radar measurements at more than 40 stations, well distributed over the globe. The height-latitude contour plots of monthly mean zonal and meridional winds for all months of the year, and of annual mean wind, amplitudes and phases of annual and semiannual harmonics of wind variations are analyzed to reveal the main features of the seasonal variation of the global wind structures in the Northern and Southern Hemispheres. Some results of comparison between the ground-based wind models and the space-based models are presented. It is shown that, with the exception of annual mean systematic bias between the zonal winds provided by the ground-based and space-based models, a good agreement between the models is observed. The possible origin of this bias is discussed.     

Lunar tidal winds in the upper atmosphere over Collm

The lunar semidiurnal tide in winds measured at around 90 km altitude has been isolated with amplitudes observed up to 4 m s^–1. There is a marked amplitude maximum in October and also a considerable phase variation with season. The average variation of phase with height indicated a vertical wavelength of more than 80 km but this, and other results, needs to be viewed in the light of the considerable averaging required to obtain statistical significance. Large year-to-year variations in both amplitude and phase were also found. Some phase comparisons with the GSWM model gave reasonable agreement but the model amplitudes above a height of 100 km were much larger than those measured. An attempt to make a comparison with the lunar geomagnetic tide did not yield a statistically significant result.     

A comparative study of winds and tidal variability in the mesosphere/lower-thermosphere region over Bulgaria and the UK

Meteor radars located in Bulgaria and the UK have been used to simultaneously measure winds in the mesosphere/lower-thermosphere region near 42.5°N, 26.6°E and 54.5°N, 3.9°W, respectively, over the period January 1991 to June 1992. The data have been used to investigate planetary waves and diurnal and semidiurnal tidal variability over the two sites. The tidal amplitudes at each site exhibit fluctuations as large as 300% on time scales from a few days to the intra-seasonal, with most of the variability being at intra-seasonal scales. Spectral and cross-wavelet analysis reveals closely related tidal variability over the two sites, indicating that the variability occurs on spatial scales large compared to the spacing between the two radars. In some, but not all, cases, periodic variability of tidal amplitudes is associated with simultaneously present planetary waves of similar period, suggesting the variability is a consequence of non-linear interaction. Calculation of the zonal wave number of a number of large amplitude planetary waves suggests that during summer 1991 the 2-day wave had a zonal wave number of 3, but that during January/February 1991 it had a zonal wave number of 4.     

Inertia-gravity waves in the troposphere and lower stratosphere associated with a jet stream exit region

Radar measurements at Aberystwyth (52.4°N, 4.1°W) of winds at tropospheric and lower stratospheric heights are shown for 12–13 March 1994 in a region of highly curved flow, downstream of the jet maximum. The perturbations of horizontal velocity have comparable amplitudes in the troposphere and lower stratosphere with downward and upward phase propagation, respectively, in these two height regions. The sense of rotation with increasing height in hodographs of horizontal perturbation velocity derived for hourly intervals show downwards propagation of energy in the troposphere and upward propagation in the lower stratosphere with vertical wavelengths of 1.7 to 2.3 km. The results indicate inertia-gravity waves propagating in a direction similar to that of the jet stream but at smaller velocities. Some of the features observed contrast with those of previous observations of inertia-gravity waves propagating transverse to the jet stream. The interpretation of the hodographs to derive wave parameters has taken account of the vertical shear of the background wind transverse to the direction of wave propagation.     

The nonlinear effects on the characteristics of gravity wave packets: dispersion and polarization relations

By analyzing the results of the numerical simulations of nonlinear propagation of three Gaussian gravity-wave packets in isothermal atmosphere individually, the nonlinear effects on the characteristics of gravity waves are studied quantitatively. The analyses show that during the nonlinear propagation of gravity wave packets the mean flows are accelerated and the vertical wavelengths show clear reduction due to nonlinearity. On the other hand, though nonlinear effects exist, the time variations of the frequencies of gravity wave packets are close to those derived from the dispersion relation and the amplitude and phase relations of wave-associated disturbance components are consistent with the predictions of the polarization relation of gravity waves. This indicates that the dispersion and polarization relations based on the linear gravity wave theory can be applied extensively in the nonlinear region.     

Planetary wave activity in the lower ionosphere during CRISTA I campaign in autumn 1994 (October-November)

On the basis of MEM spectrum analysis, the main planetary scale fluctuations formed in the lower ionosphere are studied over a period of 3–25 days during the CRISTA campaign (October-November 1994). Three dominant period bands are found: 3–5, 6–8 and 15–23 (mainly 16–18) days. For 7–8 and 16–18 day fluctuations, propagation was eastward with wave numbers K = 3 and K = 1, respectively. The magnitude of planetary wave activity in the mid-latitudes of the Northern Hemisphere during the CRISTA campaign seems to be fairly consistent with the expected undisturbed normal/climatological state of the atmosphere at altitudes of 80–100 km.     

Diurnal and semidiurnal tide in the upper middle atmosphere during the first year of simultaneous MF radar observations in northern and southern Japan (45°N and 31°N)

The climatology of mean wind, diurnal and semidiurnal tide during the first year (1996–1997) of simultaneous wind observations at Wakkanai (45.4°N, 141.7°E) and Yamagawa (31.2°N, 130.6°E) is presented. The locations of the radars allow us to describe the latitudinal dependence of the tides. Tidal amplitude and phase profiles are compared with those of the global scale wave model (GSWM). While the observed amplitude profiles of the diurnal tide agree well with the GSWM values, the observed phase profiles often indicate longer vertical wavelengths than the GSWM phase profiles. In contrast to the GSWM simulation, the observations show a strong bimodal structure of the diurnal tide, with the phase advancing about 6 hours from summer to winter.     

Penetration of the Earth’s free oscillations at 54 minute period into the atmosphere

It is known that the fundamental spheroidal mode ₀S₂ of the Earth free oscillation with a period of about 54 min forces atmospheric oscillations. We present a certain phase relationship for components of the ₀S₂ multiplet, which is based on synchronous collocated microbarograph and seismograph observations. This relationship is both the first observational manifestation of the Pekeris mode of global atmospheric oscillations with the 54 min period and a further proof of the Earths ₀S₂ mode penetrating into the atmosphere. We show that the linear non-dissipative model of steady forced oscillations in isothermal atmosphere at rest does not describe the penetration of the ₀S₂ mode into the atmosphere adequately.     

Trends in the characteristics of the annual and semiannual variations observed in the radio wave absorption in the lower ionosphere

The continuous increase in concentration of greenhouse gases in the atmosphere is expected to cool higher levels of the atmosphere. There is some direct and indirect experimental evidence of long-term trends in temperature and other parameters in the mesosphere and lower thermosphere (MLT). Here we look for long-term trends in the annual and semiannual variations of the radio wave absorption in the lower ionosphere, which corresponds to the MLT region heights. Data from central and southeastern Europe are used. A consistent tendency to a positive trend in the amplitude of the semiannual wave appears to be observed. The reality of a similar tendency in the amplitude of the annual wave is questionable in the sense that the trend in the amplitude of the annual wave is probably induced by the trend in the yearly average values of absorption. The phases of both the annual and semiannual waves display a forward tendency, i.e. shift to an earlier time in the year. A tentative interpretation of these results in terms of changes of the seasonal variation of temperature and wind at MLT heights does not contradict the trends observed in those parameters.     

VHF radar observations of gravity waves at a low latitude

Wind observations made at Gadanki (13.5°N) by using Indian MST Radar for few days in September, October, December 1995 and January, 1996 have been analyzed to study gravity wave activity in the troposphere and lower stratosphere. Horizontal wind variances have been computed for gravity waves of period (2–6) h from the power spectral density (PSD) spectrum. Exponential curves of the form e^Z/H have been fitted by least squares technique to these variance values to obtain height variations of the irregular winds upto the height of about 15 km, where Z is the height in kilometers. The value of H, the scale height, as determined from curve fitting is found to be less than the theoretical value of scale height of neutral atmosphere in this region, implying that the waves are gaining energy during their passage in the troposphere. In other words, it indicates that the sources of gravity waves are present in the troposphere. The energy densities of gravity wave fluctuations have been computed. Polynomial fits to the observed values show that wave energy density increases in the troposphere, its source region, and then decreases in the lower stratosphere.     

Spectral characteristics of spring arctic mesosphere dynamics

The spring of 1997 has represented a stable period of operation for the joint University of Tromsø/University of Saskatchewan MF radar, being between refurbishment and upgrades. We examine the horizontal winds from the February to June inclusive and also include estimates of energy dissipation rates derived from signal fading times and presented as upper limits on the turbulent energy dissipation rate, . Here we address the periodicity in the dynamics of the upper mesosphere for time scales from hours to one month. Thus, we are able to examine the changes in the spectral signature of the mesospheric dynamics during the transition from winter to summer states.     

Analysis of gravity waves in the tropical middle atmosphere over La Reunion Island (21°S, 55°E) with lidar using wavelet techniques

The capabilities of the continuous wavelet transform (CWT) and the multiresolution analysis (MRA) are presented in this work to measure vertical gravity wave characteristics. Wave properties are extracted from the first data set of Rayleigh lidar obtained between heights of 30 km and 60 km over La Reunion Island (21°S, 55°E) during the Austral winter in 1994 under subtropical conditions. The altitude-wavelength representations deduced from these methods provide information on the time and spatial evolution of the wave parameters of the observed dominant modes in vertical profiles such as the vertical wavelengths, the vertical phase speeds, the amplitudes of temperature perturbations and the distribution of wave energy. The spectra derived from measurements show the presence of localized quasi-monochromatic structures with vertical wavelengths <10 km. Three methods based on the wavelet techniques show evidence of a downward phase progression. A first climatology of the dominant modes observed during the Austral winter period reveals a dominant night activity of 2 or 3 quasi-monochromatic structures with vertical wavelengths between 1/2 km from the stratopause, 3/4 km and 6/10 km observed between heights of 30 km and 60 km. In addition, it reveals a dominant activity of modes with a vertical phase speed of –0.3 m/s and observed periods peaking at 3/4 h and 9 h. The characteristics of averaged vertical wavelengths appear to be similar to those observed during winter in the southern equatorial region and in the Northern Hemisphere at mid-latitudes.     

Effects of atmospheric oscillations on the field-aligned ion motions in the polar F-region

The field-aligned neutral oscillations in the F-region (altitudes between 165 and 275 km) were compared using data obtained simultaneously with two independent instruments: the European Incoherent Scatter (EISCAT) UHF radar and a scanning Fabry-Perot interferometer (FPI). During the night of February 8, 1997, simultaneous observations with these instruments were conducted at Tromsø, Norway. Theoretically, the field-aligned neutral wind velocity can be obtained from the field-aligned ion velocity and by diffusion and ambipolar diffusion velocities. We thus derived field-aligned neutral wind velocities from the plasma velocities in EISCAT radar data. They were compared with those observed with the FPI (=630.0 nm), which are assumed to be weighted height averages of the actual neutral wind. The weighting function is the normalized height dependent emission rate. We used two model weighting functions to derive the neutral wind from EISCAT data. One was that the neutral wind velocity observed with the FPI is velocity integrated over the entire emission layer and multiplied by the theoretical normalized emission rate. The other was that the neutral wind velocity observed with the FPI corresponds to the velocity only around an altitude where the emission rate has a peak. Differences between the two methods were identified, but not completely clarified. However, the neutral wind velocities from both instruments had peak-to-peak correspondences at oscillation periods of about 10–40 min, shorter than that for the momentum transfer from ions to neutrals, but longer than from neutrals to ions. The synchronizing motions in the neutral wind velocities suggest that the momentum transfer from neutrals to ions was thought to be dominant for the observed field-aligned oscillations rather than the transfer from ions to neutrals. It is concluded that during the observation, the plasma oscillations observed with the EISCAT radar at different altitudes in the F-region are thought to be due to the motion of neutrals.     

Non-magnetic aspect sensitive auroral echoes from the lower E region observed at 50 MHz

Backscatter from E-region irregularities was observed at aspect angles close to 90° (almost parallel to the direction of the magnetic field) using the ALOMAR SOUSY radar at Andoya/Norway. Strong electric fields and increased E-region electron temperatures simultaneously measured with the incoherent scatter facility EISCAT proved that the Farley-Buneman plasma instability was excited. In addition, strong particle precipitation was present as inferred from EISCAT electron densities indicating that the gradient drift instability may have been active, too. Backscatter at such large aspect angles was not expected and has not been observed before. The characteristics of the observed echoes, however, are in many aspects completely different from usual auroral radar results: the Doppler velocities are only of the order of 10 m/s, the half-width of the spectra is around 5 m/s, the echoes originate at altitudes well below 100 km, and they seem to be not aspect-sensitive with respect to the magnetic field direction. We, therefore, conclude that the corresponding irregularities are not caused by the mentioned instabilities and that other mechanism have to be invoked.     

Cooling of the mesosphere and lower thermosphere due to doubling of CO2

A new parameterization of infrared radiative transfer in the 15-m CO₂ band has been incorporated into the Spectral mesosphere/lower thermosphere model (SMLTM). The parameterization is applicable to calculations of heating rates above approximately 15 km for arbitrary vertical profiles of the CO₂ concentration corresponding to the surface mixing ratio in the range 150–720 ppm. The sensitivity of the mesosphere and lower thermosphere (MLT) to doubling of CO₂ has been studied. The thermal response in the MLT is mostly negative (cooling) and much stronger than in the lower atmosphere. An average cooling at the stratopause is about 14 K. It gradually decreases to approximately 8 K in the upper mesosphere and again increases to about 40–50 K in the thermosphere. The cooling and associated thermal shrinking result in a substantial density reduction in the MLT that reaches 40–45% in the thermosphere. Various radiative, chemical, and dynamical feedbacks potentially important for the thermal response in the MLT are discussed. It is noted that the results of simulations are strikingly similar to observations of long-term trends in the MLT. This suggests that during the last 3–4 decades the thermal structure in the real upper atmosphere has undergone substantial changes driven by forcing comparable with that due to doubling of CO₂.