Global estimates of gravity wave parameters from GPS radio occultation temperature data

Gravity waves (GWs) play critical roles in the global circulation and the temperature and constituent structures in the middle atmosphere. They also play significant roles in the dynamics and transport and mixing processes in the upper troposphere and lower stratosphere and can affect tropospheric weather. Despite significant advances in our understanding of GWS and their effects in different regions of the atmosphere in the past few decades, observational constraints on GW parameters including momentum flux and propagation direction are still sorely lacking. Global Positioning System (GPS) radio occultation (RO) technique provides global, all-weather, high vertical resolution temperature profiles in the stratosphere and troposphere. The unprecedentedly large number of combined temperature soundings from the Constellation Observing System for Meteorology, Ionosphere, and Climate and Challenging Minisatellite Payload GPS RO missions allows us to obtain GW perturbations by removing the gravest zonal modes using the wavelet method for each day. We extended the GW analysis method of Alexander et al. (2008) to three dimensions to estimate the complete set of GW parameters (including momentum flux and horizontal propagation direction) from the GW temperature perturbations thus derived. To demonstrate the effectiveness of the analysis, we showed global estimates of GW temperature amplitudes, vertical and horizontal wavelengths, intrinsic frequency, and vertical flux of horizontal momentum in the altitude range of 17.5-22.5 km during December 2006 to February 2007. Consistent with many previous studies, GW temperature amplitudes are a maximum in the tropics and are generally larger over land, likely reflecting convection and topography as main GW sources. GW vertical wavelengths are a minimum at equator, likely due to wave refraction, whereas GW horizontal wavelengths are generally longer in the tropics. Most of the waves captured in the analysis of the GPS data are low-intrinsic frequency inertia-GWs, and the estimated intrinsic frequencies scaled by the Coriolis parameter also show a strong maximum at equator. Enhanced wave fluxes are linked to convection, topography, and storm tracks, among others. As preliminary tests of the analysis in deriving horizontal propagation directions, we compared the GPS estimates with the corresponding estimates from the U. S. high vertical resolution radiosonde data using the conventional Stokes parameters method and we also conducted a separate analysis of the GPS data over the southern Andes in South America. We also showed the first global estimates of GW propagation directions from the GPS data. Finally, the sensitivity of the analysis to the temporal and spatial dimensions of the longitude x latitude x time cells and the uncertainties of the analysis and possible ways to reduce these uncertainties are discussed.