Prioritization of Sub-Watersheds in a Large Semi-Arid Drainage Basin (Southern Jordan) Using Morphometric Analysis, GIS, and Multivariate Statistics ()
1. Introduction
Rapid population growth(≈ 3% annually), and the persistent need for food production during the 1960s and ‘70s, caused continuous expansion of rainfed cultivation on marginal areas (average annual rainfall < 250 mm), forest and rangeland, which in turn increased the pressure on soil and water resources. The destruction of vegetation cover historically and land use abuse, coupled with the absence of efficient conservation measures, and compounded by improper farming practices accelerates soil erosion. The impact of land use/cover changes on soil erosion risk in northern Jordan was assessed [1] [2] . Cultivated land with poor conservation measures exhibit a higher rate of soil erosion and decline in soil fertility. The only conservation practice dominant is inefficient old stone terracing where rainfed “mixed farming” is practiced [3] . Further, soil erosion is becoming more serious on moderate and steep slopes which were transformed into cultivated land. Repetitive heavy rainstorms are the major significant triggering factor for extreme soil erosion, landslide activity, and floods. Severe storms with maximum daily intensity in the range of 2.1 - 6.66 mm∙h−1 are common in the highlands region [4] [5] [6] . Several case studies were carried out on watersheds prioritization in Jordan using morphometric analysis method, multivariate statistics, soil erosion susceptibility, and RUSLE soil loss modeling. For example, fourteen mini-watersheds of W. Shueib (Central Jordan) were prioritized for soil and water conservation [7] based on the method of morphometric analysis [8] [9] [10] , and soil erosion susceptibility analysis [11] using GIS. Degradation of vegetation covers including forest, and the existence of W. Shueib structure along the eastern part of the catchment, caused serious soil erosion and landslides activity. Therefore, flooding and sediment discharge into the W. Shueib reservoir have increased during heavy rainstorms. Eleven mini-watersheds (78.6% of the total) are classified in the categories of moderate, high, and very high priority. Thus, they should be considered as potential areas for preferential conservation intervention, and must be prioritized immediately for soil and water conservation practices. Moreover, prioritization was performed for thirty one third-order mini-watersheds connected to W. Kerak (Southern Jordan) using morphometric analysis and RUSLE soil loss modeling methods, and GIS techniques [12] . Based on morphometric parameters and soil loss values, and their ranking scores, five priority classes were recognized: extremely high, very high, high, moderate and low priority. Fifteen mini-watersheds (48.4% of the total) are classified as extremely high, very high, and high priority. Therefore, these mini-watersheds must be given the highest priority for soil and water conservation to ensure future sustainable agriculture. Remote sensing and GIS techniques were adopted for prioritization of the Zarqa River catchment. Forty three fourth-order sub-basins were prioritized based on morphometric and Principal Component Analysis (PCA) [13] . The optimal priority classes were recognized and the effectiveness of erosion risk parameters in watershed prioritization was assessed. Prioritization of the 43 sub-basins was conducted using ten linear and shape parameters, and three parameters which are found to be of high correlation with components 1 and 2. Then different sub-basins were prioritized by ascribing ranks based on the calculated compound parameter (Cp) using the two approaches. Comparison of the results indicated that prioritization of sub-basins based on morphometric analysis is more consistent and serves for better decision making in conservation planning as compared with the PCA approach. Seventeen sub-basins (39.5% of the total) grouped under the high priority class, and are subjected to high erosion risk. Nineteen sub-basins (44.2% of the total) are classified as moderate priority for soil conservation. The transformation of enormous areas of forest and rangeland into farming practices and the expansion of rainfed “mixed farming” (mainly cereals) accelerates soil erosion. Soil erosion eventually induced low agricultural productivity, changes in river morphology, high sediment yield and reduction in reservoir storage capacity. The future life of constructed and proposed dams in Jordan, i.e., W. Alarab [14] ; King Talal Reservoir [15] ; W .Mujib dam [16] , and W. Wala dam [17] are seriously threatened by the inflow of high sediment loads. Soil erosion loss for the catchments draining to the rift was found to be 1.328 million tons year−1, which means, 0.14 cm of the top soil is eroded annually [18] . It is obvious that the recorded soil erosion rates in different parts of the highlands in Jordan are greater than the accepted soil loss tolerance limits for the Mediterranean environment which were estimated at from 2 to 12 ton∙ha−1∙year−1 [19] [20] .
Prioritization of sub-watersheds refers to the “ranking of different sub-watersheds as the order they have to be selected for appropriate soil conservation measures adaptation” [21] . Erosion risk parameters must be calculated to prioritize sub-watersheds for soil conservation. The linear parameters possess a direct relationship with erodibility. Therefore, the highest value of the linear variables was ranked 1, the second highest value ranked 2 and so on, and the least value was rated last in rank. Furthermore, the shape parameters possess an inverse relation with erodibility, thus, the lower their values, the greater the erodibility [22] . Consequently, the lowest value of shape parameter was rated as rank 1 and the second lowest as rank 2 and so on, and the highest value was rated last in rank [23] . The Compound Parameter (Cp) was computed by adding up all the ranks of linear parameters, as well as shape parameters, and then, dividing by the number of all parameters. Following the rating of every single morphometric parameter, the ranking values for all linear and shape parameters referring to each sub-watershed are added up for each of the sub-watersheds to achieve a compound parameter (Cp) score based on the average value of these parameters. Further, the sub-watersheds having the lowest compound parameter score was designated the highest priority, the next higher value was assigned as second priority and so on [21] . Highest priority indicates the greater degree of soil erosion in that particular sub-basin, thus, it is considered a potential area for applying soil conservation measures [24] .
Several methods have been developed and elaborated for watersheds prioritization. Eleven of the twenty five studies (44% of the total) employed morphometric analysis method in prioritization [8] [10] [13] [21] [25] [26] [27] . Whereas the other 14 studies (56% of the total) adopted the morphometric analysis method combined with one or two of the following methods, such as: sediment yield index method(SYI) [8] [22] ; sediment production rate (SPR) method [28] ; USLE/RUSLE models for estimating soil loss [7] [29] [30] ; soil erosion susceptibility analysis [12] ; runoff potential method [31] ; land use/cover analysis [32] [33] ; sediment yield index(SYI)and land use/cover [23] ; weighted sum analysis technique [34] ; fuzzy analytical hierarchy process [30] [35] ; Snyder’s method of unit hydrograph, and land use/cover [36] .
In the current study, the authors employed the morphometric analysis method to:
1) Prioritize 76 fourth-order sub-watersheds for soil and water conservation based on the morphometric analysis method using GIS and RS.
2) Generate a spatial map illustrating the distribution of final priority classes representing the 76 sub-watersheds,
3) Explore the relationship of major components determined based on PCA to erosion risk parameters and the basic morphometric parameters.
4) Test the validity of the achieved priority classes by means of Discriminant Analysis, and to determine the nature of Discriminant functions in relation to the character of components,
5) Propose suitable soil and water conservation measures for the W. Mujib-Wala catchment.
Information on soil type, slope categories, and current land use/cover has been provided in order to help in suggesting appropriate soil conservation measures for sub-watersheds in each priority class. Prioritization in the present study demonstrates the role of powerful GIS, RS, and the morphometric analysis method in ranking different sub-watersheds in relation to the order in which they have to be taken for conservation measures [24] . Further, quantitative morphometric analysis of drainage networks is considered the key approach for understanding the hydro-morphological processes acting over drainage basins. Erosion risk parameters can be measured and calculated using DEM’s and Arc GIS software. Morphometric analysis of linear and shape parameters are the common approach developed and adopted to compute the compound parameter (Cp) on which prioritization of sub-watersheds is perceived [8] [10] [22] [32] . It has been argued in the recent past that erosion risk morphometric parameters are the most valid and consistent parameters in relation to erodibility and sub-watershed prioritization for conservation practice as compared to other methods of prioritization [13] .
2. Study Area
W. Mujib-Wala catchment represents the second largest watershed draining to the Rift, and discharging directly into the Dead Sea. The catchment is located between 30˚39' to 31˚ to 41'N latitude, and 35˚30' to 36˚30' East longitude (Figure 1). The watershed attains a maximum elevation of 1277 m (a.s.l) and a minimum of - 431 m (b.s.l) at the outlet of the wadi. The catchment has been subjected to progressive stages of rejuvenation due to continuous base-level changes along the Dead Sea, the lowest base level in the Jordan Rift [37] [38] . Thus, deeply dissected topography, incised drainage, and over-steepened and interrupted slopes are characteristic. Moreover, the catchment suffers from serious geo-hazards such as landslides activity, high sediment yield, severe soil erosion, and repetitive flooding [13] . Slope categories range from 0˚ - 10˚ in the east, to nearly vertical slopes bordering the canyons in the western part of the watershed, and the faulted-erosional escarpment overlooking the Dead Sea. Therefore, slope gradients have a direct impact on surface runoff and geomorphic processes over the demarcated sub-basins.
High hypsometric integral (HI) values ranging from 0.70% to 0.89% predominate in the watershed, indicating that the W. Mujib-Wala and the 76 sub-watersheds are at the youth-age stage of geomorphic development. Consequently, the sub-watersheds evince high soil erosion rates [3] , high sediment yield production [16] [17] , high possibility of flooding and landslide activity especially in the western part of the watershed [39] . Mean annual rainfall ranges from 346 mm at Madaba area, the north-western edge of the catchment, 164 mm at the W. Mujib weather station, to 331 mm and 335 mm at Qasr and Rabba town’ respectively. Daily temperatures range from a maximum of >40˚C in August (close to the Dead Sea), to a minimum of −5˚C in January close to Mazar town. Several rock units are exposed in the catchment, ranging from Cambrian sandstones to Quaternary fluvial terraces and lacustrine deposits. Rainfed cultivation is practiced on the summit surfaces and undulating terrain units receiving 250 - 350 mm of rainfall annually (Madaba, Qasr, and Rabba areas). Whereas irrigated farming from springs, Al-Haidan groundwater wells, and base flow dominated the lowlands (i.e., wadi floor terraces and flood plains) below W. Mujib and W. Wala dams. Scattered irrigated farming based on groundwater wells is characteristic of the eastern W. Mujib-Wala catchment. However, hydrological prospecting and water resources assessment are in progress at
present. The most promising sites for groundwater exploration in W. Wala were demarcated [40] . Agricultural land utilization occupies 7% of the total catchment. The western part of the catchment comprises a major component of the W. Mujib Natural Reserve. Palm trees, tamarisk trees, and oleander shrubs are observed in wadi beds.
3. Materials and Methods
3.1. Derivation of Erosion Risk Parameters
The morphometric analysis method was employed to conduct prioritization for the 76 sub-watersheds. Five linear morphometric parameters (Rb, Dd, Fs, Tr, and Lo) and five shape morphometric parameters (Rf, Bs, Re, Cc, and Re) (see Table 1) were calculated as a first step to compute the compound parameter (Cp) in order to categorize the sub-watersheds into priority classes. The linear parameters comprise the bifurcation ratio, drainage density, stream frequency, texture ratio, and length of overland flow.
The bifurcation ratio (Rb) refers to the ratio of the number of streams in lower order (Nu) to the number of the streams of the next higher order (Nu + 1), and it computed as:
, where (1)
= total number of stream segments of order -u-
= no. of segments of the next higher order [41] [42] .
High Rb values indicate high overland flow and early hydrograph peak with a high potential of possibility for flash flooding associated with severe rainstorms [43] , thus increasing soil erosion rates and sediment discharge in the main channel.
Drainage density (Dd) is introduced by Horton [41] as a significant indicator of the linear scale of landform elements in fluvially eroded topography. Dd is defined as the total length of streams in a drainage basin per unit area, or
![]()
Table 1. Morphometric characteristics of the 76 sub-watersheds
, where (2)
A is the basin area
is the total stream length
Stream frequency (Fs) is computed as the ratio of the total number of streams (Nu) of all orders in a watershed to the catchment area (A) [41] . It denotes the texture of the drainage network, and is expressed by the following equation:
, where (3)
is the total number of streams
A is the area of a drainage basin
Texture ratio (Tr) is defined as the ratio of the total number of streams of the first order (N1) to the perimeter of the drainage basin. Texture ratio is determined by:
where (4)
= the total number of streams of all orders
P = perimeter (km)
Length of overland flow (Lo) is the length of water over the ground before it gets concentrated into definite stream channels [41] , and is determined by the equation:
, where (5)
H = basin relief (Bh) m
= basin length (km)
parameter is related inversely to the average slope of the channel and is equivalent to the length of sheet flow to a large extent [41] .
The shape parameters comprise form factor, shape factor, elongation ratio, compactness coefficient, and circularity ratio.
Form factor (Rf) refers to the ratio of the area of drainage basin (A) to the square of the basin length (
) [41] [42] . It is computed according to the following:
where (6)
higher values of Rf indicate a more circular shape of a catchment, whereas small Rf values (<0.45) imply that the basin is elongated [41] . A catchment characterized by high Rf has high peak flows of shorter duration.
Shape factor (Bs) refers to the ratio of the square of the basin length to the area of the basin, or
(7)
Shape factor provides a notion towards the circular nature of the watershed. The greater the circular character, the greater is the fast response of the catchment following an intense rainstorm [44] .
Elongation ratio (Re) is denoted by
(8)
Variations in geology and climate cause prominent variation in Re values. Low values of Re indicate that the watershed is more elongated. When Re values approach 1.0, the shape of the watershed becomes circular [45] .
Compactness coefficient (Cc) parameter is developed by Gravelius [46] as a shape measure of a watershed. Cc is defined as the ratio of perimeter of a watershed to the circumference of the circle area, which is equal to the area of the watershed, or
where (9)
P = perimeter of the basin (km)
A = area of the basin (km2)
When the Cc value is 1.0, it denotes that the catchment is a perfect circle. If the value is 1.28, the basin is more square-shaped, whereas the basin is considered a very elongated one, when the Cc value > 3.0 [47] .
Circularity ratio (Rc) parameter has been developed by Miller [48] , and calculated according to:
, where (10)
(A) is the basin area, and (P) is the perimeter. If Rc is close 1.0, the shape of catchment is circular. Low, medium, and high values of Rc indicate young, mature, and old stages of geomorphic development of the watershed.
3.2. Tools, Data, and Multivariate Statistical Techniques
3.2.1. Digital Elevation Model and Morphometric Parameters
Morphometric analysis for prioritization was conducted using the ASTER Global Digital Elevation Model (TEM) of the National Aeronautic and Space Administration (NASA). The DEM is projected to the Universal Transverse Mercator (UTM) project system, Zone 36 N using Arc GIS (10.1) software. The W. Mujib-Wala, and the 76 sub-watersheds were delineated initially using topographic maps (Scale 1:50,000). Then ASTER DEM was employed to demarcate the final boundaries of the W. Mujib-Wala and the 76 sub-watersheds (Figure 2). The drainage network for the entire watershed and the sub-watershed were extracted using ASTER DEM (30 m resolution) and the Arc Hydro tool provided by Arc GIS software. Stream order was designated using the stream ordering system developed by Horton [41] and elaborated by Strahler (1964). The W. Mujib-Wala is classified as a seventh-order drainage basin (Figure 3) [42] . Fifteen morphometric
parameters were computed using DEM,
GIS
, and the standard formulae developed elsewhere [22] [41] [42] [45] [49] . Moreover, five basic parameters, five linear parameters, and five shape parameters were considered in prioritization of the 76 sub-watersheds based on morphometric analysis as illustrated earlier (Table 1). Linear parameter comprise the Bifurcation ratio (Rb), Drainage density (Dd), stream frequency (Fs), texture ratio (Tr), and length of overland flow (Lo). Shape parameters include: shape factor (Bs), circularity ratio (Rc), form factor (Rf), elongation ratio (Re), and compactness coefficient (Cc). Further, the calculated basic parameters are: basin area (A), perimeter (P), number of streams (Nu), basin length (Lb), and stream length (Lu).
3.2.2. Land Use/Cover, Soil, and Slope Categories Maps
Using ERDAS Imagine (2015, v. 15), LANDSATE 8 (July 2017), and supervised classification, the Maximum Likelihood Method of classification techniques was employed to classify Land use/cover guided by the classification system proposed by Anderson et al. [50] . The land use/cover classes recognized are: bare land; rainfed farming; irrigated farming, water body (Mujib dam and Wala dam), and irrigation pools, settlement, and major and minor roads. A soil map (Figure 4) was digitized from the National Soil and Land Use maps of Jordan [51] . A slope categories map was derived using ASTER DEM (Figure 5). A rem arkable variation exists in slopes. Slope categories of 0 - 5, 5 - 10, and 10 - 15 characterized the eastern part of the entire watershed, whereas, slope categories of 15 - 20, 20 - 30, 30 - 45, and >45 stand out in the western part of the watershed and the faulted-erosional escarpment overlooking the Dead Sea, and along the canyons and major tributaries downstream (Figure 6).
3.2.3. Multivariate Analysis
Principal Component Analysis (PCA) and Discriminant Analysis (DA) have been employed in the analysis. The PCA method is often used to analyze a large data matrix present study the 76 sub-watersheds in order to reduce the large number of variables to a small number of composite variables, or principal components (or factors) which are correlated with the erosion risk morphometric parameters. The technique is also designed to compute the correlation matrix to expose the highly loaded variables on each principal component, and how much the variability in the original parameters is explained by each principal component. The Eigen-value indicates how each of the identified components will fit with the data from all the morphometric parameters on all the PCs. Moreover, Discriminant Analysis was used to test statistically the validity of priority classes of sub-watersheds derived based on morphometric analysis, and to determine if they are significantly different from each other. This serves to guide in explaining the regional spatial variations among the fourth-order sub-basins in terms of prioritization. The development of powerful and cost-effective GIS and RS technology enables researchers to extract, measure,
compute and precisely process basic, linear, shape, and relief morphometric parameters of drainage basins . Simila rly, the availability of free access Digital Elevation Models, i.e. STRM and ASTER DEMS, with a reasonable resolution (30 m for ASTER DEM, and 90 m for STRM DEM)have enhanced quantitative analysis of drainage basin morphometry, thematic mapping, and function to expand the application of quantitative morphometry to other fields of research.
4. Results and Discussion
4.1. Morphometric Analysis
Basic parameters comprise: A, P, Nu, Lb, and Lu (Table 1). The basin area (A) is a basic morphometric parameter, and constitutes a major component in the hydrological processes [52] . Thus, it has been established by Chorley et al. [53] that the maximum discharge of flood per unit are, is inversely related to the size of the drainage basin. The total area W. Mujib-Wala is 6571.425 km2 (Table 2), and for the 76 sub-basins it varies from 12.12 km2 to 184.99 km2. The basin length (Lb) of W. Mujib-Wala is 136.8 km, and the perimeter is 564.484. The perimeters for the sub-watersheds range from 17.59 km (Sub-basin no. 23) to 153.38 km (Sub-basin no. 13). Sub-basin no. 13 is the longest in basin length (38.19 km), while Sub-basin no. 23 is the shortest (5.42 km). W. Mujib-Wala is classified as a seventh-order catchment, whereas the delineated 76 sub-catchments
![]()
Table 2. Morphometric characteristics of W. Mujib-Wala Watershed.
are of fourth order. The total stream length (Lu) of W. Mujib-Wala is 9216.34 km , and the first order streams account for 50.3% of the total stream length. The linear and shape parameters will be discussed in relation to their significance for the hydro-morphological properties of the sub-watersheds. The mean bifurcation ratio (Rbm) for W. Mujib-Wala is 4.324, and for the 76 sub-basins it ranges from 2.17 to 11.33 (Table 1). High Rbm values indicate the impact of structural distortions (faults system and dense lineaments) on drainage networks. The drainage density (Dd) for the entire W. Mujib-Wala is 1.403 which denotes moderate to well-drained catchments. An obvious variation in Dd values exists between the eastern and southern sub-basins (21 sub-basins, 28.9% of the total), and the western and the south western sub-basins (18 sub-basins, 23.7% of the total). Such variation is attributed to relatively high rainfall, and degradation of vegetation cover over the western sub-watersheds, while the southern and south eastern sub-basins are characterized by high relief, steep slopes, and influenced by the Kerak-Al-Fiha faults system and the associated dense lineaments. Thus, greater runoff, and more surface erosion is predicted. There are 39 sub-basins (51.3% of the total) with Dd values <1.4. The stream frequency (Fs) of W. Mujib-Wala is 1.209, and the highest values of Fs (1.77) are observed in sub-basin no. 6, while the lowest Fs values (1.07) are observed in sub-basins nos. 3 and 51. High Fs values are consistent with areas of a high density of lineaments. The texture ratio (Tr) of W. Mujib-Wala is 14.08, and for the 76 sub-basins, it ranges from 0.40 (sub-basin no. 61) to 1.97 (sub-basin no. 76). The values of Tr indicate that these sub-basins are of high runoff. The length of overland flow (Lo) for W. Mujib-Wala is 0.701, whereas, the Lo values for the 53 sub-basins range from 0.59 km (sub-basins no. 23) to 0.80 km (sub-basin no. 60). Further, the Rf value for W. Mujib-Wala is 0.351, and for the 76 sub-watersheds ranges from 0.09 (sub-basins no. 13), to 0.75 (sub-basin no. 47). Sixty one sub-basins have Rf values less than 0.4 (80% of the total), which implies that these sub-basins are elongated and more elongated in shape. Therefore, they are characterized as having a low peak flow of longer duration, and consequently exhibit lower probability for severe flooding [13] . The shape factor (Bs) for the entire W. Mujib-Wala is 2.849, while Bs values for the 76 sub-basins vary from 1.24 (sub-basin no. 16) to 11.59 (sub-basin 13) which indicate that elongated shape characterizes most of the sub -basins. The elongation ratio (Re) for Wadi Mujib-Wala catchment is 0.668, while it ranges from 0.33 (sub-basin no. 54), with seven sub-basins (Re > 0.95) nearly circular in shape, and sub-basin no. 16 is the only circular one. The compactness coefficient (Cc) value for W. Mujib-Wala is 3.929, whereas the Cc values for the 76 sub-basins vary from 2.30 (sub-basin no. 35) to 7.72 (sub-basin no. 13). Therefore high surface erosion dominates these sub-watersheds [41] . The circularity ratio (Rc) for the W. Mujib-Wala catchment is 0.259, whereas Rc values for the 76 sub-watersheds range from 0.10 (sub-basin no. 55) to 0.75 (sub-basin no. 35). Rc values indicate that W. Mujib-Wala and the 76 sub-basins are at the youth-age stage of geomorphic development, and most of them are elongated in shape.
4.2. Prioritization of Sub-Watersheds Based on Morphometric Analysis
Morphometric analysis of the 76 sub-watersheds of W. Mujib-Wala catchment revealed their relative properties in relation to erodibility and hydrological response. Further, the results of the morphometric analysis method clarified sub-watersheds in terms of priority classes, including higher erosion and soil loss-prone areas which required urgent suitable soil erosion control measures to protect the top soil from further erosion. Very high and high priority classes indicate the greater degree of erosion in specific sub-watersheds, the greater the need to improve available soil and water conservation practices or, the need, to install additional soil conservation structures. The “erosion risk morphometric parameters” utilized to prioritize the 76 sub-basins connected to the W. Mujib-Wala are: five linear parameters, e.g. bifurcation ratio (Rb), stream frequency (Fs), drainage density (Dd), length of overland flow (Lo), and texture ratio (Tr). Likewise, five shape parameters are employed in the process and include: shape factor (Bs), form factor (Rf), circularity ratio (Rc), compactness coefficient (Cc), and elongation ratio (Re). Based on the range of calculated compound parameter (Cp) values (Table 3), the 76 sub-watersheds were classified into four priority classes:
1) Very high priority (20 - 29)
2) High priority (30 - 39)
3) Moderate priority (40 - 49)
4) Low priority (50 - 60)
The spatial distribution of priority classes is illustrated in Figure 7 and Table 3. Out of the 76 sub-watersheds, 12 sub-basins (15.8% of the total) came under very high priority (sub-watersheds 6, 8, 10, 25, 32, 35, 47, 54, 60, 63, 64, and 67). The second category of sub-basins is assigned as high priority. It consists of 32 sub-basins (42.1% of the total) as follows: 3, 5, 11, 12, 14, 16, 19, 20, 21, 24, 26, 29, 34, 37, 38, 40, 41, 42, 52, 53, 55, 59, 62, 66, 68, 69, 71, 72, 73, 74, 75, and 76. The third class of sub-basins is designated as moderate priority class (Figure 7 and Table 3). It comprises of 25 sub-watersheds (32.9% of the total) as follows: 7, 9, 13, 15, 17, 22, 23, 27, 28, 30, 33, 36, 43, 44, 45, 48, 49, 50, 51, 56, 57, 58, 61, 65 and 70. Sub-basins of this category overlap in the eastern and central part of the W. Mujib-Wala watershed with sub-basins ranked mainly as high priority. About half of these sub-basins are located in the rainfed farming area in the western part of the catchment (Figure 8(a)), while the rest of the sub-basins are located within the degraded rangeland and bare land. Scattered irrigated agriculture is predominant based on pumping wells (Figure 8(b)). It is also practiced downstream of the W. Wala floodplain and fluvial terraces using surface water from Wala reservoir (Figure 8(c)), the extracted water from springs, and groundwater wells. Similarly, irrigated agriculture is growing rapidly down the
![]()
Table 3. Calculation of compound parameters and prioritized ranks based on morphometric analysis.
dam of the W. Mujib on the floodplain, and fluvial terraces using surface water from the reservoir and pumping wells as well. The fourth priority class of sub-basins is assigned a low priority (Figure 7 and Table 3). It consists of 7
sub-watersheds (9.2% of the total) as follows: 1, 2, 4, 18, 31, 39, and 46. Five of these sub-basins represent typical rainfed and irrigated farming areas on flat-undulating terrain in the northern part of the catchment east of Madaba, and the eastern part of the watershed. The expansion of cereals cultivation over the marginal areas and rangeland (annual rainfall varies from 200 to 250 mm) accelerate soil erosion [2] . Irrespective of scattered irrigated agriculture, the eastern part of W. Mujib-Wala forms a poor grazing land with high soil erodibility. Thus repetitive intense rainstorms over southern Jordan [54] initiate severe soil erosion where the vegetation cover is highly degraded. Additionally, degraded rangeland and bare land occupy large areas of the watershed due to the low amount of rainfall, marginality (the dominance of semiarid and arid conditions), which in turn accelerate soil erosion. All sub-watersheds ranked under very high and high priority have greater erosional potential with a high erosion risk. Consequently, they are considered potential areas for adapting soil conservation measures [21] . Rill and gully erosion, and landslide activity are common on steep slopes (15˚ - 25˚, and >30) where soft carbonate rocks of low shearing resistance are exposed. Sheet erosion is also active on gentle slopes (0˚ - 5˚) in the rainfed areas over the table lands bordering the canyons of W. Mujib-Wala downstream (i.e., the Qasr-Rabba, south east of Madaba, and Dhiban areas) where the annual rainfall ranges from 285 to 350 mm. Geological, morphological
![]()
Figure 8. Rainfed farming on gentle/undulating slopes (a); pumping wells irrigated farming (b) and irrigated farming based on the Wala reservoir (c). Source: Google Earth Bro 1/1/2017.
and anthropogenic factors account for high soil erosion rates. Among these are: the Kerak Al-Fiha fault system and the associated dens lineaments, the northwest-southeast fault system which influenced the central part of the W. Mujib-Wala, the W. Wala fault, the W. Az-Za’faran fault, W. Zerqa Ma’an fault, and the W. Al-Falij fault, and the connected dense lineaments [55] .
Morphological factors contributing to high soil erosion loss are: high relief, steep and long slopes of 10˚ - 15˚, 15˚ - 20˚, 20˚ - 30˚ and >30˚ slope categories. Sub-watersheds classified as high priority groups display high soil erosion rates and high sediment loads discharging into the Wala dam and the Mujib dam [16] [17] . These 44 sub-basins which are influenced by such processes (58% of the total) need urgent attention for implementing soil conservation practice. The long periods of resource exploitation, land use abuse, and the deterioration of vegetation cover have functioned to maximize soil erosion since the Neolithic and Chalcalithic periods. The existence of agricultural stone terraces over the rainfed highlands is an indicator of severe soil erosion at least since the Nabatean period some 3000 years ago [56] .
4.3. Proposed Soil and Water Conservation Measures
Morphometric analysis using erosion risk parameters was employed to prioritize 76 fourth-order sub-watersheds for soil and water conservation. The compound parameter values were calculated and the prioritization rating is given in Table 3. The effect of land use/cover, slope, soil, and relief were considered in order to propose appropriate soil and water conservation measures. Despite overlapping observed in the spatial distribution of sub-basins connected to each priority class, geological, morphological, and anthropogenic factors including history of land utilization (i.e. past and present abuse of land resources), were responsible for such a pattern of priority classes distribution. According to their ranking priority, sub-basin 54 with a compound parameter (Cp) value of 22.8 receives the highest priority rank, while the next in priority (Table 3) is sub-watershed 47, having a value of 23.3. Likewise, sub-basins nos. 6 and 10 rank third and fourth in the priorities with Cp values of 26.1 and 26.7 respectively.‟ Highest priority implies the greater degree of erosion potential associated in the particular area for adaptation of soil conservation measures” [21] . It is evident in this study that sub-watersheds nos. 6, 10, 25, 52, 60, and 64, for example, are designated as very high priority for soil conservation measures (Figure 7 and Table 3). Sub-basins nos. 6, 10, and 25 represent deeply dissected streams cut into Ajlune group limestone [37] , where the Chromic Camisol type developed. These sub-watersheds constitute desertified rangeland and bare land. Over-grazing, recurrent droughts, the expansion of towns and rural settlements at the expense of rainfed areas, and degradation of vegetation cover caused severe soil erosion, and high sediment yield discharge into the Wala and Mujib reservoirs. Similarly sub-watersheds 52, 60, and 64 represent undulating terrain developed on limestone, with slope categories of 0 - 5 and 5 - 10. The soil type is restricted to Yermic-Haplic Calcisol, the prime land use activity is poor grazing and scattered irrigated farming based on underground water which supports the livelihood of local people. Structural conservation measures, comprise, stone bunds, contour stone terraces, and water conservation measures; these are on both gentle and steep slopes connected to sub-watersheds classified under very high and high priority. Rangeland management is essential to protect the present degraded vegetation cover, and redeveloping of the natural vegetation by means of seeding specific sub-basins with suitable grasses, and planning for efficient range land management. Expanding tree-planting of drought-resistance species is also recommended [57] .
The high priority class of sub-watersheds occupy the northeast, east and southeastern part of the W. Mujib-Wala watershed. The main land use activity here is an overlapping of rainfed farming, rangeland, and scattered irrigated cultivation based on pumping wells. Sub-basins 3, 12, 16, 20, 34, 40, 55, 62, 68, 75 are an example (Figure 7 and Table 3). Exploitation of land resources over the last 3000 years [56] has contributed to several serious environmental problems, i.e., destruction of vegetation cover, severe soil erosion loss, continuous crop yield and livestock decline, and suffering of the local people. The most noticeable land use activity is the expansion of rainfed farming and irrigated agriculture in the west towards the marginal and rangeland in the east, and the increase in irrigated farming. Thus, a decrease in rangeland, and increase in bare land have occurred. Appropriate check dams, terracing and bunding are suitable conservation measures on sloppy land that are proposed, and efficient rangeland management is recommended to reduce soil erosion loss and sediment yields. Moderate and low priority sub-watersheds are located on the broad table-land bordering the deeply dissected wadis cut in Ajlune group limestone and basalt in the western part of the entire catchment, and the gentle-undulating terrain in the northern part. The dominant slope categories are: 0˚ - 5˚, 5˚ - 10˚, and 10˚ - 15˚. Rainfed and irrigated farming are the main land use activity. Examples of these 32 sub-basins (42% of the total area) are: 7, 17, 27, 28, 30, 48, 49 (moderate priority), and sub-basins nos. 1, 2, 4, 18, 46, (low priority) (Figure 7). Length and steepness of slopes in the moderate and low priority sub-watersheds can be modified through the construction of contour stone terraces locally termed “Mastaba” accompanied by tree-planting on slopes ranging from 3˚ to 25˚ [13] [58] . Stone terraces should be placed in long rows along the contour at various intervals depending on the length and steepness of slope [58] . Additionally, bench terraces can be constructed on slopes even steeper than 25˚, and when stones are not sufficiently available to build contours stone terraces. Such conservation techniques have been adopted by the rainfed farmers since the Nabatean period, some 3000 years ago [56] . Structural conservation measures should be integrated with technology enhancing farming practices, i.e., rotation, strip cropping, contour strip intercropping cultivation [59] , and contour plowing so as to reduce soil erosion loss, increase soil moisture at a farm, and improve crop productivity. The previous technique can be proposed for soil and water conservation in the sub-watersheds dominated by rainfed farming.
4.4. Multivariate Analysis
4.4.1. Principal Component Analysis (PCA)
Principal Component Analysis was performed to assess the effective basic and erosion risk morphometric parameters in watershed prioritization for soil and water conservation, and to illustrate the physical behavior and relationships between the basic and erosion risk parameters for a semi-arid catchment. Principal component analysis utilized as a component loading matrix which expresses numerically the level of relationship between them and the original erosion risk morphometric parameters. The weights of the original parameters in each component are called “loadings”, and each component is correlated with certain erosion risk parameters. Besides interpreting the processes that generate the observed relationship between the chosen variables, PCA also provides a simplified data matrix known as the component score (or weightings) matrix [60] . Through Principal Component Analysis, three efficient components were extracted, and the 16 basic and erosion risk parameters were reduced to these three components. The relationships between the rotated components and the original morphometric parameters are occasionally incorporated in the component loading matrix. These relationships are expressed in terms of the percentage contributed to the variance in the original data. It is also obvious that each component is more strongly correlated to some parameters that are considered more effective compared to others. Table 4 illustrates the Eigen-values, variance proportion, and cumulative proportion variance. Further, the results showed 88% of variances represented by 12 parameters which include five basic parameters, and one erosion risk parameter are highly loaded on PC1 and PC3. Most of the basic and erosion risk parameters exhibit loading values > 0.9, whereas Rc, Fs and Tr parameters display loading values < 0.9 (0.798, 0.898, and 0.703 respectively). The higher the loadings, the stronger the correlation. Three major components were resulted based on PCA analysis and accounts for 88% of the total variance explained by the basic, and erosion risk parameters and the Cp parameter. The most effective variable in PC1 to PC3 are shown by bold font in Table 5. It is clear that stream length (Lu), basin area (A), basin length (Lb), stream number (Nu), Perimeter (P), and circularity ratio (Rc) have the most variance proportions (41.85%), and turn to greatest effect on the first component. Likewise, drainage density (Dd), length of overland flow (Lo), texture ratio (Tr), Stream frequency (Fs), and the compound parameter (Cp) have the greatest effect on the second component (PC2), which comprises 25.1% of variance proportions. Additionally, form factor (Rf) and elongation ration (Re) have the greatest effect on the third component (PC3), which contributes 21.1% of variance proportions. Variable loadings show that one of the erosion risk parameters is accommodated within PC1, while six parameters including the Cp parameter are associated with PC2 and PC3 where both explain 46.2% of the total variance. The erosion risk parameters mentioned above explain more than half of the total variance, although the Rc parameter is associated with PC1. The component loadings/show that PC1 describes parameters related to the sub-watersheds drainage network (Lu, A, Lb, Nu, and P). Consequently, the PC1 is labeled as a “sub-watershed size “component. The second component corresponds to drainage density (Dd), length of overland flow (Lo), texture ratio (Tr), and stream frequency (Fs). Thus, the second component reflects the dissection/hydrographic component. The
![]()
Table 4. Total variance explained by the component.
Note: Figures in bold indicate the variables with strong loadings on the corresponding components (represent the reduced 16 variables).
third component corresponds to the form factor (Rf), and elongation ratio (Re); therefore this component refers to the sub-watershed “Shape component.”
4.4.2. Discriminant Analysis (DA)
It has been argued that prioritization based on the morphometric analysis method by utilizing the erosion risk morphometric parameters is time consuming in comparison with other approaches including the Principal Component Analysis approach, which allows for more effective parameters for prioritizing watersheds [24] . However, a prioritization for the Zarqa River has been carried out recently based on the morphometric analysis method and the PCA approach [13] . The output revealed that both methods did not produce similar results as stated earlier by Gajbhiye and Sharma [24] in a study on Shakkar River Catchment, Madhya Pradesh, India. Such disagreement in the results related to the two case studies is probably attributed to differences in physical conditions between Central India and Northern Jordan. The present study introduces a prioritization scheme for 76 sub-watersheds on the basis of the morphometric analysis method (the linear and shape parameters), and to test the achieved priority classes using Discriminant Analysis (DA) technique. The resultant priority classes displayed in Figure 7 and Table 3, are proved to be statistically valid. Consequently, the utilization of the morphometric analysis method is justified as a successful method of prioritization. Statistical validation also implies that erosion risk morphometric parameters are efficient parameters in prioritization of watersheds for soil and water conservation measures. The intention of statistical testing of sub-watersheds pertaining to the four priority classes is to test the hypothesis that there are significant differences between the four priority classes achieved earlier, and if the hypothesis can be accepted to establish a system of coordinate axis which discriminates between the four priority classes identified (1: low priority to 4: very high priority). Statistical analysis was conducted on four data matrices representing the four priority groups (i.e., 7 × 11; 25 × 11:32 × 11; and 12 × 11) with the associated ranking values (connected to linear and shape parameters), including the Cp scores. The F test of Wilks Lambda obtained is F ratio 174.9, with the degree of freedom V1 = 3 and V2 = 72. Referring to the table of percentage points of the F-distribution, with V1 = 3 and V2 = 72, it is found that at 99.9 percent of confidence, the tabulated value is 5.78, which is significantly exceeded by the computed F ratio (174.9). Subsequently, their is a great significant difference between each of the priority groups (very high, high, moderate, and low), and the four priority classes are completely separate and distinct. Moreover, 98.0 percent of the difference between the four the four priority classes is attributed to discriminant function 1 (93.7 percent) and discriminant function 2 (4.3 percent). Further, it was observed that the discriminant function 1 is positively correlated with six erosion risk morphometric parameters (the linear and shape parameters). Correlation values range from 0.413 to 0.999, and the Cp values are very strongly correlated with discriminant function 1 (0.970). By contrast, the correlation of discriminant function 2 with erosion risk parameters is relatively weak (0.17 - 0.237). The scores of each sub-basin of the priority groups (shown in Table 3, and illustrated in Figure 7) on the discriminant function 1 and 2 were plotted in Figure 9. The plot displays completely disconnected priority clusters. The present results show that
![]()
Figure 9. Results of discriminant function 1 and 2 grouping: the four priority classes are completely separated.
prioritization based on morphometric analysis is substantiated to be statistically valid, consistent and reliable, and of high capacity using GIS tools. The potential of the morphometric analysis approach as developed an elaborated earlier [8] [9] [22] is highly appreciated, and thus is, recommended for prioritization research.
5. Conclusions
Soil erosion by water has seriously threatened rainfed farming and rangeland over most of the terrain units of the W. Mujib-Wala watershed. High soil erosion rates have also increased sediment supply to the Mujib and Wala reservoirs during exceptionally heavy rainstorms, which are common in southern Jordan. The GIS-based morphometric analysis method and RS, were employed to prioritize 76 sub-watersheds, and the relationship of major components determined based on PCA to erosion risk morphometric parameters was explored. Then efficient conservation measures were suggested, especially for sub-watersheds where rainfed farming and grazing are practiced. Four priority classes were recognized: very high (12 sub-basins, 15.8% of the total), high (32 sub-basins, 42.1% of the total), moderate (25 sub-basins, 32.9% of the total, and low (7 sub-watersheds, 9.2% of the total). Moreover, 44 sub-watersheds (57.9% of the total) are ranked under very high and high priority, and are subject to high erosion risk, thus creating an urgent need for applying soil and water conservation measures so as to maintain rainfed farming and grazing sustainability. Based on Cp values and ranking priority, sub-watershed 54 with a Cp score of 22.8 receives the highest priority. The next in the priority list is sub-basin 47, having a value of 23.3 (Table 3). Likewise, sub-basins nos. 6 and 10 ranked third and fourth in the priorities, with Cp values of 26.1 and 26.7 respectively. By contrast, sub-basins nos. 15, 58 and 70 for example, are ranked as 47, 46, and 68, with Cp values 40.6, 40.4, and 49.7 respectively, and with moderate priority. Further, sub-basins nos. 1, 18, 39 are ranked as 72, 73, and with compound parameters at 51.7, 52.1, and 54.7 respectively (Table 3), and with moderate priority as well.
The W. Mujib-Wala watershed has been subjected to severe soil erosion over the last 3000 years, resulting in immense destruction of its vegetation cover. Supplementary information regarding current land use/cover, soil, slopes and topography were used as a guide in suggesting suitable soil and water conservation measures. The recommended measures were in accordance with priority ascribed in order to minimize negative impact of soil and land resources, rainfed farming, rangeland, and sedimentation in the Mujib and Wala dams. The expected advantages of expanding soil conservation measures, modernization of old soil conservation structures, enhancing farming practice, and rangeland management over sub-basins ranked as very high and high priority are manifested in the following: soil erosion loss control so as to protect soil from future erosion; minimizing sediment yield production to control sedimentation in the Mujib and Wala reservoirs; and reduced peak flow across the sub-watersheds and the entire W. Mujib-Wala catchment. Principal Component Analysis was used to reduce the original 16 basic and erosion risk morphometric parameters to three significant components which account 88% of the variance explained by the basic and erosion risk parameters. Out of ten parameters, six erosion risk variables were strongly correlated with PC1, and PC2 where both explained 46.2% of the total variance. One erosion risk parameter (circularity ratio) was highly correlated with PC1, and most of the erosion risk parameters exhibit loading values >0.9. The limitation of the morphometric analysis method of prioritization stated earlier was tested statistically using Discriminant Analysis. The results showed that the four priority classes significantly differ from each other; thus, prioritization based on the morphometric analysis approach is consistent, reliable, accepted, and of high capacity using RS and GIS technology. The present results are aimed so as to assist decision-makers in identifying priority sub-basins which need immediate adoption of appropriate conservation measures, and land management practices.