A ten year climatology of stratospheric gravity wave activity above southern Germany


Thomas Birner and Andreas Dörnbrack

DLR Oberpfaffenhofen, Institut für Physik der Atmosphäre, Wessling, Germany


FIGURES


Abstract

A climatology of stratospheric gravity wave activity is presented for the ten year period from 1990 to 1999 . The climatology is based on high vertical resolution radiosonde data of routine soundings from Munich (48° N, 12° E) and Stuttgart (49° N, 9° E), i.e. stations located close to the northern baseline of the Alps. Part of this study contributes to the cooperative 'SPARC Gravity Wave Initiative' coordinated by R.A. Vincent of the University of Adelaide, Australia.

The data are analyzed in two different stratospheric altitude ranges: a fixed one is applied for the contribution to the global climatology. A new adaptive analysis range is defined. This range takes daily variations of the tropopause height into account. Results obtained with this adaptive altitude range are contrasted with those from the fixed one. The seasonal variation of stratospheric energy density is more pronounced when analyzing with the adaptive altitude range. The use of the fixed altitude range leads to an overestimation of stratospheric energy density in summer.

In order to classify the different excitation mechanisms of gravity waves, the dependence of stratospheric wave energy density on tropospheric wind conditions is investigated. Maximum stratospheric wave energy densities are found to be connected to large tropospheric wind shears. For the investigated region the jet stream excitation seems to be the most important mechanism for the generation of stratospheric gravity waves.

Introduction

Internal gravity waves are the cause for most of the spatial and temporal fluctuations in stratospheric temperature and wind. For example, the adiabatic cooling near the crest of hydrostatic mountain waves favors the nucleation of supercooled droplets and thus the formation of polar stratospheric clouds (e.g. Wirth et al., 1999). Upward propagating gravity waves transport momentum to the middle atmosphere and are now well appreciated to play a dominant role in driving the circulation of the atmosphere (e.g. Alexander and Rosenlof., 1996). Since the present GCMs and climate models are not able to resolve the whole gravity wave spectrum, a better characterization of stratospheric gravity wave activity in terms of wave energies, periods, phase speeds, propagation directions, and spectral parameters is needed. This can be achieved by means of observational, theoretical and modeling studies.

The present work constitutes the South-German contribution to the global gravity wave climatology within SPARC, see Vincent, this issue. This climatology is based on high vertical resolution radiosonde data of about 16 groups in 13 countries, and follows the previous work of Allen and Vincent, 1995 and Vincent et al., 1997. A significant improvement for our region is reached by analyzing the data with the new adaptive altitude range. Our results help to better understand the different sources of gravity wave activity in the lower stratosphere (Nastrom and Fritts, 1992 and Fritts and Nastrom, 1992).

Measurements and Analysis Technique

For our analysis 10 years data of twice daily radiosonde launches (00 UT, 12 UT) from Munich-Oberschleißheim (MO) and Stuttgart (STU) are used. Altogether these are 14608 profiles of high vertical resolution radiosonde data. Since there are few gaps in the data (due to uncomplete or no measurements), not all of the profiles are suitable for the analysis. However, about 95% of the profiles are analyzed. The stored data of each sounding include profiles of temperature, horizontal wind speed and direction, relative humidity, and ascent rate of the balloon. In this analysis only the profiles of temperature and horizontal wind components are used. The accuracy of the temperature measurements is about 0.1 K. For the wind speed measurements it ranges from 1 to 2 ms -1 , depending on the horizontal distance to the station, since it is measured by radar tracking which leads to more inaccurate measurements for larger distances. The vertical resolution is 50  m for the temperature measurements and about 150 m for the wind speed measurements. Väisälä, RS 80 radiosondes are used and reach maximum ceilings of about 30 to 35 km in our region.

As an example, selected vertical profiles of temperature and horizontal wind speed for MO, January 1995 are shown in Figure 1. For clearness the curves are shifted from each other by 10 K and 10 ms -1 , respectively. An unperturbed background profile (T0 , u0, v0) is defined by fitting a quadratic polynomial to the measured profile in the respective altitude range (blue dashed curves). The local deviations (T', u', v' ) from this background state are considered as due to gravity waves. For this approach see e.g. Allen and Vincent, 1995; Whiteway, 1999. For our analysis two altitude ranges are applied: a fixed one as for the global climatology within SPARC (FAR = zT - zB , where zB= ZTP+ 2 km, zT = zB+ 7 km and ZTP = 11 km is the climatological zonally averaged tropopause height (Hoinka, 1998) for the investigated geographical region). The new adaptive altitude range is defined by AAR = zT - zB , where zB(t)=zTP(t) + 4 km, zT (t) = zB(t) + 11 km and zTP(t) is the time dependent tropopause height. The altitude ranges FAR and AAR are indicated in the figure by bright and dark shadings, respectively.




Figure 1: Selected vertical profiles of temperature (left) and horizontal wind speed (right) for MO, January 1995 . Subsequent curves are shifted by 10  K (left) and 10  ms -1 (right), respectively. The dashed blue lines represent the unperturbed background state and the two altitude ranges are indicated by different gray shadings.
 

The lower edge of the AAR shows very clearly the strong variability of the tropopause height in just one month. The thermal definition of the tropopause (WMO, 1957) is used to obtain its height. Within our 10 year analysis we find tropopause heights ranging from 6 to 15 km, and a mean value of about 11 km which agrees very well with the zonally averaged value of Hoinka, 1998. Thus, if one uses the FAR, there are some profiles with the tropopause lying within the analyzed range. An inversion at the tropopause (which is typical for northern midlatitudes) would then be misinterpreted as part of a wave fluctuation. Eventually, this overestimates the contribution to the calculated energy density for the respective profile. In order to avoid misinterpretation due to tropopause effects we defined the above AAR, which represents a constant range in altitude referring to the temporal varying tropopause. To ensure that there will be really no effects from the tropopause, this altitude range starts 4  km above zTP(t).

Analysis and Climatology

The interannual variability of the stratospheric total energy density is determined on the basis of the 10 year data set for MO and STU. Figure 2 shows this quantity and the ratio of kinetic to potential energy density for MO and both altitude ranges (AAR = 11 km and FAR = 7 km). The results for STU exhibit a qualitatively similar behaviour, but with smaller energy densities mainly in winter, and are not shown. The potenial and kinetic energy densities are given by Epot = 1/2(g/N)2(T'/T0)2 and Ekin = 1/2(u'2 +v'2), respectively, where we neglected the vertical part of the kinetic energy density.

Figure 2: Interannual variability of the mean stratospheric total energy density (kinetic plus potential, top) and the ratio of kinetic to potential energy density (bottom) for MO. The blue lines represent the ten year averages over each launch date, the green shaded areas indicate the range of obtained values of the ten year period, and the red lines are the monthly averaged values (all for AAR =11  km). The dashed black lines are the monthly averaged values for FAR =7  km. Note the logarithmic scales of the ordinates.

There is a strong half-yearly trend in stratospheric energy density which can be seen in Figure 2, i.e. high stratospheric energy densities in winter and lower ones in summer. Mean values for MO are 3.8  J/kg for JJA (summer), 9.2  J/kg for DJF (winter), and 6.3  J/kg for all data. For STU the mean values are 3.4  J/kg for JJA, 7.4  J/kg for DJF, and 5.4  J/kg for all data. The ratio of kinetic to potential energy density for MO (STU) is almost constant about 1.8 (1.7) during the entire year. The stratospheric energy densities derived from the FAR analysis show qualitatively the same seasonal behaviour, but they are smaller in winter and higher in summer than those derived from the AAR analysis. Power spectral analysis, which is a technique often used in the literature (e.g. VanZandt, 1982; Allen and Vincent, 1995), shows that the energy density increases with increasing vertical wavelength, Figure 3. Since larger altitude ranges can resolve longer wavelengths, the energy density should be larger for AAR =11  km than for FAR =7  km. In winter this is the case. In summer the opposite behaviour cannot be understood with the above explanation. However, this is an artifact and comes most likely from the above mentioned tropopause effect (for the FAR analysis) that occurs more often in summer, since the tropopause in summer is higher on average (Hoinka, 1998).

Figure 3: Power spectrum derived using FFT analysis of normalized temperature fluctuations (T'/T0). The dotted line indicates the power law due to saturation theory ( Smith et al., 1987). The slope t of the saturated part of the spectrum is almost constant (t = 2.9), and the mean vertical wavelength is lz =3 km (both for AAR). Note that the spectra are averaged over all data from both stations.

We now focus on the dependence of stratospheric energy density on tropospheric wind conditions, Figure 4, in order to identify the most important wave excitation mechanism north of the Alps. The tropospheric quantities are obtained between the altitudes z0+ 1 km and  zTP - 2 km, where z0 is the altitude above sea level of the individual station. In a first step we determine the energy weighted angular distribution of mean tropospheric wind direction (30° bins, insets in Figure 4). This gives the energetic significance of different tropospheric wind directions. In order to clarify the difference in stratospheric energy density for MO and STU, the energy densities are normalized by the sum of MO and STU. Obviously, the tropospheric wind directions in the westerly sector show high positive correlation with the stratospheric energy density. The energy weighted mean tropospheric wind directions derived from this are 287° for MO and 282° for STU. These are compatible with the wind directions from which in our region most of the active weather systems as fronts come from.

In a second step we investigate whether critical level filtering inhibits the vertical propagation of the gravity waves and thus reduces the stratospheric wave energy density. To see that, we consider the energy weighted distribution of the tropospheric directional shear, Figure 4. Again, the energy densities are normalized by the sum of MO and STU. This distribution indicates that about 80 % of the normalized energy densities are correlated with directional shears between 0° and 90° . Thus high stratospheric energy density is associated with small tropospheric directional shear in agreement with the results in Whiteway, 1999.

Figure 4:Energy weighted distribution of the tropospheric directional shear. Insets: energy weighted angular distributions of mean tropospheric wind directions. The dotted lines indicate the energy weighted mean tropospheric wind directions.

Finally, in a third step, we consider mean relative statistical correlations between stratospheric energy density and magnitudes derived from the tropospheric wind speed, Figure 5. They show a strong contrast between summer (correlations between -0.1 and 0.1 ) and winter (correlations between 0.2 and 0.6 ). For the winter we find positive correlation of stratospheric energy density with the surface wind speed, larger correlation with the tropospheric wind shear, and even larger correlation with both, that means the product of surface wind speed and tropospheric wind shear. Thus, since surface wind speed and tropospheric wind shear are measures for the strength of the jetstream, there is a high positive corralation in our region between jetstream strength and observed stratospheric energy density. We remark that there is a significant difference between MO and STU only for January and February referring to the correlation between stratospheric energy density and tropospheric wind shear.

Figure 5:Interannual variability of the mean relative statistical correlations . Full lines are for MO, dashed lines for STU.

Summary and Conclusions

This study constitutes the South-German contribution to the global gravity wave climatology within SPARC (see Vincent, this issue). The results derived with the newly developed adaptive analysis height range were critically compared with those of the fixed altitude range used for the global climatology. Both methods give high stratospheric energy densities occuring mostly in winter, associated with large tropospheric wind shear, small tropospheric directional shear and large surface wind speed. Lower stratospheric energies are found in summer. For the investigated region the use of the AAR analysis shows a stronger half-yearly trend in stratospheric energy density than the FAR analysis, and makes thus seasons better comparable. It avoids the overestimation of stratospheric energy density of the FAR analysis in summer. Furthermore, higher energy densities are found for MO compared to STU. This could be due to the fact that Munich lies closer to the Alps than Stuttgart and some of the MO-profiles are thereby certainly influenced by the effect of these mountains. The dependence of stratospheric energy density on tropospheric wind conditions support the hypothesis that the main source of stratospheric gravity wave activity in the investigated region is the jet stream excitation. Up to now the relative contribution of convective or orographic excitation is not clear for the investigated region. Since the Alps are very complex, 90° propagating mountain waves are possibly a part of the sources for the observed stratospheric energy densities.

For further investigations idealized studies and realistic model studies are needed. Together with this observational study they should lead to a reliable identification of the different gravity wave generation mechanisms.

Acknowledgments

We thank the DWD Oberschleißheim for the free availability of the radiosonde data. We furthermore thank Jana Freund and Winfried Beer for their help in solving the difficulties with different data formats. We are particularly grateful to Ulrich Hamann for his help in analyzing the huge amount of data.

References

Alexander, M.J. and K.H. Rosenlof, Nonstationary gravity wave forcing of the stratospheric zonal mean wind, J. Geophys. Res., 101, 23,465-23,474, 1996.

Allen, S.J. and R.A. Vincent, Gravity wave activity in the lower atmosphere: Seasonal and latitudinal variations, J. Geophys. Res., 100, 1327-1350, 1995.

Fritts, D.C. and G.D. Nastrom, Sources of mesoscale variability of gravity waves. Part II: frontal, convective, and jet stream excitation, J. Atmos. Sci., 49, 111-127, 1992.

Hoinka, K.P., Statistics of the Global Tropopause Pressure, Mon. Wea. Rev., 126, 3303-3325, 1998.

Nastrom, G.D. and D.C. Fritts, Sources of mesoscale variability of gravity waves. Part I: topographic excitation, J. Atmos. Sci., 49, 101-110, 1992.

Smith, S.A., D.C. Fritts, and T.E. VanZandt, Evidence of a saturation spectrum of atmospheric waves, J. Atmos. Sci., 44, 1404-1410, 1987.

VanZandt, T.E., A universal spectrum of buoyancy waves in the atmosphere, Geophys. Res. Lett., 9, 575-578, 1982.

Vincent, R.A., S.J. Allen, and S.D. Eckermann, Gravity wave parameters in the lower stratosphere, in Gravity wave processes: Their parametrization in global climate models, edited by K. Hamilton, NATO ASI ser. I, 50, 7-25, 1997.

Vincent, R.A., Outcomes of the SPARC gravity wave radiosonde initiative, this issue

Whiteway, J.A., Enhanced and inhibited gravity wave spectra, J. Atmos. Sci., 56, 1344-1352, 1999.

Wirth, M., A. Tsias, A. Dörnbrack, V. Weiß, K.S. Carslaw, M. Leutbecher, W. Renger, H. Volkert, and T. Peter, Model-guided Lagrangian observation and simulation of mountain polar stratospheric clouds, J. Geophys. Res., 104, 23,971-23,981, 1999.

WMO, Definition of the tropopause, WMO Bull., 6, 136, 1957.


Back to

Session 1 : Stratospheric Processes and their Role in Climate Session 2 : Stratospheric Indicators of Climate Change
Session 3 : Modelling and Diagnosis of Stratospheric Effects on Climate Session 4 : UV Observations and Modelling
AuthorData
Home Page