Gravity wave activity in the upper mesosphere over Urbana,Illinois: lidar observations and analysis of gravity wave propagation models

We analyze 375 h of Na Wind/Temperature lidar measurements of the mesopause region (≈ 80–105 km) Na density and temperature profiles on 57 nights distributed over 2 yr at Urbana, Illinois. These observations yield a high-resolution seasonal data set of gravity wave activity in the upper mesosphere. From this data, we present measurements of the Brunt-Väisälä period, the relative atmospheric density perturbations and their spectra, and the parameters of 143 quasi-monochromatic gravity waves. The direct measurement of the Brunt-Väisälä period allows accurate calculation of the horizontal velocity perturbations and vertical displacement perturbations from the density measurements. The horizontal velocity and vertical displacement vertical wave number spectrum magnitudes and indices show considerable seasonal and nightly variability. The gravity wave amplitudes, wavelengths, and observed periods exhibit systematic relationships similar to those found in previous studies, and are consistent with the MU radar measurements of intrinsic gravity wave parameters. Here, we present a detailed analysis of the observations in terms of Diffusive-Filtering Theory models of gravity wave propagation. The magnitudes of the vertical wave number spectrum, the form of the joint vertical wave number and frequency spectrum, and the systematic relationships between the monochromatic gravity wave parameters are consistent with the Diffusive-Filtering model. We compare these results with a variety of radar, lidar, and airglow observations from other sites. This observational study suggests that the complex nonlinear interactions of the gravity wave field may be modeled successfully as a diffusive damping process, where the effective diffusivity is a function of the total wave variance.     

Observational evidence of a saturated gravity wave spectrum in the mesosphere

Five vertical profiles of scalar horizontal winds have been measured at high resolution (25m) in the range from 80–95 km during the last salvo of the MAC/SINE campaign in the summer 1987 at Andenes, Northern Norway (69.3°N). Our purpose in this study is to examine the consistency of the motion spectrum with the saturated spectrum of gravity waves proposed by Smith S. A., Fritts D. C. and Van Zandt T.E., (1987, J. atmos. Sci. 44, 1404). An analysis of vertical wavenumber spectra of the five horizontal wind profiles is presented and it is found that (a) the average slope of the five vertical wavenumber spectra is −3.0 ± 0.2 for wavelengths in the range from 6.4 km to 100 m. The slope is considerably steeper than the vertical wavenumber spectra of the horizontal velocity discussed in the literature, (b) the average vertical wavenumber spectrum shows that there is excellent agreement between the observed spectrum and the saturated spectrum in both slope and amplitude, suggesting that saturation processes do indeed act to control spectral amplitudes at large wavenumbers, and (c) a dominant vertical wavelength of 6.4 km is found in the mesosphere. Taken together, our observations provide further support for the saturated spectrum theory.     

Multiple-frequency studies of the high-latitude summer mesosphere : implications for scattering processes

The characteristics of polar mesosphere summer echoes (PMSE) are studied at 53.5 and 224 MHz. Observations at 2.78 MHz, simultaneous with the ones at the other two frequencies, were carefully compared for indications of PMSE, but no obvious relation was found. Relationships between relative scattering cross-section, spectral width and vertical velocity are studied for the 224 MHz radar, and observations at 53.5 MHz are compared with those at 224 MHz. Results of aspect sensitivity measurements at 53.5 MHz are presented. The implications of these characteristics for several possible scattering mechanisms are discussed. We rule out incoherent scatter and chemically induced fluctuations from the evidence that we have. In view of the extremely low temperatures near the high-latitude mesopause in summer, we discuss several scenarios involving heavy cluster ions and charged aerosol particles.     

Dominant vertical scales of gravity waves in the middle atmosphere observed with the MU radar and rocketsondes

We have simultaneously observed wind motions in the altitude range of 5–90 km by means of the MU radar, rocketsondes and radiosondes. Dominant vertical scales of wind fluctuations due to gravity waves were 2–5 km in the lower stratosphere, about 5–15 km in the upper stratosphere and longer than 15 km in the mesosphere. The increase in the vertical scale with altitude is interpreted in terms of the saturation of upward propagating gravity waves. In the stratosphere, the observed vertical wavenumber spectra showed smaller amplitudes and more gradual slopes than the model values. Furthermore, the wind velocity variance in the stratosphere increases exponentially with an e-folding height of about 9 km, implying that the gravity waves were not fully saturated. On the other hand, the spectra in the upper stratosphere and mesosphere agreed fairly well with the model spectra. The variance in the mesosphere seems to cease increase of the wave amplitudes and agrees reasonably well with the model value.     

Numerical simulations of the saturated gravity wave spectra in the atmosphere

A two dimensional numerical model is used to compute the saturation of small scale gravity waves in the region near the critical level. The vertical wave number spectrum of horizontal velocity fluctuations in the unstable region (USR) where shear instability develops is found to be governed by wave-shear interaction and follows a theoretical saturation spectrum ~ω_b²/2m³. Wave-shear interaction is also found to be responsible for the observed fact that the variance of vertical velocity fluctuations is significantly lower than the level predicted by linear gravity wave theory. On the other hand, the corresponding spectrum in the stable region (SR) following a much shallower spectrum ~m⁻² is found to result from the combined effects of wave-wave interactions and eddy diffusion. The key step in our simulation is the separate parameterization of horizontal and vertical eddy diffusion coefficients instead of a constant molecular viscosity coefficient.     

Contributions of gravity waves to the momentum,heat and turbulent energy budget of the upper mesosphere and lower thermosphere

The role of gravity waves for the momentum and heat budget of the atmosphere between approximately 70 and 110 km height is considered. Parameterization schemes for vertical gravity wave diffusivity, generalized Rayleigh friction, viscous force, heat conduction and kinetic energy dissipation are reviewed. Eddy diffusion parameterization and its relation to the gravity wave approach is also discussed and it is shown that principal similarities exist in both concepts, especially when irregular (stochastic) contributions to the perturbations are modeled. Special attention is paid to the dissipation of perturbation kinetic energy and its contribution to the heat budget of the mesopause region. It is concluded that the amount of energy which can be attributed to the part of the gravity wave spectrum contributing to generalized Rayleigh friction above the mesopause is of the order of 10% of the total perturbation energy.     

Vertical propagation characteristics of internal gravity waves around the mesopause observed by the Arecibo UHF radar

A high resolution wind observation of the mesosphere and lower thermosphere (73–95 km) was conducted with the aid of the high power UHF Doppler radar at Arecibo (18.4°N, 66.8°W). Zonal wind velocities were continuously observed during day-time hours on 1–15 August 1980. We discuss here the observed wind fluctuations with periods of 1–4 h in the light of internal gravity waves. The phase propagation associated with these fluctuations is, on average, shown to be downward, indicating an upward energy flux. A space-time spectral analysis shows that waves with vertical wavelengths shorter than 10 km disappear around the mesopause (about 85km), while those with longer vertical wavelengths exist throughout the observational height. This result is explained in terms of wave absorption at a critical layer where the mean zonal wind has a westerly shear with height. This feature is consistent with the behavior expected for internal gravity waves around the summer mesopause in order to explain general circulation models.     

Observational evidence of local turbulence production due to gravity wave saturation in the summer mesosphere

Five foil chaff and two falling sphere rockets flown during the MAC/SINE Campaign on 15 July 1987 at Andenes, Northern Norway (69°17′N). From these rocket measurements, turbulent energy dissipation rates, vertical wind shears and Richardson numbers as functions of height were derived in the range from 82 to 92km. Turbulent energy dissipation rates generally range from 1.4 × 10⁻⁵ to 2.0 × 10⁻²W/kg and are consistent with other experiments performed at the same latitude. Strong wind shears of the order of 50–90 m/s/km are observed at various heights. Good correspondence between turbulence intensity peaks, regions of strong wind shear and low Richardson number is found. Vertical wavenumber spectra of the five scalar winds measured by the foil chaff rockets indicate that there is an excellent agreement with the saturation hypothesis, suggesting that the turbulence intensity peaks measured in this salvo are linked directly to the saturation of gravity wave motions via dynamical instabilities.     

Gravity waves from O2 nightglow during the AIDA '89 campaign III: effects of gravity wave saturation

The O₂ atmospheric (0–1) night airglow emitted within the gravity wave saturation region at ∼90–100 km can serve as a means of studying the wave activity. In this analysis, the atmospheric motions were described by a mean spectral model and an algorithm was developed to infer the wave kinetic energy density and momentum flux from variations in O₂ (0–1) airglow emission rate and rotational temperature. The method was applied to eight nights of data collected by MORTI, a mesopause oxygen rotational temperature imager, during the AIDA campaign of 1989 in Puerto Rico (18°N, 67°W). The observed r.m.s. fractional fluctuations of airglow emission rate and rotational temperature were of the order of ∼0.07–0.15 and ∼0.02–0.04, respectively, and the characteristic vertical wavelengths were estimated at ∼10 2 -20 km. The inferred r.m.s. horizontal velocities and velocity variances were found to be ∼12–25 m/s and ∼150–600 m²/s², with the majority of the horizontal velocity and its variance associated with low-frequency, large-scale wave motions. The estimated momentum fluxes, mainly contributed by high-frequency, small-scale waves, were ∼2–10 m²/s². These results are in good agreement with those obtained from other measurements using different observational methods at low and mid-latitudes.     

Simultaneous nightglow and Na lidar observations at Arecibo during the AIDA-89 campaign

Simultaneous nightglow and Na lidar observations made at Arecibo during the AIDA-89 campaign are reported from between 0000 and 0530 LT (Local Time) on 6 April 1989 and on 9 April 1989. On 9 April the observations are consistent with the presence of a large amplitude 8.5 km vertical wavelength, 4 h period wave propagating through the 85–91 km region. A lower amplitude 80 min period wave is also observed. The results imply that the O₂ atmospheric band intensity peaked near 91 km while the OH Meinel (6,2) band intensity peaked near 85 km. The OH Meinel temperature and intensity are 160° out of phase which can be explained by low eddy diffusion and high ozone densities near the mesopause. The integrated Na abundance from 87 to 89 (90–92) km correlates well with the OH Meinel (O₂ atmospheric) band intensity. On 5–6 April the OH Meinel and O₂ atmospheric band intensities are not well correlated. The OH Meinel intensity is correlated with the integrated Na abundance from 86 to 88 km. Both the Na and OH measurements reveal the presence of an approximately 1 h period wave. The OH temperature data appear to be consistent with the OH Meinel band originating near 85 km. The O₂ atmospheric band data show the presence of a 2 h period wave. The integrated Na abundance data suggest that the O₂ atmospheric band peaks between 90 and 94 km. A large sporadic Na event which occurs near 6 UT appears related to the presence of a gravity wave near 95 km. In all of the observed waves there is good agreement between the wave parameters derived separately by the optical airglow and Na lidar techniques.     

The non-linear interaction of a gravity wave with the vortical modes

The wave-wave interaction theory has been used successfully in describing one class of weakly non-linear wave phenomena. The application of this theory to the atmosphere shows the possibilities of energy and momentum transfer among three interacting gravity waves, as well as from the gravity wave to the other modes of motion. It has been found that the non-resonant interaction of a gravity wave with two vortical modes can proceed at a reasonably rapid rate. With the gravity wave viewed as the primary wave and the two vortical modes as the secondary waves, the interaction equation can be linearized and solved. The resulting analytic formula gives the growth rate of the interaction. In the absence of the Earth's rotation, the growth is limited to a threshold effect. The theory shows that whenever the horizontal air parcel velocity of a gravity wave exceeds a factor of √2 times the horizontal trace velocity of the wave, energy and momentum transfer from the gravity wave to the vortical modes can proceed. The rotation of the Earth will blur this threshold effect by making the interaction more likely to occur. Thus, through this mechanism, a gravity wave can transfer its energy and momentum to the horizontal velocity field in the vortical mode. In this sense, the small scale vortical motions would serve as the sink of both energy and momentum of a propagating gravity wave. When scales of vortical modes reach sufficiently small values, dissipation through viscosity becomes important. At this scale and smaller, the vortical modes are damped out quickly and its energy spectrum must exhibit a sharp decay.     

Equatorial spread-F and neutral atmospheric turbulence : a review and a comparative anatomy

A number of satellite and rocket plasma density spectra obtained during equatorial spread-F conditions are presented and discussed in the light of similar measurements in the neutral atmosphere. We discuss this comparison in some detail and find both distinct similarities and subtle differences. The horizontal spectral measurements show a peak at an outer scale quite similar to the scale of the undulations caused by gravity wave interactions with the ionosphere. This feature is similar to a buoyancy subrange but it is easy to show that the amplitudes of the plasma fluctuations are too large to be directly driven by the neutral atmosphere. At intermediate scales the plasma fluctuations have a one-dimensional horizontal spectrum with a power law well described by a ($\text{−5}\text{3}$) slope. Once again it can be shown that the neutral fluid cannot be similarly structured at 400 km altitude due to the high viscosity coefficient. A plasma cascade process seems to be operating but it is not at all clear how the spectrum is formed. Furthermore, vertical power spectra seem to run the gamut in spectral form with slopes (n) varying in the range from −1 to −3. So the horizontal spectra are near universal in form, while the vertical spectra are quite variable.     

Simulation of gravity wave effects under solstice conditions using a 3-D circulation model of the middle atmosphere

Vertically propagating gravity waves can transport momentum and energy from the troposphere up to the mesosphere and thus modify the circulation of the middle atmosphere. The effects of regional gravity wave sources, together with temporal changes of gravity wave activity, are studied under solstice conditions in a 3-D circulation model using a simplified parameterization scheme for the gravity momentum deposition. In this way we can reproduce the reversal of the mean zonal wind with height and very low temperatures at the summer mesopause region. Using a stochastic forcing by taking the gravity wave parameters at random, characteristic oscillations are found with periods in the planetary scale range (2, 4 and 5 days) and in the tidal range (1 day, 16 h and 12 h).     

Gravity-wave motions in the mesosphere and lower thermosphere observed at Mawson,Antarctica

We present the results of MF radar observations of mean winds and waves in the height range 78–108 km at Mawson (67°S, 63°E), Antarctica. The measurements were made in the period from 1984 to 1990. Climatologies of the prevailing zonal and meridional circulations made with a 12-day time resolution show that the mean circulation remained relatively stable over the 6 yr of observation. Climatologies of gravity-wave motions in the 1–24 h period range were also generated. These reveal that the r.m.s. amplitudes of horizontal wave motions near the mesopause (~90 km) are about 30 m s⁻¹, and that there is some anisotropy in the motions, especially at heights below 90 km. Meridional amplitudes are larger than zonal amplitudes, which suggests a preference for wave propagation in the north-south direction. Comparisons with MST radar wind observations made near the summer solstice at Poker Flat, Alaska (65°N) and at Andøya, Norway (69°N) show similarities with the Mawson observations, but the wave amplitudes and mean motions are larger in magnitude at the northern sites. This suggests hemispheric differences in wave activity that require further study.     

Interference of tidal and gravity waves in the ionosphere and an associated sporadic E-layer

Observations made on 10 July 1987 with the EISCAT UHF radar are presented. The F-region measurements of both electron density and field-aligned ion velocity show that an upward propagating gravity wave with a period of about 1 h is present. The origin of the gravity wave is probably auroral. The E-region ion velocities show a tidal wave and both upward and downward propagating gravity waves. The gravity waves have three dominant periods with a possible harmonic relationship and similar vertical wavelengths. These waves are either reflected at a single reflection level, ducted between two levels, or they are generated in a non-linear interaction between gravity and tidal waves. The E-region electron density is dominated by particle precipitation. After a short burst of more intense precipitation, a sporadic E-layer forms at 105km and then disappears 40min later. Within this time, the layer rises and falls by a few kilometres, following closely the motion of a convergent null in the velocity profile. We suggest that the formation and destruction of this layer is controlled by both the precipitation, which indirectly provides a source of metal ions through charge exchange, and the superposition of gravity waves and the tidal wave.     

Inertio-gravity waves and subtropical multiple tropopauses: Vertical wavenumber spectra of wind and temperature observed by the MU radar,radiosondes and operational rawinsonde network

We have carried out continuous observations of the tropopause region over Japan for three weeks during the Baiu (early summer rain “in Japan”) season in 1991, by using a VHF Doppler radar (the MU radar), radiosondes launched at the radar site and operational rawinsondes at five meteorological stations. Based on these observations, we try to examine the hypothesis that the multiple tropopauses and the dominant inertio-gravity waves are one and the same feature, and obtain some interesting results that are not inconsistent with this hypothesis. First, vertical wavenumber spectra and hodographs analyzed from the radar wind data in the tropopause region suggest that inertio-gravity waves with vertical wavelengths of ∼ 2 km are quasi-monochromatically dominant (with 2–3 day scale variabilities of 10–20%), and are in accordance with activities of the subtropical jet stream and mesoscale cyclone-front system activities observed by the operational network. Second, striking (potential) temperature fluctuations are detected simultaneously by the radiosondes and rawinsondes, which appear as multiple tropopauses in meridional cross-section analysis. Third, vertical wavenumber spectra analyzed from the radiosonde temperature data are consistent with the radar wind spectrum, if we assume that both wind and temperature fluctuations are mainly induced by the dominant inertio-gravity waves. Finally, we confirm that the dominant interio-gravity waves can be barely detected also from routine rawinsonde (1.5-km running-mean) wind data if the amplitude is larger than 1.5 m/s. However, the monochromatic wave structures are generally quite localized in space and time.     

Stratospheric and tropospheric wave measurements near the Andes mountains

Experimental results and interpretation of temperature, pressure and wind velocity measurements, performed with an instrumented balloon, are presented. The balloon, an open-type stratospheric one was launched from Mendoza (Argentina), near the Andes mountains. The data analysis suggests the presence of a large amplitude quasi-inertial gravity wave, with intrinsic period close to 0.5 days, and vertical wavelength of around 1.7 km just below the tropopause. The possible orographic origin of this wave is discussed. A Fourier analysis confirms the existence of this mode, simultaneously in the temperature and in the wind velocity components. A hodograph of the zonal and meridional wind components shows the expected counterclockwise sense of rotation of the horizontal velocity with increasing altitude, corresponding to a long period gravity wave, in the southern hemisphere.It is found that the vertical wind velocity variations measured by the anemometer, are mainly due to buoyancy force variations induced by the wave on the open stratospheric balloon. The vertical profile obtained by the anemometer is anticorrelated with the rate of ascent or descent of the gondola. As a consequence, the wave induced velocity is very difficult to obtain using these balloons, contrary to the case of radiosonde balloon data. The differences in the response of open stratospheric and radiosonde-type balloons to the presence of internal gravity waves may be explained by their different design and material characteristics.     

Gravity wave saturation and eddy diffusion in the middle atmosphere

A discussion is given of gravity wave saturation and its relation to eddy diffusion in the middle atmosphere. Attention is focused on the saturation process and some of its observable manifestations. It does not serve as a review of all related work. Although a theoretical point of view is taken, the emphasis is on which wave parameters need be measured to predict quantitatively the influence of gravity waves on eddy transport. The following considerations are stressed: the variation of spectra with observation time T; that eddy diffusivities are determined by velocity spectra; the anisotropic nature of diffusivity; a unified approach to saturation; an attempt to make eddy diffusivity more precise; the relationship between eddy diffusivity and wave dissipation.The subjects of ‘wave drag’ (momentum flux deposition) and heat flux need only be treated briefly, because they are related to eddy diffusivity in simple ways. Consideration is also given to two different theoretical mechanisms of wave saturation—wave induced convective instability and strong nonlinear wave interactions. The saturation theory is then used to predict a globally averaged height profile of vertical diffusivity in the middle atmosphere. This calculation shows that gravity waves are a major contributor to eddy diffusion from heights of 40–110 km, and that they are significant down to 20 km. A more detailed calculation of wave induced eddy diffusion, including latitudinal and seasonal variations, can be made if wave velocity spectra become available. The paper closes with recommendations for future research.     

Summertime stratospheric wind measurements above the South Pole

Observations of the mean wind flow and wave motions in the stratosphere at the South Pole are presented. The atmospheric motions are determined from the tracking of a high altitude, zero-pressure balloon launched from Amundsen-Scott Station during the austral summer of 1985–1986. The balloon position was precisely monitored by an optical theodolite for a large portion of the flight so that small scale motions could be resolved. The mean flow above the pole was approximately 3ms⁻¹. Atmospheric motions characteristic of internal gravity waves were observed with an intrinsic period of approximately 4.5 h and vertical and horizontal wavelengths of approximately 2.5km and 125km, respectively. The horizontal perturbation velocity of the observed waves was large compared to the mean horizontal flow velocity. The implication is that wave motions play a dominant role in the transport of stratospheric constituents in regions where the mean winds are light, such as over the South Pole during austral summer.     

Design Earthquake Ground Motion Prediction for Perth Metropolitan Area with Microtremor Measurements for Site Characterization

Perth is the largest city in Western Australia and home to three-quarters of the state's residents. In recent decades, there have been a lot of earthquake activities just east of Perth in an area known as the South-West Seismic Zone. Previous numerical results of site response analyses based on limited available geology information for PMA indicated that Perth Basin might amplify the bedrock motion by more than 10 times at some frequencies and at some sites. Hence, more detailed studies on site characterization and amplification are necessary. The microtremor method using spatial autocorrelation (SPAC) processing is a useful tool for gaining thickness and shear wave velocity (SWV) of sediments and has been adopted in many previous studies. In this study, the response spectrum of rock site corresponding to the 475-year return period for PMA is defined according to the probabilistic seismic hazard analysis (PSHA) based on the latest ground motion attenuation model of Southwest Western Australia. Site characterization in PMA is performed using two microtremor measurements, namely SPAC technique and H/V method. The clonal selection algorithm (CSA) is introduced to perform direct inversion of SPAC curves to determine the soil profiles of representative PMA sites investigated in this study. Using the simulated bedrock motion as input, the responses of the soil sites are estimated using numerical method based on the shear-wave velocity vs. depth profiles determined from the SPAC technique. The response spectrum of the earthquake ground motion on surface of each site is derived from the numerical results of the site response analysis, and compared with the respective design spectrum defined in the Australian Earthquake Loading Code. The comparison shows that the code spectra are conservative in the short period range, but may slightly underestimate the response spectrum at some long period range.