Sensitivity Study of Cloud-Resolving Convective Simulations with WRF Using Two Bulk Microphysical Parameterizations: Ice-Phase Microphysics versus Sedimentation Effects

2009 ◽  
Vol 48 (1) ◽  
pp. 61-76 ◽  
Author(s):  
Song-You Hong ◽  
Kyo-Sun Sunny Lim ◽  
Ju-Hye Kim ◽  
Jeong-Ock Jade Lim ◽  
Jimy Dudhia

Abstract This study examines the relative importance of ice-phase microphysics and sedimentation velocity for hydrometeors in bulk microphysics schemes. The two bulk microphysics schemes having the same number of prognostic water substances, the Weather Research and Forecasting (WRF) Single-Moment 6-Class Microphysics Scheme (WSM6) and the Purdue–Lin scheme (PLIN), are evaluated for a 2D idealized storm case and for a 3D heavy rainfall event over Korea. The relative importance of microphysics and sedimentation velocity for ice particles is illuminated by the additional experiments that exchange the sedimentation velocity formula for graupel in the two schemes. In a 2D idealized storm simulation test bed, it is found that, relative to the PLIN scheme, the WSM6 scheme develops the storm late with weakened intensity because of a slower sedimentation velocity for graupel. Such a weakened intensity of precipitation also appears in a 3D model framework when the WSM6 scheme is used, in conjunction with the overall distribution of the precipitation band southward toward what was observed. The major reason is found to be the ice-phase microphysics of the WSM6 and related ice-cloud–radiation feedback, rather than the smaller terminal velocity for graupel in the WSM6 than in the PLIN scheme.

2010 ◽  
Vol 138 (5) ◽  
pp. 1778-1791 ◽  
Author(s):  
Hann-Ming Henry Juang ◽  
Song-You Hong

Abstract A semi-Lagrangian advection scheme is developed for falling hydrometeors in hopes of replacing the conventional Eulerian scheme that has been widely used in the cloud microphysics scheme of numerical atmospheric models. This semi-Lagrangian scheme uses a forward advection method to determine the advection path with or without iteration, and advected mass in a two-time-level algorithm with mass conservation. Monotonicity is considered in mass-conserving interpolation between Lagrangian grids and model Eulerian grids, thus making it a positive definite advection scheme. For mass-conserving interpolation between the two grid systems, the piecewise constant method (PCM), piecewise linear method (PLM), and piecewise parabolic method (PPM) are proposed. The falling velocity at the bottom cell edge is modified to avoid unphysical deformation by scanning from the top layer to the bottom of the model, which enables the use of a large time step with reasonable accuracy. The scheme is implemented and tested in the Weather Research and Forecasting (WRF) Single-Moment 3-Class Microphysics Scheme (WSM3). In a theoretical test bed with constant terminal velocity, the proposed semi-Lagrangian algorithm shows that the higher-order interpolation scheme produces less diffusive features at maximal precipitation. Results from another idealized test bed with mass-weighted terminal velocity demonstrate that the accuracy of the proposed scheme is still satisfactory even with a time step of 120 s when the mean terminal velocity averaged at the departure and arrival points is employed. A two-dimensional (2D) squall-line test using the WSM3 scheme shows that the control run with the Eulerian advection scheme and the semi-Lagrangian run with the PCM method reveal similar results, whereas behaviors using the PLM and PPM are similar with higher-resolution features, such as mammatus-like clouds.


2010 ◽  
Vol 138 (11) ◽  
pp. 4254-4267 ◽  
Author(s):  
Andrew L. Molthan ◽  
Walter A. Petersen ◽  
Stephen W. Nesbitt ◽  
David Hudak

Abstract The Canadian CloudSat/Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) Validation Project (C3VP) was a field campaign designed to obtain aircraft, surface, and radar observations of clouds and precipitation in support of improving the simulation of snowfall and cold season precipitation, their microphysical processes represented within forecast models, and radiative properties relevant to remotely sensed retrievals. During the campaign, a midlatitude cyclone tracked along the U.S.–Canadian border on 22 January 2007, producing an extensive area of snowfall. Observations of ice crystals from this event are used to evaluate the assumptions and physical relationships for the snow category within the Goddard six-class, single-moment microphysics scheme, as implemented within the Weather Research and Forecasting (WRF) model. The WRF model forecast generally reproduced the precipitation and cloud structures sampled by radars and aircraft, permitting a comparison between C3VP observations and model snowfall characteristics. Key snowfall assumptions in the Goddard scheme are an exponential size distribution with fixed intercept and effective bulk density, and the relationship between crystal diameter and terminal velocity. Fixed values for the size distribution intercept and density did not represent the vertical variability of naturally occurring populations of aggregates, and the current diameter and fall speed relationship underestimated terminal velocities for all sizes of crystals.


Author(s):  
Jason A. Milbrandt ◽  
Hugh Morrison ◽  
Daniel T. Dawson ◽  
Marco Paukert

AbstractIn the original Predicted Particle Properties (P3) bulk microphysics scheme, all ice-phase hydrometeors are represented by one or more “free” ice categories, where the physical properties evolve smoothly through changes to the four prognostic variables (per category,) and with a 2-moment representation of the particle size distribution. As such, the spectral dispersion cannot evolve independently, which thus results in limitations in representation of ice – in particular hail – due to necessary constraints in the scheme to prevent excessive gravitational size sorting. To overcome this, P3 has been modified to include a 3-moment representation of the size distribution of each ice category through the addition of a fifth prognostic variable, the sixth moment of the size distribution.The details of the 3-moment ice parameterization in P3 are provided. The behavior of the modified scheme, with the single-ice-category configuration, is illustrated through simulations in a simple 1D kinematic model framework as well as with near large-eddy-resolving (250-m grid spacing) 3D simulations of a hail-producing supercell. It is shown that the 3-moment ice configuration controls size sorting in a physically-based way and leads to an improved capacity to simulate large, heavily-rimed ice (hail), including mean and maximum sizes and reflectivity, and thus an overall improvement in the representation of ice-phase particles in the P3 scheme.


2005 ◽  
Vol 62 (7) ◽  
pp. 2118-2135 ◽  
Author(s):  
Vlado Spiridonov ◽  
Mladjen Curic

Abstract The relative importance of various processes to sulfate production and wet deposition is examined by using a cloud-resolving model coupled with a sulfate chemistry submodel. Results using different versions of the model are then compared and principal differences with respect to their dynamics, microphysics, and chemistry are carefully discussed. The results imply that the dominant microphysical and chemical conversions of sulfate in the 3D run are nucleation, scavenging, and oxidation. Due to the lower cloud water and rainwater pH, oxidation does not contribute as significantly to the sulfate mass in the 2D run as the 3D. Sensitivity tests have revealed that in-cloud scavenging in the 2D run for continental nonpolluted and continental polluted clouds accounted for 29.4% and 31.5% of the total sulfur deposited, respectively. The 3D run shows a lower percentage contribution to sulfur deposition for about 28.2% and 29.6%. In addition, subcloud scavenging for the 2D run contributed about 32.7% and 38.2%. In-cloud oxidation in the 2D run accounted for about 24.5% to 30.4% of the total sulfur mass deposited. Subcloud oxidation contributed from 21.0% to 20.6% of the total sulfur mass removed by wet deposition. In-cloud oxidation for the 3D run shows slightly lower percentage values when compared to those from the 2D run. The relative contribution of subcloud oxidation for continental nonpolluted and polluted clouds exceeds those values in the 2D run by approximately 7% and 10%, respectively. Ignoring the ice phase and considering those types of convective clouds in the 2D run may lead to a higher value of the total sulfur mass removed by the wet deposition of about 33.9% to 39.2% for the continental nonpolluted and 36.2% to 45.6% for the continental polluted distributions relative to the base runs.


2007 ◽  
Vol 64 (4) ◽  
pp. 1068-1088 ◽  
Author(s):  
Andrew J. Heymsfield ◽  
Gerd-Jan van Zadelhoff ◽  
David P. Donovan ◽  
Frederic Fabry ◽  
Robin J. Hogan ◽  
...  

Abstract This two-part study addresses the development of reliable estimates of the mass and fall speed of single ice particles and ensembles. Part I of the study reports temperature-dependent coefficients for the mass-dimensional relationship, m = aDb, where D is particle maximum dimension. The fall velocity relationship, Vt = ADB, is developed from observations in synoptic and low-latitude, convectively generated, ice cloud layers, sampled over a wide range of temperatures using an assumed range for the exponent b. Values for a, A, and B were found that were consistent with the measured particle size distributions (PSD) and the ice water content (IWC). To refine the estimates of coefficients a and b to fit both lower and higher moments of the PSD and the associated values for A and B, Part II uses the PSD from Part I plus coincident, vertically pointing Doppler radar returns. The observations and derived coefficients are used to evaluate earlier, single-moment, bulk ice microphysical parameterization schemes as well as to develop improved, statistically based, microphysical relationships. They may be used in cloud and climate models, and to retrieve cloud properties from ground-based Doppler radar and spaceborne, conventional radar returns.


2022 ◽  
Vol 12 (3) ◽  
pp. 29-43
Author(s):  
Samarendra Karmakar ◽  
Mohan Kumar Das ◽  
Md Quamrul Hassam ◽  
Md Abdul Mannan

The diagnostic and prognostic studies of thunderstorms/squalls are very important to save live and loss of properties. The present study aims at diagnose the different tropospheric parameters, instability and synoptic conditions associated the severe thunderstorms with squalls, which occurred at different places in Bangladesh on 31 March 2019. For prognostic purposes, the severe thunderstorms occurred on 31 March 2019 have been numerically simulated. In this regard, the Weather Research and Forecasting (WRF) model is used to predict different atmospheric conditions associated with the severe storms. The study domain is selected for 9 km horizontal resolution, which almost covers the south Asian region. Numerical experiments have been conducted with the combination of WRF single-moment 6 class (WSM6) microphysics scheme with Yonsei University (YSU) PBL scheme in simulation of the squall events. Model simulated results are compared with the available observations. The observed values of CAPE at Kolkata both at 0000 and 1200 UTC were 2680.4 and 3039.9 J kg-1 respectively on 31 March 2019 and are found to be comparable with the simulated values. The area averaged actual rainfall for 24 hrs is found is 22.4 mm, which complies with the simulated rainfall of 20-25 mm for 24 hrs. Journal of Engineering Science 12(3), 2021, 29-43


2021 ◽  
Author(s):  
Antonio Ricchi ◽  
Vincenzo Mazzarella ◽  
Lorenzo Sangelantoni ◽  
Gianluca Redaelli ◽  
Rossella Ferretti

<div> <p><span>A severe weather events hit Italy on July 9-10, 2019 causing heavy damages by the falling of large-size hail. A trough from Northern Europe affected Italy and the Balkans advecting cold air on the Adriatic Sea. The intrusion of relatively cold and dry air on the Adriatic Sea, in a first stage through the "Bora jets" generated by the Dinaric Alps gave rise to a frontal structure on the ground, which rapidly moved from North to South Adriatic. The large thermal gradient (also with the sea surface), the interaction with the complex orography and the coastal zone, generated several storm structures along the eastern Italian coast. In particular, on 10 July 2019 between 8UTC and 12UTC a deep convective cell (probably a supercell) developed along the coast North of the city of Pescara, producing intense rainfall (accumulated rainfall reaching 130 mm/3h) and a violent hailstorm with hailstones larger than 10 cm in diameter. The storm quickly moved southward, evolving into a complex multicellular structure clearly visible by observing radar data. In this work the frontal dynamics and the genesis of the storm cell are investigated using the numerical model WRF (Weather Research and Forecasting system). Numerical experiments are carried out using a 1 km grid on Central Italy, initialized using the ECMWF dataset and the Sea Surface Temperature (SST) taken by MFS-CMEMS Copernicus dataset. The sensitivity study investigated both the impact of the initial conditions, the quality and the anomaly of the SST on the Adriatic basin in those days. Furthermore, in order to quantify the importance of the use of different microphysics, Planetary boundary Layer (PBL) and radiative schemes, several experiments are performed. The role of orography in the development and location of the convective cell is also investigated. Preliminary results show that initialization and SST played a fundamental role. In particular, the initialization several hours before the event, coupled with a detailed SST allows to correctly reproduce the atmospheric fields. The microphysics scheme turned out to play a key role for this event by showing a significant greater impact than the PBL, in terms of frontal genesis on both the synoptic and local scale. </span></p> </div>


2005 ◽  
Vol 44 (4) ◽  
pp. 445-466 ◽  
Author(s):  
Jerry M. Straka ◽  
Edward R. Mansell

Abstract A single-moment bulk microphysics scheme with multiple ice precipitation categories is described. It has 2 liquid hydrometeor categories (cloud droplets and rain) and 10 ice categories that are characterized by habit, size, and density—two ice crystal habits (column and plate), rimed cloud ice, snow (ice crystal aggregates), three categories of graupel with different densities and intercepts, frozen drops, small hail, and large hail. The concept of riming history is implemented for conversions among the graupel and frozen drops categories. The multiple precipitation ice categories allow a range of particle densities and fall velocities for simulating a variety of convective storms with minimal parameter tuning. The scheme is applied to two cases—an idealized continental multicell storm that demonstrates the ice precipitation process, and a small Florida maritime storm in which the warm rain process is important.


Sign in / Sign up

Export Citation Format

Share Document