Search

This study investigates the structure and evolution of a summertime convective event that occurred on 14 July 2015 over the Arabian region. We use the WRF Model with 1-km horizontal grid spacing and test three PBL parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme; the Asymmetrical Convective Model, version 2, (ACM2) scheme; and the quasi-normal scale-elimination (QNSE) scheme. Convection initiates near the Al Hajar Mountains of northern Oman at around 1100 local time (LT; 0700 UTC) and propagates northwestward. A nonorographic convective band along the west coast of the United Arab Emirates (UAE) develops after 1500 LT as a result of the convergence of cold pools with the sea breeze from the Arabian Gulf. The model simulation employing the QNSE scheme simulates the convection initiation and propagation well. Although the MYNN and ACM2 simulations show convective initiation near the Al Hajar Mountains, they fail to simulate the development of the convective band along the UAE west coast. The MYNN run simulates colder near-surface temperatures and a weaker sea breeze, whereas the ACM2 run simulates a stronger sea breeze but a drier lower troposphere. Sensitivity simulations using horizontal grid spacings of 9 and 3 km show that lower-resolution runs develop broader convective structures and weaker cold pools and horizontal wind divergence, affecting the development of convection along the west coast of the UAE. The 1-km run using the QNSE PBL scheme realistically captures the sequence of events that leads to the moist convection over the UAE and adjacent mountains.

Updrafts in wintertime cloud systems over mountainous regions can be described as fixed, mechanically driven by the terrain under a given ambient wind and stability profile (i.e., vertically propagating gravity waves tied to flow over topography), and transient, associated primarily with vertical wind shear and conditional instability within passing weather systems. This analysis quantifies the magnitude of fixed and transient updraft structures over the Payette River basin sampled during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE). Vertical motions were retrieved from Wyoming Cloud Radar measurements of radial velocity using the algorithm presented in Part I. Transient circulations were removed, and fixed orographic circulations were quantified by averaging vertical circulations along repeated cross sections over the same terrain during the campaign. Fixed orographic vertical circulations had magnitudes of 0.3–0.5 m s−1. These fixed vertical circulations were composed of a background circulation in which transient circulations were embedded. Transient vertical circulations are shown to be associated with transient wave motions, cloud-top generating cells, convection, and turbulence. Representative transient vertical circulations are illustrated, and data from rawinsondes over the Payette River basin are used to infer the relationship of the vertical circulations to shear and instability. Maximum updrafts are shown to exceed 5 m s−1 within Kelvin–Helmholtz waves, 4 m s−1 associated with transient gravity waves, 3 m s−1 in generating cells, 6 m s−1 in elevated convection, 4 m s−1 in surface-based deep convection, 5 m s−1 in boundary layer turbulence, and 9 m s−1 in shear-induced turbulence.

Vertical motions over the complex terrain of Idaho’s Payette River basin were observed by the Wyoming Cloud Radar (WCR) during 23 flights of the Wyoming King Air during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE) field campaign. The WCR measured radial velocity Vr, which includes the reflectivity-weighted terminal velocity of hydrometeors Vt, vertical air velocity w, horizontal wind contributions as a result of aircraft attitude deviations, and aircraft motion. Aircraft motion was removed through standard processing. To retrieve vertical radial velocity W, Vr was corrected using rawinsonde data and aircraft attitude measurements; w was then calculated by subtracting the mean W () at a given height along a flight leg long enough for to equal the mean reflectivity-weighted terminal velocity at that height. The accuracy of the w and retrievals were dependent on satisfying assumptions along a given flight leg that the winds at a given altitude above/below the aircraft did not vary, the vertical air motions at a given altitude sum to 0 m s−1, and at a given altitude did not vary. The uncertainty in the w retrieval associated with each assumption is evaluated. Case studies and a projectwide summary show that this methodology can provide estimates of w that closely match gust probe measurements of w at the aircraft level. Flight legs with little variation in equivalent reflectivity factor at a given height and large horizontal echo extent were associated with the least retrieval uncertainty. The greatest uncertainty occurred in regions with isolated convective turrets or at altitudes where split cloud layers were present.

In Part II, two classes of vertical motions, fixed (associated with vertically propagating gravity waves tied to flow over topography) and transient (associated primarily with vertical wind shear and conditional instability within passing weather systems), were diagnosed over the Payette River basin of Idaho during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE). This paper compares vertical motions retrieved from airborne Doppler radial velocity measurements with those from a 900-m-resolution model simulation to determine the impact of transient vertical motions on trajectories of ice particles initiated by airborne cloud seeding. An orographic forcing index, developed to compare vertical motion fields retrieved from the radar with the model, showed that fixed vertical motions were well resolved by the model while transient vertical motions were not. Particle trajectories were calculated for 75 cross-sectional pairs, each differing only by the observed and modeled vertical motion field. Wind fields and particle terminal velocities were otherwise identical in both trajectories so that the impact of transient vertical circulations on particle trajectories could be isolated. In 66.7% of flight-leg pairs, the distance traveled by particles in the model and observations differed by less than 5 km with transient features having minimal impact. In 9.3% of the pairs, model and observation trajectories landed within the ideal target seeding elevation range (>2000 m), whereas, in 77.3% of the pairs, both trajectories landed below the ideal target elevation. Particles in the observations and model descended into valleys on the mountains’ lee sides in 94.2% of cases in which particles traveled less than 37 km.

A dry-air intrusion induced by the tropopause folding split the deep cloud into two layers resulting in a shallow orographic cloud with a supercooled liquid cloud top at around −15°C and an ice cloud above it on 19 January 2017 during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE). The airborne AgI seeding of this case was simulated by the WRF Weather Modification (WRF-WxMod) Model with different configurations. Simulations at different grid spacing, driven by different reanalysis data, using different model physics were conducted to explore the ability of WRF-WxMod to capture the properties of natural and seeded clouds. The detailed model–observation comparisons show that the simulation driven by ERA5 data, using Thompson–Eidhammer microphysics with 30% of the CCN climatology, best captured the observed cloud structure and supercooled liquid water properties. The ability of the model to correctly capture the wind field was critical for successful simulation of the seeding plume locations. The seeding plume features and ice number concentrations within them from the large-eddy simulations (LES) are in better agreement with observations than non-LES runs mostly due to weaker AgI dispersion associated with the finer grid spacing. Seeding effects on precipitation amount and impacted areas from LES seeding simulations agreed well with radar-derived values. This study shows that WRF-WxMod is able to simulate and quantify observed features of natural and seeded clouds given that critical observations are available to validate the model. Observation-constrained seeding ensemble simulations are proposed to quantify the AgI seeding impacts on wintertime orographic clouds.Significance StatementRecent observational work has demonstrated that the impact of airborne glaciogenic seeding of orographic supercooled liquid clouds is detectable and can be quantified in terms of the extra ground precipitation. This study aims, for the first time, to simulate this seeding impact for one well-observed case. The stakes are high: if the model performs well in this case, then seasonal simulations can be conducted with appropriate configurations after validations against observations, to determine the impact of a seeding program on the seasonal mountain snowpack and runoff, with more fidelity than ever. High–resolution weather simulations inherently carry uncertainty. Within the envelope of this uncertainty, the model compares very well to the field observations.

The spatial distribution and magnitude of snowfall resulting from cloud seeding with silver iodide (AgI) is closely linked to atmospheric conditions, seeding operations, and dynamical, thermodynamical, and microphysical processes. Here, microphysical processes leading to ice and snow production are analyzed in orographic clouds for three cloud-seeding events, each with light or no natural precipitation and well-defined, traceable seeding lines. Airborne and ground-based radar observations are linked to in situ cloud and precipitation measurements to determine the spatiotemporal evolution of ice initiation, particle growth, and snow fallout in seeded clouds. These processes and surface snow amounts are explored as particle plumes evolve from varying amounts of AgI released, and within changing environmental conditions, including changes in liquid water content (LWC) along and downwind of the seeding track, wind speed, and shear. More AgI did not necessarily produce more liquid equivalent snowfall (LESnow). The greatest amount of LESnow, largest area covered by snowfall, and highest peak snowfall produced through seeding occurred on the day with the largest and most widespread occurrence of supercooled drizzle, highest wind shear, and greater LWC along and downwind of the seeding track. The day with the least supercooled drizzle and the lowest LWC downwind of the seeding track produced the smallest amount of LESnow through seeding. The stronger the wind was, the farther away the snowfall occurred from the seeding track.

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.

This paper reviews research conducted over the last six decades to understand and quantify the efficacy of wintertime orographic cloud seeding to increase winter snowpack and water supplies within a mountain basin. The fundamental hypothesis underlying cloud seeding as a method to enhance precipitation from wintertime orographic cloud systems is that a cloud’s natural precipitation efficiency can be enhanced by converting supercooled water to ice upstream and over a mountain range in such a manner that newly created ice particles can grow and fall to the ground as additional snow on a specified target area. The review summarizes the results of physical, statistical, and modeling studies aimed at evaluating this underlying hypothesis, with a focus on results from more recent experiments that take advantage of modern instrumentation and advanced computation capabilities. Recent advances in assessment and operations are also reviewed, and recommendations for future experiments, based on the successes and failures of experiments of the past, are given.

The Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE) project aims to study the impacts of cloud seeding on winter orographic clouds. The field campaign took place in Idaho between 7 January and 17 March 2017 and employed a comprehensive suite of instrumentation, including ground-based radars and airborne sensors, to collect in situ and remotely sensed data in and around clouds containing supercooled liquid water before and after seeding with silver iodide aerosol particles. The seeding material was released primarily by an aircraft. It was hypothesized that the dispersal of the seeding material from aircraft would produce zigzag lines of silver iodide as it dispersed downwind. In several cases, unambiguous zigzag lines of reflectivity were detected by radar, and in situ measurements within these lines have been examined to determine the microphysical response of the cloud to seeding. The measurements from SNOWIE aim to address long-standing questions about the efficacy of cloud seeding, starting with documenting the physical chain of events following seeding. The data will also be used to evaluate and improve computer modeling parameterizations, including a new cloud-seeding parameterization designed to further evaluate and quantify the impacts of cloud seeding.

The Wyoming Weather Modification Pilot Project randomized cloud seeding experiment was a crossover statistical experiment conducted over two mountain ranges in eastern Wyoming and lasted for 6 years (2008–13). The goal of the experiment was to determine if cloud seeding of orographic barriers could increase snowfall and snowpack. The experimental design included triply redundant snow gauges deployed in a target–control configuration, covariate snow gauges to account for precipitation variability, and ground-based seeding with silver iodide (AgI). The outcomes of this experiment are evaluated with the statistical–physical experiment design and with ensemble modeling. The root regression ratio (RRR) applied to 118 experimental units provided insufficient statistical evidence (p value of 0.28) to reject the null hypothesis that there was no effect from ground-based cloud seeding. Ensemble modeling estimates of the impact of ground-based seeding provide an alternate evaluation of the 6-yr experiment. The results of the model ensemble approach with and without seeding estimated a mean enhancement of precipitation of 5%, with an inner-quartile range of 3%–7%. Estimating the impact on annual precipitation over these mountain ranges requires results from another study that indicated that approximately 30% of the annual precipitation results from clouds identified as seedable within the seeding experiment. Thus the seeding impact is on the order of 1.5% of the annual precipitation, compared to 1% for the statistical–physical experiment, which was not sufficient to reject the null hypothesis. These results provide an estimate of the impact of ground-based cloud seeding in the Sierra Madre and Medicine Bow Mountains in Wyoming that accounts for uncertainties in both initial conditions and model physics.

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.

As a response to extreme water shortages in southeast Queensland, Australia, brought about by reduced rainfall and increasing population, the Queensland government decided to explore the potential for cloud seeding to enhance rainfall. The Queensland Cloud Seeding Research Program (QCSRP) was conducted in the southeast Queensland region near Brisbane during the 2008/09 wet seasons. In addition to conducting an initial exploratory, randomized (statistical) cloud seeding study, multiparameter radar measurements and in situ aircraft microphysical data were collected. This comprehensive set of observational platforms was designed to improve the physical understanding of the effects of both ambient aerosols and seeding material on precipitation formation in southeast Queensland clouds. This focus on gaining physical understanding, along with the unique combination of modern observational platforms utilized in the program, set it apart from previous cloud seeding research programs. The overarching goals of the ...

During the Queensland Cloud Seeding Research Program, the “CP2” polarimetric radar parameter differential radar reflectivity Zdr was used to examine the raindrop size evolution in both maritime and continental clouds. The focus of this paper is to examine the natural variability of the drop size distribution. The primary finding is that there are two basic raindrop size evolutions, one associated with continental air masses characterized by relatively high aerosol concentrations and long air trajectories over land and the other associated with maritime air masses with lower aerosol concentrations. The size evolution difference is during the growth stage of the radar echoes. The differential radar reflectivity in the growing continental clouds is dominated by large raindrops, whereas in the maritime clouds differential reflectivity is dominated by small raindrops and drizzle. The drop size evolution in many of the maritime air masses was very similar to those observed in the maritime air of the Caribbean Sea observed with the NCAR S-band polarimetric radar (S-Pol) during the Rain in Cumulus over the Ocean (RICO) experiment. Because the tops of the Queensland continental clouds ascended almost 2 times as fast as the maritime ones in their growth stage, both dynamical and aerosol factors may be important for the systematic difference in drop size evolution. Recommendations are advanced for future field programs to understand better the causes for the observed variability in drop size evolution. Also, considering the natural variability in drop size evolution, comments are provided on conducting and evaluating cloud seeding experiments.

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.

This second part of a two-part study examines the lightning and charge structure evolution of the 29 June 2000 tornadic supercell observed during the Severe Thunderstorm Electrification and Precipitation Study (STEPS). Data from the National Lightning Detection Network and the New Mexico Tech Lightning Mapping Array (LMA) are used to quantify the total and cloud-to-ground (CG) flash rates. Additionally, the LMA data are used to infer gross charge structure and to determine the origin locations and charge regions involved in the CG flashes. The total flash rate reached nearly 300 min−1 and was well correlated with radar-inferred updraft and graupel echo volumes. Intracloud flashes accounted for 95%–100% of the total lightning activity during any given minute. Nearly 90% of the CG flashes delivered a positive charge to ground (+CGs). The charge structure during the first 20 min of this storm consisted of a midlevel negative charge overlying lower positive charge with no evidence of an upper positive charge. The charge structure in the later (severe) phase was more complex but maintained what could be roughly described as an inverted tripole, dominated by a deep midlevel (5–9 km MSL) region of positive charge. The storm produced only two CG flashes (both positive) in the first 2 h of lightning activity, both of which occurred during a brief surge in updraft and hail production. Frequent +CG flashes began nearly coincident with dramatic increases in storm updraft, hail production, total flash rate, and the formation of an F1 tornado. The +CG flashes tended to cluster in or just downwind of the heaviest precipitation, which usually contained hail. The +CG flashes all originated between 5 and 9 km MSL, centered at 6.8 km (−10°C), and tapped LMA-inferred positive charge both in the precipitation core and (more often) in weaker reflectivity extending downwind. All but one of the −CG flashes originated from >9 km MSL and tended to strike near the precipitation core.

This is a two-part study that addresses the kinematic, microphysical, and electrical aspects of a severe storm that occurred in western Kansas on 29 June 2000 observed during the Severe Thunderstorm Electrification and Precipitation Study (STEPS) field campaign. In this first part, polarimetric and Doppler radar data are used along with a simple particle growth model to examine the evolution of the kinematic and microphysical properties of the storm from its earliest developing phase through its mature and dissipating phases. During its severe stage, the storm exhibited frequent positive cloud-to-ground lightning strikes, very large (∼5 cm) hail, and a tornado. Doppler-derived winds, radar reflectivity, and hydrometeor classifications from the polarimetric data over a nearly 4-h period are presented. It is shown that updraft velocity and vertical vorticity had to reach magnitudes of at least 10 m s−1 and 10−2 s−1 and occupy major portions of the storm before it could produce most of the observed severe storm characteristics. Furthermore, the establishment of cyclonic horizontal flow around the right flank of the updraft core was essential for hail production. Most of the largest hail grew from near millimeter-sized particles that originated in the mid- to upper-level stagnation region that resulted from obstacle-like flow of environmental air around the divergent outflow from the upper part of the updraft. These recycling embryonic particles descended around the right flank of the updraft core and reentered the updraft, intermingling with other smaller particles that had grown from cloud base along the main low-level updraft stream.

The primary objective of this study was to estimate the percentage of U.S. tornadoes that are spawned annually by squall lines and bow echoes, or quasi-linear convective systems (QLCSs). This was achieved by examining radar reflectivity images for every tornado event recorded during 1998–2000 in the contiguous United States. Based on these images, the type of storm associated with each tornado was classified as cell, QLCS, or other. Of the 3828 tornadoes in the database, 79% were produced by cells, 18% were produced by QLCSs, and the remaining 3% were produced by other storm types, primarily rainbands of landfallen tropical cyclones. Geographically, these percentages as well as those based on tornado days exhibited wide variations. For example, 50% of the tornado days in Indiana were associated with QLCSs. In an examination of other tornado attributes, statistically more weak (F1) and fewer strong (F2–F3) tornadoes were associated with QLCSs than with cells. QLCS tornadoes were more probable during the winter months than were cells. And finally, QLCS tornadoes displayed a comparatively higher and statistically significant tendency to occur during the late night/early morning hours. Further analysis revealed a disproportional decrease in F0–F1 events during this time of day, which led the authors to propose that many (perhaps as many as 12% of the total) weak QLCSs tornadoes were not reported.

This chapter focuses on the lightning and charge structure observed in supercells; a particular class of severe thunderstorms that are typically responsible for producing the most violent hailstorms and tornadoes. The chapter reviews common supercell features, the classification of supercells based on structural and visual characteristics, as well as the microphysics of supercells related to hail growth and electrification. A summary of observed flash rates, lightning polarity, and charge structure within supercells is also presented, along with descriptions of a few case studies as examples. Specific attention has been given to recent observations of anomalous positive cloud-to-ground lightning and inverted charge structures. Though this behavior is not exclusive to supercells, it has typically been observed in strong and/or severe storms that often meet supercell criteria.

Following up on an almost thirty year long history of International Cloud Modeling Workshops, that started out with a meeting in Irsee, Germany in 1985, the 8th International Cloud Modeling Workshop was held in July 2012 in Warsaw, Poland. The workshop, hosted by the Institute of Geophysics at the University of Warsaw, was organized by Szymon Malinowski and his local team of students and co-chaired by Wojciech Grabowski (NCAR/MMM) and Andreas Muhlbauer (University of Washington). International Cloud Modeling Workshops have been held traditionally every four years typically during the week before the International Conference on Clouds and Precipitation (ICCP) . Rooted in the World Meteorological Organization’s (WMO) weather modification program, the core objectives of the Cloud Modeling Workshop have been centered at the numerical modeling of clouds, cloud microphysics, and the interactions between cloud microphysics and cloud dynamics. In particular, the goal of the workshop is to provide insight into the pertinent problems of today’s state-of-the-art of cloud modeling and to identify key deficiencies in the microphysical representation of clouds in numerical models and cloud parameterizations. In recent years, the workshop has increasingly shifted the focus toward modeling the interactions between aerosols and clouds and provided case studies to investigatemore » both the effects of aerosols on clouds and precipitation as well as the impact of cloud and precipitation processes on aerosols. This time, about 60 (?) scientists from about 10 (?) different countries participated in the workshop and contributed with discussions, oral and poster presentations to the workshop’s plenary and breakout sessions. Several case leaders contributed to the workshop by setting up five observationally-based case studies covering a wide range of cloud types, namely, marine stratocumulus, mid-latitude squall lines, mid-latitude cirrus clouds, Arctic stratus and winter-time orographic clouds and precipitation. Interested readers are encouraged to visit the workshop website at http://www.atmos.washington.edu/~andreasm/workshop2012/ and browse through the list of case studies. The web page also provides a detailed list of participants and the workshop agenda. Aside from contributed oral and poster presentations during the workshop’s plenary sessions, parallel breakout sessions focused on presentations and discussions of the individual cases. A short summary and science highlights from each of the cases is presented below.« less