WRF Simulations of Extreme Rainfall over Uganda’s Lake Victoria Basin: Sensitivity to Parameterization, Model Resolution and Domain Size ()
1. Introduction
Rainfall is one of the important weather elements and its forecasting is crucial to society (Gouda et al., 2018; Mehr et al., 2019). In Uganda, occurrence of rainfall extremes has been associated with a variety of negative socio-economic impacts on human livelihood. For example, extremely heavy rainfall plays a principal role in soil erosion (Bamutaze et al., 2017), flooding (Lwasa, 2010), landslides (Mugagga et al., 2012) and transmission of waterborne diseases (Cann et al., 2013). Consequently, it causes substantial loss in rain-dependent sectors such as agriculture, public health, transport and disaster management among many others.
One direct way to tackle this problem of extreme rainfall is to use numerical weather prediction (NWP) models to forecast its occurrence. For this, the Uganda National Meteorological Authority has integrated NWP modelling in their weather forecasting service to support the generation of short and medium range forecasts and their model of choice is the Weather Research and Forecasting (WRF) model. This model’s user community is on the rise signifying the confidence that scientists have in it (Warner, 2011). Also, it can be easily set-up for operational use since its outputs can be autonomously analyzed using existing state-of-the-science tools such as the Atmospheric Model Evaluation Tool (Appel et al., 2011).
As a major action, it is important to customize the WRF model over different regions of the country. This work focused on the Lake Victoria Basin (LVB) in Uganda (Figure 1(a)). Climatological zones in this region have been documented to receive substantially heavy rainfall (Majaliwa et al., 2015), therefore it seemed logical to explore the rainfall records of this region in the quest for extremely heavy rainfall. For this, the spatial rainfall anomaly for the longest rainfall season, March to May (MAM) was calculated based on TRMM rainfall for 20 years, 1998 to 2017. The analysis revealed the MAM season of 2008 as one which experienced extremely wet conditions over the land area (Figure 1(b)). The WRF model was then used to simulate a 20-day period (17th March 2008 to 5th April 2008) of extremely heavy rainfall identified within this season.
The aspects of the customization done in this work include: 1) identifying combinations of cumulus and microphysical parameterization that generate skillful simulations of extreme rainfall over the region; 2) determining the sensitivity of the simulations to increasing the model resolution and changing the domain size. This work used 2 independent domain set-ups (Table 1).
2. Data and Methods
2.1. Data
Boundary conditions
The WRF model simulations were initialized using Final Operational Model
Figure 1. (a) Map of the Lake Victoria Basin in Uganda. The dotted rectangles are the representative regions of the basin. Region 1 is western LVB, region 2 is central LVB and region 3 is eastern LVB. (b) Spatial rainfall anomaly (mm) for March to May, 2008 based on TRMM rainfall for 20 years, 1998 to 2017.
Table 1. Domain setup for WRF model experiments.
Global Tropospheric Reanalyzes from the National Center for Environmental Prediction (NCEP FNL). This is 6-hourly data of 1˚ × 1˚ (110 km) resolution (Kalnay et al., 1996). The boundary conditions used were for 21 days, corresponding to the study period of 20 days with an additional day to cater for model spin-up.
Rainfall observations
The Tropical Rainfall Measuring Mission (TRMM) provided observations for model validation. Observed estimates of daily and 3-hourly rainfall accumulations were obtained from TRMM 3B42 version 7. This product is a 0.25˚ × 0.25˚ (27.5 km) resolution precipitation estimate obtained by applying the Goddard Profiling Algorithm to extract rainfall estimates from a variety of passive microwave sources (Huffman et al., 2007). The dataset was downloaded from the GIOVANNI data handle (https://giovanni.gsfc.nasa.gov/giovanni/) of the National Aeronautics and Space Administration (NASA). Data for a 20-day period was used for evaluating the model performance.
2.2. Methods
Model setup and experimental design
Simulation of extreme rainfall that occurred over the LVB was done using the WRF modelling system, Advanced Research WRF (ARW) core version 3.9. It is a grid point model with a terrain—following vertical coordinate near the surface which transforms to a constant pressure at the highest level. WRF performs numerical integration on a staggered Arakawa C-grid with 3rd order Runge-Kutta time integration (Skamarock et al., 2008). The model was setup with 27 vertical levels ending at the 50 hPa isobaric level. Three domains (Table 1 and Figure 2) of varying grid resolution were used with a superior 2-way domain interaction and 1:3 nesting ratio. The model was initialized using NCEP finalized reanalyzes starting 24 hours before the period of interest to allow time for model spin-up.
The domain setup was as follows; a coarse domain centered at 0.73˚ latitude, 22˚ longitude covering Africa to encompass large scale synoptic systems such as the Inter Tropical Convergence Zone and sub-tropical anti-cyclones; a nest domain d02, covering the Congo basin and the western region of the Indian Ocean to cater for moisture inflow from the Congo air mass and the Indian Ocean respectively; and a nest domain d03, covering LVB and Lake Victoria to cater for the influence of lake surface temperatures on rainfall in the basin. The ARW solver was used with a combination of physical parameterization (Table 2) coupled with static land-use data based on Moderate Resolution Imaging Spectroradiometer (MODIS) with 21 land-use categories (Friedl et al., 2002) and a special lake surface representation. The land-use was setup with varying resolution, the coarse domain used a 10-minute resolution while d02 and d03 used a 5-minute resolution.
To determine the model skill associated with the use of different cumulus – microphysical scheme combinations, the model setup was identical in all experiments, except for the cumulus and microphysical schemes (Table 2). Similar to Argent et al. (2015), the default parameterization settings of the WRF model were used as the control run.
Figure 2. Domain configuration for WRF simulations.
To determine the effect of changing grid resolution, experiments done in set-up 1 were repeated with the inner most domain set to 4 km resolution (Set-up 2). A comparison was then done between results of experiments at 12 km resolution and experiments at 4 km resolution. Furthermore, to reveal the effect of changing domain size, a comparison was done between results of domain, d03 in set-up 1 and domain, d02 in set-up 2. Both domains were set at 12 km resolution and only differed in size.
Model validation methods
The Difference variable. (Willmott et al., 1985) described in Equation (1) was used to generate a spatial difference between the WRF output and TRMM.
(1)
Continuous scores. The Root Mean Square Error (RMSE) (Chai & Draxler, 2014) was used to determine model accuracy, and the Mean Error (ME) (Mugume et al., 2016) was used to determine model bias.
(2)
(3)
In Equations (2) and (3), S is the simulated value, O is the observed value, and n is the total number of comparisons.
Categorical scores. The scores used were the Probability of Detection (POD) and False Alarm Ratio (FAR) described in Equations (4) and (5) respectively. POD and FAR values range from 0 to 1 and they are derived from a 2 × 2 contingency table (Table 3) described by Schaefer (1990).
(4)
(5)
3. Results and Discussion
3.1. WRF Simulations of Extreme Rainfall under Different Parameterization Settings
Observed and simulated spatial rainfall patterns
Figure 3 shows the rainfall distribution from TRMM observations and the WRF model runs for the 20-day period. TRMM shows that the total rainfall amount received in LVB during this period was between 100 and 280 mm although rainfall amounts between 120 to 200 mm dominated the area and central LVB received a higher rainfall amount compared to the other 2 regions. On the contrary, all model runs generated higher rainfall amounts in the higher altitude western LVB compared to the other 2 regions. This gives evidence of the terrain-sensitive nature of the WRF model which was also reported by Mugume et al. (2017) and Ntwali et al. (2016).
The model runs also generated varying rainfall distributions when different cumulus and microphysical parameterization combinations were used. Overall, the rainfall was in the range of 10 to over 320 mm. Notably, runs done with the SBU_YLin microphysical scheme (Figure 3(e), Figure 3(h), & Figure 3(k))
Figure 3. Total rainfall amount (mm) from TRMM data (a), the control run (b) and experimental runs E1 to E9 (c to k). The plots c to k are also arranged following the cumulus schemes (rows) and microphysical schemes (columns) used.
generated higher rainfall amounts over LVB as compared to those done with other microphysical schemes. This rainfall disparity caused by the microphysical scheme is because they each use different techniques for predicting hydrometeor species.
The difference maps (Figure 4) show that all model runs underestimated the rainfall in the watershed closet to the lake. Also, except for the runs done with the Grell 3D cumulus scheme, the other runs overestimated rainfall in the highland areas in the southwestern corner of the basin. Of the areas in which the model runs show close agreement, the error ranged from 25 to 50 mm and very few areas had an error below 25 mm.
Observed and simulated temporal rainfall patterns
Figure 5 shows the temporal rainfall accumulation (a, c, e) and the diurnal evolution of rainfall (b, d, f) during the 20-day period for western, central and eastern LVB respectively. The temporal rainfall accumulation indicates that for western LVB, the model runs fairly simulate the rainfall amounts while for central and eastern LVB, there was a dominant underestimation of the rainfall. A further investigation into the diurnal variation (Figure 5(b), Figure 5(d), Figure 5(f)) alludes that the underestimation happens during the nocturnal hours of the day. For central LVB, the model, this is between 18:00 hrs. to 06:00 hrs. of the following day while for eastern LVB, this happens between 18:00 hrs. and 24:00 hrs.
Figure 4. Difference in total rainfall amount (mm) between WRF model output and TRMM data for the 20-day period.
Important to note is that TRMM observations (b, d, f) show that rainfall in the basin is mostly received at night and early in the morning. Sun et al. (2014) demonstrated that during this time, the lake is warmer than the adjacent land, favoring the occurrence of a land breeze that causes the rainfall. The WRF model, however, was unable to correctly simulate this nocturnal rainfall caused by the land breeze. As explained by Argent et al. (2015), there is insufficiency in the representation of the lake surface temperature within the model. The temperature of the lake is constantly colder than the adjacent land. This hinders the reversal of winds to blow from the land to the lake (that is, failure to generate the land breeze). The lack of nocturnal rainfall also explains why the model underestimated rainfall in central and eastern LVB. Mayor and Mesquita (2015) also acknowledged a similar failure of WRF over Cuba.
Skill scores
Table 4 summarizes the skill scores and corresponding ranks of model runs done using an inner-most domain set at 12 km resolution (Set-up 1). A one-sample t-test was used to test whether the skill scores generated by model runs E1 to E9 significantly differ from those of the control run and the results show that at a 95% level of confidence, changing the microphysical and cumulus parameterization schemes from the default significantly altered the model bias (ME, p < 0.05), detection ability (POD, p < 0.05) and proportion of false alarms (FAR, p < 0.05)
Figure 5. Temporal rainfall accumulation (mm) for the 20-day period (left column) and 3-hourly rainfall accumulation (mm) (right column) for western LVB (a), (b), central LVB (c), (d) and eastern LVB (e), (f).
but not the accuracy (RMSE, p > 0.05). Overall, model runs done with the Grell 3D cumulus scheme showed closer proximity to the TRMM observations and specifically, model run E6 done with the Grell 3D cumulus scheme combined with the SBU_YLin microphysical scheme showed the most satisfactory skill and generated the best spatial representation of the rainfall.
Table 4. Ranks of the model skill scores for runs done at 12 km resolution.
3.2. Effect of Grid Resolution on WRF Model Skill in Simulating Extreme Rainfall
Effect of model resolution on the skill scores
The 10 model runs done with the inner most domain set to 12 km (Set-up 1) were repeated with the inner most domain set to 4 km (Set-up 2) and the skill scores they generated were summarized in Table 5 along with the corresponding ranks. Model runs E4 and E6 which both use the Grell 3D cumulus scheme rank higher than the rest of the model runs.
Furthermore, the skill scores of the model runs at both resolutions were summarized in Table 6 and a paired-sample t-test was used to test the effect that changing model resolution from 12 km to 4 km has on the model skill. The results show that it significantly affected the accuracy (RMSE, p < 0.05), detection ability (POD, p < 0.05) and proportion of false alarms (FAR, p < 0.05) but did not significantly alter the model bias (ME, p > 0.05). Majority of the model runs generated lower RMSE, lower POD and higher FAR at 4 km resolution compared to the 12 km resolution. This implies that increasing the model’s resolution produced a higher accuracy, lowered the detection ability by 2% or more and increased the number of false alarms by 1% to 2%. For the Control run, E1, E2, E7 and E8, increasing the model resolution was of no additional value. It caused an undesirable increase in model bias.
Effect of domain size on the skill scores
The nature of the experimental set-up allowed for a domain-size sensitivity test to be carried out. For this, the paired-sample t-test was used to compare the skill scores generated by domain, d03 set at 12 km resolution (Set-up 1) to those of domain, d02 also set at 12 km resolution (Set-up 2). These domains only differed in size. The results (Table 7) suggest that changing the size of the domain significantly altered the model’s accuracy (RMSE, p < 0.05) and proportion of false alarms (FAR, p < 0.05). Majority of the model runs generated lower RMSE in the large domain. This increase in accuracy is because synoptic features are
Table 5. Ranks of the model skill scores for runs done at 4 km resolution.
Table 6. Skill scores generated by the 12 km and 4 km model runs.
Table 7. Skill scores generated by model runs when different domain sizes of the same resolution are used. The subscripts S and L denote small and large domain respectively.
better resolved in large domains. The large domain used in this test covered an approximate area of more than 15,000,000 km2, making it suitable to capture large scale features such as the ITCZ which is responsible for rainfall in LVB during the March to May season. However, the large domain was disadvantaged by its large spatial coverage because it registered a higher number of false alarms.
4. Conclusion and Recommendations
It is difficult to objectively identify the best parameterization combination because each combination ranked differently depending on the validation metric and the model resolution being considered. However, the Grell 3D cumulus scheme combined with the SBU_YLin microphysical scheme showed higher skill than the others and therefore can be considered as more appropriate for use in numerical simulations of extreme rainfall in equatorial regions.
The WRF model showed strong performance in western LVB. Overall bias was majorly contributed to by central and eastern LVB. The model’s failure to simulate nocturnal rainfall in these 2 regions can hinder its capability to generate reliable forecasts of nocturnal thunderstorms. This can also affect its long-term simulations of the rainfall climatology over this region. Furthermore, the benefit of increasing the model grid resolution depends on the parameterization schemes used. In this study, it was mostly useful for experiments done with the Grell 3D cumulus scheme. Therefore, this scheme is suitable for high resolution numerical simulations of rainfall in equatorial regions. Also, there is added value in using a large domain in comparison to a small one.
Acknowledgements
This work was supported by the WIMEA-ICT project (https://wimea-ict.net/) with funding from the Norwegian Agency for Development Cooperation (NORAD). The NCEP FNL reanalysis data was obtained from the NCAR research data archive and the TRMM rainfall observations were obtained from the Giovanni data handle of NASA.
NOTES
#These authors made equal contribution to the work.