Evaluation of the GPM-IMERG Precipitation Product for Flood Modeling in a Semi-Arid Mountainous Basin in Morocco

: A new precipitation dataset is provided since 2014 by the Global Precipitation Measurement (GPM) satellite constellation measurements combined in the Integrated Multi-satellite Retrievals for GPM (IMERG) algorithm. This recent GPM-IMERG dataset provides potentially useful precipitation data for regions with a low density of rain gauges. The main objective of this study is to evaluate the accuracy of the near real-time product (IMERG-E) compared to observed rainfall and its suitability for hydrological modeling over a mountainous watershed in Morocco, the Ghdat located upstream the city of Marrakech. Several statistical indices have been computed and a hydrological model has been driven with IMERG-E rainfall to estimate its suitability to simulate ﬂoods during the period from 2011 to 2018. The following results were obtained: (1) Compared to the rain gauge data, satellite precipitation data overestimates rainfall amounts with a relative bias of + 35.61% (2) In terms of the precipitation detection capability, the IMERG-E performs better at reproducing the di ﬀ erent precipitation statistics at the catchment scale, rather than at the pixel scale (3) The ﬂood events can be simulated with the hydrological model using both the observed and the IMERG-E satellite precipitation data with a Nash–Sutcli ﬀ e e ﬃ ciency coe ﬃ cient of 0.58 and 0.71, respectively. The results of this study indicate that the GPM-IMERG-E precipitation estimates can be used for ﬂood modeling in semi-arid regions such as Morocco and provide a valuable alternative to ground-based precipitation measurements.


Introduction
The Mediterranean regions are often affected by deadly flood events that are considered as the most dangerous hydrometeorological risk [1]. Floods are essentially caused by intense precipitation [2,3] which can generate destructive floods that cause several deaths [4,5]. The Moroccan High Atlas is one of the most vulnerable areas to such extreme events [6]. There, in addition to the intense precipitation, steep slopes cause destructive floods [7], such as the catastrophic event of the Ourika valley in the summer of 1995 that caused several hundred deaths [8]. More recently, in autumn 2014, intense precipitation and major floods in southwestern Morocco caused 47 deaths [9,10] and in the summer of 2019 in the Anti-Atlas Mountain, a deadly flash flood caused seven deaths. Moreover, the random and Water 2020, 12, 2516 2 of 15 often unpredictable Wadi's flow in these areas requires careful study of their hydrological processes to control and better manage them.
To reduce the adverse effects of floods, several approaches are usually implemented. Some are structural, such as draining ditching [11] or building protective structures, such as dams [12]; other approaches rely on non-structural measures including hydrological modeling for flood forecasting, which has become a key tool for prevention but also to improve the alert systems design [13]. This modeling approach, benefiting from technological advances in computational tools, allows the analysis and simulation of the hydrological behaviors of streams [14][15][16]. Among these models, the rainfall-runoff models are often considered to simulate runoff [17]. They are based on a more or less detailed description of the basin's surface, topography, soil types, and land use. These variables can be explicitly included in the model or not, depending on the data availability.
The physically-based hydrological models require several physical measurements distributed over a catchment, which often requires costly monitoring strategies, and in the case of Morocco this type of data is rarely available [7,18]. Oppositely, conceptual or empirical models only require a long time series of observed data to obtain a robust model evaluation allowing for their use to issue forecasts, so they are usually preferred for operational and real time applications [19]. In Morocco, the long time series is often available only at the daily time step. However, the time of the concentration of floods can be only a few hours in most Mediterranean basins in the case of fast floods, a type of flood prominent in semi-arid regions such as Morocco. Therefore, the application of hydrological modeling is challenged with the lack of the meteorological data at a fine temporal resolution, as is the case in the majority of North African countries. To deal with this problem, satellite precipitation products could be extremely useful, particularly in ungauged catchments or in areas with a low density of hydrometeorological monitoring networks [20][21][22] or even over the land areas where the spatial coverage from surface observations is not uniform [23].
In recent years, rainfall products derived from satellites have been developed and released to the scientific community. Among these products, the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN) [24], the National Oceanic and Atmospheric Administration Climate Prediction Center morphing technique product (CMORPH) [25], the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis products (TMPA) [26] and the Global Precipitation Measurement (GPM) [27] have been released. Several inter-comparison studies have been performed to evaluate the validity of precipitation satellite estimates in many fields, such as for hydrological simulation and water resources management at different spatial and temporal scales [28][29][30][31][32][33][34][35][36].
As the successor of TRMM, the GPM mission has been launched in 2014 by the National Aeronautics and Space Administration (NASA) and the Japanese Aerospace Exploration Agency (JAXA). Compared to TRMM, GPM incorporates three important improvements: (1) greater coverage of climate zones because of the increase in the orbital inclination angle of the GPM satellite from 35 • in the TRMM satellite to 65 • , (2) liquid and solid precipitation detection based on "high-frequency" passive microwave channels (165.5 and 183.3 GHz) (3) greater sensitivity of dual-frequency precipitation radar (DPR) to light precipitation and snowfall [26,37,38]. The GPM products are divided into different categories, the early, late, and final runs corresponding to IMERG-E, IMERG-L and IMERG-F, respectively. These products have a high spatial (0.1 • × 0.1 • ) and temporal (30 min) resolution within the 60 • N-S latitude band, that could provide precipitation at the adequate spatio-temporal resolutions for flood monitoring.
In Morocco, some studies have evaluated the performance of TRMM products. For instance, ref. [39] evaluated four TMPA products (3B42: V6, V7temp, V7, RTV7) against 125 rain gauges depending on climatic zones and topography in Northern Morocco for the period 1998-2012. They concluded that the TMPA research products (3B42 V6 and V7) performed better than the real-time product (3B42 RT). Further, ref. [22] assessed five different satellite-based precipitation products (TRMM-3B42 v6 and v7, RFE 2.0, PERSIANN-CDR, CMORPH1.0 version 0.x) in north Morocco. They indicated that the TRMM-3B42 v7 product is the closest to rain gauges data and can be used for hydrological modeling. Ref. [40] evaluated The TRMM 3B42 Version 7 product by using rain gauges measurements from 19 stations in the Oum-ErBia basin from 1998 to 2011. They indicated that this satellite product provides poor estimations of rainfall at the daily time scale but improved accuracy at monthly and annual time scale.
No studies have been conducted so far to evaluate the GPM-IMERG product for flood modeling, neither in Morocco nor in Maghreb countries with similar climatic conditions. To fill this gap, the purpose of this study is to evaluate (i) the accuracy and performance of the GPM-IMERG product to reproduce observed rainfall and (ii) the capability to model flood events with a hydrological model based on Soil Conservation Service-Curve Number (SCS-CN). This type of model is chosen since it is widely used in the Mediterranean region for floods [3]. This model could serve as a tool to assist decision makers (water basin agencies, water departments) in flood prediction and management, as traditional warning systems do not provide forecasts in time for effective intervention [41].

Study Area and Hydro-Rainfall Data
The Ghdat watershed is located in the upstream of the large Tensift basin which drains the High Atlas of Marrakech ( Figure 1). Covering an area of 552 km 2 , the altitudes are ranging from 693 to 3578 m. The basin, whose average altitude is 1730 m, has a large zone located between 1000 and 2600 m representing 82% of the basin's area. The average slope over the catchment, measured from a 30 m digital terrain model, is 30.7%. modeling. Ref. [40] evaluated The TRMM 3B42 Version 7 product by using rain gauges measurements from 19 stations in the Oum-ErBia basin from 1998 to 2011. They indicated that this satellite product provides poor estimations of rainfall at the daily time scale but improved accuracy at monthly and annual time scale.
No studies have been conducted so far to evaluate the GPM-IMERG product for flood modeling, neither in Morocco nor in Maghreb countries with similar climatic conditions. To fill this gap, the purpose of this study is to evaluate (i) the accuracy and performance of the GPM-IMERG product to reproduce observed rainfall and (ii) the capability to model flood events with a hydrological model based on Soil Conservation Service-Curve Number (SCS-CN). This type of model is chosen since it is widely used in the Mediterranean region for floods [3]. This model could serve as a tool to assist decision makers (water basin agencies, water departments) in flood prediction and management, as traditional warning systems do not provide forecasts in time for effective intervention [41].

Study Area and Hydro-Rainfall Data
The Ghdat watershed is located in the upstream of the large Tensift basin which drains the High Atlas of Marrakech ( Figure 1). Covering an area of 552 km 2 , the altitudes are ranging from 693 to 3578 m. The basin, whose average altitude is 1730 m, has a large zone located between 1000 and 2600 m representing 82% of the basin's area. The average slope over the catchment, measured from a 30 m digital terrain model, is 30.7%.  Geologically, the basin consists of two large groups: on the one hand, a rigid Paleozoic and Pre-Cambrian basement upstream, dominated by schists and quartzites, and downstream, soft Permo-Triassic deposits including sandstones, clays, conglomerates and marl-limestones of the Cenozoic age [42][43][44]. These formations are generally semi-permeable and favor heavy runoff in the case of intense rain. The hydrometric station of Sidi Rahal located at an altitude of 690 m, controls the discharge of the Ghdat Wadi and also serves as a rain gauge. It receives an average of 359 mm of rainfall per year. Rainfall and discharge data cover a 50-year period from the 1967-1968 to 2016-2017 hydrological years. These data come from gauges installed by the Tensift Hydraulic Basin Agency in Marrakech. This agency installed two other rain gauges upstream of the basin: Adrar Gardouz and Tijdante at altitudes of 2285 m and 1122 m, respectively. This made the Ghdat basin instrumented by three rainfall gauges (Sidi Rahal, Adrar Gardouz and Tijdante) and a hydrometric station (Sidi Rahal) ( Table 1). The time step of the rainfall and runoff data is 10 min. The data for modeling covers a period from 2011 to 2018. A total of 17 flood events are selected from the hydro-meteorological data. The observed precipitation data at 10-min time step were converted to 30-min to be compared with the satellite data.

The GPM IMERG-E Satellite Precipitation Product
The GPM mission, as an advanced successor to the TRMM, provides new global rainfall and snowfall observation products. The GPM Core Observatory satellite was launched in February 2014 by NASA and JAXA, equipped with a Dual-frequency Precipitation Radar (DPR) and a GPM Microwave Imager (GMI). The GPM mission has different products, organized using three levels: Level-1 (IMERG-E), Level-2 (IMERG-L), and Level-3 (IMERG-F). The IMERG system is run in near-real time and provides IMERG-E and IMERG-L with 4-h and 12-h latencies, respectively. The IMERG-F is calibrated with rain gauge analyses by the Global Precipitation Climatology Center (GPCC) and then would be available after 3.5 months. These products are provided at 0.1 • × 0.1 • spatial scale between 60 • S-60 • N and at half-hourly time scale. In this study, IMERG-E half-hourly data from 2011 to 2018 were downloaded from the NASA website (https://gpm.nasa.gov/data-access/downloads/gpm). This product was chosen due to its short time latency (4-h), making it interesting for the early warning of flood events. We extracted four pixels of GPM IMERG-E that cover the Ghdat basin and we interpolated observed precipitation data with the inverse distance method to represent areal precipitation over the studied basin.

Soil Moisture Dataset
To estimate the initial soil moisture condition prior to flood events, the European Space Agency Climate Change Initiative (ESA-CCI) data is considered. ESA-CCI V04.2 is a soil moisture product (http://www.esa-soilmoisture-cci.org/) available over a 38-year period from 1978 to 2018. This database includes three products: active, passive and combined. It provides daily soil surface moisture with a spatial resolution of 0.25 degrees [31,45,46]. The product has been validated by [47], and recently validated in the High Atlas of Marrakech by [18,48] for flood modeling. The ESA-CCI product shows in the mountainous area some missing values [49] that are related to a filter that is applied to the ESA-CCI product in order to ensure the data quality [45]. The ESA-CCI soil moisture pixel covering the Ghdat basin has been extracted.

Methods
In this section we describe the methodology adopted for the statistical and hydrological evaluation of GPM IMERG-E at different time scale. First, we compare the GPM IMERG-E product extracted from January 2011 to December 2018, against the ground precipitation observations using statistical indices. Second, hydrological simulations of flood events using a hydrological model is performed with either rain gauge data or GPM IMERG-E precipitation data.

Statistical Evaluation of Satellite Precipitation Products
To evaluate the accuracy of the GPM IMERG-E product against rain gauge observations in both at basin and grid-cell scale that corresponds to the GPM pixel covering the rain gauge. Four statistical indices were adopted; The Pearson's correlation coefficient (CC) describes the agreement between the satellite precipitation and ground observations (Equation (1)). Mean error (ME) measures the average error magnitude of satellite precipitation (Equation (3)). Root-mean square error (RMSE) indicates the average error magnitude between satellite data and observed precipitation data (Equation (4)). Relative bias (BIAS) reflects the systematic bias between the SPPs and the ground observations (Equation (2)). These indices are expressed as follows: where n is the sample size; S i is the satellite precipitation estimate; G i is the gauge observation; S is the mean satellite precipitation estimate; G is the mean gauge observation. Moreover, to quantitatively evaluate the precipitation detection ability of the satellite precipitation against the rain gauge observations, we used two indices that are based on the contingency table of satellite precipitation estimates with observations. The probability of detection (POD) (Equation (5)) measures the fraction of precipitation events correctly detected by the satellite. The false alarm ratio (FAR) (Equation (6)) describes the fraction of false events detected by the satellite. (Table 2).

Hydrological Model
In this study, we used the SCS-CN model implemented in the HEC-HMS software [50]. This software includes hydrological analysis procedures such as unit hydrographs and hydrologic routing (http://www.hec.usace.army.mil/software/hec-hms/). The SCS-CN model is based on the cumulative excess of precipitation by time [50]: where Pe is the excess of precipitation at time t, P is the rain accumulation at time t and S is the maximum retention capacity, that is related to the curve number (CN). The excess precipitation is transformed into direct runoff through the concept of unit hydrograph. The modeling tool used has several transfer functions: Clark, Snyder and SCS unit hydrographs, Modclark transformation and the kinematic wave. For event modeling, Clark's unit hydrograph [50] has been widely used in the semi-arid mountainous basins which are known by a complex and variable topography [18,51]: where A t is the area of the watershed, contributing at time t, A is the total watershed area, t c is the time of concentration of the watershed.

Hydrological Model Calibration and Validation
The calibration process aims to find the optimal parameters to minimize the objective function. The calibration for each parameter is performed independently of the others, to avoid interdependence between the parameters of the production and transfer functions. The calibration is performed both automatically and manually, to assess the sensitivity of the model to small variations of the parameters. Three efficiency criterions were used to evaluate the agreement between observed and simulated hydrographs; The Nash Sutcliffe (NS) efficiency coefficient [52] (Equation (9)) and the relative error calculated based on the observed and estimated total volume runoff and peak flow (bias) (Equations (10) and (11)) where Q o,i and Q s,i are the observed and the simulated flow, respectively, and Q o is the average observed flow. V s and V o are the simulated and the observed runoff depth, respectively. The CN parameter can be identified as the initial soil moisture condition in the model. Several studies have shown that this parameter can be related to different soil moisture indicators, obtained from field measurements [53][54][55], or soil moisture models or from satellite data [56] to reproduce the CN of flood events. In the present study, the approach previously used by [18,48] in a similar basin of South Morocco is implemented: the ESA-CCI soil moisture is used to estimate the initial soil moisture condition. All the model parameters including the CN are first calibrated. Then a linear regression is established between the CN and the ESA-CCI soil moisture the day before of each flood event.
For validation, a leave-one-out procedure (or Jack Knife) is considered, with each flood event in turn removed from the sample, the model is recalibrated on the remaining events. The mean values of the model parameters obtained after recalibration are used for the validation of the event that was removed. For the CN, the value is estimated from the linear regression against ESA-CCI soil moisture.

GPM IMERG-E Product Evaluation
The GPM IMERG-E product is evaluated against the rain gauge data using six statistical indices considered (CC, ME, RMSE, BIAS, POD and FAR) at both basin scale and the individual grid-cells corresponding to rain gauges during the available period from January 2011 to December 2018. The statistical evaluation indices of the IMERG-E product are presented at the basin scale in Table 3 and at the grid-cell scale in Table 4. At the basin scale, satellite precipitation data overestimates the amount of rainfall with a relative Bias of +35.61%. Several studies found that the IMERG precipitation products overestimate the amount of rainfall; in Italy [57], in Peru and Ecuador [58], in Iran [59] and in Spain [60].
About the overall agreement between the gauge observations and satellite precipitation data, the Pearson's correlation coefficient is equal to 0.10 and 0.14 at the hourly and daily time scale, respectively. Therefore, the accuracy of IMERG-E data improved over the daily time scale.
In terms of precipitation detection capability, the near real time IMERG-E precipitation data presents satisfactory performance and can be used to detect rain events over the Ghdat basin. The same results have been found in Spain with a low bias in the IMERG-E product that allows near-real time applications [61]. The daily estimates had a higher POD equal to 0.68 and a smaller FAR equal to 0.63 compared to the hourly (POD = 0.36 and FAR = 0.86) time scale.
For the grid-cell evaluation, the comparison of IMERG-E data with the data at the Sidi Rahal rain-gauge presents higher performance compared with others pixels; in which the RMSE is lower at the hourly time scale. Concerning the agreement between observed and satellite precipitation data, Sidi Rahal also shows the higher correlation coefficients compared to the other gauges. The GPM IMERG-E presents slightly higher performance in the downstream part of the studied basin. These results may be reflected the impact of the complex topography in the upstream portions of the basin. The good agreement between observed and satellite precipitation at the Adrar station is probably due to the very good quality of the data at this station; this gauge is located at a relatively high altitude but in an open area where there is no relief effect affecting rain measurements. This gauge is also considered by the Tensift Hydraulic Basin Agency as reference station, that is also used to relay data transmission from the other gauges with a poorer network coverage.
Several studies presented similar or comparable results. In terms of elevation effects on the IMERG product performance, in the Ebro River basin in Spain the IMERG performance strongly depends on elevation [60]. Additionally, in the arid western China with high altitudes, the performance of IMERG is less satisfactory than in south and east China, indicating a significant impact of topography [62]. Moreover, the IMERG presented better performance at low latitudes than at middle latitudes over mainland China [37]. Overall, the results of this statistical evaluation indicated that the IMERG-E achieved a better performance at the basin scale at hourly and daily time scales rather than at the pixel scale.

Hydrological Model Calibration
For rain gauge-based precipitation data, The NS values range from 0.40 to 0.97 with an average of 0.77. The calibrated model parameters are shown in Table 5. The mean relative bias on peak discharge and volume for all events are −7.07% and −0.7%, respectively (Table 6). For the events simulated with the GPM IMERG-E product, the average of NS is equal to 0.74 with a relative bias on peak discharge and volume of −20.57% and 3.99%, respectively. Thus, it can be concluded that the 17 events are well simulated (Figure 2). The model was able to reproduce the peak discharge and the shape of different hydrographs for the events occurred in the studied basin between 2011 and 2018 for both rain gauge-based precipitation data and satellite precipitation data (Figure 2 and Figure 4). These results are similar to the study of [63], who found that IMERG-E gave a good performance in term of model calibration using the HEC-HMS model. The overestimation of flood volume with GPM IMERG-E is consistent with the positive bias in IMERG-E found in the previous section. Table 5 shows the calibrated model parameters with observed or IMERG-E precipitation data. The CN values are very similar between the simulations with observed or IMERG-E precipitation. For some events, the Tc and Sc parameters show a greater variation between the two different sources of precipitation data, mostly due to compensation effects during calibration due to different rainfall inputs. For rain gauge-based precipitation data, The NS values range from 0.40 to 0.97 with an average of 0.77. The calibrated model parameters are shown in Table 5. The mean relative bias on peak discharge and volume for all events are −7.07% and −0.7%, respectively (Table 6). For the events simulated with the GPM IMERG-E product, the average of NS is equal to 0.74 with a relative bias on peak discharge and volume of −20.57% and 3.99%, respectively. Thus, it can be concluded that the 17 events are well simulated (Figure 2). The model was able to reproduce the peak discharge and the shape of different hydrographs for the events occurred in the studied basin between 2011 and 2018 for both rain gauge-based precipitation data and satellite precipitation data (Figures 2 and 4). These results are similar to the study of [63], who found that IMERG-E gave a good performance in term of model calibration using the HEC-HMS model. The overestimation of flood volume with GPM IMERG-E is consistent with the positive bias in IMERG-E found in the previous section. Table 5 shows the calibrated model parameters with observed or IMERG-E precipitation data. The CN values are very similar between the simulations with observed or IMERG-E precipitation. For some events, the Tc and Sc parameters show a greater variation between the two different sources of precipitation data, mostly due to compensation effects during calibration due to different rainfall inputs.
Despite the fact that IMERG-E better reproduces observed rainfall at the daily time step rather than at the hourly time step, the hydrological model driven by IMERG-E is able to reproduce floods with a good efficiency. This could be partly explained by the fact that the hydrological model parameters are being calibrated, so the adjustment of its parameters makes it possible to compensate for the bias on the satellite precipitation. This can be seen in the Sc and Tc parameter values ( Table 5) that are quite different between observed precipitation and IMERG. We also note slightly higher average CN values with the observed than with the IMERG precipitation, due to the overestimation of IMERG compared to the observed rainfall.  Despite the fact that IMERG-E better reproduces observed rainfall at the daily time step rather than at the hourly time step, the hydrological model driven by IMERG-E is able to reproduce floods with a good efficiency. This could be partly explained by the fact that the hydrological model parameters are being calibrated, so the adjustment of its parameters makes it possible to compensate for the bias on the satellite precipitation. This can be seen in the Sc and Tc parameter values ( Table 5) that are quite different between observed precipitation and IMERG. We also note slightly higher average CN values with the observed than with the IMERG precipitation, due to the overestimation of IMERG compared to the observed rainfall.

Hydrological Model Validation with Observed and GPM IMERG-E Precipitation
The calibrated CN are correlated with the ESA-CCI soil moisture data in order to estimate the initial soil moisture condition of the model. This correlation is good with a coefficient equal to r = 0.96 when the SCS-CN model is calibrated with rain gauge data (Figure 3).

Hydrological Model Validation with Observed and GPM IMERG-E Precipitation
The calibrated CN are correlated with the ESA-CCI soil moisture data in order to estimate the initial soil moisture condition of the model. This correlation is good with a coefficient equal to r = 0.96 when the SCS-CN model is calibrated with rain gauge data ( Figure 3). For the data based on the IMERG-E satellite precipitation, the correlation between the CN and the ESA-CCI data are also strong and equal to 0.95 ( Figure 3). This correlation confirms the good relationship between the CN parameter of the production function and ancillary source of soil moisture conditions, as demonstrated in similar studies [18,55,64]. The validation will therefore be based on ESA-CCI data.
The equation based on the correlation allowed us to calculate the validation CN of each event in turn by the leave-one-out method. This method involves correlating the calibrated CN values of the sixteen events with ESA-CCI soil moisture to determine the CN of the event that remains. Table 6 presents the results of the efficiency criterions used to evaluate the agreement between observed and simulated hydrographs. For the validated events, the average Nash values is equal to For the data based on the IMERG-E satellite precipitation, the correlation between the CN and the ESA-CCI data are also strong and equal to 0.95 ( Figure 3). This correlation confirms the good relationship between the CN parameter of the production function and ancillary source of soil moisture conditions, as demonstrated in similar studies [18,55,64]. The validation will therefore be based on ESA-CCI data.
The equation based on the correlation allowed us to calculate the validation CN of each event in turn by the leave-one-out method. This method involves correlating the calibrated CN values of the sixteen events with ESA-CCI soil moisture to determine the CN of the event that remains. Table 6 presents the results of the efficiency criterions used to evaluate the agreement between observed and simulated hydrographs. For the validated events, the average Nash values is equal to 0.58. The mean relative bias on peak discharge and volume are 16.8% and 20.4%, respectively. Concerning the events modeled with satellite rainfall, the model presents slightly better performances with a NS raised to 0.71. The relative bias on peak discharge and volume are −24.2% and −0.2%, respectively.
We were able to validate 13 out of 17 events simulated with rain gauge data ( Figure 2) and 15 events for satellite precipitation data ( Figure 4). It should be noted that this validation experiment combines the uncertainties stemming from the CN/ESA-CCI soil moisture relationship but also the uncertainties on precipitation estimates of the IMERG-E product. For non-validated events from observed data (i.e., with very low negative Nash values reported as N.A. in Table 6), a detailed analysis of the results reveals that for these events, the difference in CN values between calibration and validation is three time larger than the difference in CN values for the other events. Therefore, this indicates that for some events a slight change in the CN values can have a strong negative impact on the results. This finding denotes the importance of the soil moisture data to adequately estimate the antecedent soil moisture conditions, for the validation of flood events in a semi-arid environment [18,51,65].
For the floods simulated with satellite precipitation data, the events of the 1 May 2011 and 2 May 2011 present higher values of CNs in calibration compared to other events. For these two events, the cumulative rainfall recorded which is very low (2.36 mm and 4.44 mm, respectively) and the high CN values obtained reflect the overestimation of rainfall by IMERG-E. It would be required then for these two events to increase the value of the CN in order to adequately reproduce the observed hydrograph. It should be noted also that events modelled with ground observation or satellite precipitation do not have the same performances in terms of Nash, relative bias on flow and volume. We also compared the spatial variability of rainfall fields for each event, but there is no apparent link between the hydrological model performance to reproduce floods and the spatial variability of rainfall during the event. and validation is three time larger than the difference in CN values for the other events. Therefore, this indicates that for some events a slight change in the CN values can have a strong negative impact on the results. This finding denotes the importance of the soil moisture data to adequately estimate the antecedent soil moisture conditions, for the validation of flood events in a semi-arid environment [18,51,65].

Conclusions
This study conducted a first assessment of the near real-time satellite precipitation product GPM IMERG-E in Morocco for flood modelling. In this work, an event-based hydrological model was developed for a semi-arid mountainous watershed based on two sources of precipitation (observed and estimated from IMERG-E). Following the determination of model parameters, calibration and validation tests, the model was able to reproduce flood hydrographs at a short time step, which is adapted to this type of semi-arid small basin. The main findings of this study are summarized as follows: (1) Compared to the rain gauge data, satellite precipitation data overestimates the amount of rainfall with a relative Bias of 35.61%. (2) At the hourly and daily time scales, the GPM IMERG-E presents a fair accuracy to reproduce precipitation amounts and detect precipitation occurrence, with better performances at the daily time step. (3) The comparison reveals that basin-averaged GPM IMERG-E precipitation is in better agreement with interpolated precipitation over the catchment from observed data than a direct pixel-to-station comparison. (4) The flood events can be simulated using either satellite or rain-gauge precipitation with NS criterions of 0.58 and 0.71, respectively.
To conclude, the studied satellite precipitation product presents satisfactory performance in terms of statistical evaluation and hydrological modeling. For hydrological purposes, this source of precipitation could be used in ungauged basins and in data scarce regions like Morocco for flood modelling and contribute to the development of early warning systems. Concerning the application of the IMERG product in hydrological modeling, ref. [66] previously noted that the product was suitable for heavy precipitation and [61] showed the potential added value of the IMERG-E product in near-real time applications.
There is a need for further evaluation studies of this satellite product for hydrological modeling in other basins in north Africa and other semi-arid regions with similar climate characteristics. The type of event-based models developed in the present study with satellite precipitation and soil moisture data would represent an interesting alternative to continuous models since the data requirements for this type of approach are limited and could rely on remote sensing products, with GPM IMERG-E data for precipitation and ESA-CCI soil moisture. This model could be applied to simulate floods in un-gauged basins or basins where only daily data for precipitation and river discharge are available. This modelling scheme could be a valuable tool to help decision-makers (watershed agencies, water services) and to understand the watershed behavior during heavy rainfall episodes and better manage flood risks. To this end, the development of early warning systems for flood alerts would provide viable opportunities to reduce the damage to properties and fatalities linked to floods.