Search

Cloud seeding has been used as one water management strategy to overcome the increasing demand for water despite decades of inconclusive results on the efficacy of cloud seeding. In this study snowfall accumulation from glaciogenic cloud seeding is quantified based on snow gauge and radar observations from three days in January 2017, when orographic clouds in the absent of natural precipitation were seeded with silver iodide (AgI) in the Payette basin of Idaho during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE). On each day, a seeding aircraft equipped with AgI flares flew back and forth on a straight-line flight track producing a zig-zag pattern representing two to eight lines of clouds visible through enhancements in radar reflectivity. As these seeding lines started to form precipitation, they passed over several snow gauges and through the radar observational domain. For the three cases presented here, precipitation gauges measured increases between 0.05-0.3 mm as precipitation generated by cloud seeding pass over the instruments. A variety of relationships between radar reflectivity factor and liquid equivalent snowfall rate were used to quantify snowfall within the radar observation domain. For the three cases, snowfall occurred within the radar observational domain between 25 -160 min producing a total amount of water generated by cloud seeding ranging from 123,220 to 339,540 m3 using the best-match Ze-S relationship. Uncertainties in radar reflectivity estimated snowfall are provided by considering not only the best-match Ze-S relationship but also an ensemble of Ze-S relationships based on the range of coefficients published from previous studies and then examining the percentile of snowfall estimates based on all of the Ze-S relationships within the ensemble. Considering the interquartile range and 5th/95th percentiles, uncertainties in total amount of water generated by cloud seeding can range between 20-45% compared to the best-math estimates. These results provide new insights towards understanding how cloud seeding impacts precipitation and its distribution across a region.

Supercooled large drop (SLD) icing poses a unique hazard for aircraft and has resulted in new regulations regarding aircraft certification to fly in regions of known or forecast SLD icing conditions. The new regulations define two SLD icing categories based upon the maximum supercooled liquid water drop diameter (Dmax): freezing drizzle (100–500 μm) and freezing rain (> 500 μm). Recent upgrades to U.S. operational numerical weather prediction models lay a foundation to provide more relevant aircraft icing guidance including the potential to predict explicit drop size. The primary focus of this paper is to evaluate a proposed method for estimating the maximum drop size from model forecast data to differentiate freezing drizzle from freezing rain conditions. Using in-situ cloud microphysical measurements collected in icing conditions during two field campaigns between January and March 2017, this study shows that the High-Resolution Rapid Refresh model is capable of distinguishing SLD icing categories of freezing drizzle and freezing rain using a Dmax extracted from the rain category of the microphysics output. It is shown that the extracted Dmax from the model correctly predicted the observed SLD icing category as much as 99% of the time when the HRRR accurately forecast SLD conditions; however, performance varied by the method to define Dmax and by the field campaign dataset used for verification.

Aerosol and cloud microphysical measurements were collected by a research aircraft during August 2019 over the United Arab Emirates (UAE). The majority of scientific flights targeted summertime convection along the eastern Al Hajar Mountains bordering Oman, while one flight sampled non-orographic clouds over the western UAE near the Saudi Arabian border. In this work, we study the evolution of growing cloud turrets from cloud base (9 ∘C) up to the capping inversion level (−12 ∘C) using coincident cloud particle imagery and particle size distributions from cloud cores under different forcing. Results demonstrate the active role of background dust and pollution as cloud condensation nuclei (CCN) with the onset of their deliquescence in the subcloud region. Subcloud aerosol sizes are shown to extend from submicron to 100 µm sizes, with higher concentrations of ultra-giant CCN (d>10 µm) from local sources closer to the Saudi border, compared with the eastern orographic region where smaller CCN are observed. Despite the presence of ultra-giant CCN from dust and pollution in both regions, an active collision–coalescence (C–C) process is not observed within the limited depths of warm cloud ( m). The state-of-the-art observations presented in this paper can be used to initialize modeling case studies to examine the influence of aerosols on cloud and precipitation processes in the region and to better understand the impacts of hygroscopic cloud seeding on these clouds.

In this two-part paper, influences from environmental factors on lightning in a convective storm are assessed with a model. In Part I, an electrical component is described and applied in the Aerosol–Cloud model (AC). AC treats many types of secondary (e.g., breakup in ice–ice collisions, raindrop-freezing fragmentation, rime splintering) and primary (heterogeneous, homogeneous freezing) ice initiation. AC represents lightning flashes with a statistical treatment of branching from a fractal law constrained by video imagery.The storm simulated is from the Severe Thunderstorm Electrification and Precipitation Study (STEPS; 19/20 June 2000). The simulation was validated microphysically [e.g., ice/droplet concentrations and mean sizes, liquid water content (LWC), reflectivity, surface precipitation] and dynamically (e.g., ascent) in our 2017 paper. Predicted ice concentrations (~10 L−1) agreed—to within a factor of about 2—with aircraft data at flight levels (−10° to −15°C). Here, electrical statistics of the same simulation are compared with observations. Flash rates (to within a factor of 2), triggering altitudes and polarity of flashes, and electric fields, all agree with the coincident STEPS observations.The “normal” tripole of charge structure observed during an electrical balloon sounding is reproduced by AC. It is related to reversal of polarity of noninductive charging in ice–ice collisions seen in laboratory experiments when temperature or LWC are varied. Positively charged graupel and negatively charged snow at most midlevels, charged away from the fastest updrafts, is predicted to cause the normal tripole. Total charge separated in the simulated storm is dominated by collisions involving secondary ice from fragmentation in graupel–snow collisions.

Throughout the western United States and other semiarid mountainous regions across the globe, water supplies are fed primarily through the melting of snowpack. Growing populations place higher demands on water, while warmer winters and earlier springs reduce its supply. Water managers are tantalized by the prospect of cloud seeding as a way to increase winter snowfall, thereby shifting the balance between water supply and demand. Little direct scientific evidence exists that confirms even the basic physical hypothesis upon which cloud seeding relies. The intent of glaciogenic seeding of orographic clouds is to introduce aerosol into a cloud to alter the natural development of cloud particles and enhance wintertime precipitation in a targeted region. The hypothesized chain of events begins with the introduction of silver iodide aerosol into cloud regions containing supercooled liquid water, leading to the nucleation of ice crystals, followed by ice particle growth to sizes sufficiently large such that snow falls to the ground. Despite numerous experiments spanning several decades, no direct observations of this process exist. Here, measurements from radars and aircraft-mounted cloud physics probes are presented that together show the initiation, growth, and fallout to the mountain surface of ice crystals resulting from glaciogenic seeding. These data, by themselves, do not address the question of cloud seeding efficacy, but rather form a critical set of observations necessary for such investigations. These observations are unambiguous and provide details of the physical chain of events following the introduction of glaciogenic cloud seeding aerosol into supercooled liquid orographic clouds.

In Part I of this two-part paper, a formulation was developed to treat fragmentation in ice–ice collisions. In the present Part II, the formulation is implemented in two microphysically advanced cloud models simulating a convective line observed over the U.S. high plains. One model is 2D with a spectral bin microphysics scheme. The other has a hybrid bin–two-moment bulk microphysics scheme in 3D. The case consists of cumulonimbus cells with cold cloud bases (near 0°C) in a dry troposphere. Only with breakup included in the simulation are aircraft observations of particles with maximum dimensions >0.2 mm in the storm adequately predicted by both models. In fact, breakup in ice–ice collisions is by far the most prolific process of ice initiation in the simulated clouds (95%–98% of all nonhomogeneous ice), apart from homogeneous freezing of droplets. Inclusion of breakup in the cloud-resolving model (CRM) simulations increased, by between about one and two orders of magnitude, the average concentration of ice between about 0° and −30°C. Most of the breakup is due to collisions of snow with graupel/hail. It is broadly consistent with the theoretical result in Part I about an explosive tendency for ice multiplication. Breakup in collisions of snow (crystals >~1 mm and aggregates) with denser graupel/hail was the main pathway for collisional breakup and initiated about 60%–90% of all ice particles not from homogeneous freezing, in the simulations by both models. Breakup is predicted to reduce accumulated surface precipitation in the simulated storm by about 20%–40%.

Two ground-based and two airborne cloud seeding events between 2003 and 2005 from a trace chemistry field experiment conducted by the Desert Research Institute (DRI) in the Payette river basin of Idaho are simulated by the Weather Research and Forecasting (WRF) model in Large Eddy Simulation (LES) mode with a grid spacing of 667-m using the recently developed Wintertime silver iodide ( AgI ) Seeding Parameterization (WASP, Xue et al., 2013a,b ). The WASP in WRF was specifically modified to simulate the emissions of both AgI and Indium oxide ( In 2 O 3 ) particles from ground- and aircraft-based generators and the associated interactions with hydrometeors. Validations of model results against available observations from soundings and SNOwpack TELemetry (SNOTEL) data show that WRF simulations reasonably capture the dynamics, the thermodynamics, and the precipitation patterns. The comparisons of deposited Ag and In concentrations between model results and observations, and the analyses of seeding effects and AgI seeding efficiencies show that: The simulated Ag mass concentrations inside the snow are similar to observations ranging mostly from 0 to 10 −9 kg kg −1 and have high spatial correlations with observations in ground seeding cases (∼ 80 % ) but lower in airborne seeding (∼ 40 % ). The simulated Ag / In ratios range between 60 and 600 indicating that ice nucleation of AgI particle dominates the scavenging of AgI in incorporating Ag into the snow. The overall simulated seeding effect of the four cases is 6 % increase in precipitation over the LES domain, which is similar to what was inferred from the observations for the entire experiment. The magnitude of the AgI seeding efficiency is one order higher in airborne seeding cases (∼ 10 10 kg kg −1 ) compared to ground seeding cases (∼ 10 9 kg kg −1 ).

The impacts of various assumptions of cloud properties represented within a numerical weather prediction model's radiation scheme are demonstrated. In one approach, the model assumed the radiative effective radii of cloud water, cloud ice, and snow were represented by values assigned a priori, whereas a second, “coupled” approach utilized known cloud particle assumptions in the microphysics scheme that evolved during the simulations to diagnose the radii explicitly. This led to differences in simulated infrared (IR) brightness temperatures, radiative fluxes through clouds, and resulting surface temperatures that ultimately affect model-predicted diurnally-driven convection. The combined approach of evaluating simulated versus observed IR brightness temperatures, radiation reaching the ground, and surface temperature forecasts revealed the root model biases better than evaluating any single variable. This study found that the Weather Research and Forecasting (WRF) model predicted less overall clouds than was observed, particularly in the mid‑troposphere, but that properly connecting the assumptions of particle sizes in the microphysics scheme to the radiation scheme resulted in sensible cloud-radiation indirect effects and modest improvements in simulated IR brightness temperature, amount of solar radiation reaching the ground, and surface temperature.

A new microphysics scheme has been developed based on the prediction of bulk particle properties for a single ice-phase category, in contrast to the traditional approach of separating ice into various predefined species (e.g., cloud ice, snow, and graupel). In this paper, the new predicted particle properties (P3) scheme, described in Part I of this series, is tested in three-dimensional simulations using the Weather Research and Forecasting (WRF) Model for two contrasting well-observed cases: a midlatitude squall line and winter orographic precipitation. Results are also compared with simulations using other schemes in WRF. Simulations with P3 can produce a wide variety of particle characteristics despite having only one free ice-phase category. For the squall line, it produces dense, fast-falling, hail-like ice near convective updraft cores and lower-density, slower-falling ice elsewhere. Sensitivity tests show that this is critical for simulating high precipitation rates observed along the leading edge of the storm. In contrast, schemes that represent rimed ice as graupel, with lower fall speeds than hail, produce lower peak precipitation rates and wider, less distinct, and less realistic regions of high convective reflectivity. For the orographic precipitation case, P3 produces areas of relatively fast-falling ice with characteristics of rimed snow and low- to medium-density graupel on the windward slope. This leads to less precipitation on leeward slopes and more on windward slopes compared to the other schemes that produce large amounts of snow relative to graupel (with generally the opposite for schemes with significant graupel relative to snow). Overall, the new scheme produces reasonable results for a reduced computational cost.

Four cloud-seeding cases over southern Idaho during the 2010/11 winter season have been simulated by the Weather Research and Forecasting (WRF) model using the coupled silver iodide (AgI) cloud-seeding scheme that was described in Part I. The seeding effects of both ground-based and airborne seeding as well as the impacts of model physics, seeding rates, location, timing, and cloud properties on seeding effects have been investigated. The results were compared with those from Part I and showed the following: 1) For the four cases tested in this study, control simulations driven by the Real-Time Four Dimensional Data Assimilation (RTFDDA) WRF forecast data generated more realistic atmospheric conditions and precipitation patterns than those driven by the North America Regional Reanalysis data. Sensitivity experiments therefore used the RTFDDA data. 2) Glaciogenic cloud seeding increased orographic precipitation by less than 1% over the simulation domain, including the Snake River basin, and by up to 5% over the target areas. The local values of the relative precipitation enhancement by seeding were ~20%. Most of the enhancement came from vapor depletion. 3) The seeding effect was inversely related to the natural precipitation efficiency but was positively related to seeding rates. 4) Airborne seeding is generally more efficient than ground-based seeding in terms of targeting, but its efficiency depends on local meteorological conditions. 5) The normalized seeding effects ranged from 0.4 to 1.6 under various conditions for a certain seeding event.

A silver iodide (AgI) cloud-seeding parameterization has been implemented into the Thompson microphysics scheme of the Weather Research and Forecasting model to investigate glaciogenic cloud-seeding effects. The sensitivity of the parameterization to meteorological conditions, cloud properties, and seeding rates was examined by simulating two-dimensional idealized moist flow over a bell-shaped mountain. The results verified that this parameterization can reasonably simulate the physical processes of cloud seeding with the limitations of the constant cloud droplet concentration assumed in the scheme and the two-dimensional model setup. The results showed the following: 1) Deposition was the dominant nucleation mode of AgI from simulated aircraft seeding, whereas immersion freezing was the most active mode for ground-based seeding. Deposition and condensation freezing were also important for ground-based seeding. Contact freezing was the weakest nucleation mode for both ground-based and airborne seeding. 2) Diffusion and riming on AgI-nucleated ice crystals depleted vapor and liquid water, resulting in more ice-phase precipitation on the ground for all of the seeding cases relative to the control cases. Most of the enhancement came from vapor depletion. The relative enhancement by seeding ranged from 0.3% to 429% under various conditions. 3) The maximum local AgI activation ratio was 60% under optimum conditions. Under most seeding conditions, however, this ratio was between 0.02% and 2% in orographic clouds. 4) The seeding effect was inversely related to the natural precipitation efficiency but was positively related to seeding rates. 5) Ground-based seeding enhanced precipitation on the lee side of the mountain, whereas airborne seeding from lower flight tracks enhanced precipitation on the windward side of the mountain.

Although it is generally accepted that the level of neutral buoyancy (LNB) is only a coarse estimate of updraft depth, the LNB is still used to understand and predict storm structure in both observations and modeling. This study uses case studies to quantify the variability associated with using environmental soundings to predict detrainment levels. Nine dual-Doppler convective cases were used to determine the observed level of maximum detrainment (LMD) to compare with the LNB. The LNB for each case was calculated with a variety of methods and with a variety of sources (including both observed and simulated soundings). The most representative LNB was chosen as the proximity sounding from NARR using the most unstable parcel and including ice processes. The observed cases were a mix of storm morphologies, including both supercell and multicell storms. As expected, the LMD was generally below the LNB, the mean offset for all cases being 2.2 km. However, there was a marked difference between the supercell and non-supercell cases. The two supercell cases had LMDs of 0.3 km and 0.0 km below the LNB. The remaining cases had LMDs that ranged from 4.0 km below to 1.6 km below the LNB, with a mean offset of 2.8 km below. Observations also showed that evolution of the LMD over the lifetime of the storm can be significant (e.g., >2 km altitude change in 30 min), and this time evolution is lacking from models with coarse time steps, missing significant changes in detrainment levels that may strongly impact the amount of boundary layer mass transported to the upper troposphere and lower stratosphere.

This paper describes idealized simulations of a squall line observed on 20 June 2007, in central Oklahoma. Results are compared with measurements from dual-polarization radar and surface disdrometer. The baseline model configuration qualitatively reproduces key storm features, but underpredicts precipitation rates and generally overpredicts median volume raindrop diameter. The sensitivity of model simulations to parameterization of raindrop breakup is tested under different low-level (0–2.5 km) environmental vertical wind shears. Storm characteristics exhibit considerable sensitivity to the parameterization of breakup, especially for moderate (0.0048 s−1) shear. Simulations with more efficient breakup tend to have higher domain-mean precipitation rates under both moderate and higher (0.0064 s−1) shear, despite the smaller mean drop size and hence lower mass-weighted fall speed and higher evaporation rate for a given rainwater content. In these runs, higher evaporation leads to stronger cold pools, faster propagation, larger storm size, greater updraft mass flux (but weaker convective updrafts at mid- and upper levels), and greater total condensation that compensates for the increased evaporation to give more surface precipitation. The impact of drop breakup on mass-weighted fall speed is also important and leads to a nonmonotonic response of storm characteristics (surface precipitation, cold pool strength, etc.) to changes in breakup efficiency under moderate wind shear. In contrast, the response is generally monotonic at higher wind shear. Interactions between drop breakup, convective dynamics, cold pool intensity, and low-level environmental wind shear are also described in the context of “Rotunno–Klemp–Weisman (RKW) theory,” which addresses how density currents evolve in sheared environments.

This study addresses the kinematic, microphysical, and electrical evolution of an isolated convective storm observed on 3 June 2000 during the Severe Thunderstorm Electrification and Precipitation Study field campaign. Doppler-derived vertical velocities, radar reflectivity, hydrometeor classifications from polarimetric radar, and Lightning Mapping Array (LMA) charge structures are examined over a nearly 3-h period. This storm, characterized as a low-precipitation supercell, produced modest amounts of hail, determined by fuzzy-logic hydrometeor classification as mostly small (

This study discusses radar and lightning observations of two multicellular storms observed during the Severe Thunderstorm Electrification and Precipitation Study. The Lightning Mapping Array data indicated that the charge structure of the 19 June 2000 storm was consistent with a normal polarity tripole, while the 22 June 2000 storm exhibited an overall inverted tripolar charge structure. The 19 June storm consisted of weaker convection and produced little to no hail and moderate total flash rates peaking between 80 and 120 min−1. The cells in the 22 June 2000 storm were much more vigorous, exhibited strong, broad updrafts, and produced large quantities of hail, as well as extraordinary total flash rates as high as 500 min−1. The National Lightning Detection Network (NLDN) indicated that the 19 June storm produced mostly negative cloud-to-ground (CG) lightning, while the 22 June storm produced predominantly positive CG lightning, peaking at 10 min−1 just after two cells merged. However, the Los Alamos Sferic Array indicated that many of the positive CG strokes reported by the NLDN in the 22 June storm were intracloud discharges known as narrow bipolar events. Negative CG lightning was also observed in the 22 June storm, but typically came to ground beneath an inverted dipole in the storm anvil.