Microphysics Research Articles

Optimal Configuration of Physical Process Parameterization Scheme Combination for Simulating Meteorological Variables in Weather Research and Forecasting Model: Based on Orthogonal Experimental Design and Comprehensive Evaluation Method

The weather research and forecasting (WRF) model is frequently used to investigate the meteorological field around nuclear installations. The configuration of physical process parameterization schemes in the WRF model has a significant impact on the accuracy of the simulation results. Consequently, carrying out a pre-experiment to quickly obtain the optimal combination of parameterization schemes is essential before conducting meteorological parameter research. To obtain the optimal combination of physical process parameterization schemes from the planetary boundary layer (PBL), land surface (LSF), microphysical (MP), long-wave (LW), and short-wave (SW) radiation processes of the WRF model for simulating the near-surface meteorological variables near a nuclear power plant in Sanshan Town, Fuqing City, Fujian Province, China on 4 June 2019 were observed. Orthogonal experimental design (OED), a comprehensive evaluation method based on the CRiteria Import Through Intercriteria Correlation (CRITIC) weight analysis, and comprehensive balance method were employed for the first time to conduct the research. The sensitivity of meteorological variables to physical processes was first discussed. The findings revealed that the PBL scheme configuration had a profound impact on simulating wind fields. Furthermore, the LSF scheme configuration had a significant influence on simulating near-surface temperature and relative humidity, which was much greater than that of other physical processes. In addition, the choice of the radiation scheme had a significant impact on how the temperature was distributed close to the ground and how the wind field was simulated. Furthermore, the configuration of the MP scheme was found to exert a certain influence on the simulation of relative humidity; however, it demonstrated a weak influence on other meteorological variables. Secondly, The MYNN3 scheme for PBL process, the NoahMP scheme for LSF process, the WSM5 scheme for MP process, the RRTMG scheme for LW process, and the Dudhia scheme for SW process are found to be the comprehensive optimal physical process parameterization scheme combination for simulating meteorological variables in the research area selected in this study. As evident from the findings, the use of the OED method to obtain the combinations of the optimal physical process parameterization scheme could successfully reproduce the wind field, temperature, and relative humidity in the current study. Thus, this method appears to be highly reliable and effective for use in the WRF models to explore the optimal combinations of the physical process parameterization scheme, which could provide theoretical support to quickly analyzing accurate meteorological field data for longer periods and contribute to deeply investigating the migration and diffusion behavior of airborne pollutants in the atmosphere.

Weather Research and Forecasting Model (WRF) Sensitivity to Choice of Parameterization Options over Ethiopia

Downscaling seasonal climate forecasts using regional climate models (RCMs) became an emerging area during the last decade owing to RCMs’ more comprehensive representation of the important physical processes at a finer resolution. However, it is crucial to test RCMs for the most appropriate model setup for a particular purpose over a given region through numerical experiments. Thus, this sensitivity study was aimed at identifying an optimum configuration in the Weather, Research, and Forecasting (WRF) model over Ethiopia. A total of 35 WRF simulations with different combinations of parameterization schemes for cumulus (CU), planetary boundary layer (PBL), cloud microphysics (MP), longwave (LW), and shortwave (SW) radiation were tested during the summer (June to August, JJA) season of 2002. The WRF simulations used a two-domain configuration with a 12 km nested domain covering Ethiopia. The initial and boundary forcing data for WRF were from the Climate Forecast System Reanalysis (CFSR). The simulations were compared with station and gridded observations to evaluate their ability to reproduce different aspects of JJA rainfall. An objective ranking method using an aggregate score of several statistics was used to select the best-performing model configuration. The JJA rainfall was found to be most sensitive to the choice of cumulus parameterization and least sensitive to cloud microphysics. All the simulations captured the spatial distribution of JJA rainfall with the pattern correlation coefficient (PCC) ranging from 0.89 to 0.94. However, all the simulations overestimated the JJA rainfall amount and the number of rainy days. Out of the 35 simulations, one that used the Grell CU, ACM2 PBL, LIN MP, RRTM LW, and Dudhia SW schemes performed the best in reproducing the amount and spatio-temporal distribution of JJA rainfall and was selected for downscaling the CFSv2 operational forecast.

The Improvement of Short-Term Quantitative Precipitation Forecast in Mountainous Areas by the Assimilation of Meteorological State Variables Retrieved by Multiple Doppler Radar Data

Abstract This study presents a sequential procedure formulated by combining a multiple-Doppler radar wind synthesis technique with a thermodynamic retrieval method, which can be applied to retrieve the three-dimensional wind, pressure, temperature, rainwater mixing ratio, and moisture over complex terrain. The retrieved meteorological state variables are utilized to reinitialize a high-resolution numerical model, which then carries out time integration using four different microphysical (MP) schemes, including the Goddard Cumulus Ensemble (GCE), Morrison (MOR), WRF single-moment 6-class (WSM6), and WRF double-moment 6-class (WDM6) schemes. It is found that through this procedure, the short-term quantitative precipitation forecast (QPF) skill of a numerical model over mountainous areas can be significantly improved up to 6 h. The moisture field plays a crucial role in producing the correct rainfall forecast. Since no specific microphysical scheme outperforms the others, a combination of various rainfall scenarios forecasted by different MP schemes is suggested in order to provide a stable and reliable rainfall forecast. This work also demonstrates that, with the proposed approach, radar data from only two volume scans are sufficient to improve the rainfall forecasts. This is because the unobserved meteorological state variables are instantaneously retrieved and directly used to reinitialize the model, thereby the model spinup time can be effectively shortened.

Arctic mixed-phase clouds simulated by the WRF model: Comparisons with ACLOUD radar and in situ airborne observations and sensitivity of microphysics properties

The microphysical properties of two mixed-phase clouds (MPCs), one over sea ice, and another over open ocean, are investigated using in situ measuring probes and a cloud radar installed on-board two aircraft during the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) field campaign, west of Svalbard on 17 June 2017. Comparisons between observations and different configurations of the Weather and Research Forecasting (WRF) model are conducted and reveal discrepancies in the vertical profile of the ice crystal concentrations for both MPCs due to slightly higher temperatures (by approx. 1–2 °C) in WRF than observed, affecting the ice formation. A comprehensive sensitivity analysis is carried out to address the impact of assumptions used for the microphysics (MP) and the atmospheric boundary layer (ABL) schemes on the simulated vertical properties of the MPC systems. The role of the vertical resolution of the large scale data used for model initialisation is also investigated. Results show that the model performance improves for large-scale initialisation with more numerous levels close to the surface. Moreover, modifying the ABL or MP scheme influences the altitudes of the cloud top and base as well as the liquid and ice water amounts. Our study also highlights how an increase in the cloud condensation nuclei and droplet number concentrations can impact the vertical profiles of the liquid and ice water amounts.

Sensitivity analysis of cumulus and microphysics schemes in the WRF model in simulating Extreme Rainfall Events over the hilly terrain of Nagaland

A comprehensive study is carried out to investigate the sensitivity of various combinations between cumulus (CU) and microphysics (MP) schemes using the Weather Research and Forecasting (WRF) model for simulating Extreme Rainfall Events (EREs) over the hilly terrain of Nagaland (25.1 oN – 27.2 oN; 93.2 oE– 95.3 °E). The study explores fifty-four distinct simulation combinations resulting from six CU and nine MP schemes. The simulated rainfall is compared with ground-based rain gauge network and Global Precipitation Measurement (GPM) satellite observations. The investigation focuses primarily on finding the optimal CU-MP scheme combination based on spatiotemporal precision achieved by the identified schemes. The C14 (New Simplified Arakawa Schubert) cumulus scheme is found to be the optimum individual CU scheme. Similarly, the M14 (WDM5) microphysics scheme performs as the optimum individual MP scheme. The study highlights the significance of evaluating the performance of combined schemes over individual ones. Among the fifty-four combinations, sensitivity analysis identified six top-performing combinations, out of which the C11-M8 (Multiscale Kain Fritsch - New Thompson) combination stands out to be optimum with a robust correlation coefficient (CC) of 0.75, lower standard deviation (STD) of 0.76 mm, and RMSD of 0.51 mm. A variation in lead/lag time is found to vary from (−) 3.2 to (+) 5.0 h, with a median value of (+) 1.3 h, suggesting that, on average, the model tends to predict the occurrence of peak rainfall events slightly earlier than observed one. The absolute difference in latitudes (longitudes) of peak rainfall variation ranges from 0.42o - 0.8o (0 o - 0.51o) with a median value of 0.66o (0.37o), indicating that the model could better represent the longitudinal peak rainfall positions with notable accuracy. The spatial correlation coefficient (SCC) ranges between 0.31 and 0.68 with a median value of 0.50. From the overall analysis, it is also found that double moment MP schemes (New Thompson, Morrison (2 M), WDM5, WDM6) consistently perform better than the single moment MP schemes (Kessler, Lin et al., WSM3, WSM5, WSM6), whether in combination with CU schemes or individually. Spatiotemporal variations observed across different cases underscore the pivotal role of accurate topography representation in complex terrains, emphasizing the significance of terrain data in shaping rainfall patterns.

Key ingredients in regional climate modelling for improving the representation of typhoon tracks and intensities

Abstract. There is evidence of an increased frequency of rapid intensification events of tropical cyclones (TCs) in global offshore regions. This will not only result in increased peak wind speeds but may lead to more intense heavy precipitation events, leading to flooding in coastal regions. Therefore, high impacts are expected for urban agglomerations in coastal regions such as the densely populated Pearl River Delta (PRD) in China. Regional climate models (RCMs) such as the Weather Research and Forecasting (WRF) model are state-of-the-art tools commonly applied to predict TCs. However, typhoon simulations are connected with high uncertainties due to the high number of parameterization schemes of relevant physical processes (including possible interactions between the parameterization schemes) such as cumulus (CU) and microphysics (MP), as well as other crucial model settings such as domain setup, initial times, and spectral nudging. Since previous studies mostly focus on either individual typhoon cases or individual parameterization schemes, in this study a more comprehensive analysis is provided by considering four different typhoons of different intensity categories with landfall near the PRD, i.e. Typhoon Neoguri (2008), Typhoon Hagupit (2008), Typhoon Hato (2017), and Typhoon Usagi (2013), as well as two different schemes for CU and MP, respectively. Moreover, the impact of the model initialization and the driving data is studied by using three different initial times and two spectral nudging settings. Compared with the best-track reference data, the results show that the four typhoons show some consistency. For track bias, nudging only horizontal wind has a positive effect on reducing the track distance bias; for intensity, compared with a model explicitly resolving cumulus convection, i.e. without cumulus parameterization (CuOFF; nudging potential temperature and horizontal wind; late initial time), using the Kain–Fritsch scheme (KF; nudging only horizontal wind; early initial time) configuration shows relatively lower minimum sea level pressures and higher maximum wind speeds, which means stronger typhoon intensity. Intensity shows less sensitivity to two MP schemes compared with the CuOFF, nudging, and initial time settings. Furthermore, we found that compared with the CuOFF, using the KF scheme shows a relatively larger latent heat flux and higher equivalent potential temperature, providing more energy to typhoon development and inducing stronger TCs. This study could be used as a reference to configure WRF with the model's different combinations of schemes for historical and future TC simulations and also contributes to a better understanding of the performance of principal TC structures.

Fidelity of WRF model in simulating heat wave events over India

The evaluation of Weather Research and Forecasting (WRF) model has been performed for simulating episodic Heat Wave (HW) events of 2015 and 2016 with varied horizontal resolutions of 27 km for the entire India (d01), 9 km for the North West (NW (d02)) and South East (SE (d03)) domain. Study compares the maximum temperature (Tmax) simulated by WRF model, using six different combination of parameterization schemes, with observations from the India Meteorological Department (IMD) during the HW events. Among the six experiments, Exp2 (i.e., combination of WSM6 microphysics (MP) together with radiation parameterization CAM, Yonsei (PBL), NOAH land surface and Grell-3D convective schemes) is found closest to the observations in reproducing the temperature. The model exhibits an uncertainty of ± 2 °C in maximum temperature (Tmax) for both the regions, suggesting regional temperature is influenced by the location and complex orography. Overall, statistical results reveal that the best performance is achieved with Exp2. Further, to understand the dynamics of rising HW intensity, two case studies of HW days along with influencing parameters like Tmax, RH and prevailing wind distribution have been simulated. Model simulated Tmax during 2015 reaches up to 44 °C in NW and SE part of India. In 2016, HW is more prevailing towards NW, while in SE region Tmax reaches upto 34–38 °C with high RH (60–85%). The comparative research made it abundantly evident that these episodic events are unique in terms of duration and geographical spread which can be used to assess the WRF performance for future projections of HW.

Sensitivity of Different Physical Schemes in WRF Model of a Rainfall Event in Baghdad Station

The Weather Research and Forecasting model (WRF) offers a number of physical options that let users modify it to different scales, regions, and applications. The aim of this study is to test the sensitivity of different physics schemes in the WRF model for rainfall events over Iraq. In this study, six different physics configurations of the climate version of WRF were evaluated for simulation of a rainfall event in Iraq. Possible combinations among two Planetary Boundary Layers (PBL), three Cumulus (CUM) and two Microphysics (MIC) schemes were tested. The study area is the region surrounded by the longitudes 35o E-55o E and latitudes 290o N–38o N, which typically includes the Iraq region. The WRF model is installed on a Linux platform with a 10 km grid size in the zonal and meridional directions. For the six different simulations and the process of choosing the best performing configuration for the Iraq region, the model outputs tested for a single grid point (Baghdad station) of the atmospheric parameters (temperature, pressure and total precipitation) with modeled data and ECMWF. Model outputs using statistical methods: Bias Error (BE), Mean Absolute Error (MAE) and Root Mean Square Error (RMSE). The results show All the simulations predict rainfall with values close to the actual but it was discovered that the cloud microphysics setup had the greatest impact on temperature biases, whereas the cumulus parameterization setup has the greatest impact on precipitation.

Assessment of physical schemes for WRF model in convection-permitting mode over southern Iberian Peninsula

Convection-permitting models (CPMs) enable the representation of meteorological variables at horizontal high-resolution spatial scales (≤ 4 km), where convection plays a significant role. In this regard, physical schemes need to be evaluated considering factors in the studied region such as orography and climate variability. This study investigates the sensitivity of the Weather Research and Forecasting (WRF) model as CPM to the use of different physics schemes on Andalusia, a complex orography region in the southern part of the Iberian Peninsula (IP). To do that, a set of 1-year WRF simulations was completed based on two “one-way” nested domains: the parent domain (d01) spanning the entire IP with 5 km spatial resolution and the nested domain (d02) for the region of Andalusia at 1 km of spatial resolution. 12 physic schemes were examined from combinations of microphysics (MP) schemes including THOMPSON, WRF single moment 6-class (WSM6), and WRF single moment 7-class (WSM7), and different options for the convection in d01, the Grell 3D (G3), Grell-Freitas (GF), Kain-Fritsch (KF), and deactivated cumulus parameterization (OFF). The simulated precipitation and 2-m temperature for the year 2018, characterized to be a very wet year, were compared with observational datasets from different sources to determine the optimal WRF configuration, including point-to-point and station-point comparisons at different time aggregations (from annual to hourly). In general, greater differences were shown when comparing the results of convection schemes in d01. Simulations completed with GF or OFF presented better performance compared to the reference datasets. Concerning the MP, although THOMPSON showed a better fit in high mountain areas, it generally presents a worse agreement with the reference datasets. In terms of temperature, the results were very similar and, therefore, the selection of the “best” configuration was based mainly on the precipitation results with the WSM7-GF scheme being suitable for Andalusia region.

Wind field simulation using WRF model in complex terrain: A sensitivity study with orthogonal design

Accurately simulating wind speed is of utmost importance for wind power assessments. The objective of this study is to investigate the performance of wind speed simulation in complex terrain using different parameterization schemes from the Weather Research and Forecasting (WRF) model, employing an orthogonal design methodology. Specifically, nine WRF simulations are conducted based on the orthogonal test, considering various configurations of the planetary boundary layer (PBL), microphysics (MP), and land surface (LS) options. The numerical results are then compared to actual wind data obtained from a measuring station at two different heights (50 m and 80 m) during summer and winter periods. Furthermore, range and variance analyses are employed to rank the three schemes and identify the optimal combination. Moreover, we examine the impact of each parameterization scheme on the accuracy of wind speed predictions based on the results obtained from the orthogonal simulations. Additionally, we discuss the influence of using different evaluation indices within the orthogonal test on the outcomes and analyze the WRF simulated results under the optimal scheme combination. Lastly, we conduct an uncertainty analysis of results from the optimal scheme combination. The findings reveal that both the PBL and MP schemes exhibit highly significant effects on the accuracy of wind speed predictions (significance: **), followed by the LS scheme (significance: *). The order of importance for these three options is ranked as follows: PBL > MP > LS, which is independent of the seasons. The optimal configurations vary from summer and winter periods. Specifically, the optimal scheme combination is determined to be PBL-ACM2, MP-Kessler, and LS-Noah MP in summer, while PBL-BouLac, MP-Lin, and LS-Noah MP in winter. The simulating accuracy of the wind speeds is satisfactory under this optimal combination when considering the uncertainty of on-site measurements during these two periods. These results provide valuable insights for selecting appropriate PBL, MP, and LS options (from the WRF model) for wind speed estimates and wind power development in the studied region.

Simulating Heavy Rainfall Associated with Tropical Cyclones and Atmospheric Disturbances in Thailand Using the Coupled WRF-ROMS Model—Sensitivity Analysis of Microphysics and Cumulus Parameterization Schemes

Predicting heavy rainfall events associated with Tropical Cyclones (TCs) and atmospheric disturbances in Thailand remains challenging. This study introduces a novel approach to enhance forecasting precision by utilizing the coupled Weather Research and Forecasting (WRF) and Regional Oceanic Model (ROMS), known as WRF-ROMS. We aim to identify the optimal combination of microphysics (MP) and cumulus (CU) parameterization schemes. Three CU schemes, namely, Betts-Miller-Janjic (BMJ), Grell 3D Ensemble (G3), and Kain-Fritsch (KF), along with three MP schemes, namely, Eta (ETA), Purdue Lin (LIN), and WRF Single-moment 3-class (WSM3), are selected for the sensitivity analysis. Seven instances of heavy (35.1–90.0 mm) to violent (>90.1 mm) rainfall in Thailand, occurring in 2020 and associated with tropical storms and atmospheric disturbances, are simulated using all possible combinations of the chosen physics schemes. The simulated rain intensities are compared against observations from the National Hydroinformatics Data Center. Performance was assessed using the probability of detection (POD), false alarm ratio (FAR), and critical success index (CSI) metrics. While the models performed well for light (0.1–10.0 mm) to moderate (10.1–35.0 mm) rainfall, forecasting heavy rainfall remained challenging. Certain parameter combinations showed promise, like BMJ and KF with LIN microphysics, but challenges persisted. Analyzing density distribution of daily rainfall, we found effective parameterizations for different sub-regions. Our findings emphasize the importance of tailored parameterizations for accurate rainfall prediction in Thailand. This customization can benefit water resource management, flood control, and disaster preparedness. Further research should expand datasets, focusing on significant heavy rainfall events and considering climate factors, for example, the Madden-Julian Oscillation (MJO) for extended-range forecasts, potentially contributing to sub-seasonal and seasonal (S2S) predictions.

WRF Physics Ensemble Performance Evaluation over Continental and Coastal Regions in Germany

WRF is used to simulate eight extreme precipitation events that occurred over the regions of Schleswig–Holstein and Baden–Wurttemberg in Germany. The events were chosen from the German Weather Service (DWD) catalog and exceeded the DWD’s warning level 3 (i.e., rainfall > 40 mm/h). A two-way nesting approach is used with 9 and 3 km spatial resolutions. Initial and boundary conditions are obtained from the ERA5 dataset at 0.25° × 0.25°. To model each event, thirty different parameterization configurations were used, accounting for all possible combinations of five microphysics (MP), three cumulus (CU), and two planetary boundary layer (PBL) parameterization methods, yielding a total of 240 simulations. TOPSIS multicriteria analysis technique is employed to determine the performance skill of each setup and rank them, using six categorical and five statistical metrics. Resolution increase from 9 to 3 km did not improve forecasting accuracy temporally or in intensity. According to TOPSIS ranking, when treating each event individually, the ideal parameterizations combination is spatiotemporally dependent, with certain members ranking higher. When all events are considered, the Morrison double-moment MP–Grell–Freitas CU–YSU PBL combination works best with a frequency of occurrence in the top five performing scenarios of 30%, 47.5%, and 57.5% respectively.

Optimizing simulation of summer precipitation by weather research and forecasting model over the mountainous southern Tibetan Plateau

The mountainous regions of southern Tibetan Plateau and surroundings (hereafter mountainous southern TP), mainly including Himalayas, Gangdise, Nyainqentanglha and Hengduan Mountains, are characterized by high mountains and a complex terrain, and play a crucial role in the redistribution of monsoon water vapor and consequently determine the spatial distribution of precipitation over the TP. However, numerical simulations generally have large precipitation deviation in this focus area. Therefore, optimizing simulation of summer precipitation by Weather Research and Forecasting (WRF) model with two-way nested domains was carried out. A matrix of 11 model experiments, including four cumulus convection (CU) schemes, four microphysics (MP) schemes and four planetary boundary layer (PBL) schemes, were designed. In situ observations from 44 stations, TRMM 3B42 v7 precipitation product, and statistical and Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) methods were utilized for the model evaluation. The results showed that physical parameterization schemes, Kain-Fritsch-Cumulus Potential scheme (CU), Thompson scheme (MP), Mellor-Yamada-Janjic scheme (PBL), Unified Noah land surface model, and Eta similarity surface layer scheme, greatly reduced the summer precipitation simulation error of newly High Asia Refined analysis (HAR v2) over the mountainous southern TP, with the better reproduction of the spatiotemporal distribution and vertical variability of precipitation. Compared to NCEP-FNL reanalysis data, the performance of this model setup forced by ERA5 was notably improved, especially with perspective to the vertical variability of precipitation. Our findings suggested that summer precipitation simulation was more sensitive to the CU scheme, and the choice of Kain-Fritsch-Cumulus Potential scheme (ERA5 as initial conditions) significantly improved the precipitation simulation, especially with the lowest MAE (2.00 mm/d), RMSE (2.89 mm/d) and highest Rs (0.74). Further analysis implied that accurate description of shallow convection and proper suppression of deep convection should be considered in precipitation simulations over the TP, especially in mountainous regions. This study could provide a firm foundation for the production of more realistic dynamic downscaling precipitation dataset over the TP.

Extreme Precipitation Events on the East Coast of Brazil’s Northeast: Numerical and Diagnostic Analysis

The Northeast of Brazil (NEB) is the region with the highest number of municipal decrees of emergency situation declaration caused by weather events in the period from 2013 to 2022 and with the highest rate of natural disasters per risk area. In the NEB, the city of Recife and its metropolitan region are the biggest localities with populations in risk areas. Focusing on this region, five events of extreme precipitation were chosen for simulations using the WRF model and diagnostics analyses. First, a set of configurations of the model was tested, including 11 microphysics (MPH) schemes, 9 planetary boundary layer (PBL) schemes, 5 cumulus (CUM), and 7 surface layer (SFC) schemes. Then, through diagnostic analysis, the conditional instability, the moisture supply at low levels, and the support of the medium and high levels in storm formation were verified. The model’s configurations were verified by 298 rain gauges with hourly registrations through statistical metrics such as bias, MSE, standard deviation, and Pearson’s correlation, and demonstrated that the MPH schemes of Thompson Aerosol-Aware and NSSL + CCM, ACM2, MYJ for the PBL, KFCuP for CUM, and RUC for SFC were considered the best. All the cases were better with CUM parametrizations turned on. In all cases, diagnostics analyses highlighted the strong moisture flux convergence at the low levels, the presence of wind shear on the middle layer, weak cyclonic vorticity advection at high levels, and CAPE values around 1500 J/kg, in addition to an inverse relationship between wind shear action and CAPE values. This work is part of the national strategy for monitoring, diagnosis, and modeling of information that can minimize or even prevent damage caused by severe precipitation events.

Assessment of Weather Research and Forecasting (WRF) Physical Schemes Parameterization to Predict Moderate to Extreme Rainfall in Poorly Gauged Basin

Incomplete hydro-meteorological data and insufficient rainfall gauges have caused difficulties in establishing a reliable flood forecasting system. This study attempted to adopt the remotely sensed hydro-meteorological data as an alternative to the incomplete observed rainfall data in the poorly gauged Kuantan River Basin (KRB), the main city at the east coast of Peninsula Malaysia. Performance of Weather Research and Forecasting (WRF) schemes’ combinations, including eight microphysics (MP) and six cumulus, were evaluated to determine the most suitable combination of WRF MPCU in simulating rainfall over KRB. All the obtained results were validated against observed moderate to extreme rainfall events. Among all, the combination scheme Stony Brook University and Betts–Miller–Janjic (SBUBMJ) was found to be the most suitable to capture both spatial and temporal rainfall, with average percentage error of about ±17.5% to ±25.2% for heavy and moderate rainfall. However, the estimated PE ranges of −58.1% to 68.2% resulted in uncertainty while simulating extreme rainfall events, requiring more simulation tests for the schemes’ combinations using different boundary layer conditions and domain configurations. Findings also indicate that for the region where hydro-meteorological data are limited, WRF, as an alternative approach, can be used to achieve more sustainable water resource management and reliable hydrological forecasting.

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

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.

Evaluating the Impact of Planetary Boundary Layer, Land Surface Model, and Microphysics Parameterization Schemes on Simulated GOES-16 Water Vapor Brightness Temperatures

The impact of several land surface models (LSMs) and microphysics (MP), planetary boundary layer (PBL), and surface layer schemes on the accuracy of simulated brightness temperatures (BTs) from water vapor (WV) sensitive bands was examined via comparison with observations from the GOES-16 Advanced Baseline Imager. Nine parameterization configurations were evaluated. Analysis revealed that, compared to the Thompson MP scheme, the National Severe Storms Laboratory MP scheme produced lower simulated WV BTs in the upper troposphere but higher WV BTs in the middle and lower troposphere. The configuration with the Geophysical Fluid Dynamics Laboratory MP and hybrid eddy-diffusivity mass-flux (EDMF) PBL instead of Mellor–Yamada–Nakanishi–Niino (MYNN) PBL produced higher BTs. Yet, changing the PBL from MYNN to Shin–Hong or EDMF reduced the simulated WV BTs. Changing the LSM from Noah to RUC also resulted in lower simulated WV BTs, which were further enhanced with the MYNN surface layer instead of the GFS. The location and orientation of upper-level jet streams and troughs was assessed using the location of WV gradient objects. Every configuration had an increased translation speed compared to the observations, as forecast WV gradient objects were west of the observation objects early in the forecast and then east later in the forecast.

Optimization of physical schemes in WRF model on cyclone simulations over Bay of Bengal using one-way ANOVA and Tukey\u2019s test

Evaluation of appropriate physics parameterization schemes for the Weather Research and Forecasting (WRF) model is vital for accurately forecasting tropical cyclones. Three cyclones Nargis, Titli and Fani have been chosen to investigate the combination of five cloud microphysics (MP), three cumulus convection (CC), and two planetary boundary layer (PBL) schemes of the WRF model (ver. 4.0) with ARW core with respect to track and intensity to determine an optimal combination of these physical schemes. The initial and boundary conditions for sensitivity experiments are drawn from the National Centers for Environmental Prediction (NCEP) global forecasting system (GFS) data. Simulated track and intensity of three cyclonic cases are compared with the India Meteorological Department (IMD) observations. One-way analysis of variance (ANOVA) is applied to check the significance of the data obtained from the model. Further, Tukey’s test is applied for post-hoc analysis in order to identify the cluster of treatments close to IMD observations for all three cyclones. Results are obtained through the statistical analysis; average root means square error (RMSE) of intensity throughout the cyclone period and time error at landfall with the step-by-step elimination method. Through the elimination method, the optimal scheme combination is obtained. The YSU planetary boundary layer with Kain–Fritsch cumulus convection and Ferrier microphysics scheme combination is identified as an optimal combination in this study for the forecasting of tropical cyclones over the Bay of Bengal.

Impact of meteorological parameterization schemes on CTM model simulations

Meteorology plays a key role in regional aerosol concentration and distribution. Microphysics and cloud formation processes in the atmosphere is interlinked with aerosol and their removal processes. To characterize their role in modulating the aerosols/pollutants concentration simulation by a chemical transport model (CTM, here we have taken CHIMERE model); we have performed this study with different microphysical (MP) schemes (Kessler Scheme/Lin Scheme/WRF single Moment 3-class (WSM3) scheme) and cumulus cloud parameterization (CU) schemes (Kain-Fritsch Scheme/Betts-Miller-Janjic Scheme/Grell 3D) of Weather Research and Forecasting Model (WRF). We have used the WRF model over a domain (3°S-41.8°N; 59.5°-102.5°E) with the resolution of (0.25° × 0.25°) while CHIMERE model simulations are performed over the domain (6°-37.5°N; 67°-95.5°E) with the similar resolution for a selected period of monsoon. In total nine combinations are framed with MP and CU schemes to observe the sensitivity of CTM to these schemes. Simulated results are compared with the satellite (TRMM/MODIS) and reanalysis data (MERRA-2) to appraise the model's performance with various parameterization scheme combinations. Results indicate that despite the same initial and boundary conditions and model configuration, notable differences occurred in the simulated meteorological parameters with different scheme combinations. Results suggested that CTM performed in a reliable range with cumulus scheme Betts-Miller-Janjic Scheme (BMJ) and Grell 3D scheme with microphysical parameterization scheme Purdue Lin Scheme over Indian continent. The study also suggests that a single set may not produce good results for all the parameters/pollutants; thus, we have to choose the parameterization schemes which give optimal results for all the parameters/pollutants. These results also infer that ensemble modeling could also lead to a better option than using single simulations.

A micro-genetic algorithm (GA v1.7.1a) for combinatorial optimization of physics parameterizations in the Weather Research and Forecasting model (v4.0.3) for quantitative precipitation forecast in Korea

Abstract. One of the biggest uncertainties in numerical weather predictions (NWPs) comes from treating the subgrid-scale physical processes. For more accurate regional weather and climate prediction by improving physics parameterizations, it is important to optimize a combination of physics schemes and unknown parameters in NWP models. We have developed an interface system between a micro-genetic algorithm (µ-GA) and the WRF model for the combinatorial optimization of cumulus (CU), microphysics (MP), and planetary boundary layer (PBL) schemes in terms of quantitative precipitation forecast for heavy rainfall events in Korea. The µ-GA successfully improved simulated precipitation despite the nonlinear relationship among the physics schemes. During the evolution process, MP schemes control grid-resolving-scale precipitation, while CU and PBL schemes determine subgrid-scale precipitation. This study demonstrates that the combinatorial optimization of physics schemes in the WRF model is one possible solution to enhance the forecast skill of precipitation.