The aim of this study is to build a flood simulation model for the city of Hafr Al-Batin catchment area. The model consists of a hydrologic model, a hydraulic model and pre- and post-processing tools. A geographic information system-based modelling interface was used to pre-process the terrain and rainfall data and generate input files for hydrologic and hydraulic models. Soil map data, land cover map, land use map and digital elevation model was used to delineate the physical watershed's characteristics. The runoff estimation was based on the widely known Soil Conservation Service curve number approach. The hydrologic/hydraulic model simulations simulated the runoff hydrograph corresponding to different design storms and helped to delineate the resulting flood inundation maps. The results indicate good agreement between the delineated inundation map and the hazard map developed by the municipality. They also show that the location of the city complicates the runoff response for small storms creating two distinct peaks. The results of this study can be utilized for planning purposes and in the design of flood control structures as the study has estimated the runoff corresponding to different design storms and used hydraulics and geospatial data in delineating the flood zones.

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 subcatchments 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.

Population growth and land use modification in urban areas require the use of accurate tools for rainfall-runoff modeling, especially where the topography is complex. The recent improvement in the quality and resolution of remotely sensed precipitation satisfies a major need for such tools. A physically-based, fully distributed hydrologic model and a conceptual semi-distributed model, forced by satellite rainfall estimates, were used to simulate flooding events in a very arid, rapidly urbanizing watershed in Saudi Arabia. Observed peak discharge for two flood events was used to compare hydrographs simulated by the two models, one for calibration and one for validation. To further explore the effect of watershed heterogeneity, the hydrographs produced by three implementations of the conceptual were compared against each other and against the output of the physically-based model. The results showed the ability of the distributed models to capture the effect of the complex topography and variability of land use and soils of the watershed. In general, the GSSHA model required less calibration and performed better than HEC-HMS. This study confirms that the semi-distributed HEC-HMS model cannot be used without calibration, while the GSSHA model can be the best option in the case of a lack of data. Although the two models showed good agreement at the calibration point, there were significant differences in the runoff, discharge, and infiltration values at interior points of the watershed.

Rainfall is the most important input parameter for water resource planning and hydrological studies because flood risk assessment, rainfall harvesting and runoff estimation depend on the rainfall distribution within a region. Due to practical and economic factors, it is not possible to site rainfall stations everywhere, so representative rainfall stations are sited at specific locations. Rainfall distribution is then estimated from such stations. In this study, rainfall distribution in the southwestern region of Saudi Arabia was estimated using kriging, co-kriging and inverse distance weighted (IDW) methods. Historical records of rainfall from 47 stations for the period 1965–2010 and the altitude of these stations were used. The study shows that co-kriging is a better interpolator than the kriging and IDW methods, with a better correlation between actual and estimated monthly average rainfall for the region.

The purpose of this study is to understand flooding in a 1,725-km² arid catchment in Makkah Province, Saudi Arabia, with very limited ground observations. This mountainous catchment includes densely urbanized areas built on lowlands and in the ephemeral streams of flood zones, which makes them vulnerable to flooding. A physically based, fully distributed hydrologic model was produced by three Integrated Multi-satellite Retrievals for Global Precipitation Measurement Mission (IMERG) high-resolution satellite rainfall products, together with a limited number of ground observations to simulate recent flood events. The hydrologic model was calibrated and validated using two events that occurred in 2010 and 2018, respectively. Details of the flooding in the urbanized areas were examined through simulation of a third event. Significant differences were noticed when the three IMERG rainfall products (Early, Late, and Final) were compared in terms of the spatial patterns and total half-hour accumulation for the three events. Accordingly, a method was developed to remove the biases of IMERG rainfall estimates using ground observations. The Early product was the closest to ground observations, producing peak discharges of 428 m³/s and 663.6 m³/s for calibration and validation storm events, respectively, and resulted in better estimates of the observed peak discharge even before adjustment. The Final run product significantly underestimated rainfall, as peak discharges of 204 m³/s and 123 m³/s for calibration and validation storm events, respectively, resulted in an amplified underestimation of the runoff.