The destructive effect of flash flood on properties and road networks make it one of the most dangerous natural hazards. Adopting an appropriate approach beside acquiring land cover data and related hydrological parameters is the way to predict sites sensitive to flooding  . The calculation of runoff coefficient which is the part of rainfall that turns to runoff as defined by  is the corner stone for flood estimation which affects many related activities such as flood damage determination, flood control plans and structures, water resources management and flood risk assessment  and  . There are several methods to estimate the flood parameters based on the morphometric characteristics of flood basins. Reference  categorized these methods according to their complicity into three groups simple, moderate and complex. Reference  summarized five methods for quantitative analysis of drainage basin morphometry. These methods are El-Shamy’s approach, the morphometric ranking method, Wahid’s flash flood rating method, Morphometric hazard degree assessment method and the curve number (CN) method. The curve number (CN) method was applied in this research as it requires few data, its complicity is moderate, and it is widely used for flood estimation  . CN method uses storm rainfall depth, land use and soil infiltration properties of a catchment for predicting of direct runoff depth. Remote sensing imageries are main spatial data source to produce land use and soil maps needed for flood properties estimation. Delineation of hydrological catchments, map overlay, and data analysis make GIS a powerful tool in hydrological modelling. For these reasons many studies applied GIS and remote sensing in flood calculations e.g.  -  . The use of Digital Elevation Model (DEM) facilitated the procedure of watershed delineation. The DEM provides the basic input for the delineation of the catchment area and generation of drainage networks in GIS environment. Many researches attempted to configure a relationship between basin morphometric parameters and flood to produce a flood hazard risk map. References  ,   and  ranked subbasins within a basin according to a series (more than 20 in some researches) of calculated morphometric parameters in order to produce a flood hazard risk map. The raw value for each parameter was normalized to ranged from (0 to 1) to minimize the value of total flood risk after summation of all the parameters. Reference  used only 5 morphometric parameters to produce the flood risk map. Other researchers used only one morphometric parameter to produce the flood risk map e.g.  used the runoff volume,  used flow velocity,  suggested using travel time,  acknowledged peak flow magnitude as a primary flood hazard indicator, and  used the peak flow rate as a flood hazard indicator. In this study a Flood Risk Code (FRC) was proposed to map the flood hazard zones in Jeddah. Two effective parameters were used to generate the FRC. These parameters are time of concentration (Tc) and runoff depth (Q). The Tc has functions of flow length, slope and curve number which has functions of land use and soil properties. The Q has functions of rainfall, slope and curve number. Using these two factors (Tc and Q) is exactly as using rainfall, land use, soil, slope and flow length for computing FRC. Each parameter was ranked by how risky it is on a scale of 1 to 4, with 4 being the extreme risk. The overlay analysis technique was used to generate the flood risk code map.
2. Study Area
The proposed approach was applied on Jeddah city, Saudi Arabia. Jeddah lies at the middle of west region of Saudi Arabia and represents the important commercial port. It’s surrounded by a chain of mountains at the east side from which the flood begins and flow through the city to the sea. Its population exceeds 3.4 million people live in an area of about 5460 km2. Although its weather condition is hot, it receives a rain storm from time to time. In November 2009 and January 2011, Jeddah was hit by a rainfall storm of 80 and 124 mm/day respectively as mentioned by  . These were examples of flash flood which strike occasionally and characterized by short durations and harmful results as mentioned by  . The watershed area that affected the city is located between (21˚15'N, 21˚50'N) and (39˚0'E, 39˚35'E) as shown in Figure 1.
The data set used in this study contain a SPOT satellite multispectral image of 2.5 m resolution used to generate the land use map as shown in Figure 2, ASTER DEM 30 m spatial resolution is provided on line cost free to extract slope data and stream network, soil maps used to generate the soil data as shown in Figure 3 and the rainfall data for a period of 42 years extending from January 1971 to December 2012 for Jeddah from two rain gauge stations, J134 and JMPE (airport station). The maximum annual daily rainfall for the two stations is shown in Figure 4.
The proposed GIS and RS based approach for flood hazard mapping consists of several steps. These steps are performed in ArcGIS environment using the proper tools for each. The sequences of these steps are schematically shown in Figure 5 and can be summarized as following:
Figure 1. The study site and its catchment area.
Figure 2. Land use map for the study area generated from classification of SPOT image.
Figure 3. Soil hydrological group map for the study area.
Figure 4. The maximum annual daily rainfall (mm) recorded at JMPE and J134 stations from 1971 to 2012.
Figure 5. Methodology for Flood hazard mapping based on GIS and RS.
・ Scanning and Georeferencing the soil maps. The soil polygons were digitized, and the soil types were extracted from the maps. Soil types was categorized to Hydrological Soil Groups (HSG) according to  .
・ Georeferencing, segmentation of SPOT image then the supervised classification technique was applied to convert the image to land use thematic map. The thematic map then converted to polygons.
・ Generate slope map using DEM data
・ Build the runoff curve number lookup table according to the land use map and reference  .
・ Generating curve number grid (CN Grid) using the former generated layers using HEC-GeoHMS tool.
・ Generating the geomorphological and hydrological parameters such as the drainage network, the main basins and the subbasins of each main catchment, and the longest flow path in each catchment using HEC-GeoHMS tool.
・ The rainfall intensity for a 100-years return period has been estimated as 106.3 mm and is used in the current research study. This estimation was based on the rainfall data of two rainfall gauges in Jeddah as shown in 5.1.
・ Apply flood formulas to calculate runoff depth, time of concentration and runoff volume for subbasins.
・ Reclassify the runoff depth and time of concentration values and rank them on scale of 1 to 4.
・ Integrating all data in a GIS environment to generate flood risk code and flood maps.
・ Evaluate the proposed approach using ground truth data.
5. Data Processing
5.1. Rainfall-Runoff Analysis
Among many of the probability density functions that used to estimate the rainfall depth at a specific return period, the Gumbel or Extreme Value Type I (EV1) function was recommended as it provided the best prediction  . In this function the magnitude XT of a hydrologic event may be calculated using Chow Equation (1) 
μ: the mean value of the data
σ: the standard deviation
KT: the frequency factor
T: the return period
The frequency factor can be calculated using Equation (2) 
The estimated rainfall depth for different return periods for the two rain gauge stations are shown in Table 1.
5.2. Catchment Morphometric Parameters
The curve number (CN) method was used to estimate the run off over the study area. The CN is an indicator of infiltration capacity for each soil type and is based on land use, hydrologic soil group and surface slope. The details of generating land use map through classification of SPOT image, generating the hydrologic soil group map and slope map for the study area can be found in  . The morphometric parameters are calculated at the sub-basin level. These parameters include the main basins, sub-basins, sub-basin-CN, time of concentration and the longest flow path for each sub-basin were generated in ArcGIS using HEC-GeoHMS tool.
Table 1. Estimated daily rainfall (mm) for different return periods.
5.3. Flood Hazard Parameters
The flood hazard parameters such as watershed storage, runoff depth, and volume of runoff were calculated using ArcGIS field calculation capability. The flood hazard parameters were calculated according to the CN method equations published in Technical Report 55  . Equation (3) used to calculate watershed storage, Equation (4) to calculate the runoff depth while Equation (5) used to calculate runoff volume for subbasins to demonstrate the local effect of runoff.
S: the watershed storage (mm),
Q: the direct runoff depth (mm),
P: rainfall depth for a specific return period (mm),
CN: the curve number,
VQ: Volume of runoff (m3),
A: Basin area (m2).
The time of concentration is the time required for a water drop to travel from the watershed divide (the most distant point of the watershed) to the watershed outlet. There are several methods to calculate the concentration time. The CN method was used to estimate the time concentration. Equation (6)  used to estimate the time of concentration as follows:
Tc = time of concentration, hours
L= longest flow path, m
Y = average watershed slope (percent)
CN = Curve number
6. Results & Discussion
The proposed approach of using RS and GIS for estimating flood hazard parameters was applied on Jeddah watershed affected area. The results show that there are 12 basins in Jeddah 7 of them are major and 5 are minor as shown in Figure 6. The minimum area of the major basins was 59.04 km2 with a flow path of 12.7 km while the maximum area was 555.5 km2 with a flow path of 77.4 km. Their time of concentration ranges from 5.87 to 13.80 hours. The areas of the minor basins range from 6.9 to 15.5 square kilometers, their longest flow paths range from 7.48 to 10.23 kilometers, and their time of concentration of the major basins ranges from 2.81 to 4.14 hours. Table 2 presents some geometrical and flood parameters including the runoff volume and the time of concentration of the 12 basins. A rainfall depth (P) equals 106.3 mm for a return period of 100 years as shown in Table 1 was used to calculate the flood characteristics.
The results presented in Table 2 show that the time of concentration (Tc) increased with increasing the flow length and decreasing the surface slope. The curve number and the surface slope have the great effect on time of concentration (Tc). The runoff depth (Q) is directly proportional to the curve number and inversely proportional to slope. It is affected greatly by the curve number value.
In this study a Flood Risk Code (FRC) was proposed to map the flood hazard zones in Jeddah. The two effective parameters Tc and Q were used to generate the FRC. The Tc has functions of flow length, slope and curve number which has functions to land use and soil properties. The Q has functions of rainfall, slope and curve number. Using these two factors (Tc and Q) is exactly as using rainfall, land use, soil, slope and flow length for computing FRC.
The Tc and Q were computed for all subbasins in the study area. The Tc values range from 0.26 to 5.80 hours and the Q values range from 43.7 to 94.5 mm. The Tc and Q were ranked by how risky they are on a scale of 1 to 4, with 4 being the extreme risk. The rank value increases with increasing the Q and decreasing the Tc according to the limits shown in Table 3.
The FRC is generated using overlay analysis in ArcMap 10.4.1 software. The values of FRC were classified into Minor, moderate, Major and Extreme risk and used for mapping the Flash Flood Hazard zones as shown in Figure 7.
Table 2. Geometrical and flood parameters of basins.
Table 3. Limit values of ranks.
Figure 6. Basins and main streams.
Figure 7. Flood Risk Code.
The areas at risk were obtained by overlaying the flood risk code map with the districts layer. The area and population affected by flood were summarized in Table 4 and may be used for and floodplain development control. More than 850,000 capita in 24 districts are vulnerable to extreme flood risk. About 2,320,000 in 63 districts are vulnerable to major risk. More than 560,000 in 18 districts are vulnerable to moderate flood risk. The rest of districts are vulnerable to minor flood risk or no data available about them. Jeddah area has received a rainstorm of 68 mm/day (max. value) in 21 November 2017 and the municipality of Jeddah published a map showing the road opening status. This map was used to verify the proposed approach by overlying it with the generated flood risk code map as shown in Figure 8. The figure proof the validity of the proposed approach as the closed roads due to flood (red lines) are coincide with major and extreme risk code, roads vulnerable to be closed (orange lines) are coincide with major risk code, the save roads (green lines) are coincide with moderate and minor risk codes.
Table 4. Area and population at flash flood risk.
Figure 8. Flood Risk Code verification.
This research article presented an efficient approach to accurately mapping the flood hazard areas in Jeddah city. The approach utilizes the use of remote sensing and GIS along with rainfall data to calculate the flood and basin morphometric parameters. The approach proposed uses two simple morphometric parameters (runoff depth and time of concentration) to generate a new Flood Risk Code for mapping hazard zones. By overlaying the flood risk code map with the districts layer, the hazard zones are delineated. The total runoff volume is 136.5 million m3 for a rainfall depth of 106.3 mm. The results show that around 1 million of Jeddah population live in 24 districts and an area of about 200 km2 are at extreme risk, more than 2 million of populations in 63 districts of area greater than 480 km2 are prone to major flood risk, the minority of Jeddah population are vulnerable to moderate to minor risk. The approach was verified using an actual map of flood hazard area due to 2017 rainfall storm and show its validity through the coincide of actual flood hazard classes with the obtained Flood Risk Code. The advantage of the approach is its ability to produce results using only topographical data and a satellite image, so it can be used in ungauged catchments, which is the case for most of the flood hazard regions in Saudi Arabia.