1. Introduction
Floods are one of the most devastating natural hazards in the world [
1] and their forecast is essential in flood risk reduction and disaster response decisions [
2]. Since 1995, 2.3 billion people have been affected by floods, and the death tolls caused by floods have risen in many parts of the world [
3]. In spite of the considerable efforts in flood disaster management by national and international organizations, floods still have a negative impact on terrestrial environments. In this context, the most efficient tools for flood disaster reduction are needed to provide timely emergency responses in large-scale areas. Flood detection and mapping are two important products from satellite-based remote sensing using optical and synthetic aperture radar (SAR) images. The selection of suitable sensors that are both cost-effective and efficient in the development of flood inundation maps has been a major challenge [
4]. In flood disaster reduction management, space technologies intervene in emergency responses by providing satellite-based flood detection systems such as the United Nations Platform for Space-based Information for Disaster Management and Emergency Response (UN-SPIDER) (Office for Outer Space Affairs, United Nations, Vienna, Austria) and the International Charter: Space and Major Disasters [
5]. Although analysts use visual interpretation or change detection methods [
6] for damage assessments by processing images acquired before and after flood disasters, these methods need to be improved, and space technologies can be incorporated with great potential. Flood forecasting is a complicated process [
7], so floods are very difficult to forecast quantitatively in some respects, such as intensity, the extent of flood discharge and water depth. Many countries in the world suffer from rain gauge deficiency. The use of ground stations to directly measure precipitation might increase the uncertainties due to limited spatial distribution of ground stations so that the accuracy of flood discharge estimation and mapping would be affected. To solve the problem, space technologies have been developed to measure precipitation using satellites. Due to their high temporal and spatial resolutions, satellite precipitation products (SPPs) have provided a new opportunity for flood discharge estimation.
For discharge estimation, two SPPs of the National Aeronautics and Space Administration (NASA, Washington, DC, USA) with different spatial and temporal resolutions can be used: (1) the integrated multi-satellite retrievals for global precipitation measurement (IMERG), specifically the GPM-3IMRGHH.05 product and (2) tropical rainfall measurement mission multi-satellite precipitation analysis (TMPA), the NRT version (3B42-RT).
The key part of any flood discharge estimation and forecasting system is the hydrological model [
8]. The physically-based rainfall–runoff (R–R) model that allows the transformation from precipitation to runoff is used for this purpose. Many categories of the R–R model exist [
9]. However, the challenge in this field is to improve the NRT discharge forecasting accuracy. To overcome this uncertainty, many researchers have highlighted the importance of only using in-situ data [
10].
The first objective of our study is to improve discharge forecasting accuracy by using SPPs to offer a large coverage area with more frequent data, which allows the initial condition estimation using the R–R model. For this purpose, we selected the R–R model “Research Institute for Hydrological Protection Model”, named “Continuous Lumped Hydrological Model” (MILc), which gave satisfactory results using the data from precipitation ground stations [
11]. In order to integrate SPPs into the MILc model, we developed an independent function for areal precipitation extraction based on the Theissen polygon equation. The second objective is, after discharge estimation, to evaluate the capacity to use the forecasted flood discharge from SSPs instead of ground measured discharge data for flood extent mapping. Discharge data was combined with the Digital Elevation Model (DEM) using a hydrodynamic model named the Hydrological Engineering Center River Analysis System (HEC-RAS) (HEC, Davis, CA, USA) [
12].
For the accuracy assessment, we selected SAR imagery acquired during the flood events as a benchmark, due to SAR image acquisition being independent from weather conditions and daytime. Moreover, the low values of water body backscattering in the microwave spectrum allows the detection and extraction of water bodies from other objects by using a threshold classifier.
The results of this study provide a new approach for early emergency responses that might be helpful to minimize the damage caused by floods.
3. Methodology
The early warning in flood disaster management has grown with the availability of real time data. The most crucial part of a flood forecasting system is the hydrological model. Many categories of flood forecasting methods exist [
22] with same expectations, namely the discharge hydrograph estimation, the volume of water in the stream, however, the difference concerns the modelling approach, the number of parameters and the type of used data. In our work, the R–R model based on SPPs in NRT was selected. Our approach was divided into two parts: hydrological forecasting in NRT and flood extent mapping in NRT.
The new global flood forecasting framework is described in
Figure 5:
3.1. Hydrological Forecasting in NTR
Hydrological forecasting aims to predict a flood discharge hydrograph based on SPPs. For this purpose, we selected an open source hydrological model which allows the use of SPPs. Kauffeldt et al. provide a comprehensive review of flood forecasting systems [
22], which permits an appropriate model selection. We chose the MILc R–R model developed by Brocca et al. [
9]. This model only applies for gauged watersheds, offering the possibility to validate expected results.
3.1.1. MILc Model
MILc consists of the soil–water balance model (SWB) and the rainfall–runoff model (R–R) that are used to simulate temporal soil moisture patterns and flood discharge, respectively. The two models have a linear relationship. It has been confirmed by [
17] that an accurate estimation of the antecedent wetness condition (AWC) on the hydrologic response of a watershed can provide a better accuracy for the flood hydrograph. The MILc model has been tested for flood simulation based on rain ground station data [
9,
10]. In this paper, we developed an independent function to adapt the MILc model to use SPPs.
The surface soil layer is assumed as a spatially lumped system, which can be expressed in the following water content balance equation:
where
t is time,
W(
t) is the amount of water in the soil layer,
f(
t) is the fraction of the precipitation infiltrating into the soil,
e(
t) is the evapotranspiration rate,
g(
t) is the drainage rate due to the interflow and/or the deep percolation and
Wmax is the maximum water capacity of the soil layer. The ratio
W(
t)/
Wmax represents the degree of saturation.
The R–R model is based on the soil conservation service curve number method [
23] aimed to estimate the direct runoff from precipitation excess for each sub-watershed. The respective sub-watershed drains water into the principal stream. Finally, the routing along the principal stream is estimated by the diffusive linear method.
Therefore, the precipitation excess,
εj(
t) for the element
j (j D 1, …, Nb) is given by the soil conservation service curve number formulation [
24]:
where
Rj(t) is the precipitation depth from the start of the rainstorm,
Sj is the soil potential maximum retention at the start of the rainstorm,
D is the diffusivity parameter and
λ1 is the parameter linked to the initial abstraction and assumed constant for all elements.
3.1.2. Model Adaptation to SPPs
In order to integrate SPPs with the MILc model, we developed an independent function for areal precipitation extraction, which allows us to use the global coverage SPPs in the MILc model for flood forecasting. Although many methods exist for this purpose, we selected the Thiessen polygon method [
25], which assumes that the precipitation value at a given area (A) is covered by a pixel (
i). So, the precipitation value observed at pixel
i is related only to its area. The weight of every pixel is determined by the corresponding area in the Thiessen polygon network. The Thiessen polygon equation is:
The following framework (
Figure 6) illustrates the adaption of the MILc model for SPPs.
3.1.3. Model Performance
For MILc model performance assessment, we selected the most widely used statistics on hydrological modeling, the Nash–Sutcliffe Efficiency (NSE) and the Adapted Nash–Sutcliffe Efficiency (ANSE), proposed by Nash and Sutcliffe [
26]. The ANSE, which is adapted to high flow conditions, is considered be more significant than NSE because the Ottawa River is characterized by high flow. NSE and ANSE ranges are between −∞ and 1. The United States Geological Survey (USGS, Reston, VA, USA) considers that the results are satisfactory when the NSE and ANSE values are higher than 0.5 and close to 1.
The following are the equations of NSE and ANSE:
and
where
Qobs is observed discharge and
Qsim is simulated discharge
3.2. Flood Mapping in NTR
The second part of the methodology aims to simulate a flood extent map using a two-dimensional (2D) flood model. For this, we propose a simple flood mapping approach in order to simulate flood extent based on open source data and software that would be compared with the observed benchmark flood extent. The HEC-RAS model requires geometric data that will be derived from DEM and discharge hydrograph data. Two discharge data were tested in different simulations and compared: (1) the forecasted discharge estimated by SPPs and (2) the observed discharge data by the ground station. Finally, to validate the simulated flood extent maps, we extracted the flood extent from the SAR imagery acquired during the flood event by RadarSat-2. The observed and simulated maps were compared using across topographical profiles for the accuracy assessment.
3.3. Flood Mapping Assessment
Based on the across topographical profiles, we estimated the distance in meters of boundaries between the simulated flood extent based on the forecasted discharge hydrograph using IMERG and the MILc model (S IMERG), and the simulated flood extent based on the observed discharge hydrograph using the ground station (S GDS) with the observed benchmark of the Ottawa River flood by RadarSat-2 (Obs).
Two methods are used to calculate the uncertainties. The first one is the absolute error (
A Err), which is the difference between the simulated map data and the observed map data as shown in the following equation:
The second one (
P Err) is used to calculate the error considering the pixel size (90 m) of the 2D mesh used in the HEC-RAS model for flood map simulations, based on the following equation:
3.4. Model Calibration
Calibration aims to adapt the model for the study area characteristics. The model parameters are physically based and estimated by the MILc model using the Ottawa River watershed data. Nine parameters have been estimated. These parameters and their ranges are shown in
Table 3.
5. Conclusions
It can be concluded from this paper that the estimation of discharge flood hydrographs based on IMERG products gives satisfactory accuracy. Consequently, flood forecasting based on IMERG products can play an important role in preparedness and mitigation of flood disaster management. The IMERG products might contribute to the use of space technologies more efficiently in flood risk reduction in emergency responses on global scale due to its large coverage.
The use of forecasted flood discharge hydrographs based on IMERG and the MILc model gives encouraging results on flood extent mapping. These results were evaluated and compared with simulated flood extent maps based on observed discharge hydrographs by using ground stations. The difference between the two simulations is insignificant. The two simulations were evaluated, and the results were satisfactory with the observed flood benchmark by RadarSat-2 imagery during the flood event. The use of SAR imagery as a benchmark was helpful in accuracy assessment, due to the difficulty of acquisition of benchmark data for natural phenomenon that can happen in a short time.