The diurnal cycle of the clouds extending above the tropical tropopause observed by spaceborne lidar

. The presence of clouds above the tropopause over tropical convection centers has so far been documented by spaceborne instruments that are either sun-synchronous or insensitive to thin cloud layers. Here we document, for the ﬁrst time through direct observation by spaceborne lidar, how the tropical cloud fraction evolves above the tropopause throughout the day. After conﬁrming previous studies that found such clouds most frequently above convection centers, we show that stratospheric clouds and their vertical extent above the tropopause follow a diurnal rhythm linked to convective activity. The diurnal cycle of the stratospheric clouds displays two maxima: one in the early night (19:00– 20:00 LT) and a later one (00:00–01:00 LT). Stratospheric clouds extend up to 0.5–1 km above the tropopause during nighttime, when they are the most frequent. The frequency and the vertical extent of stratospheric clouds is very limited during daytime, and when present they are found very close to the tropopause. Results are similar over the major convection centers (Africa, South America and the Warm Pool), with more clouds above land in DJF (December–January– February)


Scientific context and objectives
Low-stratospheric clouds impact the atmospheric system in several ways. First, their larger heating rate than in the clear sky (Corti et al., 2006) increases the upward mass flux and fosters the large-scale upward transport of water above the tropopause. At the hour timescale, the cloud particles penetrating the stratosphere via overshooting convection leads, on the one hand, to a direct stratospheric humidification (Schoeberl et al., 2019;Dauhut et al., 2018). On the other hand, these particles can serve as support for ice scavenging: under saturated conditions, the water vapor deposits on the particles, which grow and fall out (Corti et al., 2008), decreasing low-stratosphere humidity (Jensen et al., 2013). By all these effects the stratospheric clouds modulate the stratospheric water vapor concentrations (Iwasaki et al., 2015) and affect the overall dynamical structure near the tropopause (Corti et al., 2006), at timescales down to 1 h. This is why it is important to understand the formation and the sub-daily evolution of such clouds.
The presence of ice clouds near the tropical tropopause has long been documented by in situ measurements (e.g., Thomas et al., 2002;Jensen et al., 2013;Frey et al., 2014). Detecting occurrences of clouds extending above the tropopause by remote sensing requires documenting the vertical cloud profile with a fine resolution and a high sensitivity to optically thin clouds, which few instruments can reach. Lidar measurements are able to document such occurrences (e.g., Nee et al., 1998;Dupont et al., 2010;Gouveia et al., 2017), but for a long time they were limited to local case studies. Dessler (2009) was the first to use the cloud detections by the CALIPSO lidar (Cloud-Aerosol Lidar Infrared Pathfinder Satellite Observations) to investigate how clouds extend above the tropopause on a global scale. Pan and Munchak (2011) refined the results by using an advanced tropopause dataset. Both studies found that clouds extending into the stratosphere are frequent above seasonal deep-convection centers and rarely elsewhere, especially in midlatitudes. Both studies deplored that the fixed overpass local time of the CALIPSO dataset is far from the late af-ternoon, when land convection is at its maximum. More recently, Wang et al. (2019) documented the presence of laminar cirrus in 10 years of CALIPSO data and reported a nonnegligible cloud amount above the tropopause. Because of the sun-synchronous orbit of CALIPSO, none of these studies were able to document the diurnal cycle of the stratospheric clouds.
The diurnal evolution of the high-altitude cirrus clouds has been documented over some specific sites using groundbased lidars (Sassen et al., 2003;Dupont et al., 2010;Gouveia et al., 2017). Gouveia et al. (2017) documented the evolution of the integrated cloud fraction (no vertical distribution) over Amazonia; Sassen et al. (2003) documented the diurnal evolution of the composition of cirrus clouds over Salt Lake City; and Dupont et al. (2010) did the same over four observatories in France and in the United States. However, using ground-based lidar to document optically thin clouds extending above the tropopause is difficult for two reasons: (1) as the studies based on CALIPSO observations show, these clouds occur primarily in regions where operational ground-based sites are absent or very few (Pacific Ocean, Equatorial Africa and South America), and (2) these clouds are mainly associated with deep convection, which implies the presence of optically thick cloud systems in the troposphere beneath that will make the successful probing of optically thin clouds near the tropopause impossible in most cases due to the attenuation of lidar signals. This explains why the ground-based lidars do not document the diurnal cycle of the stratospheric clouds with a satisfying spatial and temporal coverage.
Describing the diurnal evolution of the high-altitude clouds from a global perspective becomes possible with the CATS (Cloud-Aerosol Transport System) lidar operated from the International Space Station (ISS) between February 2015and November 2017(McGill et al., 2015. Thanks to the ISS non-synchronous orbit, CATS was able to probe the vertical cloud distribution of a particular region at different times of the day (not only at 01:30 and 13:30 local time like the instruments on CALIPSO). Aggregating CATS detections over a region of interest and over enough time provides a statistical overview of the diurnal evolution of cloud vertical profiles over that region (Noel et al., 2018). Our work aims at using CATS observations to describe and better understand the diurnal evolution of the cloud fraction in the tropical stratosphere.
Finding the processes responsible for the formation of tropical stratospheric clouds proves difficult, just like with high-tropospheric clouds (Reverdy et al., 2012). Two processes have been mainly proposed. Overshooting convection can lead to the injection of ice crystals into the stratosphere (Dauhut et al., 2018;Lee et al., 2019). Stratospheric cooling triggered by gravity waves (Pfister et al., 2010) could also lead to so-called cloud in situ formation (Pan and Munchak, 2011). The ratio of stratospheric clouds that are formed in situ has not been estimated yet. The current study does not provide a further estimate, but by describing the spatiotemporal evolution of the stratospheric clouds, it highlights how important the convective activity is to drive the stratospheric cloudiness and how the twice-daily sampling by lidars on board sun-synchronous platforms can miss the highest and largest stratospheric cloud fraction over certain regions.
In this paper, we document for the first time the diurnal cycle of clouds above the tropopause in the Tropics and the extent of their penetration in the stratosphere thanks to the high vertical and temporal resolution of the cloud detection by the CATS lidar. After describing CATS cloud data and the method to retrieve the tropopause heights used to detect clouds extending in the stratosphere (Sect. 2), we present maps of stratospheric clouds and document their diurnal cycle in regions of interest (Sect. 3). We then summarize our results and conclude (Sect. 4).

CATS cloud data
Between February 2015 and November 2017, the CATS lidar reported profiles at a vertical resolution of 60 m every 350 m along-track, with an average repeat cycle of nearly 3 d (Yorks et al., 2016). CATS Level 2 Operational layer files (L2O files; Palm et al., 2016) describe altitudes where cloud layers were detected within profiles of backscatter coefficients measured at 1064 nm by the CATS lidar (Pauly et al., 2019), averaged 5 km along-track. We considered all such files over the CATS operation period and inspected each 5 km profile within. For profiles located in the Tropics (30 • S-30 • N), we inspected each atmospheric layer therein identified as a cloud layer according to the CATS layer type information. As in Noel et al. (2018), we considered layers with a feature type score above 6, to avoid any possibly mislabeled aerosol layers. We flagged the cloud layers with a top altitude above the tropopause. Since any CATS L2O layer entirely above the tropopause is labeled as an aerosol layer (like in CALIPSO; Pan and Munchak, 2011), our study will not include clouds with their base in the stratosphere. Davis et al. (2010) noted that lidars in space may miss the thinnest subvisible cirrus clouds, but with enough spatial averaging cirrus clouds with optical depths near 0.001 can be detected (Martins et al., 2011). Lidar cloud detections also suffer from a lower sensitivity in the presence of sunlight, which induces significant additional noise in the lidar signal, but climatologies are still relevant (Noel et al., 2018).

Tropopause heights
To obtain the tropopause height, we considered profiles of temperature and pressure from the ERA5 (European Centre for Medium-Range Weather Forecasts -ECMWF -reanalysis) reanalysis dataset (Albergel et al., 2018). These profiles are available every 6 h, on 37 vertical levels and on a Atmos. Chem. Phys., 20, 3921-3929, 2020 www.atmos-chem-phys.net/20/3921/2020/ 0.25 • × 0.25 • horizontal grid. Such profiles in the ERA5 reanalysis agree well with observations in the high tropical troposphere (Podglajen et al., 2014). Using these profiles, we computed the vertical lapse rate profile (as in Reichler et al., 2003) and interpolated it on a 100 m vertical grid. We then applied the World Meteorological Organization (WMO) criteria defining the presence of a tropopause -i.e., the lowest altitude at which the lapse rate falls below 2 • C km −1 , provided the lapse rate between this altitude and all higher altitudes within 2 km does not exceed 2 • C km −1 (WMO, 1957). Following the WMO definition, we also allowed for the possibility of a second tropopause if the lapse rate exceeds 3 • C km −1 at least 1 km above the first tropopause. In such a case, we started to look for another tropopause above. To limit computation overhead, we constrained the search below 22 km. Using the WMO tropopause definition further allows us to compare our results to previous efforts based on the CALIPSO database that used the same definition (Pan and Munchak, 2011).

Stratospheric cloud detection
For a given CATS 5 km profile (Sect. 2.1), we identified the ERA5 tropopause height (Sect. 2.2) closest in time and location. Given the 6 h time resolution of the ERA5 reanalysis, there is at most a 3 h difference between the observation time and the thermodynamic information used to retrieve the tropopause height. We used the cloud information contained in the 5 km profiles in two ways. First, in 2 • × 5 • latitude-longitude bins we counted how many profiles contained a cloud extending above the tropopause, compared to the total number of profiles in the bins. Aggregating such numbers observed in JJA (June-July-August) and DJF (December-January-February) over the CATS operation period produced seasonal maps of above-tropopause cloud amounts (Sect. 3.1). Second, from each CATS 5 km profile we built a vertical cloud mask, using the tropopause height as the vertical reference and considering clouds that extend above it. Within regions chosen based on the seasonal maps, we aggregated such cloud masks over the same periods as above, keeping also track of the local time of observation for the considered mask. This produced regional vertical cloud fraction profiles above the tropopause, with one profile for each local time of observation (Sect. 3.2). Figure 1 shows the fraction of CATS profiles in which a cloud is detected above the tropopause, in all DJF (top) and JJA (bottom) months of CATS operation. Figure 1 shows that clouds in the tropical stratosphere are mostly detected over continents (South America, Equatorial Africa and land masses in the South Warm Pool in DJF; Cen-tral America, Central Africa and land masses in the North Warm Pool in JJA). The cloud fraction in the lower stratosphere is largest in DJF, up to 24 % over central Amazonia and coastal areas in the South Warm Pool and up to 20 % over Equatorial Africa. It is significantly lower in JJA, up to 12 % over Africa and 16 % over the North Warm Pool, even though the lowermost stratosphere (380-420 K potential temperature) is moister in JJA than in DJF (cf. e.g., Fig. 8c in Fueglistaler et al., 2009). The seasonal variation of the stratospheric cloud fraction is neither due to changes in the tropopause height, as the tropopause is lower in JJA than in DJF, over all regions (0.2 to 0.9 km altitude difference; see Appendix A), in line with the zonally averaged tropopause heights presented by Fueglistaler et al. (2009) andRieckh et al. (2014). Several factors may contribute to this seasonal variation: the density and strength of the convective systems (Liu and Zipser, 2005), their propensity to propagate or to be stationary (Houze et al., 2015), and the activity and efficiency of the in situ formation processes (Jensen et al., 2001;Jensen and Pfister, 2004).

Stratospheric cloud distributions
The spatiotemporal distribution of the stratospheric clouds is in very good agreement with the 4-year climatology of Pan and Munchak (2011) from CALIPSO observations. The DJF distribution also matches the CALIPSO cirrus detection at 100 hPa reported by Wang et al. (2019) for January 2009 very well. We report lower cloud frequencies than , which can be explained by the fact that we investigated slightly higher altitudes. Both CATS and CALIPSO datasets find (1) significantly weaker stratospheric cloud fraction in JJA than in DJF and (2) near-zero stratospheric clouds in the subtropics. These results are also consistent with the CALIPSO cloud fractions near 16 km reported by Schoeberl et al. (2019). Since those studies consider cloud detections derived from a spaceborne lidar instrument, over several years for most, their good agreement suggests that the CATS stratospheric cloud detections at 1064 nm are as reliable as the CALIPSO ones at 532 nm. A first conclusion of our results is therefore that CATS measurements strongly support the findings of all other studies using detections of high clouds from CALIPSO data.
Our CATS results are also in very good agreement with the distributions of clouds near the tropopause from other space instruments: 2006-2007 HIRDLS (High-resolution dynamics limb sounder) observations reported by Massie et al. (2010), 2006-2014 CloudSat observations (Kim et al., 2018) and the pioneering 1989 passive Stratospheric Aerosol and Gas Experiment (SAGE) II observations (Jensen et al., 1996). Besides the specificity in the cloud detection method employed by each instrument (occultation for HIRDLS and SAGE II and radar backscattering for CloudSat), the small differences between the distributions mostly come from the year-to-year variability. Larger differences can be found with the distributions of clouds penetrating the tropical tropopause tation Radar (Liu and Zipser, 2005). The densities of overshooting systems with tops in the lower stratosphere (on which Liu and Zipser, 2005, focused rather than all stratospheric clouds) are remarkably larger in Central America and Central Africa than over the Warm Pool. Since TRMM precipitation radar reflectivities are less sensitive to thin ice particles than CATS and CALIPSO lidars, we can interpret this difference by the fact that the American and African systems, though frequently overshooting the stratosphere, produce less thin stratospheric clouds than the Asian systems or other in situ processes (like gravity wave cooling) are more efficient to produce stratospheric clouds over Asia than America and Africa.

Diurnal cycle of cloud fractions in the tropical stratosphere
In contrast with the previous studies, the CATS dataset allows us to analyze the diurnal cycle of the cloud fraction in stratosphere. The cloud fraction at a regional scale shows a consistent diurnal cycle, robust over the different regions identified in the previous section (Fig. 2). In particular and in contrast to the diurnal cycle of surface precipitation, there is no land-ocean difference. All exhibit a pronounced minimum of about 2 %-4 % during the daytime, from 07:00 to 16:00 LT. They all present a first maximum at 19:00 or 20:00 LT (early-night peak), up to 16.5 % over Equatorial Africa. For all regions except South America and the North Warm Pool, this maximum is the largest cloud fraction of the day. All regions also present a second peak (late-night peak) at 00:00 or 01:00 LT (23:00 LT for the West Pacific and 02:00 LT for Central Africa), up to 16.5 % over South America. The midnight peak over Equatorial Africa is less clear than over the other regions because of the large variations between 23:00 and 03:00 LT. The capability of a longer dataset to produce a clearer signal is to be investigated. The cirrus clouds observed over Amazonia by groundbased lidar (Gouveia et al., 2017) shows a very similar diurnal cycle: a first peak in the early night (at 18:00-19:00 LT) and a second peak later in the night (at 02:00-03:00 LT). Although Gouveia et al. (2017) do not consider the cloud above the tropopause only, their distinction between subvisible, thin and opaque cirrus clouds indicates that the opaque cirrus clouds are predominant during the early night (18:00-21:00 LT) and that the thin cirrus (and subvisible ones during the dry season) dominate during the later night (from 00:00 and 02:00 LT onward, in wet and dry seasons, respectively).
The very deep convection transports cloudy air masses beyond the tropopause via overshoots and then directly con-Atmos. Chem. Phys., 20, 3921-3929, 2020 www.atmos-chem-phys.net/20/3921/2020/ tributes to the stratospheric cloud fraction (Dauhut et al., 2016(Dauhut et al., , 2018. The diurnal cycle of the stratospheric cloud fraction observed by CATS can at the first order be explained by the diurnal cycle of very deep convection over land (Liu and Zipser, 2005), especially (i) the minimal value during daytime and (ii) the first peak in the early evening. This first peak occurs with a delay of 3 to 4 h compared to the verydeep-convection maximum. As the dataset used by Liu and Zipser (2005) is more sensitive to overshoots freshly developed into the stratosphere, this delay can be explained by the subsequent horizontal expansion of the overshoots and their spread by the winds (Dauhut et al., 2018;Lee et al., 2019). The convective generation of gravity waves, which produce transient cooling off the convective centers and in some conditions trigger cloud formation, can also contribute to the increase of the stratospheric cloud fraction after the maximum of the very deep convection and then explain the delay of the first peak and potentially the second peak. It may also explain the similar diurnal cycle over the ocean regions, either close (South Warm Pool) or remote (West Pacific) from land masses. This process remains to be investigated. Figure 3 shows how far above the tropopause the clouds extend, depending on the local time in each tropical region (Sect. 3.1). Some regions are considered in DJF, while others are considered in JJA, because the stratospheric cloud distribution changes throughout the year (Fig. 1), following the ITCZ (Intertropical Convergence Zone) position. Pat-terns appear very consistent in all the regions considered. In all regions the largest cloud fractions are found near the tropopause, with few clouds extending higher. Cloud fractions extend relatively high (up to 1 km above the tropopause) during the early night. The first peak of cloud fraction, near 19:00-20:00 LT (Fig. 2), is associated with the all-day maximum of cloud vertical extent, with clouds in 5 % of profiles reaching 1 km above the tropopause in DJF regions. During the rest of the night (after 00:00 LT) clouds are still present but extend less high (except over South America). During daytime (06:00-18:00) clouds appear very close to the tropopause. Cloud fractions are overall much smaller in JJA (maximum of 5 %-10 %; Fig. 3d, f) than in DJF (max 10 %-12 %, Fig. 3,a, b, c, e), even if the tropopause is lower in JJA than in DJF (Appendix A).
In addition to describing the evolution of the stratospheric cloud cover at hourly timescales, these observations help interpret observations with limited temporal sampling (Noel et al., 2018). The Microwave Limb Sounder (MLS), like CALIPSO and all other instruments on board platforms of the A-train (Afternoon Train), samples the atmosphere at 01:30 and 13:30 LT, providing one single night and one single day observation. Some authors (e.g., Dion et al., 2019) attempt to retrieve the diurnal cycle of the observed water contents in the tropopause region, combining MLS observations with higher-temporal-resolution observations of convective activity based on TRMM observations of precipita- tion. Dion et al. (2019) assumed an in-phase relationship between precipitation and ice water content in the upper troposphere and at the tropopause level. For the stratospheric ice water content, MLS data still provide a signal-to-noise ratio that is too low. For future investigations, our results indicate that the stratospheric cloud fraction at 13:30 LT is, whatever the region, close to the minimal value of its diurnal cycle, whereas at 01:30 LT it is more typical of the second maximum. Carminati et al. (2014) investigated, from MLS measurements between 2005 and 2012, the differences between day and night ice water contents in the upper troposphere and the tropopause level. Unlike the stratospheric cloud fraction, tropopause ice water contents are larger at 13:30 LT than at 01:30 LT over Equatorial Africa during DJF, over Central Africa during JJA and over South America during both seasons. A possible explanation to reconcile our results is that tropopause ice water content is more sensitive to fresh convective activity (very-deep-convection occurrence), whereas the stratospheric cloud cover is more sensitive to the diffusion of the injected ice in the stratosphere.

Conclusions
Our results show how clouds in the tropical stratosphere are strongly concentrated above deep-convection centers, are almost absent in subtropical regions, and are more frequent in DJF than JJA and over land than over the ocean. In addition to these results, which are consistent with most previ-ous studies, we also show that both the cloud fraction and its extension above the tropopause follow a diurnal rhythm with a maximum during the early nighttime and a near-zero minimum during daytime. During daytime, the stratospheric clouds are limited to the first 100 m above the tropopause. During nighttime, a significant average cloud fraction is found up to 1 km above the tropopause. A second maximum of stratospheric cloud fraction is observed over all regions, generally a little after midnight. These results highlight how much the evolution of stratospheric clouds can be undersampled by other spatial instruments restricted to 01:30 and 13:30 LT, which then miss for instance the first maximum and the deepest development of stratospheric clouds in the early night. The very-deep-convective activity over tropical lands drives most of this diurnal cycle and leads in particular to the minimal stratospheric cloud fraction during daytime and the second peak during nighttime, both consistent over all regions. Further investigation is necessary to describe how convection contributes to this diurnal cycle and to assess the role of other processes leading to stratospheric cloud formation like the gravity waves. Finally further research is needed to understand why the timing of this diurnal cycle is very similar over land and over the ocean.
Data availability. Primary data and scripts used in the analysis and other supplementary information that may be useful for reproducing the author's work are archived by the Max Planck Institute for Meteorology and can be obtained by contacting publications@mpimet.mpg.de.
Author contributions. TD and VN designed the data analyses, and VN carried them out. TD prepared the paper with contributions from VN and IAD.
Competing interests. The authors declare that they have no conflict of interest.
Acknowledgements. The authors would like to thank Bernard Legras (LMD) and Florian Pantillon (LA) for useful discussions on the quality of tropopause altitudes from various sources. They thank the two anonymous reviewers for their comments and suggestions, which greatly helped to improve the paper. They also thank NASA Earthdata for access to CATS data and ECMWF for access to the ERA5 reanalysis dataset. They processed CATS and ERA5 data using the storage and computing facilities of Climserv (part of the IPSL ESPRI Mesocenter) and ICARE from AERIS.
Financial support. This research has been supported by the IDEX University of Toulouse (TEASAO project).
The article processing charges for this open-access publication were covered by the Max Planck Society.
Review statement. This paper was edited by Michael Pitts and reviewed by two anonymous referees.