Determination of Potential Runoff Coefficient Using GIS and Remote Sensing

Flash floods in arid environments are a major hazard feature to human and to the infrastructure. Shortage of accurate environmental data is main reason for inaccurate prediction of flash flooding characteristics. The curve number (CN) is a hydrologic number used to describe the storm water runoff potential for drainage area. This study introduces an approach to determine runoff coefficient in Jeddah, Saudi Arabia using remote sensing and GIS. Remote sensing and geographic information system techniques were used to obtain and prepare input data for hydrologic model. The land cover map was derived using maximum likelihood classification of a SPOT image. The soil properties (texture and permeability) were derived using the soil maps published my ministry of water and agriculture in Saudi Arabia. These soil parameters were used to classify the soil map into hydrological soil groups (HSG). Using the derived information within the hydrological modelling system, the runoff depth was predicted for an assumed severe storm scenario. The advantages of the proposed approach are simplicity, less input data, one software used for all steps, and its ability to be applied for any site. The results show that the runoff depth is directly proportional to runoff coefficient and the total volume of runoff is more than 136 million cubic meters for a rainfall of 103.6 mm.


Introduction
In arid zones, there is a shortage of data needed for hydrological processes as mentioned by [1].In these zones flash flooding happened suddenly and affects R. Khalil both human and infrastructure.The effective way to reduce the damages due to flood is to predict the critical sites affected by flash flood for management planes.The part of rainfall that turns to runoff due to land use and soil hydrological parameters is defined as runoff coefficient [2].It may also be defined as the ratio between the runoff depth and the rainfall depth [3].The reasonable calculation of runoff from rainfall is the key for flood estimation [4] [5].Defining rainfall-runoff relationship leads to calculate flood properties such as runoff depth, peak discharge, runoff speed, and runoff volume.Flood properties derivation is important for planners and decision makers to avoid flood hazard effects.When flood discharge data records are missed, there are many approaches to estimate runoff depth and volume based on soil and surface properties of catchment area.These approaches may be summarized into three groups, simple, moderate and complex models as stated by [4].The curve number (CN) is one of the moderate complicity models and it is widely used for flood estimation.The CN model was developed by United States Department of Agriculture (USDA), National Resources Conservation Service (NRSC) in 1969 [6].It is an empirical model with clearly stated assumptions and few data requirements [7].CN is stable conceptual method for predicting of direct runoff depth using storm rainfall intensity, land use and soil hydrological properties of a catchment.
Remote sensing imageries considered a major source of spatial data especially for wide area or when the ground surveys are not available.They can be used to generate land use, soil and geological maps needed for flood estimation.Geographic Information System (GIS) is a powerful tool in hydrological modelling because of its capability to handle large amount of spatial and attribute data.Delineation of hydrological catchments, map overlay and analysis, which are basics of GIS software, help for derivation and aggregation of hydrologic parameters from the input data like DEM, soil map, land use map, and rainfall data.Remote sensing (RS) images and Geographic Information System (GIS) techniques have been used in flood hazard studies by many researchers, e.g.[8]- [18].
The main objective of this research article is determining the potential runoff coefficient for Jeddah, Saudi Arabia by applying CN model and using remote sensing and GIS techniques.

Study Area
Jeddah is the second main city in Saudi Arabia.It is at the middle of east coast of the Red Sea and represents the important commercial port.Its population exceeds 3.4 million people according to (Central Department of Statistics & Information 2010) with 3.5% annual rate of growth make its population exceeds 4.2 million people.Its weather condition is classified as hot.A rainfall intensity of 80 and 124 mm/day hit Jeddah in November 2009 and January 2011, respectively as mentioned by [19].Reference [20] mentioned that these occasional floods were examples of flash flood striking which characterized by short durations and destructive results.The watershed 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).

Data and Methodology
The available rainfall data for Jeddah are from two rain gauge stations, J134 which operated by the Ministry of Water and Electricity (MOWE) and 41024 (airport station) which operated by Presidency of Meteorology and Environment (JMPE) [21].The data represent the recorded observation for a period of 42 years extending from January 1971 to December 2012; no data were available for the most recent years.The maximum annual daily rainfall for the two stations is shown in (Figure 2).High resolution satellite image SPOT5 (2.5 m multispectral) was obtained from King Abdulaziz City for Science and Technology (KACST) and used generate land use map.Topography data was obtained using Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) and used to extract slope data and stream network.Soil data were obtained using Saudi general soil map, ministry of agricultural and water 1986.
The methodology used to determine the potential runoff coefficient for the study area using remote sensing and GIS is shown in the flow chart in (Figure 3).It starts with scanning and Georeferencing the soil maps using ArcMap 10.4, the soil polygons were digitized, and each soil type was assigned to a Hydrological Soil Group (HSG) according to [22].Image classification tool in ArcMap 10.4 was used to classify and convert the SPOT image to land use thematic map by applying the supervised classification technique.The DEM was used to generate slope map using surface analysis tool, and to generate stream network and basin parameters using HEC-GeoHMS.The runoff curve number lookup table was build according to the land use map and [22].All data were then overlaid in  ArcGIS environment to generate the runoff coefficient.Lastly, the flood formulas were applied to calculate the flood parameters using field calculator in Arc-Map 10.4.

Land Use Analysis
SPOT image (2.5 m) for the year 2010 was obtained from the King Abdul-Aziz City for Science and Technology (KACST).Supervised classification was performed to the image using spectral signatures collected from training samples (polygons that represent distinct sample areas of the different land cover types to be classified).After preparing the signature file, Maximum likelihood classifier then attached labels to all the image pixels according to the trained parameters to generate land cover map.The land cover map was classified into six main classes: water, vegetation, rocks, bare soil, and build up area as shown in (Figure 4).

Soil Type Analysis
The general soil map of Saudi Arabia was the source for soil data used in this

Slope Analysis
The runoff velocity and soil erosion are in direct proportional to the degree of

Potential Runoff Coefficient
As flash flood in Jeddah occurred with rainfall didn't exceed 80 mm and rainstorm duration didn't exceed 3 hours, it is important to estimate the expected flood discharges of return periods reliably for protection and future developments [23].
The rainfall data was used to predict the rainfall depth at different return periods.
Among many of the probability density functions, numerous studies recommended the Gumbel or Extreme Value Type I (EV1) function as it demonstrated the best fit in most cases and providing the best prediction of rainfall depth.The estimated rainfall depth for different return periods for the two rain gauge stations for JPME (41024) and J134 are shown in Table 1.
The flood hazard parameters such as watershed storage, runoff depth, and volume of runoff were calculated.The CN method Equation (1) was used for runoff depth calculation [22].
where Q: the direct runoff depth (mm), Figure 6.Slope map.
Table 1.Estimated daily rainfall (mm) for different return periods.P: depth of precipitation for a specific return period (mm), S: the watershed storage (mm), and can be calculated using Equation (2), CN: the curve number.
The volume of the runoff was calculated for subbasins to show the local effect of runoff and was also calculated for the major basins to show the total volume of runoff flood.The volume of the runoff can be calculated using Equation (3) where V Q : Volume of runoff (m 3 ), Q: Depth of runoff (m), A: Basin area (m 2 ).

Results
The flood characteristics were calculated using the rainfall depth (P) equals 106.3 mm for a return period of 100 years as shown in Table 1.These characteristics include area, longest flow length and the runoff volume of the main basins in Jeddah watershed area as shown in Table 2.The results show that there are 7 major and 5 minor basins in Jeddah as shown in (Figure 7).The areas of the  Refereeing to (Table 2), one can notice that the runoff depth (Q) is directly proportional to the curve number and inversely proportional to slope and the curve number has the great effect on it.
The proposed approach of using remote sensing and GIS applying the CN method for runoff coeffecient estimation has many advantages over other approaches.Firstly, it uses one software to perform all procedure steps.Secondly, only satellite image, soil maps and DEM are needed to calculate the runoff parameters.Thirdly, all needed calculations are done within the GIS environment using field calculation.Fourthly, it can be modeled using model builder so, runoff parameters estimation process can be efficient, faster, and easily performed for several return period scenarios and for any regions.

Conclusion
This research article presented an efficient approach to accurate determination of potential runoff coefficient in Jeddah city using remote sensing and GIS.The effects of land use, soil hydrological characteristics, surface slope, were considered in calculating runoff coefficient and consequently runoff depth and runoff volume.The results of the research show that the total runoff volume for a rainfall depth of 106.3 mm is 136.5 million m 3 .Results also show that the main factors affect the total flood volumes, are the basin area, and the flow length.Additionally, it has been concluded that the higher CN value and slope percent, the higher runoff and flood hazards.

Figure 1 .
Figure 1.Location of the study area.

Figure 2 .
Figure 2. Rainfall data for the two rain gauges.
study.The general soil map was scanned and georeferenced then different types of soil were digitalized in ArcMap 10.4.Each type of soil was assigned to the proper hydrological soil group that developed by US Soil Conservation Service (SCS).According to SCS, all the soil, based on their permeability and infiltration are divided to 4 groups A, B, C and D. Hydrological Soil Group of A has the lowest runoff potential while group D has the highest.The soil map was then classified into Hydrological Soil Groups (HSGs) as shown in Figure 5.

Table 2 .
Morphometric and flood parameters of basins.
9 to 15.5 square kilometers, and their longest flow paths range from 7.48 to 10.23 kilometers.The total runoff volume was calculated and found to be more than 136 million cubic meters.R. Khalil DOI: 10.4236/jgis.2017.96046760 Journal of Geographic Information System