Uncertainty Quantification of WRF Model for Rainfall Prediction over the Sichuan Basin, China

Abstract

The mesoscale Weather Research and Forecasting (WRF) model has been widely employed to forecast day-ahead rainfalls. However, the deterministic predictions from the WRF model incorporate relatively large errors due to numerical discretization, inaccuracies in initial/boundary conditions and parameterizations, etc. Among them, the uncertainties in parameterization schemes have a huge impact on the forecasting skill of rainfalls, especially over the Sichuan Basin which is located east of the Tibetan Plateau in southwestern China. To figure out the impact of various parameterization schemes and their interactions on rainfall predictions over the Sichuan Basin, the Global Forecast System data are chosen as the initial/boundary conditions for the WRF model and 48 ensemble tests have been conducted based on different combinations of four microphysical (MP) schemes, four planetary boundary layer (PBL) schemes, and three cumulus (CU) schemes, for four rainfall cases in summer. Compared to the observations obtained from the Chinese ground-based and encrypted stations, it is found that the Goddard MP scheme together with the asymmetric convective model version 2 PBL scheme outperforms other combinations. Next, as the first step to explore further improvement of the WRF physical schemes, the polynomial chaos expansion (PCE) approach is then adopted to quantify the impacts of several empirical parameters with uncertainties in the WRF Single Moment 6-class (WSM6) MP scheme, the Yonsei University (YSU) PBL scheme and the Kain-Fritsch CU scheme on WRF rainfall predictions. The PCE statistics show that the uncertainty of the scaling factor applied to ice fall velocity in the WSM6 scheme and the profile shape exponent in the YSU scheme affects more dominantly the rainfall predictions in comparison with other parameters, which sheds a light on the importance of these schemes for the rainfall predictions over the Sichuan Basin and suggests the next step to further improve the physical schemes.