A 1-km Multi-ICs/BCs and Mixed-Physics WRF Ensemble for Sea-Breeze Convection Prediction during ESCAPE: Precipitation Evaluation
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Lebo, Zachary J..
Explore the source record for details and available documents.
Aerosols serve as cloud condensation nuclei, shaping the microphysical properties of cloud droplets. Aerosol effects on convective clouds are complex and remain controversial. The debate centers around the process of aerosol-induced invigoration of deep convection, a phenomenon that could significantly affect convective cloud properties but lacks robust evidence due to methodological limitations in observational approaches and questions about the robustness of modeling studies. Resolving these discrepancies is crucial for understanding how aerosols affect the atmosphere. Here, this study examines the effects of meteorological and aerosol parameters in a weakly synoptic-driven convective environment, where the influence of aerosols may be more pronounced and observable. Daily atmospheric soundings and aerosol concentrations from several ground instruments collected during the summer of 2022 in Houston, Texas, as part of the Tracking Aerosol Convection interactions Experiment (TRACER) and Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE) field campaigns are analyzed. Statistical learning methods are applied to uncover the complex relationships between aerosols, meteorology, and convective cloud characteristics, such as cell area and echo-top height. The findings reveal that higher aerosol concentrations are associated with narrower convective cells, which we argue contradicts the idea of stronger convection with increased aerosol loading. However, once the data are clustered by the synoptic environment, the relationship between aerosol loading and convective cell area diminishes, indicating that the covariablity between synoptic-scale weather patterns, local thermodynamics, and aerosol loading makes it challenging to draw definitive conclusions about the specific impacts of aerosols on convective cloud properties.
Abstract This study investigates how entrainment’s diluting effect on cumulonimbus updraft buoyancy is affected by the temperature of the troposphere, which is expected to increase by the end of the century. A parcel model framework is constructed that allows for independent variations in the temperature ( T ), the entrainment rate ε , the free-tropospheric relative humidity (RH), and the convective available potential energy (CAPE). Using this framework, dilution of buoyancy is evaluated with T and RH independently varied and with CAPE either held constant or increased with temperature. When CAPE is held constant, buoyancy decreases as T increases, with parcels in warmer environments realizing substantially smaller fractions of their CAPE as kinetic energy (KE). This occurs because the increased moisture difference between an updraft and its surroundings at warmer temperatures drives greater updraft dilution. Similar results are found in midlatitude and tropical conditions when CAPE is increased with temperature. With the expected 6%–7% increase in CAPE per kelvin of warming, KE only increases at 2%–4% K −1 in narrow updrafts but tracks more closely with CAPE at 4%–6% in wider updrafts. Interestingly, the rate of increase in the KE with T becomes larger than that of CAPE when the later quantity increases at more than 10% K −1 . These findings emphasize the importance of considering entrainment in studies of moist convection’s response to climate change, as the entrainment-driven dilution of buoyancy may partially counteract the influence of increases in CAPE on updraft intensity. Significance Statement Cumulonimbus clouds mix air with their surrounding environment through a process called entrainment, which controls how efficiently environmental energy is converted into upward speed in thunderstorm updrafts. Our research shows that warmer temperatures will exacerbate the moisture difference between cumulonimbus updrafts and their surroundings, leading to greater mixing and less efficient conversion of environmental energy into updraft speeds. This effect should be considered in future research that investigates how climate change will affect cumulonimbus clouds.
Explore the source record for details and available documents.
Localized projections of 21st‐century hydroclimate variables obtained from downscaling Global Climate Model (GCM) output are central to informing regional impact assessments and infrastructure planning. Regional GCM biases can be significant and, for dynamical downscaling, can be addressed either before (a priori) or after (a posteriori) downscaling. However, a priori bias correction (APBC) has generally unexplored effects on climate change signals. Here we analyze dynamically downscaled solutions of CMIP6 GCMs over the Western U.S., with and without APBC, and quantify APBC's impact on climate change signals relative to other irreducible uncertainty sources. For temperature and precipitation, the uncertainty introduced by APBC is negligible compared to that arising from GCM choice or internal variability. Furthermore, APBC greatly reduces regional models' unrealistically high snow‐water‐equivalent (SWE) biases that result directly from GCM errors. We leverage this finding to encourage the dynamical downscaling community to adopt APBC as a standard operating procedure.
Polarization surrounding bias correction (BC) in creating climate projections arises from its lack of physicality. Here, we perform and analyze 18 dynamical downscaling simulations (with and without BC) to better understand the physical impacts of BC, applied before downscaling, on regional climate output across the western United States. Without BC, downscaled precipitation is systematically and unrealistically wet biased compared to a hierarchy of observationally based datasets over the 1980–2014 period due to cascading mean–state Global Climate Model (GCM) biases: (a) overly strong lower–tropospheric lapse rates (5 K/km), (b) overly cold (2 K) tropospheric temperatures, and (c) anomalous mid–tropospheric cyclonic vorticity advection. With BC, downscaled precipitation (snow) biases are virtually eliminated (halved). Identified GCM biases are common to the broader Coupled Model Intercomparison Project ensemble. Physical effects of BC on the quality of the regionalized projections, pending an evaluation of BC's distortion of the downscaled climate response, may motivate its broader application by dynamical downscalers.
Abstract. Predicting future climate change over a region of complex terrain, such as the western United States (US), remains challenging due to the low resolution of global climate models (GCMs). Yet the climate extremes of recent years in this region, such as floods, wildfires, and drought, are likely to intensify further as climate warms, underscoring the need for high-quality and high-resolution predictions. Here, we present an ensemble of dynamically downscaled simulations over the western US from 1980–2100 at 9 km grid spacing, driven by 16 latest-generation GCMs. This dataset is titled the Western US Dynamically Downscaled Dataset (WUS-D3). We describe the challenges of producing WUS-D3, including GCM selection and technical issues, and we evaluate the simulations' realism by comparing historical results to temperature and precipitation observations. The future downscaled climate change signals are shaped in physically credible ways by the regional model's more realistic coastlines and topography. (1) The mean warming signals are heavily influenced by more realistic snowpack. (2) Mean precipitation changes are often consistent with wetting on the windward side of mountain complexes, as warmer, moister air masses are uplifted orographically during precipitation events. (3) There are large fractional precipitation increases on the lee side of mountain complexes, leading to potentially significant changes in water resources and ecology in these arid landscapes. (4) Increases in precipitation extremes are generally larger than in the GCMs, driven by locally intensified atmospheric updrafts tied to sharper, more realistic gradients in topography. (5) Changes in temperature extremes are different from what is expected by a shift in mean temperature and are shaped by local atmospheric dynamics and land surface feedbacks. Because of its high resolution, comprehensiveness, and representation of relevant physical processes, this dataset presents a unique opportunity to evaluate societally relevant future changes in western US climate.
Deep convective updraft invigoration via indirect effects of increased aerosol number concentration on cloud microphysics is frequently cited as a driver of correlations between aerosol and deep convection properties. Here, we critically evaluate the theoretical, modeling, and observational evidence for warm- and cold-phase invigoration pathways. Though warm-phase invigoration is plausible and theoretically supported via lowering of the supersaturation with increased cloud droplet concentration in polluted conditions, the significance of this effect depends on substantial supersaturation changes in real-world convective clouds that have not been observed. Much of the theoretical support for cold-phase invigoration depends on unrealistic assumptions of instantaneous freezing and unloading of condensate in growing, isolated updrafts. When applying more realistic assumptions, impacts on buoyancy from enhanced latent heating via fusion in polluted conditions are largely canceled by greater condensate loading. Many foundational observational studies supporting invigoration have several fundamental methodological flaws that render their findings incorrect or highly questionable. Thus, much of the evidence for invigoration has come from numerical modeling, but different models and setups have produced a vast range of results. Furthermore, modeled aerosol impacts on deep convection are rarely tested for robustness, and microphysical biases relative to observations persist, rendering many results unreliable for application to the real world. Without clear theoretical, modeling, or observational support, and given that enervation rather than invigoration may occur for some deep convective regimes and environments, it is entirely possible that the overall impact of cold-phase invigoration is negligible. Substantial mesoscale variability of dominant thermodynamic controls on convective updraft strength coupled with substantial updraft and aerosol variability in any given event are poorly quantified by observations and present further challenges to isolating aerosol effects. Observational isolation and quantification of convective invigoration by aerosols is also complicated by limitations of available cloud condensation nuclei and updraft speed proxies, aerosol correlations with meteorological conditions, and cloud impacts on aerosols. Furthermore, many cloud processes, such as entrainment and condensate fallout, modulate updraft strength and aerosol–cloud interactions, varying with cloud life cycle and organization, but these processes remain poorly characterized. Considering these challenges, recommendations for future observational and modeling research related to aerosol invigoration of deep convection are provided.
Are the results of aerosol invigoration studies that neglect entrainment valid for diluted deep convective clouds? We address this question by applying an entraining parcel model to soundings from tropical and midlatitude convective environments, wherein pollution is assumed to increase parcel condensate retention. Invigoration of 5%–10% and <2% is possible in undiluted tropical and midlatitude parcels respectively when freezing is rapid. This occurs because the positive buoyancy contribution from freezing is larger than the negative buoyancy contribution from condensate loading, leading to positive net condensate contribution to buoyancy. However, aerosol-induced weakening is more likely when realistic entrainment rates occur because water losses from entrainment more substantially reduce the latent heating relative to the loading contribution. This leads to larger net negative buoyancy contribution from condensates in polluted than in clean entraining parcels. Our results demonstrate that accounting for entrainment is critical in conceptual models of aerosol indirect effects in deep convection.
A study of the vertical structure of postfrontal shallow clouds in the marine boundary layer over the Southern Ocean is presented here. The central question of this two-part study regards cloud phase (liquid/ice) of precipitation, and the associated growth mechanisms. In this first part, data from the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign are analyzed, starting with a 75-h case with continuous sea surface-based thermal instability, modest surface heat fluxes, an open-cellular mesoscale organization, and very few ice nucleating particles (INPs). The clouds are mostly precipitating and shallow (tops mostly around 2 km above sea level), with weak up- and downdrafts, and with cloud-top temperatures generally around –18° to –10°C. Further, the case study is extended to three other periods of postfrontal shallow clouds in MARCUS. While abundant supercooled liquid water is commonly present, an experimental cloud-phase algorithm classifies nearly two-thirds of clouds in the 0° to –5°C layer as containing ice (cloud ice, snow, or mixed phase), implying that much of the precipitation grows through cold-cloud processes. The best predictors of ice presence are cloud-top temperature, cloud depth, and INP concentration. Measures of convective activity and turbulence are found to be poor indicators of ice presence in the studied environment. The water-phase distribution in this cloud regime is explored through numerical simulations in Part II. Significance Statement: Climate models generally predict a lower albedo than observed over the Southern Ocean, and this is largely attributed to a lack of cloudiness, especially in the postfrontal cold sector of midlatitude cyclones. This in turn may be due to an excess of ice in these simulated clouds, resulting in rapid precipitation fallout and an overly brief cloud lifespan. The objective of this study is to examine whether shallow postfrontal clouds over the Southern Ocean are dominated by supercooled drops, or by snow and ice, using data collected by a U.S. Department of Energy Atmospheric Radiation Measurement Mobile Facility deployed aboard an Australian Antarctic supply vessel. We find that these clouds contain much supercooled liquid, even though cloud-top temperatures generally are around –18° to –8°C, and that about two-thirds of the clouds just above the freezing level contain ice. Much of the precipitation appears to grow through cold-cloud processes above the freezing level, rather than drizzle/rain. Updrafts and/or turbulence in convection or in cloud-top generating cells do not initiate much ice, compared to observations elsewhere in a similar temperature range. This may be attributable to the extremely low concentration of ice nucleating particles in this environment. Ultimately, the deepest clouds with the coldest cloud tops are most likely to be ice dominated.
Abstract Part I of this series presented a detailed overview of postfrontal mixed-phase clouds observed during the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign. In Part II, we focus on a multiday (23–26 February 2018) case with the aim of understanding ice production as well as model sensitivity to ice process parameterizations using the Weather Research and Forecasting (WRF) Model. The control simulation with the Predicted Particle Properties (P3) microphysics scheme underestimates the ice content and overestimates the supercooled liquid water content, contrary to the bias common in global climate models. The simulations targeted at ice production processes show negligible sensitivity to cloud droplet number concentrations. Further, neither increasing ice nuclei particle (INP) concentrations to an unrealistic level nor adjusting it to MARCUS field estimations alone guarantees more ice production in the model. However, the simulated clouds are found to be highly sensitive to the implementation of immersion freezing, the thresholding of condensation/deposition freezing initiation, and the rime splintering process. By increasing immersion freezing of cloud droplets, relaxing thresholds for condensation/deposition freezing, or removing rime splintering thresholds, the model significantly improves its performance in producing ice. The relaxation of the immersion freezing temperature threshold to the observed cloud-top temperature suggests an in-cloud seeder–feeder mechanism. The results of this work call for an increase in observations of INP, especially over the remote Southern Ocean and at relatively high temperatures, and measurements of ice particle size distributions to better constrain ice nucleating processes in models.
Explore the source record for details and available documents.
Abstract This paper presents incorporation of the microphysical piggybacking into the Weather Research and Forecasting (WRF) model. Microphysical piggybacking is to run a single simulation applying two microphysical schemes, the first scheme driving the simulation and the second piggybacking this simulated flow. “Driving the simulation” means that the simulated microphysical processes, affect the cloud buoyancy and thus force the simulated flow. In contrast, the piggybacking variables are advected by the simulated flow and undergo microphysical transformation, but they do not affect the simulated flow (like in prescribed flow—kinematic—simulations). The two sets of variables (driver and piggybacker) include temperature, water vapor mixing ratio, and all microphysical variables. We provide details of implementing piggybacking into the WRF model, illustrate its applications, and demonstrate the benefits of this methodology in two idealized three‐dimensional cases: (a) a squall line case applying two microphysics schemes, the Thompson bulk microphysics scheme and the University of Pécs/NCAR bin (UPNB) scheme. The piggybacking simulations revealed that the microphysics‐dynamics interaction plays a more important role than the pure microphysical size sorting effect in the transition zone formation. (b) A case of daytime shallow‐to‐deep convective development over land. This case uses the UPNB scheme and contrasts convection developing in environments with either pristine or polluted cloud condensation nuclei (CCN). The piggybacking results indicated that the increase of cloud cover and decrease of supersaturation are mainly associated with the microphysical effect of increasing CCN while the change of precipitation on the ground is also influenced by microphysics‐dynamics interactions.
Abstract Bulk ice-microphysical models parameterize the dynamic evolution of ice particles from advection, collection, and sedimentation through a cloud layer to the surface. Frozen hydrometeors can grow to acquire a multitude of shapes and sizes, which influence the distribution of mass within cloud systems. Aggregates, defined herein as the collection of ice particles, have a variety of formations based on initial ice particle size, shape, falling orientation, and the number of particles that collect. This work focuses on using the Ice Particle and Aggregate Simulator (IPAS) as a statistical tool to repetitively collect ice crystals of identical properties to derive bulk aggregate characteristics. A database of 9 744 000 aggregates is generated with resulting properties analyzed. After 150 single ice crystals (monomers) collect, the most extreme aggregate aspect ratio calculations asymptote toward and ϕ ca ≈ 0.50 for aggregates composed of quasi-horizontally oriented and randomly oriented monomers, respectively. The results presented are largely consistent with both a previous theoretical study and estimates derived from ground-based observations from two different geographic locations. Particle falling orientation highly influences newly formed aggregate aspect ratios from the collection of particles with extreme aspect ratios; quasi-horizontally oriented particles can produce aggregate aspect ratios an order of magnitude more extreme than randomly oriented particles but can also produce near-spherical aggregates as the number of monomers comprising the aggregate reach approximately 100. Finally, a majority of collections result in aggregates that are closer to prolate than oblate spheroids.
Abstract The Ice Particle and Aggregate Simulator (IPAS) is used to theoretically represent the aggregation process of ice crystals. Aggregates have a variety of formations based on initial ice particle size, shape, and falling orientation, all of which influence water phase partitioning. Aggregate dimensional properties and density changes are calculated for monomer–monomer (MON–MON), monomer–aggregate (MON–AGG), and aggregate–aggregate (AGG–AGG) collection to be used by ice-microphysical models for improvement in aggregation parameterizations. Aggregates are chosen from a database of 9 744 000 preformed combinations to be further collected (see Part II). AGG–AGG collection results in more extreme and a smaller range of aggregate aspect ratios than MON–AGG collection. A majority of aggregates are closer to prolate than oblate spheroids, regardless of collection type, except for quasi-horizontally oriented particles that have extreme aspect ratios to begin with. MON–AGG collection frequently results in an increase in density upon collection, whereas MON–MON and AGG–AGG collection almost always result in particle density decreases, with extreme reductions near 99% for MON–MON collection. MON–MON collection results in the greatest decreases in density but then quickly becomes unaffected by the addition of more monomers due to inherent size differences between monomers and aggregates. Finally, a holistic analysis to in situ observations of cloud particle images is presented. IPAS 2D aspect ratios surround a median value of 0.6 and closely follow that of previous studies while varying by no more than ≈12% on average from observed aggregates.
Explore the source record for details and available documents.
Explore the source record for details and available documents.