Next Article in Journal
Multi-Objective Optimization for Analysis of Changing Trade-Offs in the Nepalese Water–Energy–Food Nexus with Hydropower Development
Next Article in Special Issue
Large Differences between Glaciers 3D Surface Extents and 2D Planar Areas in Central Tianshan
Previous Article in Journal
Statistical Dependence of Pipe Breaks on Explanatory Variables
Previous Article in Special Issue
Improvements to Runoff Predictions from a Land Surface Model with a Lateral Flow Scheme Using Remote Sensing and In Situ Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physically, Fully-Distributed Hydrologic Simulations Driven by GPM Satellite Rainfall over an Urbanizing Arid Catchment in Saudi Arabia

by
Hatim O. Sharif
1,*,
Muhammad Al-Zahrani
2 and
Almoutaz El Hassan
1
1
Department of Civil and Environmental Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX 78249, USA
2
Department of Civil and Environmental Engineering, King Fahd University of Petroleum and Minerals, Water Research Group, Dhahran 31261, Saudi Arabia
*
Author to whom correspondence should be addressed.
Water 2017, 9(3), 163; https://doi.org/10.3390/w9030163
Submission received: 4 December 2016 / Accepted: 22 February 2017 / Published: 24 February 2017

Abstract

:
A physically-based, distributed-parameter hydrologic model was used to simulate a recent flood event in the city of Hafr Al Batin, Saudi Arabia to gain a better understanding of the runoff generation and spatial distribution of flooding. The city is located in a very arid catchment. Flooding of the city is influenced by the presence of three major tributaries that join the main channel in and around the heavily urbanized area. The Integrated Multi-satellite Retrievals for Global Precipitation Measurement Mission (IMERG) rainfall product was used due to lack of detailed ground observations. To overcome the heavy computational demand, the catchment was divided into three sub-catchments with a variable model grid resolution. The model was run on three sub-catchments separately, without losing hydrologic connectivity among the sub-catchments. Uncalibrated and calibrated satellite products were used producing different estimates of the predicted runoff. The runoff simulations demonstrated that 85% of the flooding was generated in the urbanized portion of the catchments for the simulated flood. Additional model simulations were performed to understand the roles of the unique channel network in the city flooding. The simulations provided insights into the best options for flood mitigation efforts. The variable model grid size approach allowed using physically-based, distributed models—such as the Gridded Surface Subsurface Hydrologic Analysis (GSSHA) model used in this study—on large basins that include urban centers that need to be modeled at very high resolutions.

1. Introduction

Storm events that result in catastrophic floods are rare but do occur in arid environments, especially in urban centers. With increasing population and urbanization, the public susceptibility and economic impact of flooding in these areas will be increasing [1]. Since these events do not occur frequently, there may not be enough pressure on decision makers to invest in the development of robust hydrometeorological observing systems or hydrologic/hydraulic flood control structures. The dramatic societal impacts of these events motivate researchers to perform studies aimed at developing science-based recommendations on best approaches to help decision makers address this issue [2]. Such information can lead to solutions that help save lives and resources and provide opportunities to harness floodwaters and turn them into a resource that can benefit the society. However, conducting hydrological studies of these events in ungauged areas is hampered by lack of adequate rainfall and physiographical data. Understanding the hydrometeorological conditions and processes that lead to destructive flood events is a first step in this process [3]. Physiographic data that enable the development of hydrologic models at reasonable resolutions over ungauged basins have become available globally in recent years. However, the accuracy of hydrologic model simulations is controlled by the accuracy of the model inputs, especially precipitation. A major challenge is the low quality and spatio-temporal resolutions of precipitation data for regions that do not have adequate ground observation networks.
High-resolution, better-quality satellite precipitation products that are increasingly becoming available can potentially lead to improvements in hydrologic modeling and forecasting, improvements as substantial as those witnessed in the United States in the 1990s when the NEXRAD radar network became operational. The potential of satellite precipitation products for various hydrometeorological applications has been reported in numerous studies [4,5,6,7]. Weather radars brought about the advantage of better spatial coverage, but satellites have even a better spatial coverage, though at lower temporal and spatial resolutions, and their field of view is not obstructed by topography. A major difference between radar and satellite precipitations products for hydrological applications is that when the NEXRAD network was deployed, the radar retrieval techniques were in place; however, it took many years to develop robust retrieval techniques after the launch of the Tropical Rainfall Measuring Mission (TRMM) satellite in late 1997. Also, systematic and random error of satellite precipitation products are typically higher than those associated with radar products [8,9]. Some of the satellite precipitation products that have been widely used in the past two decades include the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN) [10], the Climate Prediction Center (CPC) MORPHing (CMORPH) [11], the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) [12], and Global Satellite Mapping of Precipitation (GSMap) [13] products. These multi-sensor techniques that merge quality-controlled satellite products with higher resolution data from radars and rain gauges are helping improve the accuracy of satellite-based rainfall products over time [14]. In hydrological applications, however, satellite rainfall errors may either be amplified on dampened in simulated runoff at the catchment scale depending on the interaction between the spatio-temporal patterns of errors and catchment properties such as size, slope, and initial moisture conditions [15,16].
Simulation of flash floods, which are typically triggered by abrupt and intense bursts of rainfall, will benefit most from the fine resolution of the recent satellite rainfall products. For example, Anquetin et al. [17] reported that a higher resolution precipitation product was able to capture features of the precipitation system that caused the devastating 2002 flash flood in France. These features were missed by low resolution products. Several other studies demonstrated the need of high-resolution rainfall data for flash flood studies [18,19,20]. The number of hydrometeorological and climatological applications of satellite precipitation products will definitely increase with the emergence of the latest satellite product, the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG) product, with spatial and temporal resolutions of 0.1 × 0.1° and 30 min, respectively. IMERG is based on the Global Precipitation Mission (GPM), which was deployed in 2014 to consolidate and enhance precipitation measurements from a constellation of research and operational microwave sensors [21]. GPM is composed of one Core Observatory satellite, deployed by NASA and the Japan Aerospace Exploration Agency (JAXA), and carries a dual frequency radar and a multi-channel microwave imager, and about 10 partner satellites [22]. IMERG integrates the intermittent precipitation estimates from all GPM microwave sensors (high quality, low temporal resolution) with infra-red-based observations from geosynchronous satellites (lower quality, higher temporal resolution) and precipitation gauge data to produce a uniformly gridded, global, multi-sensor precipitation product. The IMERG product is designed to incorporate strengths and avoid major weaknesses of the previous multi-satellite algorithms supported by NASA: CMORPH-KF, TRMM-TMPA, and PERSIANN. This high-resolution precipitation data will significantly advance hydrological modeling and predictions worldwide, especially in ungauged and poorly gauged basins. The near-real-time IMERG Early and Late products are available within 6 h and 18 h after observations are made, respectively. After the gauge analysis is incorporated, the final satellite-gauge IMERG product becomes available, typically three or more months after the month in which the observations were made.
Since the hydrological response of a basin is very sensitive to the spatio-temporal variability in various physical attributes of soil, land use, and topography, hydrological models that consider the spatial variability are better suited for accurate flood simulation and predictions [23,24,25]. Distributed hydrologic models can also provide a detailed description of the flood hazard areas, especially in urban catchments [26]. Physically based, distributed models employ a gridded nature, which allows parameters to be constrained within certain ranges that have clear physical meanings. Recent research demonstrated that physically-based, distributed hydrologic models can potentially perform as well as—or outperform—calibrated conceptual, lumped models [24,25,26,27]. The physically-based Gridded Surface Subsurface Hydrologic Analysis (GSSHA) model is an example of a grid-based fully-distributed hydrologic models [28]. GSSHA is capable of simulating flow generated from Hortonian runoff, saturated source areas, exfiltration, and ground water discharge to streams [28]. Sharif et al. [29] successfully applied the GSSHA model to evaluate the effect of flood control structures on stream discharge in urbanized watersheds. Furl et al. [23] used the model to describe the flood hydrology of a small urbanized basin in Austin, TX. Chintalapudi et al. [30] used the GSSHA model to study the effect of land cover changes on peak discharge and runoff volumes with simulations driven by satellite rainfall products. Ogden et al. [31] compared the GSSHA distributed model to the HEC-HMS (a lumped model). Results showed that HEC-HMS failed to simulate some extreme events using standard parameters, whereas the GSSHA performed fairly well. Elhassan et al. [32] compared the simulated stream generated by the GSSHA and HEC-HMS models for different storm events. They concluded that the GSSHA simulated streamflow matched the observations much better than HEC-HMS. In addition to these studies, the model has been validated over a densely urbanized catchment in Texas and the results demonstrated the benefit of the use of 30-m model grid over urban areas [25]. Different satellite products we used as input to GSSHA in simulation of several floods over a 3000 km2 catchment in Texas [30] with satellite products of higher spatiotemporal resolutions producing the most reasonable runoff estimates. Another study over the Guadalupe River in Texas demonstrated that GSSHA was more successful in simulating events with multiple rainfall hiatuses than the HEC model [33]. An experiment using rainfall forecasts over a semi-arid urban catchment in Colorado demonstrated that GSSHA was able to produce reasonable forecasts of inundation and peak discharge for lag times of up to 70 min [34].
Recent extreme precipitation and flooding events in the Arab Peninsula led to several hydrometeorological studies. Furl et al. [35] analyzed rainfall in the southwestern region of Saudi Arabia and highlighted the lack of dense rain gauge networks. Almazroui [36] studied the TRMM rainfall data over Saudi Arabia throughout the period from 1998 to 2009. Although mixed results were obtained regarding accuracy of the TRMM estimates, the study recommended using the product to complement rainfall data from the extremely sparse rain gauge network in the country. A more recent study [37], evaluated the use of TRMM rainfall estimates for flood warning in urban areas of the country and concluded that TRMM satellite rainfall will provide some helpful information for preparation during extreme events but with low accuracy in terms of the spatiotemporal distribution of the rainfall storms. However, detailed fully-distributed hydrological analysis of runoff generation during these events is lacking, especially for urban areas that are most vulnerable to flash flooding events. A major reason for that was the lack of rainfall data at resolutions suitable for physically-based hydrologic modeling. The main objective of this study is to examine the flooding potential in the arid Wadi Al Batin catchment, Saudi Arabia, with focus on the rapidly urbanizing city of Hafr Al Batin, located near the outlet of the catchment. The city started to witness frequent flooding in the last two decades. A short-lived storm that hit the city in 2009 caused major ephemeral streams in the city to flow overbank, resulting in devastating flooding in the residential areas and significant damage to public and private properties. A flood event that occurred in October 2015 is employed as a case study in this paper. Detailed understanding of the runoff response in the urbanized part of the city of Hafr Al Batin will provide invaluable information that can help city officials to identify appropriate actions for reducing the probability of future flooding as well as to implement mitigation measures for severe storm events. GPM satellite precipitation data was used as input to a physically-based distributed hydrologic model. At model grid resolutions of 30 m to 700 m, computational requirements were reduced by dividing the basin into fully hydrologically interconnected sub-basins. The high-resolution simulations in the urbanized portion of the catchment helped identify the areas most susceptible to flooding. The relative contribution of the three streams that meet the main channel near the urban center was examined through model simulations.

2. Study Area

Hafr Al Batin is an old town in northeastern Saudi Arabia not very far from the borders with Kuwait and Iraq. The city has been witnessing a high rate of urbanization in the past few decades with a current population of over 300,000. Hafr Al Batin lies in the valley of the mostly-dry Wadi Al Batin (Figure 1) from which it takes its name and is the sole source of its groundwater supplies. Wadi Al Batin represents the now-disconnected upstream segment of an ancient large river, Wadi Al Rimmah, originates in western Saudi Arabia and empties in the Arabian Gulf [38,39]. The main channel is highly incised, a sign of a history of very frequent deluges [40], which is probably the reason for the name Hafr (incision in Arabic).
As the purpose of the hydrologic simulations of this study is to quantify the runoff generation and spatial distribution of flooding in the city of Hafr Al Batin (Figure 2), the catchment outlet was selected just northeast of the city limits such that all of the city is included in the simulation. Digital Elevation Models (DEMs) at 30-m resolution of the catchment were based on the Saudi General Directorate of Military Survey (GDMS) national DEM data. These DEMs were based on photogrammetry and ground control data [41]. ArcGIS 10.3 software was used for processing and resampling of the data to different resolutions (the grid cell size of the hydrologic model). The Watershed Modeling System (WMS) software [42] was used to delineate the stream network and sub-catchments based on the DEM at variable resolutions. The stream network based on aerial photography was superimposed on the network generated through WMS-based delineation to make adjustments when necessary.
The drainage area of the delineated catchment is 4273 km2, as seen in Figure 1, with the city of Hafr Al Batin representing only about 7.2% of the delineated area. The densely urbanized part represents about 4% of the total catchment. The catchment has a generally mild slope with elevation ranging between around 660 m above mean sea level at the western edge and 203 m at the outlet. The catchment is dominated by a combination of sand and gravel soil while the alluvial fans of the ephemeral streams (wadis) consist of weathered and fractured limestone and sandstone and permeable sediments. Soil and land use data over the catchment are shown in Figure 3. Land use/cover data were obtained from Spot 5 and Landsat remote sensing imagery. Soil and geologic data were obtained from the Saudi Geological Survey maps. Soil properties and parameters were derived from the digital soil map of the world [43]. These datasets were compiled and then processed in ArcGIS 10.3 to create GSSHA input files that represent the physical characteristics of the watershed.
The delineated catchment includes three other wadis that meet Wadi Al Batin (Figure 2) in or around the city. The North Fleaj originates northwest of the city and joins the main channel just north of the city. The Northwest Fleaj that joins Wadi Al Batin near the center of the city passes through a heavily urbanized area. The South Fleaj that originates southeast of the city and flows through sparsely urbanized areas except near its confluence with the Wadi Al Batin’s main channel in the northern part of the city.

3. Methods

3.1. Rainfall Data

The main hydrological simulations of this study were driven by satellite rainfall. The IMERG data were downloaded from the Precipitation Measurement Missions (PMM) website (http://pmm.nasa.gov/data-access/downloads/gpm). The three IMERG products, the early, late, and final were available for this event; all with spatial and temporal resolutions of 0.1 × 0.1° (approximately 11 × 11 km) and 30 min, respectively. An R-based script was used to download GPM products, convert the rainfall data from HDF5 into a gridded ASCІІ format, and prepare the rainfall input files for the GSSHA model. Figure 4 shows the 3–5 November total rainfall I the region as estimated the IMERG product.

3.2. Hydrologic Model

The physically-based, distributed-parameter Gridded Surface Subsurface Hydrologic Analysis (GSSHA) model is used to simulate recent flooding in Wadi Al Batin catchment. GSSHA is a process-based model for simulating all the hydrologic states and fluxes before, during and after storm events over each grid cell. GSSHA can accept spatially and temporally varying precipitation (e.g., from gauges, radar, satellite, or design storms). It can also ingest snowfall accumulation and simulate snow melting, abstractions due to interception, evapotranspiration, surface retention, and infiltration. The GSSHA model uses a simple two-parameter scheme to model interception of precipitation by plants. The user can select one of four infiltration methods: Green and Ampt (GA), Green and Ampt with Redistribution (GAR, [44]), multi-layered GA, and fully-implemented Richards’s equation [45]. The last two methods are best suited for modeling continuous storm event with significant hiatuses. GSSHA can simulate overland runoff routing, unsaturated zone soil moisture dynamics, saturated groundwater flow, surface sediment erosion transport and deposition, in-stream sediment transport, simplified lake storage and routing, wetland peat layer hydraulics, and overland contaminant transport and uptake [45].
The GAR method was used in simulating the flooding event of November 2015. The method, which computes inter-storm redistribution of soil water and performs multiple ponding simulations using the GA methodology, is based on the following equation:
f ( t ) = K ( S f ( θ s θ i ) F ( t ) + 1 )
where:
  • f ( t ) —potential Infiltration rate (cm/h)
  • F(t)—cumulative Infiltration (cm)
  • Sf—wetting front suction head (cm)
  • K—effective hydraulic conductivity (cm/h) = Ks/2.0
  • Ks—saturated hydraulic conductivity (cm/h)
  • θs—water content of the soil at natural saturation.
  • θi—initial soil water content.
The GAR methods uses the traditional rectangular wetting form assumption to execute variations of this equation to predict infiltration for multiple periods of ponding [44]. The configuration of GAR is illustrated in several figures and discussed in detail by Ogden and Saghafian [44]. After subtracting all abstraction, ponded flow over each grid cell is computed using the diffusive wave approximation of Saint-Venant’s equation and routed into two orthogonal directions. For grid cells adjacent to the watershed divide, only inwards flow is allowed [46]. The GSSHA model uses three numerical schemes to solve the diffusive wave equation: the Explicit, Alternative Direction Explicit (ADE), and ADE-Prediction Correction (PC) schemes. The physiographic conditions of the catchment dictate the most appropriate method. The Explicit is the fastest, simplest, and least robust method and the ADE-PC is the slowest and most robust method [46]. The ADE-PC, which is more demanding computationally, is generally recommended for watersheds of complex terrain when minimal smoothing of the DEM is warranted. The ADE scheme is used in this study. The ADE method uses the following formulas to calculate the flows.
First, inter cell flows are calculated in the x-direction by using Equation (2).
p i j N = 1 n ( d i j N ) 5 3 ( s f x N ) 1 2
Based on the flows in the x-direction, depths in each cell are calculated at the n + 1 time level by using Equation (3).
d i j N + 1 2 = d i j N + Δ t Δ x ( p i 1 , j N p i j N )
Equation (4) is used to calculate the interflows in the y-direction from each cell.
q i j N + 1 2 = 1 n ( d i j N + 1 2 ) 5 3 ( s f y N + 1 2 ) 1 2
Column depths are updated based on the interflows in the y-direction. Equation (5) is used to update the column depths.
d i j N + 1 = d i j N + 1 2 + Δ t Δ x ( q i , j 1 N + 1 2 q i j N + 1 2 )
where:
  • pij and qij are the overland flows from cell ij in the x and y directions, respectively
  • dij is the depth of water in cell ij at the Nth time level
  • sfx and sfy are the water surface slopes in the x and y directions, respectively.
  • n is Manning’s roughness coefficient
In this study, the Explicit solution scheme is used for solving the diffusive wave equation for 1-D channel routing. Overland flow that is entered into the stream is routed until it reaches the outlet. The volume of flow at each node is calculated using Equation (6).
V i N + 1 = V i N + Δ t ( q l a t N + 1 Δ x + q r e c N + 1 Δ x + Q i 1 / 2 N Q i + 1 / 2 N )
where:
  • qlat—amount of lateral flow (m2/s)
  • qrec—amount of flow exchanged between the groundwater and channel (m2/s)
Q i 1 / 2 N   a n d   Q i + 1 / 2 N are the inter cell flows in the longitudinal direction (x) computed from depths d, at the Nth time level. A simplified flowchart of GSSHA model is shown in Figure 5.

3.3. Model Setup

The GSSHA model step follows the flowchart shown in Figure 5. Al Batin catchment was divided into three sub-catchments: the upper sub-catchment of entirely barren desert with an area of 3117 km2, the middle sub-catchment with increased density of small dry channels that join the main channel (an area of approximately 847 km2), and the lower sub-catchment, which includes the city of Hafr Al Batin and its surroundings and covers an area of 309 km2. To provide the level of detail needed to simulate runoff over the three sub-catchments, three sizes of GSSHA model grid were used, 270 × 270 m2, 90 × 90 m2, and 30 × 30 m2 for the upper, middle, and lower sub-catchments, Figure 4. GSSHA model has to be run on each sub-catchment separately to allow for a variable grid size.
The hydrologic connectivity between the three adjacent sub-catchments is maintained through the downstream channel flow from an upstream sub-catchment because there is no sub-surface flow interchange among them. For every model time step, outflow from the upstream sub-catchment is added as inflow to the most upstream cell of the corresponding channel in the downstream sub-catchment. This inflow is merged with the channel flow in the same manner lateral inflow is treated at every time step—i.e., it is not just added to the downstream channel flow—as is typically done in popular semi-distributed models. This method of hydrologic connectivity (not simple routing) will ensure that dividing the catchment into sub-catchments will have minimum effect on the hydrograph and the water balance at every simulation time step. The model has to be run on the three sub-catchments in sequence or in parallel separated by one time step. A script was written to run the simulations in sequence and transfer of output from one sub-catchment to the next. The three sub-catchments are shown in Figure 6 together with the location of the outflow from one sub-catchment into the next. Since the sub-catchment delineation was not based entirely on topography, overland flow across the divide between the sub-catchments is not included. The sub-catchments were delineated such that overland flow between sub-catchments is minimal and can be neglected in the simulation.
The 30 × 30 m DEMs were used to delineate the channel network. Minor adjustments were performed after comparing the delineated network with aerial photographs of the catchment. All stream channels within the urbanized sub-catchment were modeled with irregular cross sections as well as the main channel in the other two sub-catchments. Uniform trapezoidal channels were used for the tributaries in the un-urbanized sub-catchments. Manual adjustment of stream channels using the WMS ‘smoothing’ tools were used to adjust the profiles of the stream channels to remove several regions of adverse (negative) channel slope resulting from errors in the DEM. Infiltration is simulated using the Green and Ampt with redistribution and the ADE method was selected for overland flow routing. The infiltration model parameters’ Manning roughness coefficients were taken from the GSSHA manual [46]. The initial soil moisture is estimated by running GSSHA model simulations over the few weeks prior to the event and extracting initial soil moisture for each model grid from the final map of spatially-distributed soil moisture values.

4. Results

The storm event started on 2 November 2015 at 9: 30 p.m. local time and ended at 11:30 a.m. on 3 November 2015. The spatial distribution of the total rainfall accumulations over the catchment estimated by the three GPM products, using Inverse Distance Weighted (IDW) interpolation, is shown in Figure 7. For the three products, the high amounts of rainfall fell on the urbanized portions of the catchment and immediately to the west of the city of Hafr Al Batin. Only two rain gauges were operational during the storm and they reported total accumulations (Figure 7). The Ministry of Environment, Water, and Agriculture rain gauge located inside the city recorded 63 mm while the Presidency of Meteorology and Environment rain gauge at the airport south of the city recorded 32 mm. The satellite grids collocated with the two gauges reported total accumulations of 47 mm and 23 mm for the Early IMERG product, respectively. The total storm precipitation averaged over the entire catchment for this product is 33.6–32 mm, 48 mm, and 47 mm over the upper, middle, and lower sub-catchments, respectively. The other two satellite products estimated slightly higher rainfall with the final IMERG product reporting the highest rainfall and a somewhat different spatial pattern of the storm. The total storm rainfall averaged over the entire catchment for the Late IMERG products are 37.2 (28.5, 44.7, and 43.1 mm for the three sub-catchments, respectively) and 46.2 mm for the final product (36.5, 52.9, and 51.9 mm for the three sub-catchments, respectively). The temporal distributions of rainfall the three IMERG products are shown in Figure 8.
The GSSHA model was run for the period between 12:00 p.m. local time on 30 September 2015 and 10:30 a.m. on 5 November 2015 for the three catchments. The run was started one month in advance to spin up the model in order to obtain a reasonable estimate of the initial moisture content for each model grid (see Chintalapudi et al. [30,47] for details of model spin-up). The three GPM precipitation products were used to force the GSSHA model but the results forced by the final product will be discussed since this product is the one adjusted with climatology data and ground observations. The storm lasted for a period of 14 h which allowed for infiltration of most of the rainfall except in the urbanized areas. The discharge from the upper sub-catchment into the middle sub-catchment was very small due the small amount of rainfall and high infiltration in the barren dominantly sandy soils. The runoff ratio for the event was 1% and most of the water that reached the main channel infiltrated in the channel. The middle sub-catchment generated modest amounts of runoff with a runoff ratio of 5%. Most of the flow from the middle sub-catchment into the lower one occurred through the main channel. The discharge at the watershed was primarily due to runoff generated within the urbanized sub-catchment with inflow form the upstream sub-catchments representing more than 12% of the total discharge. The runoff ratio in the urbanized area was 15% due to higher rainfall, the significant impervious fraction, and efficient drainage by the natural channel network.
The peak discharge of 127 cms at the outlet occurred at 8:00 a.m. on 3 November 2015. Peak discharges were 10, 50, and 90 min earlier for the North Fleaj, the South Fleaj, and the Northwest Fleaj tributaries, respectively. The South Fleaj contributed most to the total discharge at the outlet (about 46%) and has the highest peak among the tributaries. The contribution from the North Fleaj and North West Fleaj was very small (each contributed just about 5%). Lack of development around the South Fleaj indicates that this tributary witnesses frequent flooding, even more than the main Al Batin Wadi. Underestimation of rainfall by the Early and Late IMERG products was amplified in runoff simulation, producing much smaller outlet peak discharges of 44.2 and 66.7 cms, respectively.
The predicted outlet hydrographs when the model was forced by the three IMERG products are shown in Figure 9. As described in the introduction section, the final IMERG product is the best rainfall product. There are no streamflow observations for the event, however, GSSHA simulations highlight the fact that rainfall errors will result in higher runoff errors and give an idea about the magnitude of runoff error (in this case, difference between using calibrated and raw satellite rainfall) to be expected when the real-time product is used for flood forecasting.
GSSHA outputs also include the overland and channel flow depths at every time step. The map shown in Figure 10 illustrates the maximum inundation over the city area as a result of the storm—only depth values above 3 cm are shown. The water depth values for each grid were estimated by GSSHA forced by the final IMERG product. The map indicates that most of the flooding occurred in the highly urbanized areas and near the center of the city. The map also shows significant flooding over some major streets. There is significant flooding also over the undeveloped areas adjacent to the South Fleaj. The storm did not result in any flooding at the confluences of the North West Fleaj and South Fleaj with Wad Al Batin. It is clear that the inundation is caused by topography and neither the Wadi Al Batin nor the South Fleaj witnessed overbank flooding. These results agree with media reports and photographs released after the event that described significant street flooding without a mention of channel overflow. The map shows that there is hardly any flooding outside the urban area, reinforcing the fact that the city was built over the ancient flood plain of Wadi Al Batin, which might be the reason for the city name, and the role of urbanization in increasing the flood hazard. When GSSHA was forced by the Early and Late GPM products, the inundation maps illustrate smaller flooding depths and extent.
To understand the role of tributaries and inflow from the upper portions of the catchment in the city of Hafr Al Batin flooding, two more GSSHA simulations were performed: a simulation representing the frequent flood events in the city and another simulation representing extreme events. GSSHA was run with the 5-year, 24-h, and 100-year, 24-h storms, which are considered spatially uniform. The two storms were based on statistical analysis of the historical rainfall data from the Ministry of Environment, Water, and Agriculture rain gauge with a record starting from 1979. The 24-h rainfall accumulations for the 5-year and 100-year, 24-h storms were found to be 31 mm and 73 mm, respectively, using Log-Pearson Type III distribution (LPT III) [48]. Due to the size of the catchment, an area-reduction factor of 0.8. Was used (see [49,50]). The temporal distributions of the designed storms, which are based on Type II storms [51], are shown in Figure 11. GSSHA simulation results indicate that the upper sub-catchment contributes less than 1% of the total discharge flooding in the city for the 5-year storm. Contribution from the middle sub-catchment is also small, about 7%, for this event. The 5-year event does not generate much discharge at the outlet. The South Fleaj, Northwest Fleaj, and North Fleaj tributaries contribute about 52.2, 13.6, and 6.9% of that discharge, respectively. For the 100-year event, the contributions from the middle sub-catchment is significant. The peak discharge at the outlet of the upper sub-catchment is 5.5 cms and it contributes less than 3% of the total discharge at the outlet. The middle sub-catchment contributes 21% of the total discharge with a peak of 50 cms at the main channel. Again, South Fleaj contributed most to the total discharge of the lower sub-catchment (53%) and has the highest peak among the tributaries while the other two tributaries contributed smaller amounts −6% from the Northwest Fleaj and 5% from the North Fleaj.

5. Summary and Conclusions

In this study, the physically-based, distributed-parameter hydrologic model GSSHA, forced by Integrated Multi-satellite Retrievals for Global Precipitation Measurement Mission (IMERG) rainfall product, was used to study a recent flood event in the city of Hafr Al Batin, Saudi Arabia. Due to the large size of the catchment that encompasses the city, it was divided into three hydrologically connected sub-catchments: a lower sub-catchment that represents the urban area, an upper sub-catchment dominated by barren desert land, and a middle sub-catchment. A variable model grid size of 270 × 270 m2, 90 × 90 m2, and 30 × 30 m2 was adopted for the upper, middle, and lower sub-catchments.
The semi-real time IMERG products underestimated the event rainfall with much more pronounced underestimation of the event runoff. The IMERG product is relatively news and is currently undergoing continuous enhancement. The references cited in the introduction section of this paper provide early assessments and discussions of the validity of the product. These products are expected to improve as the IMERG algorithms get refined through more ground validations. Nonetheless, they provide valuable information that can help improve flood prediction in ungauged basins. The authors are not aware of published hydrologic applications of the IMERG products. There was no streamflow data to validate the hydrographs predicted by the hydrologic. However, the model has been validated over similar environments in many previous studies as described in the introduction section. Notwithstanding, the GSSHA model simulations forced by the final IMERG product enabled quantifying the relative contribution of the sub-catchments and the major tributaries in the urbanized area to flooding of the city. For this event, the distributed model simulations demonstrated that most of the flooding (approximately 85%) was generated in the urbanized portion of the catchments (6.8% of the total area of the simulated catchment). The contribution of the upper portion of the catchment (68% of the area) was insignificant due to it sandy soils and the limited amounts of rainfall it received. The middle sub-catchment contributed about 13% of the discharge at the outlet.
One of the tributaries, the Northwest Fleaj, that meets the main channel inside the city, does not play an important role in flooding of the city center due to the size of its drainage area and its physiographic feature. For the same reasons, the North Fleaj, which meets the main channel just north of the city, does not contribute significant discharge. However, the South Fleaj that meets the main channel just north of the city center contributes significantly to the total discharge in the main channel and therefore has to be considered in any future flood control projects and urban development in the eastern part of the city. As this tributary originates well outside the city, different types of flood control measures can be applied to reduce the discharge it contributes. Distributed model simulations demonstrate that flooding in the city is driven primarily by topography rather that overbank flow in the main channel. A well designed urban drainage network might be needed to prevent flooding of residential areas and streets. This information is important if the city officials want to implement flood mitigation measures. For example, the presence of North Fleaj will limit expansion of the urban area to the north because of increasing flood hazards unless its flow is diverted before entering the urban area. Large detention basins outside the city on the main channel and Northwest Fleaj and South Fleaj can help mitigate flooding caused by extreme event. Also, diversion of the Northwest Fleaj and South Fleaj flow outside the city can be helpful.
Additional simulations were performed to understand the roles of the unique channel network in flooding in the city of Hafr Al Batin. The barren desert upper portions of the catchment contribute to flooding in the city only when it receives significant amounts of rainfall. It contributes about 3% of the total discharge when the 100-year storm covers the entire catchment, which should have a probability of much less than 1%. No significant discharge results from the 5-year storm. The middle sub-catchment contributes 7% of the total runoff for the 5-year storm and 23% for the 100-year storm. The results demonstrate that the upper portions of the catchment do not pose a significant flood threat unless they receive exceptional amounts of rainfall. In that case, the early or even the late satellite rainfall products can be invaluable since the flood peak will arrive many hours after the precipitation peak for such a large catchment. The best approach to control flooding in the city of Hafr Al Batin is to improve storm drainage to control runoff generated within and around the urbanized area as described above and use quantitative precipitation estimates, such as IMERG, and/forecast to prepare in the case of exceptional rainfall events in the upper portions of Wadi Al Batin catchment.
The use of variable grid size resulted in significant saving of computing time. It is necessary to use a 30 × 30 m2 grid size in the urbanized portion of the catchment (about 300 km2) to make full use of the topographic resolution and include the major land surface in the simulations. However, there is no need to use such a resolution of the dominantly barren desert portion of the catchment. Employing this gridding scheme, the simulation run takes about 5% of what it would have taken had the 30 × 30 m2 grid been used for the entire catchment (more than 4000 km2). We did not perform a synchronized run of the three simulations in parallel. The computing time for such a run is effectively not more than the time to run the model on the urbanized catchment alone. In addition to significant saving of computing time, the same approach used in this study will allow merging observed discharge from upstream portions of a catchment with the simulated runoff downstream.
Using this approach, physically-based, fully-distributed models like the GSSHA can be run on large basins that include urban centers that need to be modeled at very high resolutions. Previous studies demonstrated the validity of model predictions in urban settings, including semi-arid catchments. We believe that using high-resolution calibrated satellite products and land features would result in reasonable estimates of the flood inundation. However, several issues need to be taken into account. For example, the division of sub-catchments needs special treatment to make sure that no significant over land flow is lost. Also, the inflow into a sub-catchment must have the same time step of the model simulations (typically minutes or seconds). Disaggregation of observed discharge data represents a challenge as it may introduce significant errors. In addition, inclusion of sub-surface flow from one sub-catchment to the next, if needed, would add more complexity to the hydrological connectivity.
The main limitations of the study are that the IMERG product is relatively new and has not been validated extensively and there was no streamflow data to validate the hydrographs. Yet, the authors are confident that the main conclusions regarding the spatial distribution of the flood inundation and the relative contributions of different parts of the catchment to the flood generation will hold unless the hydrological model and IMERG are of very poor quality, which they do not believe to be the case.

Acknowledgments

The first author was funded in part by the U.S. Army Research Office (Grant W912HZ-14-P-0160). This support is cordially acknowledged. The second author would like to acknowledge the support provided by the Deanship of Scientific Research (DSR) at King Fahd University of Petroleum and Minerals (KFUPM) for funding his work through project No. RG1305-1&2.

Author Contributions

Hatim O. Sharif and Muhammad Al-Zahrani designed the overall study. Almoutaz El Hassan downloaded the remote sensing products and prepared and performed the hydrologic model simulations with input from Hatim O. Sharif and Muhammad Al-Zahrani. Hatim O. Sharif conducted post-analysis of the model outputs and prepared the first draft. Muhammad Al-Zahrani reviewed and revised the manuscript. Hatim O. Sharif did the final overall proofreading of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Moel, H.; Aerts, J.C.J.H. Effect of uncertainty in land use, damage models and inundation depth on flood damage estimates. Nat. Hazards 2011, 58, 407–425. [Google Scholar] [CrossRef]
  2. Borga, M.; Stoffel, M.; Marchi, L.; Marra, F.; Jakob, M. Hydrogeomorphic response to extreme rainfall in headwater systems: Flash floods and debris flows. J. Hydrol. 2014, 518, 194–205. [Google Scholar] [CrossRef]
  3. Piaget, N.; Froidevaux, P.; Giannakaki, P.; Gierth, F.; Martius, O.; Riemer, M.; Wolf, G.; Grams, C.M. Dynamics of a local Alpine flooding event in October 2011: Moisture source and large-scale circulation. Q. J. R. Meteorol. Soc. 2015, 141, 1922–1937. [Google Scholar] [CrossRef]
  4. Prakash, S.; Mitra, A.K.; AghaKouchak, A.; Pai, D.S. Error characterization of TRMM Multisatellite Precipitation Analysis (TMPA-3B42) products over India for different seasons. J. Hydrol. 2015, 529, 1302–1312. [Google Scholar] [CrossRef]
  5. Siddique-E-Akbor, A.H.M.; Hossain, F.; Sikder, S.; Shum, C.K.; Tseng, S.; Yi, Y.; Turk, F.J.; Limaye, A. Satellite Precipitation Data–Driven Hydrological Modeling for Water Resources Management in the Ganges, Brahmaputra, and Meghna Basins. Earth Interact. 2014, 18, 1–25. [Google Scholar] [CrossRef]
  6. Wood, E.F.; Roundy, J.K.; Troy, T.J.; Van Beek, L.P.H.; Bierkens, M.F.; Blyth, E.; de Roo, A.; Döll, P.; Ek, M.; Famiglietti, J.; et al. Hyper-resolution global land surface modeling: Meeting a grand challenge for monitoring Earth’s terrestrial water. Water Resour. Res. 2011, 47, W05301. [Google Scholar] [CrossRef]
  7. Su, F.; Hong, Y.; Lettenmaier, D.P. Evaluation of TRMM Multisatellite Precipitation Analysis (TMPA) and its utility in hydrologic prediction in the La Plata basin. J. Hydrometeorol. 2008, 9, 622–640. [Google Scholar] [CrossRef]
  8. Vergara, H.J.; Hong, Y.; Gourley, J.J.; Anagnostou, E.N.; Maggioni, V.; Stampoulis, D.; Kirstetter, P.-E. Effects of Resolution of Satellite-Based Rainfall Estimates on Hydrologic Modeling Skill at Different Scales. J. Hydrometeorol. 2013, 15, 593–613. [Google Scholar] [CrossRef]
  9. Behrangi, A.; Khakbaz, B.; Jaw, T.C.; AghaKouchak, A.; Hsu, K.; Sorooshian, S. Hydrologic evaluation of satellite precipitation products over a mid-size basin. J. Hydrol. 2011, 397, 225–237. [Google Scholar] [CrossRef]
  10. Sorooshian, S.; Kuo-Lin, H.; Xiaogang, G.; Gupta, H.V.; Imam, B.; Braithwaite, D. Evaluation of PERSIANN system satellite-based estimates of tropical rainfall. Bull. Am. Meteorol. Soc. 2000, 81, 2035–2046. [Google Scholar] [CrossRef]
  11. Joyce, R.J.; Janowiak, J.E.; Arkin, P.A.; Xie, P. CMORPH: A method that produces global precipitation estimates from passive microwave and infrared data at high spatial and temporal resolution. J. Hydrometeorol. 2004, 5, 487–503. [Google Scholar] [CrossRef]
  12. 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. [Google Scholar] [CrossRef]
  13. Kubota, T.; Shige, S.; Hashizume, H.; Aonashi, K.; Takahashi, N.; Seto, S.; Takayabu, Y.N.; Ushio, T.; Nakagawa, K.; Iwanami, K.; et al. Global precipitation map using satellite-borne microwave radiometers by the GSMaP Project: Production and validation. IEEE Trans. Geosci. Remote Sens. 2007, 45, 2259–2275. [Google Scholar] [CrossRef]
  14. Stisen, S.; Sandholt, I. Evaluation of remote-sensing-based rainfall products through predictive capability in hydrologic runoff modelling. Hydrol. Process. 2010, 24, 879–891. [Google Scholar] [CrossRef]
  15. Nikolopoulos, E.I.; Anagnostou, E.N.; Borga, M. Using High-resolution satellite rainfall products to simulate a major flash flood event in Northern Italy. J. Hydrometeorol. 2013, 14, 171–185. [Google Scholar] [CrossRef]
  16. Yong, B.; Hong, Y.; Ren, L.-L.; Gourley, J.J.; Huffman, G.J.; Chen, X.; Wang, W.; Khan, S.I. Assessment of evolving TRMM-based multisatellite real-time precipitation estimation methods and their impacts on hydrologic prediction in a high latitude basin. J. Geophys. Res. 2012, 117, D09108. [Google Scholar] [CrossRef]
  17. Anquetin, S.; Yates, E.; Ducrocq, V.; Samouillan, S.; Chancibault, K.; Davolio, S.; Accadia, C.; Casaioli, M.; Mariani, S.; Ficca, G.; et al. The 8 and 9 September 2002 flash flood event in France: An intercomparison of operational and research meteorological models. Nat. Hazards Earth Syst. Sci. 2005, 5, 741–754. [Google Scholar] [CrossRef]
  18. Deng, L.; McCabe, M.F.; Stenchikov, G.; Evans, J.P.; Kucera, P.A. Simulation of flash-flood-producing storm events in Saudi Arabia using the weather research and forecasting model. J. Hydrometeor. 2015, 16, 615–630. [Google Scholar] [CrossRef]
  19. Gao, W.; Sui, C.-H. A modeling analysis of rainfall and water cycle by the cloud-resolving WRF Model over the western North Pacific. Adv. Atmos. Sci. 2013, 30, 1695–1711. [Google Scholar] [CrossRef]
  20. Xie, B.; Zhang, F. Impacts of typhoon track and island topography on the heavy rainfalls in Taiwan associated with Morakot (2009). Mon. Weather Rev. 2012, 140, 3379–3394. [Google Scholar] [CrossRef]
  21. Hou, A.Y.; Kakar, R.K.; Neeck, S.; Azarbarzin, A.A.; Kummerow, C.D.; Kojima, M.; Oki, R.; Nakamura, K.; Iguchi, T. The global precipitation measurement mission. Bull. Am. Meteorol. Soc. 2014, 95, 701–722. [Google Scholar] [CrossRef]
  22. GPM, 2016: Precipitation Measurement Missions. Global Precipitation Measurement Program, NASA. Available online: http://gpm.nasa.gov/ (accessed on 15 November 2016).
  23. Furl, C.; Sharif, H.O.; El Hassan, A.; Mazari, N.; Burtch, D.; Mullendore, G.L. Hydrometeorological Analysis of Tropical Storm Hermine and Central Texas Flash Flooding, September 2010. J. Hydrometeor. 2015, 16, 2311–2327. [Google Scholar] [CrossRef]
  24. Smith, M.B.; Koren, V.; Zhang, Z.; Zhang, Y.; Reed, S.M.; Cui, Z.; Moreda, F.; Cosgrove, B.A.; Mizukami, N.; Anderson, E.A. DMIP Participants Results of the DMIP 2 Oklahoma experiments. J. Hydrol. 2012, 418–419, 17–48. [Google Scholar] [CrossRef]
  25. Sharif, H.O.; Sparks, L.; Hassan, A.A.; Zeitler, X.W.J.; Xie, H. Application of a Distributed Hydrologic Model to the November 17, 2004 Flood of Bull Creek Watershed, Austin, Texas. J. Hydrol. Eng. 2010, 15, 651–657. [Google Scholar] [CrossRef]
  26. Sharif, H.O.; Al-Juaidi, F.H.; Al-Othman, A.; Al-Dousary, I.; Fadda, E.; Jamal-Uddeen, S.; Elhassan, A. Flood Hazards in an Urbanizing Watershed in Riyadh, Saudi Arabia. Geomat. Nat. Hazards Risk. 2014. [Google Scholar] [CrossRef]
  27. Koren, V.I.; Reed, S.; Smith, M.; Zhang, Z.; Seo, D.-J. Hydrology laboratory research modeling system (HL-RMS) of the US National Weather Service. J. Hydrol. 2004, 291, 297–318. [Google Scholar] [CrossRef]
  28. Downer, C.W.; Ogden, F.L. GSSHA: Model to simulate diverse stream flow producing processes. J. Hydrol. Eng. 2004, 9, 161–174. [Google Scholar] [CrossRef]
  29. Sharif, H.O.; Chintalapudi, S.; Elhassan, A.; Xie, H.; Zeitler, J. Physically-based hydrological modeling of the 2002 floods in San Antonio, Texas. J. Hydrol. Eng. 2013, 18, 228–236. [Google Scholar] [CrossRef]
  30. Chintalapudi, S.; Sharif, H.O.; Xie, H. Sensitivity of distributed hydrologic simulations to ground and satellite based rainfall. Water 2014, 6, 1221–1245. [Google Scholar] [CrossRef]
  31. Ogden, F.L.; Downer, C.W.; Meselhe, E. U.S. Army Corps of Engineers Gridded Surface/Subsurface Hydrologic Analysis (GSSHA) Model: Distributed-Parameter, Physically Based Watershed Simulations. In Proceedings of the World Water and Environmental Resources Congress 2003, Philadelphia, PA, USA, 23–26 June 2003; Bizier, P., DeBarry, P., Eds.; American Society of Civil Engineers: Washington, DC, USA, 2003. [Google Scholar]
  32. Elhassan, A.; Sharif, H.O.; Jackson, T.; Chintalapudi, S. Performance of a Conceptual and a Physically-Based Model in Simulating The Response of a Semi-Urbanized Watershed in San Antonio, Texas. Hydrol. Process. 2013, 27, 3394–3408. [Google Scholar] [CrossRef]
  33. Sharif, H.O.; Hassan, A.A.; Bin-Shafique, S.; Xie, H.; Zeitler, J. Hydrologic Modeling of an Extreme Flood in the Guadalupe River in Texas. J. Am. Water Resour. Assoc. 2010, 46, 881–891. [Google Scholar] [CrossRef]
  34. Sharif, H.O.; Yates, D.; Roberts, R.; Mueller, C. The use of an automated now-casting system to forecast flash floods in an urban watershed. J. Hydrometeorol. 2006, 7, 190–202. [Google Scholar] [CrossRef]
  35. Furl, C.; Sharif, H.O.; Elhassan, M.A.A.; Mazari, N. Precipitation Amount and Intensity Trends Across Southwest Saudi Arabia. J. Am. Water Resour. Assoc. 2014, 50, 74–82. [Google Scholar] [CrossRef]
  36. Almazroui, M.A.; Islam, M.N.; Jones, P.D.; Athar, H.; Rahman, M.A. Recent climate change in the Arabian Peninsula: Seasonal rainfall and temperature climatology of Saudi Arabia for 1979–2009. Atmos. Res. 2012, 111, 29–45. [Google Scholar] [CrossRef]
  37. Tekeli, A.E.; Fouli, H. Evaluation of TRMM satellite-based precipitation indexes for flood forecasting over Riyadh City, Saudi Arabia. J. Hydrol. 2016, 541 Pt A, 471–479. [Google Scholar] [CrossRef]
  38. El-Baz, F. Remote Sensing: Generating Knowledge about Groundwater. Arab Environment: Water. 2010 Report of the Arab Forum for Environment and Development (AFED). 2010. Available online: http://www.afedonline.org/Report2010/pdf/AR/prefacear.pdf. (accessed on 12 January 2017).
  39. Edgell, H.S. Aquifers of Saudi Arabia and their geological framework. Arab. J. Sci. Eng. Water Resour. Arab. Penins. Part I 1997, 22, 3–31. [Google Scholar]
  40. McClure, H.A. Ar Rub’ Al Khali. In Quaternary Period in Saudi Arabia; Springer: Vienna, Austria, 1978; Volume 1, pp. 252–263. [Google Scholar]
  41. El Hassan, I.M.; Algarni, D.A.; Dalbouh, F.M. Flood risk prediction using DEM and GIS as applied to Wijj Valley, Taif, Saudi Arabia. J. Geomat. 2014, 8, 86–89. [Google Scholar]
  42. EMRL (Environmental Modeling Research Laboratory). Watershed Modeling System (WMS) Version 10.1 Tutorial; Brigham Young University: Provo, UT, USA, 2016. [Google Scholar]
  43. FAO (Food and Agriculture Organization of the United Nations). World Reference Base for Soil Resources. World Soil Resources Report, #84; FAO: Rome, Italy, 1998; p. 88. [Google Scholar]
  44. Ogden, F.L.; Saghafian, B. Green and Ampt infiltration with redistribution. J. Irrig. Drain. Eng. 1997, 123, 386–393. [Google Scholar] [CrossRef]
  45. GSSHA Primer, 2012: GSSHA Model Description and Formulation. Available online: http://gsshawiki.com/gssha/GSSHA_Primer#Description (accessed on 18 November 2016).
  46. Gsshawiki, 2016: Information of GSSHA Model. Available online: http://gsshawiki.com/gssha/Main_Page (accessed on 18 November 2016).
  47. Chintalapudi, S.; Sharif, H.O.; Furl, C. High-resolution, Fully-distributed Event-based Hydrologic Simulations over a Large Watershed in Texas. Arab. J. Sci. Eng. 2016, in press. [Google Scholar] [CrossRef]
  48. Chow, V.T. A general formula for hydrologic frequency analysis. Trans. Am. Geophys. Union 1951, 32, 231–237. [Google Scholar] [CrossRef]
  49. World Meteorological Organization (WMO). Manual on Estimation of Probable Maximum Precipitation; Oper. Hydrol. Rep., 1, WMO Pap. 332; WMO: Geneva, Switzerland, 1986. [Google Scholar]
  50. McKay, G.A. Statistical Estimates of Precipitation Extremes for the Prairie Provinces; Canada Department of Agriculture, Prairie Farm Rehabilitation Administration (PFRA) Engineering Branch: Regina, SK, Canada, 1965. [Google Scholar]
  51. SCS (Soil Conservation Service, U.S. Department of Agriculture). A Method for Estimating Volume and Runoff in Small Watersheds, TP-149; US Department of Agriculture, Soil Conservation Service: Washington, DC, USA, 1973.
Figure 1. Location of the Wadi Al Batin catchment in Northeast Saudi Arabia.
Figure 1. Location of the Wadi Al Batin catchment in Northeast Saudi Arabia.
Water 09 00163 g001
Figure 2. Soil and land use types of the Wadi Al Batin catchment in Northeast Saudi Arabia.
Figure 2. Soil and land use types of the Wadi Al Batin catchment in Northeast Saudi Arabia.
Water 09 00163 g002
Figure 3. The city of Hafr Al Batin and Wadi Al Batin and its major tributaries.
Figure 3. The city of Hafr Al Batin and Wadi Al Batin and its major tributaries.
Water 09 00163 g003
Figure 4. IMERG rainfall totals for the 2–3 November 2015 storm over the region.
Figure 4. IMERG rainfall totals for the 2–3 November 2015 storm over the region.
Water 09 00163 g004
Figure 5. A simplified flowchart of GSSHA model.
Figure 5. A simplified flowchart of GSSHA model.
Water 09 00163 g005
Figure 6. The three sub-catchments of the Wadi Al Batin with variable model grid sizes.
Figure 6. The three sub-catchments of the Wadi Al Batin with variable model grid sizes.
Water 09 00163 g006
Figure 7. Rainfall totals for the 2–3 November 2015 storm as estimated by the three IMERG satellite products using IDW interpolation.
Figure 7. Rainfall totals for the 2–3 November 2015 storm as estimated by the three IMERG satellite products using IDW interpolation.
Water 09 00163 g007
Figure 8. 30-min hyetographs of the the three IMERG products for the 2–3 November 2015 storm.
Figure 8. 30-min hyetographs of the the three IMERG products for the 2–3 November 2015 storm.
Water 09 00163 g008
Figure 9. Outlet discharge predicted by GSSHA model for 2–3 November 2015 storm as estimated by the three GPM satellite products.
Figure 9. Outlet discharge predicted by GSSHA model for 2–3 November 2015 storm as estimated by the three GPM satellite products.
Water 09 00163 g009
Figure 10. Flood inundation in the city of Hafr Al Batin caused by the 2–3 November 2015 storm as estimated by the GSSHA model.
Figure 10. Flood inundation in the city of Hafr Al Batin caused by the 2–3 November 2015 storm as estimated by the GSSHA model.
Water 09 00163 g010
Figure 11. Hourly hyeographs of the design storms used for addition GSSHA model simulations.
Figure 11. Hourly hyeographs of the design storms used for addition GSSHA model simulations.
Water 09 00163 g011

Share and Cite

MDPI and ACS Style

Sharif, H.O.; Al-Zahrani, M.; Hassan, A.E. Physically, Fully-Distributed Hydrologic Simulations Driven by GPM Satellite Rainfall over an Urbanizing Arid Catchment in Saudi Arabia. Water 2017, 9, 163. https://doi.org/10.3390/w9030163

AMA Style

Sharif HO, Al-Zahrani M, Hassan AE. Physically, Fully-Distributed Hydrologic Simulations Driven by GPM Satellite Rainfall over an Urbanizing Arid Catchment in Saudi Arabia. Water. 2017; 9(3):163. https://doi.org/10.3390/w9030163

Chicago/Turabian Style

Sharif, Hatim O., Muhammad Al-Zahrani, and Almoutaz El Hassan. 2017. "Physically, Fully-Distributed Hydrologic Simulations Driven by GPM Satellite Rainfall over an Urbanizing Arid Catchment in Saudi Arabia" Water 9, no. 3: 163. https://doi.org/10.3390/w9030163

APA Style

Sharif, H. O., Al-Zahrani, M., & Hassan, A. E. (2017). Physically, Fully-Distributed Hydrologic Simulations Driven by GPM Satellite Rainfall over an Urbanizing Arid Catchment in Saudi Arabia. Water, 9(3), 163. https://doi.org/10.3390/w9030163

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop