Next Article in Journal
Morphology of Traveling Wave Disturbances Recorded in Eastern Siberia in 630 nm Atomic Oxygen Emission
Next Article in Special Issue
Distinguishing Convective-Transition Moisture-Temperature Relationships with a Constellation of Polarimetric Radio Occultation Observations in and near Convection
Previous Article in Journal
Effect of Selected Meteorological Variables on Full Flowering of Some Forest Herbs in the Western Carpathians
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relationships between Extratropical Precipitation Systems and UTLS Temperatures and Tropopause Height from GPM and GPS-RO

1
COSMIC Program Office, University Corporation for Atmospheric Research, Boulder, CO 80301, USA
2
Department of Physical and Environmental Sciences, Texas A&M University–Corpus Christi, Corpus Christi, TX 78412, USA
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(2), 196; https://doi.org/10.3390/atmos13020196
Submission received: 5 November 2021 / Revised: 15 January 2022 / Accepted: 24 January 2022 / Published: 26 January 2022
(This article belongs to the Special Issue Advances in GNSS Radio Occultation Technique and Applications)

Abstract

:
This study characterizes the relationship between extratropical precipitation systems to changes in upper troposphere and lower stratosphere (UTLS) temperature and tropopause height within different environments. Precipitation features (PFs) observed by the Global Precipitation Measurement (GPM) satellite are collocated with GPS radio occultation (RO) temperature profiles from 2014 to 2017 and classified as non-deep stratospheric intrusion (non-DSI; related to convective instability) or deep stratospheric intrusion (DSI; related to strong dynamic effects on the tropopause). Non-DSI PFs introduce warming (up to 1 K) in the upper troposphere, transitioning to strong cooling (up to −3.5 K) around the lapse rate tropopause (LRT), and back to warming (up to 2.5 K, particularly over the ocean) in the lower stratosphere. UTLS temperature anomalies for DSI events are driven predominantly by large scale dynamics, with major cooling (up to −6 K) observed from the mid-troposphere to the LRT, which transitions to strong warming (up to 4 K) in the lower stratosphere. Small and deep non-DSI PFs typically result in a lower LRT (up to 0.4 km), whereas large but weaker PFs lead to a higher LRT with similar magnitudes. DSI events are associated with larger LRT height decreases, with anomalies of almost −2 km near the deepest PFs. These results suggest intricate relationships between precipitation systems and the UTLS temperature structure. Importantly, non-DSI PF temperature anomalies show patterns similar to tropical convection, which provides unification of previous tropical research with extratropical barotropic convective impacts to UTLS temperatures.

1. Introduction

The upper troposphere and lower stratosphere (UTLS) is a coupling region in the atmosphere, which is distinct in radiation, dynamics, chemistry, and microphysics [1,2]. A strong connectivity amongst these various processes makes the UTLS highly susceptible to climate change [3] and has attracted much research attention in recent decades. The UTLS is generally defined as the region ±5 km around the tropopause [4], and this boundary plays a crucial role in UTLS variability [5,6]. A fundamental characteristic of the tropopause is the change in static stability (temperature lapse rate) across the interface [7]. The WMO definition of the tropopause is based on lapse rate criteria (lapse rate tropopause or LRT) and is defined as “the lowest level at which the lapse rate decreases to 2 K km−1 or less, provided that the average lapse rate between this level and all higher levels within 2 km does not exceed 2 K km−1” [8]. The thermal tropopause provides a convenient way of defining the tropopause and is most widely used. There are also other ways of defining the tropopause, such as the potential vorticity-based dynamical tropopause [9] or the ozone-based chemical tropopause [10]. Regardless of tropopause definition, it has been suggested that changes in tropopause height can serve as a useful indicator of climate change [3,11].
Stratosphere-troposphere exchange (STE) across the tropopause is an important bidirectional process influencing the chemistry of the UTLS [5]. Many recent studies have focused on understanding the role convection plays in STE throughout the tropics [12,13,14] since radiative-convective balance is the dominant physical process in the region [15]. Deep convection influences climate processes by regulating stratospheric water vapor through direct convective injection, by enhancing thin cirrus cloud presence, and by modulating the ozone budget in the upper troposphere [10,16,17]. However, understanding STE processes due to deep convection requires accurate temperature observations in the UTLS [18], so understanding the role convection plays in the heat budget of the tropical UTLS has been a prominent research topic over the past few decades. Warm anomalies have been observed throughout the mid-to-upper troposphere using radiosonde data [16] and climate models [1]. Cool anomalies have been observed in the tropical tropopause layer (TTL) using AIRS and radiosonde data [19] along with GPS radio occultation (RO) data [20,21,22,23,24]. The resulting temperature anomalies varied greatly among these studies, and the magnitude of these anomalies depended on surface type (land/ocean) and intensity of the deep convection [1,16,19,20,21,22,23,24].
On the other hand, the relationship of extratropical precipitation to UTLS temperature and tropopause height changes has attracted much less attention because of the significant challenge of quantifying the relative contributions of a variety of dynamical processes associated with extratropical STE [25]. These processes include tropopause folds near the subtropical and polar jets due to baroclinic wave dynamics [26,27], cutoff lows [28], gravity wave breaking [29], and mesoscale convective complexes [30]. For example, significant STE occurs due to turbulent mixing during baroclinic tropopause folds [31], and these folds can even reach nearly as low as the boundary layer [32]. Additionally, simulations have shown that a large fraction of STE near folds occurs due to precipitation below the depressions in the tropopause [33]. The diversity of weather systems that occur in the extratropics can also provide considerable challenges. Similar to the tropics, precipitation can form due to convective instability (such as from single cell thunderstorms and mesoscale convective complexes), but it also commonly occurs due to large-scale dynamics (such as within extratropical cyclones). Since these precipitation systems form in vastly different environments, characterizing their relationship to the UTLS thermodynamic structure is not as straightforward relative to the tropics. However, it is critical to study both synoptic and mesoscale extratropical precipitation systems and their relationship to UTLS temperatures and the tropopause to facilitate a better understanding of extratropical STE.
Space-borne radar observations from the Tropical Rainfall Measuring Mission (TRMM) have been useful in identifying the vertical extent of convection, but were restricted to the tropics and subtropics. Over the extratropics, most passive microwave and infrared satellite data lack the detailed vertical structure of storms and their associated vertical transport. Highly sensitive radar and lidar measurements from CloudSat and CALIPSO have aided the observations of vertical cloud structure in the extratropics. However, both satellites on the A-train orbit only capture storms at fixed local times of 1:30 P.M. and 1:30 A.M., which misses the peak occurrence of land convection in the extratropics [34]. Moreover, the small swath of both satellites limits the spatial sampling of individual storms. These gaps have been filled by the Global Precipitation Measurement (GPM) mission, which was launched in February 2014. The space-borne radar onboard the GPM satellite, along with high orbit inclination, extends the vertical scan of precipitation into the high latitudes [35]. Recent research utilizing GPM has enhanced our understanding of extratropical precipitation, including identifying the global distribution of storms with large sizes [36] and overshooting convection [37]. To better understand storm structure, additional information from nearby atmospheric thermodynamic profiles have historically been obtained from radiosonde soundings with limited spatial and temporal sampling or meteorological analyses/reanalyses with coarse vertical resolution and large uncertainty near the convection. GPS-RO soundings have filled these gaps by offering global observations of UTLS temperatures with a high vertical resolution in all-weather conditions [38].
The main goal of this research is to better understand the relationship between extratropical UTLS thermodynamic structure variations and the precipitation features (PFs) underneath. Specifically, extratropical PFs under two environments (e.g., Deep Stratospheric Intrusion (DSI) and non-Deep Stratospheric Intrusion (non-DSI)) are classified using potential vorticity and investigated separately because PF development is similar in nature within each environment, as non-DSI PFs form mainly due to convective instability, whereas PFs near DSI events typically form due to baroclinic instability.
To achieve this main goal, we have three major objectives: (1) To classify and establish the UTLS temperature and tropopause height anomaly pattern for non-DSI PFs. Since PFs within the non-DSI group are similar to tropical convection, we hypothesize that the UTLS anomalies should look relatively similar to previous results shown in the tropics [1,16,19,20,21,22,23,24]. Additionally, we want to quantify how PF depth, size, surface type, and season of occurrence can impact the sign and magnitude of UTLS temperature and tropopause height anomalies. Many studies [1,16,19,20,21,22,23,24] have attributed differences in the tropical convective temperature signal to those various factors, but this has yet to be analyzed in-depth for the extratropics. (2) To identify DSI events and determine the UTLS temperature and tropopause height anomalies that occur during these events. Most of the PFs that occur during these events are very different from tropical convection and are typically associated with extratropical cyclones or strong frontal boundaries. Thus, we hypothesize that the anomaly patterns should look different from tropical convection and will be strongly influenced by synoptic-scale extratropical dynamics. (3) To determine where and when the two populations of PFs occur and what their properties are, as this knowledge will aid in providing physical explanations for the observed anomalies.
In this study, extratropical PFs are identified using GPM, while the temperature structure within these PFs is provided by collocated high-resolution GPS-RO soundings. The structure of this paper is as follows: Section 2 introduces the GPM, GPS, and ERA-Interim data used in this study; Section 3 describes the methodology used, including how the background profiles are generated and how the PFs are separated into two groups; Section 4 provides sampling for PFs in both environments along with PF characteristics; Section 5 presents the key results of the study, including UTLS temperature and tropopause height anomalies within each environment; Lastly, a discussion, including the main conclusions and study limitations, is provided in Section 6.

2. Data Description

2.1. GPM Precipitation Feature Product

The GPM mission is an international network of satellites that provide global observations of precipitation. Building upon the many successes of TRMM, which focused primarily on quantifying the three-dimensional distribution of moderate-to-heavy rain throughout the tropics [39], the advanced radar/radiometer system onboard the GPM Core Observatory extends the measurement range to include light precipitation and snow. The GPM Core Observatory is equipped with the first spaceborne dual-frequency phased array precipitation radar (DPR), which operates at the Ku and Ka bands (13 and 35 GHz, respectively), and a conical-scanning multichannel (10–183 GHz) microwave imager (GMI) [35]. The 65° orbit inclination allows for observations into the high latitudes where much of the precipitation has a lighter intensity. GPM horizontal resolution is 5 km × 5 km and Ku-band vertical resolution is 250 m. In this study, GPM radar PFs are obtained from 2014 to 2017. The GPM PF database uses an algorithm similar to the TRMM PF database originally developed at the University of Utah [40]. The PFs are defined by grouping contiguous areas with nonzero near-surface precipitation using the Ku-band radar [41], and the location (lat/lon) of the PF is the centroid of a best-fit ellipse. This analysis method condenses the original pixel-level measurements into the properties of events, which greatly increases the efficiency of searching and sorting the observed historical events [40]. Maximum echo-top heights for calculating PF depth are obtained using the GPM KuPR, which has a minimum detectable reflectivity near 12 dBZ [42]. For some PFs with large areas, the PF centroid and maximum echo-top height location could be quite different (most of these PFs are elongated fronts with continuous precipitation). Therefore, PFs with >2° lat/lon difference between the centroid and maximum echo-top height are excluded from this study, which removes 1132 PFs (0.16%) prior to collocation.

2.2. GPS Radio Occultation Data

GPS-RO soundings are obtained from three missions for this study: the joint US-Taiwan six-satellite FORMOSA Satellite Series No. 3/Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) mission [38], the German TerraSAR-X satellite [43], and the GRACE-B satellite [44]. The COSMIC constellation provided over 1000 soundings per day with relatively homogeneous sampling coverage around the globe when GPM was launched in early 2014. However, the number of daily soundings decreased to roughly 250 per day at the end of 2017. To increase the sampling numbers, additional GPS-RO observations from the TerraSAR-X and GRACE-B satellites are obtained. Both satellites provided ~250 soundings/day and comprise ~25% of the total collocations with GPM. The reprocessed level-2 RO profiles for all three missions are obtained from the COSMIC Data Analysis and Archive Center (CDAAC) at the University Corporation for Atmospheric Research (UCAR). The profiles are quality controlled by excluding the ones with “bad” flags (such as if the observation bending angles exceed the climatology by a specific threshold). We use the “atmPrf” product, which provides refractivity and dry temperature (Tdry) from usually near the surface up to ~60 km. At microwave frequencies, the atmospheric refractivity N is related to atmospheric pressure P, temperature T, and water vapor partial pressure e [45]:
N = k 1 P T + k 2 e T 2
where k1 is 77.6 K hPa−1 and k2 is 3.73 × 105 K2 hPa−1. The dry temperature is derived from the refractivity Equation (1) by neglecting atmospheric moisture [46], the second term within (1), such that:
T d r y = k 1 P d r y N
where Pdry is the dry pressure derived through hydrostatic integration by again neglecting the moisture within (1). In this study, dry temperature is used because it can be treated as an independent satellite retrieval, whereas the real temperature retrieval (“wetPrf”) relies on a priori moisture information from the European Centre for Medium-Range Weather Forecasts (ECMWF) low-resolution analysis. Note, however, the dry temperature retrieval is nearly identical to the real temperature in the UTLS region, as moisture is generally negligible when temperatures are lower than 250 K [47]. Thus, only sounding altitudes with a Tdry less than 250 K were used in this study to minimize the impact of water vapor on the temperature [24]. Therefore, the RO dry temperature in the UTLS region is denoted as RO temperature as it represents the real atmospheric temperature. The retrieved profiles are reported as a function of geometric height above mean sea level and the location of each profile nearest to the surface is used. The vertical resolution of RO soundings varies from 0.2 km in the lower troposphere to 1.4 km in the upper stratosphere [38] with an average of 0.5 km in the UTLS [48], while the horizontal footprint of tangent point trajectories is ~200 km [38]. Consistency and structural uncertainty of RO temperatures are similar between the various missions and uncertainties have been demonstrated to be smallest in the UTLS region between 8–25 km [49].
The monthly zonal mean (using 2.5° grids) global temperatures, temperature standard deviation, and LRT height are shown in Figure 1 using COSMIC profiles from 2006–2017. Throughout the tropics, LRT heights are consistently high year-round (~15 to 17 km) and UTLS temperature standard deviation remains low. A transition region occurs in the midlatitudes around 30° along with a sharp “tropopause break” near the subtropical jet, when the tropical tropopause extends to higher latitudes and overlies the extratropical lower stratospheric air. This transition region migrates with season. In the wintertime, it moves equatorward, becomes sharper and, consequently, the temperature standard deviation increases, whereas in the summer, it moves poleward and the tropopause height decreases more gradually. Several temperature standard deviation “hotspots” are evident throughout the UTLS, such as during the Northern Hemisphere (NH) summer monsoon (JJA) around ~30° N when deep convection frequently occurs. In the high latitudes, the tropopause remains low year-round (~8 to 12 km) with a large temperature standard deviation that reaches a maximum during the winter. The extremely large standard deviation throughout the polar lower stratosphere is due to variability along the edges of the polar vortex. The Arctic vortex reaches its highest variability during the NH winter (DJF), while the Antarctic vortex displays its largest variability in the Southern Hemisphere (SH) spring (SON) since it is stronger, larger, and longer-lasting [50]. The large standard deviation of UTLS temperatures observed throughout the mid-to-high latitudes is of particular interest in this study. In the tropics, diabatic heating from frequent deep convection is the dominant process, yet temperature standard deviation remains relatively small year-round. However, in the extratropics, both moist and dry dynamics play a role in the development of precipitation and contribute to the large temperature standard deviation. Thus, proper characterization of extratropical environments is paramount in accurately determining how PFs relate to their environment and impact UTLS temperature variation.

2.3. ERA-Interim Data

ERA-Interim is a global reanalysis produced by the ECMWF. This gridded dataset covers 1979–2019 in six-hour intervals [51]. Pressure-level potential vorticity data are used, which have a horizontal resolution of 0.75° latitude × 0.75° longitude with 37 vertical levels from 1000 to 1 hPa. Roughly six to seven levels are provided within the UTLS.

3. Methodology

To analyze the relationship between PFs and the thermodynamic structure of the UTLS, GPS-RO temperature profiles are collocated with PFs observed by GPM from 2014 to 2017 throughout the extratropics (20–65°). The collocated GPS profile must be within 3 h after PF occurrence and within a 300 km radius of the PF. Note that the 3 h and 300 km collocation criteria are chosen to provide a good balance between robust sampling and observation of the direct impact to the PF’s surrounding environment. In addition, for each RO temperature profile, the nearest-neighbor ERA-I potential vorticity profile is obtained.
Throughout the extratropics, considerable seasonal and meridional variation of the local tropopause height is observed [52], and the top of the troposphere coincides closely to the level above which convection rarely penetrates. This is because the upper bound of tropospheric clouds is strongly constrained by radiative cooling from water vapor and the depth of tropospheric mixing, and these physical linkages have been applied to both tropical and extratropical convection [53]. Therefore, to account for the large variation of tropopause heights, we introduce the PF relative depth (RD), which measures the fractional PF depth (e.g., the maximum echo-top height) relative to the local LRT height:
RD = M a x i m u m   H e i g h t   o f   D e t e c t a b l e   R a d a r   E c h o e s L a p s e   R a t e   T r o p o p a u s e   H e i g h t
We use a relative definition of depth because the absolute difference does not take into account the “size” of the quantities involved (e.g., the absolute difference of an echo-top height of 4 km and LRT height of 9 km is the same as an echo-top height of 12 km and LRT height of 17 km). Since we focus on PFs that reach near/into the UTLS, only PFs with a maximum echo-top height of at least 50% of the LRT height (e.g., RD ≥ 0.5) are studied.

3.1. PF and GPS Temperature Profile Classification Using Potential Vorticity

Over the tropics, PFs typically result from convective instability. However, over the extratropics, PFs can develop not only from tropical-like convective instability, but also baroclinic instability. Since the environments these PFs develop in can be drastically different, it is necessary to devise a method of classifying the PFs that form in these two different environments. In the extratropics, potential vorticity (PV) is a tracer-like variable based on the conservation of thermodynamic properties and momentum. It has long been used to identify the dynamical tropopause, as the stratosphere is characterized by significantly higher values of PV [9,54,55]. The extratropical tropopause is remarkably close to the ±2-PVU surface, where PVU denotes the standard potential vorticity unit (1 PVU = 10−6 m2 s−1 K kg−1) [5]. PV on various isentropic surfaces has often been used to identify the intrusion of stratospheric air into the troposphere. Namely, the intersection of the 320 K isentrope to the 2-PVU surface has been used to identify moderate to deep stratospheric intrusions or tropopause folds, which frequently occur near extratropical cyclones during a strong frontal passage from late fall through early spring [25,56,57,58]. It has also been suggested that the intrusion of stratospheric air may impact precipitation development, as the convective environment is often enhanced when drier stratospheric air descends to the middle levels [59,60].
In this study, we use PV on the 320 K isentrope to classify the observed PFs into two categories where the trigger for PF development is similar in nature within the environment (baroclinic instability vs. convective instability). All GPS profiles with PV of at least 2-PVU on the 320 K isentrope (identified from collocated ERA-I PV) are labeled DSI as the associated PFs coincide with strong dynamic effects on the tropopause through the aforementioned stratospheric intrusions. Most of the PFs within this classification occur due to large baroclinic instability and are found within extratropical cyclones, cut-off lows, or along strong frontal boundaries. These PFs occur most frequently in fall, winter, and spring throughout the mid/high latitudes. On the other hand, profiles with PV of less than 2-PVU at 320 K are labeled non-DSI as the associated PFs are more likely to be related to convective instability. The PFs within this group mainly include a variety of weather phenomena that occur most frequently from mid-spring through early fall, such as single cell thunderstorms, supercells, squall lines, mesoscale convective systems, and tropical cyclones. After classifying the GPS profiles, non-DSI and DSI median monthly background temperatures/tropopause heights are derived using profiles within 2.5° latitude × 5° longitude grids. Note that RO profiles from 2006–2017 are used to derive the background profiles and the grid sizes were chosen to allow for enough RO profiles in each grid cell. Additionally, a minimum of 30 profiles are required in a grid cell for deriving the median background profile for either environment, and the median statistical method was chosen to reduce the impact of outliers on the background profiles. Since we use monthly gridding, the presence of both environments can exist within one grid during that month (e.g., as different synoptic patterns propagate over the midlatitudes during winter). To reduce the uncertainty of the PF classification solely based on PV values, the simple least-squares method is used to determine the PF’s classification by computing the median absolute deviation (MAD) of the PF temperature profile and the two median background profiles within ±5 km around the PF LRT height. The MAD was used in this scenario since it is more resistant to outliers than the standard deviation. Then, the background profile with a smaller MAD is the classification assigned to the PF, which makes the PF classification more robust.
Figure 2 displays the seasonal GPS-RO sampling for both non-DSI (blue contours) and DSI (red contours) background profiles along with gridded median LRT height differences (LRTnon-DSI–LRTDSI, color shaded) for the transition regions where both types of PFs coexist. The number of GPS profiles in each grid generally ranges from 60 to 120 in areas with one type of PF, but can be significantly reduced in areas where both environments occur (e.g., the midlatitudes). The transition region shifts equatorward as winter approaches and moves poleward in the summer. The corresponding LRT height differences between the median non-DSI and DSI backgrounds also shows large seasonal variation, especially in the NH due to larger land/ocean contrast. The largest differences are typically observed over the oceans. In January (Figure 2a), height differences of 2.5 to 3.5 km are seen over the Kuroshio and Gulf Stream warm ocean currents. This can be attributed to consistent changes in air mass type and large baroclinicity along the coasts. On the other hand, the smallest variability occurs over the continental interior as conditions are frequently cold and dry. The transition region shifts much farther north in July (Figure 2c), with few DSI profiles and small LRT height differences. In general, the transition region coincides well with the latitude band showing the largest UTLS temperature variation (between 40–60° from Figure 1), which further confirms the importance of separating the two distinct PF categories using PV for an extratropical UTLS study.

3.2. Additional Quality Control

In addition to the large tropopause height standard deviation within 40–60°, there can be large intraseasonal LRT height variations throughout the subtropics (~25–35°), especially during the winter to early spring months in both hemispheres where standard deviations often exceed 2 km. For example, studies have shown a bimodal distribution of subtropical tropopause heights, with one mode above 15 km (characteristic of the tropical tropopause) and another one below 13 km (typical of the extratropical tropopause) that has been shown to be related to the formation of multiple tropopauses [61,62]. Pan et al. [63] has also shown the occurrence of multiple tropopauses extending to nearly 60° N. Therefore, in order to minimize the impact that intraseasonal LRT height variation could have on monthly median GPS background temperature profiles and anomaly calculations, PFs that occur within any grids with monthly LRT height standard deviations exceeding 2 km are excluded from this study (which removes ~6.9% of PFs).
Note that there also may be a PF-type ambiguity within larger DSI PFs. For example, many large-size DSI PFs form along elongated cold fronts in the wintertime with a large meridional and small zonal extent. If the pixels identified by the GPM KuPR are contiguous along the length of the front, then the front would be considered a single PF. However, there may be some deeper convective features near the equatorward-end of the PF, which may be some distance away from the stratospheric intrusion. While steps are taken to minimize this limitation (such as filtering out PFs with a distance greater than 2° between the PF center and the location of PF maximum height, see Section 2.1), some stronger convective features with more non-DSI characteristics may remain within large DSI PFs.

3.3. UTLS Temperature Anomaly and Tropopause Height/Temperature Anomaly Calculations

In this study, each individual PF temperature anomaly profile (T′) is derived by subtracting the respective non-DSI or DSI gridded median GPS background temperature (TBG) profile from the GPS PF temperature profile (TPF):
T ( z ) = T PF ( z ) T B G ( z )
where height z is from within ±5 km of the PF’s LRT height to account for the large seasonal and meridional variations in tropopause height. Figure 3 shows two examples of GPS temperature profiles collocated with PFs (red) and their gridded median non-DSI and DSI background profiles (black). Figure 3a shows a relatively large-size (20,229 km2) overshooting non-DSI PF (RD = 1.07). Overshooting (or tropopause-penetrating) PFs are generally considered some of the most intense convective systems on Earth and have attracted much attention with regards to STE [64]. The temperature anomaly profile (Figure 3c) is calculated by subtracting the non-DSI background profile (solid black line) from the PF profile. The temperature anomalies for each PF are centered on the PF’s LRT height and obtained for a ±5 km window around this height. Then, the anomalies are grouped according to the different PF characteristics (PF size and RD) and averaged to generate a mean anomaly profile for each group. The classification of the PF temperature profile from Figure 3a is relatively straightforward. However, there are some instances (e.g., Figure 3b) that can be somewhat ambiguous. While the upper tropospheric section of the PF profile is certainly “non-DSI” in nature, the tropopause/lower stratosphere is rather uncertain. Thus, the temperature anomaly profile from Figure 3d looks quite different than Figure 3c due to the PF classification uncertainty.
Similarly, LRT height/temperature anomalies are also computed by subtracting the median background from the PF tropopause height/temperature:
LRT H = LRT H PF LRT H B G
LRT T = LRT T PF LRT T B G
The PFs are again categorized by their area and RD to determine the relationship between PF characteristics and anomaly magnitude.
It is worth noting the challenge in classifying PFs in the extratropics for a composite study such as this. For example, when using PV, any potential temperature surface chosen (e.g., 320 K) will intersect the tropopause at different latitudes throughout the year [25]. However, the joint distribution of LRT height/temperature anomalies derived from collocated GPS profiles and PV values for extratropical PFs clearly shows two distinct PF categories separated by the ±2-PVU threshold at 320 K (Figure 4). Note that positive PVU values indicate PFs within the NH and negative PVU values indicate PFs within the SH. Non-DSI PFs have most samples clustered to minor LRT anomalies, although a few LRT height/temperature anomalies can reach ±5 km/±20 K. However, there is a distinct change in the anomalies observed around ±2-PVU. The majority of DSI events display negative LRT height anomalies along with a corresponding LRT temperature increase. The anomalies are generally larger in the SH than the NH. Additionally, the anomaly magnitudes tend to become larger as PVU values also become larger. These clear differences in tropopause anomalies support how the PV method can successfully separate the PFs into two general groups that have distinct opposite influences on the tropopause.

4. Sampling and PF Characteristics

The extratropics in both the NH and SH (20–65°) are chosen for analysis in this paper. There are 832,858 PFs observed by GPM throughout the extratropics with 10,736 non-DSI PFs and 11,239 DSI PFs with an RD > 0.5 collocated with a GPS-RO temperature profile. Figure 5 displays the distribution of GPS-collocated PFs with a RD > 0.5 for both non-DSI PFs (a,b) and DSI PFs (c,d) categorized by size and RD. Larger-size non-DSI PFs (over 50,000 km2) are most frequently observed over the oceanic regions (Figure 5a). These larger systems are a combination of mesoscale convective complexes, extratropical cyclones that have progressed from mature and into the occlusion stage of development, and tropical cyclones. In contrast, the deepest (overshooting, RD > 1.0) PFs are mainly found over land, such as North America, Europe, and Argentina (Figure 5b). These deepest PFs are primarily summertime features such as squall lines and supercells fueled by warm land surface temperatures and large convective instability. However, some deep PFs also form over warm ocean currents such as the Gulf Stream, some of which are likely tropical cyclones. There are also many shallower non-DSI PFs over both land and ocean, which are mostly comprised of single-cell thunderstorms and warm-frontal precipitation. Also, note that there are few non-DSI PFs in the subtropics off the continental west coasts due to the general lack of precipitation in these regions. On the other hand, DSI PFs (Figure 5c,d) are most often observed in the mid-and-high latitudes, with few occurrences in the subtropics. Most DSI PFs occur over the oceans due to contrasting air masses (baroclinic instability), with the cold and dry continental air masses interacting with the warmer and wetter oceanic air masses. The majority of these PFs include precipitation associated with extratropical cyclones, along frontal boundaries, or within cutoff lows. Several hotspots stand out, such as the wintertime storm tracks over the Southern Ocean, the Icelandic Low, and the Aleutian Low. No consistent pattern is shown for different PF size or RD categories. However, as expected, many more DSI PFs have large sizes, but are rather shallow since extratropical cyclones typically have much weaker uplift than summertime convection.
Additionally, Figure 6 shows an example of an overshooting convective system (RD = 1.13) observed by GPM in the early morning of 28 July 2015 over South Dakota in the United States. This system is a strong bow-echo, with near-surface reflectivity reaching up to 60 dBZ (Figure 6a) and maximum echo-top heights near 17 km (Figure 6b). A cross-section of the PF is also displayed, with very low brightness temperatures and high rain rates (50 mm h−1) occurring near the core of the PF (Figure 6e). The collocated ERA-I PV is near or slightly above 1-PVU, indicating a non-DSI PF (Figure 6c). Finally, the GPS temperature profile near the maximum echo-top height (44.2° N, 101° W) is shown along with its corresponding background profile, and the LRT height and maximum echo-top height are also displayed (Figure 6d). Strong upper tropospheric warming along with cooling around the tropopause are observed for this overshooting PF.

5. Results

5.1. Relationship of PF Characteristics to UTLS Temperature Anomalies

This section focuses on extratropical PF characteristics and UTLS temperature anomalies within ±5 km of the PF LRT height. GPS-RO temperature anomalies near both non-DSI and DSI PFs are separated by their surface properties (land or ocean) and further divided into subgroups based on PF size (area) and depth (RD). For example, PFs reaching close to the LRT are grouped together within the “0.8–1.0” RD classification, while PFs reaching above the LRT (overshooting) are grouped within the “>1.0” classification. Finally, PFs are grouped by their month of occurrence to quantify any seasonal differences in anomalies.

5.1.1. Non-DSI PFs

Figure 7 displays the mean GPS-RO temperature anomalies near non-DSI PFs over the extratropics (20–65°). For all PFs classified by various RDs (Figure 7a), warm anomalies are observed for the shallower PFs within roughly 1.5 to 5 km below the LRT, ranging from 0.1 to 0.8 K. Interestingly, the strongest warming (~0.8 K) occurs for the shallower PFs (RD < 0.6), whereas either slight cooling or no anomalies are observed for deeper PFs reaching near or above the tropopause (RD > 1). We speculate that this is because the more intense deep convection typically occurs in the earlier part of a PF’s life cycle [65]. As a result, the cumulative influence of latent heat release that occurs around the convective core may not be sufficient to affect the surrounding environment yet [66,67]. Above, the upper tropospheric warming, a strong layer of cool anomalies centered on the PF LRT is observed, ranging from roughly −2.5 to −3 K. Both the shallowest and deepest PFs show comparable magnitudes of strong cooling. Above the cooling layer, warming of 0.5 to 1 K is generally observed with the largest magnitudes for the deepest overshooting PFs. There are major differences in anomaly magnitude between land and oceanic PFs (Figure 7b,c), with much more variation evident for oceanic PFs. For example, oceanic PFs display a wide range of anomalies below the tropopause layer, with moderate warming for the shallower RDs and strong cooling for the deeper RDs. However, warming is observed for all land PF RD categories. Additionally, the largest anomaly magnitudes occur near oceanic PFs, as the strongest tropopause-level cooling, as well as warming below and above the tropopause, is observed in these groups.
Similar temperature anomaly patterns are seen for all PFs with different sizes (Figure 7d). An increase in anomaly magnitude from smaller to larger PF size is observed throughout the vertical extent of the profiles. A warming layer is displayed up to 1.5 km below the LRT for all subgroups, with the largest PFs exhibiting the strongest warming (up to 0.8 K). The strongest cooling at the LRT also occurs for the largest PFs (roughly −3 K). Above the cooling layer, a much smaller range of anomalies is displayed, with minor warming (~0.3 K) observed for all PFs. Differences between land and ocean PF size subgroups are relatively minor (Figure 7e,f), highlighted by slightly stronger warming in both the upper troposphere and lower stratosphere for oceanic PFs. Larger anomaly variation is again observed among different size oceanic PFs, but with much smaller variation compared to the RD dependency.
The UTLS temperature anomalies near extratropical non-DSI PFs show similar patterns to those near tropical deep convection [24], albeit with slightly larger magnitudes. The warm anomalies observed throughout the mid-to-upper troposphere result from latent heat release in-and-near the convective clouds [1,18]. Additionally, some shallower non-DSI PFs are associated with warm frontal precipitation in the extratropics, and latent heat also is released during the formation of different hydrometeor species induced by warm conveyor belts [68,69]. The cool anomalies that occur near the convective tropopause have been characterized in a few different ways. Holloway and Neelin [19] showed that as gravity waves spread convective warming through the free troposphere, hydrostatic pressure gradients will extend above the heating, causing divergence, ascent, and adiabatic cooling aloft. The convective “cold-top” should be thought of as an inherent part of quasi-equilibrium temperature adjustment. Khaykin et al. [21] also discussed two cooling mechanisms: the first being nonmigrating tides generated by convective diabatic heating; and the second being the systematic injection and turbulent mixing of adiabatically cooled air by cross-tropopause updrafts. Some previous studies have shown that warming occurs in the lower stratosphere, which is also observed in this study. For example, Chae et al. [70] observed warm anomalies at these altitudes and downward motions above the convection. They suggested that strongly divergent flow and turbulent mixing near the cloud tops would mechanically drag the air just above the cloud outward, which would pull comparatively warmer lower stratospheric air down from above the clouds. It has also been shown that fine-scale features in the lower stratospheric static stability field are reminiscent of the planetary wave response to the lower stratospheric diabatic heating associated with convection [71].
Figure 8 displays the monthly mean temperature anomalies for non-DSI PFs within the extratropics in both hemispheres. The extratropics are separated into two bands to determine any latitudinal differences (20–40° and 40–65°). Throughout each latitude band, large differences in the pattern and magnitude of anomalies are observed as the seasons change. Within 20–40° N (Figure 8a), the magnitude of upper tropospheric warming remains similar in each month, although the sharpness of the transition from warm to cool anomalies differs. However, tropopause-level cooling is weakest from late spring through early fall and strongest during the wintertime. The magnitude of summertime tropopause-level cooling (between 1.5–2 K) is similar to results shown in previous research throughout the tropics [24]. Considerable seasonal variation is displayed for the lower stratospheric warming. Similar results are also observed within 20–40° S (Figure 8c). For example, the weakest tropopause-level cooling is again observed during the summer, whereas the strongest occurs during the winter. In contrast, seasonal magnitude changes generally show the opposite patterns in the high latitudes for both hemispheres. Within 40–65° (Figure 8b,d), the strongest tropopause-level cooling is seen during the summer whereas the weakest occurs during the winter. Upper tropospheric warming is slightly stronger in the high latitudes, with the strongest warming occurring during austral summer and boreal winter. We attribute the seasonal magnitude changes in the two latitude bands to seasonal variations in PF characteristics (not shown). Within 20–40° N, the wintertime tends to have more frequent shallow/larger-size PFs relative to the summertime, when more frequent deep convective cores with smaller sizes are observed. However, within 40–65° N, the shallower convection that occurs during the summer has larger areas than the wintertime PFs of similar depth. The SH PFs also display similar characteristics albeit to a smaller degree as the majority of PFs occur over the ocean. This confirms the results from Figure 7, where PFs with larger sizes produced the strongest tropopause-level cooling.

5.1.2. DSI PFs

Figure 9 displays the mean GPS-RO temperature anomalies near DSI PFs within the extratropics (20–65°), with major differences in the anomaly vertical structure observed compared to non-DSI PFs. Near PFs classified by various RDs (Figure 9a), cooling is observed from 5 km below the LRT to roughly 0.5–1 km above the LRT. The cooling noticed throughout the upper troposphere cannot be interpreted by latent heat release from precipitation systems, as it is likely dominantly driven by larger scale dynamic mixing of baroclinic air masses. Above the LRT, temperature anomalies rapidly transition to strong warming. Maximum warming increases with larger RD (from 2 to 4 K) and occurs at progressively higher altitudes. This warming is largely contributed to by descent in the stratosphere due to large scale mixing associated with extratropical cyclones. Major temperature anomaly differences are evident between land and ocean PFs (Figure 9b,c). Throughout the extratropics, DSI events occur more frequently over the oceans due to the increase in baroclinicity over the warm ocean currents. This stronger baroclinicity is reflected in anomaly magnitudes that are much larger near the oceanic DSI events. On the other hand, the weaker baroclinicity for land DSI events means smaller temperature anomalies are observed. Differences in anomaly structure between land and oceanic PFs with different sizes (Figure 9e,f) are less pronounced, with oceanic PFs again displaying larger anomaly variation between different sizes.
It is worth noting the UTLS temperature anomalies (and tropopause anomalies shown later) observed near DSI PFs are strongly influenced by extratropical dynamics rather than just the PFs themselves (unlike non-DSI PFs). Extratropical cyclones obtain their energy from sharp horizontal temperature and moisture differences between tropical and polar air (e.g., highly baroclinic environments). We speculate that the differences in the magnitude of temperature anomaly variation as DSI PF area and RD increase are related to the maturing and dissipation process in these extratropical cyclones. During the early stages of extratropical cyclone development, surface frontogenesis is accompanied by strong differential temperature advection of cold and warm air. The PFs are often driven by large thermodynamic instability, with a stronger intensity, reaching higher altitudes, and having relatively smaller sizes. Thus, UTLS temperature anomalies are largest at this stage due to the stronger baroclinicity. As the cyclone matures, the system becomes occluded, thereby growing in size, but weakening in strength. The continued dynamic mixing among the two contrasting air masses would then be displayed as weaker cooling in the upper troposphere (Figure 9a,d). After the cyclone has mixed away the temperature contrasts across the front, there is no longer any potential energy available to keep the cyclone going so it dissipates. On the other hand, the larger, more mature systems have ample time to spread latent heat release to the surrounding upper tropospheric environment (evidenced in the upper tropospheric warming displayed only near the largest PFs). However, these speculations would need further validation through additional observations and model simulations in the future. Thermal advection also plays a role in the observed temperature anomalies for DSI events, but quantifying this role would be far from straightforward. For example, strong frontal boundaries in highly baroclinic environments have a mixture of warm and cold air advection. These boundaries would need to be analyzed separately, their direction of motion verified, and the location of the GPS-RO profile determined relative to this motion (as the collocation can be within a 300 km radius of the PF centroid) to better appreciate how thermal advection impacts the results.
DSI cases are associated with tropopause folding, where deep frontogenetic forcing can occur commonly near and upstream of precipitation events in baroclinic environments [72]. A tropopause fold is an extrusion of stratospheric air within an upper-tropospheric baroclinic zone that slopes downward from the normal tropopause level to the middle or lower troposphere [73]. The fold is associated with a substantial lowering of tropopause heights and has been identified to be a contributor to upper-level frontogenesis and the rapid development of surface storms [72,74]. In summary, the DSI PFs themselves are not the main cause of the UTLS temperature anomalies observed and are more directly related to the nearby dynamical processes.
Figure 10 displays the monthly mean temperature anomalies for DSI PFs within 40–65° in both hemispheres. For DSI events, even larger seasonal differences are observed in the pattern and magnitude of anomalies as the seasons change. In the NH, wintertime events display cool anomalies throughout the entirety of the upper troposphere, with maximum cooling of −4 K occurring near the tropopause. Additionally, lower stratospheric warming is also at a maximum in the winter (~3 K). In contrast, summertime events display weak warming in the upper troposphere (~0.5 K) along with relatively weaker cooling (less than −2 K) near the tropopause and warming (~1 K) in the lower stratosphere. Anomaly magnitudes during the transition seasons generally remain in between summer and winter magnitudes. Similar results are also shown for the SH. Cool anomalies are observed throughout the entirety of the upper troposphere in all seasons, but with stronger summertime cooling near the tropopause (−2 to −3 K). The larger-amplitude wintertime anomalies are likely due to stronger temperature contrasts fueling the development of intense midlatitude cyclones, which cause deeper tropopause folding and significant stratosphere-troposphere exchange. These contrasts weaken significantly in the summer, leading to reduced dynamic UTLS mixing, which results in anomaly profiles that look more similar to non-DSI PFs. Additionally, the minor differences in anomaly structure displayed between the NH and SH are in contrast to previous results from Biondi et al. [75], which observed considerable structural differences in anomalies between NH and SH extratropical cyclones.

5.2. Relationship of PFs to Tropopause Height/Temperature Anomalies

In this section, we focus on the changes that occur specifically to the thermal tropopause near PFs. Two-dimensional histograms are constructed relating the RD and size of the PF to the magnitude of tropopause height and temperature anomalies observed. The PFs are binned using intervals of 0.1 for RD and bin sizes increase logarithmically from 100 km2 to ~120,000 km2 for area. Contours of the number of PFs within each bin are also plotted (excluding the bins with fewer than ten profiles).
Figure 11 displays the LRT height and temperature anomalies for non-DSI and DSI PFs within the extratropics (20−65°). For non-DSI PFs, moderate LRT height decreases are observed (0.1 to 0.4 km), on average, for PFs with a high RD (Figure 11a). In contrast, considerable LRT height increases (0.2 to 0.6 km) are observed for PFs with a low RD, and magnitudes increase as PF size increases. The corresponding temperature anomalies (Figure 11c) show a similar pattern, but with the opposite sign compared to the height anomalies. On the other hand, DSI events are consistently associated with negative LRT height anomalies (Figure 11b) along with corresponding positive LRT temperature anomalies (Figure 11d). Much larger anomaly magnitudes are observed relative to non-DSI PFs. Anomaly magnitude is larger near PFs with a stronger RD, which likely occurs in the earlier stages of cyclone development. The largest LRT height anomalies (−1.5 to −2.0 km) and temperature anomalies (6 to 8 K) are observed near the highest RDs, whereas the smallest LRT height anomalies (−0.25 to −0.5 km) and temperature anomalies (0.1 to 2 K) are observed near the largest PFs with a shallow RD, when cyclones are likely in the occlusion stage.
Figure 12 displays the observed LRT temperature/height anomalies (solid contours) along with the temperature anomalies at the PF LRT relative to the background profile at the same altitude (the anomalies at the “zero line” from Figure 7 and Figure 9; shaded contours), as well as individual cases representing typical RO profiles near each type of PF. Non-DSI PFs throughout the extratropics (Figure 12a) display a wide range of tropopause anomalies, with LRT temperature anomalies up to ±12 K and height anomalies up to ±2.5 km. However, most of these samples are clustered near the middle with only a slight trend toward the lower-right quadrant. Most PFs are typically found in the upper-left and lower-right quadrants, which indicates that the tropopause is most commonly raised/cooled or lowered/warmed. Even though the PF tropopause displays considerable variation relative to the background tropopause, strong cooling (shaded contours) is almost always observed for the PF LRT relative to the background temperature at the same altitude (displayed in Figure 7). The strongest cooling (between −8 to −10 K) occurs when the PF LRT is pushed much higher (e.g., 2 km), but notable cooling is also observed (between −2 to −6 K) even when the PF LRT is at a lower altitude than normal. This is explained by a commonly-observed type of profile seen in the lower-right quadrant and is shown in Figure 12c, where a PF-collocated GPS profile (red) is compared to its background profile (black). In this instance, the PF LRT is almost 1 km lower and 3 K warmer than the background LRT. However, note that the PF LRT is roughly 2 K colder than the background temperature at this altitude. Near the top of the PF (blue dotted line), the lapse rate increases compared to the background and a relative minimum temperature is observed above the convective cloud. Above this relative minimum, an inversion commonly occurs, and then the temperature progresses towards another relative minimum at a slightly higher altitude. This often results in a double tropopause structure, as a new, lower LRT forms near the cloud top and the second (higher) LRT remains either near or a few kilometers above the climatological LRT, similar to previous results shown by Biondi et al. [76]. While this type of profile was observed throughout the entire extratropics, it was most commonly seen in the subtropics. Additional research into the formation of these “lower” tropopauses near the convective top is recommended, especially for determining how long this type of UTLS environment persists after PF occurrence.
Similarly, DSI PFs throughout the extratropics (Figure 12b) display a wide range of tropopause anomalies. The amplitude of maximum anomalies nearly doubles for DSI events, with LRT temperature anomalies up to 20 K and LRT height anomalies near −4 km. Most samples are within the lower-right quadrant, indicating the PF LRT typically becomes much lower and warmer. Nevertheless, similar to non-DSI PFs, strong cooling is again observed at the PF LRT relative to the background temperature at the same altitude. However, the strongest cooling now occurs when the PF LRT height decreases the most. A typical GPS temperature profile for a high-RD DSI PF is displayed in Figure 12d. The PF temperatures are much colder than the background environment throughout the troposphere. The cold UTLS air descends deeply into the troposphere along a tropopause fold, which is associated with typical baroclinic features such as strong cold fronts and extratropical cyclones [31]. The PF LRT shows a much sharper temperature transition, with a height about 2 km lower (at ~7.5 km) and temperature ~3 K warmer than normal. The tropopause inversion layer, which is a region of enhanced static stability above the extratropical tropopause associated with a narrow-scale temperature inversion [77], occurs between ~7.5 to 10.5 km and becomes somewhat stronger. Finally, temperatures remain above normal throughout the lower stratosphere near the PF. This type of profile is the most common profile observed near DSI PFs, although there are some instances when the PF LRT can be pushed higher and become colder.

6. Discussion and Conclusions

In this paper, the relationship between extratropical precipitation and the thermodynamic structure of the UTLS in both hemispheres (20–65°) was studied by collocating precipitation features observed by the GPM precipitation radar with high vertical resolution GPS-RO temperature soundings from 2014 to 2017. PFs were classified into two different categories (non-DSI and DSI) using PV. UTLS temperature and tropopause height/temperature anomalies were calculated and evaluated according to PF characteristics (e.g., PF sizes, depths, surface types, and seasonality).
For non-DSI PFs, warm anomalies ranging from 0.1 to 1 K were observed from roughly 5 to 1.5 km below the LRT due to latent heat release in-and-near the convective clouds. This transitioned to a layer of strong cooling ranging from −2 to −3 K centered around the PF LRT height, which occurs due to the injection of adiabatically cooled air near the tops of the PFs. Beginning about 1 km above the PF LRT, minor-to-moderate warming ranging from 0.2 to 1 K was observed. Furthermore, a consistent progression in anomaly magnitude was observed for upper tropospheric warming as well as cooling around the PF LRT as PF size increased. In contrast, the RD displayed mixed results on anomaly magnitude, as the strongest warming in the upper troposphere and cooling near the PF LRT generally occurred for PFs with the shallowest RD. We speculate this may be related to the life cycle of a PF, as intense overshooting convection typically occurs in the earlier stages of convective development and enough time may not have passed to spread the warming and cooling to the environment surrounding the PF.
On the other hand, DSI events displayed distinctly different patterns as the UTLS temperature response was dominated by large-scale dynamics, with strong cooling ranging from −2 to −5 K (except near large-size PFs) observed throughout the mid/upper troposphere and maximum values occurring slightly below the PF LRT. The anomalies transitioned to moderate/strong warming above the LRT, ranging from 2 to 4 K. The strongest upper tropospheric cooling and lower stratospheric warming occurred near the deepest PFs, whereas the weakest magnitudes were observed near the largest size PFs, which we attributed to properties of PFs under the different stages of extratropical cyclone development. The temperature anomaly pattern observed during DSI events likely occurred due to strong horizontal dynamic mixing among contrasting tropical and polar air masses and contributes to the large temperature standard deviation observed throughout the upper troposphere within the mid/high latitudes. Additionally, larger magnitudes and more variation in the temperature anomalies was evident for oceanic DSI events, which was expected as baroclinicity is typically strongest near the coasts.
In addition, seasonal differences in anomalies for both PF categories were explored. For non-DSI PFs between 20–40° in both hemispheres, tropopause-level cooling was weakest from late spring through early fall and strongest during the wintertime. However, within the higher latitudes (40–65°), seasonal magnitude changes showed the opposite patterns, as the strongest tropopause-level cooling was seen during the summer, whereas the weakest occurred during the winter. This was attributed to latitudinal seasonal variations of PF characteristics. For DSI events within 40–65° in both hemispheres, wintertime PFs displayed cool anomalies throughout the entirety of the upper troposphere and the strongest near-tropopause cooling of any season. In contrast, minor warming was observed in the upper troposphere along with weaker near-tropopause cooling for summertime events. The larger-amplitude wintertime anomalies likely occurred due to stronger temperature contrasts fueling the development of intense midlatitude cyclones.
Finally, the relationship between extratropical PFs and LRT height/temperature was quantified relative to PF area and RD. For non-DSI PFs, minor decreases in LRT height (−0.1 to −0.4 km) and increases in LRT temperature (0.1 to 1.0 K) were generally observed for PFs with a deeper RD, whereas PFs with a shallower RD and large area displayed large LRT height increases (0.2 to 1 km) and temperature decreases (−0.1 to −2 K). Overall, non-DSI PFs throughout the extratropics displayed a wide range of tropopause anomalies, with LRT temperature anomalies up to ±12 K and height anomalies up to ±2.5 km. Two types of UTLS temperature profiles were commonly observed for non-DSI PFs: (1) the PFs either push the tropopause higher and it becomes colder; or (2) a double tropopause structure develops, with a new, lower tropopause forming near the top of the convective cloud beneath the higher climatological LRT. DSI events were associated with much larger tropopause anomalies. On average, LRT heights were 0.25 to 2 km lower and LRT temperatures were 2 to 8 K warmer than normal during significant deep stratospheric intrusion events, indicative of subsidence and tropopause folding. Additionally, LRT height decreases of almost −4 km and temperature increases near 20 K were observed. There was one common profile type observed near DSI PFs, as the LRT displayed a much stronger than normal temperature transition and a robust tropopause inversion layer was seen.
It is worth noting that there are some limitations in this study. First, GPM only provides a snapshot of the PF at the time of observation so it cannot capture the stage of PF development (although it can be inferred to some extent). This additional information would be useful for providing a more in-depth explanation of the differences in the temperature anomaly magnitudes under different PF sizes and depths. Second, UTLS temperature and tropopause anomalies would ideally be determined relative to the larger neighboring environment of each PF. Unfortunately, this is not possible with GPS-RO sounding density, so a composite study is unable to describe the full influence of convection on its environment. Third, the number of GPS-RO soundings available in the extratropics has steadily declined in recent years, as satellite missions such as COSMIC have reached the end of their life cycle. The recently launched COSMIC-2 mission only provides soundings within the tropics and subtropics, so additional GPS-RO satellites are needed to provide necessary coverage throughout the extratropics. However, despite these limitations, these results suggest intricate relationships between different types of precipitation systems and the UTLS temperature structure and tropopause throughout the extratropics. Importantly, non-DSI PF temperature anomaly patterns and magnitude agreed well with previous research on tropical convection, which we believe provides a “unification” of tropical and extratropical deep convection impacts to UTLS temperatures.
This study enhances understanding of extratropical precipitation and UTLS temperature responses, and it is hopeful that this will lead to the continued improvement of precipitation system representation in weather models. Additional studies on the role that PFs play in moisture transport in the extratropical UTLS would also be beneficial in understanding extratropical STE and will likely be a topic of future research.

Author Contributions

Conceptualization, B.R.J., F.X. and C.L.; Methodology, B.R.J., F.X. and C.L.; Software, B.R.J.; Validation, B.R.J.; Formal analysis, B.R.J., F.X. and C.L.; Investigation, B.R.J.; Resources, F.X.; Data curation, B.R.J.; Writing—original draft preparation, B.R.J.; Writing—review and editing, B.R.J., F.X. and C.L.; Visualization, B.R.J.; Supervision, F.X.; Funding acquisition, F.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation, grant number 2054356; and by the National Aeronautics and Space Administration, grant numbers NNX14AK17G and NNX15AQ17G.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All datasets used in this study are freely and openly available. Thanks to the Precipitation Processing System (PPS) team at NASA Goddard Space Flight Center, Greenbelt, MD, for GPM data processing assistance. GPM precipitation feature database is available at http://atmos.tamucc.edu/trmm/ (accessed on 10 January 2019). GPS-RO data was obtained from CDAAC at UCAR: https://data.cosmic.ucar.edu/gnss-ro/ (accessed on 10 January 2019). ERA-Interim data was obtained from ECMWF: https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=pl/ (accessed on 10 January 2019).

Acknowledgments

The authors would like to thank Bill Randel, Thomas Winning, and Kevin Nelson for helpful discussions on the project.

Conflicts of Interest

The authors declare no conflict of interest and the funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Gettelman, A.; Birner, T. Insights into Tropical Tropopause Layer Processes Using Global Models. J. Geophys. Res. 2007, 112, D23104. [Google Scholar] [CrossRef] [Green Version]
  2. Peevey, T.R.; Gille, J.C.; Homeyer, C.R.; Manney, G.L. The Double Tropopause and Its Dynamical Relationship to the Tropopause Inversion Layer in Storm Track Regions. J. Geophys. Res. 2014, 119, 10,194–10,212. [Google Scholar] [CrossRef] [Green Version]
  3. Birner, T. Fine-Scale Structure of the Extratropical Tropopause Region. J. Geophys. Res. 2006, 111, D04104. [Google Scholar] [CrossRef] [Green Version]
  4. Gettelman, A.; Pan, L.L.; Randel, W.J.; Hoor, P.; Birner, T.; Hegglin, M.I. The Extratropical Upper Troposphere and Lower Stratosphere. Rev. Geophys. 2011, 49, RG3003. [Google Scholar] [CrossRef] [Green Version]
  5. Holton, J.R.; Haynes, P.H.; Mcintyre, M.E.; Douglass, A.R.; Rood, R.B.; Pfister, L. Stratosphere-Troposphere Exchange. Rev. Geophys. 1995, 95, 403–439. [Google Scholar] [CrossRef]
  6. Haynes, P.; Scinocca, J.; Greenslade, M. Formation and Maintenance of the Extratropical Tropopause by Baroclinic Eddies. Geophys. Res. Lett. 2001, 28, 4179–4182. [Google Scholar] [CrossRef] [Green Version]
  7. Randel, W.J.; Seidel, D.J.; Pan, L.L. Observational Characteristics of Double Tropopauses. J. Geophys. Res. 2007, 112, D07309. [Google Scholar] [CrossRef] [Green Version]
  8. World Meteorological Organization (WMO). Atmospheric Ozone 1985; WMO: Geneva, Switzerland, 1986; pp. 1–1095. [Google Scholar]
  9. Hoskins, B.J.; McIntyre, M.E.; Robertson, A.W. On the Use and Significance of Isentropic Potential Vorticity Maps. Q. J. R. Meteorol. Soc. 1985, 111, 877–946. [Google Scholar] [CrossRef]
  10. Pan, L.L.; Randel, W.J.; Gary, B.L.; Mahoney, M.J.; Hintsa, E.J. Definitions and Sharpness of the Extratropical Tropopause: A Trace Gas Perspective. J. Geophys. Res. 2004, 109, D23103. [Google Scholar] [CrossRef]
  11. Santer, B.D.; Sausen, R.; Wigley, T.M.L.; Boyle, J.S.; AchutaRao, K.; Doutriaux, C.; Hansen, J.E.; Meehl, G.A.; Roeckner, E.; Ruedy, R.; et al. Behavior of Tropopause Height and Atmospheric Temperature in Models, Reanalyses, and Observations: Decadal Changes. J. Geophys. Res. Atmos. 2003, 108, 4002. [Google Scholar] [CrossRef] [Green Version]
  12. Grosvenor, D.P.; Choularton, T.W.; Coe, H.; Held, G. A Study of the Effect of Overshooting Deep Convection on the Water Content of the TTL and Lower Stratosphere from Cloud Resolving Model Simulations. Atmos. Chem. Phys. 2007, 7, 4977–5002. [Google Scholar] [CrossRef] [Green Version]
  13. Fueglistaler, S.; Dessler, A.E.; Dunkerton, T.J.; Folkins, I.; Fu, Q.; Ote, P.W. Tropical Tropopause Layer. Rev. Geophys. 2009, 47, RG1004. [Google Scholar] [CrossRef]
  14. Hassim, M.E.E.; Lane, T.P. A Model Study on the Influence of Overshooting Convection on TTL Water Vapour. Atmos. Chem. Phys. 2010, 10, 9833–9849. [Google Scholar] [CrossRef] [Green Version]
  15. Held, I.M. On the Height of the Tropopause and the Static Stability of the Troposphere. J. Atmos. Sci. 1982, 39, 412–417. [Google Scholar] [CrossRef] [Green Version]
  16. Sherwood, S.; Horinouchi, T.; Zeleznik, H. Convective Impact on Temperatures Observed near the Tropical Tropopause. J. Atmos. Sci. 2003, 60, 1847–1857. [Google Scholar] [CrossRef] [Green Version]
  17. Kim, H.; Dessler, A.E. Observations of Convective Cooling in the Tropical Tropopause Layer in AIRS Data. Atmos. Chem. Phys. Discuss. 2004, 4, 7615–7629. [Google Scholar] [CrossRef] [Green Version]
  18. Randel, W.J.; Wu, F.; Rios, W. Thermal Variability of the Tropical Tropopause Region Derived from GPS/MET Observations. J. Geophys. Res. 2003, 108, 4024. [Google Scholar] [CrossRef]
  19. Holloway, C.; Neelin, J.D. The Convective Cold Top and Quasi Equilibrium. J. Atmos. Sci. 2007, 64, 1467–1488. [Google Scholar] [CrossRef]
  20. Paulik, L.C.; Birner, T. Quantifying the Deep Convective Temperature Signal within the Tropical Tropopause Layer (TTL). Atmos. Chem. Phys. 2012, 12, 12183–12195. [Google Scholar] [CrossRef] [Green Version]
  21. Khaykin, S.M.; Pommereau, J.-P.; Hauchecorne, a. Impact of Land Convection on Temperature Diurnal Variation in the Tropical Lower Stratosphere Inferred from COSMIC GPS Radio Occultations. Atmos. Chem. Phys. 2013, 13, 6391–6402. [Google Scholar] [CrossRef] [Green Version]
  22. Xian, T.; Fu, Y. Characteristics of Tropopause-Penetrating Convection Determined by TRMM and COSMIC GPS Radio Occultation Measurements. J. Geophys. Res. Atmos. 2015, 120, 7006–7024. [Google Scholar] [CrossRef] [Green Version]
  23. Kim, J.; Randel, W.J.; Birner, T. Convectively Driven Tropopause-Level Cooling and Its Influences on Stratospheric Moisture. J. Geophys. Res. Atmos. 2018, 123, 590–606. [Google Scholar] [CrossRef]
  24. Johnston, B.R.; Xie, F.; Liu, C. The Effects of Deep Convection on Regional Temperature Structure in the Tropical Upper Troposphere and Lower Stratosphere. J. Geophys. Res. Atmos. 2018, 123, 1585–1603. [Google Scholar] [CrossRef]
  25. Sprenger, M.; Wernli, H.; Bourqui, M. Identification and ERA-15 Climatology of Potential Vorticity Streamers and Cutoffs near the Extratropical Tropopause. J. Atmos. Sci. 2007, 64, 1569–1586. [Google Scholar] [CrossRef]
  26. Danielsen, E.F. Stratospheric-Tropospheric Exchange Based on Radioactivity, Ozone and Potential Vorticity. J. Atmos. Sci. 1968, 25, 502–518. [Google Scholar] [CrossRef]
  27. Baray, J.L.; Daniel, V.; Ancellet, G.; Legras, B. Planetary-Scale Tropopause Folds in the Southern Subtropics. Geophys. Res. Lett. 2000, 27, 353–356. [Google Scholar] [CrossRef] [Green Version]
  28. Price, J.D.; Vaughan, G. The Potential for Stratosphere-Troposphere Exchange in Cut-off Low Systems. Q. J. R. Meteorol. Soc. 1993, 119, 343–365. [Google Scholar] [CrossRef]
  29. Lamarque, J.-F.; Langford, A.O.; Proffitt, M.H. Cross-Tropopause Mixing of Ozone through Gravity Wave Breaking: Observation and Modeling. J. Geophys. Res. Atmos. 1996, 101, 22969–22976. [Google Scholar] [CrossRef]
  30. Stenchikov, G.; Dickerson, R.; Pickering, K.; Ellis, W.; Doddridge, B.; Kondragunta, S.; Poulida, O.; Scala, J.; Tao, W.K. Stratosphere-Troposphere Exchange in a Midlatitude Mesoscale Convective Complex: 2. Numerical Simulations. J. Geophys. Res. Atmos. 1996, 101, 6837–6851. [Google Scholar] [CrossRef]
  31. Reid, H.J.; Vaughan, G. Convective Mixing in a Tropopause Fold. Q. J. R. Meteorol. Soc. 2004, 130, 1195–1212. [Google Scholar] [CrossRef] [Green Version]
  32. Vaughan, G.; Price, J.D.; Howells, A. Transport into the Troposphere in a Tropopause Fold. Q. J. R. Meteorol. Soc. 1994, 120, 1085–1103. [Google Scholar] [CrossRef]
  33. Gray, S.L. A Case Study of Stratosphere to Troposphere Transport: The Role of Convective Transport and the Sensitivity to Model Resolution. J. Geophys. Res. 2003, 108, 4590. [Google Scholar] [CrossRef]
  34. Homeyer, C.R.; Pan, L.L.; Dorsi, S.W.; Avallone, L.M.; Weinheimer, A.J.; O’Brien, A.S.; DiGangi, J.P.; Zondlo, M.A.; Ryerson, T.B.; Diskin, G.S.; et al. Convective Transport of Water Vapor into the Lower Stratosphere Observed during Double Tropopause Events. J. Geophys. Res. Atmos. 2014, 119, 10941–10958. [Google Scholar] [CrossRef]
  35. Hou, A.Y.; Kakar, R.K.; Neeck, S.; Azarbarzin, A.A.; Kummerow, C.D.; Kojima, M.; Oki, R.; Nakamura, K.; Iguchi, T. The Global Precipitation Measurement Mission. Bull. Am. Meteorol. Soc. 2014, 95, 701–722. [Google Scholar] [CrossRef]
  36. Liu, C.; Zipser, E.J. The Global Distribution of Largest, Deepest, and Most Intense Precipitation Systems. Geophys. Res. Lett. 2015, 42, 3591–3595. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, N.; Liu, C. Global Distribution of Deep Convection Reaching Tropopause in 1 Year GPM Observations. J. Geophys. Res. Atmos. 2016, 121, 3824–3842. [Google Scholar] [CrossRef] [Green Version]
  38. Anthes, R.A.; Ector, D.; Hunt, D.C.; Kuo, Y.-H.; Rocken, C.; Schreiner, W.S.; Sokolovskiy, S.V.; Syndergaard, S.; Wee, T.-K.; Zeng, Z.; et al. The COSMIC/FORMOSAT-3 Mission: Early Results. Bull. Am. Meteorol. Soc. 2008, 89, 313–333. [Google Scholar] [CrossRef]
  39. Kummerow, C.; Barnes, W.; Kozu, T.; Shiue, J.; Simpson, J. The Tropical Rainfall Measuring Mission (TRMM) Sensor Package. J. Atmos. Ocean. Technol. 1998, 15, 809–818. [Google Scholar] [CrossRef]
  40. Liu, C.; Zipser, E.J.; Cecil, D.J.; Nesbitt, S.W.; Sherwood, S. A Cloud and Precipitation Feature Database from Nine Years of TRMM Observations. J. Appl. Meteorol. Climatol. 2008, 47, 2712–2728. [Google Scholar] [CrossRef]
  41. Seto, S.; Iguchi, T.; Oki, T. The Basic Performance of a Precipitation Retrieval Algorithm for the Global Precipitation Measurement Mission’s Single/Dual-Frequency Radar Measurements. IEEE Trans. Geosci. Remote Sens. 2013, 51, 5239–5251. [Google Scholar] [CrossRef] [Green Version]
  42. Hamada, A.; Takayabu, Y.N. Improvements in Detection of Light Precipitation with the Global Precipitation Measurement Dual-Frequency Precipitation Radar (GPM DPR). J. Atmos. Ocean. Technol. 2016, 33, 653–667. [Google Scholar] [CrossRef]
  43. Beyerle, G.; Grunwaldt, L.; Heise, S.; Köhler, W.; König, R.; Michalak, G.; Rothacher, M.; Schmidt, T.; Wickert, J.; Tapley, B.D.; et al. First Results from the GPS Atmosphere Sounding Experiment TOR Aboard the TerraSAR-X Satellite. Atmos. Chem. Phys. 2011, 11, 6687–6699. [Google Scholar] [CrossRef] [Green Version]
  44. Beyerle, G.; Schmidt, T.; Michalak, G.; Heise, S.; Wickert, J.; Reigber, C. GPS Radio Occultation with GRACE: Atmospheric Profiling Utilizing the Zero Difference Technique. Geophys. Res. Lett. 2005, 32, L13806. [Google Scholar] [CrossRef] [Green Version]
  45. Smith, E.K.; Weintraub, S. The Constants in the Equation for Atmospheric Refractive Index at Radio Frequencies. Proc. IRE 1953, 41, 1035–1037. [Google Scholar] [CrossRef] [Green Version]
  46. Foelsche, U.; Borsche, M.; Steiner, A.K.; Gobiet, A.; Pirscher, B.; Kirchengast, G.; Wickert, J.; Schmidt, T. Observing Upper Troposphere-Lower Stratosphere Climate with Radio Occultation Data from the CHAMP Satellite. Clim. Dyn. 2008, 31, 49–65. [Google Scholar] [CrossRef] [Green Version]
  47. Kursinski, E.R.; Hajj, G.A.; Schofield, J.T.; Linfield, R.P.; Hardy, K.R. Observing Earth’s Atmosphere with Radio Occultation Measurements Using the Global Positioning System. J. Geophys. Res. Atmos. 1997, 102, 23429–23465. [Google Scholar] [CrossRef]
  48. Zeng, Z.; Sokolovskiy, S.; Schreiner, W.S.; Hunt, D. Representation of Vertical Atmospheric Structures by Radio Occultation Observations in the Upper Troposphere and Lower Stratosphere: Comparison to High-Resolution Radiosonde Profiles. J. Atmos. Ocean. Technol. 2019, 36, 655–670. [Google Scholar] [CrossRef]
  49. Steiner, A.K.; Ladstädter, F.; Ao, C.O.; Gleisner, H.; Ho, S.P.; Hunt, D.; Schmidt, T.; Foelsche, U.; Kirchengast, G.; Kuo, Y.H.; et al. Consistency and Structural Uncertainty of Multi-Mission GPS Radio Occultation Records. Atmos. Meas. Tech. 2020, 13, 2547–2575. [Google Scholar] [CrossRef]
  50. Labitzke, K.; Kunze, M. Variability in the Stratosphere: The Sun and the QBO. In Climate and Weather of the Sun-Earth System (CAWSES); TERRAPUB: Tokyo, Japan, 2009; pp. 257–278. [Google Scholar]
  51. Dee, D.P.; Uppala, S.M.; Simmons, A.J.; Berrisford, P.; Poli, P.; Kobayashi, S.; Andrae, U.; Balmaseda, M.A.; Balsamo, G.; Bauer, P.; et al. The ERA-Interim Reanalysis: Configuration and Performance of the Data Assimilation System. Q. J. R. Meteorol. Soc. 2011, 137, 553–597. [Google Scholar] [CrossRef]
  52. Zhang, X.; Gao, P.; Xu, X. Variations of the Tropopause Over Different Latitude Bands Observed Using COSMIC Radio Occultation Bending Angles. IEEE Trans. Geosci. Remote Sens. 2014, 52, 2339–2349. [Google Scholar] [CrossRef]
  53. Thompson, D.W.J.; Bony, S.; Li, Y. Thermodynamic Constraint on the Depth of the Global Tropospheric Circulation. Proc. Natl. Acad. Sci. USA 2017, 114, 8181–8186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. McIntyre, M.E.; Norton, W.A. Potential Vorticity Inversion on a Hemisphere. J. Atmos. Sci. 2000, 57, 1214–1235. [Google Scholar] [CrossRef]
  55. Kunz, A.; Konopka, P.; Müller, R.; Pan, L.L. Dynamical Tropopause Based on Isentropic Potential Vorticity Gradients. J. Geophys. Res. Atmos. 2011, 116, D01110. [Google Scholar] [CrossRef] [Green Version]
  56. Chen, P. Isentropic Cross-Tropopause Mass Exchange in the Extratropics. J. Geophys. Res. 1995, 100, 16661–16673. [Google Scholar] [CrossRef]
  57. Appenzeller, C.; Davies, H.C.; Norton, W.A. Fragmentation of Stratospheric Intrusions. J. Geophys. Res. Atmos. 1996, 101, 1435–1456. [Google Scholar] [CrossRef]
  58. Škerlak, B.; Sprenger, M.; Pfahl, S.; Tyrlis, E.; Wernli, H. Tropopause Folds in ERA-Interim: Global Climatology and Relation to Extreme Weather Events. J. Geophys. Res. Atmos. 2015, 120, 4860–4877. [Google Scholar] [CrossRef]
  59. Funatsu, B.M.; Waugh, D.W. Connections between Potential Vorticity Intrusions and Convection in the Eastern Tropical Pacific. J. Atmos. Sci. 2008, 65, 987–1002. [Google Scholar] [CrossRef]
  60. Goering, M.A.; Gallus, W.A.; Olsen, M.A.; Stanford, J.L. Role of Stratospheric Air in a Severe Weather Event: Analysis of Potential Vorticity and Total Ozone. J. Geophys. Res. 2001, 106, 11813–11823. [Google Scholar] [CrossRef]
  61. Seidel, D.J.; Randel, W.J. Recent Widening of the Tropical Belt: Evidence from Tropopause Observations. J. Geophys. Res. 2007, 112, D20113. [Google Scholar] [CrossRef] [Green Version]
  62. Johnston, B.; Xie, F. Characterizing Extratropical Tropopause Bimodality and Its Relationship to the Occurrence of Double Tropopauses Using COSMIC GPS Radio Occultation Observations. Remote Sens. 2020, 12, 1109. [Google Scholar] [CrossRef] [Green Version]
  63. Pan, L.L.; Randel, W.J.; Gille, J.C.; Hall, W.D.; Nardi, B.; Massie, S.; Yudin, V.; Khosravi, R.; Konopka, P.; Tarasick, D. Tropospheric Intrusions Associated with the Secondary Tropopause. J. Geophys. Res. 2009, 114, D10302. [Google Scholar] [CrossRef] [Green Version]
  64. Liu, N.; Liu, C.; Hayden, L. Climatology and Detection of Overshooting Convection From 4 Years of GPM Precipitation Radar and Passive Microwave Observations. J. Geophys. Res. Atmos. 2020, 125, e2019JD032003. [Google Scholar] [CrossRef]
  65. Machado, L.A.T.; Rossow, W.B.; Guedes, R.L.; Walker, A.W. Life Cycle Variations of Mesoscale Convective Systems over the Americas. Mon. Weather. Rev. 1998, 126, 1630–1654. [Google Scholar] [CrossRef]
  66. Houze, R.A. Observed Structure of Mesoscale Convective Systems and Implications for Large-Scale Heating. Q. J. R. Meteorol. Soc. 1989, 115, 425–461. [Google Scholar] [CrossRef]
  67. Schumacher, C.; Houze, R.A.; Kraucunas, I. The Tropical Dynamical Response to Latent Heating Estimates Derived from the TRMM Precipitation Radar. J. Atmos. Sci. 2004, 61, 1341–1358. [Google Scholar] [CrossRef] [Green Version]
  68. Joos, H.; Wernli, H. Influence of Microphysical Processes on the Potential Vorticity Development in a Warm Conveyor Belt: A Case-Study with the Limited-Area Model COSMO. Q. J. R. Meteorol. Soc. 2012, 138, 407–418. [Google Scholar] [CrossRef]
  69. Schemm, S.; Wernli, H. The Linkage between the Warm and the Cold Conveyor Belts in an Idealized Extratropical Cyclone. J. Atmos. Sci. 2014, 71, 1443–1459. [Google Scholar] [CrossRef] [Green Version]
  70. Chae, J.H.; Wu, D.L.; Read, W.G.; Sherwood, S.C. The Role of Tropical Deep Convective Clouds on Temperature, Water Vapor, and Dehydration in the Tropical Tropopause Layer (TTL). Atmos. Chem. Phys. 2011, 11, 3811–3821. [Google Scholar] [CrossRef] [Green Version]
  71. Grise, K.M.; Thompson, D.W.J.; Birner, T. A Global Survey of Static Stability in the Stratosphere and Upper Troposphere. J. Clim. 2010, 23, 2275–2292. [Google Scholar] [CrossRef]
  72. Uccellini, L.W. The Possible Influence of Upstream Upper-Level Baroclinic Processes on the Development of the QE II Storm. Mon. Weather. Rev. 1986, 114, 1019–1027. [Google Scholar] [CrossRef]
  73. Reed, R.J. A Study of a Characteristic Type of Upper-Level Frontogenesis. J. Meteorol. 1955, 12, 226–237. [Google Scholar] [CrossRef]
  74. Bosart, L.F.; Lin, S.C. A Diagnostic Analysis of the Presidents’ Day Storm of February 1979. Mon. Weather. Rev. 1984, 112, 2148–2177. [Google Scholar] [CrossRef] [Green Version]
  75. Biondi, R.; Steiner, A.K.; Kirchengast, G.; Rieckh, T. Characterization of Thermal Structure and Conditions for Overshooting of Tropical and Extratropical Cyclones with GPS Radio Occultation. Atmos. Chem. Phys. 2015, 15, 5181–5193. [Google Scholar] [CrossRef] [Green Version]
  76. Biondi, R.; Randel, W.J.; Ho, S.-P.; Neubert, T.; Syndergaard, S. Thermal Structure of Intense Convective Clouds Derived from GPS Radio Occultations. Atmos. Chem. Phys. 2012, 12, 5309–5318. [Google Scholar] [CrossRef] [Green Version]
  77. Randel, W.J.; Wu, F. The Polar Summer Tropopause Inversion Layer. J. Atmos. Sci. 2010, 67, 2572–2581. [Google Scholar] [CrossRef]
Figure 1. Climatology of UTLS monthly zonal mean temperatures (solid contours, in K) and temperature standard deviation (color contours, in K) between 80° N and 80° S derived from COSMIC GPS-RO profiles from 2006 to 2017. Lapse-rate tropopause height is shown with a solid white contour. Temperatures greater than 250 K are removed. Additional temperature contours are added between 205–225 K to better display UTLS structure.
Figure 1. Climatology of UTLS monthly zonal mean temperatures (solid contours, in K) and temperature standard deviation (color contours, in K) between 80° N and 80° S derived from COSMIC GPS-RO profiles from 2006 to 2017. Lapse-rate tropopause height is shown with a solid white contour. Temperatures greater than 250 K are removed. Additional temperature contours are added between 205–225 K to better display UTLS structure.
Atmosphere 13 00196 g001
Figure 2. Number of GPS-RO profiles for non-DSI (solid blue contours) and DSI (solid red contours) backgrounds from 2006 to 2017 for (a) January, (b) April, (c) July, and (d) October. Shaded contours display LRT height differences (km) between non-DSI and DSI background profiles, provided the sampling is at least 10 profiles within each grid.
Figure 2. Number of GPS-RO profiles for non-DSI (solid blue contours) and DSI (solid red contours) backgrounds from 2006 to 2017 for (a) January, (b) April, (c) July, and (d) October. Shaded contours display LRT height differences (km) between non-DSI and DSI background profiles, provided the sampling is at least 10 profiles within each grid.
Atmosphere 13 00196 g002
Figure 3. (a,b) Two overshooting non-DSI PFs with collocated GPS temperature profiles (solid red) along with their associated non-DSI (solid black) and DSI (dotted green) background profiles. Blue asterisks and the associated dashed lines indicate the height and temperature of the background and PF LRT. Dotted blue line indicates the PF maximum echo-top height. Blue-shaded area shows the region ±5 km around the PF’s LRT height for anomaly calculation. Temperature anomaly profiles are shown on the bottom (c,d) with the “zero” height corresponding to the PF’s LRT height.
Figure 3. (a,b) Two overshooting non-DSI PFs with collocated GPS temperature profiles (solid red) along with their associated non-DSI (solid black) and DSI (dotted green) background profiles. Blue asterisks and the associated dashed lines indicate the height and temperature of the background and PF LRT. Dotted blue line indicates the PF maximum echo-top height. Blue-shaded area shows the region ±5 km around the PF’s LRT height for anomaly calculation. Temperature anomaly profiles are shown on the bottom (c,d) with the “zero” height corresponding to the PF’s LRT height.
Atmosphere 13 00196 g003
Figure 4. Number of PF-collocated GPS-RO temperature profiles comparing the potential vorticity values (PVU) at the 320 K isentrope with the lapse rate tropopause height (a, km) and temperature (b, K) anomalies that occurred near the PF. Solid red lines (±2-PVU) separate non-DSI and DSI PFs.
Figure 4. Number of PF-collocated GPS-RO temperature profiles comparing the potential vorticity values (PVU) at the 320 K isentrope with the lapse rate tropopause height (a, km) and temperature (b, K) anomalies that occurred near the PF. Solid red lines (±2-PVU) separate non-DSI and DSI PFs.
Atmosphere 13 00196 g004
Figure 5. Distribution of extratropical PFs (RD > 0.5) observed by GPM from 2014 to 2017 for (a) non-DSI PFs categorized by size, (b) non-DSI PFs categorized by RD, (c) DSI PFs categorized by size, and (d) DSI PFs categorized by RD.
Figure 5. Distribution of extratropical PFs (RD > 0.5) observed by GPM from 2014 to 2017 for (a) non-DSI PFs categorized by size, (b) non-DSI PFs categorized by RD, (c) DSI PFs categorized by size, and (d) DSI PFs categorized by RD.
Atmosphere 13 00196 g005
Figure 6. Overshooting convection observed by GPM over South Dakota on the morning of 28 July 2015. (a) Precipitation radar (PR) Ku-band near-surface reflectivity (dBZ). (b) PR Ku-band maximum echo-top heights (km). (c) ERA-Interim PV at 320 K potential temperature over the GPM PR swath. (d) PF-collocated GPS-RO temperature profile (red) located at 44.2° N, 101° W (denoted by the asterisk in (ac)) and the corresponding background profile (black), along with the PF lapse rate tropopause height (red) and maximum echo-top height (blue). (e) Cross-section of PR Ku-band reflectivity (dBZ) through the center of the PF (indicated by the solid black line in (a,b)). Red line is the GMI 37 GHz polarization-corrected temperature (PCT) (K); blue line is the GMI 89 GHz PCT (K); and black line is the rain rate (10 mm h−1).
Figure 6. Overshooting convection observed by GPM over South Dakota on the morning of 28 July 2015. (a) Precipitation radar (PR) Ku-band near-surface reflectivity (dBZ). (b) PR Ku-band maximum echo-top heights (km). (c) ERA-Interim PV at 320 K potential temperature over the GPM PR swath. (d) PF-collocated GPS-RO temperature profile (red) located at 44.2° N, 101° W (denoted by the asterisk in (ac)) and the corresponding background profile (black), along with the PF lapse rate tropopause height (red) and maximum echo-top height (blue). (e) Cross-section of PR Ku-band reflectivity (dBZ) through the center of the PF (indicated by the solid black line in (a,b)). Red line is the GMI 37 GHz polarization-corrected temperature (PCT) (K); blue line is the GMI 89 GHz PCT (K); and black line is the rain rate (10 mm h−1).
Atmosphere 13 00196 g006
Figure 7. GPS-RO temperature anomalies (K) near non-DSI PFs classified by relative depth (top) and area (bottom, in km2) for (a,d) all PFs, (b,e) land PFs, and (c,f) oceanic PFs. The number of PFs within each category is also displayed after each label.
Figure 7. GPS-RO temperature anomalies (K) near non-DSI PFs classified by relative depth (top) and area (bottom, in km2) for (a,d) all PFs, (b,e) land PFs, and (c,f) oceanic PFs. The number of PFs within each category is also displayed after each label.
Atmosphere 13 00196 g007
Figure 8. Monthly mean GPS-RO temperature anomaly profiles (K) near non-DSI PFs over (a) 20–40° N, (b) 40–65° N, (c) 20–40° S, and (d) 40–65° S. The number of PFs within each category is also displayed after each label.
Figure 8. Monthly mean GPS-RO temperature anomaly profiles (K) near non-DSI PFs over (a) 20–40° N, (b) 40–65° N, (c) 20–40° S, and (d) 40–65° S. The number of PFs within each category is also displayed after each label.
Atmosphere 13 00196 g008
Figure 9. GPS-RO temperature anomalies (K) near DSI PFs classified by relative depth (top) and area (bottom, in km2) for (a,d) all PFs, (b,e) land PFs, and (c,f) oceanic PFs. The number of PFs within each category is also displayed after each label.
Figure 9. GPS-RO temperature anomalies (K) near DSI PFs classified by relative depth (top) and area (bottom, in km2) for (a,d) all PFs, (b,e) land PFs, and (c,f) oceanic PFs. The number of PFs within each category is also displayed after each label.
Atmosphere 13 00196 g009
Figure 10. Monthly mean GPS-RO temperature anomaly profiles (K) near DSI PFs over (a) 40−65° N and (b) 40−65° S. The number of PFs within each category is also displayed after each label.
Figure 10. Monthly mean GPS-RO temperature anomaly profiles (K) near DSI PFs over (a) 40−65° N and (b) 40−65° S. The number of PFs within each category is also displayed after each label.
Atmosphere 13 00196 g010
Figure 11. Two-dimensional histogram (solid contours) and joint distribution of lapse rate tropopause height and temperature anomalies (shaded contours) with corresponding PF size and relative depth for non-DSI (a,c) and DSI (b,d) PFs. Note the color bar scale differences in each panel.
Figure 11. Two-dimensional histogram (solid contours) and joint distribution of lapse rate tropopause height and temperature anomalies (shaded contours) with corresponding PF size and relative depth for non-DSI (a,c) and DSI (b,d) PFs. Note the color bar scale differences in each panel.
Atmosphere 13 00196 g011
Figure 12. Two-dimensional histogram of the PF LRT temperature (K) and height (km) anomalies (solid contours) and the joint distribution of the temperature anomalies (K) at the PF LRT from Figure 7 and Figure 9 (shaded contours) for (a) non-DSI and (b) DSI PFs. Note the different axis ranges in each panel. Also shown are examples of common PF-collocated GPS temperature profiles (red) for (c) non-DSI (5 July 2015) and (d) DSI (2 May 2016) PFs along with their associated background profiles (black). Dashed lines indicate the tropopause height and temperature for both profiles and dotted lines indicates maximum echo-top height.
Figure 12. Two-dimensional histogram of the PF LRT temperature (K) and height (km) anomalies (solid contours) and the joint distribution of the temperature anomalies (K) at the PF LRT from Figure 7 and Figure 9 (shaded contours) for (a) non-DSI and (b) DSI PFs. Note the different axis ranges in each panel. Also shown are examples of common PF-collocated GPS temperature profiles (red) for (c) non-DSI (5 July 2015) and (d) DSI (2 May 2016) PFs along with their associated background profiles (black). Dashed lines indicate the tropopause height and temperature for both profiles and dotted lines indicates maximum echo-top height.
Atmosphere 13 00196 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Johnston, B.R.; Xie, F.; Liu, C. Relationships between Extratropical Precipitation Systems and UTLS Temperatures and Tropopause Height from GPM and GPS-RO. Atmosphere 2022, 13, 196. https://doi.org/10.3390/atmos13020196

AMA Style

Johnston BR, Xie F, Liu C. Relationships between Extratropical Precipitation Systems and UTLS Temperatures and Tropopause Height from GPM and GPS-RO. Atmosphere. 2022; 13(2):196. https://doi.org/10.3390/atmos13020196

Chicago/Turabian Style

Johnston, Benjamin R., Feiqin Xie, and Chuntao Liu. 2022. "Relationships between Extratropical Precipitation Systems and UTLS Temperatures and Tropopause Height from GPM and GPS-RO" Atmosphere 13, no. 2: 196. https://doi.org/10.3390/atmos13020196

APA Style

Johnston, B. R., Xie, F., & Liu, C. (2022). Relationships between Extratropical Precipitation Systems and UTLS Temperatures and Tropopause Height from GPM and GPS-RO. Atmosphere, 13(2), 196. https://doi.org/10.3390/atmos13020196

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop