water Article Potential of the Coupled WRF/WRF-Hydro Modeling System for Flood Forecasting in the Ouémé River (West Africa) Gandomè Mayeul Leger Davy Quenum 1,2,3,* , Joël Arnault 3, Nana Ama Browne Klutse 1,4, Zhenyu Zhang 3, Harald Kunstmann 3 and Philip G. Oguntunde 5 1 African Institute of Mathematical Sciences (AIMS), Sector Remera, Kigali 20093, Rwanda; nklutse@ug.edu.gh 2 National Institute of Water (NIW), University of Abomey-Calavi, Godomey, Cotonou 01 BP 4521, Benin 3 Institute for Meteorology and Climate Research, Karlsruhe Institute of Technology, Campus Alpine, 82467 Garmisch-Partenkirchen, Germany; joel.arnault@kit.edu (J.A.); zhenyu.zhang@kit.edu (Z.Z.); harald.kunstmann@kit.edu (H.K.) 4 Department of Physics, University of Ghana, Legon P.O. Box LG 63, Ghana 5 Department of Agricultural and Environmental Engineering, Federal University of Technology, Akure 340110, Nigeria; poguntunde@yahoo.com * Correspondence: malgdaq2000@yahoo.fr; Tel.: +229-9606-2003 Abstract: Since the beginning of the 2000s, most of the West-African countries, particularly Benin, have experienced an increased frequency of extreme flood events. In this study, we focus on the case of the Ouémé river basin in Benin. To investigate flood events in this basin for early warning, the coupled atmosphere–hydrology model system WRF-Hydro is used, and analyzed for the period 2008–2010. Such a coupled model allows exploration of the contribution of atmospheric components into the flood event, and its ability to simulate and predict accurate streamflow. The potential of  WRF-Hydro to correctly simulate streamflow in the Ouémé river basin is assessed by forcing the  model with operational analysis datasets from the European Centre for Medium-Range Weather Citation: Quenum, G.M.L.D.; Forecasts (ECMWF). Atmospheric and land surface processes are resolved at a spatial resolution Arnault, J.; Klutse, N.A.B.; Zhang, Z.; of 5 km. The additional surface and subsurface water flow routing are computed at a resolution of Kunstmann, H.; Oguntunde, P.G. Potential of the Coupled 500 m. Key parameters of the hydrological module of WRF-Hydro are calibrated offline and tested WRF/WRF-Hydro Modeling System online with the coupled WRF-Hydro. The uncertainty of atmospheric modeling on coupled results for Flood Forecasting in the Ouémé is assessed with the stochastic kinetic energy backscatter scheme (SKEBS). WRF-Hydro is able to River (West Africa). Water 2022, 14, simulate the discharge in the Ouémé river in offline and fully coupled modes with a Kling–Gupta 1192. https://doi.org/10.3390/ efficiency (KGE) around 0.70 and 0.76, respectively. In the fully coupled mode, the model captures w14081192 the flood event that occurred in 2010. A stochastic perturbation ensemble of ten members for three Academic Editor: Renato Morbidelli rain seasons shows that the coupled model performance in terms of KGE ranges from 0.14 to 0.79. Additionally, an assessment of the soil moisture has been developed. This ability to realistically Received: 1 March 2022 reproduce observed discharge in the Ouémé river basin demonstrates the potential of the coupled Accepted: 1 April 2022 WRF-Hydro modeling system for future flood forecasting applications. Published: 8 April 2022 Publisher’s Note: MDPI stays neutral Keywords: fully coupled WRF-Hydro modeling; flood forecasting; stochastic kinetic energy backscatter with regard to jurisdictional claims in scheme (SKEBS); Ouémé river basin published maps and institutional affil- iations. 1. Introduction In its fifth report, the Intergovernmental Panel on Climate Change (IPCC) stresses the Copyright: © 2022 by the authors. Licensee MDPI, Basel, Switzerland. increment of the number of extreme weather events for the 21st century due to climate This article is an open access article change [1]. Tropical countries of West Africa are threatened by climatic hazards, such as distributed under the terms and droughts, floods, high winds, the elevation of the sea level, etc. Droughts and floods are the conditions of the Creative Commons most important in terms of damages and impacts. Countries in West Africa such as Benin, Attribution (CC BY) license (https:// Burkina Faso, Cote d’Ivoire, Niger, Senegal, and Togo have suffered from catastrophic creativecommons.org/licenses/by/ floods [2] with severe consequences including loss of life, property, and damage of property. 4.0/). Floods happen in geomorphologically low-lying areas, such as rivers, basins, lakes, etc., Water 2022, 14, 1192. https://doi.org/10.3390/w14081192 https://www.mdpi.com/journal/water Water 2022, 14, 1192 2 of 23 and occur often with heavy rain associated with severe thunderstorms, hurricanes, tropical storms, or meltwater from ice or snow flowing over ice sheets or snowfields [3]. Ouémé river, the main watershed in the Benin Republic, is used for irrigation, hydroelectricity, transport, and water resource management. Its utility led to the construction of hydrological structures; however, this construction exposes the residents on the riverfronts to flooding risks. In 2010, Ouémé river has recorded an overflow driving the country (Benin) to experience dramatic flooding, which affected 680,000 people, leading to 43 deaths [4] over the length of the Ouémé river (at Bétérou, Savè, and Bonou outlets).. To prevent these damages during the rainy season, the ideal solution proposed by Benin is the improvement of national and decentralized capacities to better manage and plan for extreme weather events and climate-related risks through the development of an early warning system (EWS) and better climate information (IC). Therefore, citizens are earlier informed about the level of the risks based on color patterns, which depend on the water heights at the level of the water gauges. The flooding water heights limit is defined based on historical rainfall episode statistics. In addition, EWS warns when there is a risk that the rise of the water level upward of the river could spread downstream. Such destructive floods make it necessary to improve hydrological forecasts to reduce the vulnerability of regional and local communities. Identifying the drivers of change in flood regimes in West African watersheds is a complex task due to the heterogeneity of the region and the changing hydrological functioning of watersheds as a result of human activity. Ref. [5] found that there is a significant trend both in terms of flooding magnitude and frequency in West Africa, with two main patterns: the Sahelian showed increasing flood trends, whilst some Sudanian areas presented decreasing flood trends. Since IPCC [1] announced and highlighted an increase in the frequency of extreme climate events, and also the potentially significant increase of their magnitudes due to climate change, identification of accurate flood forecasting methods is required. The introduction of flood forecasting and early warning systems, the building of population awareness and preparedness, urban planning, discouraging human settlements in flood-prone areas, and the development of local institutional capacities, as well as forecasting methods development, should be appropriate actions to plan. To determine an accurate flood forecasting method, several studies applied hydrologi- cal modeling based on conceptual, semi-distributed, physical-based, data-driven or ma- chine leaning methods. The conceptual methods adopted in some previous studies [2,6–8] have been proved to be useful for hydrological modeling options because most are success- ful for rainfall-runoff simulation and flood modelling. The data-driven approach, based on the empirical relationship amongst datasets, which is usually driven through statistical constraints [9–11], has been investigated; despite their advantages (no need of physical equations and parameters or catchment characteristics), they present some limitations such as static models which cannot evaluate the change (e.g., land-use/land-cover change). With the increasing development in the machine leaning field, numerous methodologies are applied in operational flooding forecasting [12,13]. Thus, some studies [14–16] have proven that the machine learning-based model, performs better when it is compared with the physical-based, conceptual, and data-driven models. These methodologies do not take into account the influence of both atmosphere and land surface complexity in the performance of the forecasting. The operational global weather forecast centers routinely provide relatively coarse precipitation forecasts with resolutions of 16–27 km (e.g., ref. [17]). Because these kinds of forecasts could not provide necessary details of complex, intense precipitation structures led by the mesoscale orography, land-surface heterogeneities, and land-water contrasts, ref. [18,19] used the Weather Research and Forecast (WRF) model to provide high-resolution precipitation forecasts at a 1.3–4 km horizontal resolution. They found that the WRF model was able to provide precipitation forecast both in terms of amount and spatial distribution. Additionally, ref. [17] used the output from the WRF model (over the upper- Jordan River basin) as forcing to run the Hydrological Model for Karst Environment Water 2022, 14, 1192 3 of 23 (HYMKE). Similar experiments were set in different areas [20–23], and the authors showed that the precipitation estimates with WRF are able to reproduce the spatial distribution of precipitation, but may underestimate the magnitude of the heavy precipitation events when compared to rain gauges or global observation datasets. Refs. [24–30] have shown the benefits of choosing coupled atmospheric-land-surface models for estimating temperature and precipitation in different areas. Ref. [31] investigated and evaluated simulations with a fully coupled atmospheric–hydrological model (WRF- HMS) and an uncoupled model for various meteorological variables, and demonstrated the better performance of the fully coupled model against the uncoupled. Another study led by [32] compared a one-way forced implementation of the WRF-Hydro system to a fully two-way coupled instance of WRF/WRF-Hydro in order to assess the impact of two-way coupling on simulated precipitation and streamflow. They found that the two setups performed well for the precipitation, but the matching to observed data was higher for the two-way coupled WRF/WRF-Hydro simulation in terms of statistical performance criteria. Assessing the advantages and limitations of one-way vs. two-way coupled modeling systems for flood prediction over the Ayalon basin (Israel), ref. [17] used both a hydrological model (Hydrological Engineering Center-Hydrological Modeling System, HEC-HMS) and the WRF-Hydro modeling system. The models were forced by observed, interpolated precipitation from rain-gauges within the Ayalon basin, and with modeled precipitation from the WRF atmospheric model. Comparing simulations with the one-way coupled WRF model and the two-way coupled WRF/WRF-Hydro modeling system, they found that the use of two-way atmospheric-hydrological coupling has the potential to enhance precipitation simulation and, therefore, hydrological forecasts for early flood warning applications. It is important to acknowledge here that the model uncertainty in simulating precipitation, as well as streamflow, can be relatively high. Model uncertainty can, for example, be evaluated with a model ensemble (e.g., [33,34]). Ref. [35] developed the stochastic kinetic energy backscatter scheme (SKEBS) within the WRF model in order to generate such an ensemble. The present study attempts to develop a model system that is able to simulate flood events in the Ouémé river basin (Benin, West Africa). According to successful examples in the literature, we use the fully coupled atmospheric-hydrological modeling system WRF/WRF-Hydro for this purpose. Section 2 gives the characteristics of the study area, different sources of the dataset used, and the modeling approaches. Results and discussion are provided in Section 3, whilst Section 4 is dedicated to a summary and conclusion. 2. Materials and Methods 2.1. Study Area and Observation Datasets The study area is in West Africa, located between latitudes 0◦ N and 18◦ N and longitude 7◦ W and 12◦ E, as displayed in Figure 1. This also constitutes the setup of the WRF domain. The region is bordered in the south by the Gulf of Guinea, and in the north by Mali and Niger. Nigeria’s highlands form the eastern boundary, while Mauritania, Mali, and Ghana form the western limit. The annual mean temperature is about 18 ◦C, but the monthly mean can be more than 30 ◦C over the southern part of the Sahara. The rainfall pattern over this region is mostly affected by ocean currents and local features such as topography. In terms of climatic zones, West Africa is characterized by mainly three different regions: the first region covers the Sahel and is characterized as a semi-arid zone located from western Senegal to eastern Sudan between 12◦ N and 20◦ N. The second zone is the Sudano-Sahelian region, and the third zone comprises the Guinea coast, which is characterized by a bimodal mode driven by the inter-tropical divergence (ITD). The basin is located in the Benin Republic. Waatteerr 22002222,, 1144,, x1 1F9O2R PEER REVIEW 4 4oof f2233 FFiigguurree 11. .Weeaatthheerr Reesseeaarrcch Forrecastt ((WRF)) and WRF/WWRRFF-H-Hyyddrroo ssimimuulalatteedd ddoommaainin, ,aanndd ssttuuddiieedd ccaattcchhmeennttss ((SSaavvèè aanndd BBééttéérroouu)).. TThhee RReeppuubblliicc ooff BBeenniinn iiss iinn tthee iintteerr--ttrropiiccaall zonee ((bettweeeen 06°◦1100′ ′N aanndd 1122°◦2255′ ′NN),) , whhiicchh hhaass aa wweet tanandd drdyr ytrtorpoipciacla cllicmlimatae t[e2][.2 I]t. cIotnctoanintasi tnhset rhieverrisv e(Srsav(Sè aavnèda BnédtéBroétué)r oun ) wonhiwchh wiceh fwoceufso icnu tshiins sthtuisdyst (usdeey F(sigeeurFei g1u). rTeh1e) .OTuhéemOé ucéamtchémcaetncth amt ethnet aStavthèe (rSeasvpè. B(réetsép- . rBoéut:é irnonue:ri-ncnatecrh-cmatecnhtm toe nStavtoè)S oavuètl)eot ucotlveetrcso avne rasraena aorf e2a4o.8f0204 k.8m002 (krmes2p(. r1e0s.p4.7150 k.4m725).k Imt i2s) . lIotciastelodc bateetwd ebeentw 7°e5e8n′–71◦05°81′2–′1 N0◦ a1n2d′ N 1°a3n5′d–31°◦0355′ ′E–,3 a◦n0d5′ rEep, arensdenretsp 4re7s%en otfs t4h7e% whofolteh eOwuéhmolée rOivueér m(réefr.i v[3e6r])(.r eTfh. e[3 s6e]a).soTnhse csoerarseosnpsoncodr rtoes tphoen pdetroiotdhse opfe drioomdsinoafndceo mofi nthaen cweeotf ttrhoepiwcaelt ctoronptiincaenl ctaoln atiinr emntaaslsaeisr. mTahses esse.aTsohneasle adsiostnraibl udtiisotrni bouft iroaninofaf lrla finolflaollwfos ltlhowe sdtirheecdtiiornec otifo nthoef ItThDe IaTnDd avnadrievsa arilemsoasltm porsotpporrotipoonratilolyn walliythw diitshtadnicseta fnrocemf rtohme ctohaestc.o Tahste.rTefhoerree,f Boéreté, rBoéuté hraosu ah uasnaimuondimalo pdraelcpipreitcaiptiiotant iroengirmegei m(Me a(My atoy OtocOtocbteorb)e, rw),hwilhsitl stthteh seosuotuhtehrenr npparatr toof fththee SSaavvèè ccaattcchhmmeenntt hhaass aa ttrraannssiittiioonnaall rreeggimimee( A(Apprirlila nadndO Octcotboebr,erw, iwthitahs ah osrhtodrrt ydpryer pioedriiondA inu gAuust-). gTuhset)a. vTehrea gaevearnanguea alnraniunafal lrlabineftawlle ebnet1w9e6e0n– 21090670i–s21020070 ism 1m20a0t mthme Baét ttéhroe uBératéinrofaul lrsatiantfiaolnl , satantdio1n1,0 a0nmd m110a0t Smamvè a. tT Shaevflè.o Twhed yflnoawm dicynisacmhiacr iasc ctehrairzaecdtebryizaedh bigyh ad hisigchh adrigsechdaurrgien gduthr-e irnagin tyhes eraasionny. sTehaesomna. xTihmeu mmaflxoimwubmet wfloewen bMetawyeaennd MSeapyt eamndb eSreopvteemr tbheerp oevrieord th1e9 6p0e–r2i0o0d7 is in the o 3 31960–2007r idse irno fth2e7 0omrde/r soaft 2B7é0t émro3u/sa antd B4é8t0érmou/ asnadt t4h8e0 Smav3/ès oautt lte3h te. FSraovmè Nouotvleetm. Fbreormto NMoavye,malbmeor sttoa Mll tahye, railvmerosstd rayll uthpew riitvhearsv edrrayg euspo wf liothw aflvoewrasgoefs 5omf lo/ws aftloSwavs èoan 3 f 5d m2 3m◦ ◦ /s a /ts SaatvBèé aténrdo u2. mTh3/es aant nBuétaélrmoue.a Tnhtee manpneuraatlu mreeraann tgeemispebreattwueree nra2n4geC isa nbedtw33eeCn .24 °C and 33 °C. Our assessment is focused on the year 2010, known as the year when the Ouémé river expeOriuenr caesdsedssrmamenatt iics flfoocoudsiendg oans tahcea yseeasrt u20d1y0. ,R kenfo. w[2n] sahso twhee ydetahra wt thheenm thaex Oimuuémmév railvueers eoxfpdeirsiecnhcaergde dr 3 3 r eacmoradtiecd fldooudriinngg tahse a pcearsieo sdtu1d98y9. –R2e0f0. [92i]s shleosws tehda tnha1t4 t0h0em ma/xsimatuSmav vèa,luanesd o6f5 d0ismch/arsgaet rBeéctoérrdoeud. dAunrainlygs itsheo fpsetraiotido n19d8a9t–a2f0o0r9 tihse lepsesr tihoadn1 1946000– 2m030/7s asht Soawveèd, atnhdat 6t5h0e mpe3/ask ast oBfédtéisrcohua. rAgenat Savè (3alysis of rsetaspti.oBn édtéartoa ufo) ra rtehea bpoeurito9d1 019(6re0s–p2.00477 0s)h, o1w06e7d( rthesapt .th56e 0p),eaaknsd o1f2 d00is(crheasrpg.e6 4a0t )Smavè/s, respectively, for a 5-, 10-, and 20-year return period.and precipitation st a(trieosnp.d Baétatéurosued) airne tahbisousttu 9d1y0 w(reesrep.c 4o7ll0e)c, t1e0d67o v(erersSpa. v5 T è60 h a) e n, d dan is Bd c h1a2r0geétéro0u (resp 3outle.t s64o0f)t hme O/su, éremspéercivtievreblyas, ifno.rT ah 5e-t,h 1r0e-e,- hanoudr 2ly0-syaetealrl irteeteusrtnim paetreisoodf. tThheeT rdoipscichaalrRgae inanfadll pMreecaispuitraintigonM sistasitoionn( TdRatMa Mus,e3dB 4in2 tvh7isd setruivdeyd wdeariley caotl0le.2c5te◦dh oovriezro Snatavlèr easnodlu Btiéotnér,o1u99 o8u–tnleeatsr - opfr ethseen Ot; ureéfm. [é3 7r]i)vdear tabsaesti,na. nTdhteh ethdraeiely-hColuimrlya tesaHtealzliaterd essGtimroautpesI noffr atrheed TPrreocpiipciatal tRioaninwfaitlhl MSteaatisounrsin(gC HMIiRssPiSo;nc h(TirRpMs vM2.,0 3aBt402. 0v57◦ dhoerriivzoednt adlarielys oaltu t0i.o2n5;° 1h9o8r1i–znoenatra-lp rreessoenlut;tiroenf., [13989])8–is nuesaerd-pfroersmenotd; reelfe. v[3a7lu])a tdioanta.set, and the daily Climate Hazards Group Infrared Precipita- tion with Stations (CHIRPS; chirps v2.0 at 0.05° horizontal resolution; 1981–near-present; ref. [38]) is used for model evaluation. Water 2022, 14, 1192 5 of 23 2.2. Experimental Design of the Weather Research and Forecasting (WRF) and Coupled WRF/WRF-Hydro Model The advanced Weather Research and Forecasting (WRF, [39,40]), developed by the U.S. National Center for Atmospheric Research (NCAR) version 3.7.1, is used for both the WRF-only and fully coupled WRF/WRF-Hydro modeling in the study. The fully coupled WRF/WRF-Hydro, hereafter referred to as WRF-H, is a non-hydrostatic, mesoscale numerical weather prediction and atmospheric simulation system. Table 1 summaries the different physics schemes and experimental details. The model setup uses one domain at 5 km spatial resolution, covering the area 7◦ W–12◦ E, 0◦–18◦ N, and 400 × 400 grid points, with 30 s as numerical simulation time step. The vertical structure of the domain consists of 50 levels, up to a 10 hPa pressure top. The WRF-Hydro model is an extension package of a hydrological component to WRF. This supplementary model allows the lateral representation of the hydrological condition at the land surface [41]. The option of land use categories “Moderate Resolution Imaging Spectroradiometer (MODIS, 20 classes; ref. [42])” is selected. The Noah LSM model ref. [25] is used as the column land surface physics model. Table 1. Experimental details of the atmosphere model, WRF-only and WRF-H. Subject Option Reference Driving data European Centre for Medium-Range WeatherForecasts (ECMWF) operational analysis ECMWF Horizontal resolution 5 km Horizontal grid 400 × 400 Integration time step 30 s Projection resolution Mercator Vertical discretization 50 layers Output interval 24 h for WRF, 30 days for WRF/WRF-Hydro Simulation period 1st January 2008–31st December 2010 Pression top 10 hPa Microphysics scheme Single Moment Microphysics class 5 (WSM5) [43] Longwave radiation Rapid Radiative Transfer Model (RRTM) [44] Shortwave radiation Dudhia [45] Planetary boundary lager Asymmetric Convection Model (ACM2) [46,47] Land use Moderate Resolution Imaging Spectroradiometer(MODIS) [42] Land surface scheme Noah Land Surface Model (LSM) [25] For purposes of hydrometeorological simulations with WRF-H, the WRF domain is additionally coupled with routing processes at 500 m resolution for 4000 × 4000 grid points in the east-west and north-south directions. The fully coupled mode simulations are performed for three years, from January 2008 to December 2010, with January–February 2008 as the spin-up period. The driving data is the operational analysis dataset from the European Centre for Medium-Range Weather Forecasts (ECMWF), which provides the initial and lateral boundary conditions at 0.125◦ horizontal resolution and 6 h time intervals. Both WRF-H and WRF-only components of the coupled modeling system share the same physics parameterizations in Table 1. The WRF is the atmospheric model which solves the equations of atmospheric motions. The atmospheric state variables are the wind components, pressure, temperature, humidity and hydrometeor mixing ratios, etc. The lower boundary of the atmospheric variables is forced with a land surface model. In the case of WRF-Hydro, a complex land surface model, with an explicit description of river streamflow, is employed. The uncoupled WRF-Hydro Water 2022, 14, 1192 6 of 23 model consists of a variety of parameters (e.g., ref. [48]), which usually require calibration. Since the aim of the research is to evaluate the performance of WRF-H to simulate discharge, and therefore analyze its predicting skills about floods, the calibration is performed based on discharge at the Savè catchment outlet (Figure 1). 2.3. Calibration of WRF-Hydro in Offline Mode For calibrating the model WRF-Hydro 4.0, we focus on selected sensitive parameters highlighted in previous studies, such as REFKDT, SLOPE, RETDEPRTFAC, OVROUGHRT- FAC, and MannN. Applying a stepwise approach, following previous WRF-Hydro stud- ies [17,32,49,50], we first focus on the parameters controlling the total water volume, namely infiltration factor (REFKDT) and surface retention depth (RETDEPRT, see Table 2). It is noted that REFKDT is a tunable parameter that significantly impacts surface infiltration, and hence the partitioning of total runoff into the surface and subsurface runoff; increasing REFKDT decreases surface runoff. Since there is not a historical range to estimate these parameters over the interested domain, the study calculates them from 0.1 to 10, with 0.1 increments. The second step of the calibration is to evaluate the coefficient governing deep drainage (SLOPE); the same method used in the case of REFKDT and RETDEPRTFAC for selecting the optimum value is applied, by testing values from 0.1 to 1.0, with 0.1 as the increment. The adjustment of the roughness parameter, which controls the overland flow, is performed from the default value to the optimum one. Table 2. Validation of the performance of simulated WRF-Hydro and observed discharge at Savè, with parameters such as the infiltration-runoff parameter REFKDT, the retention factor RETDEPRTFAC, the SLOPE, the overland flow roughness scaling factor OVROUGHTFAC, and the Manning’s roughness coefficients, MannN. Experiments in bold show the optimum parameters based on the objective criteria (Corr, NSE, and KGE). REFKDT Range 0.1 0.5 0.8 1.0 1.5 2.0 3.0 3.5 4.5 NSE 0.29 0.33 0.38 0.46 0.52 0.34 0.20 0.17 0.11 Corr 0.64 0.67 0.63 0.63 0.58 0.60 0.61 0.64 0.60 KGE 0.38 0.39 0.41 0.47 0.49 0.32 0.09 0.07 0.07 RETDEPRTFAC Range 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 NSE 0.50 0.52 0.49 0.49 0.49 0.49 0.49 0.49 0.49 Corr 0.58 0.58 0.59 0.59 0.59 0.59 0.59 0.58 0.58 KGE 0.48 0.49 0.49 0.49 0.48 0.46 0.47 0.46 0.46 SLOPE Range 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 NSE 0.52 0.61 0.60 0.50 0.37 0.20 0.12 −0.03 −0.23 Corr 0.58 0.62 0.65 0.65 0.66 0.61 0.63 0.55 0.61 KGE 0.49 0.56 0.52 0.39 0.33 0.24 0.19 0.11 0.04 OVROUGHRTFAC Range 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 NSE 0.60 0.60 0.60 0.61 0.63 0.62 0.62 0.61 0.61 Corr 0.67 0.65 0.61 0.60 0.66 0.66 0.64 0.64 0.62 KGE 0.56 0.56 0.57 0.60 0.60 0.60 0.59 0.56 0.56 MannN Range 0.2 0.4 0.6 0.8 1.0 1.4 1.6 1.8 2.0 NSE 0.42 0.44 0.50 0.56 0.61 0.62 0.66 0.62 0.60 Corr 0.66 0.63 0.73 0.68 0.70 0.72 0.67 0.69 0.67 KGE 0.36 0.37 0.41 0.48 0.50 0.57 0.63 0.60 0.61 Water 2022, 14, 1192 7 of 23 In addition, sensitivity tests are done on the surface (OVROUGHRTFAC) and chan- nel roughness parameter (MannN), which controls the shape of the hydrograph. The efficiency criteria—Kling–Gupta efficiency (KGE), Correlation coefficient (Corr), Nash- Sutcliffe Efficiency (NSE) and BIAS are used to evaluate the model performance within the calibration process. √[ n( xy)− ( x)( y)Corr ∑ ]∑= [ ∑ ] (1) n ∑ x2 − (∑ x)2 n ∑ y2 − (∑ y)2  ( )  ∑n 2 obs sim i=1 Yi −Yi NSE 1− =  ( )  (2) ∑n Y obs −Y mean 2√ i=1 i i( ) ( ) − − 2 σ 2 2 KGE = 1 (r 1) + sim( − µ 1 + sim − 1 (3) σobs ) µobs ∑n Y sim −Y obsi=1 i i BIAS = × 100 (4) ∑n obsi=1 Yi where x is the observations and y is the simulations; Y obsi is the ith observation for the constituent being evaluated, Y simi is the ith simulated value for the constituent being evalu- ated, Y meani is the mean of observed data for the constituent being evaluated, and n is the total number of observations [51]. For the KGE, r is the linear correlation between observa- tions and simulations, σobs is the standard deviation in observations, σsim is the standard deviation in simulations, µsim is the simulation mean, and µobs is the observation mean. In order to harmonize the uncoupled and coupled setups, the uncoupled simulations use the same time step as the WRF-only and WRF-H simulations (30 s). The calibration of the model is performed using hourly dataset input, and the focus is on the performance skill in reproducing daily discharge in the sub-catchments. One-year calibration is considered sufficient to evaluate the basic parameter sensitivities [32]. WRF-Hydro is calibrated on period P1 (2008) and validated on period P2 (2009–2010), where P1 and P2 are the shared periods containing the whole study period, named P (2008–2010). 2.4. Evaluation of Model Uncertainty with the Stochastic Kinetic Energy Backscatter Scheme The stochastic kinetic energy backscatter scheme (SKEBS; [35,52,53]), which primarily acts on the dynamical tendencies at the lateral boundaries, is activated into WRF-H for the fully coupled simulation (WRF-H-SKEBS). The SKEBS technique provides several advantages over perturbation techniques that only perturb the initial state. The method aims to represent model uncertainties associated with scale interactions that take place in the real atmosphere but are absent in a truncated numerical model [54]. SKEBS perturbs the model fields by adding random and amplitude perturbations (noise) to the horizontal wind and potential temperature tendency equations at the lateral boundaries for each time step [55]. An ensemble of 10 members using the WRF-H-SKEBS model is generated for each rain season of the period P (2008–2010). 3. Results and Discussion 3.1. Evaluation of WRF-Only Precipitation The WRF model (no coupled, referee as WRF-only) is run over the domain displayed in Figure 1 in order to generate atmospheric input data for the calibration of the supplement model WRF-Hydro. Figure 2 shows the comparison between the simulated WRF-only and observed dataset. The weekly WRF precipitation for the Savè catchment is relatively close to that derived from CHIRPS (Figure 2a) and TRMM (Figure 2b), with the mean coefficient Water 2022, 14, x FOR PEER REVIEW 8 of 23 3. Results and Discussion 3.1. Evaluation of WRF-Only Precipitation The WRF model (no coupled, referee as WRF-only) is run over the domain displayed in Figure 1 in order to generate atmospheric input data for the calibration of the supple- ment model WRF-Hydro. Figure 2 shows the comparison between the simulated WRF- only and observed dataset. The weekly WRF precipitation for the Savè catchment is rela- tively close to that derived from CHIRPS (Figure 2a) and TRMM (Figure 2b), with the mean coefficient of determination (R2) equal to 0.64, and 0.59, respectively. The agreement between the two observed datasets (CHIRPS and TRMM) is about 0.87 for R2. Since the investigation is on the potential of WRF-H (model WRF-Hydro coupled with WRF) for flood prediction, the available WRF-only precipitation at the highest spa- tiotemporal resolution is used to force the uncoupled WRF-Hydro model. In particular, the hourly output of WRF at 5 km spatial resolution is used as meteorological forcing data, containing necessary variables such as incoming shortwave radiation (W/m2), incoming longwave radiation (W/m2), specific humidity (kg/kg), air temperature (K), surface pres- sure (Pa), u and v components of near-surface wind (m/s), and liquid water precipitation rate (mm/s). The meteorological forcing data needed by the Noah LSM (land surface hy- drological modeling system) are prepared as hourly gridded data. The Noah LSM static Water 2022, 14, 1192 data (topography, land cover, soil type) are too coarse for a WRF-Hydro applica8toiof 2n3. Therefore, datasets from the Shuttle Elevation Derivatives at Multiple Scales (Hy- droSHEDS) database [56], high-resolution topography and channel network) are consid- oerfedde tteor amccinuartaitoenly( Rro2u) teeq wuaaltteor a0c.6r4o,sas nthde0 l.a5n9d, rsecsappeec tthivreoluy.gThh oevaegrlraenedm, esnutbbsuetrwfaeceen, othr echtwano- onbesl eflrovwed. datasets (CHIRPS and TRMM) is about 0.87 for R2. Fiigure 2.. SSccaattteter rpploltost sofo sfimsimuluatleadte Wd RWFR-oFn-olyn lvys.v osb. seorbvsedrv (ead) th(ae) CthliemCatleim HaatezaHrdasz GarrdosupG rIonufrparIend- fPrraerceidpPitraetcioipni twatitohn SwtaitthioSntsa t(ioCnHsI(RCPHSI)R, PaSn)d, a(nbd) (tbh)e thTeroTproicpailc aRlaRinaifnaflla lMl Meaesausurirningg MMisissioionn ((TRMM) weeeeklly ccuummuulalatitviveep prerceicpiiptaittaiotinonfo fropre prieordio2d0 0280–0280–1200,1(0c,) c(co)m cpoamripsaornisboentw beeetwn oeebnse orvbesder,vaend, (adn)da n(adl)- yasniasloyfsisse ovfe nse-dveany-fidlatyer feidltedraeidly dparileyc ipprietactiipointaftoiornW foRrF W-onRlFy-,oCnHlyI,R CPHS,IRanPdS,T aRnMd MTR.MM. 3.2. CSainlibceratthioeni navneds Etivgaltuioatnioins ofn WthReFp-oHtyednrtoia Ol offliWneR F-H (model WRF-Hydro coupled with WRF)Tfhoer Nfloooadh pLrSeMdi cctoiuonp,lethde waivthai WlabRleF-WHRydF-roon elvyaplureactiepsi tsautriofancea twthaetehr idgehpestht scphaatniogtee mh(pmo)- raasl ar ersaotelu otvioenr itsimues eodf tionffioltrrcaetitohne euxncceosus p[2le5d]: WRF-Hydro model. In particular, the hourly output of WRF at 5 km spatial resolution is used as meteorological forcing data, containing necessary variables such as incoming shortwave radiation (W/m2), incoming longwave radiation (W/m2), specific humidity (kg/kg), air temperature (K), surface pressure (Pa), u and v components of near-surface wind (m/s), and liquid water precipitation rate (mm/s). The meteorological forcing data needed by the Noah LSM (land surface hydrological mod- eling system) are prepared as hourly gridded data. The Noah LSM static data (topography, land cover, soil type) are too coarse for a WRF-Hydro application. Therefore, datasets from the Shuttle Elevation Derivatives at Multiple Scales (HydroSHEDS) database [56], high-resolution topography and channel network) are considered to accurately route water across the landscape through overland, subsurface, or channel flow. 3.2. Calibration and Evaluation of WRF-Hydro Offline The Noah LSM coupled with WRF-Hydro evaluates surface water depth change h(m) as a rate over time ofinfiltratio[n excess [25]: ][ ( )] ∂h ∂P  ∑[ 4 i=1 ∆Zi(θs − θi) 1][− exp −(k Ks δtd Kre f 86400 )]= 1− (5) ∂t ∂t P 4d + ∑i=1 ∆Zi(θs − θi) 1− exp −k Ks δt  Kre f 86400 where Pd(m) is the precipitation that is not intercepted by the canopy; ∆Zi(m) is the soil layer i in Noah LSM, θi is the volumetric water content (soil moisture) of the soil layer i, θs is the saturated soil moisture, and the saturated hydraulic conductivity is Ks (m/s); both θs Water 2022, 14, 1192 9 of 23 and Ks varied in terms of the soil texture. The reference saturated hydraulic conductivity Kre f is equal to 2× 10−6(m/s). The calibration of the model (WRF-Hydro) is performed in offline mode for the Savè catchment with parameters REFKDT, SLOPE, RETDEPRTFAC, OVROUGHRTFAC, and MannN, using both KGE and Corr as efficiency criteria. Efficiency criteria results of the optimization of those parameters are listed in Table 2. The first parameter evaluated is one of the parameters controlling the total water volume (the REFKDT), whose values range from 0.1 to 10.0, with 0.1 increments. The default simulation (REFKDT = 3.0) realized had shown an underestimation of the observed streamflow because the REFKDT controls the infiltration of the surface water, and its value shall be reduced to disable many infiltrations. As illustrated with [49] for the case of the Sissili in West Africa, and [48] in Kenya (East Africa), we found that the model discharge performance is highly sensitive to parameter REFKDT. For the case in study, the REFKDT = 1.5 performs better than using the default value (3.0) with statistics (NSE = 0.52, KGE = 0.49, and Corr = 0.58). The second parameter assessed is the RETDEPRTFAC, which adjusts the initial retention depth in the model. Scaling factors (RETDEPRTFAC) between 0 and 10 with 1.0 increments show that the modeled discharge remains the same regardless of the change in RETDEPRTFAC, as the scores are very close to each other. As mentioned by [41], increases in the RETDEPRTFAC on channel pixels can encourage more local infiltration near the river channel leading to wetter soils. Therefore, the default RETDEPRTFAC = 1.0 seems to provide better results and will be considered such for the next calibration steps. In the case of the present study, the uncoupled model WRF-Hydro is also very sensitive to the parameter SLOPE. The optimized value of the parameter SLOPE is 0.2. The parameters controlling the hydrograph shape are also investigated. As illustrated by [50], the surface roughness (OVROUGHRT) plays an important role in transmitting infiltration excess water to channel networks and is calibrated in WRF-Hydro using a scaling factor (OVROUGHRTFAC) between 0.2 and 1.0, with 0.1 increments. Considering the correlation coefficient (Corr), the NSE and the KGE statistics, and the hydrograph shape match between simulated and observed hydrographs at Savè, the scaling factor value of OVROUGHRTFAC = 0.6 is judged as the best to fit the simulated hydrograph to the observed hydrograph. Furthermore, the calibrated Manning’s coefficient (MannN) for the river channel routing used are set as 1.75 for stream order 1, 1.70 for stream order 2, 1.65 for stream order 3, 1.60 for stream order 4, and 1.55 for stream order 5. It can be seen that the Water 2022, 14, x FOR PEER REVIEW observed discharge hydrograph at Savè is reasonably well reproduced, with KGE = 0.6310, of 23 and Corr = 0.67, for the calibration period March-December 2008 (Figure 3a). FigFuirgeu r3e. O3.bsOerbvsedrv aednda nsidmsuimlauteladt e(dun(cuonucopulepdle Wd RWFR-FH-Hydyrdor)o )dadialyily hyhdydroroggrraapphhss aatt SSaavvèè,, aanndd catch- mecnatt-cahvmereangt-eadv edragileyd pdraeilcyipirteactiipointa tdioenridverdiv ferdofmro mWWRFR-Fo-nolnyly: :(a(a) )ccaalliibrationp peeriroiodd2 020080; 8(b; )(bva) lvi-alida- tiond apteiorniopde r2i0o0d92–020091–02.0 10. The robustness of the calibrated model in the uncoupled module is also evaluated over the Bétérou catchment (an inner-domain of the Savè basin). Figure 4 shows the best- fit performance of the model in simulating the shape of the observed hydrograph at Bété- rou. It reproduces the peaks of discharge relatively well, and the trend of the shape of the observed hydrograph is also well reproduced. The performances of the model are re- vealed by the statistics KGE = 0.74, Corr = 0.85, and NSE = 0.68 for the calibration period over the Bétérou catchment for the uncoupled module simulation. Figure 4. Evaluation of the uncoupled calibrated WRF-Hydro model over the Bétérou catchment during the calibrated year (2008). 3.3. Evaluation of the Coupled Model WRF-H The calibrated (optimized) model parameters are used for the online WRF-Hydro model (referred to as WRF-H), to assess the performance of the calibrated WRF-H to sim- ulate discharge and precipitation in the study area. Water 2022, 14, x FOR PEER REVIEW 10 of 23 Water 2022, 14, 1192 10 of 23 The calibrated model is evaluated offline for the period P2: 2009–2010 (Figure 3b). The above-mentioned efficiency criteria allow us to evaluate the performance of the model. It can be noticed that it simulates the trend and peaks of the observed discharge relatively well, even slightly better in comparison to the calibration period, with model efficiencies KGE of 0.86 and Corr of 0.87. This enhanced performance for the validation period is related to theFmiguucrhe 3h.i gOhbesrerdviesdc haanrdg seimpeualakteind 2(u0n1c0o, ui.pel.,edth We flRoFo-Hdyindgroy) edaari,lyw hhyicdhroigsrraeplahtsi vate lSyavè, and catch- well reprodumceednt-bayvetrhaegemd oddaielyl epvreacliupaittaetdiono fdfleirnivee. dG frloomba WllyR, Ff-oornltyh:e (as)i mcauliblartaitoionn ppeerriioodd 2P008; (b) valida- (2008–2010), WtioRnH p-eHrioydd r2o00in9–o2f0fl1i0n. e mode is able to simulate discharge, with KGE and Corr equal to 0.70 and 0.74. The robustneTshseo frothbuesctanleibssr aotef dthme ocadleibl rinatethde muondcoeul pinl etdhem uondcuoluepilseadl smooedvualleu aiste adlso evaluated over the Bétéroovuerc athtceh Bméetnérto(aun ciantcnhemr-deonmt (aain ionfntehre-dSoamvèabina soinf )t.hFei gSauvreè 4bashsionw). sFtihgeurbee s4t -sfihtows the best- performance foitf ptheerfmoromdaenl cine soifm thuela mtinogdetlh iens shiampuelaotfinthge tohbes sehravpede ohfy tdhreo ogbrasperhvaetdB héytédrrooug.raph at Bété- It reproducesrotuh.e Itp reeapkrsodofucdeiss cthhaer pgeearkesl aotfi vdeislychwareglle, raenladtitvheelyt rweneldl, oanf dth tehes htraepnedo off tthhee shape of the observed hydorbosgerravpehdi shayldsorowgreallprhe pirso adlusoce dw.eTllh erepperrofdorumceadn.c eTshoef ptheerfmoromdaenl caeres roefv ethaele dmodel are re- by the statistivcesaKleGdE b=y 0t.h7e4 ,sCtaotirsrti=cs0 .K85G,Ea n=d 0N.7S4E, C=o0rr.6 =8 0fo.8r5t,h aencda lNibSraEt i=o n0.6p8e rfioord thove ecratlhiberation period Bétérou catchomveern tthfoe rBtéhteéruonuc ocautpchlemdemnot dfourl ethseim unuclaotuiopnle. d module simulation. Figure 4. EvalFuiagtuiorne o4f. Ethvealuunactoiounp olefd thcea liubnrcaoteudplWedR cFa-lHibyrdatreodm WodReFl-Hovyedrroth me Bodéteélr ouvecra thche mBeéntétrou catchment during the calidburartiendg ytheea rca(2li0b0r8a)t.ed year (2008). 3.3. Evaluation3.3of. EthveaCluoautpiolned oMf tohde eCl oWuRplFed-H Model WRF-H The calibrateTdh(eo pcatilmibrizaetedd) m(oopdtieml ipzaerda)m meotedresl apraeruamseedtefrosr athree uosneldin efoWr tRhFe -Honyldinreo WRF-Hydro model (referrmedodteol (arsefWerRreFd- Hto) ,atso WaRssFe-sHs)t, htoe apsesrefsosr tmhae npceerfoofrmthaenccea loibf rtahtee dcaWlibRraFt-eHd WtoRF-H to sim- simulate discuhlaartgee dainscdhparrgecei panitda tpiorneciinpitthaetisotnu diny tahree as.tudy area. 3.3.1. Precipitati on Simulations The agreement skills of WRF-H are evaluated for precipitation both temporally and spatially for the rainy season of the period 2008–2010. Figure 5 exposes a comparison between the weekly precipitations from WRF-only, WRF-H, and observed datasets. The coefficient of determination between the WRF-H and WRF-only, is equal to 0.88 (Figure 5a). This shows clearly that WRF-only and WRF-H simulate precipitation differently, which has been already illustrated by [17,32,57,58]. Figure 5b,c compare WRF-H with CHIRPS and TRMM, and it illustrates a good agreement between these datasets. The slightly Water 2022, 14, x FOR PEER REVIEW 11 of 23 3.3.1. Precipitation Simulations The agreement skills of WRF-H are evaluated for precipitation both temporally and spatially for the rainy season of the period 2008–2010. Figure 5 exposes a comparison be- tween the weekly precipitations from WRF-only, WRF-H, and observed datasets. The co- efficient of determination between the WRF-H and WRF-only, is equal to 0.88 (Figure 5a). Water 2022, 14, 1192 This shows clearly that WRF-only and WRF-H simulate precipitation differently, 1w1hofic2h3 has been already illustrated by [17,32,57,58]. Figure 5b,c compare WRF-H with CHIRPS and TRMM, and it illustrates a good agreement between these datasets. The slightly better agreement of CHIRPS (compared to TRMM) with WRF-only (Figure 2b) and WRF-H (Fig- buertet e5rba) gcroeuelmd ebnet eoxfpClaHinIRedP Sby(c tohme phaigrehd retosoTluRtMioMn o) fw bioththW CRHFI-RoPnSly a(nFdig WurReF2-bH) apnrdecWipRitaF-- Htio(nF. iAgusr see5enb) pcroeuvlidoubselye (xFpiglauinree d2db)y, CthHeIRhPigSh arneds oTluRtMioMn opfrebcoitphitaCtHioInRsP’ Shaavned thWe RsaFm-He ptrreencidp,i otantliyo nd.ifAfesrs iene nthperire vmioaugsnliytu(Fdiegsu; rtehe2 dW),RCFH-oInRlPyS daenpdicTtsR wMeMll tphree coibpsiteartvioatniso’nhsa fvoer tthhee swahmoelet rseenadso, no.n Ilnyddeiefdfe,r thine tchoemirpmaraisgonnit oufd Feisg;uthreesW 2 RanFd-o 5n illyludsetrpaitcetss twhealtl WthReFo-bHs eprevraftoiromnss fsolirgthhtelyw bhetotleer stehaasno Wn.RInFd-oeneldy, itnh teecrmoms pofa rwiseoenkloyf pFriegcuipreitsa2tioannsd. O5 vilelurasltlr,a Fteigsutrhea 5t dW sRhFo-wHs pthearfto WrmRsFs-lHig hftollylobwetter than WRF-only in terms of weekly precipitations. Overall, Figushows that WRF-H fos llwowelsl wine tlleirnmtse romf sthoef tthreentrde,n dth,eth We WRFR-Fo-nolnyl.y A. Adddditiitoionnaalllyly, ,FFiigguur r re e e 55d5dd eennhhaanncceess tthheessee rreessuullttss wwiitthh CCoorrrr == 00..6688 bbeettwweeeenn WWRRFF--HH aanndd CCHHIIRRPPSS,, aaggaaiinnsstt CCoorrrr == 00..5599 bbeettwweeeenn WWRRFF--HH aanndd TTRRMMMM. D. eDsepsiptei ttehet hgeoogdo opderfpoerrmfoarnmcea nocf emoofdmelso,d thelesy, pthreeysenptr essoemnet sinocmeretiitnucdeer twituhdeerewvehre trheevye ratrhe eaypaprleiedap. pRleife.d [.59R]e sfh. o[5w9e]ds htohwate tdhet hhaitghth veahriigabhilvitayr ioafb iplirtey- ocifppitraetciiopnit iant iWonesint AWferisctaA rfersiucaltsr efsroumlts af rloamrgea ulanrcgeertuanincteyr tiani nWtyRFin sWimRuFlastiimonusl.a Ttihoins su.nTcheirs- utanincetyrt aisi nintyveisstiignavteesdt iigna dteedtaiinl ind ethtaei lcainset hofe WcaRsFe-oHf iWn SReFc-tHionin 3.S4e bcyti omno3d.i4fybinygm booduinfydianrgy bcoonudnidtiaornysc wonitdhi taio sntoscwhiatshtiac spteorcthuarsbtaictiponer. turbation. FFiigguurree 55.. EEvvaalluuaattiioonno offw weeekelkylys ismimulualtaetdedpr pecriepciitpaittiaotniosnwsi twhiothb soebrvseerdvdedat adsaettassaettsth aet Stahveè Scaavtcèh cmatecnht-: (ma)eWnt:R (Fa-)H WvRsF. -WHR vFs-.o WnlRy,F(-bo)nWly,R (Fb-)H WvRs.FC-HH vIRs.P CS,H(cIR) WPSR, (Fc-)H WvRs.FT-HR MvsM. T, RanMdM(d, )aannda (ldy)s iasnoaflysesvise onf- dseavyefinl-tdeareyd fidltaeirleydp dreaciilpy iptarteicoinpiftoartiWonR fFo-rH W, WRFR-FH-o, WnlyR,FC-HonIRlyP, SC,HanIRdPTSR, ManMd .TRMM. TThhee moonntthhllyy ttrreenndd ooff pprreecciippiittaattiioonn dduurriinngg tthhee Weesstt AAffrriiccaann Moonnssoooonn ((WAAM)) ppeerriioodd ((JJuunnee––SSeepptteemmbbeerr:: JJJJAASS)) hhaass bbeeeenn iinnvveessttiiggaatteedd aatt tthhee OOuuéémméé––SSaavvèè oouuttlleett ((SSaavvèè)).. IItt sshhoowwss tthhaatt tthhee mmoonntthhllyyc cuummuulalatitviveet rternenddf oflolollwows so bosbesrevravtaiotinosnssh sahpaepbeu btuvta vriaersieins min amgnaigtunditeu,daes, ialslu isllturastteradteind Finig Fuirgeusr2esa n2 dan5d. F5i.g Fuirgeu6res h6 oswhoswthse thmea mxiamxuimmummo mntohnltyhclyu mcuumlautilvateivraei rnafianll- afatlSl aavt èSaivnèA inu gAuusgt,uastn, danwdh wilheicleo mcopmaprianrgintgh tehteh trhereeep aptatettrenrsn,si, tits hshoowwsst thhaatt tthhee mmoonntthhllyy ccuummuullaattiivvee pprreecciippiittaattiioonn rreeccoorrddss dduurriinngg tthhee ttwwoo llaasstt mmoonntthhss ((AAuugguusstt aanndd SSeepptteemmbbeerr)) iinn 2010 are higher than both 2008 and 2009, expressing the heavy rainfall recorded in this particular year. To observe more detail of the distribution of precipitation over the focusing domain for either WRF-only and WRF-H, as well as CHIRPS, Figure 7 presents the spatial distribution of precipitation during the WAM (JJAS) of the flooding year (2010). Water 2022, 14, x FOR PEER REVIEW 12 of 23 2010 are higher than both 2008 and 2009, expressing the heavy rainfall recorded in this Water 2022, 14, 1192 particular year. To observe more detail of the distribution of precipitation ov1e2r othf e23 focus- ing domain for either WRF-only and WRF-H, as well as CHIRPS, Figure 7 presents the spatial distribution of precipitation during the WAM (JJAS) of the flooding year (2010). Water 2022, 14, x FOR PEER REVIEW Figure 6. PrecipitaFtiiognurse t6r.e PnrdeciipnitaJtuionnes –trSeenpdt einm Jubneer–S(eJJpAteSm)baert (tJJhAeS)S atv thèe cSaatvcèh cmatechnmt efnotr fotrh teh 1ep3 p eoerrfi io2d3 2008– 2008–2010. 2010. Figure 7a,b depict the WRF-only simulations and WRF-H precipitations, respec- tively, while the difference between the two models in Figure 7c shows that WRF-H either underestimated or overestimated simulation in comparison to WRF-only, depending on the location (as illustrated [31]). Figure 7d depicts the spatial distribution of CHIRPS over the focusing domain; its color bar is as in Figure 7a,b. The mean precipitation for the do- main is about 864 mm for WRF-only (Figure 7a), and 947 mm for WRF-H (Figure 7b). This means that WRF-H increases the simulated precipitation from WRF-only by about 1%. The observed precipitation in this domain is about 817 mm (i.e., less than simulated pre- cipitation from both WRF-only and WRF-H). Similar results are obtained for the Savè catchment, with seasonal spatial-averaged precipitation of about 1049 mm for WRF-H, 998 mm for WRF-only, and 977 mm for CHIRPS. As the model performed globally well precipitations both for Savè and Bétérou at temporal scale and spatial distribution, the following analysis focused on those two basins of the Oueme river by investigating the temporal cumulative variability of precipitation. Therefore, Figures 8 and S1 present, respectively, the cumulative precipitations derived from WRF-H, and the satellite datasets CHIRPS and TRMM at Beterou and Savè. They attempt to assess WRF-H behavior about the cumulative precipitation over the study pe- riod (2008–2010). A very good reproductivity of the observed dataset by WRF-H is no- ticed, and also the clear rainy season period is well captured by the model. It can also be noticed that the total seasonal precipitation is based on the rainfall from April to October over the basins. The cumulative seasonal amount (between April and October) is indica- tive of the average of annual precipitation for a given year in the region. Figure 7. JJJJAS prreciipiitations for the ffllooding year 2010 in the Savè catchment (red contour), and the Bétérou catchment (purplle contour):: ((a)) WRF--onlly siimullattiionss;; ((b)) couplleed WRF--Hydrroo ((WRF--H)) ssiimuullaattiioonnss;; ((cc)) ddiiffffeerreenncceeb beetwtweeeennW WRRFF-H-Ha nadndW WRFRF(W (WRFR-FH-Hm minuinsuWs WRFR);Fa)n; adn(dd )(dC)H CIRHPIRS PpSre pcirpei-- ctaiptiiotant.ioTnh.e Tchoelo crobloarr boafrF oigf uFriegu5rde i5sdu sise dusaesda acso am cmomonmcoonlo crobloarr bfoarr Ffoigr uFriegu5are,b 5,da.,b,d. According to Figures 8 and S1, we can make the assumption that there is no or neg- ligible rainfall recorded in the dry season. At Bétérou (Figure 8), the model underestimates the cumulative precipitation compared to TRMM and CHIRPS, albeit showing a good capture of seasonal variability over the years. The respective total three-yearly amounts are WRF-H = 3921 mm, TRMM = 4077 mm, and CHIRPS = 4000 mm, which confirms the consistency with the results discussed above that the WRF-H precipitation simulations are closer to the observed CHIRPS than TRMM. The bias of the cumulative precipitation is between −8.82% and 0.00% in regard to CHIRPS, and between −8.16% and 2.13% for TRMM. At Savè (with WRF-H = 3883 mm), conversely to analysis at Bétérou, the cumula- tive precipitation of the model fits well with TRMM (that value is about 3985 mm) rather than CHIRPS (value equal to 4020 mm). The error related to the evaluation of the cumu- lative precipitation is between −5.82% and 2.18% for TRMM, and about −7.00% and 0.00% for CHIRPS. Water 2022, 14, 1192 13 of 23 Figure 7a,b depict the WRF-only simulations and WRF-H precipitations, respectively, while the difference between the two models in Figure 7c shows that WRF-H either un- derestimated or overestimated simulation in comparison to WRF-only, depending on the location (as illustrated [31]). Figure 7d depicts the spatial distribution of CHIRPS over the focusing domain; its color bar is as in Figure 7a,b. The mean precipitation for the domain is about 864 mm for WRF-only (Figure 7a), and 947 mm for WRF-H (Figure 7b). This means that WRF-H increases the simulated precipitation from WRF-only by about 1%. The observed precipitation in this domain is about 817 mm (i.e., less than simulated precipitation from both WRF-only and WRF-H). Similar results are obtained for the Savè catchment, with seasonal spatial-averaged precipitation of about 1049 mm for WRF-H, 998 mm for WRF-only, and 977 mm for CHIRPS. As the model performed globally well precipitations both for Savè and Bétérou at temporal scale and spatial distribution, the following analysis focused on those two basins of the Oueme river by investigating the temporal cumulative variability of precipitation. Therefore, Figure 8 and Figure S1 present, respectively, the cumulative precipitations derived from WRF-H, and the satellite datasets CHIRPS and TRMM at Beterou and Savè. They attempt to assess WRF-H behavior about the cumulative precipitation over the study period (2008–2010). A very good reproductivity of the observed dataset by WRF-H is noticed, and also the clear rainy season period is well captured by the model. It can also be Water 2022, 14, x FOR PEER REVIEnWo ticed that the total seasonal precipitation is based on the rainfall from April to Oct1o4b oefr 23 over the basins. The cumulative seasonal amount (between April and October) is indicative of the average of annual precipitation for a given year in the region. FigFuigreur8e. 8T.h Tehceu cmuumlautliavteivteo ttaoltaplr epcriepciitpaittiaotniodne drievreivdefdr ofmromTR TMRMMMan adnCdH CIHRPIRSP, Sa,n adnsdim siumlautleadteidn in WRWFR-HF-oHv oerveBré tBéértoéuroduu driunrginthg ethpee rpioerdio2d0 0280–0280–1200.10. 3.3A.2c. cDoirsdchinagrgteo SFimiguulraeti8onasn d Figure S1, we can make the assumption that there is no or negligible rainfall recorded in the dry season. At Bétérou (Figure 8), the model undereDstiismchaatergset hreesucultms autl aStaivvèe aprree dciipspitlaatyioend icno Fmigpuarree d9at,o shToRwMinMg tahned dCaiHlyI RtimPSe, saelrbieeist of shsoiwmiunlgataedg o(gordeecnap) atunrde oobfsseeravseodn (arlevda) rsitarbeialmity doisvcehratrhgeesy,e aanrds. rTehlaeterdes WpeRcFti-vHe ptorteaclipthitraeteio- n ye(abrlluyea)m fooru tnhtes apreeriWodR F20-H08=–23091201. mm, TRMM = 4077 mm, and CHIRPS = 4000 mm, which confirms the consistency with the results discussed above that the WRF-H precipitation simulations are closer to the observed CHIRPS than TRMM. The bias of the cumulative precipitation is between −8.82% and 0.00% in regard to CHIRPS, and between −8.16% and 2.13% for TRMM. At Savè (with WRF-H = 3883 mm), conversely to analysis at Bétérou, the cumulative precipitation of the model fits well with TRMM (that value is about 3985 mm) rather than CHIRPS (value equal to 4020 mm). The error related to the evaluation of the cumulative precipitation is between −5.82% and 2.18% for TRMM, and about −7.00% and 0.00% for CHIRPS. Figure 9. Observed and simulated (fully coupled WRF-Hydro) hydrographs and derived precipita- tion from WRF-Hydro: (a) full year at Savè, and (b) at Bétérou (Savè’s inner catchment). Water 2022, 14, x FOR PEER REVIEW 14 of 23 Water 2022, 14, 1192 14 of 23 Figure 8. The cumulative total precipitation derived from TRMM and CHIRPS, and simulated in WRF-H over Bétérou during the period 2008–2010. 33.3.3.2.2. .D DisiscchhaarrggeeS Simimuulalatitoionnss DDisiscchhaargrgeer eresusultlstsa attS Saavvèèa arered disispplalayyeeddi ninF Figiguurere9 9aa, ,s shhoowwininggt htheed daailiylyt itmimees seerireiesso off ssimimuulalatetedd( g(grereeenn) )a annddo obbsseerrvveedd( r(reedd))s strtreeaammd disiscchhaargrgeess, ,a annddr reelalateteddW WRRFF-H-Hp prreeccipipitiatatitoionn (b(bluluee) )f oforrt htheep peeriroiodd2 2000088––22001100. . FFigiguurere9 .9.O Obbsesrevrvededa nadndsi smimulualtaetded(f u(flulyllcyo cuopulepdleWd RWFR-HF-yHdryod)rhoy) dhryodgrroagprhaspahnsd adnedr idveerdivperdec pipreitcaitpioitna- frtoiomn WfroRmF -WHyRdFr-oH: y(ad)rofu: l(lay) efuarlla yteSaarv aèt, Saanvdè(, ba)nadt (Bbé)t éart oBuét(éSraovuè ’(sSianvnèe’sr icnantcehr mcaetncth)m. ent). A good agreement can be seen between the observed and the simulated hydrographs, and an approximate good representation of the peaks of discharge, as well as hydrograph shapes, as quantified by the performance measures KGE and Corr, equal to 0.76 and 0.84, respectively. This better performance, in comparison to the offline simulation, could be explained by the time step of the meteorological data in the fully coupled mode, which is 30 s and not hourly as in offline mode. Since one objective of the study is to evaluate the performance of WRF-H to simulate the discharge, and therefore to predict potential floods, here we focus on the ability of the model in reproducing only the rainy seasons. We obtain for Savè, from Figure 9, a KGE equal to 0.22, 0.64, and 0.80 for the rainy seasons of 2008, 2009, and 2010, respectively, which gives solid information about the model’s simulation skills. It is noted that the model has a better performance in 2010. The robustness of the calibrated WRF-H over Savè is evaluated in a second catchment, i.e., in Bétérou (Savè’s inner-catchment). Figure 9b shows that WRF-H here reproduces well the discharge trend, as well as the peaks, so that WRF-H can also be used successfully for this inner-catchment. Table S1 illustrates the discharges peaks obtain for the two basins during the three years. An evaluation from a recent work of ref. [2] over the same study domain revealed that the averages of annual rainfall between 1960–2007 are 1205 mm at the Bétérou rainfall station Water 2022, 14, 1192 15 of 23 and 1098 mm at Savè. The dynamics of the flow are characterized by a high discharge during the rainy season. The maximum flow between May and September over the period 1960–2007 is in the order of 240 m3/s to 740 m3/s at Bétérou, and 1000 m3/s to 1750 m3/s at the Savè outlet. From November to May, almost all the rivers dry up, and the averages of low flows are about 5 m3/s at Savè and 2 m3/s at Bétérou. WRF-H is able to capture the flood event which occurred in September–October 2010 over Savè, as well as over Bétérou. In particular, although the predicted highest discharge peak occurs earlier than in the observation at Savè and Bétérou. The second “weak” peak in 2010, which could amplify damage intensities of the flood in the study area, is also well reproduced. According to results from Figure 9, this second peak noticed for the discharge, should result from the highest precipitation simulated and observed in September 2010. The first important peak at Savè in 2010 is also reflected from the highest simulated precipitation of August 2010. Figures 10 and 11, respectively, over Bétérou and Savè show, on the one hand, the simulated and observed discharge at daily time step for the whole WRF-H evaluation period (2008–2010: Figure 10a,b; Figure 11a,b), and the flooding period (2010: Figure 10c; Figure 11c) based on linear regression with values of line 1:1 and coefficient of determination (R2); and on the other hand, the cumulative discharge of simulated and observed datasets. Figure 11a indicates that, at Bétérou, the WRF-H underestimates the low flows, whilst in contrast, at Savè (Figure 11a), it extracts the low flow relatively well. This underestimation of the low flow at Bétérou and Savè is related to the fact that in early August at the beginning of the rainy season over the sub-basins, the model started producing streamflow (2008 and Water 2022, 14, x FOR PEER REVIEW 2009 on Figure 10a,b). Figure 9a justifies that the underestimation of the low flow at16S aofv è23 (Figure 11) is due to the lags noticed between the simulated and observed “peaks”. FFigiguurree1 100. .S Sccaatttteerrp pllootts shhoowwininggc coommppaarrisisoonno offs simimuulalatteedda annddo obbsseerrvveeddd disiscchhaarrggeess( a(a))f oforrt htheep peerrioiodd 22000088––22001100,, aanndd ((bb))f oforrt htheefl foloooddininggy yeeaarr2 2001100, ,a anndd( c(c))t htheec ucummuulalatitviveeo offs ismimuulalatetedda annddo obbsseerrvveedd ddisiscchhaarrggeea attB Béététérorouu. . Figure 11. Scatter plot showing comparison of simulated and observed discharges (a) for the period 2008–2010, and (b) for the flooding year 2010, and (c) the cumulative of simulated and observed discharge at Savè. Water 2022, 14, x FOR PEER REVIEW 16 of 23 Figure 10. Scatter plot showing comparison of simulated and observed discharges (a) for the period Water 2022, 14, 1192 2008–2010, and (b) for the flooding year 2010, and (c) the cumulative of simul1a6teofd2 3and observed discharge at Bétérou. FigureF i1g1u.r Sec1a1.ttSecra tptelroptl osthsohwowiningg ccoommpparaisroisnoonf soimf usilamteudlaantdedob asenrdve odbdsisecrhvaergde sd(ias)cfhoratrhgeepse r(iao)d for the period 2008–22001080–, 2a01n0d, a(nbd) (bfo) rfo trhtehe ffllooooddiinnggy eyaera2r0 1200, 1an0d, a(cn)dth e(cc)u mthuel actiuvme ouf lsaimtivulea toedf asnimd oublsaetrevded and observed dischadrigsceh aartg Seaavt Sèa. vè. In general, WRF-H captures well the variability of the seasonal discharge for Savè and Bétérou, but overestimates the cumulative total discharge (Figures 10c and 11c). These figures also indicate that the discharge is recorded from July to November, with a high increase in August–September. The total volume of discharge observed at Bétérou during the period 2008–2010 is 86,045 m3 against 100,182 m3 for the WRF-H simulation, with a rise of 16% as a difference. At Savè, the recorded total volume of water is about 186,475 m3, while the corresponding simulated discharge is approximately 220,367 m3, overestimating the observed discharge up to 18%. 3.4. Evaluation of the Soil Water Content Figures 12 and S2 show the daily averaged soil moisture extracted for soil depth between 0.5–2 cm [60,61] from the Climate Change Initiative (CCI) of the European Space Agency (ESA: https://www.esa.int/ESA (12 May 2018).; thereafter referred to as θCCI), respectively, for Bétérou and Savè, and the corresponding volumetric soil water content in the first Noah soil layer between 0 and 10 cm (referred to as θWRF−H). The θWRF−H is relatively high at the beginning of the simulation (January 2008: about 0.3 m3·m−3), but experiences a decreasing trend from the first months of the simulation, and reached approximatively the same value in February as those of other years (February 2009 and February 2010). This suggests that there is an excess of θWRF−H both at Bétérou and Savè at the initial time of the simulation, and that a two-month spin-up period appears to be sufficient for soil moisture in the first Noah LSM soil layer [33]. The θCCI values are globally lower than WRF-H simulations either for Bétérou or Savè during the rainy season and simulate reasonably well at the onset end the end of the season. However, despite the important bias, WRF-H correlated well with the observed. Water 2022, 14, x FOR PEER REVIEW 17 of 23 In general, WRF-H captures well the variability of the seasonal discharge for Savè and Bétérou, but overestimates the cumulative total discharge (Figures 10c and 11c). These figures also indicate that the discharge is recorded from July to November, with a high increase in August–September. The total volume of discharge observed at Bétérou during the period 2008–2010 is 86,045 m3 against 100,182 m3 for the WRF-H simulation, with a rise of 16% as a difference. At Savè, the recorded total volume of water is about 186,475 m3, while the corresponding simulated discharge is approximately 220,367 m3, overesti- mating the observed discharge up to 18%. 3.4. Evaluation of the Soil Water Content Figures 12 and S2 show the daily averaged soil moisture extracted for soil depth be- tween 0.5–2 cm [60,61] from the Climate Change Initiative (CCI) of the European Space Agency (ESA: https://www.esa.int/ESA (12 May 2018).; thereafter referred to as 𝜃 ), re- spectively, for Bétérou and Savè, and the corresponding volumetric soil water content in the first Noah soil layer between 0 and 10 cm (referred to as 𝜃 ). The ∙𝜃 is rela-tively high at the beginning of the simulation (January 2008: about 0.3 m3 m−3), but expe- riences a decreasing trend from the first months of the simulation, and reached approxi- matively the same value in February as those of𝜃 other y ears (February 2009 and February 2010). This suggests that there is an excess of both at Bétérou and Savè at the initial time of the simulation, and that a two-month spin-up period appears to be sufficient for soil moisture in the first Noah LSM soil layer [33]. The 𝜃 values are globally lower Water 2022, 14, 1192 than WRF-H simulations either for Bétérou or Savè during the rainy season and simu1la7toef 23 reasonably well at the onset end the end of the season. However, despite the important bias, WRF-H correlated well with the observed. FigFuigruer1e2 1.2.E Evvaalluuaattiioonn ooff tthhee sosoili ml moiostiustrue rbeetbweetwene tehne tshime usilmateudla wteitdh wWiRthF-WH RofF t-hHe foirfstth Neofiarhs Lt SNMo ah LSsMoil sloaiylelra (yferorm(f r0o tmo 100 tcom1)0 ancmd )thaen ddaitlhye advearialygea fvreormag CeCfIr o(fmromCC 0I.5( ftroo 2m cm0.5) otvoer2 tchme )Béotvéerorut he catchment. Bétérou catchment. FuFruthrtehremrmoroer,eF, iFgiguurere1s3 1a3n adndF iSg3u prereSs3enptr tehsee ndtatihlye sdoaili lmyosiositlumreo iins ttuhree foinurt hNeofaohu rLSNMo ah LSsMoil sloaiylelrasy oefr sthoef tW𝜃heRFW-HR Fs-imHusliamtiounla, trieosnp,ercetsivpeelcyt,i vfoerl yt,hfeo rbathsiensb aBséitnésroBué taénrdo uSaavnèd. ASasv è.Asobosbesrevrevde dfofro rthteh (with soil layer between 0–10 cm) where only a two-month spin-up is enoughe tθo1 WsiRmHu−lHate( wthiteh ssooili lmlaoyisetrubree,t wtheee dnif0fe–r1e0nccme b)ewtwheeeren othnely twa otw yoea-mrs oins th splianr-guep foisr ethneo suogihl latoyesri m10u–l2a0t0e cthme. sFooirl imnsotaisntcuer,e f,otrh tehde isfofeilr emnociestbuertew veaelunet hofe tthwe oseyceoanrds is lar(g10e–f4o0r ctmhe) asnodil tlhaey tehrir1d0 –(4200–01c0m0 c.mF)o sroiinl slataynerc,e t,hfeo srtathbeilistyo iils mreoaicshteudr eafvtearlu theeo sfimthuelasteicoonn d Water 2022, 14, x FOR PEER REVIEW (10o–f 4m0ocrme t)haannd otnhee ytehairr.d T(h4e0r–e1fo0r0ec, mfor) tshoiisl slatuydeyr, ctahsee,s wtaeb inlieteydi smroeraec thimede afoftre srptihne-uspim (muloarteio n18 of 23 oftmhaonr ea tyheaanr)o tnoe aydeeaqru.aTtehleyr esfimoruel,aftoer tthhei ssositlu mdyoicsatusere, wbeetwneeeedn m10o–r2e00ti mcme. for spin-up (more than a year) to adequately simulate the soil moisture between 10–200 cm. Figure 13. Daily average (five-day filtered) time series of soil water content (SWC) of the four Noah FigLuSrMe 1s3o.i lDlaayileyrs aovfetrhaegWe (RfFiv-He-dsiamyu flialttieornedat) Btiémtéero sue.ries of soil water content (SWC) of the four Noah LSM soil layers of the WRF-H simulation at Bétérou. 3.5. Evaluation of Uncertainty of WRF-H In order to evaluate the forecasting uncertainties of WRF-H (analyzed in Section 3.3), a stochastic kinetic energy backscatter scheme (SKEBS: [21,34,36]) is used and activated in WRF-H; it is referred to as WRF-Hydro-SKEBS. The purpose here is that the SKEBS ap- proach adds random perturbations with prescribed spatial and temporal decorrelations. In particular, SKEBS produces perturbation into the lateral boundary conditions. The am- plitude of the stochastic perturbations is chosen as the default in the WRF model. There- fore, we generated an ensemble of ten members for the evaluation of the uncertainties. Both stochastic physics and initial condition perturbations into WRF-Hydro-SKEBS result in an ensemble spread for the three rainy seasons over both Savè and Bétérou. Figure 14 shows that WRF-Hydro-SKEBS has a relatively large impact on precipitation and dis- charge results in the study region. The WRF-Hydro-SKEBS ensembles depicted well the discharge shape, but showed either an overprediction or underestimation of the discharge intensity. The assessment of the forecast uncertainty showed that the high variation of discharge amounts is associated with a speedy variation of the forecasted precipitation ensemble as demonstrated [62]. The ensemble also results in a large range of simulated discharge performance, as can be seen in Table S2. As in the case of the current study where the forecasting is performed under the convective-scale (5 km horizontal resolu- tion), [63,64] concluded that SKEBS (at convective-scale) provides a balance between in- creased ensemble spread and forecast bias change. This demonstrates the sensitivity of WRF-H to lateral boundary perturbations and confirms the uncertainty of the model re- garding discharge and precipitation simulations, which is of the uttermost importance for flood forecasting. Water 2022, 14, 1192 18 of 23 3.5. Evaluation of Uncertainty of WRF-H In order to evaluate the forecasting uncertainties of WRF-H (analyzed in Section 3.3), a stochastic kinetic energy backscatter scheme (SKEBS: [21,34,36]) is used and activated in WRF-H; it is referred to as WRF-Hydro-SKEBS. The purpose here is that the SKEBS approach adds random perturbations with prescribed spatial and temporal decorrelations. In particular, SKEBS produces perturbation into the lateral boundary conditions. The amplitude of the stochastic perturbations is chosen as the default in the WRF model. Therefore, we generated an ensemble of ten members for the evaluation of the uncertainties. Both stochastic physics and initial condition perturbations into WRF-Hydro-SKEBS result in an ensemble spread for the three rainy seasons over both Savè and Bétérou. Figure 14 shows that WRF-Hydro-SKEBS has a relatively large impact on precipitation and discharge results in the study region. The WRF-Hydro-SKEBS ensembles depicted well the discharge shape, but showed either an overprediction or underestimation of the discharge intensity. The assessment of the forecast uncertainty showed that the high variation of discharge amounts is associated with a speedy variation of the forecasted precipitation ensemble as demonstrated [62]. The ensemble also results in a large range of simulated discharge performance, as can be seen in Table S2. As in the case of the current study where the forecasting is performed under the convective-scale (5 km horizontal resolution), [63,64] concluded that SKEBS (at convective-scale) provides a balance between increased ensemble Water 2022, 14, x FOR PEER REVIEW spread and forecast bias change. This demonstrates the sensitivity of WRF-H to latera1l 9 of 23 boundary perturbations and confirms the uncertainty of the model regarding discharge and precipitation simulations, which is of the uttermost importance for flood forecasting. Figure 14. Cont. Figure 14. Ensemble (WRF-Hydro-SKEBS) of simulated hydrographs and precipitations at Savè and Bétérou. 4. Summary and Conclusions The present study explored the abilities of the fully coupled WRF-Hydro modeling system to simulate discharge and precipitation in the Ouémé river in West Africa. The main objective of the study was to determine an accurate flood forecasting tool that ex- hibits the interaction between atmosphere and land surface models, since the feedback of these interactions impact the precipitation dynamic [62,65]. Therefore, the supplementary hydrological extension of WRF model has been calibrated in offline mode for one year, and tested for two years using hourly outputs from WRF simulations. Optimized param- eters from the calibration were used to perform the fully coupled WRF-Hydro model, which was used to investigate the forecasting performance skills of WRF over the Ouémé. We have found that the coupled modeling system WRF-H performed slightly better re- garding the precipitation episodes than the WRF-only model (as already illustrated by [66] over the interested period (2008–2010)). The evaluation of simulated precipitation showed its good performance skills, and provides confirmation about the uncertainty of WRF-H to simulate precipitation [43,46]. Water 2022, 14, x FOR PEER REVIEW 19 of 23 Water 2022, 14, 1192 19 of 23 FigFuigreu r1e41. 4E.nEsenmsebmleb l(eW(WRFR-FH-Hydydroro-S-SKKEEBBSS)) ooff ssiimuullaatteedd hhyyddrrooggrraapphhssa annddp prerceicpipitaittaiotinosnast aSta Svaèvè and BétaénrdouB.é térou. 4. Summary and Conclusions 4. Summary and Conclusions The present study explored the abilities of the fully coupled WRF-Hydro modeling sysTtehme ptoressimenutl astteudiys cehxaprgloeraendd thpere acibpiiltiatiteiosn oifn ththee fOululéym céouripvleerdi nWWReFs-tHAyfrdicrao. mThoedeling sysmteamin otob jesicmtivuelaotfet hdeisscthuadrygwe asntdo pderetecrimpiitnaetiaonn aicnc utrhaet eOfluoéomd éfo rievcaesrt ing Wtoeosltt hAaftreicxa- . The mahiinb iotsbtjheectiinvtee roafc ttiohne bsetutwdeye nwaatsm toos pdheetreermanidnela nand saucrcfuacreatme ofdloeolsd, sfionrceectahsetifneegd tboaockl tohfat ex- hibthitess tehien tienrtaecrtiaocntsioinm bpaectwt teheenp raetcmipoitsaptihoenred yannadm liacn[d62 s,6u5r]f.aTchee mrefoodree,lsth, esisnucpep tlheme feenetdarbyack of thehsyed irnotloegraiccatlioenxtse nimsiopnacotf tWheR Fprmeocidpeilthaatisobne ednycnaalimbriact e[6d2i,n65o]f.fl Tinheemreofdoerefo, rthoen esuypeaprl,eamndentary hytdersoteldogfoicratlw eoxyteenarssiouns inogf WhoRurFl ymooudtpeul ths afrso mbeWenR cFasliimbrualtaetdio nins. oOfpfltiinmei zmedodpear faomr eotenres year, from the calibration were used to perform the fully coupled WRF-Hydro model, which was anuds teedstteodin fvoers ttwigoat eyethaersf ourseicnagst hinoguprelyrf oorumtpauntcse fsrkoilmls Wof RWFR sFimovuelratthioenOsu. Oémpéti.mWiezehdav pearam- etefrosu nfrdotmha tththee ccaoluibprlaedtiomno dweelirneg usysestde mtoW pReFrf-oHrmpe rtfhoerm fuedllysl icgohutlpylbeedt teWr rReFg-aHrdyindgroth me odel, whpircehc iwpiatast uiosnedep tios oindevsesthtiagnatteh ethWe RfoFr-eocnalystminogd peler(afosramlreaandcye islkluilsltsr aotfe Wd bRyF[ 6o6v]eorv tehret Oheuémé. Wein htearveset efdoupnerdio tdh(a2t0 t0h8e– 2c0o1u0)p).led modeling system WRF-H performed slightly better re- gardingT htehe vparlueactipiointaotifosnim eupliastoedepsr etchiapnit atthioen WshRoFw-eodnliyts mgooodepl e(rafso ramlraenacdeys kiillllus,starnadted by [66p]r oovviedre tshceo innfitremreasttieodn apbeoruiotdth (e20u0n8c–er2t0a1in0t)y). of WRF-H to simulate precipitation [43,46]. WRTFh-eH eavlsaolushatoiwoned oaf gsoimoduplaetrefdor mpraenccieptiotastiimonu lsahteodwisecdh airtsg eg, owoidth paeKrfGoErmeqaunaclet osk0.i7ll6s, and profovridthees pcoernifoidrm20a0ti8o–n2 0a1b0o. uTth tehero ubnusctenretassinotfyW ofR WF-HRFh-aHs tboe esnimasusleastsee dpraetcBipéittéartoiuo,na [n43,46]. inner-catchment of the Ouémé river at Savè, where it provided a good agreement with respect to observed discharge, with a KGE equal to 0.66. Additionally, WRF-H was able to capture the flood event which occurred in 2010 over both Savè and Bétérou. Indeed, in the WRF-Hydro simulation in the fully coupled mode, the atmospheric and hydrological processes are simulated in a consistent way, which enhances the confidence in the results. The uncertainty of predictability skills of WRF-H with respect to discharge in the Ouémé river at Savè was performed with an ensemble of ten members using a random perturbation scheme. Results demonstrated the large sensitivity of simulated discharge to perturbations introduced into the atmosphere. In summary, WRF-H is considered as a suitable model for evaluating discharge prediction uncertainties in the Ouémé river, and we encourage the implementation of the model for further basins in West Africa. Further studies are needed to assess how such WRF-H predicted discharges uncertainties could be used to improve an EWS in West-African regions. Water 2022, 14, 1192 20 of 23 Fully coupled WRF-H allows simulation of streamflow from global meteorological data. As such, coupled WRF-H can easily be applied to any region of the world, which is an advantage of the method used. Nevertheless, the application of WRF-H to another region may require a different calibration to optimize the quality of the modeled streamflow results over the concerned area. Additionally, the global forcing dataset to drive the model must be available and provided. Moreover, the model has some shortcomings; the human impacts (factors) on streamflow, such as irrigation and dam managements, are not yet considered in WRF-Hydro, which we see as a future and necessary model improvement for flood risk management purposes. Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/w14081192/s1, Table S1: Recapitulation of yearly highest discharge values both for simulated WRF-H and station dataset, followed by KGE during the rainy season. Table S2: Comparison of KGE between the deterministic model WRF-H, and WRF-Hydro-SKEBS for Savè (resp. Bétérou). Figure S1: cumulative total precipitation derived from TRMM and CHIRPS, and simulated in WRF-H over Savè during the period 2008–2010. Figure S2: Evaluation of the soil moisture between the simulated with WRF-H of the first Noah LSM soil layer (from 0 to 10 cm) and the daily average from CCI (from 0.5 to 2 cm) over the Savè catchment. Figure S3: Daily average time series of soil water content (SWC) of the four Noah LSM soil layers of the WRF-H simulation at Savè. Author Contributions: Conceptualization, G.M.L.D.Q., methodology, G.M.L.D.Q. and J.A.; software, G.M.L.D.Q. and J.A.; validation: G.M.L.D.Q., J.A., N.A.B.K., H.K. and P.G.O., formal analysis, G.M.L.D.Q. and J.A.; investigation, G.M.L.D.Q.; resources, G.M.L.D.Q.; writing—original draft preparation, G.M.L.D.Q.; writing—review and editing, G.M.L.D.Q., J.A., N.A.B.K., Z.Z., H.K. and P.G.O.; visualization G.M.L.D.Q. and J.A.; supervision, J.A, N.A.B.K., H.K. and P.G.O.; project administration, J.A., N.A.B.K. and H.K.; funding acquisition, G.M.L.D.Q. and N.A.B.K. All authors have read and agreed to the published version of the manuscript. Funding: We are grateful to the Deutscher Akademischer Austauschdienst (DAAD), which supports G.M.L.D.Q within the framework of the climapAfrica programme (German Academic Exchange Service: 57556650), and via the CONCERT and FURIFLOOD projects. This work is also supported by a grant from the Government of Canada, provided through Global Affairs Canada, www.international. gc.ca (accessed on 1 January 2021), and the International Development Research Centre, www.idrc.ca (Number: 108246-001). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: ECMWF: https://www.ecmwf.int/en/forecasts/dataset/operational- archive (accessed on 14 April 2018); TRMM: http://daac.gsfc.nasa.gov/precipitation/TRMM_ README/TRMM_3B43_readme.shtml (accessed on 12 May 2018); CHIRPS: https://data.chc.ucsb. edu/products/CHIRPS-2.0/africa_daily/ (accessed on 12 May 2018). Acknowledgments: We would like to thank Ilse Hamann and the German Climate Computing Center (DKRZ) for providing the computing facilities. The Center for High-Performance Computing (CHPC, Cape town, South Africa) also provided computing facilities for G.M.L.D.Q. for some parts of the present study. We would also like to acknowledge the European Center for Medium-Range Weather Forecasts (ECMWF), for providing the operational analysis dataset and products. The TRMM data used in this study were acquired as part of the Tropical Rainfall Measuring Mission (TRMM), and the observational station data obtained from the National Meteorological Agency of Benin (METEO-BENIN). Conflicts of Interest: The authors declare no conflict of interest. References 1. IPCC. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team, Pachauri, R.K., Meyer, L.A., Eds.; IPCC: Geneva, Switzerland, 2014; p. 151. 2. Hounkpè, J.; Diekkrüger, B.; Badou, D.F.; Afouda, A.A. Non-Stationary Flood Frequency Analysis in the Ouémé River Basin, Benin Republic. Hydrology 2015, 2, 210–229. [CrossRef] Water 2022, 14, 1192 21 of 23 3. Ryu, Y.; Lim, Y.-J.; Ji, H.-S.; Park, H.-H.; Chang, E.-C.; Kim, B.-J. Applying a coupled hydrometeorological simulation system to flash flood forecasting over the Korean Peninsula. Asia Pacific J. Atmos. Sci. 2017, 53, 421–430. [CrossRef] 4. UNHCR. 2010. Le Monde Afrique. Available online: https://www.lemonde.fr/afrique/article/2010/10/22/680-000-personnes- touchees-par-les-inondations-au-benin_1429957_3212.html (accessed on 13 December 2018). 5. Nka, B.N.; Oudin, L.; Karambiri, H.; Paturel, J.E.; Ribstein, P. Trends in West African floods: A comparative analysis with rainfall and vegetation indices. Hydrol. Earth Syst. Sci. Discuss. 2015, 12, 5083–5121. [CrossRef] 6. Mohammed, K. Comparative Study of Performance for Real-Time Flash Flood Forecasting in the Upper Meghna Basin. Master’s Thesis, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh, 2017. 7. Arnold, J.G.; Srinivasan, R.; Muttiah, R.S.; Williams, J.R. Large area hydrologic modeling and assessment part I: Model develop- ment. JAWRA J. Am. Water Resour. Assoc. 1998, 34, 73–89. [CrossRef] 8. Arnold, J.G.; Fohrer, N. SWAT2000: Current capabilities and research opportunities in applied watershed modelling. Hydrol. Process. 2005, 19, 563–572. [CrossRef] 9. Anderson, M.G.; Burt, T.P. Modelling strategies. In Hydrological Forecasting; Anderson, M.G., Burt, T.P., Eds.; John Wiley & Sons: Suffolk, UK, 1985; pp. 1–13. 10. Wilby, R.L. Contemporary hydrology: Towards holistic environmental science. In Hydrological Modeling in Practice; Watts, G., Ed.; John Wiley & Sons: Chichester, UK, 1997. 11. Kodja, D.J.; Quenum, G.M.L.D.; Mahe, G.; Paturel, J.-E.; Boko, M. Indicators of Hydroclimatic Extreme Hazards to Flooding in the Ouémé Watershed in Bénin Investigation on Climate Extreme Events in Africa View Project 4th FRIEND/IAHS International Conference on the Hydrology of African Large River Basins View Project. Available online: https://www.researchgate.net/ publication/344340892 (accessed on 21 June 2021). 12. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep residual learning for image recognition. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 770–778. 13. Devlin, J.; Chang, M.W.; Lee, K.; Toutanova, K. BERT: Pre-Training of Deep Bidirectional Transformers for Language Understand- ing. In Proceedings of the 2019 Conference of the North American Chapter of the Association for Computational Linguistics: 1; 2019; pp. 4171–4186, EM-DAT. Available online: https://www.emdat.be/ (accessed on 10 October 2021). 14. Agudelo-Otálora, L.M.; Moscoso-Barrera, W.D.; Paipa-Galeano, L.A.; Mesa-Sciarrotta, C. Comparison of physical models and artificial intelligence for prediction of flood levels. Tecnol. Cienc. Del Agua 2018, 9, 209–235. [CrossRef] 15. Mosavi, A.; Ozturk, P.; Chau, K.-W. Flood Prediction Using Machine Learning Models: Literature Review. Water 2018, 10, 1536. [CrossRef] 16. Kratzert, F.; Klotz, D.; Shalev, G.; Klambauer, G.; Hochreiter, S.; Nearing, G. Towards learning universal, regional, and local hydrological behaviors via machine learning applied to large-sample datasets. Hydrol. Earth Syst. Sci. 2019, 23, 5089–5110. [CrossRef] 17. Givati, A.; Gochis, D.; Rummler, T.; Kunstmann, H. Comparing One-Way and Two-Way Coupled Hydrometeorological Forecasting Systems for Flood Forecasting in the Mediterranean Region. Hydrology 2016, 3, 19. [CrossRef] 18. Givati, A.; Lynn, B.; Liu, Y.; Rimmer, A. Using the WRF Model in an Operational Streamflow Forecast System for the Jordan River. J. Appl. Meteorol. Clim. 2011, 51, 285–299. [CrossRef] 19. Fiori, E.; Comellas, A.; Molini, L.; Rebora, N.; Siccardi, F.; Gochis, D.; Tanelli, S.; Parodi, A. Analysis and hindcast simulations of an extreme rainfall event in the Mediterranean area: The Genoa 2011 case. Atmos. Res. 2013, 138, 13–29. [CrossRef] 20. Kerandi, N.M.; Laux, P.; Arnault, J.; Kunstmann, H. Performance of the WRF model to simulate the seasonal and interannual variability of hydrometeorological variables in East Africa: A case study for the Tana River basin in Kenya. Arch. Meteorol. Geophys. Bioclimatol. Ser. B 2016, 130, 401–418. [CrossRef] 21. Igri, P.M.; Tanessong, R.S.; Vondou, D.A.; Panda, J.; Garba, A.; Mkankam, F.K.; Kamga, A. Assessing the performance of WRF model in predicting high-impact weather conditions over Central and Western Africa: An ensemble-based approach. Nat. Hazards 2018, 93, 1565–1587. [CrossRef] 22. Liu, J.; Bray, M.; Han, D. Sensitivity of the Weather Research and Forecasting (WRF) model to downscaling ratios and storm types in rainfall simulation. Hydrol. Process. 2011, 26, 3012–3031. [CrossRef] 23. Uribe, I.M.; Lugo-Morin, D.R. Performance of the WRF model with different physical parameterizations in the precipitation simulation of the state of Puebla. Atmósfera 2020, 33, 357–383. [CrossRef] 24. Bouilloud, L.; Chancibault, K.; Vincendon, B.; Ducrocq, V.; Habets, F.; Saulnier, G.-M.; Anquetin, S.; Martin, E.; Noilhan, J. Coupling the ISBA Land Surface Model and the TOPMODEL Hydrological Model for Mediterranean Flash-Flood Forecasting: Description, Calibration, and Validation. J. Hydrometeorol. 2010, 11, 315–333. [CrossRef] 25. Chen, F.; Dudhia, J. Coupling an advanced land surface–hydrology model with the Penn State–NCAR MM5 modeling system. Part I: Model implementation and sensitivity. Mon. Weather. Rev. 2001, 129, 569–585. [CrossRef] 26. Seuffert, G.; Gross, P.; Simmer, C.; Wood, E.F. The Influence of Hydrologic Modeling on the Predicted Local Weather: Two-Way Coupling of a Mesoscale Weather Prediction Model and a Land Surface Hydrologic Model. J. Hydrometeorol. 2002, 3, 505–523. [CrossRef] 27. Jasper, K.; Gurtz, J.; Lang, H. Advanced flood forecasting in Alpine watersheds by coupling meteorological observations and forecasts with a distributed hydrological model. J. Hydrol. 2002, 267, 40–52. [CrossRef] Water 2022, 14, 1192 22 of 23 28. Marty, R.; Zin, I.; Obled, C. Sensitivity of hydrological ensemble forecasts to different sources and temporal resolutions of probabilistic quantitative precipitation forecasts: Flash flood case studies in the Cévennes-Vivarais region (Southern France). Hydrol. Process. 2012, 27, 33–44. [CrossRef] 29. Moreno, H.A.; Vivoni, E.R.; Gochis, D.J. Limits to Flood Forecasting in the Colorado Front Range for Two Summer Convection Periods Using Radar Nowcasting and a Distributed Hydrologic Model. J. Hydrometeorol. 2013, 14, 1075–1097. [CrossRef] 30. Zabel, F.; Mauser, W. 2-way coupling the hydrological land surface model PROMET with the regional climate model MM5. Hydrol. Earth Syst. Sci. 2013, 17, 1705–1714. [CrossRef] 31. Wagner, S.; Fersch, B.; Yuan, F.; Yu, Z.; Kunstmann, H. Fully coupled atmospheric-hydrological modeling at regional and long-term scales: Development, application, and analysis of WRF-HMS. Water Resour. Res. 2016, 52, 3187–3211. [CrossRef] 32. Senatore, A.; Mendicino, G.; Gochis, D.J.; Yu, W.; Yates, D.N.; Kunstmann, H. Fully coupled atmosphere-hydrology simulations for the central Mediterranean: Impact of enhanced hydrological parameterization for short and long time scales. J. Adv. Model. Earth Syst. 2015, 7, 1693–1715. [CrossRef] 33. Errico, R.M.; Langland, R.; Baumhefner, D.P. The workshop in atmospheric predictability. Bull. Am. Meteorol. Soc. 2002, 83, 1341–1343. 34. Zhang, Z.; Arnault, J.; Laux, P.; Ma, N.; Wei, J.; Shang, S.; Kunstmann, H. Convection-permitting fully coupled WRF-Hydro ensemble simulations in high mountain environment: Impact of boundary layer- and lateral flow parameterizations on land– atmosphere interactions. Clim. Dyn. 2021, 1–22. [CrossRef] 35. Berner, J.; Shutts, G.J.; Leutbecher, M.; Palmer, T. A Spectral Stochastic Kinetic Energy Backscatter Scheme and Its Impact on Flow-Dependent Predictability in the ECMWF Ensemble Prediction System. J. Atmos. Sci. 2009, 66, 603–626. [CrossRef] 36. Le Barbé, L.; Alé, G.; Texier, H. Les Ressources en Eaux Superficielles de la République du Bénin; Monographies Hydrologiques; ORSTOM: Paris, France, 1993; 540p, ISBN 2-7099-1168-X. ISSN 0335-6906. 37. Huffman, G.J.; Bolvin, D.T.; Nelkin, E.J.; Wolff, D.B.; Adler, R.F.; Gu, G.; Hong, Y.; Bowman, K.P.; Stocker, E.F. The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-Global, Multiyear, Combined-Sensor Precipitation Estimates at Fine Scales. J. Hydrometeorol. 2007, 8, 38–55. [CrossRef] 38. Funk, C.; Peterson, P.; Landsfeld, M.; Pedreros, D.; Verdin, J.; Shukla, S.; Husak, G.; Rowland, J.; Harrison, L.; Hoell, A.; et al. The climate hazards infrared precipitation with stations—A new environmental record for monitoring extremes. Sci. Data 2015, 2, 150066. [CrossRef] 39. Skamarock, W.C.; Klemp, J.B.; Dudhia, J.; Gill, D.O.; Barker, D.M.; Wang, W.; Powers, J.G. A Description of the Advanced Research WRF Version 2; Tech. Note NCAR/TN-4681STR; National Center For Atmospheric Research: Boulder, CO, USA, 2005; 101p. 40. Skamarock, W.C.; Klemp, J.B. A time-split nonhydrostatic atmospheric model for weather research and forecasting applications. J. Comput. Phys. 2008, 227, 3465–3485. [CrossRef] 41. Gochis, D.; Mc Creight, J.; Yu, W.; Dugger, A.; Sampson, K.; Yates, D.; Wood, A.; Clark, M.; Rasmussen, R. Multi-Scale Water Cycle Predictions Using the Community WRF-Hydro Modeling System; NCAR: Boulder, CO, USA, 2015; Available online: http://www.ral.ucar.edu/projects/wrf_hydro/ (accessed on 1 March 2022). 42. Friedl, M.A.; McIver, D.K.; Hodges, J.C.F.; Zhang, X.Y.; Muchoney, D.; Strahler, A.H.; Woodcock, C.E.; Gopal, S.; Schneider, A.; Cooper, A.; et al. Global land cover mapping from MODIS: Algorithms and early results. Remote Sens. Environ. 2002, 83, 287–302. [CrossRef] 43. Hong, S.Y.; Dudhia, J.; Chen, S.H. A revised approach to ice microphysical processes for the bulk parameterization of clouds and precipitation. Mon. Weather. Rev. 2004, 132, 103–120. [CrossRef] 44. Mlawer, E.J.; Taubman, S.J.; Brown, P.D.; Iacono, M.J.; Clough, S.A. Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res. Atmos. 1997, 102, 16663–16682. [CrossRef] 45. Dudhia, J. Numerical study of convection observed during the winter monsoon experiment using a mesoscale two-dimensional model. J. Atmos. Sci. 1989, 46, 3077–3107. [CrossRef] 46. Pleim, J.E. A Combined Local and Nonlocal Closure Model for the Atmospheric Boundary Layer. Part I: Model Description and Testing. J. Appl. Meteorol. Climatol. 2007, 46, 1383–1395. [CrossRef] 47. Pleim, J.E. A Combined Local and Nonlocal Closure Model for the Atmospheric Boundary Layer. Part II: Application and Evaluation in a Mesoscale Meteorological Model. J. Appl. Meteorol. Clim. 2007, 46, 1396–1409. [CrossRef] 48. Kerandi, N.; Arnault, J.; Laux, P.; Wagner, S.; Kitheka, J.; Kunstmann, H. Joint atmospheric-terrestrial water balances for East Africa: A WRF-Hydro case study for the upper Tana River basin. Arch. Meteorol. Geophys. Bioclimatol. Ser. B 2017, 131, 1337–1355. [CrossRef] 49. Arnault, J.; Wagner, S.; Rummler, T.; Fersch, B.; Bliefernicht, J.; Andresen, S.; Kunstmann, H. Role of Runoff–Infiltration Partitioning and Resolved Overland Flow on Land–Atmosphere Feedbacks: A Case Study with the WRF-Hydro Coupled Modeling System for West Africa. J. Hydrometeorol. 2016, 17, 1489–1516. [CrossRef] 50. Yucel, I.; Onen, A.; Yilmaz, K.K.; Gochis, D. Calibration and evaluation of a flood forecasting system: Utility of numerical weather prediction model, data assimilation and satellite-based rainfall. J. Hydrol. 2015, 523, 49–66. [CrossRef] 51. Moriasi, D.N.; Arnold, J.G.; van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [CrossRef] 52. Berner, J.; Fossell, K.R.; Ha, S.Y.; Hacker, J.P.; Snyder, C. Increasing the skill of probabilistic forecasts: Understanding performance improvements from model-error representations. Mon. Weather. Rev. 2015, 143, 1295–1320. [CrossRef] Water 2022, 14, 1192 23 of 23 53. Shutts, G. A kinetic energy backscatter algorithm for use in ensemble prediction systems. Q. J. R. Meteorol. Soc. 2005, 131, 3079–3102. [CrossRef] 54. Leutbecher, M.; Lock, S.; Ollinaho, P.; Lang, S.; Balsamo, G.; Bechtold, P.; Bonavita, M.; Christensen, H.M.; Diamantakis, M.; Dutra, E.; et al. Stochastic representations of model uncertainties at ECMWF: State of the art and future vision. Q. J. R. Meteorol. Soc. 2017, 143, 2315–2339. [CrossRef] 55. Judt, F.; Chen, S.S. Predictability and Dynamics of Tropical Cyclone Rapid Intensification Deduced from High-Resolution Stochastic Ensembles. Mon. Weather Rev. 2016, 144, 4395–4420. [CrossRef] 56. Lehner, B.; Verdin, K.; Jarvis, A. New Global Hydrography Derived From Spaceborne Elevation Data. Eos Trans. Am. Geophys. Union 2008, 89, 93–94. [CrossRef] 57. Naabil, E.; Lamptey, B.; Arnault, J.; Olufayo, A.; Kunstmann, H. Water resources management using the WRF-Hydro modelling system: Case-study of the Tono dam in West Africa. J. Hydrol. Reg. Stud. 2017, 12, 196–209. [CrossRef] 58. Zhang, Z.; Arnault, J.; Wagner, S.; Laux, P.; Kunstmann, H. Impact of Lateral Terrestrial Water Flow on Land-Atmosphere Interactions in the Heihe River Basin in China: Fully Coupled Modeling and Precipitation Recycling Analysis. J. Geophys. Res. Atmos. 2019, 124, 8401–8423. [CrossRef] 59. Klein, C.; Heinzeller, D.; Bliefernicht, J.; Kunstmann, H. Variability of West African monsoon patterns generated by a WRF multi-physics ensemble. Clim. Dyn. 2015, 45, 2733–2755. [CrossRef] 60. Dorigo, W.A.; Gruber, A.; De Jeu, R.A.M.; Wagner, W.; Stacke, T.; Loew, A.; Albergel, C.; Brocca, L.; Chung, D.; Parinussa, R.M.; et al. Evaluation of the ESA CCI soil moisture product using ground-based observations. Remote Sens. Environ. 2015, 162, 380–395. [CrossRef] 61. Wagner, W.; Dorigo, W.; de Jeu, R.; Fernandez, D.; Benveniste, J.; Haas, E.; Ertl, M. Fusion of active and passive microwave observations to create an essential climate variable data record on soil moisture. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, 7, 315–321. [CrossRef] 62. Arnault, J.; Fersch, B.; Rummler, T.; Zhang, Z.; Quenum, G.M.; Wei, J.; Graf, M.; Laux, P.; Kunstmann, H. Lateral terrestrial water flow contribution to summer precipitation at continental scale–A comparison between Europe and West Africa with WRF-Hydro-tag ensembles. Hydrol. Processes 2021, 35, e14183. [CrossRef] 63. Romine, G.S.; Schwartz, C.S.; Berner, J.; Fossell, K.R.; Snyder, C.; Anderson, J.L.; Weisman, M.L. Representing Forecast Error in a Convection-Permitting Ensemble System. Mon. Weather Rev. 2014, 142, 4519–4541. [CrossRef] 64. Duda, J.D.; Wang, X.; Kong, F.; Xue, M.; Berner, J. Impact of a Stochastic Kinetic Energy Backscatter Scheme on Warm Season Convection-Allowing Ensemble Forecasts. Mon. Weather Rev. 2016, 144, 1887–1908. [CrossRef] 65. Baisya, H.; Pattnaik, S.; Rajesh, P.V. Land surface-precipitation feedback analysis for a landfalling monsoon depression in the Indian region. J. Adv. Model. Earth Syst. 2017, 9, 712–726. [CrossRef] 66. Wang, W.; Liu, J.; Li, C.; Liu, Y.; Yu, F.; Yu, E. An Evaluation Study of the Fully Coupled WRF/WRF-Hydro Modeling System for Simulation of Storm Events with Different Rainfall Evenness in Space and Time. Water 2020, 12, 1209. [CrossRef]