Numerical evaluation of WRF microphysics parameterizations for a major rainfall event over an arid region

In this study, the Weather Research and Forecasting (WRF) model version 3.2 is used to examine the impact of precipitating ice and especially snow-graupel partitioning in the simulation of a heavy rainfall event over Chalkidiki peninsula in Northern Greece. This major precipitation event, associated with a case of cyclogenesis over the Aegean Sea, occurred on the 8th of October 2006 causing severe flooding and damage. Two widely used microphysical parameterizations, the Purdue Lin (PLIN) and WRF Single-Moment 6-class scheme (WSM6) are compared with available raingauge measurements over the complex topography of Chalkidiki. To further investigate the importance of snow and graupel relative mass content and the treatment of precipitating ice sedimentation velocity, two older versions of the WSM6 scheme were compiled and run with the current model. The verification results indicate that all simulations were found to match raingauge data more closely over the eastern mountainous Chalkidiki peninsula where maximum accumulations were observed. In other stations all schemes overestimate 24h accumulated rainfall except a station situated at the western part of the peninsula, where none of the simulations was able to reproduce observed rainfall. Graupel dominance in PLIN generates rapid precipitation fallout at the point of maximum predicted 24h accumulation. Similar behavior is shown in WSM6 from WRF version 2, but with significant less rainfall. Increasing snow amounts aloft, due to the unified treatment of precipitating ice in WSM6 from WRF version 3, modifies rain dynamics which decrease rainfall rates, but increases 24h accumulations. A sensitivity experiment where PLIN is used with snow accretion by graupel turned off, indicated that this process seems to be the most important factor controlling the differences in surface precipitation between PLIN and WSM6 from WRF version 3, determining the spatial and temporal distribution of this heavy precipitation event. The results also revealed that snow overestimation can lead to high rainfall accumulations, even though rain is more evenly distributed over the 24h period, deteriorating precipitation forecast.

A new method for driving a One-Dimensional Stratiform Cold (1DSC) cloud model with Weather Research and Forecasting (WRF) model outputs was developed by conducting numerical experiments for a typical large-scale stratiform rainfall event that took place on 4–5 July 2004 in Changchun, China. Sensitivity test results suggested that, with hydrometeor profiles extracted from the WRF outputs as the initial input, and with continuous updating of soundings and vertical velocities (including downdraft) derived from the WRF model, the new WRF-driven 1DSC modeling system (WRF-1DSC) was able to successfully reproduce both the generation and dissipation processes of the precipitation event. The simulated rainfall intensity showed a time-lag behind that observed, which could have been caused by simulation errors of soundings, vertical velocities and hydrometeor profiles in the WRF output. Taking into consideration the simulated and observed movement path of the precipitation system, a nearby grid point was found to possess more accurate environmental fields in terms of their similarity to those observed in Changchun Station. Using profiles from this nearby grid point, WRF-1DSC was able to reproduce a realistic precipitation pattern. This study demonstrates that 1D cloud-seeding models do indeed have the potential to predict realistic precipitation patterns when properly driven by accurate atmospheric profiles derived from a regional shortrange forecasting system. This opens a novel and important approach to developing an ensemble-based rain enhancement prediction and operation system under a probabilistic framework concept.

Abstract. Regions with high ice water content (HIWC), composed of mainly small ice crystals, frequently occur over convective clouds in the tropics. Such regions can have median mass diameters (MMDs) <300 µm and equivalent radar reflectivities <20 dBZ. To explore formation mechanisms for these HIWCs, high-resolution simulations of tropical convective clouds observed on 26 May 2015 during the High Altitude Ice Crystals – High Ice Water Content (HAIC-HIWC) international field campaign based out of Cayenne, French Guiana, are conducted using the Weather Research and Forecasting (WRF) model with four different bulk microphysics schemes: the WRF single‐moment 6‐class microphysics scheme (WSM6), the Morrison scheme, and the Predicted Particle Properties (P3) scheme with one- and two-ice options. The simulations are evaluated against data from airborne radar and multiple cloud microphysics probes installed on the French Falcon 20 and Canadian National Research Council (NRC) Convair 580 sampling clouds at different heights. WRF simulations with different microphysics schemes generally reproduce the vertical profiles of temperature, dew-point temperature, and winds during this event compared with radiosonde data, and the coverage and evolution of this tropical convective system compared to satellite retrievals. All of the simulations overestimate the intensity and spatial extent of radar reflectivity by over 30 % above the melting layer compared to the airborne X-band radar reflectivity data. They also miss the peak of the observed ice number distribution function for 0.1<Dmax<1 mm. Even though the P3 scheme has a very different approach representing ice, it does not produce greatly different total condensed water content or better comparison to other observations in this tropical convective system. Mixed-phase microphysical processes at −10 ∘C are associated with the overprediction of liquid water content in the simulations with the Morrison and P3 schemes. The ice water content at −10 ∘C increases mainly due to the collection of liquid water by ice particles, which does not increase ice particle number but increases the mass/size of ice particles and contributes to greater simulated radar reflectivity.

Long‐term variations in the frequencies of different rainfall events during summer monsoon season (June through September) are analysed over a period of 38‐years (1980–2017) using the simulations of weather research and forecasting (WRF) model evaluated against the India Meteorological Department (IMD). Further, we have divided the study period into two sub‐periods such as before (1980–1998), and after the mega ENSO event (1999–2017). Our analysis suggests the WRF exhibits good skill in capturing the increased trends in the frequency of heavy and very heavy (extreme) rainfall events during the recent epoch (1999–2017).We have assessed the skill of WRF in capturing the observed probability distributions of daily rain rates over different sub‐regions of India. Our results depict that the downscaled products of WRF captured well the extreme rainfall as observed in IMD, signifying that this dynamically downscaled model with continuous 36‐hr re‐initialization makes it suitable to simulate extreme rainfall events compared to light rainfall event. However, our analysis reports that there is a slight overestimation in WRF rainfall and fails in capturing the light rainfall events compared to IMD. The overestimation of rainfall by WRF possibly be due to increased pressure gradients and enhances the convection in the upper streams and associated rainfall. Our study emphasizes that the variability in summer monsoon rainfall during the recent period is due to the changes in the monsoon circulations, particularly monsoon low‐level jet that leads to alter the moisture budget in terms of moisture convergence and transport. Moreover, the variability in frequencies of rainfall events are analysed for strong La Niña/El Niño years and observed an increase in extremes over major part of India during strong La Niña except over parts of Central India and NEI where strong El Niño dominates; which is well captured by the WRF model.

A two‐moment bulk microphysical scheme has been implemented into the Weather Research and Forecasting (WRF) model to investigate the aerosol‐cloud interaction. The microphysical scheme calculates the mass mixing ratios and number concentrations of aerosols and five types of hydrometeors and accounts for various cloud processes including warm and mixed phase microphysics. The representation of the aerosol size distribution is evaluated, showing that the three‐moment modal method produces results better in agreement with the sectional approach than the two‐moment modal method for variable supersaturation conditions in clouds. The effects of aerosols on cloud processes are investigated using the two‐moment bulk microphysical scheme in a convective cumulus cloud event occurring on 24 August 2000 in Houston, Texas. The modeled evolution of the distribution of radar reflectivity in the y‐z section, the cell lifetime, and averaged accumulated precipitation with the aerosol concentration under the polluted urban condition are qualitatively consistent with the measurements. Sensitivity simulations are initialized using a set of aerosol profiles with the number concentrations ranging from 200 to 50,000 cm−3 and mass ranging from 1 to 10 μg m−3 at the surface level. The response of precipitation to the increase of aerosol concentrations is nonmonotonic, because of the complicated interaction between cloud microphysics and dynamics. The precipitation increases with aerosol concentrations from clean maritime to continental background conditions, but is considerably reduced and completely suppressed under highly polluted conditions, indicating that the aerosol concentration exhibits distinct effects on the precipitation efficiency under different aerosol conditions. The maximal cloud cover, core updraft, and maximal vertical velocity exhibit similar responses as precipitation. Comparison is made to evaluate the effects of different autoconversion parameterizations and bulk microphysical schemes on cloud properties. Because of its broad application in numerical weather prediction, implementation of the two‐moment microphysical scheme to the WRF model will greatly facilitate assessment of aerosol‐cloud interaction from individual cumulus to mesoscale convective systems.

This study evaluates the performance of several cloud microphysics parameterizations in simulating surface precipitation for two snowstorm cases during the International Collaborative Experiment held at the PyeongChang 2018 Olympics and Winter Paralympic Games (ICE-POP 2018) field campaign. We compared four different schemes in the Weather Research and Forecasting (WRF) Model, namely the double-moment 6-class (WDM6), the WRF single-moment 6-class (WSM6), and Thompson and Morrison parameterizations. Both WSM6 and WDM6 overestimated the precipitation amount for the shallow precipitation system because of the substantial amount of cloud ice, mostly generated by the deposition process. The simulated precipitation amount and distribution for the deep precipitation system showed no noticeable differences in the different cloud microphysics parameterizations. However, the simulated hydrometeor type at the surface using WSM6 and WDM6 showed good agreement with observations for all cases. The accuracy of the mean mass-weighted terminal velocity of cloud ice applied in WSM6 and WDM6 is ±20%. The number concentration of cloud ice and the ice microphysics processes are newly retrieved with 1.2 times increased . For the shallow snowstorm, the precipitation amount was reduced by approximately 8% because of the inefficient deposition and its effects on the subsequent ice microphysical processes, such as the accretion of cloud ice by snow and the conversion from cloud ice to snow.

This study uses the Weather Research and Forecasting (WRF) Model to investigate the performance of hail parameterizations of the WRF double-moment 7-class (WDM7), aerosol-aware Thompson (AAT), and National Taiwan University triple-moment (NTU3M) bulk microphysics schemes (BMSs) on a real case of a hailstorm initiated in Shandong Province, China. The maximum hail size is particularly evaluated because it is crucial to hail severity prediction, along with areal coverage and intensity of the 24-h solid precipitation during the simulations. Compared with the radar-derived maximum hail size, the objective analysis shows that the NTU3M scheme has the best score in the forecast skill of hail-fall coverage and size, while two BMSs with single-moment rimed ice species overestimate hail diameters aloft but underpredict the coverage at the surface. A deeper investigation suggests that the derived size tendencies from the three BMSs are comparable to the benchmark solutions from the detailed hailstone growth and melting models. The NTU3M scheme displays the most consistent size tendency of the maximum diameter with the benchmark solution in the growth processes. The behaviors of melted diameter by parameterizations are highly related to the treatments of number concentration, which are consistent with the predicted hail severity and coverage. Finally, the sensitivity study shows that increasing the model resolution does not improve the forecast of the maximum hail size, given the biases in the hail mass budget equations and the parameterization of particle size distribution, with single-moment rimed ice species of the AAT scheme. Significance Statement Improving hail-forecasting skill, including the size, severity, and the spatial and temporal coverage of hail fall, has become an important subject for numerical weather prediction models as the model resolution increases. The objective of this study is to investigate the fundamental differences in hail parameterizations of three bulk microphysics schemes that lead to differences in the prediction of severe hail events and the spatial coverage of hail fall, hopefully providing insights into hail prediction with a regional numerical weather prediction model in the future.

This study is intended to simulate and better understand mean and extreme precipitation over the Loess Plateau (LP) in China using the Weather Research and Forecasting (WRF) model. We performed a long-term (1980–2011) simulation with the WRF model at a 10 km horizontal resolution, forced by the ERA-Interim reanalysis. A series of sensitivity experiments were conducted to investigate the influence of resolution and physical parameterization schemes on the simulation ability of the WRF model over the LP. The region has a dominant semi-arid climate and is adversely affected by extreme precipitation. Results show that the WRF model produced better simulations at a 10 km resolution for the LP with complex terrain than it did at coarser resolutions. WRF simulations for precipitation over the LP were most sensitive to cumulus schemes, moderately sensitive to planetary boundary layer schemes, and less sensitive to microphysics schemes. The WRF model not only adequately captured the spatial distribution of precipitation in the LP but also reproduced its magnitude and variability at different time scales. Moreover, the WRF model captured 73% of the observed extreme precipitation events with an intensity of more than 12 mm day− 1, which cause substantial soil erosion in the LP. Although it was unable to capture 27% of the observed events, the WRF model reasonably reproduced ERA-Interim extreme precipitation and its dynamic processes. These analyses strongly suggest that the biased extreme precipitation generated with the WRF model most likely resulted from inaccurate ERA-Interim forcing data. This study provides a better understanding of extreme precipitation simulations as well as insight into predictions of soil erosion in the LP.

Numerical weather prediction (NWP) models such as the Weather Research and Forecasting (WRF) model are increasingly used over the Indian region to forecast extreme rainfall events. However, studies which explore the application of high-resolution rainfall simulations obtained from the WRF model in urban hydrology are limited. In this paper, the utility of a model coupling framework to predict urban floods is explored through the case study of Bangalore city in India. This framework is used to simulate multiple extreme events that occurred over the city for the monsoons of years 2020 and 2021. To address the uncertainty from the WRF model, a 12-member convection permitting ensemble is used. Model configurations using Kain Fritsch and WSM6 parameterization schemes could simulate the spatial and temporal pattern of the selected event. The city is easily flooded with rainfall events above a threshold of 60 mm/day and to capture the response of the urban catchment, the Personal Computer Storm Water Management Model (PCSWMM) is used in this study. Flood forecasts are created using the outputs from the WRF ensemble and the Global Forecasting System (GFS). The high temporal and spatial resolution of the rainfall forecasts (<4 km at 15-min intervals), has proved critical in reproducing the urban flood event. The flood forecasts created using the WRF ensemble indicate that flooding and water levels are comparable to the observed whereas the GFS underestimates these to a large extent. Thus, the coupled WRF–PCSWMM modelling framework is found effective in forecasting flood events over an Indian city.

Atmospheric wet delay caused by Precipitable Water Vapor (PWV) significantly impacts the performance of many geodetic surveying systems such as Global Navigation Satellite System (GNSS). In this study, we use wet delay corrections forecast by the Weather Research and Forecasting (WRF) model to enhance GNSS Precise Point Positioning (PPP) during two observation periods with two different weather conditions, that is, period 1: March 01 to 14, 2020 (average PWV: 23.5 kg/m2) and period 2: June 02 to 15, 2020 (flooding weather with average PWV: 55.6 kg/m2), over the South China. PWV data from 277 to 263 GNSS stations are assimilated into WRF model to enhance the WRF water vapor forecasting capability for period 1 and period 2, respectively. Wet delay corrections from two different WRF configurations, that is, WRF no data assimilation and WRF with assimilation of GNSS PWV, are used to augment the PPP. Totally, eight WRF‐enhanced PPP schemes are tested. The results show that WRF‐enhanced PPP schemes generally have a better positioning performance in the up component than traditional PPP. After using WRF wet delay corrections, for static mode, the vertical positioning accuracy is improved by 14.6% and 33.7% for period 1 and period 2, respectively. The corresponding convergence time are reduced by 41.8% and 25.0% for period 1 and period 2, respectively. For kinematic mode, the positioning accuracy improvements in the up component reach 13.8% and 19.0% for period 1 and period 2, respectively. The kinematic PPP convergence time is reduced by up to 8.2% for period 1.

Impact of different reanalysis data on WRF dynamical downscaling over China

In this study, we conduct a 15‐year long (2008–2022) convection‐permitting simulation over the Tibetan Plateau (TP) using the Weather Research and Forecasting (WRF) model to evaluate its added value compared to ERA5 in characterizing subdaily precipitation and its extremes over the TP. Our findings indicate that the overestimation of precipitation amounts caused by overestimated precipitation frequency in ERA5 is mitigated in WRF mainly through a reduction in light precipitation occurrences. This improvement is likely linked to increased temperatures and decreased relative humidity in the lower troposphere. Moreover, the much earlier precipitation peak in ERA5 is corrected in WRF, a result attributed to differences in related net water vapor transport onto the TP. Additionally, WRF provides a more accurate estimation of the percentage of short‐ and long‐duration precipitation events with only slight underestimations and overestimations, respectively, whereas ERA5 exhibits significantly larger errors. The relatively larger temperature changes in WRF, indicating increased solar heating, are likely responsible for the higher occurrence of short‐duration precipitation events. In contrast, the anticyclonic anomalies and northeastern wind anomalies in WRF reduce the persistence of long‐duration precipitation. Importantly, WRF also produces more realistic timings of precipitation with different durations compared to ERA5. For extreme precipitation, short‐duration extreme events dominate accounting for 50% of extreme precipitation. These events are better reproduced by WRF than by ERA5. Although WRF outperforms ERA5 in estimating the proportion and contribution of extreme events to total precipitation and in capturing their triggering timing, a slight underestimation persists in the simulated proportion and contribution.

The fire weather of south-east Australia from 1985 to 2009 has been simulated using the Weather Research and Forecasting (WRF) model. The US National Oceanic and Atmospheric Administration Centers for Environmental Prediction and National Center for Atmospheric Research reanalysis supplied the lateral boundary conditions and initial conditions. The model simulated climate and the reanalysis were evaluated against station-based observations of the McArthur Forest Fire Danger Index (FFDI) using probability density function skill scores, annual cumulative FFDI and days per year with FFDI above 50. WRF simulated the main features of the FFDI distribution and its spatial variation, with an overall positive bias. Errors in average FFDI were caused mostly by errors in the ability of WRF to simulate relative humidity. In contrast, errors in extreme FFDI values were driven mainly by WRF errors in wind speed simulation. However, in both cases the quality of the observed data is difficult to ascertain. WRF run with 50-km grid spacing did not consistently improve upon the reanalysis statistics. Decreasing the grid spacing to 10 km led to fire weather that was generally closer to observations than the reanalysis across the full range of evaluation metrics used here. This suggests it is a very useful tool for modelling fire weather over the entire landscape of south-east Australia.

Artificial intelligence (AI) models have demonstrated advancements in computational efficiency and forecast accuracy relative to the Numerical Weather Prediction (NWP), but they are unable to fully represent high-dimensional atmospheric dynamics. Thus, some AI-NWP coupled frameworks have been proposed, such as integrating AI-driven boundary conditions with numerical models to leverage the strengths of both approaches. However, in this coupled framework, ensemble forecasts and associated error propagation and energy dynamics remain under-explored. In this study, an AI-NWP coupled system that also uses the stochastic kinetic energy backscatter scheme (SKEBS) to generate ensemble forecasts is established. Ensemble simulations of Typhoon Yutu (2018) are carried out with the Weather Research and Forecasting (WRF) model employing Pangu-Weather and FuXi forecast data as boundary forcing. The results show that the ensemble WRF_Pangu (WRF_FuXi) improved Yutu’s track forecast by 67% (50%) compared to the traditional physics-based WRF_GFS (Global Forecast System), and reduced its intensity underestimation by about 67% relative to their AI global counterparts. Nonetheless, WRF_FuXi and WRF_Pangu exhibited limited ensemble spread and linear error growth, reflecting deterministic tendencies. Comparison of global and regional experiments show that Pangu-Weather is more physically constrained and thus better aligned with the WRF model for regional applications, while the adaptation of FuXi to the regional model is less robust. Spectral analysis revealed that AI-derived boundaries introduced excessive small-scale energy and underestimated larger-scale energy. The regional model WRF acted as a “conveyor belt”, propagating additive small-scale energy upscale, ultimately overwhelming the stochastic perturbations for ensemble generation. These findings underscore the need to incorporate more physical features into the AI-derived boundary conditions for ensemble forecasting.

Assimilation of accurate external satellite-based water vapor data has been demonstrated to be able to improve the Numerical Weather Prediction (NWP) model’s forecasting performance. In this study, two types of near-infrared Precipitable Water Vapor (PWV) data from the Ocean and Land Colour Instrument (OLCI) of the Sentinel-3A and -3B satellites are assimilated into the Weather Research and Forecasting (WRF) model to enhance the weather forecasting performance over the South China area. One type is the raw PWV products obtained directly from the Sentinle-3A and -3B satellites; the other type is the Sentinel-3A and -3B PWV calibrated by PWV from Global Navigation Satellite System (GNSS). The Back Propagation Neural Network (BPNN) model is used in the PWV calibration. Correspondingly, two WRF assimilation schemes are adopted to assimilate the Sentinel-3 raw PWV and Sentinel-3 calibrated PWV. Additionally, the WRF without data assimilation scheme is used as the background run to investigate the impact of data assimilation. Three WRF experimental schemes are performed for two periods, i.e., February 2020 (dry period) and July 2020 (wet period). The forecasting performance of these three WRF schemes are evaluated by observations from GNSS, radiosonde, and surface meteorological stations. The comparison results with GNSS PWV show that assimilation of Sentienl-3 calibrated PWV can improve PWV forecasting accuracy by 5.6% for February period for the first 6 h after data assimilation. Evaluated by radiosonde profiles, the assimilation of Sentinel-3 both raw and calibrated PWV can improve WRF humidity and temperature forecasting accuracy over the July period. The accuracy improvement can be up to 9.2% (at 2 km to 3 km of altitude) for water vapor mixing ratio and 6.0% for temperature (at 4 km to 5 km of altitude). The WRF rainfall evaluation result has been validated by rainfall observations from meteorological stations and it shows that the assimilation of Sentinel-3 PWV can improve the rainfall forecasting performance over the July period. The rainfall forecasting success rates of three schemes of assimilating no data, Sentinel-3 raw PWV, and Sentinel-3 calibrated PWV are 69.4%, 69.6%, and 71.1%, respectively, over the July 2020 in the South China region.