skip to main content

Title: Improving Simulations of Warm Rain in a Bulk Microphysics Scheme

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

Significance Statement

Although the accurate simulation of warm rain is critical to forecasting the hydrology of coastal areas and windward slopes, many warm rain parameterizations underpredict precipitation in these locations. This study introduces and evaluates modifications to the Thompson–Eidhammer microphysics parameterization scheme that significantly improve the accuracy of rainfall prediction in those regions.

more » « less
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Monthly Weather Review
Medium: X Size: p. 169-185
["p. 169-185"]
Sponsoring Org:
National Science Foundation
More Like this
  1. The OLYMPEX field campaign, which took place around the Olympic Mountains of Washington State during winter 2015/16, provided data for evaluating the simulated microphysics and precipitation over and near that barrier. Using OLYMPEX observations, this paper assesses precipitation and associated microphysics in the WRF-ARW model over the U.S. Pacific Northwest. Model precipitation from the University of Washington real-time WRF forecast system during the OLYMPEX field program (November 2015–February 2016) and an extended period (2008–18) showed persistent underprediction of precipitation, reaching 100 mm yr−1over the windward side of the coastal terrain. Increasing horizontal resolution does not substantially reduce this underprediction. Evaluating surface disdrometer observations during the 2015/16 OLYMPEX winter, it was found that the operational University of Washington WRF modeling system using Thompson microphysics poorly simulated the rain drop size distribution over a windward coastal valley. Although liquid water content was represented realistically, drop diameters were overpredicted, and, consequently, the rain drop distribution intercept parameter was underpredicted. During two heavy precipitation periods, WRF realistically simulated environmental conditions, including wind speed, thermodynamic structures, integrated moisture transport, and melting levels. Several microphysical parameterization schemes were tested in addition to the Thompson scheme, with each exhibiting similar biases for these two events. We show that the parameterization of aerosols over the coastal Northwest offered only minor improvement.

    more » « less
  2. Abstract

    Radar retrievals of drop size distribution (DSD) parameters are developed and evaluated over the mountainous Olympic Peninsula of Washington State. The observations used to develop retrievals were collected during the 2015/16 Olympic Mountain Experiment (OLYMPEX) and included the NASA S-band dual-polarimetric (NPOL) radar and a collection of second-generation Particle Size and Velocity (PARSIVEL2) disdrometers over the windward slopes of the barrier. Nonlinear and random forest regressions are applied to the PARSIVEL2 data to develop retrievals for median volume diameter, liquid water content, and rain rate. Improvement in DSD retrieval accuracy, defined by the mean error of the retrieval relative to PARSIVEL2 observations, was achieved when using the random forest model when compared with nonlinear regression. Evaluation of disdrometer observations and the retrievals from NPOL indicate that the radar retrievals can accurately reproduce observed DSDs in this region, including the common wintertime regime of small but numerous raindrops that is important there. NPOL retrievals during the OLYMPEX period are further evaluated using two-dimensional video disdrometers (2DVD) and vertically pointing Micro Rain Radars. Results indicate that radar retrievals using random forests may be skillful in capturing DSD characteristics in the lowest portions of the atmosphere.

    more » « less
  3. null (Ed.)
    The Remote sensing of Electrification, Lightning, And Meso-scale/micro-scale Processes with Adaptive Ground Observations (RELAMPAGO) and the Cloud, Aerosol, and Complex Terrain Interactions Experiment Proposal (CACTI) field campaigns provided an unprecedented thirteen-disdrometer dataset in Central Argentina during the Intensive (IOP, 15 November to 15 December 2018) and Extended (EOP, 15 October 2018 to 30 April 2019) Observational Periods. The drop size distribution (DSD) parameters and their variability were analyzed across the region of interest, which was divided into three subregions characterized by the differing proximity to the Sierras de Córdoba (SDC), in order to assess the impact of complex terrain on the DSD parameters. A rigorous quality control of the data was first performed. The frequency distributions of DSD-derived parameters were analyzed, including the normalized intercept parameter (logNw), the mean volume diameter (D0), the mean mass diameter (Dm), the shape parameter (μ), the liquid water content (LWC), and the rain rate (R). The region closest to the SDC presented higher values of logNw, lower D0, and higher μ, while the opposite occurred in the farthest region, i.e., the concentration of small drops decreased while the concentration of bigger drops increased with the distance to the east of the SDC. Furthermore, the region closest to the SDC showed a bimodal distribution of D0: the lower values of D0 were associated with higher values of logNw and were found more frequently during the afternoon, while the higher D0 were associated with lower logNw and occurred more frequently during the night. The data were analyzed in comparison to the statistical analysis of Dolan et al. 2018 and sorted according to the classification proposed in the cited study. The logNw-D0 and LWC-D0 two-dimensional distributions allowed further discussion around the applicability of other mid-latitude and global precipitation classification schemes (startiform/convection) in the region of interest. Finally, three precipitation case studies were analyzed with supporting polarimetric radar data in order to relate the DSD characteristics to the precipitation type and the microphysical processes involved in each case. 
    more » « less
  4. Abstract

    Rainfall prediction by weather forecasting models is strongly dependent on the microphysical parametrization being utilized within the model. As forecasting models have become more advanced, they are more commonly using double‐moment bulk microphysical parametrizations. While these double‐moment schemes are more sophisticated and require fewera prioriparameters than single‐moment parametrizations, a number of parameter values must still be fixed for quantities that are not Prognosed or diagnosed. Two such parameters, the width of the rain drop size distribution and the choice of collection efficiencies between liquid hydrometeors, are examined here. Simulations of deep convective storms were performed in which the collection efficiency dataset and thea prioriwidth of the rain drop size distribution (RSD) were individually and simultaneously modified. Analysis of the results show that thea prioriwidth of the RSD was a larger control on the total accumulated precipitation (a change of up to 75% over the typical values tested in this article) than the choice of collection efficiency dataset used (a change of up to 10%). Changing the collection efficiency dataset produces most of the impacts on precipitation rates through changes in the warm rain process rates. On the other hand, the decrease in precipitation with narrowing RSDs occurs in association with the following processes: (a) decreased rain production due to increased evaporation, (b) decreased rain production due to decreased ice melting, and (c) slower raindrop fall speed which leads to longer residency times and changes in rain self‐collection. These results add to the growing body of work showing that the representation of hydrometeor size distributions is critically important, and suggests that more work should be done to better represent the width of the RSD in models, including further development of triple‐moment and bin schemes.

    more » « less
  5. This study evaluates moist physics in the Weather Research and Forecasting (WRF) Model using observations collected during the Olympic Mountains Experiment (OLYMPEX) field campaign by the Global Precipitation Measurement (GPM) satellite, including data from the GPM Microwave Imager (GMI) and Dual-Frequency Precipitation Radar (DPR) instruments. Even though WRF using Thompson et al. microphysics was able to realistically simulate water vapor concentrations approaching the barrier, there was underprediction of cloud water content and rain rates offshore and over western slopes of terrain. We showed that underprediction of rain rate occurred when cloud water was underpredicted, establishing a connection between cloud water and rain-rate deficits. Evaluations of vertical hydrometeor mixing ratio profiles indicated that WRF produced too little cloud water and rainwater content, particularly below 2.5 km, with excessive snow above this altitude. Simulated mixing ratio profiles were less influenced by coastal proximity or midlatitude storm sector than were GMI profiles. Evaluations of different synoptic storm sectors suggested that postfrontal storm sectors were simulated most realistically, while warm sectors had the largest errors. DPR observations confirm the underprediction of rain rates noted by GMI, with no dependence on whether rain occurs over land or water. Finally, WRF underpredicted radar reflectivity below 2 km and overpredicted above 2 km, consistent with GMI vertical mixing ratio profiles.

    more » « less