A Survey of Precipitation-Induced Atmospheric Cold Pools over Oceans and Their Interactions with the Larger-Scale Environment

Pools of air cooled by partial rain evaporation span up to several hundreds of kilometers in nature and typically last less than 1 day, ultimately losing their identity to the large-scale flow. These fundamentally differ in character from the radiatively-driven dry pools defining convective aggregation. Advancement in remote sensing and in computer capabilities has promoted exploration of how precipitation-induced cold pool processes modify the convective spectrum and life cycle. This contribution surveys current understanding of such cold pools over the tropical and subtropical oceans. In shallow convection with low rain rates, the cold pools moisten, preserving the near-surface equivalent potential temperature or increasing it if the surface moisture fluxes cannot ventilate beyond the new surface layer; both conditions indicate downdraft origin air from within the boundary layer. When rain rates exceed ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document} 2 mm h-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}, convective-scale downdrafts can bring down drier air of lower equivalent potential temperature from above the boundary layer. The resulting density currents facilitate the lifting of locally thermodynamically favorable air and can impose an arc-shaped mesoscale cloud organization. This organization allows clouds capable of reaching 4–5 km within otherwise dry environments. These are more commonly observed in the northern hemisphere trade wind regime, where the flow to the intertropical convergence zone is unimpeded by the equator. Their near-surface air properties share much with those shown from cold pools sampled in the equatorial Indian Ocean. Cold pools are most effective at influencing the mesoscale organization when the atmosphere is moist in the lower free troposphere and dry above, suggesting an optimal range of water vapor paths. Outstanding questions on the relationship between cold pools, their accompanying moisture distribution and cloud cover are detailed further. Near-surface water vapor rings are documented in one model inside but near the cold pool edge; these are not consistent with observations, but do improve with smaller horizontal grid spacings.


Introduction
Cold pools are defined by the American Meteorological Society (AMS) as ''a region or pool of relatively cold air surrounded by warmer air,'' or, ''any large-scale mass of cold air'' (Glickman 2000). In one example, nighttime radiative cooling of a land surface can create a near-surface cold pool of air. Over the open ocean, far from land, the cooling of near-surface air is arguably exclusively caused by the partial evaporation of precipitation within subsaturated air. This precipitation can be liquid or solid, though in subtropical and tropical regions the precipitation is entirely liquid near the surface. The precipitationinduced downdrafts introduce denser air underneath warmer, lighter environmental air, giving rise to a density current. The increased surface pressure establishes a horizontal pressure gradient force that drives the cold pool air outward, establishing an outflow boundary known as a gust front. The air inside the cold pool is less able to support surfacebased buoyancy-driven convection. The gust front itself behaves as a material surface, and air moving over the gust front can encourage further secondary convection through mechanical lifting and anomalous buoyancy. Such convectively produced atmospheric cold pools are typically observed to span 10-200 km in diameter, and to last for less than a day, after which they lose their individual identity to the large-scale or synoptic flow. This range of sizes and lifetimes identifies them as mesoscale features, although cold pools are often embedded in and can help define the edges of larger synoptic systems (e.g., squall lines or hurricane rainbands).
Although the latent cooling from evaporation of precipitation is necessary for creating a cold pool, not all precipitation leads to well-defined cold pools at the bottom of the atmosphere. Precipitative downdrafts are driven by both condensate loading and evaporation, and, if the resulting downdraft is not strong enough to reach the surface, elevated patches of moist, cool air remain. These most likely will not fit the AMS cold pool definition, in that they are unlikely to be completely surrounded by warmer air. One example, common in stratocumulus clouds because of their relatively smaller precipitation drop sizes (e.g., Wood 2005b), is light precipitation (drizzle) evaporating below cloud base (e.g., Wood 2005a). The local cooling can serve to destabilize the subcloud layer and reinforce the updrafts in stratocumulus clouds (e.g., Feingold et al. 1996) and can also be visualized in lower-level scud clouds. In the shallow cumulus regime, such an elevated enhanced moisture layer may aid future convection (e.g., Li et al. 2014;Schlemmer and Hohenegger 2014). As drops reach larger sizes through collision-coalescence (e.g., Baker et al. 2009) and their corresponding size-dependent fall velocities increase, drops can reach the ground with only partial evaporation.
The convective organization established by surface-based precipitation-induced cold pools differs fundamentally from the larger-scale convective self-aggregation discussed elsewhere in this volume (e.g., Holloway et al. 2017;Wing et al. 2017). In convective selfaggregation, radiative subsidence produces an expansive ([ 1000 km), convectively suppressed, dry region (a ''dry pool'') that is growing in time, neighboring a moist, deep convective region. Although an analogy between the radiative dry pool and evaporationdriven cold pools is often made, there are important differences. In the simulation shown in Fig. 1a and b, the driest precipitable water (PW) percentile corresponds to the ''dry pool,'' and the moistest PW percentile to the convective regime capable of supporting cold pools. A dry region is seen to form and expand, confining the deep convection to the remaining moist area. Convective self-aggregation is most clearly exhibited within radiative-convective equilibrium models of tropical convection. Observations in nature of convective self-aggregation remain elusive ) and difficult to attribute to specific feedbacks, but one study examining observed aggregation depicts properties similar to modeled self-aggregation, including large-scale drying (Tobin et al. 2012).
Convective organization leads to an enhancement of moisture gradients, as moist regions become moister and dry regions become drier, evident in Fig. 1. The mesoscale convective cold pools within the convectively self-aggregated regions can act to oppose the aggregation, by transporting water vapor from moist to dry areas, reflecting the divergence resulting from strong downdrafts in the subcloud layer in the model moist cold pools (Fig. 1e). In nature, the most dramatic visual example of this is arguably cold-pool-containing squall lines emanating from the African monsoon into the Sahara desert (e.g., Flamant et al. 2009;Trzeciak et al. 2017). Over the tropical ocean, cold pools within  ) also shows the low-level winds, which are seen to converge into the moist aggregate. Vertical profiles of c temperature, d water vapor and e vertical velocities in the radiatively subsiding dry pool, identified by PW 20th percentile, and in the precipitation-induced cold pools, identified by the near-surface (first atmospheric level z ¼ 37 m) temperature T air 20th percentile, and precipitation [ 0 at the end of the simulation. An inset in the first 2 km added to c to indicate the inversion in the dry pool at the top of the subcloud layer Surv Geophys (2017) 38:1283-1305 1285 simulated squall lines help broaden the precipitating intertropical convergence zone poleward (Nolan et al. 2016). Precipitatively-generated cold pools also counteract convective self-aggregation through suppression of local convection by divergence of the nearsurface air (Fig. 1e) (Jeevanjee and Romps 2013;Muller and Bony 2015). In theory, such mesoscale subsidence within a precipitation cold pool could morph into a new, larger radiatively induced ''dry pool,'' although this model behavior has not yet been witnessed, to our knowledge (see also Held et al. 1993;Muller and Held 2012). Cold pools act to disperse convection. Atmospheric cold pools over oceans are receiving attention for their ability to reorganize the mesoscale cloud distributions and potentially facilitate the transitions from high-to low-albedo shallow cloud cover. In the Tropics, the expansion of the spatial and height distribution of convection by cold pools may facilitate transitions from shallow to deep convective regimes, and therefore the eastward propagation of the Madden-Julian oscillation (MJO) into moistening environments (Rowe and Houze 2015;Feng et al. 2015;Ruppert and Johnson 2015;Schlemmer and Hohenegger 2016;Hannah et al. 2016;Ciesielski et al. 2017).
These foci on cold pool impacts justify the timeliness of a survey of cold pool characteristics and processes in different convective regimes. New observational capabilities and strategies were employed during the Rain in Cumulus over Ocean campaign (RICO; Rauber et al. 2007) and Dynamics of the MJO campaign (DYNAMO; Yoneyama et al. 2013) and will be during the upcoming Elucidating the role of clouds-circulation coupling in climate campaign (EUREC 4 A; Bony et al. 2017). These advance early observational studies stymied by imprecise data collocation and measurement misunderstanding (e.g., Warner et al. 1979;LeMone 1980) that had difficulty perceiving the larger-scale mesoscale organization. Measurements either came from one point, e.g., a ship (Addis et al. 1984;Young et al. 1995;Saxen and Rutledge 1998) or from aircraft (e.g., Zipser 1977;Kingsmill and Houze 1999), convoluting space and time. Modeling capabilities are improving also, either in their spatial grid spacing (Romps and Jeevanjee 2016;Grant and van den Heever 2016) or domain size Hohenegger 2014, 2016), and at times both (Seifert and Heus 2013;Vogel et al. 2016). The recent years have also seen advances in the formulation of cold pool parameterizations (e.g., Qian et al. 1998;Rozbicki et al. 1999;Rio et al. 2009Rio et al. , 2013Hohenegger and Bretherton 2011;Del Genio et al. 2015;Pantillon et al. 2015) and their coupling with convective schemes. These provide an avenue for rectifying a model diurnal cycle that is too closely linked to the solar cycle in models (Rio et al. 2009;Schlemmer and Hohenegger 2014).
The focus in this survey is over ocean in the subtropical and tropical latitudes. Cold pools are extremely important over land for initiating convection, where their depth and gust fronts are often substantial (e.g., Bryan and Parker 2010), surface sensible fluxes are significant and aerosol effects are more pronounced (e.g., van den Heever 2015, 2016;Schlemmer and Hohenegger 2016). Cold pools are also important for squall lines and tornadogenesis (e.g., Markowski and Richardson 2014). The omission of landbased cold pools is merely to keep the scope of this particular survey tractable. The focus on the oceanic regions equatorward of $ 30 N also emphasizes near-surface latent cooling induced by the partial evaporation of the liquid phase, as opposed to cooling by ice melting or sublimation. Cold pool behavior and significance vary with the depth of the originating convection, its degree of organization and relationship to the large-scale environment (trade wind vs. equatorial). This is reflected in the structure of the survey, which begins with shallow convection and moves toward deeper convection. The subtropical and trade wind regions are defined by climatologically steady winds imposing a background wind shear, extending to Barbados in the examples provided and cited literature. The tropical ocean examples and cited literature are primarily near-equatorial, where wind shear can be ignored at times, and convection is readily able to span the full free troposphere.

Cold Pools from Boundary Layers not Exceeding km Altitude
The strong inversion capping subtropical stratocumulus clouds maintains cloud top heights at $ 1.5 km or less (Zuidema et al. 2009), yet drizzle is ubiquitous (Leon et al. 2008). Some of the precipitating clouds occupying slightly deeper boundary layers (e.g., Mechem et al. 2012) are capable of downdrafting air that can reach the surface and develop cold pools (e.g., Fig. 2), also documented in Savic-Jovcic and  and Wood et al. (2011). Aircraft measurements within stratocumulus cold pools indicate that near the surface the equivalent potential temperature (h e ) tends to increase, rather than maintain a constant value or decrease (Zanten and Stevens 2005;Savic-Jovcic and Stevens 2008;Terai and Wood 2013). 1 The increase in h e is attributed to an accumulation of the surface fluxes underneath a cold pool capping stratification. Cold pool depth is difficult to observe. Inferences made from pressure increases suggest altitudes of $ 300 m (Terai and Wood 2013, and references therein). Further aloft, evaporation at constant moist static energy still brings air closer to saturation, increasing its susceptibility to convection. Wind speed convergence of the convectively favorable air near the cold pool edges helps perpetuate open-celled organization (Feingold et al. 2010), ventilating the accumulated surface fluxes and mixing air through the boundary layer. Thus, alterations to both the thermodynamic vertical structure and surface dynamics provide mechanisms by which cold pools can contribute to the longevity of precipitation within stratocumulus regions, in both closedand open-celled organizations. Precipitation is most pronounced pre-dawn (e.g., Burleyson et al. 2013), suggesting that cold pools are also most effective at night.
Radar observations indicate that the speed of advection of cold pools is similar to that of the cloud layer (Wilbanks et al. 2015). This suggests that cold pools are better thought of as tracers or artifacts of stratocumulus precipitation as opposed to drivers. This is also concluded from simulations relying on fixed cloud droplet numbers (Zhou et al. 2017) and is consistent with simulations showing little influence on cloud organization from the inhomogeneization of surface fluxes by cold pools (Kazil et al. 2014). Precipitation is necessary for transitions from closed-to open-cell stratocumulus cloud organizations (e.g., Savic-Jovcic and Stevens 2008; Xue et al. 2008;Wang and Feingold 2009;Feingold et al. 2010;Wood et al. 2011), thereby implicating cold pools in mesoscale organization transitions, but only indirectly. Aerosol number concentrations must also deplete sufficiently for open-celled or cumuli organization to form (Wood et al. 2011;Zhou et al. 2017;Yamaguchi et al. 2015). That said, aerosol budgets are unlikely to be influenced by cold pools per se, as the wind increases at cold pool fronts are shortlived (Terai and Wood 2013), except perhaps in strongly aerosol-depleted conditions (Kazil et al. 2014). It is also noteworthy that, while an important motivation for studying shallow mesoscale transitions is their influence on cloud fraction and the planetary albedo, a clear relationship between precipitation and cloud cover is not apparent in Fig. 2, as the precipitation feeds an upper stratiform layer below the trade wind inversion. Nuijens et al. (2015) clarify that it is this cloud layer, rather than a lower cloud layer at the lifting condensation level, that is most variable.
As air advects from the stratocumulus regions to warmer sea surface temperatures, the altitude of the trade wind capping inversion rises only slowly (Schubert 1995), but the weakening temperature inversion does begin to permit deeper convection. Early experiments studying the suppressed trade wind/tropical convective environment, namely the Atlantic Trade wind Experiment (ATEX; Augstein et al. 1973) and the Barbados Oceanographic and Meteorological Experiment (BOMEX; Nitta and Esbensen 1974), estimated the fractional area occupied by actively ascending cumuli to be a negligible 0.02 (Augstein et al. 1973). The prevailing cloud organization was isolated cumuli or cloud lines aligned with the mean wind that precipitated very little if at all (LeMone and Pennell 1976; Nair et al. 1998), suggesting that cold pools can likely be ignored for this environment (e.g., Albrecht 1993;Siebesma et al. 2003). This supported a paradigm begun with Riehl et al. (1951) in which precipitation within the boundary layer does not change the moist static energy, but rather the cooling introduced by evaporating rain is energetically balanced by an increase in moisture.

Cold Pools from Convection Reaching the Mid-Troposphere
The diffusing trade wind inversion strength does allow some deeper clouds to develop, however, and along with them, cold pools. Larger drop sizes, encouraged by a stronger collision-coalescence process, allow precipitation to return approximately one-third of the surface evaporation to the ocean in the trade wind cumulus region (Snodgrass et al. 2009). Visible in Figs  deeper clouds spawn downdrafts capable of bringing down air that is drier than the nearsurface air, contributing to a lowering of the equivalent potential temperature (Zuidema et al. 2012). The presence of less convectively favorable air near the surface, spreading out as a density current, explains an organization of mesoscale cloud arcs circumscribing mostly cloud-free regions. The small areal coverage of the deeper clouds can thus alter a much larger area of near-surface air. More quantitatively, the satellite-derived cloud fraction producing rain rates exceeding 1 mm h À1 is a mere 0.02 for wintertime Caribbean cumuli (Snodgrass et al. 2009), over an area of approximately $ 10 4 km 2 scanned by a   coincident with the Next-generation Airborne Remote Sensing for Validation (NARVAL) aircraft campaign . The wind is flowing from right to left. The two largest cold pools span approximately 100 km and are better defined on the downwind side, particularly the left side of the image. The higher cloud tops within the cold pool centers reach approximately 4 km w in d Fig. 4 Aircraft view of a cold pool taken on August 25, 2016, southeast of Barbados. Noteworthy are the cloud lines to the left of the image, with the nearby cold pool convection organized in a circle, a portion of which is oriented perpendicular to the cloud lines in the left-hand side. Isolated convection reaching a higher altitude is detraining, most likely into a layer of increased stability precipitation radar, also consistent with estimates from large-eddy-scale simulations (Neggers et al. 2002). The diurnal cycle is weak (Snodgrass et al. 2009).
A further notable feature of Figs. 3 and 4 is their siting at approximately 10 N of the intertropical convergence zone (ITCZ), after a multi-day equatorward advection of boundary layer air away from the stratocumulus regions. At what point along the journey shallow convection can deepen substantially remains underexplored. Space-based lidar suggests the northern hemisphere ocean basins are more conducive of 'deeper' shallow convection ( Fig. 1 of Medeiros et al. 2010), perhaps because northern hemisphere surface parcels do not advect over cooler equatorial waters in their journey to the ITCZ. The environmental conditions supporting shallow convection that can become deep enough to support coherent downdrafts bringing down air of lower h e from above the cloudy boundary layer is also not known and may depend on the history of the air parcel as well.
Liquid-only clouds reaching 4 km produce cold pools that are very similar to those produced by tropical deep convection extending throughout the depth of the troposphere, shown next.

Cold Pools from Deep Tropical Convection
The mean properties of $ 300 cold pools composited from conventional surface meteorological datasets from Gan Island (0.6 S, 73.1 E) and the Research Vessel Roger Revelle located $ 700 km to the east, at 80.5 E are shown in Fig. 5. A cold pool is identified through a temperature drop of 0.5 K, applied to a 5-min time series that has been smoothed using a Haar wavelet filter, combined with the requirement of no rain within the hour prior. The temperature drops are normalized to their mean time span of 20 min.  The cold pool front is defined through a temperature drop of 0.5 K applied to the 5-min time series smoothed using a Haar wavelet filter. The front is normalized to a 20-min time interval, corresponding to the mean frontal passage time. a, g Air temperature (black) and water vapor mixing ratio (red); b, h water vapor mixing ratio (red) and equivalent potential temperature (black); c, i surface wind speed (black) and surface pressure (red); d, j latent heat fluxes (LHF, black) and the Bowen ratio (surface heat fluxes/latent heat fluxes; red); e, k rain rate (black) and column water vapor (red, e) only); f liquid water path sample is drawn from the full tropical convective spectrum. The composite-mean includes a near-surface drying of $ 1 g kg À1 and decrease in the equivalent potential temperature of 2 K. The mean wind speed increases by $ 2 m s À1 , but only for 15-30 min. The pressure increase corresponds to a mean cold pool depth of approximately 300 m, is that stronger temperature drops are accompanied by stronger decreases in the water vapor mixing ratio (and h e ), indicating either higher origin heights for the downdrafts or more coherent structures less modified by environmental mixing (de Szoeke et al. 2017). The mean near-surface relative humidity of the cold pools indicate subsaturation (not shown), suggesting the downdrafts rarely if ever maintain saturation, despite being initially saturated. All of these traits are shared with cold pools in the northeast Atlantic trade wind regime (Zuidema et al. 2012).
Cold pools are much more likely to occur underneath obscuring cirrus shields (e.g., de Szoeke et al. 2017), but one example of a tropical cold pool visible from space is shown in Figs. 6 and 7, coinciding with the DYNAMO field campaign (Yoneyama et al. 2013). A vertically pointing cloud radar and meteorological station operated within the range of a scanning precipitation radar, with a range-height indicator scan dedicated to the precipitation radar azimuth encompassing the cloud radar (Feng et al. 2014), allowing a precise collocation of cold pool convective features with its surface features. The cold pool spanning $ 120 km is outlined on the 0615 UTC image, reaching the surface site at approximately 0915 UTC, where the vertically pointing radar reveals a cloud depth of 12 km in places lacking wind shear. The visible image shows other cold pools that are less obvious in the radar image. The surface meteorological time series shows a water vapor path increasing from 5.1 to 5.6 cm prior to the cold pool, whose edge has a temperature drop of $ 5 K, a more quickly recovering drop in the vapor mixing ratio, and a short, almost unidentifiable, wind speed increase. The relative humidity at altitudes above 5 km was 50% (not shown). Stratiform precipitation occurring later in the day prevented recovery of the near-surface temperature and maintained a near-surface relative humidity of 85-90%. Earlier observational studies focused on longer-lasting cold pools also selected from mesoscale systems with stratiform precipitation (e.g., Young et al. 1995), and longer-lasting surface fluxes changes are documented (Saxen and Rutledge 1998). For these, the altered surface fluxes, with their higher Bowen ratio, may be enough to influence the mean. This is one difference from cold pools in the trade wind regime. Deep tropical convection with stratiform precipitation can also sustain mesoscale downdrafts of warmer air. Early aircraft measurements concluded mesoscale downdrafts are too warm to reach the surface (Zipser 1977), but this is contradicted by Kilpatrick and Xie (2015), who relied on satellite scatterometer data combined with surface buoys.
Other influences on tropical cold pool characteristics, besides the level of mesoscale organization, include the atmospheric moisture distribution and the amount of mixing with environmental air during the downdraft. Since tropical deep convection easily attains cloud top heights of 8 km and upwards (Fig. 8), the downdraft air can in theory originate from a higher altitude. The modification of near-surface air properties is indeed more pronounced when the cloud top heights of the parent convection are higher (Fig. 8). Most studies point to an origin altitude for the downdraft air of 2 km or less, however (Betts 1976;Betts and Dias 1979;Torri and Kuang 2016a;de Szoeke et al. 2017;Schiro and Neelin 2017). That higher clouds are associated with stronger cold pools may reflect a correlation between downdraft width and clouds that are wider as a result of organization, discouraging environmental dilution of the downdraft air (Schlemmer and Hohenegger 2014; Schiro and . The low-altitude origin discounts contributions to the convective downdrafts (treated separately from stratiform precipitation) by the melting or sublimation of ice particles (Srivastava 1987).
Convection able to impinge upon relatively dry mid-tropospheres should in theory also produce stronger cold pools, by allowing more evaporation into the downdraft while simultaneously transporting drier air of lower h e to the surface (e.g., Chen et al. 2016). Many examples of drier atmospheres coexisting with more isolated convection are documented within shallow-to-deep transition studies (Feng et al. 2014;Ruppert and Johnson 2015), as well as with linearly organized convection such as squall lines (Takemi and Satomura 2000;Mapes et al. 2017;Schiro and Neelin 2017). A correspondence between stronger cold pools and drier mid-tropospheres has not yet been robustly identified in observations, but is consistent with published emphases (see also Takemi et al. 2004). One important compensation is the condensate loading within the downdrafts, as more condensate will decrease the surface parcel buoyancy (James and Markowski 2010; Torri and Kuang 2016a). Tropical cold pools are more frequent during times of lower outgoing longwave radiation, reflecting more organized deep convection with larger cloud covers, cooler sea surface temperatures (SSTs), and cooler near-surface layers that are closer to saturation (de Szoeke et al. 2017), and more fully saturated atmospheres. The number of observed cold pools only varies weakly with time of day (de Szoeke et al. 2017), though satellite measurements suggest mesoscale downdrafts able to reach the surface occur 8-12 h after the peak rainfall (Kilpatrick and Xie 2015). The boundary layer is cooled and moistened over a large area by the evaporation of stratiform precipitation, reducing the buoyancy of individual surface parcels. Nevertheless, cold pools and the surface-originating convection they reflect clearly do occur (de Szoeke et al. 2017), perhaps more than originally thought at the time of Houze and Betts (1981). This is also evident in space-based cloud radar observations of cumulus congestus/cumulonimbus occurring underneath upper-level stratiform cloud (e.g., Riley et al. 2011). The moisture distribution is influenced more by layer-lifting and less by surface-based buoyancy (see, e.g., in this issue, Mapes et al. (2017), dating back to at least Houze and Betts (1981), so that it cannot be argued that cold pools are important for redistributing moisture.
A connection between deep convective cold pools and those in the trade wind regions may then be their relationship to the column water vapor path. In the trade wind regions, column water vapor paths of 4.5-5.0 cm correspond to moister free tropospheres capable of supporting deeper convection and more cold pools (Zuidema et al. 2012). The DYNAMO water vapor path frequency distribution contains a plateau at 5.8-5.9 cm, and a mean of 5.1 cm , suggesting the slightly drier conditions in which cold pools can be more influential are close in value to those for the trade wind regions. This is consistent with the subcloud moisture field remaining critical for convective initiation (Kingsmill and Houze 1999;Seifert and Heus 2013), particularly at cloud base level (Takemi and Satomura 2000), as well as anomalous moisture in the lower free troposphere (Sherwood et al. 2010), since column water vapor paths between 4.5 and 5.5 cm allow drier mid-tropospheres to coexist with deep, moist lower free tropospheres ). One long-standing question in the Tropics, to which cold pools are relevant, is what maintains its boundary layer. The powerful Arakawa and Schubert (1974) parameterization assumes a quasi-equilibrium moisture closure in which surface moisture fluxes are balanced by the clear-sky entrainment of air drier than the saturation surface vapor mixing ratio. This is questioned within Raymond (1995), which concluded that surface fluxes must be balanced by convective downdrafts rather than clear-sky entrainment, because downdrafts provide a larger reduction in near surface h e . In the more comprehensive observational assessment of de Szoeke et al. (2017) based on DYNAMO soundings, downdrafts contribute to 20-30% of the boundary layer moisture and temperature budget, on days with diminished sea surface temperature (SST), and less otherwise. This suggests the effects of downdrafts are secondary, if not negligible, to those from clear-air entrainment through the boundary layer top, on the boundary layer h e . Thayer-Calder and Randall (2015) in an updated simulation also conclude clear-air entrainment predominantly balances surface evaporation. Torri and Kuang (2016a) also find a small contribution from convective downdrafts to the net flux of moist static energy into the boundary layer, and a much larger contribution from turbulent mixing across the boundary layer top, using a Lagrangian particle tracking analysis. These recent, independent studies all indicate a secondary role for downdrafts on the large-scale boundary layer moist static energy budget.

Remaining Questions
This survey motivates two remaining, intertwined questions.

The Relationship of Trade Wind Cold Pools to Cloud Cover
An outstanding question with shallow convection remains the still poorly known relationship between convection and cloud cover in the trade wind region, where the cloud cover is important for the planetary albedo. Modeling simulations, despite their many advances (e.g., Seifert and Heus 2013), remain inconclusive (van Zanten et al. 2011), to a substantial degree because many microphysical parameterizations are not optimal for trade wind cumuli representations (Li et al. 2015). While the cold pool itself discourages further surface-based convection, detrainment elsewhere in the troposphere increases cloud cover. In Fig. 4, the deeper convection is detraining at a higher altitude, affecting the overall cloud cover, arguably also visible in Fig. 3. Such detrainment can occur into layers of increased stability, and although the altitude of the higher stable layer is not clear in Fig. 2, the 0 C level, which occurs at $ 500 hPa, is enriched with such layers (Zuidema 1998;Stevens et al. 2017).
This undertaking also means critically assessing microphysical parameterizations and their interaction with representations of environmental mixing. Experimentation with two popular microphysical schemes has revealed significant differences in cloud fraction ensuing from cold pools that are similarly simulated (Li et al. 2015), and relevant observations should also be acquired to assess results from arguably better-suited schemes (Seifert and Heus 2013). The sensitivity of the convective downdraft to condensate loading and mixing above the boundary layer is also important. In addition, this means further evaluating the effects of low-and high-wind shear, as wind shear can interact dynamically with the convective structure (Li et al. 2014). Related to all these questions is also the timescale of the boundary layer recovery within a cold pool, and how often cold pools are able to fully recover before being impinged upon by a subsequent cold pool; observations suggest cold pools occur in clumps (e.g., Fig. 2). If so, and embedded in larger-scale moisture envelopes, those in tandem can alter the atmospheric longwave radiative cooling on longer time and larger spatial scales than of just an individual cold pool.
As computational capabilities continue to improve, the trade-off between spatial grid spacing and domain size in the modeling of cold pools can be expected to diminish. The enabled sophisticated examinations of cold pools with their larger-scale environment include examining the relative importance of colliding cold pools as opposed to single pools. The further development of tracking algorithms (Torri et al. 2015;Feng et al. 2015;Böing 2016; Drager and van den Heever 2017) within such more complex simulations will allow for a better understanding of collisions and their impact of the cold pool life cycle. Such efforts will also further parameterization development.

Thermodynamic Secondary Initiation Processes
Models, through their ability to capture the full four-dimensional fields (Fig. 9), highlight an important distinction between dynamical and thermodynamical influences on secondary convection. Early ideas for triggering new convection focused on dynamics, in which the cold density currents interacted with the prevailing near-surface wind shear to encourage upright updrafts favorable for convection (Rotunno et al. 1988). The trade wind regions in particular experience steady winds, with a maximum wind speed at cloud base decreasing both above and below. Dynamical forcing of boundary layer air is stronger when the cold pool gust front is aligned with the mean wind at the downwind side of a trade wind region cold pool, and the vorticity interactions help explain the more pronounced convection on the downwind side of a cold pool (e.g., in Figs. 3 and 4, further explored in Li et al. 2014). Stronger dynamical forcing is also related to higher downdraft heights in Jeevanjee and Romps (2015).
In tropical oceanic regions lacking rotation and strong wind shear (Fig. 5), a separate paradigm emphasizing the role of water vapor in initiating new convection has also been influential (Tompkins 2001). In the study of Tompkins (2001), cloud-resolving simulations using a doubly periodic domain of 90 km, at a horizontal grid spacing of 350 m, and with Fig. 9 A snapshot of the nearsurface equivalent potential temperature distribution from Torri and Kuang (2016b) no imposed mean wind vertical structure, accumulate water vapor within the full outer onethird of the cold pool, by 0.25 g kg À1 in the mean. The moistening remained inside the gust front edge. The anomalous moisture was attributed to the evaporation of precipitation into temperature-recovered air before the downdraft of air from a higher-altitude source could reach the surface, and a straightforward calculation of how much rain could evaporate into the subcloud layer is consistent with the 0.25 g kg À1 increase evident in the composite mean (Langhans and Romps 2015). In this way, rain evaporation not only moistens the air, but it also raises the air's h e (which is otherwise conserved with phase changes), creating a positive feedback in triggering further convection. Surface fluxes (and winds) were discounted because these would be diminishing at the cold pool edge, as also evident in the observational composite (Fig. 5).
More recent modeling experimentation with 'single-bubble' convection further emphasizes the existence and importance of enhanced near-surface moisture resulting from the parent convection, but attributes the source of the moisture primarily to surface fluxes (Langhans and Romps 2015;Romps and Jeevanjee 2016). Rain evaporation nevertheless remains important by contributing the anomalous moisture capable of reaching the cloud base (Torri and Kuang 2016b). Similar to Tompkins (2001), the moisture in both these studies remains inside the cold pool edge.
The observational composites in Fig. 5 do not show widespread enhancements of water vapor mixing ratio q v interior to the cold pool. A slight observed enhancement in q v , of 0.25 g kg À1 , instead, appears at the cold pool edge, and prior to the wind increase. The increase in temperature of $ 0.1 K could conceivably reflect how the temperature drop is identified, but that would not explain the q v increase. Such prior moisture increases were also noted in the tethered balloon measurements analyzed within Addis et al. (1984) and associated with upward vertical velocities. Their presence prior to the wind speed increases points to a convergence of the surface wind, and with it, moisture. The slight temperature increase could also be explained in this way.
The question then arises why models produce water vapor rings (e.g., Langhans and Romps 2015). Model cold pools can be composited similarly to observations, using time series at model grid points in the lowest model layer. Three examples, all based on the System for Atmospheric Modeling (SAM; Khairoutdinov and Randall 2006), are shown in Figs. 10 and 11. The composites shown in Fig. 10 incorporate large-scale forcings driven DYNAMO observations (Wang et al. 2015) and utilize a one-km horizontal grid spacing. Figure 11 show composites based on the Torri and Kuang (2016b) simulations, with an additional simulation performed using a smaller grid spacing (80 vs. 250 m), but otherwise identical. All three simulations lack wind shear.
Similarly to the Tompkins (2001) composites, a water vapor mixing ratio increase is seen to occur within the cold pool frontal edge, before a decrease after the frontal passage. This occurs within all three simulations. The simulations in Fig. 11 do place the water vapor increase closer to the front, but are themselves not sensitive to the horizontal grid spacing. The model q v increases exceed those observed, also evident in Fig. 3 of Feng et al. (2015). These comparisons suggest models may overproduce water vapor rings. Postulated explanations can include an under-entrainment of drier air into the gust front edge or excessive evaporation, reflecting turbulent mixing and microphysical processes that are difficult to represent accurately, or a surface flux-wind feedback that is too strong. In Fig. 11, the simulation with the smaller grid spacing produces surface fluxes that are slightly lower than in the other simulation (102.7 vs. 100.6 W m À2 for the latent heat fluxes, and 9.58 vs. 9.48 W m À2 for the sensible heat fluxes), and may be due to different gust front velocities, while evaporation rates are slightly higher when the grid spacing is smaller. More work is required to confidently explain the differences beyond the scope of this current contribution, with this example primarily introduced to highlight the need for a benchmark dataset of both observations and model simulations.
In addition, model simulations of cold pools generally make trade-offs between domain size and grid spacing. 2 Studies that use larger domain sizes to represent cold pool mesoscale organization are more likely to attribute anomalous moisture sources for secondary convection outside of the cold pool (Li et al. 2014;Schlemmer and Hohenegger 2016). This is also consistent with the observed inhomogeneity of the secondary convection around the mesoscale arc (Figs. 2, 4) locations, where a dynamical lifting of a surface parcel is most able to access moisture, will be the most likely to see convection thrive (see also Torri et al. 2015). Further advantages of larger domains are the ability to capture colliding cold pools (e.g., Fig. 9), and, when incorporating the large-scale forcing at the boundary as opposed to using doubly periodic domains, the ability to capture cold pool asymmetries (Li et al. 2014(Li et al. , 2015. Wind shear aloft can modify the cloud's geometry, thus limiting cloud deepening (Zuidema et al. 2012;Li et al. 2014) but can also increase evaporation (Schlemmer and Hohenegger 2014), as well as allow the precipitation to fall outside of the main updraft core, and has bearing on the overall cloud fraction. (e) Fig. 10 Composites of a, g: air temperature (black) and water vapor mixing ratio (red); b, h water vapor mixing ratio (red) and equivalent potential temperature (black); c, i surface wind speed (black) and surface pressure (red); d, j latent heat fluxes (LHF, black) and the Bowen ratio (surface to latent heat fluxes; red); e, k rain rate (black) and column water vapor (red, e) only); f, l liquid water path, based on SAM simulations.
The left a-f panel applies to Gan island from two time periods combined, October 8-16 and December 18-25, 2011. The right g-l panel applies to R/V Revelle simulations spanning October 2 to November 1, 2011. The doubly periodic simulations at a one-km horizontal grid spacing incorporate daily large-scale forcings developed for the DYNAMO time period from the campaign observations (Wang et al. 2015). 5-min model output was evaluated at 16 grid points within a 256 km by 256 km spatial domain The dominating remaining observational challenge arguably remains the need to better resolve the four-dimensional humidity field. The idea that the updraft air feeding convection tends to have a high h e (or equivalently, moist static energy) is not controversial, either observationally (Kingsmill and Houze 1999) or theoretically (Emanuel et al. 1994). Consensus about the source of the anomalous moisture and equivalent potential energy (h e ) for the secondary convection that is triggered by cold pools has not yet been reached. The distribution of moisture near cloud base remains poorly known, along with the vertical structure of cold pools. This is unlikely to be solved soon via remote sensing (see, e.g., Mapes et al. 2017;Zhang et al. 2017), with what is possibly the best passive remote sensing technique, using spectral infrared (Blumberg et al. 2015) still primarily limited to clear skies. Differential absorption and Raman lidars can profile moisture  but are obstructed by cloud, although extrapolation from points nearby could provide insight, and in combination with Doppler wind lidar might be able to resolve moist updrafts. The measurement of the most important state variables can also be readily done as the standard package on a small research aircraft and now also on drones and unmanned aerial vehicles (e.g., http://vandenheever.atmos.colostate.edu/vdhpage/c3loud-ex/index. php), or using small, light balloons as pseudo-Lagrangian drifters (https://sites.psu.edu/ pmarkowski/2017/06/07/markowski-and-richardson-fly-dozens-of-airborne-probes-into-threesevere-storms-in-oklahoma-and-kansas/). This would shed further light on the relationship of the thermodynamics to the cold pool dynamics as well as conceivably interaction with surface fluxes (Ross et al. 2004;Gentine et al. 2016;Grant and van den Heever 2016). Future challenges also remain on the modeling front. For example, the influence of entrainment at the gust front remains an open question. Work has been done to determine the entrainment in gravity currents (see e.g., Hacker et al. 1996;Hallworth et al. 1996;Fragoso et al. 2013), but mostly in idealized scenarios, which raises questions on the applicability of these results to real-world cold pools. To address this question using numerical models is a challenging task as the simulations required for this purpose would have to be conducted in a large enough domain, and, at the same time, with a spatial grid spacing capable of properly representing turbulent mixing at the gust front (see Grant and van den Heever (2016) for such an attempt). Such efforts will also improve parameterizations of the boundary layer turbulent transports towards representing clear-air entrainment and up/downdrafts correctly (see also Tompkins and Semie 2017).