New insights into Rossby wave packet properties in the extratropical UTLS using GNSS radio occultations

The present study describes Rossby wave packet (RWP) properties in the upper-troposphere and lower-stratosphere (UTLS) with the use of Global Navigation Satellite System radio occultation (GNSS-RO) measurements. This global study covering both hemisphere’s extratropics is the first to tackle medium and synoptic-scale waves with GNSS-RO. We use one decade of GNSS-RO temperature and pressure data from the CHAMP, COSMIC, GRACE, Metop-A, Metop-B, SAC-C and TerraSAR-X 5 missions; combining them into one gridded dataset for the years 2007-2016. Our approach to extract RWP anomalies and their envelope uses Fourier and Hilbert transforms over longitude without preor post-processing the data. Our study is purely based on observations, only using ERA-Interim winds to provide information about the background wind regimes. The RWP structures that we obtain in the UTLS agree well with theory and earlier studies, in terms of coherent phase/group propagation, zonal scale and distribution over latitudes. Furthermore, we show that RWP pressure anomalies maximize around 10 the tropopause, while RWP temperature anomalies maximize right above tropopause height with a contrasting minimum right below. RWP activity follows the zonal-mean tropopause during all seasons. RWP anomalies in the lower stratosphere are dynamically coupled to the upper troposphere. They are part of the same system with a quasi-barotropic structure across the UTLS. RWP activity often reaches up to 20 km height and occasionally higher, defying the Charney-Drazin criterion. We note enhanced amplitude and upward propagation of RWP activity during 15 sudden stratospheric warmings. We provide observational support for improvements in RWP diagnostics and wave trend analysis in models and reanalyses. Wave quantities follow the tropopause, and diagnosing them on fixed pressure levels (which the tropopause does not follow) can lead to aliasing. Our novel approach analysing GNSS-RO pressure anomalies provides wave signals with better continuity and coherence across the UTLS and the stratosphere, compared to temperature anomalies. Thus, RWP vertical propagation is 20 much easier to analyse with pressure data. 1 https://doi.org/10.5194/acp-2020-124 Preprint. Discussion started: 4 March 2020 c © Author(s) 2020. CC BY 4.0 License.

Abstract. The present study describes Rossby wave packet (RWP) properties in the upper troposphere and lower stratosphere (UTLS) with the use of Global Navigation Satellite System radio occultation (GNSS-RO) measurements. This global study covering both hemispheres' extratropics is the first to tackle medium-and synoptic-scale waves with GNSS-RO. We use 1 decade of GNSS-RO temperature and pressure data from the CHAMP, COSMIC, GRACE, Metop-A, Metop-B, SAC-C and TerraSAR-X missions, combining them into one gridded dataset for the years 2007-2016. Our approach to extract RWP anomalies and their envelope uses Fourier and Hilbert transforms over longitude without pre-or post-processing the data. Our study is purely based on observations, only using ERA-Interim winds to provide information about the background wind regimes.
The RWP structures that we obtain in the UTLS agree well with theory and earlier studies, in terms of coherent phase or group propagation, zonal scale and distribution over latitudes. Furthermore, we show that RWP pressure anomalies maximize around the tropopause, while RWP temperature anomalies maximize right above the tropopause height with a contrasting minimum right below. RWP activity follows the zonal-mean tropopause during all seasons.
RWP anomalies in the lower stratosphere are dynamically coupled to the upper troposphere. They are part of the same system with a quasi-barotropic structure across the UTLS. RWP activity often reaches up to 20 km height and occasionally higher, defying the Charney-Drazin criterion. We note enhanced amplitude and upward propagation of RWP activity during sudden stratospheric warmings.
We provide observational support for improvements in RWP diagnostics and wave trend analysis in models and re-analyses. Wave quantities follow the tropopause, and diagnosing them on fixed pressure levels (which the tropopause does not follow) can lead to aliasing. Our novel approach analyzing GNSS-RO pressure anomalies provides wave signals with better continuity and coherence across the UTLS and the stratosphere, compared to temperature anomalies. Thus, RWP vertical propagation is much easier to analyze with pressure data. therefore an intermediate range of w4-15 or similar is generally used to study RWPs (e.g., Zimin et al., 2003;Glatt and Wirth, 2014;Wolf and Wirth, 2017;Quinting and Vitart, 2019). The extent in longitude of RWPs, or the size of the envelope that modulates the wave's amplitude, also shows high variability: in the Southern Hemisphere (SH) a more global scale is observed, while in the Northern Hemisphere (NH) the waves tend to be more localized, due to the extent and zonality of the waveguide (e.g., Lau, 1979;Randel and Stanford, 1985). Nevertheless, either behavior can occur in both hemispheres at specific times, and even when RWPs reach a global scale, the amplitude of the wave is rarely constant along all longitudes. Global or local wave modes show no apparent differences in their dynamics (Randel and Stanford, 1985). For consistency we will use the term RWP throughout the paper to refer to wave activity, since it simultaneously covers the medium and synoptic scales.
RWPs are the main driver of extratropical weather and climate, also at larger spatial and timescales. Variations in the zonal-mean westerlies of the order of 2-3 weeks were found to originate from energy conversion between zonalmean flow and medium-scale waves in the SH (Webster and Keller, 1975;Randel and Stanford, 1985). Eddy feedbacks are responsible for the increased timescale of annular modes in the extratropical winter troposphere which are coupled with the stratospheric polar vortex behavior Hartmann, 2001, 2003), making RWPs an important contributor to stratosphere-troposphere coupling (Kidston et al., 2015).
The theory for lower-upper level PV anomaly coamplification (Hoskins et al., 1985), idealized life cycles (Gall, 1976;Simmons and Hoskins, 1978;Thorncroft et al., 1993), and conceptual and case studies (Shapiro and Keyser, 1990) have all highlighted the importance of tropopause processes in the formation and evolution of baroclinic waves for decades, but research has only recently started to focus on the role of the stratosphere. Williams and Colucci (2010) showed that RWP properties in the upper troposphere depended on lower-stratospheric conditions, namely the strength of the polar vortex. New case studies about extreme cyclones are pointing to a significant role of stratospheric conditions in their development (Odell et al., 2013;Tao et al., 2017a, b), and a more general study about severe European windstorms found that in 20 of 60 cases the stratospheric contribution during their deepening phases was over 10 % (Pirret et al., 2017). Nevertheless, despite their interaction with the stratosphere, it is generally assumed that RWPs cannot propagate upward due to the typical wind regimes in the stratosphere (Charney and Drazin, 1961).
There are several issues concerning RWP representation in forecast models. Atmospheric models present a general bias in their simulated waveguides, which are not as sharp as observed (Gray et al., 2014;Giannakaki and Martius, 2016). The horizontal and vertical structures of forecast errors tend to spread together with the RWP's envelope (Dirren et al., 2003;Hakim, 2005;Sellwood et al., 2008;Zheng et al., 2013), and near-tropopause dynamics in particular are the ones dominating error growth (Baumgart et al., 2018). Apart from baroclinic energy generation, other important processes that affect PV distribution in RWPs near the tropopause (and thereby forecast error growth) are latent-heat release and longwave radiative cooling (e.g., Martínez-Alvarado et al., 2016;Teubler and Riemer, 2016). The organized structure of RWPs and their large traveling distances, compared to individual troughs and ridges or cyclones and anticyclones, offers the opportunity of improving and extending the range of weather forecasts with a better understanding and modeling of their dynamics (Grazzini and Vitart, 2015).
Climate change can have opposing influences on the storm tracks (Shaw et al., 2016) or the polar front jet (Hall et al., 2015) position. As an example, Arctic amplification will decrease the lower troposphere's meridional temperature gradient and baroclinicity, tending to shift the storm track equatorward, while stratospheric radiative cooling by CO 2 will increase near-tropopause baroclinicity towards the pole, counterbalancing the shift in the lower troposphere. This is one reason why climate projections show sensitivity to the inclusion of the stratosphere in models (Scaife et al., 2012), and RWP dynamics play a central role in establishing where the projected storm tracks and polar front jets will shift to. Models used for climate projections, e.g., the Coupled Model Intercomparison Project Phase 5 (CMIP5), are known to have biases in their storm tracks (Chang et al., 2012;Hall et al., 2015).
Recognizing the importance of RWPs for weather forecasting, climate change projections and stratospheretroposphere interactions, the scientific community would greatly benefit from increased observational knowledge of RWP structures and behavior across the upper troposphere and lower stratosphere (UTLS), which is the main goal of our study in an area lacking research thus far. The recent availability of Global Navigation Satellite System radio occultation (GNSS-RO) measurements (e.g., Kursinski et al., 1997;Anthes, 2011) enables the study of RWP temperature and pressure anomalies with global coverage across the UTLS at unprecedented high vertical resolution. The exclusive use of 1 decade of GNSS-RO observations (2007 will avoid any limitation from reanalysis or model output in the form of biases or lack of vertical resolution within the UTLS. The present paper intends to describe general and purely observational RWP properties across the UTLS. GNSS-RO observations have been used to study atmospheric waves in many recent studies. In the equatorial regions, Kelvin waves have been studied the most so far (e.g., Randel and Wu, 2005;Flannaghan and Fueglistaler, 2013;Scherllin-Pirscher et al., 2017), with Zeng et al. (2012) and Kim and Son (2012) also studying the Madden-Julian oscillation (MJO; Madden and Julian, 1994). Alexander et al. (2008) andPilch Kedzierski et al. (2016a) had a wider approach and considered all equatorial wave types. In the ex-tratropics, studies using GNSS-RO measurements have concentrated mainly on the extraction of gravity wave parameters (e.g., Tsuda et al., 2000;Wang and Alexander, 2010;Kohma and Sato, 2011;Tsuda, 2014;Schmidt et al., 2016). Madhavi et al. (2015) studied the 2 d wave in the upper stratosphere and lower mesosphere in both the NH and SH extratropics. Alexander and Shepherd (2010) studied planetary wave activity in the Arctic and Antarctic lower and mid-stratosphere, while Shepherd and Tsuda (2008) studied planetary waves in the SH polar summer at 30 km height. Pilch Kedzierski et al. (2017) used GNSS-RO to study extratropical tropopause modulation by waves, although only separating wave activity by timescale and propagation direction, without describing wave properties. Our study is a first attempt to describe RWP properties in the extratropics with the use of GNSS-RO.
This paper is organized as follows: Sect. 2 will introduce the GNSS-RO missions used for our study and the analysis methods applied, Sect. 3 will present how the zonally averaged RWP activity is distributed in the UTLS globally, and Sect. 4 will concentrate on RWP zonal and vertical structures. Both Sects. 3 and 4 will start with case studies to introduce how RWPs evolve across the UTLS, moving on to climatological statistics with a focus on general and common RWP properties in the extratropics of both hemispheres. Section 5 will discuss some implications of our results, and Sect. 6 will summarize our main findings.

Datasets
We make use of GNSS-RO measurements from the following satellite missions: CHAMP (Wickert et al., 2001), COSMIC (Anthes et al., 2008), GRACE (Beyerle et al., 2005), Metop-A (von Engeln et al., 2011), the successive Metop-B, SAC-C (Hajj et al., 2004) and TerraSAR-X (Beyerle et al., 2011). All data are reprocessed or post-processed occultation profiles with moisture information ("wetPrf" product) from the COSMIC Data Analysis and Archive Center (CDAAC, https: //cdaac-www.cosmic.ucar.edu/cdaac/products.html, last access: 15 July 2020), for the years 2007-2016. GNSS-RO profiles of temperature and pressure are provided interpolated on a 100 m vertical grid between the surface and 40 km height, although their effective physical resolution varies from ∼ 1 km in regions of constant stratification down to the order of 100 m where stratification gradients occur, such as the tropopause or the top of the boundary layer (Kursinski et al., 1997;Gorbunov et al., 2004).
The wetPrf product is the best suited for our study for several reasons. In regions where the atmosphere is dry (roughly above 10 km height), the wetPrf and the dry temperature "atmPrf" profiles basically coincide (Alexander et al., 2014;Danzer et al., 2014). At lower levels, water vapor is increas-ingly influential and dry temperatures get colder than the real temperature (e.g., at 8 km height in the extratropics, the difference is already of the order of 0.5-2 K and increasing downwards in Danzer et al., 2014). Our wave analyses start at 6 km height, and we note that the extratropical tropopause can vary between 7 and 12 km height over one latitude band on the same day (see Sect. 4.1 and 4.2). Water vapor content within the RWP's troughs and ridges can be very different at a given vertical level and its effect on the GNSS-RO profiles, especially within the lower part of the UTLS, cannot be neglected if reliable wave amplitudes are desired. To account for water vapor effects on the retrieved GNSS-RO profile, a background state (typically analyses) is assimilated using one-dimensional variational (1D-Var) technique (Healy and Eyre, 2000;Poli et al., 2002). Although this procedure combines the GNSS-RO measurement with the "firstguess" background, it has been shown that between 9 and 22 km height the GNSS-RO measurement dominates the retrieved vertical profile (Healy and Eyre, 2000), whereas at lower levels there are temperature improvements from the background towards radiosonde measurements where available (Poli et al., 2002).
Different studies have shown the consistency, mission independence and good precision among different GNSS-RO satellite missions as well as compared to radiosondes (Hajj et al., 2004;Wickert et al., 2009;Schreiner et al., 2011;Anthes et al., 2008;Ho et al., 2017). The advantage of GNSS-RO profiles over radiosondes relies on their global coverage, weather independence and higher sampling density. The time period of our analysis spans 2007-2016, and thanks to the simultaneous use of several missions, the total number of GNSS-RO profiles is fairly stable with around 2500 profiles per day.
All GNSS-RO profiles undergo an in-depth quality check (see Sect. 2.2) to eliminate erroneous or unphysical profiles, which, if present often enough, can trigger false signals when applying space-time filters. After merging all GNSS-RO data into one gridded dataset (see Sect. 2.3), we find no discontinuities or artifacts arising from the presence of GNSS-RO from different missions, due to the self-calibration and consistency of each satellite instrument and the post(re)processing at the same center (CDAAC). The RWP temperature anomalies we find after filtering our gridded dataset (see Sect. 2.4) are 1-2 orders of magnitude higher than the spread of large-scale and long-term temperature differences across different GNSS-RO processing centers, depending on the height range (Steiner et al., 2020).
The GNSS-RO measurements are complemented with zonal-mean zonal wind data from the ERA-Interim reanalysis  for the same period 2007-2016, in order to provide a context about background wind regimes. The programming language "R" (R Core Team, 2015) is used to perform the quality control (Sect. 2.2), gridding (Sect. 2.3) and wave analysis (Sect. 2.4) from the GNSS-RO data.

Quality control for GNSS-RO profiles
The quality control for the GNSS-RO measurements from all satellite missions is performed in two steps.
1. The first step is intended as a general screening for profiles whose temperatures or tropopause stability fall outside global climatological values. Temperature profiles with values < −100 or > 50 • C are excluded, as are those with T < −90 • C above 35 km height, since the coldest values within the winter polar vortex occur well below 35 km. The tropopause height (TP z ) is defined following the World Meteorological Organization lapse-rate criterion (WMO, 1957). Profiles where the tropopause cannot be found, or those with tropopauses that are unreasonably warm (> −45 • C), are excluded. This is close to the mean tropopause temperature in polar summer (−50 • C), but we do not find discontinuities in the time availability of GNSS-RO profiles arising from this criterion. Static stability is calculated as the Brunt-Väisälä frequency squared, N 2 = (g/ ) · (∂ /∂z), where g is the gravitational acceleration, the potential temperature and ∂z its vertical derivative. The N 2 maximum above TP z (up to 3 km) is computed, and profiles with tropospheric values (N 2 < 3 × 10 −4 s −2 ) or values that are too high (N 2 > 100 × 10 −4 s −2 ) are excluded. These conditions are based on Son et al. (2011) for TP z temperature climatologies from GNSS-RO measurements and Pilch Kedzierski et al. (2016a) for maximum N 2 climatologies, whose highest values are found in the tropics. From the initial 10 053 153 profiles available for 2007-2016 from all GNSS-RO satellite missions, this first step keeps 9 369 092 (93.2 %) of them.
2. After step 1 there remains a significant number of GNSS-RO profiles whose stratospheric temperatures markedly stand out when compared to nearby occultations in space and time, although passing the climatological criteria. On a global 5 • by 5 • longitude-latitude grid, for each day throughout 2007-2016 and at every grid point, GNSS-RO profiles within ±3 d, ±10 • longitude and ±5 • latitude are selected. A distribution of their 30 km temperature is computed, and a mean temperature profile is calculated from those GNSS-RO profiles that fall between the 0.2 and 0.8 quantiles, which avoids the influence of possibly erroneous profiles. The integrated squared temperature difference from the mean profile between 20 and 40 km height is calculated for each occultation: A i = 40 km 20 km (T i ) 2 . The 0.5 quantile of A represents the half of the selected GNSS-RO profiles that are closer to the mean profile. The profiles with A i exceeding 20 times the 0.5 quantile (that fall far out of the distribution) are excluded.
Step 2 eliminates GNSS-RO profiles with unrealistic stratospheric temperatures for their season and location, and even in extreme situations with stark temperature contrasts such as sudden stratospheric warmings this criterion is not met.
After the quality control is carried out, 9 215 804 profiles (91.7 % of the initial GNSS-RO dataset) are used.

Gridding
The method to grid GNSS-RO profiles in our study is a refinement of Pilch Kedzierski et al. (2017) and Pilch Kedzierski et al. (2016a), with better horizontal resolution and with ground-based instead of tropopause-based averaging. Both were developed after the method by Randel and Wu (2005). GNSS-RO profiles of temperature and pressure for the years 2007-2016 are gridded daily on a 5 • by 5 • longitude-latitude grid, from 85 • S to 85 • N. The height range analyzed is between 6 and 40 km, although because the focus of this paper is the UTLS region, most of the presented results will use a lower lid. The amount of GNSS-RO profiles that penetrate deeper than 6 km diminishes at lower altitudes. Having fewer data available to grid forces more interpolation and less reliability of filtered signals at lower levels. The GNSS-RO profiles that fall within the grid point area are averaged as follows: where λ is longitude, φ is latitude, z is height and t is time. The weight w i is a Gaussian function determined by each GNSS-RO profile's distance from the grid center: with D x and D y as the distances in degrees longitude and degrees latitude from the grid's center and D t as the time distance in hours from 12:00 UTC, divided by the grid's half size in all dimensions: 2.5 • longitude, 2.5 • latitude and 12 h. Typically 2-3 GNSS-RO profiles of the same day are selected within the grid's area for averaging, with the following exceptions. In a 28 % of grid points the allowed distance to search for GNSS-RO profiles needs to be expanded to ±5 • longitude and ±3.5 • latitude in order to avoid gaps (in 9 % of cases it expands further to ±10 • longitude and ±5 • latitude), without changing the weighting function w i in any case. The remaining 1 % of gaps are filled by averaging adjacent grid points (±1 in longitude, then ±1 in time). The exception's percentages presented here belong to the Northern Hemisphere from 15 • to 85 • latitude and are nearly identical in the Southern Hemisphere. The gridded tropopause height TP z (λ, φ, t) is calculated with the same weighting of all profiles' tropopauses.
Overall the real resolution of the gridded dataset is slightly coarser than 5 • longitude-latitude and may be viewed as an interpolation to a certain degree. Nonetheless this setting resolves the horizontal scale of the RWPs analyzed in our study very well (see Sects. 3 and 4).

Wave analysis
Given that the main goal of this paper is to describe properties of RWPs in the UTLS in a very general manner, we intend to keep the wave analysis as simple as possible. In most analyses a fast Fourier transform (FFT) in longitude is used to extract the wave anomalies from the gridded GNSS-RO temperature and pressure data, either for individual wavenumbers (Sect. 3) or using an intermediate range typical of RWPs (Sect. 4). For RWP envelope reconstruction, the Hilbert transform (Zimin et al., 2003) is applied to RWP pressure anomalies without any further pre-or postprocessing of the RWP anomalies or their envelope (Sect. 4). A couple of exceptions that add a degree of complexity to the analysis are disclosed below.
FFT in longitude and time. To avoid including stationary waves when showing the climatological distribution of RWP activity over latitudes, we use a two-dimensional FFT in longitude and time (Schreck, 2009)  Vertical scale analysis. RWP anomalies are extracted by taking zonal wavenumbers 4-8 with an FFT in the longitude dimension, obtaining daily longitude-height snapshots. In a second step, at every longitude grid an additional FFT is performed in the vertical dimension on the profile of RWP anomalies between 6 and 36 km height, in order to obtain the power spectrum of the different vertical wavelengths. The power spectra are then averaged for midlatitudes in Fig. 10 in Sect. 4.2.
Zonality of RWP envelopes. An FFT in longitude is used to obtain RWP anomalies (w4-8), on which the Hilbert transform is applied to reconstruct the RWP envelope. In a third step, another FFT in longitude is applied to the reconstructed envelope to obtain its zonal wavenumber power spectrum. The average wavenumber spectra of RWP envelopes at midlatitudes are shown in Fig. 12 in Sect. 4.2.
Equivalence of RWP pressure anomalies to geostrophic meridional winds. The RWP pressure anomalies in our study are defined relative to the zonal-mean pressure at each level: let us name it P rel = P /P . This notation is very convenient, since P rel is proportional to geostrophic meridional winds as explained next.
Under geostrophic balance, f v g = 1 ρ · δP δx , with f the Coriolis parameter, v g geostrophic meridional wind, ρ air density, P pressure and δ/δx its partial derivative in the longitude dimension. Knowing ρ = P /RT where R is the specific gas constant for air and T absolute temperature, Note that both T and P in Eq. (4) are in absolute terms, with zonal variations being relatively small and therefore influencing v g much less than the δP δx term. It follows that wave amplitudes of v g and P rel should remain proportional over a given height range. Our results throughout Sects. 3 and 4 indicate that P rel decreases away from the tropopause, which agrees with v g resulting from jet stream undulations whose wind maxima are located nearby.

Distribution of wave activity in the extratropical UTLS
Section 3 concentrates on the analysis of wave activity as the amplitude of individual harmonics, which only allows us to study RWPs from a zonally integrated perspective. Section 3.1 will introduce examples of wave activity behavior from gridded temperature and pressure GNSS-RO data for 1 year and one latitude band: July 2008 till June 2009 at 50 • N. Section 3.2 will generalize those results with climatological statistics for the 2007-2016 period throughout the whole extratropics. Results are presented this way so the reader gets an overall impression of where wave activity is located as well as how it evolves on a day-to-day basis. Section 4 will analyze RWP zonal structures and their horizontal propagation.

Time-height section examples
We begin by analyzing RWP activity in time-height sections. Figure 1 shows the evolution of the amplitudes of wavenumbers 4 to 8 during 2008-2009 at 50 • N, filtered from gridded GNSS-RO temperature fields. The 50 • N latitude represents the one with most wave activity in the NH extratropics. The highest amplitudes (orange and red shading) are found directly above the zonal-mean lapse-rate tropopause (magenta line) and in the upper troposphere, with a stark minimum in temperature signals in between (blue shading). This vertical distribution of wave activity in terms of temperature tightly follows the zonal-mean tropopause over time. Wavenumbers 4 and 5 reach amplitudes of 6-7 K, with higher wavenumbers showing lower amplitudes: 3-4 K for w8. The activity of all wavenumbers in general seems to be confined within westerly winds (solid grey lines): very little penetrates into the summer easterlies (dashed white lines). Interestingly, stratospheric zonal-mean zonal winds of 10-20 m s −1 are no impediment for wave activity to propagate beyond 20 km height, which can be seen most clearly for wavenumbers 4 to 6: they show amplitudes in excess of 4 K around 20-26 km height before and during a sudden stratospheric warming (SSW) with a central date on 24 January 2009, whose wind anomalies at 50 • N can be seen until April in the lower stratosphere, with a patch of easterly winds (thick white lines) starting in February 2009 at 26 km. All wavenumbers show temperature signals of 2-3 K (yellow shading) reaching 18-20 km height very often. The observed GNSS-RO temperature signals of wavenumbers 4-8 reaching that high into the stratosphere in Fig. 1 is in con-tradiction with the Charney-Drazin criterion (Charney and Drazin, 1961).
In Fig. 1 it may also be noticed that the temperature signals in the lower stratosphere (LS) and the upper troposphere (UT) tend to amplify and dissipate simultaneously, apart from showing very similar amplitudes. For wavenumbers 6 and higher this happens throughout the whole year; for w5 and w4 it happens less during the summer months. The vertical propagation of the temperature signals is difficult to discern in Fig. 1. There are hints of rapid upward propagation in the stratosphere which are better visible in the analysis of pressure signals, which will be discussed later in this section.
It should pointed out that the wave signals shown in Fig. 1 mostly belong to traveling waves: if the same figure is done with signals filtered also in time (eastward-propagating with 2-20 d periods), it looks very similar (not shown). We proceed to repeat the analysis from Fig. 1 on GNSS-RO gridded pressure fields. The filtered pressure anomalies are displayed relative to the mean pressure of their vertical level (in percentage units) to make the comparison between tropospheric and stratospheric signals fairer, given the exponential decrease in pressure with height. Also, this notation is proportional to wave amplitude in terms of geostrophic meridional winds as explained in Sect. 2.4. Figure 2 shows time-height sections of the wave pressure amplitudes (in percent of the zonal-mean pressure) of wavenumbers 4 to 8, also for 50 • N and 2008-2009. The same bursts of wave activity can be observed as in Fig. 1, but the most striking difference is that the pressure anomalies in Fig. 2 maximize around the tropopause height, following it constantly throughout the year for all studied wavenumbers. Wavenumbers 4 and 5 reach amplitudes of 3 %-4 % (dark red and brown shading), diminishing towards higher wavenumbers: ∼ 1.5 % for w8.
In Fig. 2 the wave pressure anomalies appear to form around the tropopause height and radiate outward vertically, diminishing their amplitude away from the tropopause. Wavenumbers 4 and 5 (and w6 in smaller amounts) show frequent upward propagation in the winter stratosphere well beyond 20 km height. Although generally decaying with height, in some cases higher wave amplitudes (red shading) are retained even at 26 km height: w4 at the beginning of December 2008 and mid-January 2009 or w5 in January 2009. Wave activity in terms of pressure of w7 and w8 tends to be confined closer to the tropopause, with very little activity reaching beyond 18 km height in Fig. 2. Compared to Fig. 1, the wave pressure anomalies from Fig. 2 show a much betterdefined continuity in their upward propagation in the stratosphere, making them much easier to track in time-height sections. For example, w4 and w5 temperature anomalies in January 2009 in Fig. 1 show no attachment between the lower and middle stratosphere (red patches are separated by blue regions of very low temperature amplitudes). Meanwhile, the corresponding pressure anomalies in Fig. 2 show clear continuity and upward propagation from the lower to the middle stratosphere. This can be discerned in the w4 bursts in the beginning of December 2008, mid-March 2009 and mid-April 2009, with the red and yellow tracks being tilted upward in the positive time direction in Fig. 2. In January 2009, the vertical propagation of several w5 bursts up to 26 km height is observed coinciding with the 2009 SSW. In this case, w4 vertical propagation seems to be very rapid in January 2009 (the red track appears quasi-vertical), maintaining very high wave amplitude throughout the whole lower and middle stratosphere. Interestingly w4, w5 and w6 to some degree show some re-amplification of their pressure anomalies near the zero-wind line (thick solid white line in Fig. 2), which is propagating downward during the 2009 SSW. This should not be interpreted as downward propagation of the RWP pressure anomalies, only that they encounter their critical level for upward propagation there, amplifying before breaking.
Overall the temperature wave structures in Fig. 1 appear much noisier and with more discontinuities compared to the same analysis done with pressure data in Fig. 2, which provides a much cleaner view of wave activity location and propagation. One explanation for this could be that the pressure data show the dynamical structure of the wave which remains together as one system. The thermal structure of the same wave is subject to more variations due to heat fluxes, meridional temperature gradients, energy conversion, etc., which feed the wave but need not be proportional to its dynamical structure locally. To our knowledge, our study is the first to analyze wave activity in terms of pressure with GNSS-RO data: only temperature or parameters derived from it have been used to study waves thus far. From the results presented in Fig. 2, the benefits of analyzing pressure GNSS-RO data for wave analysis are noticeable.
From Figs. 1 and 2 it can be concluded that the separate temperature anomalies in the UT and LS and their simultaneous amplification are part of the same pressure system that goes across the UTLS and occasionally propagates far up into the stratosphere, more often for wavenumbers 4 and 5. The wave pressure anomalies in Fig. 2 that maximize around the zonal-mean tropopause height can be considered equivalent to the upper-level PV anomalies (Hoskins et al., 1985) which then couple with surface circulation. Our results are in agreement with Teubler and Riemer (2016), who showed PV anomalies maximizing around the tropopause height. Since our gridded GNSS-RO dataset includes altitudes of 6 km and higher, the wave behavior in the lower troposphere cannot be diagnosed to study its coupling with the near-tropopause anomalies, but this will be the subject of future work in combination with other datasets.
Another topic of high interest for future work is the behavior of RWP activity, namely w4 and w5, and the possible influence on SSWs since in Fig. 2   with the troposphere: Karpechko et al. (2017) showed that northern annular mode (NAM) anomalies at the 150 hPa level define which SSWs affect the troposphere. RWP activity is known to force low-frequency variability of the zonal-mean circulation (Webster and Keller, 1975;Randel and Stanford, 1985).
Given the large amounts of RWP activity we diagnosed in the LS (and stratosphere during SSWs; see Fig. 2), RWP interaction with the LS mean flow through which plane-tary waves propagate should be expected: this would affect LS NAM conditions that control SSW downward propagation (Karpechko et al., 2017) and indirectly affect the phase speeds of upward-propagating planetary waves (Domeisen et al., 2018) by Doppler shifting. One may even speculate RWPs could directly affect planetary wave phase speeds and amplitudes in the LS if some wave-wave interaction and energy transfer or positive interference happened between the two. In either case, we state that RWP activity in the LS is potentially relevant for both onset and downward propagation stages of SSWs, opening an interesting field of research. Figures 1 and 2 have shown an example for 1 year (July 2008 till June 2009) and one latitude band: 50 • N. Section 3.2 will generalize the results shown above, giving climatological statistics for all extratropical latitudes from a decade (2007-2016) of GNSS-RO data.

Climatological distribution of wave activity over latitude and height
How is the wave activity of individual wavenumbers distributed over the whole extratropical latitude range? Taking the zonal-mean tropopause level as a reference, since wave activity in terms of pressure follows it and maximizes there (Fig. 2), we calculate the mean amplitude of each wavenumber's filtered pressure signal. We select the eastward-propagating components with 2-20 d periods, in order to isolate traveling waves throughout the NH extratropics. The resulting wave spectrum is shown in Fig. 3a for the whole NH gridded GNSS-RO dataset (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016) and separated into winter and summer climatologies in Fig. 3b and c.
The spectrum from Fig. 3a has a similar shape to the spectrum from Wolf and Wirth (2017) (their Fig. 6). Although the parameter used by Wolf and Wirth (2017) is 300 hPa meridional winds and the time filtering is slightly different (30 d high-pass), their spectrum also shows the two relative maxima of wave activity: one at midlatitudes around w6 and the other at polar latitudes and lower wavenumbers, similarly to our Fig. 3a. Traveling waves tend to have a similar range of wavelengths at all latitudes, and the spectrum shape results from their resulting zonal wavenumber Fourier composition at the different latitudes.
The wave spectrum for NH winter (Fig. 3b) is more elongated towards lower latitudes, reaching 30 • N. Meanwhile in the NH summer spectrum (Fig. 3c), there is very little wave activity south of 45 • N, which agrees well with the seasonality of the jet stream position. Interestingly, transient wave activity in the Arctic region (poleward of 70 • N, w1-3) seems to be very active all year round at the zonal-mean tropopause level.
We take the time series of the amplitude of wave temperature anomalies 1.5 km above and 3 km below the zonal-mean tropopause at each latitude band and separate wavenumber as a measure of LS and UT thermal wave activity. These heights were chosen to avoid the minimum in wave activity in terms of temperature around and right below the tropopause in Fig. 1. The correlation of the LS and UT time series will indicate the degree of their simultaneous amplification, which is shown in Fig. 3d-f. After seeing Figs. 1 and 2, one would expect a tendency of the upper and lower rows in Fig. 3 to match, but this is not the case: higher values of the measure for simultaneous amplification are shifted towards higher wavenumbers at all seasons. Furthermore, the NH summer (Fig. 3f) shows markedly higher correlations than the NH winter season (Fig. 3e). The results shown in Fig. 3d-f indicate that the simultaneous amplification of wave temperature signals between the UT and LS is done more towards the synoptic scale and the upper end of the medium scale of the wave spectrum, any season. This is an unexpected finding. One possibility of explaining it would be that the temperature signal requires meridional advection to be formed, and meridional winds have been shown to have a trough-ridge asymmetry Wirth, 2015, 2017) due to their semi-geostrophic nature, with the strongest meridional winds found closer to the trough's center. Therefore, due to non-geostrophic processes, meridional winds and temperature could typically have slightly shorter wavelengths compared to their corresponding pressure anomalies' wavelength. Proving this would require the use of other datasets beyond GNSS-RO; thus this is beyond the scope of our study.
We repeat the same analysis from Fig. 3 with RWP activity latitudinal distribution and the UT-LS simultaneous amplification measure in the Southern Hemisphere (Fig. 4). The wave spectra, as well as the UT-LS simultaneous amplification measure, are qualitatively similar in the NH (Fig. 3) and SH (Fig. 4). The SH shows maximized wave activity at w4 and ∼ 55 • S throughout the year, with markedly higher values in terms of mean relative pressure amplitude: SH maxima reach ∼ 1.3 % in Fig. 4a-c vs. ∼ 0.7 % in the NH in Fig. 3a-c. SH midlatitudes show increased eastward-propagating wave activity overall, maximized at lower wavenumbers compared to the NH midlatitudes, in agreement with previous studies (see Randel and Stanford, 1985, and references therein). This interhemispheric difference in the wave spectrum can be explained by the absence of mountains in the SH: a model study by Hayashi and Golder (1983) noted a similar change in the wave spectrum when NH mountains were removed.
The maximum of SH transient wave activity shows little seasonal variation in the latitude where it is located (Fig. 4b  and c), although in winter it shows a broader latitudinal extent, indicating larger meridional advection scales. In contrast to the NH, the SH polar latitudes show very little wave activity during summer months (Fig. 4c). As in the NH, the simultaneous amplification of SH wave temperature signals in the UT and LS is maximized at higher wavenumbers and during the summer months ( Fig. 4d-f). The SH distribution of the UT-LS correlation factor, relative to the distribution of transient wave activity, is very similar to the NH both in magnitude and shape.
We proceed to show a climatology of the vertical distribution of RWP activity in both hemispheres. We noted earlier in Figs. 1 and 2 that the temperature signal had a relative minimum near and right below the zonal-mean tropopause, and a maximum located right above it, while the pressure signal maximized and followed the zonal-mean tropopause. We take the standard deviation (SD) of zonal wavenumbers 4-8 daily temperature and pressure anomalies and average the daily SD value relative to the zonal-mean tropopause level at both NH and SH midlatitudes. We select w8 as the upper limit because we note very little wave activity at higher wavenumbers from Figs. 3 and 4. The resulting mean vertical profiles of RWP activity are shown in Fig. 5.
The vertical distribution of the RWP temperature signals in Fig. 5a is very similar in both NH and SH midlatitudes, with tropospheric and LS relative maxima above ∼ 3 K and a minimum close to and below the zonal-mean tropopause level of ∼ 1.5 K. The pressure signals in Fig. 5b maximize around and below zonal-mean tropopause height, with the SH reaching slightly higher values (∼ 1.8 %) than the NH (∼ 1.5 %), but they are otherwise alike in shape. The higher RWP activ-ity in the SH was already noted comparing the wave spectra from Figs. 3a-c and 4a-c. The climatology in Fig. 5 shows the expected vertical distributions from the 50 • N examples from Figs. 1 and 2 and shows very little interhemispheric differences. Our climatological profile of RWP temperature signals (w4-8) qualitatively agrees very well with the seasonal averages for individual wavenumbers (w4-7) by Randel and Stanford (1985) (their Fig. 4) and also with the vertical distribution of sensible heat fluxes in Blackmon and White (1982) (their Figs. 4-5). Our climatological profile of RWP pressure signals shows very good agreement with the vertical distribution of momentum fluxes in Blackmon and White (1982) (their Figs. 7-8) as well as the seasonal vertical profiles of w5 geopotential height amplitudes in the SH in Hirooka et al. (1988) (their Fig. 8a). Also, early model experiments predicted RWP activity in terms of geopotential or eddy kinetic energy to maximize around the tropopause (Gall, 1976;Simmons and Hoskins, 1978).

Zonal and vertical structures of RWPs
A classical RWP example with eastward phase propagation and faster group speed will be presented in Sect. 4.1 with the use of Hovmöller diagrams and longitude-height snapshots of RWP anomalies from the gridded GNSS-RO data. Recent studies have pointed out that under specific resonance conditions, RWPs can become stationary and form large-scale teleconnection patterns in midlatitudes, which have been linked to the extreme heat waves of 2003, 2010 and 2018 (Petoukhov et al., 2013;Kornhuber et al., 2017Kornhuber et al., , 2019. To test whether RWP properties in the UTLS are different under these conditions, we will repeat the same analysis from Sect. 4.1 during the Russian heat wave in summer 2010, which will be presented in Sect. 4.2. After these two casestudies, in Sect. 4.3 we will present climatological statistics of the vertical scale of RWP anomalies, the relation of UT and LS zonal temperature structures, and the zonal scale of RWP envelopes for both NH and SH.

Classical RWP example
RWP activity was shown to follow the zonal-mean tropopause in Figs. 1 and 2; therefore we take this level to extract RWP pressure anomalies from the gridded GNSS-RO data. This case study was selected to show a classical RWP with eastward phase propagation and a high-amplitude envelope of faster speed, which is the case at the end of  Zimin et al., 2003;Wolf and Wirth, 2015). We select pressure anomalies belonging to wavenumbers 5-8 in order to avoid stationary waves with w1-4 present at the time. The evolution of the RWP pressure anomalies is shown in the Hovmöller diagram in Fig. 6a, and the RWP envelope after applying the Hilbert transform (Zimin et al., 2003) to the w5-8 anomalies is shown in Fig. 6b.
In Fig. 6b the formation of an RWP around 27 February 2009 over the Pacific Ocean (∼ 150 • W) can be observed, with the envelope propagating eastward until 7 March 2009 when it covers the Atlantic Ocean (60-0 • W). The eastward phase propagation of the RWP (individual blue troughs and red ridges) in Fig. 6a can be compared to the markedly faster movement of the envelope (green and dark colors) in Fig. 6b. This RWP reaches amplitudes exceeding 4 %. By 7 March 2009, a new RWP is forming again over the Pacific Ocean, this one with an even faster envelope and circumnavigating the globe by 12 March 2009, and even reaching global scale for a couple of days although having lower amplitudes than the first RWP, of around 2 %-3 %. It also has to be noted that low-amplitude fluctuations of ∼ 1 % within the w5-8 range are present almost constantly (e.g., light yellow shading in Fig. 6b). Recent studies using meridional winds and the Hilbert transform have shown the benefits of refining the methodology with semi-geostrophic coordinate adjustment Wirth, 2015, 2017) and even using wave activity flux (Takaya and Nakamura, 2001) for RWP diagnostics and optimizing the envelope's shape. Our approach to describe the phase and envelope propagation of an RWP in Fig. 6, with the only use of relative pressure anomalies from gridded GNSS-RO and the Hilbert transform over a limited wavenumber range (5-8), is quite simplistic in comparison but enough to get a good qualitative view. We concentrate now on the first RWP in Fig. 6 (27 February to 7 March 2009), exploring how its pressure and temperature structures evolve in the UTLS. Longitude-height snapshots of this RWP are presented at 2 d intervals in Fig. 7.
The RWP pressure anomalies in Fig. 7 (top row) are centered around the tropopause level, while the temperature anomalies (bottom row) have a clear separation between the UT and LS and a noisier appearance. This is in agreement with the results from the zonal-mean perspective from Figs. 1 and 2. The RWP pressure anomalies have a near-global coverage near the tropopause at almost all times, despite having segments of very low amplitudes: e.g., ∼ 60 • E in 1 March 2009 and ∼ 120 • E in 3 March 2009 in Fig. 7, top row. The zonal extent of the RWP pressure anomalies decreases away from the tropopause: e.g., compare pressure anomalies at 10 and 16 km in Fig. 7. The choice of a specific vertical level to extract the RWP's envelope could affect the extent of the diagnosed RWP (see Discussion,Sect. 5).
Looking at the individual phases of the RWP pressure anomalies, the same troughs and ridges can be recognized going from 6 km up to ∼ 20 km height in Fig. 7 (top row). A slight westward tilt with height can be spotted in the stratosphere, which is typical of Rossby waves, although the RWP pressure structures in the UTLS can be considered almost barotropic. From the RWP pressure anomalies in Fig. 7 (top row) it can be concluded that RWPs form a direct dynamical connection between the UT and the stratosphere.
The RWP temperature anomalies in Fig. 7 (bottom row), apart from the break around tropopause level, show an outof-phase behavior between the UT and LS: troughs (blue phases in Fig. 7, top row) correspond to negative temperature anomalies in the troposphere and positive temperature anomalies in the stratosphere. This is expected from the meridional advection associated with the barotropic pressure anomalies and the inversion of the meridional temperature gradient between the UT and LS (Volkmar Wirth, personal communication, 2019). However, the meridional temperature gradients in the lowermost stratosphere are of low magnitude, so it is not clear why the RWP temperature anomalies maximize so close to the tropopause and not higher up. Pilch Kedzierski et al. (2017) also reported wave amplitudes maximizing within a short height range close to the tropopause.
The longitude-height structures of RWP pressure and temperature anomalies in the UTLS from Fig. 7 show a high resemblance with those reported by Hakim (2005) (e.g., see their Fig. 7), although they used the leading empirical orthogonal functions (EOFs) of meridional winds and temperature. These structures are strongly related to forecast error propagation (Hakim, 2005), with wind errors maximizing near the tropopause, while temperature errors maximize near the surface. It is possible to collocate forecast error directly with observed RWP structures like those we show in Fig. 7, another interesting topic for a future study.
Recent case studies about extreme Arctic cyclones (Tao et al., 2017a, b) highlighted the influence of strong positive PV anomalies in the LS with a warm core in the deepening of the surface cyclones, which would be consistent with the presence of an RWP with UTLS pressure and temperature structures such as those depicted in Fig. 7, a high-amplitude trough in this case. Odell et al. (2013) noted large stratospheric geopotential height tendencies during the 1993 Braer storm, the deepest extratropical low on record with 914 hPa minimum surface pressure. A study of 60 severe European windstorms (Pirret et al., 2017) found that in 20 cases the stratosphere contributed 10 % or more to their deepening. The pressure tendency equation used in Pirret et al. (2017) used the 100 hPa level (∼ 16 km) as the upper boundary, with the stratospheric contribution being dependent on geopotential height tendencies at this boundary (their dPhi term). The presence of a high-amplitude RWP in the UTLS can significantly affect geopotential anomalies at the 100 hPa level, as our results show (Fig. 7, top row; also Fig. 2 indicates pressure anomalies reaching this level very often). Therefore we hypothesize that 100 hPa geopotential height tendencies (used to calculate the stratospheric contribution to surface cyclone deepening) could be very sensitive to the presence of RWP structures in the UTLS.

The 2010 Russian heat wave case
Recent studies have shown that in NH summer and with a specific setting of the waveguide(s), resonant conditions appear and drive high-amplitude waves which become quasistationary, with wavenumbers in the w6-8 range (Petoukhov et al., 2013;Kornhuber et al., 2017Kornhuber et al., , 2019. This process has been linked to the 2003, 2010, 2015 and 2018 heat waves in Europe and Russia. Is the appearance of RWPs in the UTLS any different under these special conditions? Throughout this section we will seek an answer by repeating the diagnostics of Sect. 4.1 for the summer 2010 Russian heat wave. Figure 8a shows the evolution of RWP pressure anomalies (w4-8) during summer 2010 and Fig. 8b their corresponding envelope calculated with the Hilbert transform. It can be observed in Fig. 8b that the RWP envelope goes around the globe one and a half times between 10 June and 10 July. Eastward phase propagation is generally visible throughout this period, except for the time when the RWP is around 0 • longitude, where the phases become nearly stationary (Fig. 8a). Between 10 July and 10 August, a trough (blue) sets in around 0 • and barely moves, with ridges constantly present to its west and east. The ridge around ∼ 30-45 • E is present nearby Moscow during this time. It is not until 20 August when a trough passes over Moscow's longitude.
Rather than one RWP expanding and becoming stationary, the envelope evolution in Fig. 8b suggests recurrent Atlantic RWPs with eastward group propagation terminating over 0-60 • E, with near-zero phase speeds over this region between 10 July and 10 August. The amplitudes of the multiple RWPs in Fig. 8 are not higher than those of the typical winter case in Fig. 6, meaning that RWP anomalies in the UTLS need not be very strong to cause extreme events on the surface. The summer 2010 case is unusual for the low phase speeds around 0 • longitude and the recurrence of RWPs coming into the 0-60 • E region, giving the event a much longer duration.
In Fig. 9 we present longitude-height snapshots of the summer 2010 RWPs. Given the longer event duration, we take a total of eight snapshots at 4 d intervals. RWP pressure anomalies in the UTLS (1st and 3rd rows) are maximized near the tropopause level and are quasi-barotropic, while the RWP temperature anomalies (second and fourth rows) maximize above the tropopause with a clear break between the UT and LS signals. No qualitative differences are found in the UTLS structures of RWPs between the typical winter case (Fig. 7) and the recurrent RWPs with stationary phases during the summer 2010 Russian heat wave (Fig. 9). A similar analysis of longitude-height snapshots during the 2015 European heat wave leads to the same qualitative conclusions (not shown) in terms of RWP pressure and temperature structures in the UTLS.
We conclude that the general behavior of RWPs across the UTLS, dynamically and thermally, is the same in normal winter conditions (Figs. 6 and 7 in Sect. 4.1) and summer under quasi-resonant waveguide setting (Figs. 8 and 9 in Sect. 4.2). The only differing factors during the summer 2010 heat wave are the low RWP phase speeds and their recurrence, although their UTLS appearance is alike to typical winter RWPs.

Vertical scale analysis
Figures 7 and 9 showed the difference between the RWP appearance in terms of pressure and temperature anomalies: whereas pressure anomalies have a long vertical wavelength, temperature anomalies show a break around the tropopause and opposite phases in the UT and LS. Therefore it is expected that the vertical scale of RWP temperature anomalies is significantly shorter than that of RWP pressure anomalies, especially in the UTLS. To quantify this, we take all longitude-height snapshots of RWP (w4-8) anomalies between 40 and 60 • latitude, computing Fourier power spectra in the vertical direction between 6 and 36 km height. The resulting average Fourier power spectra for both hemispheres are shown in Fig. 10a. Temperature (solid blue and dashed cyan lines) and pressure (red and dashed orange lines) vertical power spectra Figure 9. Longitude-height snapshots of (fist and third rows) pressure and (second and fourth rows) temperature anomalies of the RWP from Fig. 8 at 55 • N. The grey line denotes the gridded lapse-rate tropopause height.
for NH and SH show very little interhemispheric differences in Fig. 10a. Both RWP temperature and pressure anomalies (w4-8) have increasing power towards longer vertical wavelengths. The pressure power spectra show the steepest increase in power for vertical wavelengths > 10 km, with a still noticeable amount of power between 5 and 10 km vertical wavelengths. Temperature vertical power spectra show the steepest increase in power at wavelengths > 5 km and a noticeable amount of power between 3 and 5 km wavelengths. Also, the integrated power of 1.5-3 km vertical wavelengths might be of significance in the temperature spectrum. Unlike pressure, the temperature structures of RWPs have a relatively large amount of power at vertical wavelengths 3-10 km in their Fourier spectra, and those are the vertical wavelengths where the power ratio between temperature and pressure is highest in Fig. 10b.
We repeated Fig. 10 selecting the same zonal wavenumbers (4-8) and additionally filtering in the time dimension to obtain 2-20 d eastward-propagating waves. This makes sure that neither stationary waves nor gravity waves (that oscillate at frequencies higher than the inertial frequency) are present. The result (see Fig. S4) is a near-identical spectrum.
In an atmosphere model, the vertical resolution (or the separation between model levels) needed to resolve wave temperature structures is ∼6 times less than the wavelength. Most global climate models (GCMs) have vertical resolu- tions coarser than 1 km in the UTLS, meaning that they will struggle to represent vertical wavelengths shorter than 6-8 km. From the temperature vertical power spectra in Fig. 10a one could anticipate that RWP temperature structures in the UTLS such as Figs. 7 and 9 may be partially underrepresented by GCMs. In fact, GCMs, reanalyses and forecast models are known to produce temperature gradients in the tropopause region that are too smooth Hegglin et al., 2010;Birner et al., 2006;Pilch Kedzierski et al., 2016b). The tendency of forecast models to smooth the tropopause with lead time (Gray et al., 2014) is partially compensated for by diabatic and parameterized processes in the modeled RWPs (Saffin et al., 2017).
A better understanding of RWP diabatic processes offers the chance of forecast improvement; therefore a variety of frameworks to analyze RWP evolution have been developed, mainly using reanalysis or forecast model output (e.g., Hoskins et al., 1985;Orlanski and Katzfey, 1991;Nielsen-Gammon and Lefevre, 1996;Chang, 2001;Gray et al., 2014;Teubler and Riemer, 2016). The vertical structure of such processes is of high importance for RWP evolution, and analyses such as our Figs. 1,5,7,9 or 10 would enable a comparison of model or reanalysis output directly with GNSS-RO observations. See Sect. 5 for a more detailed discussion.

UT and LS out-of-phase temperature behavior
Figures 7 and 9 indicate that the RWP thermal structures in the UT and LS have constantly opposing phases. If we take the zonal structures of RWP temperature anomalies 1.5 km above and 3 km below the zonal-mean tropopause as proxies for the LS and UT, respectively, in the case of Figs. 7 and 9 the correlation factors between the LS and UT would be close to R = −1. We test the prevalence of this behavior by creating probability density functions (PDFs) of the UT-LS correlation factor R against the wave's amplitude (as the standard deviation of the wave pressure anomalies at zonalmean tropopause level) for different wavenumber ranges and latitude bands in Fig. 11.
The PDF for midlatitude RWPs (w4-8, Fig. 11a) shows a general tendency for UT-LS anticorrelation of the temperature structures: higher densities (yellow-red colors) are increasingly packed towards R = −1 the higher the RWP pressure amplitude gets. The RWP temperature structures with opposite phases between the UT and LS are thus a typical behavior with very few exceptions as one can guess from the near-zero probability densities away from R = −1 in Fig. 11a.
Midlatitude planetary waves (w1-2, Fig. 11b) show a much more disperse distribution resulting in lower densities (grey-black colors) covering the whole R range, evenly distributed especially at higher wave amplitudes. In Fig. 11b, the relative concentration towards R = −1 at lower wave amplitudes could be the result of the w1-2 components being part of an RWP, as explained next. One may consider an RWP with carrier w4 (e.g., the idealized example in Wolf and Wirth (2015); their Fig. 2d) and its zonal wavenumber Fourier spectrum: the RWP wavenumber spectrum peaks at w4, with decreasing power of the wavenumbers around it, and w1-2 are low-amplitude contributors to the RWP zonal structure the same way as w6-7. Higher-amplitude planetary waves, which cannot be part of RWPs, do not show such a tendency of the PDF to maximize near R = −1.
In the subtropics ( Fig. 11c and d), both RWPs and planetary waves have lower amplitudes than their midlatitude counterparts. This is expected since Figs. 3 and 4 show very little transient (2-20 d period) wave activity there, and any wave amplitude of importance in Fig. 11c and d, which has not been filtered in the time dimension, should therefore come from stationary waves. The PDF for subtropical RWPs (w4-8, Fig. 11c) is concentrated at low amplitudes and shows very little shift towards R = −1, being spread throughout all R values and looking nothing like the midlatitude PDF.
Planetary waves in the subtropics show a PDF packed near R = −1 (w1-2, Fig. 11d). As mentioned above, higher amplitudes of w1-2 belong to stationary signals, which translate into large-scale and long-lasting temperature anomalies of opposite sign in the subtropical UT and LS. A plausible explanation for Fig. 11d is represented by ENSO-related temperature anomalies in the UT and LS: they extend to the sub- tropics, are long-lasting and their zonal scale fits well with a combination of zonal wavenumbers 1-2. Indeed Domeisen et al. (2019), their Fig. 5, showed opposite-signed UT and LS temperature regressions to the ENSO index. UT warm anomalies relate to increased convective latent-heat release and LS cold anomalies from the resulting enhanced largescale upwelling. These anomalies would combine into the anticorrelation shown in Fig. 11d.
We repeated the analysis of Fig. 11 on the SH, finding nearly identical results (see Fig. S5 in the Supplement).

Zonality of RWP envelopes: hemispheric comparison
Early studies noticed the tendency of RWP envelopes to have a near-global scale in the SH, while the NH shows a more zonally confined behavior (see Lau, 1979;Randel and Stanford, 1985, and references therein). Hemispheric comparisons by Chang (1999) and Souders et al. (2014b) with more advanced methods to define RWPs support these earlier findings. A Fourier analysis of the RWP envelope describes its zonal scale: the w0 component of the envelope corresponds to the global (or zonal-mean) extent, while w1 and higher components modulate how the envelope's amplitude varies regionally. Here we want to find out whether the SH has higher amounts of all envelopes' Fourier components or only of the longest ones. A comparison of the wavenumber power spectra of RWP envelopes from NH and SH midlatitudes is shown in Fig. 12. A surprising similarity of the RWP envelope's wavenumber spectra for both hemispheres can be observed in Fig. 12 (solid black and dashed grey lines; also note the logarithmic y scale). There is a steady exponential increase in power from w10 to w5, appearing as a line of constant slope in Fig. 12 between values from 0.5 to 2. The increase in power from w4 to w1 is well above exponential, with a jump-like shift to values from 10 to ∼ 30. The w0 component has the highest power (∼ 150), which is expected because an RWP envelope is by definition always positive. It can be concluded that at tropopause level, w0-4 components markedly dominate the RWP envelope spectrum, with higher wavenumbers contributing marginally. This is observational support for previous studies which filtered out higher wavenumber components of the RWP envelope in their methodology: e.g., Wolf and Wirth (2015) eliminated the envelope's components >w6 to avoid undesired small-scale wiggles in the envelope's shape.
To highlight where the differences between hemispheres are, we add the NH/SH power ratio for each wavenumber as the dashed blue line in Fig. 12. It can be observed that for w3-10, the NH/SH ratio is very close to 1: it ranges from 0.95 to 1.03. The NH/SH power ratio becomes lower at larger scales: 0.92 for w2, 0.88 for w1 and 0.85, the biggest hemispheric difference, for w0. This means that the longest wavenumber components (w0 and w1) of RWP envelopes in the NH climatologically have 85 % and 88 % of the power found in the SH spectrum. Our results in Fig. 12, with the SH having more w0 and w1 in the RWP envelopes, are consistent with previous studies (Lau, 1979;Randel and Stanford, 1985;Chang, 1999;Souders et al., 2014b) in that the zonal extent of RWPs in the SH tends to be more global than in the NH. In addition, we show that the amount of zonal variability of the envelope's amplitudes at a sub-planetary scale (w2 and higher) is basically the same in both hemispheres.

Discussion
The following points outline the importance of following the tropopause.

For diagnosing RWP properties
Throughout Sects. 3 and 4 our results depict RWP activity being centered around the tropopause height, following it at all times. At midlatitudes, the tropopause in the NH varies from ∼ 220-300 hPa in winter to ∼ 175-250 hPa in summer, while in the SH it varies from a fairly meridionally uniform ∼ 250 hPa in winter to ∼ 200-280 hPa in summer (Son et al., 2011). This has important implications for RWP diagnostics, which are usually performed with meridional winds at a fixed pressure level, which is not equidistant from the tropopause over latitude and season. The most widely used pressure level for RWP diagnostics is 300 hPa, either for showcasing RWP envelope reconstruction or tracking methods (Zimin et al., 2003(Zimin et al., , 2006Souders et al., 2014a;Wirth, 2015, 2017), producing RWP climatologies or composites (Chang, 1999;Chang and Yu, 1999;Williams and Colucci, 2010;Souders et al., 2014b), or studying their predictability (Quinting and Vitart, 2019). The 250 hPa level was also used by Glatt and Wirth (2014) and Grazzini and Vitart (2015) to study RWP properties and predictability in forecast models.
Generally a fixed threshold for RWP envelope amplitude is defined to determine whether an RWP is present and delimit the RWP's horizontal extent. Having the longitudeheight snapshots of RWP anomalies from Figs. 7 and 9 in mind (relative pressure RWP amplitudes being proportional to geostrophic meridional wind amplitudes, as explained in Sect. 2.4), the diagnosed RWP amplitude and extent would increase the closer the chosen vertical level is to the zonalmean tropopause. For a given pressure level, its distance from the tropopause level will vary seasonally and over latitude, introducing an aliasing factor if RWP properties are compared between summer and winter or NH and SH. Some-times this issue could even affect the number of diagnosed RWPs: if the pressure level is further from the tropopause than usual, some part of an RWP may show meridional winds below the threshold, which leads to fragmentation into two or more RWPs; or a low-amplitude RWP may not be detected at all, while at tropopause level the threshold is still exceeded. Meridional wind on the 2 PVU surface (the so-called dynamical tropopause) is an available reanalysis product that is a perfect candidate to avoid the above-discussed aliasing effect on RWP properties. Some methods make use of a threshold relative to the zonally averaged RWP envelope amplitude (e.g., Glatt and Wirth, 2014;Grazzini and Vitart, 2015;Quinting and Vitart, 2019). While solving the issue of pressure level distance to the tropopause, in this case an even bigger problem arises: the presence of stronger wave activity over a near-global longitude range would increase the threshold substantially, reducing the diagnosed RWP size (Wolf and Wirth, 2017). Also, in the case of a zonally localized RWP growing into a highamplitude and global wave due to resonant conditions (e.g., Petoukhov et al., 2013), a relative threshold would by contrast yield a shrinking RWP.

For RWP energy budgets and fluxes
RWP eddy kinetic energy (EKE) budgets are usually presented as vertical integrals, from the surface up to the upper boundary at 100 hpa (Orlanski and Katzfey, 1991;Chang, 2000Chang, , 2001 or 50 hPa (Orlanski and Sheldon, 1993). Additionally, Chang (2001) showed volume-averaged composites of RWP wave activity fluxes for the UTLS (400-200 hPa), noting that the baroclinic growth -barotropic decay paradigm of RWP life cycles has to be treated with caution when interpreting zonally localized waves. PV inversion diagnostics by, e.g., Nielsen-Gammon and Lefevre (1996) use 1000 and 100 hPa as boundaries for the PV tendency equation, presenting different terms and compositing at a fixed pressure level (e.g., near 300 hPa).
Our results suggest that compositing energy budget terms and fluxes relative to the tropopause height, or even making separate budgets for the UT and LS, would add precision to the magnitude, location and role of each term. Winddependent terms would maximize near the tropopause, since our RWP pressure anomalies are proportional to geostrophic winds. The temperature-dependent terms would show high sensitivity to their location relative to the tropopause, since our RWP temperature signals maximize right above the tropopause and become very low right below it (see Figs. 1, 2, 5, 7 and 9). Averaging on pressure levels may lead to smoothing of the vertical structures of specific EKE budget or PV tendency terms. Tropopause-based averaging retains UTLS gradients as introduced by Birner et al. (2002) and Birner (2006) for radiosonde temperature structures.
Recent publications are already applying a tropopauserelative framework to study PV tendencies (e.g., Cavallo and Hakim, 2009;Saffin et al., 2017), compositing for cyclones and troughs and ridges. Such an approach combined with a refined method to localize RWP envelopes (see discussion above) could improve the understanding of the forcings involved in RWP life cycles and their interaction with the background flow, e.g., in 3-D composites similar to those presented by Chang (2001).

For diagnosing wave trends
Quantifying the waviness of the extratropical westerly flow is challenging (see review by Coumou et al., 2018): each measure may have a different physical meaning and different degrees of complexity for systematic implementation on reanalysis or model data. Waviness trends are inconsistent among different methodologies (e.g., Barnes, 2013;Coumou et al., 2015;Cattiaux et al., 2016). These studies have in common that they perform their analyses on fixed pressure levels (e.g., 500 and 250 hPa for wave amplitudes - Barnes, 2013 andCattiaux et al., 2016 -or 850-250 hPa integrals for EKE trends - Coumou et al., 2015). Thus, they do not consider tropopause height biases in the models used for climate projections, which can be different for each model, potentially adding spread to the results. Neither do these studies account for modeled or reanalyzed tropopause height trends (e.g., tropospheric expansion leads to an increase in extratropical tropopause height over time with climate change). Interannual variability of the extratropical tropopause height alone would introduce aliasing effects over the diagnosed waviness on a fixed pressure level. Wave amplitudes diagnosed on the 2 PVU surface, as our results show, would always be equidistant from the location of the maximum in wave activity in the jet stream, thus providing the fairest trend measure and model intercomparison in terms of wave amplitudes.
The issue of tropopause trends was already noted for the attribution of stratospheric residual circulation trends by Oberländer-Hayn et al. (2016). We add that, since RWPs show a quasi-barotropic structure throughout the UTLS (See Figs. 7 and 9), RWP phase speeds in principle are not subject to the aliasing effect of varying tropopause distance from the diagnosed pressure level, so phase speed trends like those studied by Coumou et al. (2015) are not affected by this issue.

Conclusions
Our study is a first attempt to describe RWP properties globally with the sole use of GNSS-RO observations, focusing on both hemispheres' extratropical UTLS. Observational knowledge about RWPs such as the one presented here is much needed for comparison and interpretation of climate model and reanalysis output, especially with the current scientific effort to better understand storm track and jet stream responses to climate change. Our results are relevant for stratosphere-troposphere coupling: they indicate systematic RWP activity propagation into the stratosphere, which markedly increases during SSWs, in addition to the wellknown RWP driving of tropospheric weather. We summarize our main findings below.
1. RWP activity follows the zonal-mean tropopause level and maximizes around it during all seasons (Figs. 1, 2, 5, 7 and 9). RWP pressure anomalies tend to be centered at the tropopause, with decreasing amplitude away from it. RWP temperature anomalies maximize right above the tropopause, with a contrasting minimum right below the tropopause.
2. We note frequent RWP propagation into the stratosphere, sometimes beyond 20 km height in the NH, which mostly manifests itself for wavenumbers 4 and 5 ( Figs. 1 and 2). Enhanced vertical propagation of RWP activity coincided with the SSWs of 2009, 2010 and 2013 (see Figs. S1, S2 and S3). We will explore this further in an upcoming study, given the importance of the lower stratosphere in modulating the onset and downward propagation of SSWs.
3. Since RWP activity constantly follows the tropopause, we note that the use of fixed pressure levels for RWP diagnostics or wave trends (e.g., typically 300 hPa) may induce some aliasing in the resulting quantities, since it is not equidistant from the tropopause over time and latitude. RWPs would lose envelope extent and/or amplitude the further the pressure level is from the tropopause. We suggest using the 2 PVU surface to avoid this (see discussions in Sect. 5).
4. In addition to gridded temperature, our study is the first to analyze GNSS-RO (relative) pressure anomalies, which are conveniently proportional to geostrophic meridional wind wave amplitudes (see Sect. 2.4). The filtered RWP pressure signals show much better continuity (i.e., no breaks in their vertical structure) and less noisiness in the UTLS and the stratosphere, compared to temperature wave signals (Figs. 1, 2, 7 and 9). This new approach shows many benefits for studying extratropical wave propagation with GNSS-RO data, especially in the vertical direction.
5. The dynamical and thermal appearance of RWPs in the UTLS is generally similar in different seasons and waveguide settings: no qualitative differences were observed after comparing a typical winter RWP with RWP activity during the 2010 Russian heat wave where resonant conditions were present (see Sects. 4.1 and 4.2).
6. Overall, RWPs in the SH show a preference for lower carrier wavenumbers compared to the NH (Figs. 3 and  4). In terms of envelope properties, the SH shows a higher amount of the longest scales compared to the NH