Received 7 March 2016; accepted 26 April 2016; published 29 April 2016
The vulnerable and stressful conditions of arid southern Jordan, including Wadi Yutum watershed and the lower Wadi Araba, impose careful site selection for urban and rural land use development to assure sustainability  -  . The long intervals between high-magnitude, low frequency flood events, and the lack of hydrological data have often led the planners to ignore, or pay little attention to the fluvial processes and the associated geomorphic hazards, i.e., flash floods, gully erosion and abundant sediment supply. Over the last four decades, there have been a notable number of applied geomorphic studies carried out in arid areas worldwide, thus enhancing geomorphic techniques which can be adopted to avoid destructive hazards, and to improve levels of awareness regarding floods geomorphology and the role of applied geomorphic surveys and information in planning processes   -  . Past experience from arid coastal cities such as Aqaba and other comparable cities, i.e. Eilat at the opposite corner of the Gulf of Aqaba   , Gulf of Suez, Red Sea, and Sinai Peninsula, Egypt  -  indicates that flash floods are considered the most recurrent disaster recorded over the last three decades (1981-2011), not only in drylands of Jordan and Egypt, but in the entire Middle East  . Unplanned urban growth has frequently resulted in heavy disturbance of land surface and natural drainage. In the light of poor perception of floods hydrology and the associated fluvial processes, and the prominent lack of adequate and efficient hydraulic structures, the Aqaba area was exposed to additional serious floods initiated by Wadi Yutum, and other wadis descending from the granite mountains overlooking lower Wadi Araba and the Gulf of Aqaba. Recent flood monitoring  demonstrates that “a small flood occurred every few years, more rarely a large one, and perhaps a truly catastrophic flood once in a few decades”. However, the hydro-geomorphic characteristics of flash floods in the Aqaba area are poorly understood due to the lack of long-term monitoring of fully instrumented arid watersheds and paleoflood studies as compared with Eilat on the opposite NW corner of the Gulf of Aqaba. Furthermore, long-term hydro-meteorological data are not available. Therefore, due to difficulties in predicting flash floods, planning adjustment to flooding risk, emergency and evacuation plans are of high priority, where Aqaba city has been almost disrupted by a flash flood of low or medium magnitude every two years in the last past decade   . Accordingly, the spatial distribution of land use categories, the expanding of residential areas across wadi floors, economic and life losses are expected to increase due to continuous exposure of assets and population to risk. An increase in frequency and intensity of heavy rainfall, and probably certain effects of climatic change  -  , will cause an increase in flash floods occurrences. The residential areas, the airport, the industrial zone, and the Aqaba-Amman highway (the most vital highway in Jordan) across Wadi Yutum, are all located in flood-prone areas on the apex, middle and toe parts of the alluvial fans of Wadi Yutum and other wadis to the north and south of Aqaba  . Such conditions are a particularly challenging task at present and will remain so in the future   . Unfortunately, all previous plans for the Aqaba area  -  suffer from a prominent lack of appreciation of the morpho-hydrological setting of the bajada at the lower W. Araba including the W. Yutum watershed, and other wadis north and south of W. Yutum. Limited or even poor perception of flash flood disasters among engineers, planners, managers, developers and the inhabitants as a whole is obvious in the Aqaba area. As Schick  stated over four decades ago, perhaps because of the “inexorable monotonous dryness”, the dry land urban population does not expect and is usually not prepared for flood hazards. Often, urban expansion and development is conducted in complete ignorance of the potential threat of flash flooding and the associated sediment discharge. Nevertheless, scattered hydrological studies/re- ports on flash floods have been carried out for W. Yutum   -   -  . These studies confirm that the Aqaba area has been exposed to recurrent flash floods since the 1950s, where low magnitude (5 - 7 year return periods), medium magnitude (20 - 25 year return periods), and high magnitude (50 year return periods) floods have been recorded. It is postulated that the northern parts of the Aqaba area are the terrain system most vulnerable to flash flooding hazard. Flash floods always have a significant negative impact on the city of Aqaba. It causes immense damage to lives and the infrastructure (i.e. residential areas, highway, roads, hydraulic structures) and coastal pollution. Therefore, since the city and the port of Aqaba are considered crucial to the Jordanian economy, the protection of the city/port from repetitive flooding is important to maintain future urban and economic development  . However, several physical factors contribute to the flooding problem, and influence flash flood severity in dry land urban centers. Among these factors are: morphological setting and site characteristics, topography and slope, drainage properties, available hydraulic structures, land cover, hydrology and climate (i.e. rainfall intensity and duration, water loss (evaporation and infiltration), drainage networks and orders, environmental processes and human interventions   . It is not practical to conduct a flood hazard evaluation for large areas/watersheds from time to time based on field monitoring and measurements. Alternatively, remote sensing and GIS techniques coupled with geomorphometric and hydrologic parameters provide a rapid, low cost and efficient approach to determine areas of potential flooding, levels of risk, basins and sub- basins with similar susceptibility to flash flooding.
Watershed characteristics have been established in studies worldwide using conventional geomorphometric analysis methods. Geometric properties of drainage networks in relation to typology of stream networks, and quantitative analysis of linear, areal, shape and relief aspects of drainage basins were intensively addressed earlier  -  . Moreover, morphometric properties of drainage networks have been utilized to predict hydromorphic processes such as peaks of floods, sediment yields, and erosion rates  . The relationship between watershed morphometry and flooding impact has been explored thoroughly by Patton  . At present, delimitation of drainage basins, and calculation of morphometric parameters can be achieved using remote sensing tools DEMs, GIS software, and mathematical formulas elaborated by Strahler  -   -  . Flooding susceptibility assessment for dryland watersheds was executed recently for various watersheds along the Red Sea Coast and the Gulf of Suez, Egypt, which are environmentally comparable to the Gulf of Aqaba, W. Yutum watershed and other wadis of lower wadi Araba    -    . Evaluation of periodic flash floods based on El-Shamy’s model was carried out on dryland watersheds in the southern Red Sea coast of Egypt  . Furthermore, analysis of drainage basin morphometry was performed to assess flood prone areas utilizing El-Shamy’s approach in western Saudi Arabia  . Similarly, Abdel-Lattif and Sherief  evaluated flash floods for two wadis (W. Sudr and W. Wardan, Gulf of Suez, Egypt) based on morphometric analysis and El-Shamy’s approach. It is obvious that previous investigations employed either El-Shamy’s model to assess flash floods risk or, combined this with elaborated morphometric analysis methods developed and employed successfully for flash flood risk estimation. By contrast, Youssef et al.  developed a GIS morphometric analysis approach to estimate levels of flash flood risk for sub-watersheds within the main watershed. The results were compared with those achieved through the application of El-Shamy’s model. It can be concluded that both methods are rational for flash flood risk assessment. El-Maghraby et al.  employed morphometric analysis to evaluate the flash flood hazard of W. Al Hamd sub-basins, western Saudi Arbia. Nine morphometric parameters which have a direct influence on flooding were selected and analyzed in relation to flash floods. These parameters are: watershed area (A), drainage density (Dd), stream frequency (fs), shape index (Si%), slope index (Sw), relief ratio (Rf), ruggedness number (Rn), texture ratio (Rt), and weighted mean bifurcation ratio (WMRb). A hazard scale number starting with 1 (lowest) to 5 (highest) has been assigned to all parameters. Calculation and spatial distribution of hazard degrees have been discussed in detail based on Davis  .
2. Description of the Study Area
2.1. Geology and Geomorphology
Wadi Yutum and the Aqaba area are located in southern Jordan, east of the Gulf of Aqaba and lower Wadi Araba. The study area lies between 29.15 N to 30.01 N latitude and 34.58 E to 35.52 E longitude (Figure 1(a) and Figure 1(b)). The Aqaba area comprises the northern edge of the Arabian plate east of the lower W. Araba fault and the Gulf of Aqaba  . The exposed rocks sequence consists mainly of igneous rocks (pre Cambrian)which represent the granite horst overlooking lower W. Araba and the Gulf of Aqaba (Figure 2)   . The granite horst/inselberg landscape west of Quweria-Ram area is composed of aplite-granite, granodiorite, and quartzdiorite with an age ranging between 585 + 8 Ma  . Granodiorite and aplite granite are exposed at the base of sandstone inselberg (the erosional residuals) west of Jebel Rum to form a minor horizontal bench, or plinth at the igneous/sandstone contact  . Slopes underlain by these rocks range from 30˚ to 40˚, while slopes overlying sandstones are almost near-vertical. To the west of Quweira, remnants of large strips of the pre-sandstone peneplain were preserved. The granite is highly weathered, and extensively intruded by dykes. Due to variations in the composition of granite   , dyke type and density  , and the rate of weathering  topographic expression is varied, ranging from rugged mountainous terrain to undulating ridges and subrounded morphology  .
The important landforms of the central part of W. Yutum watershed are the “cliffed sandstone inselbergs”, which rise 800 mm above the playas and sand flats of the desert floor. To the north-east, the watershed is bounded by Ras En Naqb escarpment (200 - 350 m high), which is formed of Cretaceous limestone, shales, marls and Kurnub sandstones (Lower Cretaceous). It is believed that these strata originally extended over what is now the inselberg landscape although the original southern terminus is unknown. However, on the west side of the Dead Sea Rift, remnants of these strata are found clear down to the southern extreme of the Sinai Peninsula  . Landslides are active on Ras En Naqb escarpment following recurrent intense and heavy rainstorms,
Figure 1. (a) The study area; (b) Google Earth Pro Image (2013) for W. Yutum watershed.
Figure 2. Geology of W. Yutum watershed.
and occasionally have threatened the Amman-Aqaba vital highway. Oversteepening of slopes due to highway/ road construction, fluvial erosion and differential weathering encourage rock falls, rock slides, and slumping from heavily jointed limestone beds. Percolation of water from relatively abundant precipitation over Ras En Naqb (140 mm annually), towards impervious marl, shale, and clay of the nodular limestone unit reduce the shear strength of the materials, thus activating landslides. At the base of the escarpment, pure white Ordovician sandstone with scattered quartz granules and pebbles is exposed (Figure 2). The sandstone strata dip gently to the east; consequently, younger sandstone units are exposed in succession on the top of inselbergs as one proceeds to the east. Although the total thickness of sandstone units in S-W Jordan is about 900 m, no complete section can be seen at one site. The sandstones are fairly weakly cemented, a property which aids in their weathering and erosion along joints and fissures  .
Three types of inselbergs at the central part of W. Yutum around Quweira exist (Figure 3). These can be differentiated based on slope morphology and slope categories  . The first type is of relatively gentle slopes (up to 40˚ - 45˚), with numerous steps a few meters high and benches covered by debris. The second type is steep-side inselbergs, with cliffs towering above the base by as much as 110 - 150 m, often with a very sharp piedmont angle, little talus, and rather flattish summit surfaces. The slopes are subjected to cavernous weathering and tafoni formation. The third type is a dome-like hill, with steeply sloping sides, and convex in profile. It is associated with iron-depleted beds of pure white sandstone, and it gives rise to locally more rounded hillslope morphology within the steep-sided inselbergs with sharp breaks, and steeped slope morphology.
The bajada of the lower W. Araba, is composed of several large coalescing alluvial fans deposited by W. Yutum and other wadis (to the north and south of it), at the foot of the granite mountains. The mean slope recorded ranges between 4˚ and 5˚. The fan surface area is found to be proportional to the drainage basin and because of the greater sediment inputs deposited at the mouth of the wadi. Alluvial fans of the bajada revealed periods of heavy rains or flash floods with a resultant debris flow and stream flood deposition  . The alluvial fans of the lower W. Araba were terminated with Qa’ed Dafiya (2 km width and 12 km length) at the western margin of the bajada.The southernmost sector of the W. Yutum huge fan directs its distributaries towards Aqaba City.
To the west of Ras En Naqb and the playa system, and east of Quweira, the W. Yutum major course (although almost always dry) cuts through the granite horst and drains into lower Wadi Araba just north of Aqaba. The wadi is separated from the playas by sand sheets. Apparently, with the dominance of the hyper-arid climate, there is rarely enough water collected in the playas to drain out into the wadi. Yet in a gross sense, Wadi Yutum is the ultimate drain of the entire inselberg region north of the Saudi border drainage divide, while the playas in this regard serve as intermediate base levels  .
Figure 3. Slope categories.
Structurally, W. Yutum watershed is characterized by normal faults and vertical joint sets. Lloyd  distinguished two major groups of faults: one trending roughly N-S (or more specifically 5˚ to 20˚E of N), the other trending NW-SE. The Quweira fault with a large vertical down throw of up to 600 m (Figure 2), and a fault running down of W. Marsad, with vertical offset of 400 m exemplify the former group of faults. These two adjacent faults form a graben, and similar features occur among smaller displacement faults further east. NW- trending faults give rise to having minor offsets, as stated by Heimbach and Meiser  based on geoelectrical resistivity survey. The most prominent NW-trending fault is the Hiswa fault with a downthrow of 300 m. It runs along a wadi/playa system southeast of Qa’ Disa. Several graben-forming normal faults in the Quweira area are hidden beneath unconsolidated sand sheets   . Lloyd  also reported the presence of two major joint directions in the area, corresponding roughly with the two main fault directions. Numerous faults were also mapped by Bender  along part of the Ras En Naqb escarpment and belong to the NW-trending category. The alluvial fan of W. Yutum endangers directly a quarter of the unplanned dwelling zone. The steep overall slope of the fan (4˚ - 5˚) facilitates flood waters to reach any point along an 8 km circumference of the fan sector, and together with smaller fan adjoining to the north, creates flooding hazards for the Aqaba airport with its 3 km long runway. The high relative relief and bare granite facilitate the background for extremely destructive flooding in the downstream urbanized area developed between the granite horst front and the inland Qa’s and sabkhas with a large damage potential  .
2.2. Climatic Conditions and Flooding Problems
The climate of W. Yutum is hyper-arid in the Aqaba-Quweira, and Ras en Naqb areas, with a slight variation between these areas in terms of rainfall amounts and temperature. The climate is hot, dry and arid in the Aqaba area (lower W. Yutum) and Quweira (middle W. Yutum), and cold, dry and arid on the crest of Ras En Naqb (upper W. Yutum). Due to variation in topography, the average annual rainfall varies from 37 mm at Aqaba, to 65 mm at Quweira town (800 m a.s.l.), to 140 mm at Ras En Naqb village (1573 m a.s.l., Figure 4). It has been concluded by Lloyd  that tributaries of W. Yutum watershed have about five to seven runoff producing storms per year, which last no more than a few hours, and are capable of producing serious flash floods. Seventy percent of rainfall events occur between December and January. Furthermore, single storm rainfall totals can be very high and may far exceed the mean annual rainfall; thus, rainfall intensities are very high. Rainfall is often of the convective type with some effect of orography   . The convective rains usually occur during warm
Figure 4. Elevation (m).
months such as October, March and April, and are characterized by intense rainfall of short duration similar to that which occurred, for example, in 1991, 1993, 1994, 2006, 2010, 2012, 2013, 2014, and 2015 respectively ( mainly based on field observations, local news papers, and published work, i.e. reference  ). A localized pattern of rainfall is also characteristic of W. Yutum catchment, where rainfall might affect limited areas only. Sixty percent of the total rainfall comes from spotty rain, where spottiness of rainfall is clearly pronounced in the fall and spring thunderstorms  , which occur at an annual average of four events/year in southern Jordan  . Most of the annual rainfall may fall during one event and concentrated occasionally in but a few hours, or less. Accordingly, rainfall intensity is predicted to be high and it may reach 120 mm/h for 10 min  .
Since the 1950’s, the Aqaba weather station has recorded great events of maximum rainfall intensities over 24 hrs. Examples are: 34 mm (1953), 27 mm (1955), 36.3 mm (1963), 32.8 mm (1966), 65 mm (1975), 42.5 mm (1980). Such figures show that high rainfall intensity occasionally resulted in severe flash flooding associated with large sediment load. During the storm of 11 March 1966 for example, Aqaba recorded 37 mm, Ras En Naqb 71.2 mm, Ma’an 39 mm, and the highlands west of Ma’an received 60 mm fall in 4 hrs, with an intensity of 15 mm, Ras En Naqb 140 mm and Ma’an 40 mm. The flash flood of the 11 March 1966 rainstorm is categorized as a 50-year return period. The rainfall recorded during the thunderstorm of May 2014, was high in southern Jordan generally (in Tafila for example 80 mm, in Shawbak 45.6 mm, Wadi Musa and Petra 32.8 mm, and in Ma’an 24 mm), all of them are located north of the W. Yutum watershed. By contrast, Aqaba City received only 3.6 mm in 20 min, thus rainfall intensity was high and recorded 13 mm∙hr−1 in 10 min. It can be concluded that Aqaba City, the Amman-Aqaba highway, the towns and villages, and irrigated farms across W. Yutum catchment have been subjected to flash floods several times. The frequency of flash floods has increased recently (almost every year in the Aqaba area and lower W. Yutum) and this has caused serious damage to properties, infrastructure, loss of life, residential areas, roads, and brought hardship for the residents of Aqaba.
In the present investigation, two methods of flash flood risk assessment have been employed. The first method is elaborated by El-Shamy  for assessing the susceptibility of sub-basins to flash flooding. The second method was formulated to estimate the degree of flooding risk based on morphometric analysis  . Two flood hazard risk map layers were generated for the seventeen sub-basins of W. Yutum based on each method. A third map was then compiled by integrating the result of flood hazard risk mapping achieved from both methods in order to assess the correlation if any, between the produced maps, and to explore the common hazardous flooding categories that may be found between the sub-basins.
The objective of this research is to:
1) Conduct morphometric analysis for the seventeen sub-basins within Wadi Yutum,
2) Generate flood hazard maps based on El-Shamy’s approach and the morphometric ranking method, to demarcate the sub-basins with low, intermediate, and high flooding susceptibility,
3) Identify the sub-basins which are expected to endanger the Aqaba area, the Amman-Aqaba highway, the Bedouin-rural settlements, and the irrigated farm across W. Yutum watershed.
4) Provide significant information which can help in planning adjustment with flood hazards, and to construct efficient flood prevention structures  -  .
3. Materials and Methods
Topographic maps with a scale 1:50,000 (20 m contour interval) were provided by the Royal Jordanian National Geographic Centre (RJNGC). Adigital elevation model (DEM) is utilized to acquire and calculate morphometric parameters for the entire W. Yutum watershed and the seventeen sub-basins. The ASTER DEM (30 m spatial resolution) is provided on line for free to all users, and available in GeoTIFF format with geographical coordinates. It has a mesh size of 1 arc-second posting  . The elevation error of the DEM is found to be +7.4 m for forest land cover, and −0.7 m for bare land (such as Ras En Naqb, or the present study area), while the horizontal error is reported to be East/West shift of −0.13 arc-seconds and North-South shift of 0.19 arc-seconds when compared to the 10-m mesh DEM produced by the Geographical Survey Institute (GSI) of Japan  . It can be concluded that ASTER DEM is utilized worldwide in morphometric analysis research. The Arc GIS ve.10.1 software along with the available tools (mainly the Spatial Analyst extension) was employed to derive the drainage networks, mophometric parameters, and to generate different related maps (i.e. slope categories, elevation, etc.). The ArcHydro tool was also used to delineate the watersheds and stream networks of the study area. The extraction process passes through different steps discussed in detail in  . Geological maps of scale 1:250,000   , scale 1:100,000   , scale 1:500,000  and 1:50,000  , and LANDSAT ETM+ (30 m spatial resolution, dated 2010) were employed to:
1) distinguish various igneous rock, sandstones, carbonate rocks and Quaternary materials, faults and major joints affected the study area,
2) outline landform features of W. Yutum watershed and to verify drainage networks derived from ASTER DEM, and
3) understand the terrain units in the bajada of lower W. Araba, inselberg landscape and Ras en Naqb escarpment in the study watershed. Additional geomorphic and hydrographic information were achieved during field observation.
These data are geo-referenced to the UTM coordinate system (WGS 1984, Zone 36˚N) based on 1:50,000 topographic maps using Arc GIS. Areas affected by recent flash floods and expected to be exposed to future flooding were field inspected. Arc GIS 10.1 and ASTER DEM were employed to generate and digitize drainage networks for W. Yutum watershed and the seventeen sub-basins (Figure 5 and Figure 6) which impose serious hazards for the residents, infrastructure facilities, the Aqaba-Amman highway, roads, and economic activities in area. The stream order maps were compiled from flow direction map for each watershed using Stream Order tool. The stream ordering system used was based on Strahler’s method     . The entire W. Yutum watershed was found to be of seventh order.
4. Flash Floods Analysis Methods
Morphometric analysis for W. Yutum catchment and the seventeen sub-basing have been carried out within a GIS environment. The morphometric parameters which elucidate their impact on flash flooding behavior were calculated for the seventeen sub-basins using Arc GIS software, and the mathematical equations elaborated by Strahler  -  and illustrated in Table 1. The hydro-geomorphometric relations were employed to assess the degree of flooding risk based on the two methods noted above.
4.1. Morphometric Ranking Method
The first method was labeled the morphometric ranking method  . Two groups of morphometric parameters were used in the analysis and consist of 17 parameters. Three of these parameters are related to drainage network
Figure 5. W. Yutum sub-basins.
Table 1. Calculation of morphometric parameters.
Figure 6. Stream order of the 17 sub-basins.
aspect, six parameters depict basin geometry aspect, three represent drainage texture aspect, and five describe the relief characteristics aspect. Group I accommodates ten morphometric parameters which have a direct proportion to the degree of risk. Consequently, the higher the parameter value, the higher the risk degree. This group comprises the following parameters:
sub-basin area (A), circularity ratio (Rc), sub-basin relief (Bh), relief ratio (Rf), drainage density (Dd), stream frequency (Fs), stream order (U), total stream number (Nu), total stream length (Lu), and sub-basin slope (degrees) (Sw).
Group II comprises seven parameters which are inversely proportional to the risk. Subsequently, the higher the parameter value, the lower the degree of risk. This group consists of the following parameters: elongation ratio (Re), shape factor (Bs), compactness coefficient (Cc), ruggedness number (Rn), hypsometric integral (Hi), geometry number (Gn), and length of overland flow (Lo). The bifurcation ratio (Rb) has also been calculated in order to be employed in El-Shamy’s approach (or the second method), along with drainage density (Dd) and stream frequency parameters (Fs). Based on the computed values of the seventeen morphometric parameters, it is possible to calculate the flash flood risk according to the following steps:
With reference to the relationship between the values of the morphometric parameters and the risk of flash floods, the evaluation of each parameter has been carried out using a simple statistical approach elaborated and employed by Pradhan  and Youssef et al.  . Each morphometric parameter is classified into five rank categories. Rank 1 refers to low risk, and rank 5 denotes high risk. The equation applied is (Max-Min)/5. Then the values for each parameter are categorized into five intervals. For group I parameters, the higher the parameter value, the higher the degree of risk, thus the higher the rank. By contrast, for the group II of parameters, the higher the parameter value, the lower the degree of risk, hence, the lower the rank. Then the morphometric parameter values for each sub-basin were summed and normalized in order to obtain a morphometric number for each sub-basin. Based upon these summated numbers, the sub-basins were classified in terms of their risk of occurrence of flash floods.
4.2. El-Shamy’s Approach
The second approach for assessing susceptibility of sub-basins to flash flooding risk is that developed by El-Shamy  . He designates a simple morphometric method to estimate the flash flood risk levels and the degree of hazardousness for each sub-basin. Two different approaches were elaborated to determine the hazardous sub-watersheds. The first is based on the relationship between bifurcation ration (Rb) and drainage density (Dd), whereas the second approach utilized the relationship between bifurcation ratio (Rb) and stream frequency (Fs). Drainage density (Dd) refers to topographic dissection, runoff potential, infiltration capacity of surface materials, climate, and land cover of the watershed. In this regard, low values of Dd indicate optimal conditions for infiltration, thus decreasing runoff potential, while, high stream frequency (Fs) represents impermeable sub-surface materials, poor vegetation cover, high relief, and low infiltration capacity, hence increasing runoff potential   . Applying this relationship separately to each sub-basin will provide reasonable information on estimation of flooding risk and recharge potential. In order to assess such relations, three morphometric parameters (Rb, Dd, and Fs) must be calculated for each sub-basin  .
The resultant illustrations for Dd Rb and Fs Rb have to be plotted graphically, where each illustration contains two curves, dividing the area into three fields. These fields are described as follows:
- Field A refers to low possibility for floods, and high recharge potential.
- Field B represents sub-basins with intermediate possibility for floods and moderate potential for groundwater recharge.
- Field C indicates high possibility for floods and low recharge property. If a sub-watershed has two different fields, then the appropriate classification plot has been selected.
However, an integration of the results was achieved based on the morphometric ranking method with El- Shamy’s model. Such a procedure was executed through superimposition of flood risk map produced according to the morphometric ranking method, with the flood risk map generated through El-Shamy’s approach. This task would make it possible to identify: (i) the common sub-watersheds falling under each category of flooding risk, (ii) determine the most risky sub-watersheds in terms of flash flooding, and (iii) recognize the most vulnerable sites/sub-basins against flooding, and the section of Amman-Aqaba highway (across W. Yutum)which may be affected by flood damage in the future.
5. Results and Discussion
A methodology for flash floods assessment has been employed to understand the nature and extent of flooding risk within W. Yutum watershed. The potential hazard was analyzed, and areas vulnerable to flooding, and expected to cause heavy damage to the residents and their livelihood were delineated. Since the 1950s, W. Yutum and Aqaba City were exposed to flash floods of different magnitudes, as a result of climatic conditions, human intervention and engineering activities, and recurrent heavy rainfall storms. In order to assess the hazardous effect of flash floods on seventeen sub-basins related to W. Yutum, hydro-morphometric parameters were utilized in the analysis. The degree of flooding risk characterizing each sub-basin was established using two different methods: (a) El-Shamy’s approach for flooding hazard assessment based on three morphometric parameters;and (b) the morphometric ranking method elaborated in     . The second method is considered a semi-quantitative measure determined on seventeen morphometric parameters.
5.1. Flash Floods Risk Assessment: El-Shamy’s Approach
In line with this approach, the relationship between bifurcation ratio (Rb) and drainage density (Dd), and then the relationship between Rb and stream frequency (Fs) have been utilized. The bifurcation ratios for the seventeen sub-basins range between 3.13 and 7.16, whereas, the mean bifurcation ratio for the entire W. Yutum is 4.6, which suggests a remarkable influence of structural control on sub-basins nos. 4, 5, 7, 9, 10, 12, 13, 14, 15, and 17. Little influence was observed on sub-basins nos. 3, 11 and 16 (Table 2). Furthermor, low values were achieved for the inselbergs―sand sheet areas. In general, these range between 1.29 and 2.66. Such figures indicate high infiltration, thereby decreasing runoff potential. By contrast, relatively high stream frequency values characterized the granite watersheds, which denotes impermeable sub-surface materials, sparse vegetation, high-rugged terrain, and low infiltration capacity.
Sub-basins vulnerable and exposed to flash flooding were demarcated and assessed in order to determine sub-basins of low, intermediate and high susceptibility to flooding based on the relationship of the two morphometric parameters (Rb vs. Dd and Rb vs. Fs); the final flood hazard map was then generated with the help of Arc GIS tool. Based on the relationship between Rb and Dd, sub-basins nos. 1 (the eastern part of Ras En Naqb
Table 2. Morphometric characteristics of W. Yutum sub-basins.
(A) area, (P) perimeter, (Lb) length, (Bh) relief, (Rr) relief ratio, (D) slope (degrees), (Rn) ruggedness number, (Hi) hypsometric integral, (Gn )geometry number, (u) stream order, (Nu) total stream number, (Lu) total stream length, (Rb) bifurcation ratio, (Dd) drainage density, (Fs) stream frequency, (Lo) length of overland flow, (Re) elongation ratio, (Rc) circularity ratio, (Sf) shape factor, (Cc) compactness coefficient.
escarpment), 3 (Q’a Disa-El Ghal), 8 (Yutum El-Umran), 11 and 16 (Figure 7 and Figure 8) are characterized as having high susceptibility to flooding. Sub-basins nos. 2, 4 (middle part of Ras En Naqb escarpment), 6, 10 (W. Marsad and W. Rumman), 12, 13, and 14 (northwestern part of the granite horst) represent the category of intermediate flooding susceptibility. By contrast, sub-basins nos. 5, 7, 9 (inselberg landscape), 15 and 17 are of low susceptibility to flooding. Similarly, and based on the relationship between Rb and Fs (Figure 9), sub-basins nos. 1, 3, 8, 11 and 16 are considered of high susceptibility to flooding, while, sub-basins 2, 4, 6 and 14 are of intermediate susceptibility to flooding. Furthermore, sub-basins 5, 7, 9, 10, 12, 13, 15 and 17 are ranked as having low susceptibility to flooding.
Field observations reveal that the section of the old Amman-Aqaba highway was routed initially along W. Yutum bed at the lower catchment, and has been damaged several times due to flash floods. Following 1973, this section of the highway was re-located five meters above a fluvial terrace to protect the highway against flooding. However, a hydrological investigation  reported that a flood of 50-year return period (resembling the Ma’an flood of 1966) or greater is likely to damage this section of the highway, assuming that surface water will rise 0.3 - 0.5 m above the present road surface, with relative high velocities. The most hazardous sub-basin found is sub-watershed no. 8 (Yutum El-Umran), located completely on granite, and representing an area of 770.95 km2 (17.5% of the total area) of the W. Yutum watershed. At the lower W. Yutum, a water gap has been
Figure 7. (a) & (b) Flooding susceptibility: El-Shamy’s approach (Rb vs. Dd).
formed and has diverted the main W. Yutum watershed towards the Gulf of Aqaba. The water gap starts where the present main wadi course changes its southerly direction by 90˚ to follow a western course through a prominent gorge cut 400 m in a granite saddle. For the last few kilometers above the recent fan apex, the wadi floor narrows to less than 100 m between the bases of the granite rock walls on both sides. This water gap has been established by a westerly flowing stream which continued its headward erosion until it captured the interior basin just a few kilometers below the entrance of the main tributary of Yutum El-Umran. During the erosional phase following wadi diversion, the pirate stream extended its course rapidly in the easily eroded alluvium to form the main channel of the present Wadi Yutum across the axis of the watershed. Successive phases of incision may be attributed either to rift tectonics and subsidence of the Rift valley floor, or may be correlated with sea level fluctuations as the new wadi reached the Gulf   . Moreover, it is also obvious, that the most ha-
Figure 8. (a) & (b) Flooding susceptibility: El-Shamy’s approach (Rb vs. Fs).
zardous basin located on sedimentary carbonate rocks is sub-basin no. 1 (area: 488.69 km2, 11% of the total area). By contrast, sub-basins nos. 4 and 13 are of intermediate susceptibility to flash flooding. However, the relative abundance of rainfall on the Ras En Naqb area and the presence of heavily jointed, fissured weak lithology, facilitate landsliding activity along the section of the highway extending between Ras En Naqb town and north of Quweira. Thereby, repetitive landslides have disrupted the highway occasionally.
It can be concluded that similar and consistent results of flood hazard susceptibility can be achieved through El-Shamy’s approach for sub-basins categorized based on the relationship between Rb vs. Dd, and the relationship between Rb vs. Fs (Figure 7 and Figure 8). For example, sub-basins nos. 1, 3, 8, 11, and 16 are categorized in both morphometric relationships under the class of high flooding susceptibility. Moreover, sub-basins nos. 2, 4, 6, and 14 are ranked under the intermediate flooding susceptibility and sub-basins 5, 7, 9, 15, and 17 are classified under low flooding susceptibility.
Figure 9. Flooding susceptibility: morphometric ranking method.
5.2. Flash Floods Risk Assessment: Morphometric Ranking Method
The morphometric ranking method for determining the flash flood hazard and the degree of risk for the seventeen sub-basins, was applied to execute the required morphometric analysis. Table 3 illustrates the results of this method as expressed by ranking score for different morphometric parameters along with the total score values for each sub-basin. To generate the flood risk map using GIS, the total score values for the seventeen sub-basins were grouped into three categories of flooding susceptibility.
1) Low flooding susceptibility (30 - 40)
2) Intermediate flooding susceptibility (41 - 49)
3) High flooding susceptibility (50 - 60)
With respect to the seventeen sub-basins of W. Yutum watershed, sub-basins nos. 3, 4, 10, 12, 14, and 15 (35.3% of the total) have the lowest overall values, therefore representing the low flooding susceptibility rank (Table 3 and Figure 9). Whereas, sub-basins nos. 1, 5, 7, 9, 11, 13, 16 and 17 (47% of the total) have intermediate score values, thus characterized by moderate flooding susceptibility. By contrast, sub-basins nos. 2, 6, and 8 (17.7% of the total) have the highest overall score values, and thus represent the most dangerous sub-basins with high flooding susceptibility. In this regard two-thirds of W. Yutum sub-basins are expected to suffer moderate and high susceptibility to flooding. The results also indicated that the main hazardous and threatening sub-basins are located on the basement complex (granite), i.e. the lower part of W. Yutum and its major tributary Yutum El-Umran, and sub-basins descending from Ras en Naqb escarpment. Here, flash floods are mostly contributed by these sub-basins. Other sub-basins which terminated to the playas of the inselberg landscape are often of low susceptibility to flooding. Consequently, the Amman-Aqaba highway and the alluvial fan of W. Yutum which at present accommodate a considerable residential area, the Aqaba International Estate, the King Hussein International Airport, and the Northern Light Industrial and Logistic Area are seriously threatened by flooding. Thus, the protection of the city from repetitive flooding is essential to maintain future urban and economic development.
An integration of the results achieved based on the morphometric ranking method and El-Shamy’s approach was conducted through superimposition of the layers/maps produced. Such a procedure makes it possible to recognize the common sub-basins falling under each category of flooding susceptibility. The correlation shows that sub-basin no. 8 is the common sub-basin and ranked under high flooding susceptibility (Table 4, Figure 10 and Figure 11). Similarly, sub-basin no.13 is a common sub-basin and classified under moderate flooding sus-
Figure 10. Flooding susceptibility: Shamy’s approach (Rb vs. Dd) & morphometric ranking method.
Figure 11. Flooding susceptibility: Shamy’s approach (Rb vs. Fs) & morphometric ranking method.
Table 3. Ranks score for different morphometric parameters and the total score value.
Table 4. Classification of sub-basins flooding susceptibility based on the two methods employed in flash floods assessment.
*Correlation between sub-basins under the two methods of floods risk assessment. +Correlation between sub-basins under the two morphometric relationships proposed by El-Shamy for flood risk assessment.
ceptibility based on both the morphometric ranking method and El-Shamy’s method (Rb vs. Dd). Moreover, three sub-basins, nos. 10, 12 and 15, are also common, and ranked under the category of low flooding susceptibility. Sub-basin no. 6 is categorized under high flooding susceptibility based on the morphometric ranking method, and is classified under intermediate flooding susceptibility based on El-Shamy’s approach. Equally, two sub-basins nos.9 and 17 are ranked under intermediate flooding susceptibility based on the morphometric ranking method, and of low flooding susceptibility based on El-Shamy’s approach. In parallel, sub-basins nos. 10 and 12 are ranked under low susceptibility to flooding based on the morphometric ranking method, and El- Shamy’s approach (Rb vs. Fs), and of intermediate flooding category based on El-Shamy’s approach (Rb vs. Dd). It can be demonstrated that 47% (8 sub-basins) of the total number of sub-basins can be categorized under high flooding susceptibility based on the morphometric ranking method, and El-Shamy’s approach. Thus they are considered the highly hazardous sub-basins within W. Yutum watershed. GIS tool and field observations enable assessment of flooding risk, and demarcation of flood-prone sub-basins as displayed on the drainage basin hazard map.
In the present investigation, GIS based hydro-morphometric analysis, field and satellite data were employed to estimate flash flood risk, and to compile a spatial potential flood hazard map which illustrated flash flood prone―areas across W. Yutum watershed. The generated maps/information were meant to help planners and decision makers understand the spatial distribution of flood hazard conditions, in order to formulate flood protection measures to minimize the negative impacts of repetitive flooding on the residents and infrastructure of Aqaba city. Floods risk analysis was performed using two morphometric analysis methods. With the aid of GIS, it was possible to delineate the sub-basins of W. Yutum which might be affected by floods of low, moderated and high susceptibility. Results of flash flood risk estimation based on the relationship between Rb and Dd (El-Shamy’s approach) showed that sub-basins nos. 1, 3, 8, 11 and 16 were characterized as being of high susceptibility to flooding, whereas, sub-basins nos. 2, 4, 6, 10, 12, 13 and 14 represented the category of moderate flooding susceptibility. By contrast, sub-basins nos. 5, 7, 9, 15 and 17 are classified under the category of low flooding susceptibility. Similarly, and based on the relationship between Rb and Fs, sub-basins nos. 1, 3, 8, 11 and 16 are categorized as of high susceptibility to flooding, while, sub-basins nos. 2, 4, 6 and 14 are of intermediate susceptibility to flooding. Furthermore, sub-basins 5, 7, 9, 10, 12, 13, 15 and 17 are considered to be of low susceptibility to flooding.
Based on the morphometric ranking method, sub-basins nos. 3, 4, 10, 12, 14, and 15 (35% of the total) have the lowest overall score values, and thus represent the low flooding susceptibility rank. Sub-basins nos. 1, 5, 7, 9, 11, 13, 16, and 17 (47% of the total) have intermediate score values, and thus are characterized by moderate flooding susceptibility. By contrast, sub-basins nos. 2, 6, and 8 (17.7% of the total) have the highest overall score values and consequently, represent the most dangerous sub-basins with high flooding susceptibility. In this regard, two-thirds of W. Yutum sub-basins are expected to suffer moderate and high flooding susceptibility.
However, the superimposition of the thematic maps produced based on the morphometric ranking method, and El-Shamy’s approach, showed that sub-basin no. 8 is the common sub-basin and ranked under high flooding susceptibility. Similarly, sub-basin no. 13 is a common sub-basin and classified under moderate flooding susceptibility based on both the morphometric ranking method and El-Shamy’s approach (Rb vs. Dd). Moreover, three sub-basins nos.10, 12, and 15 are also common, and ranked under the category of low flooding susceptibility. Sub-basin no.6 is categorized under high flooding liability based on the morphometric ranking method, and is classified under moderate flooding potentiality based on El-Shamy’s approach. In parallel, two sub-basins nos. 9 and 17 are ranked under moderate flooding susceptibility based on morphometric ranking method, and of low flooding susceptibility based on El-Shamy’s approach. Equally, sub-basins nos. 10 and 12 are ranked under low flooding susceptibility based on the morphometric ranking method and El-Shamy’s method (Rb vs. Fs), and of moderate flooding probability based on El-Shamy’s approach (Rb vs. Dd). It can be concluded that 47% (8 sub-basins) of the total can be categorized under high flooding potentiality based on the morphometric ranking method, and El-Shamy’s method. Consequently, they are considered the highly hazardous sub-basins within W. Yutum watershed. It is obvious that the most hazardous sub-basin found is Yutum El-Umran (no. 8)which is developed in the basement rocks (granite)with an area of 770.95 km2 (17.5% of the total area). Furthermore, sub-basin no. 1, which is located on sedimentary carbonate rocks, and descending from the Ras En Naqb escarpment, and lower W. Yutum are considered hazardous as well. Flash floods here are mostly attributed to these sub-basins. As a result, the Amman-Aqaba highway and the alluvial fan of W. Yutum which accommodates substantial residential areas, and major infrastructure facilities are seriously exposed to flash floods. Therefore, the protection of the city from disastrous flooding is essential to maintain future sustainable development.
Repetitive flash floods can be of great advantage if appropriate dams constructed at the narrowest sections of W. Yutum and Yutum El-Umran to retain water to reduce flow velocity and to enhance recharge in the alluvial fan of W. Yutum. Dams can be constructed by utilizing available local weathering products in the form of loose materials and boulders and deposits of fluvial terrace of W. Yutum. An earlier field survey of sediment sources across areas surrounding the back road of Aqaba revealed the occurrence of large boulder fields to the west of section A of the road, a few kilometers south of W. Yutum bridge. Dimensions of each dam depend on local geological and geomorphic conditions of the selected site. High potential flash floods in W. Yutum and Yutum El-Umran are expected to hold a high potential recharge of ground water towards the huge alluvial fan of W. Yutum. The residential areas and other urban services are encroaching at present across the fan surface. Although the toe part of the fan close to Qa el Dafiya is often inundated following flash floods, palm-tree plantations could be planted, and a deep well can be installed at the apex of the fan to provide water for the light industrial area north of Aqaba City. It is also possible to carry out an afforestation project in the Ras En Naqb escarpment based on a water harvesting system suitable for the dryland areas. Furthermore, a warning system is recommended for the Aqaba city to alarm officials and residents before the occurrence of a flash flood disaster.
 Farhan, Y. (1989) Urban Planning and the Evaluation of Flood Hazard in Arid Lands of Jordan: Aqaba Case Study. In: Farhan, Y., Beheiry, S. and Abu Satat, M., Eds., Geomorphological Studies on Southern Jordan, Publication of the University of Jordan, 181-206. (In Arabic)
 Schick, A., Grodek, T. and Wolman, M.G. (1999) Hydrologic Process and Geomorphic Constraints on Urbanization of Alluvial Fan Slopes. Geomorphology, 31, 325-335.
 Cooke, R.U., Goudie, A.S. and Doornkamp, J.C. (1978) Middle East-Review and Bibliography of Geomorphological Contributions. Quarterly Journal of Engineering Geology and Hydrogeology, 11, 9-18.
 Doornkamp, J.C., Brunsden, D., Jones, D.K.C., Cooke, R.U. and Bush, P. (1979) Rapid Geomorphological Assessment for Engineering. Quarterly Journal of Engineering Geology and Hydrogeology, 12, 189-209.
 Potter, K.W. and Frevert, D.K., Eds. (2010) Proceedings of Watershed Management 2010: Innovations in Watershed Management under Land Use and Climate Change. Madison, 23-27 August 2010, American Society of Civil Engineers.
 Peng, S.H., Shieh, M. and Fan, S. (2010) Potential Hazard Map for Disaster Prevention Using GIS Base Linear Combination Approach and Analytic Hierarchy Method. Journal of Geographic Information System, 4, 403-411.
 Grodek, T., Lekack, J., and Schick, A. (2000) Urbanizing Aluvial Fans as Flood-Conveying and Flood-Reducing Systems: Lessons from the October 1997 Eilat Flood. In: Hassan, M., Ed., The Hydrology-Geomorphology Interface: Rainfall, Floods, Sedimentation, Land Use, International Association of Hydrological Sciences Publication, No. 261, 229-250.
 Youssef, A.M., Pradhan, B., and Hassan, A.M. (2011) Flash Flood Risk Estimation Along the St. Katherine Road, Southern Sinai, Egypt Using GIS Based Morphometry and Satellite Imagery. Environmental Earth Sciences, 62, 611-623.
 Arnous, O., Aboulela, A. and Green, R. (2011) Geo-Environmental Hazards Assessment of the North Western Gulf of Suez, Egypt. Journal of Coastal Conservation, 15, 37-50.
 Abdel-Lattif, A. and Sherief, Y. (2012) Morphometric Analysis and Flash Floods of Wadi Sudr and Wadi Wardan, Gulf of Suez, Egypt: Using Digital Elevation Model. Arabian Journal of Geosciences, 5, 181-195.
 Abdalla, F., El-Shamy, I., Bamousa, A., Mansour, A., Mohamed, A. and Tahoon, M. (2014) Flash Floods and Ground- water Recharge Potentials in Arid Land Alluvial Basins, Southern Red Sea Coast, Egypt. International Journal of Geosciences, 5, 971-982.
 Abduladheem, A., Elmewafey, M., Beshr, A. and Elnaggar, A.A. (2015) Using GIS Based Morphometry Estimation of Flood Hazard Impacts on Desert Roads in South Sinai, Egypt. International Journal of Scientific and Engineering Research, 6, 1593-1599.
 Redwan, M., Abdel Moneim, A., Abo Amra, M. and Youssef, A. (2015) Watershed Characteristics and Flood Hazard, Assessment of the Area West of Sohag, Upper Egypt. International Journal of Environmental Engineering Research, 4, 18-26.
 Elmoustafa, A. and Mohamed, M. (2013) Flash Flood Risk Assessment Using Morphological Parameters in Sinai Peninsula. Open Journal of Modern Hydrology, 3, 122-129.
 Murphy, B. (2010) Hydrologic Analysis for a Hyper-Arid Region in the Middle East. In: Potter, J. and Frevert, D., Eds., Proceedings of Watershed Management 2010: Innovations in Watershed Management under Land Use and Climate Change, 23-27 August 2010, Madison, 1238-1247.
 Farhan, Y. and Anbar, A. (2014) Fragile Landscape: Impact and Consequences of May 2014 Flash-Flood Disaster in the Aqaba Area, Southern Jordan. Research Journal of Environmental and Earth Sciences, 6, 451-465.
 Hamdi, M.R., Abu Allaban, M., Al Shayeb, A., Jaber, M. and Momani, N.M. (2009) Climatic Change in Jordan. American Journal of Environmental Sciences, 5, 58-68.
 Al-Qinna, M., I., Hammouri, N.A., Obeidat, M.M, and Ahmad, F.Y. (2011) Drought Analysis in Jordan under Current and Future Climates. Climatic Change, 106, 421-440.
 Schick, A. (1971) A Desert Flood: Physical Characteristics, Effect on Man, Geomorphic Significance, Human Adaptation: A Case Study in the Southern Arava Watershed. Jerusalem Studies in Geography, 2, 91-155.
 Horton, R. (1945) Erosional Development of Streams and Their Drainage Basins, Hydrological Approach to Quantitative Morphology. Geological Society of America Bulletin, 56, 275-370.
 Miller, V. (1953) A Quantitative Geomorphic Study of Drainage Basin Characteristics in the Clinch Mountain Area, Virginia and Tennessee. Project NR 389-402, Technical Report 3, Department of Geology, ONR, Columbia University, New York.
 Schumm, S. (1956) Evolution of Drainage Systems and Slopes in Badlands at Pert Amboy, New Jersey. Geological Society of America Bulletin, 67, 597-646.
 Morisawa, M.E. (1959) Relation of Quantitative Geomorphology to Stream Flow in Representative Watersheds of the Appalachian Plateau Province. Project NR 389-042, Technical Report 20, Department of Geology, ONR, Columbia University, Geography Branch, New York.
 Kamp, U., Bolch, T. and Olsenholler, J. (2005) Geomorphometry of Cerro Sillajhauy (Andes, Chile/Bolivia): Comparison of Digital Elevation Models (DEMs) from ASTER Remote Sensing Data and Contour Maps. Geocarto International, 20, 23-33.
 Ozdemir, H. and Bird, D. (2009) Evaluation of Morphometric Parameters of Drainage Networks Derived from Topographic Maps and DEM in Point of Floods. Environmental Geology, 56, 1405-1415.
 Prasannakumar, V., Vijith, H. and Geetha, N. (2013) Terrain Evaluation through the Assessment of Geomorphometric Parameters Using DEM and GIS: Case Study of Two Major Sub-Watersheds in Attapady, South India. Arabian Journal of Geosciences, 6, 1141-1151.
 Singh, P., Thakur, J. and Singh, U. (2013) Morphometric Analysis of Morar River Basin, Madhya Pradesh, India, Using Remote Sensing and GIS Techniques. Environmental Earth Sciences, 68, 1967-1977.
 Markose, V., Dinesh, A. and Jayappa, K. (2014) Quantitative Analysis of Morphometric Parameters of Kali River Basin, Southern India, Using Bearin Azimuth and Drainage (bAd) Calculator and GIS. Environmental Earth Sciences, 72, 2887-2903.
 Youssef A.M., Pradhan, B., Gaberr, A.F.D. and Buchroithner, M.F. (2009) Geomorphological Hazard Analysis along the Egyptian Red Sea Coast between Safaga and Quseir. Natural Hazard and Earth System Science, 9, 751-766.
 Hassan, O. (2000) Salient Geoenvironmental Parameters of Ras Malaab-Abu Zenima Area, Gulf of Suez, Egypt, with an Emphasis on Flash Flood Potential and Mitigative Measures. Egyptian Journal of Remote Sensing and Space Science, 3, 37-58.
 El Maghraby, M., Masoud, M. and Niyazi, B. (2014) Assessment of Surface Runoff in Arid, Data Scarce Regions; An Approach Applied to Wadi Al-Hamd, Al Madinah Al Munawarah, Saudi Arabia. Life Science Journal, 11, 271-289.
 Natural Resources Authority (1976) Medium Grained Alkali Granites of the Aqaba-Quweira Area for Use as a Source of Feldspar of Ceramics. Unpublished Report Geological Survey and Bureau of Mines, Natural Resources Authority, Amman.
 Van den Boom, G. and Lahloub, M. (1964) Geological and Petrological Investigations of the Igneous Rocks in the Area of Quweira, Southern Jordan. Unpublished Report, Natural Resources Authority, Amman.
 Abed, A.M. (1985) On the Supposed Precambrian Paleostructure along the Dead Sea Rift, Jordan. Journal of the Geological Society of London, 142, 527-531.
 Osborn, G. and Duford, J.M. (1981) Geomorphological Processes in the Inselberg Region of South-Western Jordan. Palestine Exploration Quarterly, 112, 1-17.
 Borga, M. Boscolo, P., Zanon, F. and Sangati, M. (2007) Hydrometeorolgical Analysis of the 29 August 2003 Flash Flood in the Eastern Italian Alps. Journal of Hydrometeorology, 8, 1049-1067.
 Rusjan, S., Kobold, M. and Mikos, M. (2009) Characteristics of the Extreme Rainfall Event and Consequent Flash Floods in W Slovenia in September 2007. Natural Hazards and Earth Systems, 9, 947-956.
 Roca, M. and Davison, M. (2010) Two Dimensional Model Analysis of Flash-Flood Processes: Application to the Boscastle Event. Journal of Flood Risk Management, 3, 63-71.
 Arghius, V., Ozunu, A., Samara, I. and Rosian, G. (2014) Results of the Post-Flood Disaster Investigations in the Transylvanian Depression (Romania) during the Last Decade (2001-2010). Natural Hazards and Earth System Science, 14, 535-544.
 ASTER GDEM Validation Ream. “ASTER Global Digital Elevation Model Version 2 Summary of Validation Results”.
https://www.jspacesystems.or.jp/ersdac/GDEM/ver2Validation/Summary GDEM2 validation report_ final.pdf
 Strahler, A. (1958) Dimensional Analysis Applied to Fluvially Eroded Landforms. Geological Society of America Bulletin, 69, 279-299.