A method for the parameterization of ice-phase microphysics is proposed and used to develop a new bulk microphysics scheme. All ice-phase particles are represented by several physical properties that ...evolve freely in time and space. The scheme prognoses four ice mixing ratio variables, total mass, rime mass, rime volume, and number, allowing 4 degrees of freedom for representing the particle properties using a single category. This approach represents a significant departure from traditional microphysics schemes in which ice-phase hydrometeors are partitioned into various predefined categories (e.g., cloud ice, snow, and graupel) with prescribed characteristics. The liquid-phase component of the new scheme uses a standard two-moment, two-category approach. The proposed method and a complete description of the new predicted particle properties (P3) scheme are provided. Results from idealized model simulations of a two-dimensional squall line are presented that illustrate overall behavior of the scheme. Despite its use of a single ice-phase category, the scheme simulates a realistically wide range of particle characteristics in different regions of the squall line, consistent with observed ice particles in real squall lines. Sensitivity tests show that both the prediction of the rime mass fraction and the rime density are important for the simulation of the squall-line structure and precipitation.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
A new two-moment stratiform cloud microphysics scheme in a general circulation model is described. Prognostic variables include cloud droplet and cloud ice mass mixing ratios and number ...concentrations. The scheme treats several microphysical processes, including hydrometeor collection, condensation/evaporation, freezing, melting, and sedimentation. The activation of droplets on aerosol is physically based and coupled to a subgrid vertical velocity. Unique aspects of the scheme, relative to existing two-moment schemes developed for general circulation models, are the diagnostic treatment of rain and snow number concentration and mixing ratio and the explicit treatment of subgrid cloud water variability for calculation of the microphysical process rates.
Numerical aspects of the scheme are described in detail using idealized one-dimensional offline tests of the microphysics. Sensitivity of the scheme to time step, vertical resolution, and numerical method for diagnostic precipitation is investigated over a range of conditions. It is found that, in general, two substeps are required for numerical stability and reasonably small time truncation errors using a time step of 20 min; however, substepping is only required for the precipitation microphysical processes rather than the entire scheme. A new numerical approach for the diagnostic rain and snow produces reasonable results compared to a benchmark simulation, especially at low vertical resolution. Part II of this study details results of the scheme in single-column and global simulations, including comparison with observations.
Celotno besedilo
Dostopno za:
BFBNIB, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Idealized simulations of the 15 May 2009 squall line from the Second Verification of the Origins of Rotation in Tornadoes Experiment (VORTEX2) are evaluated in this study. Four different ...microphysical setups are used, with either single-moment (1M) or double-moment (2M) microphysics, and either hail or graupel as the dense (rimed) ice species. Three different horizontal grid spacings are used: Delta x = 4, 1, or 0.25 km (with identical vertical grids). Overall, results show that simulated squall lines are sensitive to both microphysical setup and horizontal resolution, although some quantities (i.e., surface rainfall) are more sensitive to Delta x in this study. Simulations with larger Delta x are slower to develop, produce more precipitation, and have higher cloud tops, all of which are attributable to larger convective cells that do not entrain midlevel air. The highest-resolution simulations have substantially more cloud water evaporation, which is partly attributable to the development of resolved turbulence. For a given Delta x, the 1M simulations produce less rain, more intense cold pools, and do not have trailing stratiform precipitation at the surface, owing to excessive rainwater evaporation. The simulations with graupel as the dense ice species have unrealistically wide convective regions. Comparison against analyses from VORTEX2 data shows that the 2M setup with hail and Delta x = 0.25 km produces the most realistic simulation because (i) this simulation produces realistic distributions of reflectivity associated with convective, transition, and trailing stratiform regions, (ii) the cold pool properties are reasonably close to analyses from VORTEX2, and (iii) relative humidity in the cold pool is closest to observations.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
Idealized three-dimensional supercell simulations were performed using the two-moment bulk microphysics schemes of Morrison and Milbrandt–Yau in the Weather Research and Forecasting (WRF) model. ...Despite general similarities in these schemes, the simulations were found to produce distinct differences in storm structure, precipitation, and cold pool strength. In particular, the Morrison scheme produced much higher surface precipitation rates and a stronger cold pool, especially in the early stages of storm development. A series of sensitivity experiments was conducted to identify the primary differences between the two schemes that resulted in the large discrepancies in the simulations.
Different approaches in treating graupel and hail were found to be responsible for many of the key differences between the baseline simulations. The inclusion of hail in the baseline simulation using the Milbrant–Yau scheme with two rimed-ice categories (graupel and hail) had little impact, and therefore resulted in a much different storm than the baseline run with the single-category (hail) Morrison scheme. With graupel as the choice of the single rimed-ice category, the simulated storms had considerably more frozen condensate in the anvil region, a weaker cold pool, and reduced surface precipitation compared to the runs with only hail, whose higher terminal fall velocity inhibited lofting. The cold pool strength was also found to be sensitive to the parameterization of raindrop breakup, particularly for the Morrison scheme, because of the effects on the drop size distributions and the corresponding evaporative cooling rates. The use of a more aggressive implicit treatment of drop breakup in the baseline Morrison scheme, by limiting the mean–mass raindrop diameter to a maximum of 0.9 mm, opposed the tendency of this scheme to otherwise produce large mean drop sizes and a weaker cold pool compared to the hail-only run using the Milbrandt–Yau scheme.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
Hugh Morrison argues that children's support of Protestant missionary activity since the early 1800s has been an educational movement rather than a financial one and outlines how it has shaped minds ...and bodies for the sake of God, empire and nation.
A state‐of‐the‐art Lagrangian microphysics scheme is used in a large‐eddy simulation to investigate the stratocumulus transition from closed to open cell structure. Processes controlling ...precipitation development, which is a key to the transition, are analyzed by leveraging unique benefits of Lagrangian microphysics, particularly the ability to track computational drops in the flow. Sufficient time is needed for coalescence growth of cloud drops to drizzle within the updraft‐downdraft cycle of large eddies. This favors broad drop size distributions (DSDs) and drizzle growth in downdrafts, where drops are typically much older than in updrafts. During the closed cell stage, mean cloud drop radius is too small, and the DSDs are too narrow, so that the timescale for coalescence is much longer than the large eddy turnover time and drizzle growth is limited. The closed‐to‐open cell transition occurs when these timescales become comparable and the precipitation flux increases sharply.
Plain Language Summary
Low‐level stratocumulus clouds critically impact Earth's energy budget by reflecting incoming solar radiation. Accurate representation of these shallow clouds in climate models requires information on fine‐scale cloud properties and processes (e.g., cloud and drizzle drop sizes and their growth) due to significant variability driven by turbulent flow. Sufficient drizzle formation in these clouds changes their flow structure, shifting clouds from a uniform layer to a layer with an open cell structure. Using a state‐of‐the‐art particle‐based cloud model, this study showed that cloud droplet sizes and their coalescence growth change during this process in a way that sharply enhances precipitation. Once the clouds change to an open structure, cloud and drizzle drop sizes and growth differ in upward and downward moving branches of large turbulent eddies, with more drizzle formation and variation in drop sizes in downward flow regions. Information on the drop size variability and growth presented here could lead to better representation of clouds in climate models and measurements using ground and satellite‐based instruments.
Key Points
The first Lagrangian microphysics simulation of the stratocumulus transition from closed to open cell structure is documented
During the transition, the rain rate increases sharply as the coalescence timescale decreases relative to the large eddy turnover
Drop size distributions in open cell stratocumulus are broader in downdrafts than updrafts from coalescence, evaporation and drop mixing
Abstract
Simulations of a squall line observed on 20 May 2011 during the Midlatitude Continental Convective Clouds Experiment (MC3E) using 750- and 250-m horizontal grid spacing are performed. The ...higher-resolution simulation has less upshear-tilted deep convection and a more elevated rear inflow jet than the coarser-resolution simulation in better agreement with radar observations. A stronger cold pool eventually develops in the 250-m run; however, the more elevated rear inflow counteracts the cold pool circulation to produce more upright convective cores relative to the 750-m run. The differing structure in the 750-m run produces excessive midlevel front-to-rear detrainment, reinforcing excessive latent cooling and rear inflow descent at the rear of the stratiform region in a positive feedback. The contrasting mesoscale circulations are connected to early stage deep convective draft differences in the two simulations. Convective downdraft condensate mass, latent cooling, and downward motion all increase with downdraft area similarly in both simulations. However, the 750-m run has a relatively greater number of wide and fewer narrow downdrafts than the 250-m run averaged to the same 750-m grid, a consequence of downdrafts being under-resolved in the 750-m run. Under-resolved downdrafts in the 750-m run are associated with under-resolved updrafts and transport mid–upper-level zonal momentum downward to low levels too efficiently in the early stage deep convection. These results imply that under-resolved convective drafts in simulations may vertically transport air too efficiently and too far vertically, potentially biasing buoyancy and momentum distributions that impact mesoscale convective system evolution.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
The rapid glaciation of tropical cumulus clouds has been an enigma and has been debated in the literature for over 60 years. Possible mechanisms responsible for the rapid freezing have been ...postulated, but until now direct evidence has been lacking. Recent high-speed photography of electrostatically suspended supercooled drops in the laboratory has shown that freezing events produce small secondary ice particles. Aircraft observations from the Ice in Clouds Experiment-Tropical (ICE-T), strongly suggest that the drop-freezing secondary ice production mechanism is operating in strong, tropical cumulus updraft cores. The result is the production of small ice particles colliding with large supercooled drops (hundreds of microns up to millimeters in diameter), producing a cascading process that results in rapid glaciation of water drops in the updraft. The process was analyzed from data collected using state-of-the-art cloud particle probes during 54 Learjet penetrations of strong cumulus updraft cores over open ocean in a temperature range from 5 degree to -20 degree C. Repeated Learjet penetrations of an updraft core containing 3-5 g m super(-3) supercooled liquid showed an order-of-magnitude decrease in liquid mass concentration 3 min later at an elevation 1-1.5 km higher in the cloud. The aircraft observations were simulated using a one-dimensional cloud model with explicit bin microphysics. The model was initialized with drop and ice particle size distributions observed prior to rapid glaciation. Simulations show that the model can explain the observed rapid glaciation by the drop-freezing secondary ice production process and subsequent riming, which results when large supercooled drops collide with ice particles.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
9.
Modeling Condensation in Deep Convection Grabowski, Wojciech W.; Morrison, Hugh
Journal of the atmospheric sciences,
07/2017, Letnik:
74, Številka:
7
Journal Article
Recenzirano
Odprti dostop
Abstract
Cloud-scale models apply two drastically different methods to represent condensation of water vapor to form and grow cloud droplets. Maintenance of water saturation inside liquid clouds is ...assumed in the computationally efficient saturation adjustment approach used in most bulk microphysics schemes. When super- or subsaturations are allowed, condensation/evaporation can be calculated using the predicted saturation ratio and (either predicted or prescribed) mean droplet radius and concentration. The study investigates differences between simulations of deep unorganized convection applying a saturation adjustment condensation scheme (SADJ) and a scheme with supersaturation prediction (SPRE). A double-moment microphysics scheme with CCN activation parameterized as a function of the local vertical velocity is applied to compare cloud fields simulated applying SPRE and SADJ. Clean CCN conditions are assumed to demonstrate upper limits of the SPRE and SADJ difference. Microphysical piggybacking is used to extract the impacts with confidence. Results show a significant impact on deep convection dynamics, with SADJ featuring more cloud buoyancy and thus stronger updrafts. This leads to around a 3% increase of the surface rain accumulation in SADJ. Upper-tropospheric anvil cloud fractions are much larger in SPRE than in SADJ because of the higher ice concentrations and thus longer residence times of anvil particles in SPRE, as demonstrated by sensitivity tests. Higher ice concentrations in SPRE come from significantly larger ice supersaturations in strong convective updrafts that feature water supersaturations of several percent.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract
A novel bulk microphysics scheme that predicts the evolution of ice properties, including aspect ratio (shape), mass, number, size, and density is described, tested, and demonstrated. The ...scheme is named the Ice-Spheroids Habit Model with Aspect-Ratio Evolution (ISHMAEL). Ice is modeled as spheroids and is nucleated as one of two species depending on nucleation temperature. Microphysical process rates determine how shape and other ice properties evolve. A third aggregate species is also employed, diversifying ice properties in the model. Tests of ice shape evolution during vapor growth and riming are verified against wind tunnel data, revealing that the model captures habit-dependent riming and its effect on fall speed. Lagrangian parcel studies demonstrate that the bulk model captures ice property evolution during riming and melting compared with a bin model. Finally, the capabilities of ISHMAEL are shown in a 2D kinematic framework with a simple updraft. A direct result of predicting ice shape evolution is that various states of ice from unrimed to lightly rimed to densely rimed can be modeled without converting ice mass between predefined ice categories (e.g., snow and graupel). This leads to a different spatial precipitation distribution compared with the traditional method of separating snow and graupel and converting between the two categories, because ice in ISHMAEL sorts in physical space based on the amount of rime, which controls the thickness and therefore fall speed. Predicting these various states of rimed ice leads to a reduction in vapor growth rate and an increase in riming rate in a simple updraft compared with the traditional approach.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK