Global Estimates of Gravity Wave Parameters in the UTLS from GPS Temperature Retrievals

Global Estimates of Gravity Wave Parameters in the UTLS from GPS Temperature Retrievals

Global Estimates of Gravity Wave Parameters in the UTLS from GPS Temperature Retrievals

Ling Wang and M. Joan Alexander

NorthWest Research Associates, CoRA Div.

Ling Wang and M. Joan Alexander

NorthWest Research Associates, CoRA Div.

Overview of TalkOverview of Talk

Introduction COSMIC and CHAMP GPS temperature

data New method to derive gravity wave

parameters and some preliminary results Uncertainties of the new GW analysis Summary

Introduction COSMIC and CHAMP GPS temperature

data New method to derive gravity wave

parameters and some preliminary results Uncertainties of the new GW analysis Summary

2/192/19

Introduction (1)Introduction (1) GWs play significant roles in the dynamics and

transport and mixing processes in the UTLS and they can also affect the formation of cirrus clouds.

At present, observational constraints on GW parameters including momentum flux & propagation direction are sorely needed to improve the modeling of GW effects on various processes and phenomena in the UTLS (and elsewhere).

GWs play significant roles in the dynamics and transport and mixing processes in the UTLS and they can also affect the formation of cirrus clouds.

At present, observational constraints on GW parameters including momentum flux & propagation direction are sorely needed to improve the modeling of GW effects on various processes and phenomena in the UTLS (and elsewhere).

3/193/19

Introduction (2)Introduction (2) Hodograph analysis and Stokes parameters

method (e.g., Eckermann and Vincent, 1989, Pure Appl. Geophys) have been used to derive GW propagation direction and other wave parameters.

These methods require, however, both temperature & wind information, which is NOT available for most satellite data.

Hodograph analysis and Stokes parameters method (e.g., Eckermann and Vincent, 1989, Pure Appl. Geophys) have been used to derive GW propagation direction and other wave parameters.

These methods require, however, both temperature & wind information, which is NOT available for most satellite data.

4/194/19

Introduction (3)Introduction (3) Ern et al. (2004, JGR) use phase differences between adjacent

T profiles to get global estimates of GW horizontal wavelength and momentum flux from the CRISTA satellite data.

Ern et al. (2004, JGR) use phase differences between adjacent T profiles to get global estimates of GW horizontal wavelength and momentum flux from the CRISTA satellite data.

5/195/19Horizontal wavelength and momentum flux at 25 km altitude in August 2007 (Ern et al., 2004).Horizontal wavelength and momentum flux at 25 km altitude in August 2007 (Ern et al., 2004).

Introduction (4)Introduction (4)

Alexander et al. (2008, JGR) use a similar approach to derive GW horizontal structures from the HIRDLS satellite data.

Alexander et al. (2008, JGR) use a similar approach to derive GW horizontal structures from the HIRDLS satellite data.

6/196/19

Momentum flux and horizontal wavenumber between 20-30 km altitude in May 2006 (Alexander et al., 2008).

Momentum flux and horizontal wavenumber between 20-30 km altitude in May 2006 (Alexander et al., 2008).

Introduction (5)Introduction (5)

In most GW analyses of satellite data so far, GW propagation directions are not routinely derived.

This study introduces a new method to estimate the complete set of GW parameters (including horizontal propagation direction) from the GPS RO temperature data alone.

In most GW analyses of satellite data so far, GW propagation directions are not routinely derived.

This study introduces a new method to estimate the complete set of GW parameters (including horizontal propagation direction) from the GPS RO temperature data alone.

7/197/19

Data (1)Data (1)

(from Anthes et al., 2008, BAMS)(from Anthes et al., 2008, BAMS)

GPS RO technique: phase delay of GPS signals --> bending angle of radio wave --> atmospheric refractivity --> temperature profiles

GPS RO technique: phase delay of GPS signals --> bending angle of radio wave --> atmospheric refractivity --> temperature profiles

8/198/19

COSMIC (Constellation Observing System for Meteorology Ionosphere and Climate/Formosa Satellite 3) GPS RO data

CHAMP (Challenging Minisatellite Payload) GPS RO data

COSMIC (Constellation Observing System for Meteorology Ionosphere and Climate/Formosa Satellite 3) GPS RO data

CHAMP (Challenging Minisatellite Payload) GPS RO data

Data (2)Data (2)

Map showing a typical daily COSMIC/CHAMP data coverage.

Unlike most other satellite data, GPS data do not have regular orbits.

COSMIC: 2006-present; ~ 1500 daily profiles

CHAMP: 2001-2008; ~ 150 daily profiles

We use dry T profiles processed by UCAR CDAAC.

Vertical resolution: ~ 1 km Accuracy: sub-Kelvin Temperature data cover

troposphere and stratosphere.

Unlike most other satellite data, GPS data do not have regular orbits.

COSMIC: 2006-present; ~ 1500 daily profiles

CHAMP: 2001-2008; ~ 150 daily profiles

We use dry T profiles processed by UCAR CDAAC.

Vertical resolution: ~ 1 km Accuracy: sub-Kelvin Temperature data cover

troposphere and stratosphere.

9/199/19

Gravity Wave Analysis (1)Gravity Wave Analysis (1)GW Temperature PerturbationGW Temperature Perturbation

One example of GPS temperature, background T determined from One example of GPS temperature, background T determined from our procedure and the resulting GW temperature perturbation profile.our procedure and the resulting GW temperature perturbation profile.

Since large-scale wave motions (e.g., Kelvin waves) can have similar vertical wavelength as GWs, we extract GW perturbations by removing zonal modes 0-6.

Since large-scale wave motions (e.g., Kelvin waves) can have similar vertical wavelength as GWs, we extract GW perturbations by removing zonal modes 0-6.

10/1910/19

We derive the dominant GW T amp & vertical wavelength at each height from wavelet analysis.

GW amplitudes generally decrease with increasing latitudes.

Vertical wavelengths also display strong latitudinal variability.

We derive the dominant GW T amp & vertical wavelength at each height from wavelet analysis.

GW amplitudes generally decrease with increasing latitudes.

Vertical wavelengths also display strong latitudinal variability.

Gravity Wave Analysis (2)T Amplitude &

Vertical WavelengthT Amplitude &

Vertical Wavelength

11/1911/19Dominant T amp & vertical wavelength averaged between 17.5-22.5 km during December 2006-February 2007.

Gravity Wave Analysis (3)Gravity Wave Analysis (3)Horizontal Wave VectorHorizontal Wave Vector

If we neglect time variation, at the same z & for the same vertical wavelength, phase differences among adjacent soundings are related with horizontal wave vector (k, l) and distances among the soundings as follows:

We can determine phase differences from cross-wavelet analysis. If there are at least 3 soundings in one space-time cell (e.g., 15 deg by 15 deg by 6 hours), we can solve for an over-determined linear problem for (k, l) using LSF.

n soundings --> [n!/(n-2)!/2] equations, e.g., 4 soundings --> 6 equations to solve for (k, l)

12/1912/19

Gravity Wave Analysis (4)Gravity Wave Analysis (4)Horizontal WavelengthHorizontal Wavelength

We use (15 deg X 15 deg X 6 hours) cells across the globe.

GWs are generally longer at lower latitudes, being consistent with many previous studies.

13/1913/19averaged between 17.5-22.5 km during December 2006-February 2007

Gravity Wave Analysis (5)Gravity Wave Analysis (5)

Most of the waves captured in our analysis are low intrinsic frequency inertia-GWs.

The latitudinal variation of intrinsic frequency is consistent with radiosonde analysis (Wang, Geller, and Alexander, 2005, JAS).

<-- GW dispersion relation

14/1914/19averaged between 17.5-22.5 km during December 2006-February 2007

Gravity Wave Analysis (6)Gravity Wave Analysis (6)With T’, k, l, m, & intrinsic frequency derived, the momentum flux components are evaluated using the left formula.

The magnitudes of MF are reasonable comparing with previous studies.

The effects of both convection and topography on GW excitations are clearly seen.

15/1915/19averaged between 17.5-22.5 km during December 2006-February 2007

Gravity Wave Analysis (7)Gravity Wave Analysis (7)

We compare GPS results with radiosonde analysis using the Stokes parameters method.

GWs generally propagate eastward in the tropics and westward in mid- and high-latitudes, being opposite to the prevailing background winds.

The two analyses are not dealing with exactly the same part of GW spectrum due to the details of the two analyses.

16/1916/1918-24.9 km, December 2006-February 2007

Gravity Wave Analysis (8)Gravity Wave Analysis (8)Topography is one of the major GW sources and mountain waves are expected to propagate upwind and orthogonal to the orientation of topography.

Results largely consistent with what is expected from the orientation of the southern Andes and the prevailing westerlies

17/1917/19

17.5-22.5 km,June-August 2007

Uncertainties of AnalysisUncertainties of Analysis We assume the existence of a single & dominant

GW mode in each (lon, lat, time) cell at each height

- Can reduce uncertainties by reducing the size of each cell and setting amplitude thresholds for dominant modes

- & by including treatment for secondary modes (if there are any)

We neglect time variation of GW modes - Can solve this by using the following complete formula to solve

for k, l, &k, l, & ground-based frequency simultaneously (which puts additional requirement for data coverage though):

We assume the existence of a single & dominant GW mode in each (lon, lat, time) cell at each height

- Can reduce uncertainties by reducing the size of each cell and setting amplitude thresholds for dominant modes

- & by including treatment for secondary modes (if there are any)

We neglect time variation of GW modes - Can solve this by using the following complete formula to solve

for k, l, &k, l, & ground-based frequency simultaneously (which puts additional requirement for data coverage though):

18/1918/19

SummarySummary We introduce a new method based on the cross-

wavelet analysis to estimate the complete set of GW parameters from the GPS RO temperature data.

We present some preliminary results in the UTLS which demonstrate the effectiveness of the new method.

The uncertainties of this new GW analysis and possible ways to deal with them are discussed.

We introduce a new method based on the cross-wavelet analysis to estimate the complete set of GW parameters from the GPS RO temperature data.

We present some preliminary results in the UTLS which demonstrate the effectiveness of the new method.

The uncertainties of this new GW analysis and possible ways to deal with them are discussed.

19/1919/19