Statistics of gravity waves in the lower stratosphere over Beijing based on high vertical resolution radiosonde

Atmospheric gravity waves, with small to medium scales, prevail in the atmosphere and have global ef- fects. Many researches show that gravity waves are the main source that causes the variation of wind and temperature field in the stratosphere, and that the break-up of upward propagating gravity waves is the dominant sources of small scale turbulent and mixing processes in the middle atmosphere. Theories and ob- servations indicate that the redistribution of momen- tum, caused by the generati…     

The structure of the natural stratosphere and the impact of chlorofluoromethanes on the ozone layer investigated in a 2-D time dependent model

two-dimensional time dependent model of the stratosphere incorporating the major interactions between radiative-photochemical and dynamical processes is described. The main prognostic equations considered are the thermodynamic equation and the general conservation equation for the minor chemical constituents representing the odd oxygen (O _x =O+('D)+O₃), odd hydrogen (HO _x =HO+HO₂), N₂O, odd nitrogen (NO _x =NO+NO₂+HNO₃), CF₂Cl₂, CFCl₃ and odd chlorine (Cl _x =Cl+ClO+HCl). The zonal wind and mean meridional circulations are determined diagnostically by the integration of the thermal wind equation and the stream function equation in the meridional plane espectively. The large scale eddy processes are parameterized in terms of zonal mean quantities using the generalized diffusion formulation on a sloping surface. The radiative heating and cooling and the hotochemical sources and sinks are incorporated in a form which allows for the major interactions among the minor trace constituents, temperature and mean circulation.Two integrations consisting of natural stratosphere and a stratosphere contaminated by the chlorofluoromethanes through lower boundary fluxes are carried out for 23 model years by changing the declination of the sun every day and using 6-hour time step. The model simulations of temperature, mean circulation, ozone, HO _x , N₂O and NO _x in the meridional plane for the normal stratosphere, show satisfactory agreement with the available observations. Based on the results of second integration it is found that the injection of chlorofluoromethanes in the atmosphere at the estimated current production rates can lead to significant changes in the meridional distribution of ozone, temperature and NO _x in the middle and upper stratosphere. The results also indicate that the percentage total ozone depletion increases from tropics to high latitudes and from summer to winter high latitudes. Also discussed are the results of additional experiments incorporating the reaction of HO₂ with NO and the reactions involving ClNO₃.     

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.     

Mean meridional circulations of the stratosphere and mesosphere

Observations from the Nimbus 6 pressure modulator radiometer (PMR) have been used to estimate monthly mean planetary wave fluxes of heat and momentum in the stratosphere and mesosphere. While the eddy heat fluxes play an important role in the mean meridional circulation of the winter stratosphere they are shown to be less important in the upper mesosphere. Incorporation of the observed momentum fluxes into the Oxford two-dimensional circulation model has shown that they are incapable of providing the momentum transport necessary to balance the zonal flow accelerations induced by the mean meridional motion. Other unspecified transfer processes represented by Rayleigh frictional damping of the zonal fow are shown to dominate. In contrast the observed fluxes in the stratosphere achieve the necessary redistribution of momentum. Moreover their interannual variability profoundly influences the stratospheric circulation, as demonstrated in the model by the use of two different annual sets of observed momentum fluxes. The desirability of calculating the planetary wave behaviour within the model is indicated.     

Quasigeostrophic planetary waves in a two-layer ocean with one-dimensional periodic bottom topography

Although the study of topographic effects on the Rossby waves in a stratified ocean has a long history, the wave property over a periodic bottom topography whose lateral scale is comparable to the wavelength is still not clear. The present paper treats this problem in a two-layer ocean with one-dimensional periodic bottom topography by a simple numerical method, in which no restriction on the wavelength and/or the horizontal scale of the topography is required. The dispersion diagram is obtained for a wavenumber range of [?π/L _b , π/L _b ], where L _b is the periodic length of the topography. When the topographic?β?is not negligible compared to the planetary β, the Rossby wave solutions around the wavenumbers which satisfy the resonant condition among the waves and topography disappear and separate into an infinite number of discrete modes. For convenience, each mode is numbered in order of frequency. As topographic height is increased, the high frequency barotropic Rossby wave (mode 1) becomes a topographic mode which can exist even on the f plane, and the highfrequency baroclinic mode (mode 2) becomes a surface intensified mode. Behaviors of low frequency modes are somewhat complicated. When the topographic amplitude is small, the low frequency baroclinic modes tend to be bottom trapped and the low frequency barotropic modes tend to be surface intensified. As topographic amplitude further increases, the relation between the mode number and vertical structure changes. This change can be attributed to the increase of the frequency of the topographic mode with the topographic amplitude.     

Effects of wave-wave and wave-mean flow interactions on the evolution of a baroclinic wave

Abstract

The weakly nonlinear evolution of a free baroclinic wave in the presence of slightly supercritical, vertically sheared zonal flow and a forced stationary wave field that consists of a single zonal scale and an arbitrary number of meridional harmonics is examined within the context of the conventional two-layer model. The presence of the (planetary-scale) stationary wave introduces zonal variations in the supercriticality and is shown to alter the growth rate and asymptotic equilibrium of the (synoptic-scale) baroclinic wave via two distinct mechanisms: The first is due to the direct interaction of the stationary wave with the shorter synoptic wave (wave-wave mechanism), and the second is due to the interaction of the synoptic wave with that portion of the mean field that is corrected by the zonally rectified stationary wave fluxes (wave-mean mechanism). These mechanisms can oppose or augment each other depending on the amplitude and spatial structure of the stationary wave field. If the stationary wave field is confined primarily to the upper (lower) layer and consists of only the gravest cross-stream mode, conditions are favorable (unfavorable) for nonzero equilibrium of the free wave.

In addition to the time dependent heat flux generated by baroclinic growth of the free wave, its interaction with a stationary wave field consisting of two or more meridional harmonics generates time dependent heat fluxes that vary with period of the free wave. However, if the stationary wave field contains several meridional harmonics of sufficiently large amplitude, the free baroclinic wave is destroyed.     

Finite-amplitude mountain waves in a compressible atmosphere including the effects of dissipation

Abstract

Adiabatic, two-dimensional, steady-state finite-amplitude, hydrostatic gravity waves produced by flow over a ridge are considered. Nonlinear self advection steepens the wave until the streamlines attain a vertical slope at a critical height z_c. The height z_c , where this occurs, depends on the ridge crest height and adiabatic expansion of the atmosphere. Dissipation is introduced in order to balance nonlinear self advection, and to maintain a marginal state above z_c. The approach is to assume that the wave is inviscid except in a thin layer, small compared to a vertical wavelength, where dissipation cannot be neglected. The solutions in each region are matched to obtain a continuous solution for the streamline displacement δ. Solutions are presented for different values of the nondimensional dissipation parameter β. Eddy viscosity coefficients and the thickness of the dissipative layer are expressed as functions of β, and their magnitudes are compared to other theoretical evaluations and to values inferred from radar measurements of the stratosphere.

The Fourier spectrum of the solution for z ≫ z_c is shown to decay exponentially at large vertical wave numbers n. In comparison, a spectral decay law n ^?-8/3 characterizes the marginal state of the wave at z = z_c .     

Instability of planetary waves in a high vertical resolution model

Abstract

A high vertical resolution model is used to examine the instability of a baroclinic zonal flow and a finite amplitude topographically forced wave. Two families of unstable modes are found, consisting of zonally propagating most unstable modes, and stationary unstable modes. The former have time scale and spatial structure similar to baroclinic synoptic disturbances, but are localized in space due to interaction with the zonally asymmetric forcing. These modes transport heat efficiently in both the zonal and meridional directions. The second family of stationary unstable modes has characteristics of modes of low frequency variability of the atmosphere. They have time scales of 10 days and longer, and are of planetary scale with an equivalent barotropic vertical structure. The horizontal structure resembles blocking flows. They are maintained by available potential energy of the basic wave, and have large zonal heat fluxes. The results for both families of modes are interpreted in terms of an interaction between forcing and baroclinic instability to create favoured regions for eddy development. Applications to baroclinic planetary waves are also considered.     

Climatological Ertel’s potential-vorticity flux and mean meridional circulation in the extratropical troposphere – lower stratosphere

We investigated to what extent the isentropic, non-geostrophic formulation of zonally averaged circulation derived for stratospheric conditions is applicable to climatological transport in the extratropical troposphere and lower stratosphere. The study is based on 10 years of daily data of ECMWF analysis and on the ECHAM3 climate model of the German Climate Computing Centre. The main result is a scalar isentropic mixing coefficient, K_yy, and a mean meridional transport circulation consistently derived from the same data base. For both data sources, isentropic mean meridional circulation is derived from horizontal mass flow rate for 4 representative months. Alternatively, a mean meridional circulation is calculated from total diabatic heating rates of the ECHAM3 model. It is shown that only the latter is in good agreement with the ECMWF mean meridional circulation. Isentropic analysis also comprises the seasonal cycle of the climatological meridional gradient and flux of Ertels potential vorticity (PV). Application of Tungs flux-gradient relation yields that for all seasons K_yy is positive in height-latitude regions where statistical significance is reached. Large K_yy values, marking regions of more efficient mixing, have been found in the subtropical vertical band of weak westerly wind and in mid-latitudes in regions of upward-propagating baroclinic wave activity in the middle and upper troposphere. Based on the ECMWF data and results of baroclinic-wave behaviour, strong indications are presented that positive zonally averaged PV flux polewards of the jet core in the NH is strengthened by stationary waves and nonlinear effects. Reduced eddy transport is apparent in winter and spring slightly below the subtropical tropopause jet. The seasonal cycle of K_yy from ECHAM3 data is to a great extent in agreement with the result based on ECMWF analysis. In the model, reduced interannual variability enlarges the height-latitude range where sign of K_yy is significant.     

On the normal mode instability of modons and Wu-Verkley waves

Abstract

The normal mode instability of steady Wu-Verkley (1993) wave and modons by Verkley (1984, 1987, 1990) and Neven (1992) is considered. All these flows are solutions to the vorticity equation governing the motion of an ideal incompressible fluid on a rotating sphere. A conservation law for infinitesimal perturbations to each solution is derived and used to obtain a necessary condition for its exponential instability. By these conditions, Fjörtoft's (1953) average spectral number of the amplitude of an unstable mode must be equal to a specific number that depends on the degree of the solution in its inner and outer regions as well as on spectral distribution of the mode energy in these regions. Some properties of the conditions for different types of modons are discussed. The maximum growth (and decay) rate of the modes is estimated, and the orthogonality of the amplitude of each unstable, decaying, or non-stationary mode to the basic solution is shown in the energy inner product.

The new instability conditions confine the unstable disturbances of the WV wave and modon to a hypersurface in the perturbation space and allow interpretation of their energy structure. They are also useful both in estimating the maximum growth rate of unstable modes and in testing the numerical algorithms designed for the linear stability study.     

On Rossby wave propagation in a meridionally stratified channel

Rossby wave propagation in the presence of a nonseparable Brunt-Väisälä frequency,N(y,z), and the associated geostrophic zonal flow,U(y,z), is examined in this paper. The usual quasi-geostrophic potential vorticity equation only includes vertical variations in Brunt-Väisälä frequency (i.e.N(z)). We derive a linearised quasi-geostrophic potential vorticity equation which explicitly includesN(y, z), where variations inN may occur on the internal Rossby radius length scale. A mixed layer distribution that monotonically deepens in the poleward direction leads to a nonseparableN(y,z). The resulting meridional pressure gradient is balanced by an eastward zonal geostrophic flow.By assuming mixed layer depth changes occur slowly, relative to a typical horizontal wavelength of a Rossby wave, a local analysis is presented. The Rossby wave is found to have a strongly modulated meridional wavenumber,l, with amplitude proportional to |l|^–1/2. To elucidate whether the modulations of the Rossby wave are caused by the horizontal variations inN orU we also consider the cases where eitherN orU vary horizontally. Mixed layer depth changes lead to largestl where the mixed layer is deepest, whereasl is reduced in magnitude whereU is nonzero. When bothU(y,z) andN(y,z) are present, the two effects compete with one another, the outcome determined by the size of |c|/U _max, wherec is the Rossby wave phase speed. Finally, the slowly varying assumption required for the analytical approach is removed by employing a numerical model. The numerical model is suitable for studying Rossby wave propagation in a rectangular zonal channel with generalN(y, z) andU(y, z).     

Nonhydrostatic Atmospheric Normal Modes on Beta-Planes

--To facilitate the understanding of nonhydrostatic effect in global and regional nonhydrostatic models, the normal modes of a nonhydrostatic, stratified, and compressible atmosphere are studied using Cartesian coordinates on midlatitude and equatorial #-planes. The dynamical equations without forcing and dissipation are linearized around the basic state at rest, and solved by using the method of separation of variables. An eigenvalue-eigenfunction problem is formulated, consisting of the horizontal and vertical structure equations with suitable boundary conditions. The wave frequency and the separation parameter, referred to as "equivalent height," appear in both the horizontal and vertical characteristic equations as a coupled problem, unlike the hydrostatic case. Therefore, the nonhydrostatic equivalent height depends not only on the vertical modal scale, as in the hydrostatic case, but also on the zonal and meridional modal scales. Numerical resu lts on the dispersion relations are presented for an isothermal atmosphere. Three kinds of normal modes, namely acoustic, gravity, and Rossby modes, are solved and compared with the corresponding global solutions. Nonhydrostatic effects are studied in terms of normal modes in a wide range of wavelengths from small to planetary scales. It is demonstrated that Rossby modes are hardly affected by nonhydrostatic effects regardless of wavelengths. However, nonhydrostatic effects on gravity modes become significant for smaller horizontal and deeper vertical scales of motion. The equivalent height plays a particularly important role in evaluating nonhydrostatic effects of normal modes on the equatorial #-plane, because the equivalent height appears in the scaling of meridional distance variable of the eigenfunctions. The implementation of nonhydrostatic normal mode analysis on high-resolution numerical modeling is also discussed.     

An evaluation of the classical and extended Rossby wave theories in explaining spectral estimates of the first few baroclinic modes in the South Pacific Ocean

Previous literature has suggested that multiple peaks in sea level anomalies (SLA) detected by two-dimensional Fourier Transform (2D-FT) analysis are spectral components of multiple propagating signals, which may correspond to different baroclinic Rossby wave modes. We test this hypothesis in the South Pacific Ocean by applying a 2D-FT analysis to the long Rossby wave signal determined from filtered TOPEX/Poseidon and European Remote Sensing-1/2 satellite altimeter derived SLA. The first four baroclinic mode dispersion curves for the classical linear wave theory and the Killworth and Blundell extended theory are used to determine the spectral signature and energy contributions of each mode. South of 17°S, the first two extended theory modes explain up to 60% more of the variance in the observed power spectral energy than their classical linear theory counterparts. We find that Rossby wave modes 2–3 contribute to the total Rossby wave energy in the SLA data. The second mode contributes significantly over most of the basin. The third mode is also evident in some localized regions of the South Pacific but may be ignored at the large scale. Examination of a selection of case study sites suggests that bathymetric effects may dominate at longer wavelengths or permit higher order mode solutions, but mean flow tends to be the more influential factor in the extended theory. We discuss the regional variations in frequency and wave number characteristics of the extended theory modes across the South Pacific basin.     

Further ozone transport calculations and the spring maximum in ozone amount

Summary Calculations of the covariance between ozone amounts and meridional wind in the lower stratosphere are presented for all stations in the northern hemisphere for the IGY-IGC. Northward ozone transport occurs by large-scale quasi-horizontal transient and standing eddies and the transport is a maximum early in the year. It is suggested that the transport is governed by the exchange of energy between the troposphere and stratosphere and data are presented on the energy transformations within the lower stratosphere and the transfer of energy into the region which support this suggestion. The vertical flux of energy is also calculated from tropospheric data and its seasonal changes are seen to be in the correct phase to explain the spring maximum in ozone amount.The research reported in this article was sponsored by the Atomic Energy Commission under Contract AT (30-1)2241.     

Stratigraphic filtering of seismic waves: The influence of mode conversion multiples onP-wave andS-wave phase and attenuation

We demonstrate how multiples, generated at the interfaces of plane parallel beds, modify the propagation characteristics of an originally coherent seismic wave. For waves propagating at an angle to the bedding plane we find that theSV andP-waves couple so that neither is a pure mode. TheSH-wave, while modified in its propagation characteristics by multiples, remains a pure mode. The coupling ofSV-multiples into the quasi-P-mode appears weaker than the coupling ofP-wave multiples into the quasi-SV mode; at least this is so for the two simple cases of (a) density fluctuations only and (b) correlatedV _p andV _s fluctuations which conserve Poisson's ratio.We also find that the coupling is sensitive to both the angle of propagation and frequency. In addition there is a cut-off angle forP-wave multiples influencing the quasi-SV mode. Propagation angles larger than the cut-off permit theP-multiples to modify the phase of the quasi-SV mode, but not its effective attenuation. No such cut-off effect is found for SV-multiples influencing the quasi-P mode, whose angle-dependent and frequency-dependent phase distortion and effective attenuation are influenced both byP-wave multiples andSV-multiples.In view of the mathematical complexity of the expressions describing the phase, and effective attenuation of modes when allowance is made forP-andS-wave multiples, we strongly advocate numerical coding of the major mathematical formulae. By so doing a systematic study can be undertaken of the frequency and offset dependence of seismic waves as a function of seismic source input and power spectral behavior of the fluctuations in density and elastic constants of beds. It is our opinion that the full mathematical expressions are too involved to permit an analytic, systematic investigation to be given of the phase and attenuation of seismic waves with any degree of sophistication or generality.     

On the plausible linkage of thermospheric meridional winds with the equatorial spread F

Some of the characteristic features of thermospheric meridional winds during equinoctial period, associated with equatorial spread F (ESF) and their possible role in the triggering of ESF are presented through case studies of observational events under different geophysical conditions that essentially control the post-sunset F-layer height (h′F) rise. The present study reveals that the polarity and magnitude of the meridional winds become significant with the equatorward wind being present when the h′F is below a critical height for the instability to get triggered. The distinctly different characteristic features of the meridional winds during ESF and non-ESF events are presented and discussed.     

Vertical and off-diagonal eddy diffusivities in the Northern Hemispheric stratosphere and lower mesosphere

Meteorological rocket soundings, launched between 1968–74 at six locations representative of low, middle, and high latitudes in the northern hemisphere, are employed to determine the vertical, meridional and off-diagonal components of the eddy diffusivity in the northern hemispheric statosphere and lower mesosphere.It is shown that the distribution of the vertical and meridional components of the eddy diffusivity are similar in the northern hemisphere, although the magnitude of the former is 10⁷ smaller than that of the latter; the magnitude of the off-diagonal eddy diffusivity is about 10³ smaller than that of the meridional eddy diffusivity. In the troposphere, a maximum eddy diffusivity occurs in the mid-latitude at about 7 km above the mean sea level for both the summer and winter seasons. In the stratosphere, a maximum eddy diffusivity occurs in the mid-latitude at about 33 km in the winter, but no maximum in the summer.Paper presented at the World Meterological Organization Technical Conference on Global Observations of Atmospheric Pollution Relative to Climate, Boulder, Colorado, 20–24 August 1979.     

Lagrangian motion of air parcels in the stratosphere in the presence of planetary waves

Development of thoughts on tracer transport mechanisms in the stratosphere which lead to new approaches to two-dimensional modeling of the tracer problem is reviewed.Three-dimensional motions of individual air parcels affected by a planetary wave are investigated theoretically, treating a steady, upward propagating wave in a uniform flow. It is shown that trajectories of air parcels are of elliptical form when projected onto the meridional plane and that they have no mean meridional or vertical motion, even though the usual zonal Eulerian-mean vertical motion exists. The origin of the difference between the mean air parcel motion and the Eulerian-mean motion is discussed.On the basis of the knowledge of air parcel motion, two approaches to two-dimensional modeling are considered. The generalized Lagrangian mean motion (quasi-zonal weighted mean taken over a meandering material tube), recently introduced by Andrews and McIntyre, is identical with the mean motion of an air parcel in a steady state. Such a mean meridional circulation may be used for advecting a tracer in the meridional plane in a two-dimensional model. The transport effect is represented solely by the advection and an eddy transport does not appear in this scheme, to a first approximation.The finding that trajectories of air parcels are elliptical necessitates a reexamination of the Reed-German eddy diffusivity currently used in two-dimensional chemical-dynamical models. By applying a mixing length type hypothesis, we derive an eddy diffusivity formula for use in Eulerian-mean calculations, which, in the case of a conservative tracer is dominated by an anti-symmetric tensor. The eddy transport due to this anti-symmetric tensor diffusivity is of advective type (not diffusive) and has the effect of taking the Stoke drift effect into account, when used in the usual Eulerian-mean formulation.     

The meridional mode of global tropospheric–stratospheric circulation

The spatio–temporal variations of major meridional modes are studied by using the ECMWF and NCEP/NCAR reanalyzed geopotential data between 70 and 10 hPa during 1979 and 2001. The variance contribution rates from the first and second modes are 56–69% and 14–22%, respectively, for ECMWF, but 76–85% and 9–10% for NCEP/NCAR. The climatic trend coefficients are positive and negative in the troposphere and stratosphere, respectively. The reversal is remarkably correlated with the AO/NAM, AAO/SAM and polar vortex, suggesting their important role in connecting the mid-low and upper circulation and the interaction between SH and NH.     

Quasi-free-decay magnetic modes in planetary cores

A new category of hydromagnetic waves in a rotating conducting fluid within a spherical shell geometry is investigated. These quasi-free-decay magnetic modes are based on particular solutions of the induction equation where the magnetic diffusion plays the central role. These solutions, normally only decaying with time, become propagative owing to the combined action of the background magnetic field and the rotation. The amplitude and sign of their azimuthal phase drift strongly depend on morphology and magnitude of the background magnetic field. The validity domain of these quasi-free-decay (QFD)-modes is related to the Elsasser number and is written as Λ???1. It follows that these modes dissipate quickly before propagating out. This restriction falls when the above criterion is no longer fulfilled (Λ?～?1), the corresponding modes evolving towards distorted QFD-modes. A systematic study of these QFD-modes is made in the limit of small Elsasser number (Λ???1), for the different symmetries allowed. Application to the Earth's and other planetary cores is then examined for an Elsasser number up to Λ?≈?O(1), in relation to the geomagnetic secular variation and the frozen-flux approximation.