Search

Cloud microphysical processes are an important facet of atmospheric modeling, as they can control the initiation and rates of snowfall. Thus, parameterizations of these processes have important implications for modeling seasonal snow accumulation. We conduct experiments with the Weather Research and Forecasting (WRF V4.3.3) Model using three different microphysics parameterizations, including a sophisticated new scheme (ISHMAEL). Simulations are con-ducted for two cold seasons (2018 and 2019) centered on the Colorado Rockies’;750-km2 East River watershed. Precipitation efficiencies are quantified using a drying-ratio mass budget approach and point evaluations are performed against three NRCS SNOTEL stations. Precipitation and meteorological outputs from each are used to force a land surface model (Noah-MP) so that peak snow accumulation can be compared against airborne snow lidar products. We find that micro-physical parameterization choice alone has a modest impact on total precipitation on the order of 63% watershed-wide, and as high as 15% for certain regions, similar to other studies comparing the same parameterizations. Precipitation biases evaluated against SNOTEL are 15% 6 13%. WRF Noah-MP configurations produced snow water equivalents with good correlations with airborne lidar products at a 1-km spatial resolution: Pearson’s r values of 0.9, RMSEs between 8 and 17 cm, and percent biases of 3%–15%. Noah-MP with precipitation from the PRISM geostatistical precipitation product leads to a peak SWE underestimation of 32% in both years examined, and a weaker spatial correlation than the WRF configurations. We fall short of identifying a clearly superior microphysical parameterization but conclude that snow lidar is a valuable nontraditional indicator of model performance.

Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement: Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months. [ABSTRACT FROM AUTHOR]

Different turbulent entrainment‐mixing mechanisms between clouds and environment are essential to cloud‐related processes; however, accurate representation of entrainment‐mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment‐mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the "black box" scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment‐mixing. Specifically, number concentration after entrainment‐mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment‐mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment‐mixing is negatively correlated with liquid water content affected by entrainment‐mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment‐mixing in weather/climate models. Plain Language Summary: Entrainment‐mixing processes occurring between clouds and environmental air have significant effects on cloud‐climate feedback, precipitation, and radiative transfer. Accurately representing these processes has been challenging with previously proposed parameterizations. Machine learning (ML) excels at identifying complex nonlinear relationships and avoids the inherent limitations of conventional parameterizations. Thus, we explore the use of ML to acquire the number concentration and cloud droplet spectral dispersion affected by entrainment‐mixing processes. Simulation data from Explicit Mixing Parcel Model are employed to train, validate, and test the four MLs, including Light Gradient Boosting Machine (LGB), eXtreme Gradient Boosting, Random Forest and Support Vector Regression. After evaluation, the LGB is shown to obtain the most accurate microphysics among the four ML schemes and traditional parameterizations. Additionally, the interpreting method for peeking inside the "black box" reveals that a smaller number concentration and a larger spectral dispersion are indicative of the more inhomogeneous entrainment‐mixing. The relative correlations between predicted microphysics and various features align with our expectations and the underlying physical principles. Sensitivity tests further confirm that incorporating more features produces a more robust and efficient prediction. Overall, this study affirms the reliability and applicability of ML to develop/replace subgrid parameterizations in actual weather/climate models. Key Points: Machine learning schemes are trained to predict cloud microphysical properties affected by turbulent entrainment‐mixing processesThe proposed Light Gradient Boosting Machine provides more accurate microphysical properties compared with traditional parameterization schemesThe partial dependencies of microphysics on features prove a robust alignment with physical mechanisms through the interpreting method [ABSTRACT FROM AUTHOR]

Numerous Earth system models exhibit excessive aerosol effective forcing at the top of the atmosphere (TOA), including the Department of Energy's Energy Exascale Earth System Model (E3SM). Here, in the context of the E3SM version 3 effort, the predicted particle property (P3) stratiform cloud microphysics scheme and an enhanced deep convection parameterization suite (ZM_plus) are implemented into E3SM. The ZM_plus includes a convective cloud microphysics scheme, a multi‐scale coherent structure parameterization for mesoscale convective systems, and a revised cloud base mass flux formulation considering impacts of the large‐scale environment. The P3 scheme improved cloud and radiation particularly over the Northern Hemisphere and the frequency of heavy precipitation over the tropics, and the ZM_plus improved clouds in the tropics. P3 decreases aerosol effective forcing by 0.15 W m−2, while the ZM_plus increases it by 0.27 W m−2, resulting from excessive direct (0.31 W m−2) and indirect forcing (−1.79 W m−2). The excessive aerosol forcings are due to aerosol overestimation associated with insufficient aerosol wet removal. By improving the physical treatments in the aerosol wet removal, we effectively mitigate anthropogenic aerosol overestimation and thus attenuate direct (0.09 W m−2) and indirect aerosol forcing (−1.52 W m−2). Adjustment to primary organic matter hygroscopicity reduces direct and indirect forcing to more reasonable values: −0.13 W m−2 and −1.31 W m−2, respectively. On climatology, improved aerosol treatments mitigate overestimation of aerosol optical depth. Plain Language Summary: The Energy Exascale Earth System Model (E3SM) exhibits strong direct and indirect aerosol forcings, after advanced stratiform and convective cloud treatments are implemented. This study identifies the primary cause of these excessive aerosol forcings as the significant overestimation of anthropogenic aerosols due to insufficient removal of aerosols by precipitation. To address this issue, we made aerosol wet removal representation more physical, which effectively reduced the overestimation of aerosols, bringing both direct and indirect forcings to the expected ranges. Furthermore, cloud and aerosol climatology are notably improved as the result of these developments in cloud and aerosol treatments. Key Points: The new cloud microphysics scheme P3 improves simulations of cloud properties in the NH and aerosol forcingInsufficient wet removal leads to an overly strong aerosol forcing after incorporating the enhanced deep convection parameterization suiteThe aerosol wet removal improvements lead to more reasonable aerosol forcing and improved aerosol climatology in Energy Exascale Earth System Model [ABSTRACT FROM AUTHOR]

Snowpack melting is a crucial water resource for local ecosystems, agriculture, and hydropower in the Intermountain West of the United States. Glaciogenic seeding, a method widely used in mountain regions to enhance precipitation, has been subject to numerous field studies aiming to understand and validate this mechanism. However, investigating precipitation distribution and amounts in mountainous areas is complicated due to the intricate interplay of synoptic circulation patterns and local complex topography. These interactions significantly influence microphysical processes, ultimately affecting the amount and distribution of surface precipitation. To address these challenges, this study leverages Weather Research and Forecasting (WRF) Model simulations, providing high-resolution (900 m), hourly data, spanning the Payette region of Idaho from January to March 2017. We applied the self-organizing map approach to categorize the most representative synoptic circulation patterns and conducted a multiscale analysis to explore their associated environmental conditions and microphysical processes, aiming to assess the cloud seeding potential. The analysis identified four primary synoptic patterns: cold zonal flow (CZF), cold southwesterly flow (CSWF), warm zonal flow (WZF), and warm southwesterly flow (WSWF), constituting 21.3%, 23.1%, 30.0%, and 25.5%, respectively. CSWF and WSWF demonstrated efficiency in generating natural precipitation. These patterns were characterized by abundant supercooled liquid water (SLW) and ice particles, facilitating cloud droplet growth through seeder–feeder processes. On the other hand, CZF exhibited the least SLW and limited potential for cloud seeding, while WZF displayed a lower ice water content but substantial SLW in the diffusion/dendritic growth layer, suggesting a favorable scenario for cloud seeding. Significance Statement: Understanding snowfall amounts and distribution in the mountains and how it is linked to topography, synoptic flow, and microphysical processes will help in the development of effective strategies for cloud seeding operations, managing runoff, reservoir, and mitigating flood risks, garnering substantial interest from stakeholders and the government agencies. [ABSTRACT FROM AUTHOR]

Rain formation is a critical factor governing the lifecycle and radiative forcing of clouds and therefore it is a key element of weather and climate. Cloud microphysics-turbulence interactions occur across a wide range of scales and are challenging to represent in atmospheric models with limited resolution. Based on past experiments and idealized numerical simulations, it has been postulated that cloud turbulence accelerates rain formation by enhancing drop collision-coalescence. We provide substantial evidence for significant impacts of turbulence on the evolution of cloud droplet size distributions and rain formation by comparing high-resolution observations of cumulus congestus clouds with state-of-the-art large-eddy simulations coupled with a Lagrangian particlebased microphysics scheme. Turbulent coalescence must be included in the model to accurately represent the observed drop size distributions, especially for drizzle drop sizes at lower heights in the cloud. Turbulence causes earlier rain formation and greater rain accumulation compared to simulations with gravitational coalescence only. The observed rain size distribution tail just above cloud base follows a power law scaling that deviates from theoretical scalings considering either a purely gravitation collision kernel or a turbulent kernel neglecting droplet inertial effects, providing additional evidence for turbulent coalescence in clouds. In contrast, large aerosols acting as cloud condensation nuclei ("giant CCN") do not significantly impact rain formation owing to their long timescale to reach equilibrium wet size relative to the lifetime of rising cumulus thermals. Overall, turbulent drop coalescence exerts a dominant influence on rain initiation in warm cumulus clouds, with limited impacts of giant CCN. [ABSTRACT FROM AUTHOR]

Numerical simulations of turbulent moist Rayleigh–Bénard convection driving CCN activation and droplet growth in the laboratory Pi chamber are discussed. Supersaturation fluctuations come from isobaric mixing of warm and humid air rising from the lower boundary with colder air featuring lower water vapor concentrations descending from the upper boundary. Lagrangian particle–based microphysics is used to represent the growth of haze CCN and cloud droplets with kinetic, solute, and surface tension effects included. Dry CCN spectra in the range between 2- and 200-nm radii from field observations are considered. Increasing the total CCN concentration from pristine to polluted conditions results in an increase in the droplet concentration and reduction in the mean droplet radius and spectral width. These are in agreement with Pi chamber observations and numerical simulations, as well as with numerous past studies of CCN cloud-base activation in natural clouds. The key result is that a relatively small fraction of the available CCN is activated in the Pi chamber fluctuating supersaturations, from about a half in the pristine case to only a 10th in the polluted case. The activation fraction as a function of the dry CCN radius is similar in all simulations, close to zero at the CCN small end, increasing to a maximum at CCN radius around 50 nm, and decreasing to close to zero at the large CCN end. This is explained as too small supersaturations to activate small CCN as in natural clouds and insufficient time to allow large CCN reaching the critical radius. Significance Statement: Impact of turbulence on the formation and growth of cloud droplets is an important cloud physics problem. Laboratory experiments in the Michigan Technological University cloud chamber provide key insights into this problem. Numerical simulations of cloud chamber processes discussed in this paper complement laboratory experiments by providing insights difficult or impossible to obtain in the laboratory. The study contrasts the formation and growth of cloud droplets in the laboratory cloud chamber with processes taking place in natural clouds. The differences documented in this paper pose questions concerning the impact of turbulence on the formation and growth of cloud droplets as a result of interactions of clouds with their environment. [ABSTRACT FROM AUTHOR]

The dual-polarization radar characteristics of severe storms are commonly used as indicators to estimate the size and intensity of deep convective updrafts. In this study, we track rapid fluctuations in updraft intensity and size by objectively identifying polarimetric fingerprints such as ZDR and KDP columns, which serve as proxies for mixed-phase updraft strength. We quantify the volume of ZDR and KDP columns to evaluate their utility in diagnosing temporal variability in lightning flash characteristics. Specifically, we analyze three severe storms that developed in environments with low-to-moderate instability and strong 0–6-km wind shear in northern Alabama during the 2016–17 VORTEX-Southeast field campaign. In these three cases (a tornadic supercell embedded in stratiform precipitation, a nontornadic supercell, and a supercell embedded within a quasi-linear convective system), we find that the volume of the KDP columns exhibits a stronger correlation with the total flash rate. The higher covariability of the KDP column volume with the total flash rate suggests that the overall electrification and precipitation microphysics were dominated by cold cloud processes. The lower covariability with the ZDR column volume indicates the presence of nonsteady updrafts or a less prominent role of warm rain processes in graupel growth and subsequent electrification. Furthermore, we observe that the majority of cloud-to-ground (CG) lightning strikes a carried negative charge to the ground. In contrast to findings from a tornadic supercell over the Great Plains, lightning flash initiations in the Alabama storms primarily occurred outside the footprint of the ZDR and KDP column objects. Significance Statement: This study quantifies the correlation between mixed-phase updraft intensity and total lightning flash rate in three severe storms in northern Alabama. In the absence of direct updraft velocity measurements, we use polarimetric signatures, such as ZDR and KDP columns, as proxies for updraft strength. Our analysis of polarimetric radar and lightning mapping array data reveals that the lightning flash rate is more highly correlated with the KDP column volume than with the ZDR column volume in all three storms examined. This contrasts with previous findings in storms over the central Great Plains, where the ZDR column volume showed higher covariability with flash rate. Interestingly, lightning initiation in the Alabama storms mainly occurred outside the ZDR and KDP column areas, contrary to previous findings. [ABSTRACT FROM AUTHOR]

Machine learning (ML) has rapidly transitioned from a niche activity to a mainstream tool for environmental research applications including atmospheric science cloud microphysics studies. Two recently developed cloud particle probes measure the light scattered in the near forward direction and save digital images of the scattering light. Scattering pattern images collected by the Particle Phase Discriminator mark 2, Karlsruhe edition (PPD-2K), and the Small Ice Detector, version 3 (SID-3) provide valuable information for particle shape and size characterization. Since different particle shapes have distinctly different light scattering characteristics, the images are ideally suited for ML. Here, results of a ML project to characterize ice particle shapes sampled by the PPD-2K in ice fog and diamond dust during a 3-yr project in Fairbanks, Alaska. About 2.15 million light-scattering pattern images were collected during 3 years of measurements with the PPD-2K. Visual Geometry Group (VGG) convolutional neural network (CNN) was trained to categorize light-scattering patterns into eight categories. Initial training images (120 each category) were selected by human visual examination of data, and the training dataset was augmented using an automated iterative method for image identification of further images which were all visually inspected by a human. Results were well correlated to similar categories identified from previously developed classification algorithms. ML identifies characteristics not included in automated analysis such as sublimation. Of the 2.15 million images analyzed, 1.3% were categorized as spherical (liquid), 43.5% were categorized as having rough surfaces, 15.3% were pristine, 16.3% were categorized as sublimating, and the remaining 23.6% did not fit into any of those categories (irregular or saturated). Significance Statement: The shapes and sizes of cloud particles can be extremely important for understanding the conditions that exist in the cloud. In this study, we show that more information about cloud particle characteristics can be identified by using machine learning than by traditional means. To demonstrate this, data are analyzed from a 3-yr study of ice fog and diamond dust events in Fairbanks, Alaska. The Particle Phase Discriminator instrument collected 2.15 million light-scattering pattern images of cloud particles during ground-based measurements. Neither traditional techniques nor machine learning were able to identify all categories, but a combination of both techniques led to a more complete view of ice particle shapes. [ABSTRACT FROM AUTHOR]

The article provides an overview of the use of Global Navigation Satellite System (GNSS) polarimetric radio occultation (PRO) measurements for weather prediction and applications. It summarizes a workshop that brought together participants from various countries to discuss the status and potential of PRO. The article highlights the benefits of PRO in improving weather forecasts and validating cloud microphysical processes. It also discusses ongoing experiments and the potential use of PRO observations in weather forecasting, climate modeling, and evaluating satellite precipitation products. The article also mentions the consideration of future constellations of small satellites for PRO measurements. [Extracted from the article]

The growth and breakup processes of raindrops within a cloud influence the rain intensity and the sizes of raindrops on the surface. The Doppler velocity spectrum acquired by a vertically pointing radar (VPR) contains information on atmospheric turbulence and the size classification of falling hydrometeors. In this study, the four types of Convective Cells (CC) during precipitation events with more than 700 mm of precipitation in southern China are described. The characteristics of four types of CCs correspond to the isolated convection, the early stage, the mature stage, and the decline stage of organizational convection, in that order. Microphysical analysis using retrieval of vertical air motion (Vair) and raindrop evolution in clouds from Doppler velocity spectra collected by C-band VPR revealed the growth and breakup of falling raindrops with dynamic impact. Larger raindrops appear in the early stages and are accompanied by ice particles, which are impacted by the falling path's downdraft. Raindrop aggregation, which is primarily related to the alternation of updraft and downdraft, accounts for the mature stage's high efficiency of surface rainfall. The CCs in the decline stage originate from the shallow uplift in the weak and broad downdraft under conditions of enough water vapor. The updraft dominates the stage of isolated convection. Observations of convective cells could be more accurately represented in model evaluations. [ABSTRACT FROM AUTHOR]

Understanding how entrainment and mixing shape the cloud droplet size distribution (DSD) is crucial for understanding the optical properties and precipitation efficiency of clouds. Different mixing scenarios, mainly homogeneous and inhomogeneous, shape the DSD in a distinct way and alter the cloud's impact on climate. However, the prevalence of these mixing scenarios and how they vary in space and time is still uncertain, as underlying processes are commonly unresolved by conventional numerical models. To overcome this challenge, we employ the L3 model, which considers supersaturation fluctuations and turbulent mixing down to the finest relevant lengthscales, making it possible to represent different mixing scenarios realistically. We investigate the spatial and temporal evolution of mixing scenarios over the life cycle of shallow cumulus clouds for varying boundary layer humidities and aerosol concentrations. Our findings suggest homogeneous mixing is generally predominant in cumulus clouds, while different mixing scenarios occur concurrently in the same cloud. Notably, inhomogeneous mixing increases over the cloud life cycle across all analyzed cases. The mean and standard deviation of supersaturation are found to be the most capable indicators of this evolution, providing a comprehensive insight into the characteristics of mixing scenarios. Finally, we show inhomogeneous mixing is more prevalent in drier boundary layers and for higher aerosol concentrations, underscoring the need for a more comprehensive investigation of how these mixing dynamics evolve in a changing climate. Plain Language Summary: Clouds play a crucial role in Earth's climate system by influencing the radiation balance and moisture transfer. When clouds mix with their environment, which occurs during a process known as entrainment and mixing, the number and size of cloud droplets change, affecting cloud optical properties and their ability to precipitate. This mixing can happen in two major ways: either the cloud and environmental air are rapidly mixed (homogeneous mixing) or slowly mixed (inhomogeneous mixing). However, accurately representing these effects, particularly inhomogeneous mixing, is challenging in nearly all atmospheric models, including weather or climate models, as essential scales and processes are unresolved. By applying a high‐resolution model that resolves all essential scales and processes, we investigate how mixing scenarios evolve over the cloud life cycle. We show that while homogeneous mixing prevails, older (mature or dissipating) clouds tend to mix more inhomogeneously. Moreover, clouds in drier and more polluted environments are shown to favor inhomogeneous mixing. In these environments, inhomogeneous mixing becomes the predominant mixing scenario for older clouds, emphasizing the need for further investigations on the impact of different mixing scenarios in a changing climate. Key Points: The evolution of predominant mixing scenarios during the shallow cumulus life cycle is investigatedWhile homogeneous mixing is prevalent, inhomogeneous mixing gradually increases as the cloud agesInhomogeneous mixing is more favored for higher aerosol concentrations and in drier boundary layers [ABSTRACT FROM AUTHOR]

The determination of whether a cloud will evolve into a thunderstorm is beneficial for understanding thunderstorm formation and important for ensuring the safety of society. However, a clear understanding of the microphysics in clouds for the occurrence of lightning activity has not been attained. Vast field observations and laboratory experiments indicate that graupel, which is rimed ice, is a vital hydrometeor for lightning generation, and is the foundation of riming electrification. In this study, polarimetric radar and lightning observations are used to compare the ice microphysics associated with graupel between 57 isolated thunderstorms and 39 isolated non-thunderstorms, and the differences in radar parameters are quantified. Our results for the occurrence of lightning activity in clouds showed the following results: 1) the maximum difference in graupel volume on the −10 °C isotherm height between thunderstorms and non-thunderstorms reached approximately 7.6 km³; 2) the graupel particles approached spherical shapes with a mean ZDR value of 0.3 dB, which likely indicated heavily rimed graupel was present; and 3) 98.2 % of thunderstorms were equipped with the ZDR column, and the mean depth was ~2.5 km³. Our study deepens our understanding of lighting physics and thunderstorm formation. [ABSTRACT FROM AUTHOR]

This paper describes the rainfall and microphysical structure of precipitation associated with Tropical Cyclone Ockhi using polarimetric Doppler weather radar (PDWR) products. The study reports the statistical analysis of precipitation types of tropical cyclone cloud systems over the north Indian Ocean, by combining the observations of the PDWR and disdrometer for the first time. There are few studies that mention initial DWR observations of TC over this low-latitude region below 23.5°N. We tried to further carry out a statistical analysis of precipitating clouds in a cyclone from its depression stage to severe cyclonic stage. This study tries to classify and quantify the contribution of convective and stratiform rain to the total TC rainfall. Precipitating clouds have been classified into convective and stratiform based on reflectivity measurements. The vertical profiles (VPRs) of the radar reflectivity Z and the differential reflectivity ZDR obtained for the stratiform and the convective events are compared. A study of the VPR of the convective events reveals that the Z and ZDR parameters tend to increase as the raindrops descend toward the ground owing to enhanced collision–coalescence processes. The VPR of stratiform rain shows signatures of the bright band (BB). The drop size distribution (DSD) parameters and rainfall rate pertaining to the two different precipitation regimes are estimated from the radar data and have been compared. The Joss–Waldvogel Disdrometer measurements have been used to derive DSD parameters and polarimetric rain-rate relation. [ABSTRACT FROM AUTHOR]

Previous work found that cold pools in ordinary convection are more sensitive to the microphysics scheme when the lifting condensation level (LCL) is higher owing to a greater evaporation potential, which magnifies microphysical uncertainties. In the current study, we explore whether the same reasoning can be applied to supercellular cold pools. To do this, four perturbed-microphysics ensembles are run, with each using an environment with a different LCL. Similar to ordinary convection, the sensitivity of supercellular cold pools to the microphysics increases with higher LCLs, though the physical reasoning for this increase in sensitivity differs from a previous study. Using buoyancy budgets along parcel trajectories that terminate in the cold pool, we find that negative buoyancy generated by microphysical cooling is partially countered by a decrease in environmental potential temperatures as the parcel descends. This partial erosion of negative buoyancy as parcels descend is most pronounced in the low-LCL storms, which have steeper vertical profiles of environmental potential temperature in the lower atmosphere. When this erosion is accounted for, the strength of the strongest cold pools in the low-LCL ensemble is reduced, resulting in a narrower distribution of cold pool strengths. This narrower distribution is indicative of reduced sensitivity to the microphysics. These results suggest that supercell behavior and supercell hazards (e.g., tornadoes) may be more predictable in low-LCL environments. Significance Statement: Thunderstorms typically produce "pools" of cold air beneath them owing in part to the evaporation of rain and melting of ice produced by the storm. Past work has found that in computer simulations of thunderstorms, the cold pools that form beneath thunderstorms are sensitive to how rain and ice are modeled in the simulation. In this study, we show that in the strongest thunderstorms that are capable of producing tornadoes, this sensitivity is reduced when the humidity in the lowest few kilometers above the surface is increased. Exploring why the sensitivity is reduced when the humidity increases provides a deeper understanding of the relationship between humidity and cold pool strength, which is important for severe storm forecasting. [ABSTRACT FROM AUTHOR]

As an essential moisture transport channel over the Tibetan Plateau (TP), accurate simulation of precipitation over the southeastern TP (SETP) is crucial for downstream moisture transport and water resource assessment. To improve SETP's precipitation simulation accuracy, a Turbulent Orographic Form Drag (TOFD) scheme, which represents the TOFD effect of subgrid orography on model grids, is employed. The impact of TOFD on cloud and precipitation during the summer monsoon over the SETP is investigated from dynamical and cloud microphysical perspectives. Through comparison with ground observations and satellite precipitation data, TOFD reduces the simulated precipitation bias. Despite reducing moisture transport, TOFD enhances low-level convergence and high-level divergence by influencing the low-level wind, thereby strengthening updraft and cloud and precipitation processes before the slope of SETP. TOFD also increases total precipitation before the slope by promoting heavy precipitation and reducing weak precipitation. However, TOFD diminishes the updraft and cloud and precipitation processes over the slope of SETP by influencing low-level wind. The reduced moisture transport by TOFD leads to decreased precipitation over the slope. Additionally, TOFD reduces total precipitation over the slope by promoting weak precipitation and reducing heavy precipitation. Furthermore, the impact mechanism of TOFD on precipitation depends on the wind speed, which is weakened during weak precipitation period with reduced wind speed and enhanced during heavy precipitation period with increased wind speed. These findings enhance the understanding of the complex terrain's influence on cloud and precipitation and provide a basis for improving cloud and precipitation simulation over the SETP. [ABSTRACT FROM AUTHOR]

Cloud seeding is considered a practical but unproved method to enhance precipitation or suppress hail, due to insufficient knowledge of ice formation and evolution after seeding clouds with ice nucleating particles. This study investigates the size effects on the immersion freezing of aerosol produced from commercial silver iodide (AgI) containing flares at mixed‐phase cloud temperatures from 243 to 267 K. Flare‐generated aerosol exhibited comparable ice nucleation ability (INA) to pure AgI particles in the size range of 200 and 400 nm. Non‐AgI impurities reduced the INA of flare‐generated particles ≤90 nm, which is lower than pure AgI particles ≤40 nm. The critical mass ice‐active site density of the generated aerosols (critical‐nm) was derived, indicating the minimum mass of AgI particles required for efficient ice nucleation. The new parameterization to predict critical‐nm can serve as a reference to optimize the effectiveness of cloud‐seeding materials for practical use. Plain Language Summary: Ice‐forming aerosol is commonly added to clouds, expecting precipitation enhancement via promotion of ice production. In this work, silver iodide (AgI) containing aerosol was generated from commercial cloud‐seeding products under different wind speed conditions. Its ice‐forming ability was studied at mixed‐phase cloud temperatures. The lower size limit for effective ice‐forming ability of the cloud‐seeding particles (90 nm) is higher than that of pure AgI particles (40 nm). The non‐AgI components produced by cloud‐seeding products are hypothesized to decrease the ice‐forming ability of smaller particles, as the mass fraction of ice‐nucleating AgI decreases. To estimate the minimum mass of AgI in a particle required for efficient ice nucleation under cloud‐seeding relevant conditions, we derived the critical ice‐activated mass fraction of the generated aerosols. These findings provide valuable insights into the optimization of cloud‐seeding practices for enhanced precipitation. Key Points: Silver iodide (AgI) containing cloud‐seeding aerosols exhibit comparable ice‐forming abilities to pure AgI at sizes of 200 and 400 nmNon‐AgI impurities produced from flare burning decrease the ice nucleation ability of particles smaller than 90 nmA new parameterization is presented to estimate the minimum mass of AgI particles required to maximize glaciogenic seeding [ABSTRACT FROM AUTHOR]

A heavy orographic precipitation event associated with the postfrontal sector of a midlatitude cyclone over the Australian Snowy Mountains (ASM) was analyzed using field observations and numerical simulations. This event, observed during a 2018 intensive field campaign, was of particular interest as three distinct precipitation episodes were identified within a prolonged postfrontal period. Deep mixed‐phase clouds (MPCs) characterized by cold cloud‐top temperatures (colder than −30°C) and the presence of updrafts extending 3.5–4.5 km above the boundary‐layer height, produced the three enhanced precipitation events over the windward slopes of the ASM. The presence of conditional instabilities and deep updrafts were also found in the sounding and Doppler velocity observations respectively, while the cloud radar observations show the deep MPCs with cloud tops reaching to 6–7 km a.s.l. Orographic convection invigoration was found to be the main mechanism producing the precipitation enhancement over the windward slopes and higher terrain. Using the Weather Research/Forecasting model, we analyzed the rates of microphysical processes to explicitly account for the enhancement of precipitation formation processes in these MPCs. This analysis showed that the precipitation formation processes were further enhanced through depositional and riming growth of ice‐phase hydrometeors during the three precipitation events. Deposition is simulated at higher levels (above the −15°C level) and most likely enabled by deep convective updrafts through the midtroposphere, whereas riming is stronger at lower levels (below −10°C level) due to the persistent production of feeder supercooled liquid water clouds sustained by the orographic lifting. [ABSTRACT FROM AUTHOR]

Abstract Different turbulent entrainment‐mixing mechanisms between clouds and environment are essential to cloud‐related processes; however, accurate representation of entrainment‐mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment‐mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment‐mixing. Specifically, number concentration after entrainment‐mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment‐mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment‐mixing is negatively correlated with liquid water content affected by entrainment‐mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment‐mixing in weather/climate models.

Abstract Numerous Earth system models exhibit excessive aerosol effective forcing at the top of the atmosphere (TOA), including the Department of Energy's Energy Exascale Earth System Model (E3SM). Here, in the context of the E3SM version 3 effort, the predicted particle property (P3) stratiform cloud microphysics scheme and an enhanced deep convection parameterization suite (ZM_plus) are implemented into E3SM. The ZM_plus includes a convective cloud microphysics scheme, a multi‐scale coherent structure parameterization for mesoscale convective systems, and a revised cloud base mass flux formulation considering impacts of the large‐scale environment. The P3 scheme improved cloud and radiation particularly over the Northern Hemisphere and the frequency of heavy precipitation over the tropics, and the ZM_plus improved clouds in the tropics. P3 decreases aerosol effective forcing by 0.15 W m−2, while the ZM_plus increases it by 0.27 W m−2, resulting from excessive direct (0.31 W m−2) and indirect forcing (−1.79 W m−2). The excessive aerosol forcings are due to aerosol overestimation associated with insufficient aerosol wet removal. By improving the physical treatments in the aerosol wet removal, we effectively mitigate anthropogenic aerosol overestimation and thus attenuate direct (0.09 W m−2) and indirect aerosol forcing (−1.52 W m−2). Adjustment to primary organic matter hygroscopicity reduces direct and indirect forcing to more reasonable values: −0.13 W m−2 and −1.31 W m−2, respectively. On climatology, improved aerosol treatments mitigate overestimation of aerosol optical depth.

Bin and bulk schemes are the two primary methods to parameterize cloud microphysical processes. This study attempts to reveal how their structural differences (size-resolved vs. moment-resolved) manifest in terms of cloud and precipitation properties. We use a bulk scheme, the Arbitrary Moment Predictor (AMP), which uses process parameterizations identical to those in a bin scheme but predicts only moments of the size distribution like a bulk scheme. As such, differences between simulations using AMP's bin scheme and simulations using AMP itself must come from their structural differences. In one-dimensional kinematic simulations, the overall difference between AMP (bulk) and bin schemes is found to be small. Full-microphysics AMP and bin simulations have similar mean liquid water path (mean percent difference

Anemophilous (wind‐driven) pollen is one type of primary biological aerosol particle, which can rupture under high humidity conditions and form smaller sub‐pollen particles (SPPs). Both pollen and SPPs can reach the upper troposphere under convective conditions, acting as cloud condensation nuclei (CCN) and ice nucleating particles (INPs), thus influencing cloud formation and precipitation. However, the impacts of these biological aerosols on cold cloud formation and local climate remain unclear as there are large uncertainties on their emission flux and ice nucleating abilities. Here, we incorporate pollen emission and rupture processes in the Weather Research and Forecasting Model with Chemistry (WRF‐Chem) simulations and update the Morrison microphysics scheme within WRF‐Chem using aerosol‐aware INP parameterizations to account for pollen in addition to other anthropogenic and biogenic aerosol. INP parameterizations for pollen and SPP are derived from laboratory experiments. When including pollen rupture rates as observed in a series of chamber studies, SPP concentrations increase, leading to an increase of cloud ice and water by up to 50% and potentially extending the duration of the convective system. Among all simulated hydrometeors, graupel and raindrops exhibit the largest enhancements from the inclusion of SPPs, with intensifying precipitation at the backside of the convective system and a greater spatial extent. Sensitivity simulations indicate that SPPs have a greater effect on cloud microphysical processes than whole pollen grains, and further observational evidence is needed to constrain these processes. Plain Language Summary: Pollen grains are emitted by vegetation and carry the reproductive material of plants. Under certain conditions (e.g., high humidity, lightning), pollen grains can rupture to form smaller sub‐pollen particles (SPPs). Both intact pollen grain and SPPs have been observed to influence the formation of cloud liquid droplets and ice crystals, potentially impacting cloud properties and precipitation. In this study, we include pollen emission and pollen‐cloud interactions into a meteorology‐chemistry model and investigate the effects of pollen grains and SPPs on a strong thunderstorm event associated with a mesoscale convective system. We find that SPPs have a greater impact than whole pollen grains on the convective system properties by increasing cloud ice and water content, facilitating the transition of cloud droplets and ice into graupel, and ultimately intensifying precipitation. These changes alter the impact of the convective system, leading to an extension of the storm width and lengthening of its duration. Key Points: Sub‐pollen particles (SPPs) enhance cloud microphysical processes more than intact pollen grainsSPPs influence the spatial extent and vertical structure of convective systemsSimulated SPP impacts are intensified when incorporating laboratory measured pollen rupture rates [ABSTRACT FROM AUTHOR]

Lightning stroke strength, characterized by energy and peak currents, over the Tibetan Plateau (TP), is investigated by utilizing datasets from the World Wide Lightning Location Network and the Chinese Cloud-to-Ground Lightning Location System during 2016–2019. Focused on the south-central (SC) and southeast (SE) of the TP, it reveals that SE-TP experiences strokes with larger average energy and peak currents. Strong strokes (energy ≥ 100 kJ or peak currents ≥ |100| kA), exhibiting bimodal distribution in winter and summer, are more frequent and have larger average values over the SE-TP than the SC-TP, with diurnal distribution indicating peaks in energy and positive strokes in the middle of the night and negative strokes peaking in the morning. Utilizing the ECMWF/ERA-5 and MERRA-2 reanalysis, we find that stronger strokes correlate with thinner charge zone depths and larger CIWCFs but stable warm cloud depths and zero-degree levels over the SC-TP. Over the SE-TP, stronger strokes are associated with smaller CIWCFs and show turning points for warm cloud depths and zero-degree levels. Thicker charge zone depths correlate with stronger negative strokes but weaker positive strokes. Generating strokes of similar strength over the SC-TP requires larger CIWCFs, thinner warm cloud depths, and lower zero-degree levels than over the SE-TP. [ABSTRACT FROM AUTHOR]

The Sichuan-Xizang Railway is an important part of the railway network in China, and geological disasters, such as mountain floods and landslides, frequently occur in this region. Precipitation is an important cause of these disasters; therefore, accurate simulation of the precipitation in this region is highly important. In this study, the descriptions for uncertain processes in the cloud microphysics scheme are improved; these processes include cloud droplet activation, cloud-rain autoconversion, rain accretion by cloud droplets, and the entrainment-mixing process. In the default scheme, the cloud water content of different sizes corresponds to the same cloud droplet concentration, which is inconsistent with the actual content; this results in excessive cloud droplet size, unreasonable related conversion rates of microphysical process (such as cloud-rain autoconversion), and an overestimation of precipitation. Our new scheme overcomes the problem of excessive cloud droplet size. The processes of cloudrain autoconversion and rain accretion by cloud droplets are similar to the stochastic collection equation, and the mixing mechanism of cloud droplets is more consistent with that occurred during the actual physical process in the cloud. Based on the new and old schemes, multiple precipitation processes in the flood season of 2021 along the Sichuan-Xizang Railway are simulated, and the results are evaluated using ground observations and satellite data. Compared to the default scheme, the new scheme is more suitable for the simulation of cloud physics, reducing the simulation deviation of the liquid water path and droplet radius from 2 times to less than 1 time and significantly alleviating the overestimation of precipitation intensity and range of precipitation center. The average root-mean-square error is reduced by 22%. Our results can provide a scientific reference for improving precipitation forecasting and disaster prevention in this region. [ABSTRACT FROM AUTHOR]

Raindrop formation processes in warm clouds mainly consist of condensation and collision–coalescence of small cloud droplets. Once raindrops form, they can continue growing through collection of cloud droplets and self-collection. In this study, we develop novel emulators to represent raindrop formation as a function of various physical or background environmental conditions by using a sophisticated aerosol–cloud model containing 300 droplet size bins and machine learning methods. The emulators are then implemented in two microphysics schemes in the Weather Research and Forecasting Model and tested in two idealized cases. The simulations of shallow convection with the emulators show a clear enhancement of raindrop formation compared to the original simulations, regardless of the scheme in which they were embedded. On the other hand, the simulations of deep convection show a more complex response to the implementation of the emulators, in terms of the changes in the amount of rainfall, due to the larger number of microphysical processes involved in the cloud system (i.e., ice-phase processes). Our results suggest the potential of emulators to replace the conventional parameterizations, which may allow us to improve the representation of physical processes at an affordable computational expense. Significance Statement: Formation of raindrops marks a critical stage in cloud evolution. Accurate representations of raindrop formation processes require detailed calculations of cloud droplet growth processes. These calculations are often not affordable in weather and climate models as they are computationally expensive due to their complex dependence on cloud droplet size distributions and dynamical conditions. As a result, simplified parameterizations are more frequently used. In our study we trained machine learning models to learn raindrop formation rates from detailed calculations of cloud droplet evolutions in 1000 parcel-model simulations. The implementation of the developed models or the emulators in a weather forecasting model shows a change in the total rainfall and cloud characteristics, indicating the potential improvement of cloud representations in models if these emulators replace the conventional parameterizations. [ABSTRACT FROM AUTHOR]

NASA's Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) field campaign gathered data using "satellite-simulating" (albeit with higher-resolution data than satellites currently provide) and in situ aircraft to study snowstorms, with an emphasis on banding. This study used three IMPACTS microwave instruments—two passive and one active—chosen for their sensitivity to precipitation microphysics. The 10–37-GHz passive frequencies were well suited for detecting light precipitation and differentiating rain intensities over water. The 85–183-GHz frequencies were more sensitive to cloud ice, with higher cloud tops manifesting as lower brightness temperatures, but this did not necessarily correspond well to near-surface precipitation. Over land, retrieving precipitation information from radiometer data is more difficult, requiring increased reliance on radar to assess storm structure. A dual-frequency ratio (DFR) derived from the radar's Ku- and Ka-band frequencies provided greater insight into storm microphysics than reflectivity alone. Areas likely to contain mixed-phase precipitation (often the melting layer/bright band) generally had the highest DFR, and high-altitude regions likely to contain ice usually had the lowest DFR. The DFR of rain columns increased toward the ground, and snowbands appeared as high-DFR anomalies. Significance Statement: Winter precipitation was studied using three airborne microwave sensors. Two were passive radiometers covering a broad range of frequencies, while the other was a two-frequency radar. The radiometers did a good job of characterizing the horizontal structure of winter storms when they were over water, but struggled to provide detailed information about winter storms when they were over land. The radar was able to provide vertically resolved details of storm structure over land or water, but only provided information at nadir, so horizontal structure was less well described. The combined use of all three instruments compensated for individual deficiencies, and was very effective at characterizing overall winter storm structure. [ABSTRACT FROM AUTHOR]

In the last decade, substantial improvements have been achieved in quantitative satellite precipitation estimates, which are essential for a wide range of applications. In this study, we evaluated the performance of Integrated Multi-satellitE Retrievals for GPM (IMERG V06B) at the sub-daily and daily scales. Ten years of half-hourly precipitation records aggregated at different sub-daily periods were evaluated over a region in the Western Mediterranean. The analysis at the half-hourly scale examined the contribution of passive microwave (PMW) and infrared (IR) sources in IMERG estimates, as well as the relationship between various microphysical cloud properties using Cloud Microphysics (CMIC–NWC SAF) data. The results show the following: (1) a marked tendency to underestimate precipitation compared to rain gauges which increases with rainfall intensity and temporal resolution, (2) a weaker negative bias for retrievals with PMW data, (3) an increased bias when filling PMW gaps by including IR information, and (4) an improved performance in the presence of precipitating ice clouds compared to warm and mixed-phase clouds. This work contributes to the understanding of the factors affecting satellite estimates of extreme precipitation. Their relationship with the microphysical characteristics of clouds generates added value for further downstream applications and users' decision making. [ABSTRACT FROM AUTHOR]

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column. Significance Statement: Secondary ice production (SIP) remains a key challenge in our understanding of boundary layer mixed-phase clouds. Here, we use sensitivity experiments performed with the ICON model at the cloud-resolving scale to explore potential interactions between primary nucleation, SIP, and mesoscale cloud organization. We simulate an Arctic single-layer cold air outbreak stratocumulus deck that was sampled during the M-PACE campaign. We find that once established, SIP alone is sufficient to maintain the mixed-phase cloud state until the end of the simulation. Our sensitivity analysis also shows that numerically more efficient treatments of immersion freezing are statistically invariant from simulations with a full prognostic INP budget. [ABSTRACT FROM AUTHOR]

Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow ("ice crystal process"), forming "cold" precipitation. Equally, cloud droplets, when large enough, coalesce to form "warm" precipitation by the "warm rain process." Warm rain could be supercooled and freeze as "warm" graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, "cold ice" (crystals, snow), and "warm rain" (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system's precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement: The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and "wins" the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space. [ABSTRACT FROM AUTHOR]

Copyright of Plateau Meteorology is the property of Plateau Meteorology Editorial Office and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Supercooled liquid clouds are ubiquitous over the Southern Ocean (SO), even to temperatures below −20°C, and comprise a large fraction of the marine boundary layer (MBL) clouds. Earth system models and reanalysis products have struggled to reproduce the observed cloud phase distribution and occurrence of cloud ice in the region. Recent simulations found the microphysical representation of ice nucleation and growth has a large impact on these properties, however, measurements of SO ice nucleating particles (INPs) to validate simulations are sparse. This study presents measurements of INPs from simultaneous aircraft and ship campaigns conducted over the SO in austral summer 2018, which include the first in situ observations in and above cloud in the region. Our results confirm recent observations that INP concentrations are uniformly lower than measurements made in the late 1960s. While INP concentrations below and above cloud are similar, higher ice nucleation efficiency above cloud supports model simulations that the dominant INP composition varies with height. Model parameterizations based solely on aerosol properties capture the mean relationship between INP concentration and temperature but not the observed variability, which is likely related to the only modest correlations observed between INPs and environmental or aerosol metrics. Including wind speed in addition to activation temperature in a marine INP parameterization reduces bias but does not explain the large range of observed INP concentrations. Direct and indirect inference of marine INP size suggests MBL INPs, at least during Austral summer, are dominated by particles with diameters smaller than 500 nm. Plain Language Summary: Although Antarctica is remote, the continent and the Southern Ocean (SO) that surrounds it play a fundamental role in shaping regional and global climate. The clouds in this region are unique, with less ice and more liquid water present at low temperatures than in other areas. This is likely related to very low concentrations of rare aerosol particles called ice nucleating particles (INPs), which cause liquid water droplets in clouds to freeze. Largely due to a lack of observations, SO clouds are poorly represented in global models, and the interactions between aerosol particles and clouds are one of the largest remaining uncertainties. This study presents results of INP measurements from several recent field campaigns over the SO, including the first observations within and above clouds in the region. Our results suggest different types of particles are present below and above clouds, which have varying ability to nucleate ice. They also highlight the need for additional measurements of INP composition and size, which are key variables needed to improve model simulations. Key Points: First vertically resolved measurements of ice nucleating particles (INPs) over the Southern Ocean, including in‐cloud observationsCorrelation between normalized INP concentrations and wind speed suggests marine active site density is variableHigher ice nucleation efficiency observed above cloud, consistent with an increasing influence of mineral dust with height [ABSTRACT FROM AUTHOR]

Recent studies from the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE) demonstrated definitive radar evidence of seeding signatures in winter orographic clouds during three intensive operation periods (IOPs) where the background signal from natural precipitation was weak and a radar signal attributable to seeding could be identified as traceable seeding lines. Except for the three IOPs where seeding was detected, background natural snowfall was present during seeding operations and no clear seeding signatures were detected. This paper provides a quantitative analysis to assess if orographic cloud seeding effects are detectable using radar when background precipitation is present. We show that a 5-dB change in equivalent reflectivity factor Ze is required to stand out against background natural Ze variability. This analysis considers four radar wavelengths, a range of background ice water contents (IWC) from 0.012 to 1.214 g m−3, and additional IWC introduced by seeding ranging from 0.012 to 0.486 g m−3. The upper-limit values of seeded IWC are based on measurements of IWC from the Nevzorov probe employed on the University of Wyoming King Air aircraft during SNOWIE. This analysis implies that seeding effects will be undetectable using radar within background snowfall unless the background IWC is small, and the seeding effects are large. It therefore remains uncertain whether seeding had no effect on cloud microstructure, and therefore produced no signature on radar, or whether seeding did have an effect, but that effect was undetectable against the background reflectivity associated with naturally produced precipitation. Significance Statement: Operational glaciogenic seeding programs targeting wintertime orographic clouds are funded by a range of stakeholders to increase snowpack. Glaciogenic seeding signatures have been observed by radar when natural background snowfall is weak but never when heavy background precipitation was present. This analysis quantitatively shows that seeding effects will be undetectable using radar reflectivity under conditions of background snowfall unless the background snowfall is weak, and the seeding effects are large. It therefore remains uncertain whether seeding had no effect on cloud microstructure, and therefore produced no signature on radar, or whether seeding did have an effect, but that effect was undetectable against the background reflectivity associated with naturally produced precipitation. Alternative assessment methods such as trace element analysis in snow, aircraft measurements, precipitation measurements, and modeling should be used to determine the efficacy of orographic cloud seeding when heavy background precipitation is present. [ABSTRACT FROM AUTHOR]

This study aims to characterize the shapes and fall speeds of ice pellets formed in various atmospheric conditions and to investigate the possibility to use a laser-optical disdrometer to distinguish between ice pellets and other types of precipitation. To do so, four ice pellet events were documented using manual observations, macrophotography, and laser-optical disdrometer data. First, various ice pellet fall speeds and shapes, including spherical, bulged, fractured, and irregular particles, were associated with distinct atmospheric conditions. A higher fraction of bulged and fractured ice pellets was observed when solid precipitation was completely melted aloft while more irregular particles were observed during partial melting. These characteristics affected the diameter–fall speed relations measured. Second, the measurements of particles' fall speed and diameter show that ice pellets could be differentiated from rain or freezing rain. Ice pellets larger than 1.5 mm tend to fall > 0.5 m s−1 slower than raindrops of the same size. In addition, the fall speed of a small fraction of ice pellets was < 2 m s−1 regardless of their size, as compared with a fall speed > 3 m s−1 for ice pellets with diameter > 1.5 mm. Video analysis suggests that these slower particles could be ice pellets passing through the laser-optical disdrometer after colliding with the head of the instrument. Overall, these findings contribute to a better understanding of the microphysics of ice pellets and their measurement using a laser-optical disdrometer. Significance Statement: Ice pellets are challenging to forecast and to detect automatically. In this study, we documented the fall speed and physical characteristics of ice pellets during various atmospheric conditions using a combination of a laser-optical disdrometer, manual observations, and macrophotography images. Relationships were found between the shape and fall speed of ice pellets. These findings could be used to refine the parameterization of ice pellets in atmospheric models and, consequently, improve the forecast of impactful winter precipitation types such as freezing rain. Furthermore, they will also help to physically interpret laser-optical disdrometer data during ice pellets and freezing rain. [ABSTRACT FROM AUTHOR]

High‐resolution measurements from the Airborne Cloud‐Turbulence Observation System (ACTOS) during the Azores Stratocumulus Measurements of Radiation, Turbulence and Aerosols (ACORES) campaign are analysed for an investigation of the vertical profiles of microphysical properties and entrainment velocity (We) in marine stratocumulus clouds. The vertical profiles show the transition from the cloudy layer to free troposphere with nearly linear profiles of total water mixing ratio, liquid water potential temperature and virtual potential temperature, but the thickness of entrainment interfacial layer varies significantly. Sharp transitions of cloud microphysical and optical properties within a single horizontal flight leg are found in one stratocumulus cloud system. They seem to be related to the local environmental conditions, such as the wind shear and turbulent dissipation rate. We values estimated by three methods show consistent tendencies in general and are clearly related to the local environmental conditions, such as vertical shear of the horizontal wind and turbulence intensity. However, the magnitudes of We values differ by up to two orders of magnitude depending on the methods, which suggests that the estimation of We from in situ measurements is still a challenge. Analysis of the microphysical response to entrainment suggests that inhomogeneous mixing occurs dominantly. On the other hand, the analysis results for the clouds under more humid conditions indicate a higher likelihood of homogeneous mixing. [ABSTRACT FROM AUTHOR]

An accurate forecast of Global Horizontal solar Irradiance (GHI) and Direct Normal Irradiance (DNI) in cloudy conditions remains a major challenge in the solar energy industry. This study focuses on the impact of cloud microphysics on GHI and its partition into DNI and Diffuse Horizontal Irradiance (DHI) using the Weather Research and Forecasting model specifically designed for solar radiation applications (WRF-Solar) and seven microphysical schemes. Three stratocumulus (Sc) and five shallow cumulus (Cu) cases are simulated and evaluated against measurements at the US Department of Energy's Atmospheric Radiation Measurement (ARM) user facility, Southern Great Plains (SGP) site. Results show that different microphysical schemes lead to spreads in simulated solar irradiance components up to 75% and 350% from their ensemble means in the Cu and Sc cases, respectively. The Cu cases have smaller microphysical sensitivity due to a limited cloud fraction and smaller domain-averaged cloud water mixing ratio compared to Sc cases. Cloud properties also influence the partition of GHI into DNI and DHI, and the model simulates better GHI than DNI and DHI due to a non-physical error compensation between DNI and DHI. The microphysical schemes that produce more accurate liquid water paths and effective radii of cloud droplets have a better overall performance. [ABSTRACT FROM AUTHOR]

The supersaturation equation for a vertically moving adiabatic cloud parcel is analyzed. The effects of turbulent updrafts are incorporated in the shape of a stochastic Lagrangian model, with spatial and time correlations expressed in terms of turbulent kinetic energy. Using the Fokker–Planck equation, the steady-state probability distributions of supersaturation are analytically computed for a number of approximations involving the time-scale separation between updraft fluctuations and phase relaxation, and droplet or ice particle size fluctuations. While the analytical results are presented in general for single-phase clouds, the calculated distributions are used to compute mixed-phase cloud properties—mixed fraction and mean liquid water content in an initially icy cloud—and are argued to be useful for generalizing and constructing new parameterization schemes. Significance Statement: Supersaturation is the fuel for the development of clouds in the atmosphere. In this paper, our goal is to better understand the supersaturation budget of clouds embedded in a turbulent environment by analyzing the basic equations of cloud microphysics. It is found that the turbulent characteristics of an air parcel substantially affect the cloud's supersaturation budget and hence its life cycle. This is also shown in the context of mixed-phase clouds where, depending on the turbulent regime, different liquid-to-ice ratios are found. Consequently, the theoretical approach of this paper is crucial to develop tools to parameterize small-scale atmospheric features, like clouds, into global circulation models to improve climate projections for the future. [ABSTRACT FROM AUTHOR]

Water vapor supersaturation in clouds is a random variable that drives activation and growth of cloud droplets. The Pi Convection–Cloud Chamber generates a turbulent cloud with a microphysical steady state that can be varied from clean to polluted by adjusting the aerosol injection rate. The supersaturation distribution and its moments, e.g., mean and variance, are investigated for varying cloud microphysical conditions. High-speed and collocated Eulerian measurements of temperature and water vapor concentration are combined to obtain the temporally resolved supersaturation distribution. This allows quantification of the contributions of variances and covariances between water vapor and temperature. Results are consistent with expectations for a convection chamber, with strong correlation between water vapor and temperature; departures from ideal behavior can be explained as resulting from dry regions on the warm boundary, analogous to entrainment. The saturation ratio distribution is measured under conditions that show monotonic increase of liquid water content and decrease of mean droplet diameter with increasing aerosol injection rate. The change in liquid water content is proportional to the change in water vapor concentration between no-cloud and cloudy conditions. Variability in the supersaturation remains even after cloud droplets are formed, and no significant buffering is observed. Results are interpreted in terms of a cloud microphysical Damköhler number (Da), under conditions corresponding to Da ≲ 1 , i.e., the slow-microphysics regime. This implies that clouds with very clean regions, such that Da ≲ 1 is satisfied, will experience supersaturation fluctuations without them being buffered by cloud droplet growth. Significance Statement: The saturation ratio (humidity) in clouds controls the growth rate and formation of cloud droplets. When air in a turbulent cloud mixes, the humidity varies in space and time throughout the cloud. This is important because it means cloud droplets experience different growth histories, thereby resulting in broader size distributions. It is often assumed that growth and evaporation of cloud droplets buffers out some of the humidity variations. Measuring these variations has been difficult, especially in the field. The purpose of this study is to measure the saturation ratio distribution in clouds with a range of conditions. We measure the in-cloud saturation ratio using a convection cloud chamber with clean to polluted cloud properties. We found in clouds with low concentrations of droplets that the variations in the saturation ratio are not suppressed. [ABSTRACT FROM AUTHOR]

In simulations of Arctic mixed-phase clouds, cloud persistence and the liquid water path (LWP) are sensitive to ice particle number concentrations. Here, we explore sensitivities of cloud microphysical properties to the dominant ice particle shape (dendrites, plates, columns, or spheres) using the SCALE-AMPS large-eddy simulation model. AMPS is a bin microphysics scheme that predicts particle shapes based on the inherent growth ratio (IGR) of spheroids, which determines vapor depositional growth rates along the a and c axes, and the rimed and aggregate mass fractions. We examine the impacts of various IGR values on simulations of clouds observed during the M-PACE and SHEBA experiments. Under M-PACE (SHEBA) conditions, LWP varies between 49 (1.1) and 230 (6.7) g m−2, and the ice water path (IWP) varies between 3 (0.03) and 40 (0.12) g m−2, depending on the ice shape. The lowest LWP and the highest IWP are obtained when columnar particles dominate because their low terminal velocities and large capacitance and collisional area result in large vapor deposition and riming rates, whereas the highest LWP and lowest IWP are obtained when spherical particles dominate because their vapor deposition and riming rates are low. Because ice particle shape significantly influences simulated Arctic mixed-phase clouds, reliable simulations require accurately estimated IGR values under various atmospheric conditions. Finally, comparisons between the simulation results and observations show that the size distribution larger than 2000 μm is better reproduced when the increase in rimed mass that causes ice particles to become spherical is suppressed. Significance Statement: Atmospheric models have difficulties in reproducing Arctic mixed-phase clouds because of uncertainties in the parameterization of microphysical processes. This is the first study to use a large-eddy simulation model implemented with a habit-predicting bin microphysics scheme to demonstrate the important role of ice particle shape on the microphysical properties of both heavy-riming and no-riming mixed-phase clouds. We found the vapor deposition and riming rates to be greatly influenced by ice particle shape. By comparing the ice particle size distribution, mass–diameter relationship, and area ratio between simulation results and observations, we show that a hexagonal ice shape model and a riming model that simply converts ice crystals to graupel may not accurately reproduce actual heavy-riming clouds. [ABSTRACT FROM AUTHOR]

This study uses a new, unique dataset created by combining multi-Doppler radar wind and reflectivity analysis, diabatic Lagrangian analysis (DLA) retrievals of temperature and water substance, and a complex hail trajectory model to create millions of numerically simulated hail trajectories in the Kingfisher, Oklahoma, supercell on 29 May 2012. The DLA output variables are used to obtain a realistic, 4D depiction of the storm's thermal and hydrometeor structure as required input to the detailed hail growth trajectory model. Hail embryos are initialized in the hail growth module every 3 min of the radar analysis period (2251–0000 UTC) to produce over 2.7 million hail trajectories. A spatial integration technique considering all trajectories is used to identify locations within the supercell where melted particles and subsevere and severe hailstones reside in their lowest and highest concentrations. It is found that hailstones are more likely to reside for longer periods closer to the downshear updraft within the midlevel mesocyclone in a region of decelerated midlevel mesocyclonic horizontal flow, termed the downshear deceleration zone (DDZ). Additionally, clusters of trajectories are analyzed using a trajectory clustering method. Trajectory clusters show there are many trajectory pathways that result in hailstones ≥ 4.5 cm, including trajectories that begin upshear of the updraft away from ideal growth conditions and trajectories that grow within the DDZ. There are also trajectory clusters with similar shapes that experience widely different environmental and hailstone characteristics along the trajectory. Significance Statement: The purpose of this study is to understand how hail grew in a thunderstorm that was observed by numerous instruments. The observations were input into a hail trajectory model to simulate hail growth. We found a part of the storm near the updraft where hailstones could remain aloft longer and therefore grow larger. Most modeled severe hailstones were found in the storm in this region. However, we also found that there are many different pathways hailstones can take to become large, although there are still some common characteristics among the pathways. [ABSTRACT FROM AUTHOR]

Current bulk microphysical parameterization schemes underpredict precipitation intensities and drop size distributions (DSDs) during warm rain periods, particularly upwind of coastal terrain. To help address this deficiency, this study introduces a set of modifications, called RCON, to the liquid-phase (warm rain) parameterization currently used in the Thompson–Eidhammer microphysical parameterization scheme. RCON introduces several model modifications, motivated by evaluating simulations from a bin scheme, which together result in more accurate precipitation simulations during periods of warm rain. Among the most significant changes are 1) the use of a wider cloud water DSD of lognormal shape instead of the gamma DSD used by the Thompson–Eidhammer parameterization and 2) enhancement of the cloud-to-rain autoconversion parameterization. Evaluation of RCON is performed for two warm rain events and an extended period during the Olympic Mountains Experiment (OLYMPEX) field campaign of winter 2015/16. We show that RCON modifications produce more realistic precipitation distributions and rain DSDs than the default Thompson–Eidhammer configuration. For the multimonth OLYMPEX period, we show that rain rates, rainwater mixing ratios, and raindrop number concentrations were increased relative to the Thompson–Eidhammer microphysical parameterization, while concurrently decreasing raindrop diameters in liquid-phase clouds. These changes are consistent with an increase in simulated warm rain. Finally, real-time evaluation of the scheme from August 2021 to August 2022 demonstrated improved precipitation prediction over coastal areas of the Pacific Northwest. Significance Statement: Although the accurate simulation of warm rain is critical to forecasting the hydrology of coastal areas and windward slopes, many warm rain parameterizations underpredict precipitation in these locations. This study introduces and evaluates modifications to the Thompson–Eidhammer microphysics parameterization scheme that significantly improve the accuracy of rainfall prediction in those regions. [ABSTRACT FROM AUTHOR]

Abstract Cloud seeding is considered a practical but unproved method to enhance precipitation or suppress hail, due to insufficient knowledge of ice formation and evolution after seeding clouds with ice nucleating particles. This study investigates the size effects on the immersion freezing of aerosol produced from commercial silver iodide (AgI) containing flares at mixed‐phase cloud temperatures from 243 to 267 K. Flare‐generated aerosol exhibited comparable ice nucleation ability (INA) to pure AgI particles in the size range of 200 and 400 nm. Non‐AgI impurities reduced the INA of flare‐generated particles ≤90 nm, which is lower than pure AgI particles ≤40 nm. The critical mass ice‐active site density of the generated aerosols (critical‐nm) was derived, indicating the minimum mass of AgI particles required for efficient ice nucleation. The new parameterization to predict critical‐nm can serve as a reference to optimize the effectiveness of cloud‐seeding materials for practical use.

In this study, we developed a global cloud system‐resolving model (GCSRM) incorporating stable water isotopes (NICAM‐WISO). Using a single‐moment cloud microphysics scheme, we applied the new model to conduct a current climate simulation at a horizontal resolution of 56 km. NICAM‐WISO simulated the seasonal means of the atmospheric hydrological cycle, as well as water isotopic ratios of precipitation and vapor. The model captured the general features of precipitation isotope effects and its spatial correlations were comparable to those of other isotope‐incorporated global atmospheric models. The model showed better spatial correlation between simulated and observed values for a fine‐horizontal‐resolution (14 km) simulation compared to coarse‐horizontal‐resolution (56 km) simulation. However, the model had isotopic biases in tropical mid‐troposphere, ocean, and cold continental regions. A comparison of stable water isotopes between the simulation and observations offered clues for improving the model. For example, in the tropical mid‐troposphere, we found a negative bias in the mixing ratio and isotopic ratio of water vapor. Our analysis using satellite retrievals revealed that these underestimations were caused by weak mixing with the boundary layer vapor and low raindrop evaporation with a small evaporation fraction. The underestimations indicated weak shallow convective mixing in the model, inducing negative bias in the mixing ratio and isotopic ratio of the mid‐tropospheric vapor. These biases were also seen in the fine horizontal‐resolution simulation. Furthermore, we conducted several km‐scale atmospheric isotope circulation simulations using NICAM‐WISO. We expect that global‐scale fine‐horizontal‐resolution simulations using isotope‐incorporated GCSRMs will improve our understanding of the atmospheric hydrological cycle. Plain Language Summary: Stable water isotopes (i.e., H218O and 1H2H16O) are useful tracers that can be used to improve our understanding of the atmospheric hydrological cycle. Isotopic ratios change during water phase changes, providing a unique constraint. In this study, we developed a new global cloud system‐resolving model equipped with stable water isotopes (NICAM‐WISO), and used this model to perform climate simulations. The model simulated the seasonal means of the atmospheric hydrological cycle, and captured geophysical and meteorological isotope effects. The performance of the model was comparable with that of other isotope‐incorporated global atmospheric models. We found that using fine‐horizontal‐resolution led to stronger spatial correlations than using coarse‐horizontal‐resolution. However, the model showed some isotopic bias. For example, in the tropical mid‐troposphere, there was negative bias in the mixing ratio and isotopic ratio of water vapor. Satellite retrieval analysis revealed that this negative bias was caused by weak mixing with boundary layer vapor and low raindrop evaporation, with a small evaporation fraction. Our results demonstrated that stable water isotopes are useful for identifying model biases associated with atmospheric humidity processes. Further fine‐horizontal‐resolution simulations using NICAM‐WISO will improve our understanding of the atmospheric hydrological cycle. Key Points: We developed a global cloud‐system‐resolving model equipped with stable water isotopes (NICAM‐WISO)Simulations performed using the model captured geophysical and meteorological isotope effects, but with some isotopic biasSimulations using NICAM‐WISO with a fine horizontal resolution will improve our understanding of the atmospheric hydrological cycle [ABSTRACT FROM AUTHOR]

Results of Russian studies on cloud physics, precipitation, and weather modification in 2019–2022 are presented based on a survey prepared for the Russian National Report on Meteorology and Atmospheric Sciences to the 28th General Assembly of the International Union of Geodesy and Geophysics. Results concerning general issues of the observation and modeling of clouds and precipitation, including convective clouds; issues of studying the microphysical and optical characteristics of clouds; and the weather modification are discussed. [ABSTRACT FROM AUTHOR]

The search for an explanation of correlations of weather and climate with solar activity has uncovered some subtle effects of the varying solar inputs, which include irradiance, energetic particles, electric fields induced by the solar wind, and the modulation of the galactic cosmic ray flux. Relevant phenomena have been identified in stratospheric chemistry and dynamics and in aerosols, atmospheric electricity, and clouds. The effects are small and appear intermittently in statistical analyses of observations. Correlations and theory permit responses to several solar inputs effective on the decadal time scales, making it difficult to distinguish them observationally. No mechanism for tropospheric decadal change can be considered to be definitely established, but for the day-to-day time scale there is clear observational evidence for clouds responding to solar wind-induced current flow (JZ) in the global electric circuit, due to changes in ionospheric potential, and independently from Forbush decreases of the cosmic ray flux. The responses to JZ also correlate separately with ionospheric potential changes due to changes in day-to-day thunderstorm activity. The identified mechanism is for electric charge effects on in-cloud scavenging, and this implies a decadal response, in view of decadal changes in JZ. However, a complete and quantitative model of the inferred electrically induced changes in cloud microphysics is not yet available. The failures, successes, and controversies of the search illustrate the somewhat chaotic process of scientific discovery. [ABSTRACT FROM AUTHOR]

Measurements from the South Dakota School of Mines and Technology T-28 hail-penetrating aircraft are analyzed using recently developed data processing techniques with the goals of identifying where the large hail is found relative to vertical motion and improving the detection of hail microphysical properties from radar. Hail particle size distributions (PSD) and environmental conditions (temperature, relative humidity, liquid water content, air vertical velocity) were digitally collected by the T28 between 1995 and 2003 and synthesized by Detwiler et al. The PSD were forward modeled by Cecchini et al. to simulate the radar reflectivity of the PSD at multiple radar wavelengths. The T-28 penetrated temperatures primarily between 08 and 2108C. The largest hailstones were sampled near the updraft/downdraft interface. Liquid water contents were highest in the updraft cores, whereas total (liquid 1 frozen) water contents were highest near the updraft/downdraft interface. The fitted properties of the PSD (intercept and slope) are directly related to each other but do not show any dependence on the region of the hailstorm where sampled. The PSD measurements and the radar reflectivity calculations at multiple radar wavelengths facilitated the development of relationships between the PSD bulk properties-hail kinetic energy and kinetic energy flux-and the radar reflectivity. Rather than using the oft-assumed sphericity and solid ice physical properties, actual measurements of hail properties are used in the analysis. Results from the maximum estimated size of hail (MESH) and vertical integrated liquid water (VIL) algorithms are evaluated based on this analysis. SIGNIFICANCE STATEMENT: Hailstorms in the United States have caused over $10 billion in damage for each of the last 14 years, according to insurance industry estimates (Heymsfield and Giammanco 2020). Algorithms have been developed to identify the presence and size of hail from radar. Numerical simulations of hailstorms have improved significantly since the 1970s, and further improvements will provide better resolution and more accurate estimates of the sizes of hailstones falling to the ground. Measurements of the properties of hailstones-their mass and terminal velocities-have been improved in recent years but in general are not incorporated in the algorithms developed for radar estimates of hail sizes or for the properties of hail used in the model simulations. This study synthesizes in situ aircraft data, computed radar backscatter cross sections, together with recent estimates of the physical characteristics of hailstones to improve the representation of hail in numerical models and quantitative assessment hail properties in storms using weather radar. [ABSTRACT FROM AUTHOR]

A new version of the Weather Research and Forecasting (WRF) double-moment 6-class (WDM6) microphysics scheme was developed based on the existing WDM6 scheme by predicting snow and graupel number concentrations. The new WDM6 scheme was tested for summer rainfall and winter snowfall cases to evaluate the effects of prognostic number concentration of snow and graupel on the simulated precipitation. The number concentration of snow decreases at the upper layers and the one of graupel also decreases at all layers in the new WDM6 scheme compared to the diagnosed ones in the original WDM6 scheme. Rain number concentration is remarkably reduced in the new WDM6 scheme due to the newly added and modified sink processes. Therefore, the new scheme produces a larger size of raindrops with a reduced number concentration than the original scheme, which hinders raindrop evaporation and produces more surface rain. Even though the enhanced surface rainfall in the new scheme deteriorates the bias score, the new scheme improves the statistical skill of the equitable threat score and probability of detection in most cases. These scores all improved for warm-type summer cases in the new scheme. The new scheme also shows more comparable features to the observation for the probability density functions of simulated liquid equivalent precipitation rates by alleviating the overprediction problem of precipitation frequencies belonging to heavy precipitation categories. Therefore, the new scheme improves the precipitation forecast for warm-type summer cases, which occur most frequently during the summer season over the Korean Peninsula. [ABSTRACT FROM AUTHOR]

Secondary ice productions (SIPs) can produce ice crystals with a number concentration much higher than that of ice nucleating particles in mixed-phase clouds and therefore influence cloud glaciation and precipitation. For midlatitude continental mesoscale convective systems (MCSs), how SIPs affect the microphysical properties and precipitation is still not clear. There are few studies of SIPs in midlatitude continental MCSs. This study investigates the roles of three SIPs (rime splintering, freezing drop shattering, and ice-ice collisional breakup) on a squall line case in North China on 18 August 2020 using the WRF model with a modified Morrison double-moment bulk microphysical scheme. Including SIPs, especially ice-ice collisional breakup, in the model simulations markedly improves the simulated convective area and convective precipitation rate of the squall line, while slightly improving the area and precipitation of the stratiform region. Within the mixed-phase layer in both the convective and stratiform regions of the squall line, ice-ice collisional breakup is the dominant process to generate ice crystals. In contrast, rime splintering generates an order of magnitude fewer ice crystals than ice-ice collisional breakup, while freezing drop shattering plays a negligible role due to the lack of large drops. Ice multiplication through ice-ice collisional breakup and rime splintering produces numerous snowflakes and graupel. This leads to enhanced depositional growth and weaker riming, which in turn weakens rime splintering. It is recommended to add SIP parameterization to the model. [ABSTRACT FROM AUTHOR]

Stratocumulus occur in closed- or open-cell states, which tend to be associated with high or low cloud cover and the absence or presence of precipitation, respectively. Thus, the transition between these states has substantial implications for the role of this cloud type in Earth's radiation budget. In this study, we analyze transitions between these states using an ensemble of 127 large-eddy simulations, covering a wide range of conditions. Our analysis is focused on the behavior of these clouds in a cloud fraction (fc) scene albedo (A) phase space, which has been shown in previous studies to be a useful framework for interpreting system behavior. For the transition from closed to open cells, we find that precipitation creates narrower clouds and scavenges cloud droplets for all fc. However, precipitation decreases the cloud depth for fc > 0.8 only, causing a rapid decrease in A. For fc < 0.8, the cloud depth actually increases due to mesoscale organization of the cloud field. As the cloud deepening balances the effects of cloud droplet scavenging in terms of influence on A, changes in A are determined by the decreasing fc only, causing a linear decrease in A for fc < 0.8. For the transition from open to closed cells, we find that longwave radiative cooling drives the cloud development, with cloud widening dominating for fc < 0.5. For fc > 0.5, clouds begin to deepen gradually due to the decreasing efficiency of lateral expansion. The smooth switch between cloud widening and deepening leads to a more gentle change in A compared to the transitions under precipitating conditions. Significance Statement: By reflecting a substantial fraction of solar shortwave radiation back to space, shallow clouds constitute a major cooling agent in Earth's radiation budget. To constrain this effect, a profound understanding of cloud cover and cloud albedo is necessary. In this study, we analyze the processes that drive the variability in these cloud properties in stratocumulus clouds, a very common cloud type covering approximately 20% of the globe. For these clouds, we show that changes from low to high or high to low cloud cover are different due to the underlying cloud micro- and macrophysics, elucidating this crucial aspect of aerosol–cloud–climate interactions. [ABSTRACT FROM AUTHOR]

Blue corona discharges are often observed at the top of thunderclouds. They are bursts of streamers, but the cloud conditions that enable them are not well known. Here we present observations by the Atmosphere‐Space Interactions Monitor (ASIM) of 92 discharges during its ∼1 min pass over tropical cyclone Fani in the Bay of Bengal from 20:10:55 to 20:12:05 UTC on 30 April 2019. The discharges were observed in convective cells forming in the rainbands of the cyclone where Convective Available Potential Energy reached ∼6,000 J kg−1. The Cloud‐Aerosol Lidar and Infrared Pathfinder Satellite Observation satellite passed over one of the cells ∼12 min after ASIM from 20:23:58 to 20:24:14 UTC. It measured the cloud microphysics related to the discharges and indicated they occurred in a convection region with the cloud top overshooting for over 20 min. The updraft lifted ice particles to lower stratospheric altitudes and formed the gullwing‐shaped cirrus. The discharges are found at an average altitude of ∼16 km where the cloud environment contained ∼2 × 107 m−3 ice particles with ∼50 µm radius, resulting in a photon mean free path of ∼3 m. Around 20% of the blue corona discharges coincide with Narrow Bipolar Events indentified from the Earth Networks Total Lightning Network. Our observations suggest that the overshooting cloud top formed by deep convection and a surge in lightning activity facilitated conditions for the blue corona discharges. This work provides the first‐ever estimate of important microphysical parameters related to blue corona discharges based on data measurements, establishing a reference for future empirical and theoretical studies. Plain Language Summary: Blue corona discharges are bursts of streamer discharges often observed at the top of thunderclouds, but the conditions in the clouds that generate them are not well understood. In this study, we discuss observations of blue corona discharges in convective cells detected by the Atmosphere‐Space Interactions Monitor on the International Space Station as it passed over cyclone Fani in the Bay of Bengal. For the first time, we have observations of the cloud particle characteristics at the cloud tops taken shortly afterward by the Cloud‐Aerosol Lidar and Infrared Pathfinder Satellite Observation satellite, which are related to the blue corona discharges. The observations indicate that the blue corona discharges may be Narrow Bipolar Events associated with strong convection in cloud cells reaching into the stratosphere. The cloud microphysical parameters related to blue corona discharges are important for future empirical studies and for theoretical models and simulations. Key Points: First combined observations of optical activity, cloud microphysics and lightning related to blue corona discharges in a tropical cycloneThe source depth of the blue corona discharges inside the cloud estimated from optical and radio observations agree wellBlue corona discharges are related to a strong updraft, high‐density large ice particles and an increase in the total lightning rate [ABSTRACT FROM AUTHOR]

The National Center for Atmospheric Research released a global atmosphere model named Community Atmosphere Model version 5.0 (CAM5), which aimed to provide a global climate simulation for meteorological research. Among them, the cloud microphysics scheme is extremely time-consuming, so developing efficient parallel algorithms faces large-scale and chronic simulation challenges. Due to the wide application of GPU in the fields of science and engineering and the NVIDIA's mature and stable CUDA platform, we ported the code to GPU to accelerate computing. In this paper, by analyzing the parallelism of CAM5 cloud microphysical schemes (CAM5 CMS) in different dimensions, corresponding GPU-based one-dimensional (1D) and two-dimensional (2D) parallel acceleration algorithms are proposed. Among them, the 2D parallel algorithm exploits finer-grained parallelism. In addition, we present a data transfer optimization method between the CPU and GPU to further improve the overall performance. Finally, GPU version of the CAM5 CMS (GPU-CMS) was implemented. The GPU-CMS can obtain a speedup of 141.69 × on a single NVIDIA A100 GPU with I/O transfer. In the case without I/O transfer, compared to the baseline performance on a single Intel Xeon E5-2680 CPU core, the 2D acceleration algorithm obtained a speedup of 48.75 × , 280.11 × , and 507.18 × on a single NVIDIA K20, P100, and A100 GPU, respectively. [ABSTRACT FROM AUTHOR]