Cloud-resolving Model Simulations Research Articles

Fitting Cumulus Cloud Size Distributions From Idealized Cloud Resolving Model Simulations

AbstractWhereas it is now widely accepted that cumulus cloud sizes are power‐law distributed, characteristic exponents reported in the literature vary greatly, generally taking values between 1 and >3. Although these differences might be explained by variations in environmental conditions or physical processes organizing the cloud ensembles, the use of improper fitting methods may also introduce large biases. To address this issue, we propose to use a combination of maximum likelihood estimation and goodness‐of‐fit tests to provide more robust power‐law fits while systematically identifying the size range over which these fits are valid. The procedure is applied to cloud size distributions extracted from two idealized high‐resolution simulations displaying different organization characteristics. Overall, power‐laws are found to be outperformed by alternative distributions in almost all situations. When clouds are identified based on a condensed water path threshold, using power‐laws with an exponential cutoff yields the best results as it provides superior fits in the tail of the cloud size distributions. For clouds identified using a combination of water content and updraft velocity thresholds in the free troposphere, no substantial improvement over pure power‐laws can be found when considering more complex two‐parameter distributions. In this context however, exponential distributions provide results that are as good as, if not better than power‐laws. Finally, it is demonstrated that the emergence of scale free behaviors in cloud size distributions is related to exponentially distributed cloud cores merging as they are brought closer to each other by underlying organizing mechanisms.

What Controls the Mesoscale Variations in Water Isotopic Composition Within Tropical Cyclones and Squall Lines? Cloud Resolving Model Simulations in Radiative‐Convective Equilibrium

AbstractWater isotopes are tracers of convective processes and are often used as proxies for past precipitation. These applications require a better understanding of the impact of convective processes on the isotopic composition of water vapor and precipitation. One way to advance this understanding is to analyze the isotopic mesoscale variations during organized convective systems such as tropical cyclones or squall lines. The goal of this study is to understand these isotopic mesoscale variations with particular attention to isotopic signals in near‐surface vapor and precipitation that may be present in observations and in paleoclimate proxies. With this aim, we run cloud resolving model simulations in radiative‐convective equilibrium in which rotation or wind shear is added, allowing us to simulate tropical cyclones or squall lines. The simulations capture the robust aspects of mesoscale isotopic variations in observed tropical cyclones and squall lines. We interpret these variations using a simple water budget model for the sub‐cloud layer of different parts of the domain. We find that rain evaporation and rain‐vapor diffusive exchanges are the main drivers of isotopic depletion within tropical cyclones and squall lines. Horizontal advection spreads isotopic anomalies, thus reshaping the mesoscale isotopic pattern. This study contributes to our understanding of mesoscale isotopic variability and provides physical arguments supporting the interpretation of paleoclimate isotopic archives in tropical regions in terms of past cyclonic activity.

How Accurately Can Warm Rain Realistically Be Retrieved with Satellite Sensors? Part II: Horizontal and Vertical Heterogeneities

Abstract There are many sources of uncertainty in satellite precipitation retrievals of warm rain. In this paper, the second of a two-part study, we focus on uncertainties related to spatial heterogeneity and surface clutter. A cloud-resolving model simulation of warm, shallow clouds is used to simulate satellite observations from three theoretical satellite architectures—one similar to the Global Precipitation Measurement Core Observatory, one similar to CloudSat, and one similar to the planned Atmosphere Observing System (AOS). Rain rates are then retrieved using a common optimal estimation framework. For this case, retrieval biases due to nonuniform beamfilling are very large, with retrieved rain rates negatively (low) biased by as much as 40%–50% (depending on satellite architecture) at 5 km horizontal resolution. Surface clutter also acts to negatively bias retrieved rain rates. Combining all sources of uncertainty, the theoretical AOS satellite is found to outperform CloudSat in terms of retrieved surface rain rate, with a bias of −19% as compared with −28%, a reduced spread of retrieval errors, and an additional 17.5% of cases falling within desired uncertainty limits. The results speak to the need for additional high-resolution modeling simulations of warm rain so as to better characterize the uncertainties in satellite precipitation retrievals.

Ice-nucleating particles near two major dust source regions

Abstract. Mineral dust and sea spray aerosol represent important sources of ice-nucleating particles (INPs), the minor fraction of aerosol particles able to trigger cloud ice crystal formation and, consequently, influence multiple climate-relevant cloud properties including lifetime, radiative properties and precipitation initiation efficiency. Mineral dust is considered the dominant INP source in many parts of the world due to its ice nucleation efficiency and its sheer abundance, with global emission rates of up to 4700 Tg a−1. However, INPs emitted from the ocean surface in sea spray aerosol frequently dominate INP populations in remote marine environments, including parts of the Southern Ocean where cloud-resolving model simulations have demonstrated that cloud radiative properties are likely strongly controlled by INPs. Here we report INP concentrations measured in aerosol and seawater samples during Air Quality and Climate Change in the Arabian Basin (AQABA), a shipborne campaign that spanned the Red Sea, Gulf of Aden, Arabian Sea, Arabian Gulf and part of the Mediterranean. In aerosol samples collected within a few hundred kilometers of the first and second ranked sources of dust globally, the Sahara and Arabian Peninsula, INP concentrations ranged from 0.2 to 11 L−1 at −20 ∘C with observed ice-active surface site densities (ns) 1–3 orders of magnitude below levels predicted by mineral dust INP parameterizations. Over half of the samples (at least 14 of 26) were collected during dust storms with average dust mass concentrations between 150 and 490 µg m−3 (PM10), as simulated by the Modern-Era Retrospective analysis for Research and Application, version 2 (MERRA-2). The impacts of heat and peroxide treatments indicate that organics dominated the observed ice nucleation (IN) activity at temperatures ≥ −15 ∘C with proteinaceous (heat-labile) INPs frequently observed at high freezing temperatures >−10 ∘C. INP concentrations in seawater samples ranged between 3 and 46 mL−1 at −19 ∘C, demonstrating the relatively low INP source potential of seawater in the region as compared to seawater from multiple other regions reported previously. Overall, our results demonstrate that despite proximity to the Sahara and the Arabian Peninsula and the dominance of mineral dust in the aerosol sampled, existing mineral dust parameterizations alone would not skillfully represent the near-surface ns in the observed temperature regime (−6 to −25 ∘C). Future efforts to develop or improve representations of dust INPs at modest supercooling (≥-15 ∘C) would benefit from a characterization of the specific organic species associated with dust INPs. More generally, an improved understanding of the organic species associated with increased IN activity and their variability across dust source regions would directly inform efforts to determine whether ns-based parameterizations are appropriate for faithful representation of dust INPs in this sensitive temperature regime, whether region-specific parameterizations are required, or whether an alternative to the ns approach is necessary.

Dependence of Convective Cloud Properties and Their Transport on Cloud Fraction and GCM Resolution Diagnosed from a Cloud-Resolving Model Simulation

The scale-aware convective parameterization for high resolution global climate models must satisfy the requirement that the parameterized subgrid convective transport diminishes as the model resolution increases to convection-resolving resolutions. A major assumption in current scale-aware convection schemes is that the differences between convective cloud properties and their environmental counterparts are independent of cloud fraction. This study examines convective cloud vertical velocity, moist static energy (MSE), moisture, and the vertical eddy transport of MSE and moisture for different averaging subdomain sizes and fractional convective cloudiness using a cloud resolving model simulation of a midlatitude mesoscale convective system. Results show that convective cloud fraction, mass flux, and vertical transport of MSE and moisture increase with decreasing subdomain size. The differences between convective cloud properties in both updrafts and downdrafts and their environment depend on both cloud fraction and the averaging subdomain size. For a given subdomain size, the differences increase with cloud fraction, in contrast to the assumption used in current scale-aware convection parameterization schemes. A consequence of this is that the parameterized convective eddy transport reaches maximum at a higher cloud fraction than believed in previous studies. This has implications on how fast the subgrid convective transport should diminish as GCM resolution increases.

Investigation of aerosol effects on a deep convection case in Southeast China with budget analysis of hydrometeor mass and latent heat

It remains an open question which microphysical processes determine the sign of aerosol effects on deep convective rain. This lack of knowledge prevents us from understanding the uncertainties of the aerosol effects on deep convection. We try to answer this question through a cloud-resolving model (CRM) simulation on a summer storm that occurred in Southeast China, by analyzing both column hydrometeor mass and latent heat budget and quantifying the contribution of each microphysical process to the rain variations with aerosol. Our results show that warm-phase microphysical processes have stronger contributions (52% −99%) to the rain variations than ice-phase microphysical processes (1%−48%) when gradually increasing cloud condensation nuclei (CCN) concentration from 125 cm−3 (clean environment) to 16,000 cm−3 (highly polluted environment). Condensation and evaporation of cloud droplets are identified as two crucial microphysical processes, both of which are enhanced with increased aerosol loading. The competition between them (i.e., net condensation) essentially determines the sign of aerosol effects on surface rain. We also explore aerosol effects on updraft velocity from the perspective of buoyancy and latent heat budget, it reveals that the largely enhanced convective updrafts below 7 km are attributed to the increase of condensation heating. While the competition between enhanced condensation heating and suppressed deposition heating restricts the aerosol-induced convective invigoration above 7 km. This study highlights the dominant role of net condensation in the aerosol-imposed invigoration of deep convection under highly polluted environments, suggesting that an explicit treatment of vapor diffusion processes is necessary to simulate aerosol effects on deep convection more accurately.

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

AbstractWe design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

Assessing free tropospheric quasi-equilibrium for different GCM resolutions using a cloud-resolving model simulation of tropical convection

This study examines the free-tropospheric quasi-equilibrium at different global climate model (GCM) resolutions using the simulation of tropical convection by a cloud-resolving model during the Tropical Western Pacific International Cloud Experiment. The simulated dynamic and thermodynamic fields within the model domain are averaged over subdomains of different sizes equivalent to different GCM resolutions. These coarse-grained fields are then used to compute CAPE and its change with time, and their relationships with simulated convection. Results show that CAPE change with time is controlled predominantly by variations of thermodynamic properties in the planetary boundary layer for all subdomain sizes ranging from 64 to 4 km. Lag correlation analysis shows that CAPE generation by the free-tropospheric dynamical advection (dCAPEls) leads convective precipitation but is in phase with convective mass flux at 600 mb and 500 mb vertical velocity for all subdomain sizes. However, the correlation coefficients and regression slopes decrease as the subdomain size decreases for subdomain sizes smaller than 16 km. This is probably due to increased randomness of convection and more scale-dependence of the relationships when the subdomain size reaches the grey zone. By examining the sensitivity of the relationships of convection with dCAPEls to temporal scales in different subdomain size, it shows that the quasi-equilibrium between dCAPEls and convection holds well for timescales of 30 min or longer at all subdomain sizes. These results suggest that the free tropospheric quasi-equilibrium assumption may still be useable even for GCM resolutions in the grey zone.

Moisture transport enhanced by the nocturnal low-level jet in association with the passage of a monsoon depression over the Indian subcontinent

Characteristics of moisture transport associated with a monsoon depression (MD) formed over the Bay of Bengal and later moved westward over central India are presented. The investigation is carried out mainly using the Weather Research and Forecast (WRF) cloud-resolving model simulations and the vertically integrated water vapor transport (IVT) to assess the path and features of the moisture transport. The model performance evaluation is carried out using micrometeorological and microwave radiometer observations from the Cloud Aerosol Interaction and Precipitation Enhancement Experiment (CAIPEEX). The results indicate that the strengthened LLJ during the nocturnal periods helps to enhance the moisture convergence in the lower troposphere and thus the rainfall, both before the formation of MD as well as during its passage over land. The moisture budget over a selected region along the path of the MD gave further insights on the moisture convergence, associated precipitation and related features during the study period. The analyses using IVT and LLJ strength are illustrated to be useful in identifying areas of intense convection.

Diurnal Differences in Tropical Maritime Anvil Cloud Evolution

AbstractSatellite observations of tropical maritime convection indicate an afternoon maximum in anvil cloud fraction that cannot be explained by the diurnal cycle of deep convection peaking at night. We use idealized cloud-resolving model simulations of single anvil cloud evolution pathways, initialized at different times of the day, to show that tropical anvil clouds formed during the day are more widespread and longer lasting than those formed at night. This diurnal difference is caused by shortwave radiative heating, which lofts and spreads anvil clouds via a mesoscale circulation that is largely absent at night, when a different, longwave-driven circulation dominates. The nighttime circulation entrains dry environmental air that erodes cloud top and shortens anvil lifetime. Increased ice nucleation in more turbulent nighttime conditions supported by the longwave cloud-top cooling and cloud-base heating dipole cannot compensate for the effect of diurnal shortwave radiative heating. Radiative–convective equilibrium simulations with a realistic diurnal cycle of insolation confirm the crucial role of shortwave heating in lofting and sustaining anvil clouds. The shortwave-driven mesoscale ascent leads to daytime anvils with larger ice crystal size, number concentration, and water content at cloud top than their nighttime counterparts.Significance StatementDeep convective activity and rainfall peak at night over the tropical oceans. However, anvil clouds that originate from the tops of deep convective clouds reach their largest extent in the afternoon hours. We study the underlying physical mechanisms that lead to this discrepancy by simulating the evolution of anvil clouds with a high-resolution model. We find that the absorption of sunlight by ice crystals lofts and spreads the daytime anvil clouds over a larger area, increasing their lifetime, changing their properties, and thus influencing their impact on climate. Our findings show that it is important not only to simulate the correct onset of deep convection but also to correctly represent anvil cloud evolution for skillful simulations of the tropical energy balance.

Influence of Upper-Troposphere Stratification and Cloud-Radiation Interaction on Convective Overshoots in the Tropical Tropopause Layer

AbstractIt is still debated whether radiative heating observed in the tropical tropopause layer (TTL) is balanced primarily by cooling from convective overshoots, as in an entrainment layer, or by adiabatic cooling from large-scale eddy-driven upwelling. In this study, three-dimensional cloud-resolving model simulations of radiative-convective equilibrium were carried out with three different cloud microphysics schemes and 1-km horizontal resolution. We demonstrate that overshooting cooling in the TTL can be strongly modulated by upper-troposphere stratification. Two of the schemes produce a hard-landing scenario in which convective overshoots reach the TTL with frequent large vertical velocity leading to strong overshooting cooling (~ −0.2 K day-1). The third scheme produces a soft-landing scenario in which convective overshoots rarely reach the TTL with large vertical velocity and produce little overshooting cooling (~ −0.03 K day-1). The difference between the two scenarios is attributed to changes in the upper-troposphere stratification related to different atmospheric cloud radiative effects (ACRE). The microphysics scheme that produces the soft-landing scenario has much stronger ACRE in the upper troposphere leading to a ~3K warmer and more stable layer which acts as a buffer zone to slow down the convective updrafts. The stratification mechanism suggests the possibility for the ozone variation or eddy-driven upwelling in the TTL to modulate convective overshoots. We further test the sensitivity of overshooting cooling to changes in model resolution by increasing the horizontal resolution to 100 m. The corresponding change of overshooting cooling is much smaller compared with the difference between the hard-landing and soft-landing scenarios.

Transport and chemistry of isoprene and its oxidation products in deep convective clouds

Deep convective clouds can transport trace gases from the planetary boundary layer into the upper troposphere where subsequent chemistry may impact aerosol particle formation and growth. In this modelling study, we investigate processes that affect isoprene and its oxidation products injected into the upper troposphere by an isolated deep convective cloud in the Amazon. We run a photochemical box model with coupled cloud microphysics along hundreds of individual air parcel trajectories sampled from a cloud-resolving model simulation of a convective event. The box model simulates gas-phase chemical reactions, gas scavenging by liquid and ice hydrometeors, and turbulent dilution inside a deep convective cloud. The results illustrate the potential importance of gas uptake to anvil ice in regulating the intensity of the isoprene oxidation and associated low volatility organic vapour concentrations in the outflow. Isoprene transport and fate also depends on the abundance of lightning-generated nitrogen oxide radicals (NOx = NO + NO2). If gas uptake on ice is efficient and lightning activity is low, around 30% of the boundary layer isoprene will survive to the cloud outflow after approximately one hour of transport, while all the low volatile oxidation products will be scavenged by the cloud hydrometeors. If lightning NOx is abundant and gas uptake by ice is inefficient, then all isoprene will be oxidised during transport or in the immediate outflow region, while several low volatility isoprene oxidation products will have elevated concentrations in the cloud outflow. Reducing uncertainties associated with the uptake of vapours on ice hydrometeors, especially HO2 and oxygenated organics, is essential to improve predictions of isoprene and its oxidation products in deep convective outflows and their potential contribution to new particle formation and growth.

On the Semidiurnal Variation in Surface Rainfall Rate over the Tropics in a Global Cloud-Resolving Model Simulation and Satellite Observations

On the Semidiurnal Variation in Surface Rainfall Rate over the Tropics in a Global Cloud-Resolving Model Simulation and Satellite Observations

Formation and maintenance of subsiding shells around non‐precipitating and precipitating cumulus clouds

AbstractA shell of subsiding air is generally known to develop around cumulus clouds and shield them from their environment. We seek here to improve our understanding of such shells by (a) revealing the detailed vertical and horizontal structure of shells surrounding both shallow and deeper clouds, and (b) identifying the mechanisms responsible for in‐shell subsidence generation and maintenance. To that end, a high‐resolution Cloud Resolving Model simulation of the shallow‐to‐deep convection transition over a tropical land surface is analysed with an emphasis on the cloud's near environment. Shells surrounding shallow and deep clouds are found to possess surprisingly similar characteristics. However important differences are observed near cloud top where the deepest clouds are associated with stronger subsidence and broader shells. In the convective outflow region, stronger in‐shell subsidence coincides with strong buoyancy reversal, but also with strong pressure gradients naturally generated by cloud‐top vortex dynamics. A more delicate balance between various processes takes place below, and in‐shell subsidence is only barely sustained as buoyancy reversal is largely compensated by pressure gradients. Finally, while evaporation is clearly the main source of buoyancy reversal everywhere around cloud edges, it is also shown that the downward transport of warmer air from aloft through the subsiding shells may compensate for evaporative cooling to slowly bring in‐shell buoyancy to a near‐neutral state. Overall, while it cannot be denied that evaporative cooling and buoyancy reversal play important roles in generating and sustaining in‐shell subsidence, the present results also emphasise that mechanical forcing at cloud top and downward transport within the shells should not be overlooked.

The NASA‐JAXA Global Precipitation Measurement mission – part II: New frontiers in precipitation science

The <scp>NASA‐JAXA</scp> Global Precipitation Measurement mission – part <scp>II</scp>: New frontiers in precipitation science

The Impact of Warm and Moist Airmass Perturbations on Arctic Mixed-Phase Stratocumulus

AbstractThe Arctic is known to be particularly sensitive to climate change. This Arctic amplification has partially been attributed to poleward atmospheric heat transport in the form of airmass intrusions. Locally, such airmass intrusions can introduce moisture and temperature perturbations. The effect of airmass perturbations on boundary layer and cloud changes and their impact on the surface radiative balance has received increased attention, especially over sea ice with regard to sea ice melt. Utilizing cloud-resolving model simulations, this study addresses the impact of airmass perturbations occurring at different altitudes on stratocumulus clouds for open-ocean conditions. It is shown that warm and moist airmass perturbations substantially affect the boundary layer and cloud properties, even for the relatively moist environmental conditions over the open ocean. The cloud response is driven by temperature inversion adjustments and strongly depends on the perturbation height. Boundary layer perturbations weaken and raise the inversion, which destabilizes the lower troposphere and involves a transition from stratocumulus to cumulus clouds. In contrast, perturbations occurring in the lower free troposphere lead to a lowering but strengthening of the temperature inversion, with no impact on cloud fraction. In simulations where free-tropospheric specific humidity is further increased, multilayer mixed-phase clouds form. Regarding energy balance changes, substantial surface longwave cooling arises out of the stratocumulus break-up simulated for boundary layer perturbations. Meanwhile, the net surface longwave warming increases resulting from thicker clouds for airmass perturbations occurring in the lower free troposphere.

Mind the gap – Part 2: Improving quantitative estimates of cloud and rain water path in oceanic warm rain using spaceborne radars

Abstract. The intrinsic small spatial scales and low-reflectivity structure of oceanic warm precipitating clouds suggest that millimeter spaceborne radars are best suited to providing quantitative estimates of cloud and rain liquid water paths (LWPs). This assertion is based on their smaller horizontal footprint; high sensitivities; and a wide dynamic range of path-integrated attenuations associated with warm-rain cells across the millimeter wavelength spectrum, with diverse spectral responses to rain and cloud partitioning. State-of-the-art single-frequency radar profiling algorithms of warm rain seem to be inadequate because of their dependence on uncertain assumptions about the rain–cloud partitioning and because of the rain microphysics. Here, high-resolution cloud-resolving model simulations for the Rain in Cumulus over the Ocean field study and a spaceborne forward radar simulator are exploited to assess the potential of existing and future spaceborne radar systems for quantitative warm-rain microphysical retrievals. Specifically, the detrimental effects of nonuniform beam filling on estimates of path-integrated attenuation (PIA), the added value of brightness temperature (TB) derived adopting radiometric radar modes, and the performances of multifrequency PIA and/or TB combinations when retrieving liquid water paths partitioned into cloud (c-LWPs) and rain (r-LWPs) are assessed. Results show that (1) Ka- and W-band TB values add useful constraints and are effective at lower LWPs than the same-frequency PIAs; (2) matched-beam combined TB values and PIAs from single-frequency or multifrequency radars can significantly narrow down uncertainties in retrieved cloud and rain liquid water paths; and (3) the configuration including PIAs, TB values and near-surface reflectivities for the Ka-band–W-band pairs in our synthetic retrieval can achieve an RMSE of better than 30 % for c-LWPs and r-LWPs exceeding 100 g m−2.

Potential Vorticity Mixing and Rapid Intensification in the Numerically Simulated Supertyphoon Haiyan (2013)

Abstract The inner-core dynamics of Supertyphoon Haiyan (2013) undergoing rapid intensification (RI) are studied with a 2-km-resolution cloud-resolving model simulation. The potential vorticity (PV) field in the simulated storm reveals an elliptical and polygonal-shaped eyewall at the low and middle levels during RI onset. The PV budget analysis confirms the importance of PV mixing at this stage, that is, the asymmetric transport of diabatically generated PV to the storm center from the eyewall and the ejection of PV filaments outside the eyewall. We employ a piecewise PV inversion (PPVI) and an omega equation to interpret the model results in balanced dynamics. The omega equation diagnosis suggests eye dynamical warming is associated with the PV mixing. The PPVI indicates that PV mixing accounts for about 50% of the central pressure fall during RI onset. The decrease of central pressure enhances the boundary layer (BL) inflow. The BL inflow leads to contraction of the radius of the maximum tangential wind (RMW) and the formation of a symmetric convective PV tower inside the RMW. The eye in the later stage of the RI is warmed by the subsidence associated with the convective PV towers. The results suggest that the pressure change associated with PV mixing, the increase of the symmetric BL radial inflow, and the development of a symmetric convective PV tower are the essential collaborating dynamics for RI. An experiment with 500-m resolution shows that the convergence of BL inflow can lead to an updraft magnitude of 20 m s−1 and to a convective PV tower with a peak value of 200 PVU (1 PVU = 10−6 K kg−1 m2 s−1).

Polarimetric Radar Characteristics of Simulated and Observed Intense Convective Cores for a Midlatitude Continental and Tropical Maritime Environment

AbstractThis study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivityZdrand larger specific differential phaseKdpabove the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with theZdrandKdpobservations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describingZdrandKdpof ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

Tobac 1.2: towards a flexible framework for tracking and analysis of clouds in diverse datasets

Abstract. We introduce tobac (Tracking and Object-Based Analysis of Clouds), a newly developed framework for tracking and analysing individual clouds in different types of datasets, such as cloud-resolving model simulations and geostationary satellite retrievals. The software has been designed to be used flexibly with any two- or three-dimensional time-varying input. The application of high-level data formats, such as Iris cubes or xarray arrays, for input and output allows for convenient use of metadata in the tracking analysis and visualisation. Comprehensive analysis routines are provided to derive properties like cloud lifetimes or statistics of cloud properties along with tools to visualise the results in a convenient way. The application of tobac is presented in two examples. We first track and analyse scattered deep convective cells based on maximum vertical velocity and the three-dimensional condensate mixing ratio field in cloud-resolving model simulations. We also investigate the performance of the tracking algorithm for different choices of time resolution of the model output. In the second application, we show how the framework can be used to effectively combine information from two different types of datasets by simultaneously tracking convective clouds in model simulations and in geostationary satellite images based on outgoing longwave radiation. The tobac framework provides a flexible new way to include the evolution of the characteristics of individual clouds in a range of important analyses like model intercomparison studies or model assessment based on observational data.