Land degradation in Jordan is induced mainly by runoff soil erosion which was active prehistorically and historically in the rainfed highlands region. Geoarchaeological studies conducted on Paleo-soil erosion in the alluvial deposits of W.Wala in historic time, suggest that historical soil erosion, agricultural intensification, and agricultural terraces, were prevalent over the Levantine highland since the Iron Age  . Further, destruction of vegetation cover, intensive exploitation of land resources, and land use abuse entailed high soil erosion rates since the Neolithic and Chalcolithic periods. The presence of field soil erosion features and old agricultural terraces, indicate that the highlands of Jordan have experienced severe soil erosion at least since the Nabatean period, 3000 years ago  . Since the 1960s, quantitative and qualitative studies have been performed on soil erosion and conservation in Jordan  -  . Several studies also reported the occurrence of repetitive heavy rainstorms which induced excessive soil erosion, landslide activity and flooding    . The recoded maximum rainfall intensity in 24 hrs ranges between 100 mm and 150 mm    , and the maximum daily intensity (mm/hr−1) varies from 2 and 6.6 mm/hr−1  . Exceptional daily rainfall intensity recorded was 15 mm/hr−1 during the 1960s. Several surveys were performed on soil erosion and conservation for Wadi Ziqlab (Northern Jordan)  , Wadi Shueib and Wadi Kufrein (Central Jordan)  , and Wadi Hasa (Southern Jordan)  . In situ field measurements on soil erosion were carried out across different bioclimatic regions of Jordan, i.e., sub-humid Mediterranean (Salt, Kufranja, and Jerash areas), semi-arid (Muwaqar and Shuwbak-W. Musa-Petra area), and arid (Azraq area) using different methods/techniques, i.e., splash, sheet, and runoff erosion  -  . The high estimated soil erosion loss for Wadi Kufranja  , and Wadi Kerak  watersheds using the RUSLE model, reveals that continuous deterioration of top productive soil. Thus, high soil erosion rates seriously endanger the rainfed agricultural sustainability, and the present and future proposed dams in the highlands of Jordan. The predicted average annual sediment yield for Wadi Wala and Wadi Mujib using SWAT model, and the Zarqa River using AGNPS erosion model is also considered a critical threat to W.Wala and W.Mujib reservoirs, and the KTD of the Zarqa River  -  . Annual and seasonal rainfall erosivity for Jordan have been estimated based on three regression models, and using annual and seasonal precipitation data over a long record (40 - 53 years) pertaining to 40 weather stations  . Spatial distribution of rainfall erosivity in Jordan is controlled largely by morphological (relief and slope) and climatic factors. Lowest R-Values are found in arid areas (annual rainfall is below 100 mm), whereas the highest R-values are found in the northern highlands (annual rainfall approaches 650 mm). The correlation between annual erosivity and the annual and seasonal precipitation (mm) exhibits a very strong relationship (R range from 0.964 to 1.0, and all correlations are significant at 0.01 level, 2 tailed test). Prioritization of sub-watersheds based on morphometric analysis method using GIS and RS techniques, are found to be efficient for large areas (watershed, or sub-watershed scale), and it is fundamental to provide proper soil conservation plan. Adaptation of conservation measures with respect to a priority approach will minimize soil erosion loss, thus increasing moisture in soil profile, and in turn reducing the impact of drought and the possibility of flooding. Finally, prioritization is important for positioning suitable sites for water harvesting, and designing efficient water harvesting structures  -  .
In the light of striking high soil erosion rates and sediment yield loads, distinct sub-basins exhibit potential areas to promote conservation activity, and must be prioritized immediately for soil and water conservation practices to maintain future agricultural sustainability. Quantitative morphometric analysis for watersheds was utilized recently as a successful approach to prioritize sub-basins for soil and water conservation measures        -  . Morphometric erosion risk parameters represented by linear and shape variables must be computed to prioritize sub-basins for soil conservation. In the recent past, it has been argued that linear parameters maintain a direct relationship with erodibility. Therefore, the highest value ranked as 2 and so on. Alternatively, the shape parameters have an inverse relation with erodibility, thus, the lower their values, the greater the erodibility     . Consequently, the lowest value of shape parameter was rated as rank 1 and the second lowest as rank 2 and so on. Compound parameter (Cp) was computed by adding up all the ranks of linear variables, as well as shape variables, and then, dividing by the number of all parameters.
Following the rating of every single morphometric variable, the ranking values for all linear and shape variables referring to each sub-basin are added up for each of the sub-basins to achieve a compound parameter (Cp) score based on the average value of these variables. Furthermore, the sub-basins having the lowest compound parameter score were designated the highest priority, the next higher value was assigned as second priority and so on  . 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  . Several methods of analysis were employed to perform watershed prioritization. Among these methods is: morphometric analysis based on linear and shape parameters    . Moreover, morphometric analysis combined with sediment yield index (SYI) model, and sediment product rate (SPR) was also utilized  . Other studies employed morphometric analysis, land use/land cover variables, and sediment yield index  . Snyder’s synthetic unit hydrograph method has incorporated in the analysis by Singh and Singh  . Several prioritization studies combined morphometric analysis with the USLE and RUSLE models for soil loss estimation, or soil erosion susceptibility assessment    . Others adopted morphometric indices and Fuzzy Analytical Hierarchy, including the Saaty Analytical Hierarchy Process   . Morphometric analysis of linear parameters (bifurcation ratio, drainage density (km/km2), texture ratio, stream frequency, and overland flow (Km)); and shape parameters (elongation ratio, circularity ratio, shape factor, form factor, and compactness coefficient) are the common parameters employed to calculate the compound parameter (Cp) on which prioritization of sub-watersheds is perceived. It has been reported  that the proposed linear and shape morphometric parameters for prioritization are consistent in relation to erodibility. However, the limitations argued regarding morphometric analysis method in prioritization were examined in the present study, and the validity of the priority classes was tested statistically using Discriminant Analysis. It is intended in this study to:
1) prioritize the 23 sub-watersheds of W.Wala using GIS-based morphometric analysis to identify sub-basins of high potential for applying soil and water conservation measures,
2) generate spatial maps displaying the final priority classes by means of GIS tools,
3) propose appropriate soil and water conservation measures for W.Wala in order to minimize soil erosion rates and sediment yield production,
4) test the validity of achieved priority classes (sub-watershed groups) through morphometric analysis by means of Discriminant Analysis (DA).
W.Wala is considered as a promising agricultural watershed, and water resources development, morphometric analysis and the retained information are significant for future soil water conservation planning and water resources management.
2. Study Area
The W.Wala catchment comprises the upper part of W. Mujib-Wala watershed. The catchment has a triangular shape with an area of 2063.6 km2. It is located at 35˚65' to 36˚30'E longitude, and 31˚55' to 31˚90'N latitude (Figure 1). Elevation varies from −327 m (b.s.l) at the confluent point where W. Wala merges with W. Mujib (3 km before the wadi system discharge into the Dead Sea) to 1007 m (a.s.l) northwest of the watershed (Figure 2). Successive lowering of the Dead Sea base level as a result of uplifting of the eastern shoulder of the Dead Sea Rift, and renewed downward movement of the Rift (during late Tertiary and Quaternary tectonics)  caused progressive rejuvenation, river incision and down cutting of W.Wala. Consequently, undulation to moderate slopes (0˚ - 5˚), (5˚ - 10˚), and (10˚ - 15˚) dominate the eastern part of the catchment, whereas, steep slopes (>35) and highly dissected topography characterize the western part (Figure 3). High HI values (70.0% - 89.0%) predominate in the catchment, indicating that W.Wala watershed and the 23 sub-basins are at the youth-age stage of geomorphic evolution. Accordingly, the sub-watersheds have a high possibility for
Figure 1. Location of the study area.
Figure 2. DEM of W.Wala.
flooding, high soil erosion rates, high sediment load, and landsliding activity. Cretaceous carbonate rocks outcrop in most of the W.Wala watershed. The massive limestone unit of the Turonian age is the oldest rock exposed. The average thickness of this unit is 67 m.
The massive chalk limestone unit is the youngest, and ranges in age from early Paleocene the Middle Eocene  . Basalt flows of the Pleistocene age are exposed in the upper W.Wala. Quaternary deposits of fluviatile and lacustrine gravels cover part of the deep tributaries of the catchment. The most productive aquifers in the watershed are those restricted to the Amman-Wadi As Sir Limestone formation (Upper Cretaceous) termed A7-B1  . In 2003, a reservoir was constructed on W.Wala with a capacity of 10 MCM, and the future plan is to
Figure 3. Stream order.
raise the height of the reservoir to store 26 MCM. The purpose of the reservoir is to recharge the groundwater, to provide water for springs and pumping wells, and to practice irrigated agriculture on a larger scale.
Hydrological prospecting and assessment of water resources is in progress at present. The most promising sites for groundwater exploration were demarcated  . The climate is classified as dry-Mediterranean, with relatively cold winters and hot summers. The canyon shape tributaries close to the Dead Sea are arid. Average annual rainfall ranges from 346 mm at Madaba (a few kilometers to the northwest of the catchment) to 282 mm at Dhiban, and 266 mm at the W.Wala weather station. Rainfall is concentrated in October and March. Large seasonal variations in temperature are evident, where daily temperatures range from a maximum ≈ 40˚C in August to a minimum of −5˚C in January. Rainfed cultivation of cereals (wheat and barley) is practiced in 38% of the total area of the catchment, and open rangeland constitutes 47% of the total area. Some 7% of the total area is urban  .
3. Materials and Methods
3.1. Extraction of Morphometric Parameters
The bifurcation ratio (Rb) is defined as the ratio of streams number of a given order to the number of the streams of the next higher order  , and it is computed by:
= total number of stream segments of order “u”,
= no. of segments of the next higher order.
Bifurcation ratio is elaborated as an index of relief and dissection. Rb values of drainage basins vary from 2 for flat/undulating or rolling terrain, to 6 for catchment where drainage network is distorted heavily by geological structure    . High Rb values indicate high overland flow and an early hydrograph peak with a high potential of susceptibility to flash flooding during exceptionally heavy rainstorms  , which in turn increase soil erosion rates and sediment discharge in the main channel.
Drainage density (Dd) is computed as the total length of streams in a drainage basin per unit area    , or
= the basin area,
= is the total stream length.
Dd is a measure of topographic dissection and runoff potential of the drainage basin. High Dd value denotes high runoff, a quick stream response, and in turn, a low infiltration rate. Whereas low drainage density of a catchment implies low runoff and high infiltration  .
Stream frequency (Fs) is defined as the ratio of the total number of streams (Nu) of all orders in a catchment to the watershed area (A)  . It is expressed by the following equation:
Fs are positively correlated with Dd values of a watershed. Low Fs values indicates a relatively low infiltration rate of surface water, and thus the ground water potential is relatively low  . High stream frequency denotes more infiltration, and thus high groundwater potential.
Texture ratio (Tr) Tr is refers to the ratio of the total number of streams of the first order (N1) to the perimeter of the drainage basin. It is considered a highly significant parameter in drainage basin morphometry. Tr variable relies on slope materials and lithology, infiltration capacity and topographic relief  . Texture ratio is determined by:
= the total no. of streams of all orders,
P = perimeter (km).
Length of overland flow (Lo) is determined by the equation:
, where (5)
H = basin relief (Bh) m,
= basin length (km).
Lo is the length of water over the ground before it is concentrated into definite stream channels  . Lo parameter is one of the most significant independent variables affecting both hydrologic and hydrographic development of drainage basins. This parameter is related inversely to the average slope of the channel and is equivalent to a large extent to the length of sheet flow  .
Form factor (Rf) is calculated according to Horton  :
Rf refers to the ratio of the area of drainage basin to the square of the basin length  . Higher values of form factor denote a more circular shape of a catchment, while smaller Rf values (<0.45) imply that the basin is elongated in shape.
Shape factor (Bs) represent the ratio of the square of the basin length to the area of the basin, or
Shape factor provides a conception regarding the circular character of the catchment. The greater the circular character, the greater the fast response of watersheds following an intense heavy rainstorm  .
Elongation ratio (Re) is defined by Schumm  according to the following equation:
Low values of Re denote that the catchment is more elongated. Where the Re value approaches 1.0, the shape of the watershed becomes a circular shape  .
Compactness coefficient (Cc) is developed by Gravelius  , and defined as the ratio of perimeter of a watershed to circumference of the circle area, which is equal to the area of the watershed. Cc is computed according to the following equation:
P = perimeter of the basin (km),
A = area of the basin (km2).
When Cc value is 1, it indicates that the watershed is a perfect circle. If the Cc value is 1.28, the basin is more square shaped, while the basin is considered highly elongated, when the Cc value > 3.0  .
The circularity ratio (Rc) of a catchment is calculated according to:
, where, (10)
(A) is the basin area, and (P) is the perimeter  . If Rc is close to 1, the shape of a catchment is circular. Low, medium, and high values of Rc indicate young, mature, and old stages of geomorphic development of the catchment respectively.
3.2. Tools, Data Used, and Statistical Techniques
Morphometric analysis for prioritization of the 23 fourth-order sub-basins related to W.Wala watershed was conducting using topo sheets, ASTER DEM and Arc GIS software. Topographic maps of scale 1:50,000 were acquired from the Royal Jordanian National Geographic Centre (Amman). They were scanned, geo-referenced, and converted to a zone 36N projection system using Arc GIS 10.1, and the associated packages. The W.Wala watershed and the 23 sub-basins were delineated using the topo sheets. Then ASTER DEM (30 m resolution) was utilized to extract the drainage networks using the Arc Hydro tool. Stream order was designated using the stream ordering system developed by Horton  and elaborated by Strahler  . The W.Wala watershed is classified as a sixth-order drainage basin (Figure 3). Fifteen morphometric parameters were calculated using DEM, GIS, and the mathematical equations developed elsewhere    . Five basic parameters were derived to illustrate the morphometric characteristics of the 23 sub-basins, Area (A), Basin length (Lb), Perimeter (P), Stream order (u), and Stream Length (Lu). Furthermore, five linear parameters, and five shape parameters were considered in prioritization of the 23 sub-basins based on morphometric analysis as exemplified earlier. These are: Bifurcation ratio (Rb), Drainage density (Dd), Stream frequency (Fs), Texture ratio (Tr), and Length of overland flow (Lo). Moreover, the five shape parameters used in this research are: Shape factor (Bs), Form factor (Rf), Circularity ratio (Rc), Elongation ratio (Re), and compactness coefficient (Cc). All morphometric parameters were measured directly from the DEM using GIS software. A land use/land cover map (Figure 4) was compiled using ERDAS IMAGINE 2015 (v.15), LANDSAT 8 (July 2017), and supervised classification. The Maximum likelihood Method of classification techniques was employed to classify land use/cover, based on the classification system proposed by Anderson et al.  . A soil map (Figure 5) was digitized from the 1994 National soil survey maps and reports pertaining to the National soil and land use maps  . A slope categories map was generated using ASTER DEM (Figure 5). A remarkable variation exists in slope categories. Slope categories of 0˚ - 5˚, 5˚ - 10˚, and 10˚ - 15˚ dominated the eastern part of
Figure 4. Land use/land cover.
Figure 5. Soils of W.Wala.
the catchment. Whereas, slope categories 15˚ - 20˚, 20˚ - 30˚, and >30˚ stand out in the western part, with the presence of cliffs (>45˚) along the canyon tributaries downstream. The development of an efficient and cost, effective GIS and remote sensing techniques enables researchers to extract, measure, compute and process precisely basic, linear, shape, and relief morphometric parameters of drainage basins. Moreover, the availability of free access Digital Elevation Models (i.e., STRM and ASTER DEMs) with a reasonable resolution (90 m and 30 m respectively) have improved the quantitative analysis approach of drainage networks, and morphometric mapping, thus, expanding the application of morphometric analysis to other fields of research. Discriminant Analysis (DA) was employed to test statistically the validity of priority groups of sub-basins generated based on morphometric analysis, and to determine if they are significantly different from each other, and to help in explaining regional spatial difference among the fourth-order sub-watersheds in term of prioritization. Through this technique, each sub-watershed can be related to a set of discriminant functions by calculating its “scores” on those functions. The mean value of these scores, on each discriminant function, can be calculated for each priority class. The mean values are then utilized to test if the priority classes are significantly different from each other, or in other words, if the resultant priority classes are accepted or not.
4. Results and Discussion
4.1. Morphometric Analysis
4.1.1. Basic Parameters
The basin area (A) is a principal component in hydrological processes  . In this regard, Chorley et al.  argued that the maximum discharge of flood per unit area is inversely related to the size of the drainage basin. The total area of Wadi Wala is 2063.6 km2, and for the 23 sub-watersheds, it ranges from 18 km2 to 184.99 km2. The basin length (Lb) corresponds to the maximum length of the watershed and sub-basins measured parallel to the main drainage line. The length of the W.Wala basin is 88.8 km, whereas the lengths of the sub-basins vary from 5.423 to 38.186 km. The perimeter (P) of W.Wala is 403.9 km, and the perimeters of the sub-basins vary from 5.423 to 38.186 (Table 1). Sub-basins no.23 represents the shortest, but with the longest perimeter, while sub-watershed no. 13 is the longest, but with the highest perimeter. In terms of area, sub-basin 12 is the largest, and sub-basin 10 is the smallest. Nevertheless, the greater area sub-basins and the longest are located close to the northwestern, southeastern, and the eastern borders of the main catchment. Whereas the shortest in length and smallest in area are located in the rejuvenated belt due to the dominance of steep slopes and topographic dissection. The W.Wala watershed is classified as a sixth-order basin, while all the delineated 23 sub-basins are of fourth-order. Stream length (Lu) is measured from the origin of a stream to the drainage divide. The total stream length of W.Wala is 2839.2 km, and the first-order streams represent 49.9% of the total stream length. The following linear and shape morphometric parameters will be discussed with reference to their significance to morphological and hydrological properties, and implications for soil and water conservation.
4.1.2. Linear Parameters
The morphometric parameters are: bifurcation ratio (Rb), drainage density (Dd), stream frequency (Fs), texture ratio (Tr), and length of overland flow (Lo).
Bifurcation ratio (Rb)
The bifurcation ratio is developed by Horton  as a morphological index for relief and dissection. The Rb value for W.Wala is 4.5, and for the 23 sub-basins varies from 2 to 7 (Table 1), with an average of 4.55. High Rb values demonstrate that most of the sub-basins are crucially affected by tectonics and structural disturbances caused by the W.Wala fault, the W. Az-Za’faran fault, the W. ZerqaMa’an fault, and the W. Al-Falij fault, and the associated dense lineaments  .
Drainage density (Dd)
Drainage density value for W.Wala watershed is 1.385, and for the 23 sub-watersheds varies from 1.19 to 1.77 which indicates moderate to well-drained watersheds. Slight variation is seen in Dd values between the upper sub-watersheds (x = 1.43), and the lower sub-watershed (x = 1.32)  . Higher drainage density in the upper reaches is probably attributed to the availability of rainfall, high relief, steep slopes,
Table 1. Morphometric characteristics of the 23 sub-basins.
and overgrazed areas which resulted in greater runoff and thus, more erosion and surface drainage lines  .
Stream frequency (Fs)
Fs values are positively correlated with drainage density of all sub-basins. Therefore, any increase in stream population led to an increase in Dd value  . Low Dd values indicate a low infiltration rate of surface water, thus, suggesting relatively low groundwater potential  . The lowest value of Fs is observed in sub-basin 7 (Fs = 1.064), whereas the highest value is observed in sub-basin 6 (Fs = 1.771). High Fs values are characteristic of areas with high density of lineaments associated with major faults, i.e., the W.Wala and W. Rumeil faults  .
Texture ratio (Tr)
Texture ratio is considered an important factor in drainage basin morphometry. The Tr for W.Wala is 4.8, and for the 23 sub-basins, it ranges from 1.36 (sub-basin no. 2) to 2.8 for sub-basin no. 10. Tr values indicate that the watersheds are of high runoff.
Length of overland flow (Lo)
The length of overland flow is one of the most significant independent parameters affecting the hydrographic and hydrologic development of a drainage basin  . The length of overland flow for the 23 sub-basins varies from 0.594 (sub-basins no. 23) to 0.791 for sub-basin no. 11.
4.1.3. Shape Parameter
Form factor (Rf)
Form factor is a dimensionless property and used as a quantitative expression of the shape of watersheds  . Catchments with high Rf values experience high peak flow of short duration. By contrast, an elongated catchment with low form factor has a low peak flows of longer duration. The Rf value for W.Wala is 0.286, and for the 23 sub-basins it ranges from 0.086 (sub-basin no. 13), to 0.809 (sub-basin no. 16). Most of sub-watersheds have Rf values varying from 0.3 to 0.5, which indicates that these sub-basins are elongated and more elongated, shape is associated with low peak flow of longer duration, and thus, of less probability for severe flooding  .
Shape factor (Bs)
The shape factor for the entire W.Wala is 3.821, whereas, the 23 sub-basins exhibit a range from 1.24 (sub-basin no. 16) to 11.58 (sub-basin no. 13), which implies that elongated shapes dominate most of the sub-basins.
Elongation ratio (Re)
Re values close to 1.0 are characteristic of watersheds with low and simple relief, while values range from 0.6 to 0.8 representing catchments with high rugged relief and steep slopes. The elongation ratio for W.Wala catchment is 0.577, while it ranges from 0.331 (sub-basin no. 13) to 1.0, where sub-basin no. 16 is the only circular basin.
Compactness coefficient (Cc)
The Cc parameter is independent of size of the catchment and dependent mainly on slope. Low i values imply more elongation and high erosion   . The Cc value for W.Wala watershed is 2.51, whereas the Cc values for the 23 sub-basins vary from 2.85 (sub-watershed no. 9) and 7.72 for sub-basin no. 17, which reflects high surface erosion.
Circularity ratio (Rc)
Circularity ratio is the most useful shape morphometric measure in correlation with stream discharge. Rc is influenced by geology, morphology (relief and slope), climate, land use/land cover of the catchment  . The circularity ratio for W.Wala is 0.159, whereas, the Rc values for the 23 sub-basins range from 0.126 (sub-basin no. 12) to 0.492 (sub-basin no. 23). These values imply that W.Wala and the sub-basins are at the youth-age stage of geomorphic development, and they are elongated in shape.
4.2. Prioritization of Sub-Basins Based on Morphometric Analysis and Proposed Conservation Measures
Recently, morphometric analysis was employed broadly to assess watersheds of different size (sub-watersheds, mini-watersheds, and micro watersheds) for soil and water conservation          . The linear and shape morphometric parameters termed “erosion risk parameters”   were utilized for prioritization of the 23 sub-basins connected with W.Wala, and designated 1 - 23 (Figure 6). The linear parameters include(a): bifurcation ratio (Rb), drainage density (Dd), stream frequency (Fs), texture ratio (Tr), and length of overland flow (Lo); and the shape parameters include (6): form factor (Rf), shape factor (Bs), elongation ratio (Re), compactness coefficient (Cc), and circularity ratio (Rc). Based on the range of computed compound parameter (Cp) values (Table 2), the 23 sub-basins of W.Wala were classified into three priority groups:
1) High priority (3 - 4),
2) Moderate Priority (4 - 5),
3) Low priority (5 - 6).
Using the Arc GIS tool, the spatial distribution of the 23 sub-basins groups based on the compound parameter values (Cp) was generated. Figure 7 illustrates the priority classes of sub-basins for soil and water conservation.
Prioritization was successful in classifying the 23 sub-watersheds into three groups of high, moderate, and low priority for soil conservation measures. The three groups were determined mainly by land use/land cover, and historical exploitation of soil resources. Other factors such as morphology (relief and slope),
Figure 6. Fourth-order sub-basins.
Figure 7. Priority of the 23 sub-basins based on morphometric analysis.
and geomorphic agents, i.e., rejuvenation are also influential, although its impact varies spatially from west to east over the W.Wala catchment. The long periods of human intervention, land use abuse, and the destruction of vegetation cover was decisive in maximizing soil erosion over the sub-basins which traditionally have been utilized for rainfed farming, or rangeland transformed to rainfed cultivation since the 1950s in the eastern part of the watershed. Spatially, sub-watersheds are characterized by high priority and are found in the rainfed areas, which occupy a substantial part of the watershed, over-grazed range land, bare land, and the western rejuvenated belt of the watershed.
Accordingly, the highest priority reveals the great degree of erosion risk associated with particular sub-basins, and thus, represents a potential area for adaptation soil and water conservation measures  . Out of the 23 sub-basins of W.Wala, 6 sub-basins (26% of the total) (sub-basins nos. 6, 8, 10, 11, 19, 21) are classified as high priority (Figure 7). Eight sub-basins (34.8% of total) (sub-basins nos. 2, 5, 12, 13, 14, 15, 16, 29) come under moderate priority (Figure 7) for soil and water conservation. It is obvious that 62% of the sub-basins are grouped under high and moderate priority for soil and water conservation. Eleven of these sub-basins are located at the eastern-southeastern part of W.Wala watershed, and three sub-basins are located at the northwestern part of the watershed (Figure 7). Their distribution is found in accordance with the land use/land cover pattern, especially overgrazed rangeland, and bare lands where the vegetation cover has been highly degraded (Figure 4).
These sub-basins are characterized by highest relief over the entire watershed (800 > 900 m a.s.l), where slope categories of 0 - 5, 5 - 10, and 10 - 15 are
Table 2. Calculation of compound parameters and prioritized ranks based on morphometric analysis.
predominant. The destruction of vegetation cover, caused severe soil erosion since the Neolithic and Chalcolithic periods. Indicators of soil erosion and the existence of agricultural stone terraces, imply that the highlands of rainfed farming (including W.Wala) experienced severe soil erosion at least since the Nabatean period, some 3000 years ago  . High relief ratio (Rr) has also been recorded for several sub-basins pertaining to high and moderate priority for soil and water conservation (i.e., sub-basins nos. 6, 10, 11, and 16) which can be considered an essential factor in the assessment of soil erosion and sediment load in W.Wala watershed. Furthermore, high hypsometric integral (HI) values for sub-basins (e.g. nos. 2, 5, 6, 8, 12) classified as high and moderate priority vary from 70% to 89%, which indicate that these drainage basins are influenced by rejuvenation processes, and are considered at the youth-age stage of geomorphic development. Consequently, surface erosion rates are considerably high  .
Nine sub-basins (39.2% of the total) (sub-basins nos. 1, 3, 4, 7, 9, 17, 18, 22, and 23) are assigned to be of low priority for soil and water conservation. Three sub-watersheds out of nine consist mainly of wide and flat, undulating table lands overlooking the deep canyons downstream which occupied the rejuvenated belt. Other sub-watersheds (nos. 9, 7, 4, 3, and 1) extend to the northern/northeast part of W.Wala watershed.
The nine sub-watersheds constitute a major part of the rainfed agricultural region of Jordan, where the cultivation of wheat and barley is practiced over these sub-watersheds. Here, the average annual rainfall varies from 250 mm to 350 mm. Scattered irrigated agriculture based on groundwater wells is practiced over the sub-basins categorized as being of high and moderate priority for soil and water conservation. Irrigated farming is also practiced downstream of W.Wala (Al-Haidan area) based on extracting the water from the springs and Al-Haidan groundwater wells (Figure 8(a)). High soil erodibility of Calcaric Cambisol and Eurtic Cambisols soil type (Figure 5), steep slopes (14˚ - 20˚), poor conservation measures and degraded vegetation cover accelerates soil erosion, and exposed the sub-watersheds to high rates of soil erosion loss. However, the transformation of vast areas since the 1990s from rangeland into farming practice, and the expansion of rainfed “mixed farming” towards the marginal areas (Figure 8(b)), accelerate soil erosion and have increase the susceptibility of soil erosion over these sub-watersheds although, they rank under low priority for soil and water conservation. The cultivated land with poor conservation measures exhibits a higher rate of soil erosion and decline in soil fertility. Consequently, all sub-basins ranked under high and moderate priority are vulnerable to soil erosion, thus, they should be prioritized for conservation. Even over the sub-basins with mixed rainfed farming areas, traditional up and down slope tillage without conservation measures (i.e., contouring and terracing) is common in W.Wala. Thus, priority must be given to the protection of sub-watersheds exposed to high soil erosion rates, in order to reduce runoff coefficient, increase infiltration rates and soil moisture in soil profile. Reduction in soil erosion loss has been verified recently by Al-Sheriadeh and Al-Hamdan  where they concluded that in the years 1987-1990; the estimated sediments yields were less than the actual yields in the Kink Talal Reservoir of the Zerqa River. They attributed this to the implementation of a soil conservation program launched by the government of Jordan, in cooperation with the German government. Further, Al-Alawi and, Abujamous  estimated the average annual sediment loss in Central Jordan (Salt area) at 78 ton∙ha−1∙year−1 before the construction of stone terraces. Twenty years later following construction and tree planting, the estimated average soil loss have been reduced considerably to an average of 33 ton∙ha−1∙year−1 (35.7%). Such encouraging results emphasize the need to execute intensive soil conservation programs over the highland to reduce soil erosion rates, in accordance with enhancing practiced soil conservation techniques, and to modify RUSLE parameters (i.e., C, P, and LS factors which are considered key factors in soil
Figure 8. Irrigated farming (W.Wala flood plain) (a), rainfed cultivation on the summits, and pumping well irrigated farming (b). Source: Google Earth bro 1/1/2017.
erosion). It has been argued that RUSLE parameters can be modified significantly through altering farmers environmental practices  with the support of local governments. Slope length and steepness factor (LS), the conservation practice factor (P), and the cover and management factor (C), can be improved noticeably to reduce soil erosion loss and to conserve water in soil on the farm to maintain crop productivity. The LS factor can be modified mainly through the construction of what is locally termed mastaba, or contour stone terraces associated with tree planting on slopes ranging from 5˚ - 25˚. Stone terraces should be aligned in long rows following the contours at various intervals depending on length and steepness of slope. On steep slopes (>25˚) with a shortage of stones, bench terraces are more appropriate   . Terrace farming in the W.Wala highlands has been practiced since some 3000 years ago  . The construction of check dams across the gullies (1 - 3 m of depth) can decrease runoff coefficient and soil erosion rates under different land utilization  . Structural soil conservation measures can be applied on gentle slopes (0˚ - 5˚, 5˚ - 10˚), and moderate to steep slopes (10˚ - 15˚, 15˚ - 20˚) characterize the sub-watersheds of high priority for soil conservation measures (sub-basins nos. 6, 8, 19, and 21), and sub-basins with moderate priority (sub-basins nos. 2, 20, 14, 15, and 16) which are considered of high soil erosion rates, and soil loss-prone areas. It is more appropriate to integrate structural conservation measures with technology improving farming practice (i.e., rotation and contour ploughing) of rainfed cultivation to reduce soil erosion, and to enhance soil moisture in the soil profile. Such recommendations can be applied almost to all sub-basins categorized as of high, moderate, and low priority. The expansion of rainfed farming on rangeland and steep slopes (>20˚) since the 1950s, has stimulated an increase in studies showing that the presence of rock fragments on the surface were highly effective in reducing soil erosion loss. The presence of stone cobbles and fragments of different size at the soil surface at 5% to 15% coverage, caused a momentous reduction in runoff by an average of 17% and 30% respectively. The corresponding reductions in soil loss for both stone streamlets were as large as 35% and 53% respectively  . If plantation of trees is accompanied by stone coverage practice, it will decrease soil loss, and increase infiltration rate and soil moisture. Rehabilitation of vegetation cover through planting trees of drought-resistance species is highly recommended to improve the grazing potential of the marginal areas of the entire watershed.
4.3. Validation of Priority Classes: Discriminant Analysis (DA)
The validity of the three priority classes determined through morphometric analysis, and represent the 23 sub-basins was tested using DA. The intention is to test the hypothesis that there is a significant difference between the priority classes obtained earlier, and if this hypothesis is accepted, to establish a system of coordinate axes which discriminate between the recognized priority groups. With reference to the three classes, it is obvious that there is a significant difference between low priority class, moderate class, and high priority class, with a noticeable concentration of high and moderate priority classes in the eastern part of W.Wala watershed. That is; in accordance with the deteriorated rangeland, bare land, and degraded vegetation cover as clarified in this part of the watershed earlier. Soil resources have been exhausted here over the last 3000 years, with prominent shortage of conservation measures. Statistical testing using Discriminant Analysis was performed on a data matrix representing the three priority groups (i.e. 6 × 11, 8 × 11, and 9 × 11) with the associated ranking values and including the Cp scores. The F test of Wilks Lambda obtained is F ratio is 45.7 with the degrees of freedom V1 = 2 and V2 = 20. Referring to the table of percentage points of the F-distribution, with V1 = 2 and V2 = 20, it is found that at 99.9 percent of confidence, the tabulated value is 9.95, which is significantly exceeded by the computed F ratio value (45.7). Consequently, there is a high significant difference between each of the priority classes (high, moderate, and low), and the three priority groups are highly separated and distinct. Moreover, 100 percent of the difference between the three priority groups is attributed to discriminant function 1 (81.7 percent) and discriminant function 2 (18.3 percent). It was also revealed that discriminant function 1is positively correlated with Cp, Tr, Dd, Lo, Fs, and Cc parameters. By contrast discriminant function 2 is positively correlated with Bs, and Re parameters. The scores of each sub-basin of the three priority groups on the discriminant functions 1 and 2 were plotted in Figure 9. The plot displays highly separated and distinct priority groups. With reference to the present findings, it can be concluded that prioritization based on morphometric analysis is proven to be statistically valid, consistent and reliable, and of high capacity using the GIS platform. The potential of morphometric analysis method as elaborated by the pioneers, i.e.      is proven to be capable for prioritization research.
Soil erosion by water has seriously threatened rainfed farming over most of the highlands of W.Wala catchment. Further, soil erosion has also increased sediment supply to the W.Wala dam especially during exceptionally heavy rainstorms.
Figure 9. Results of discriminant functions 1 and 2 grouping: the three priority classes are completely separated and distinct.
Consequently, all sub-watersheds ranked under high and moderate priority in the present study, should be prioritized for soil conservation measures in order to maintain future farming sustainability. Integrated GIS and remote sensing, and morphometric analysis of linear and shape parameters, were employed to prioritize the 23 sub-basins, and then to prescribe efficient conservation measures especially for sub-basins subjected to rainfed farming.
The results of prioritization based on Cp values, indicate that sub-watershed no. 10 has been ranked 1 with the lowest compound parameter (Cp) at 3.2; while sub-basin no. 6 is ranked as the second with compound parameter at 3.5, and sub-basin no. 21 is ranked third, and all of them are with high priority (Figure 7). By contrast, sub-basins nos. 12, 5, and 15 are ranked as fourth, fifth and sixth with compound parameters at 4.0, 4.2, and 4.3 respectively and with moderate priority. All sub-basins with compound parameters equal or greater than 5 are ranked as low priority (Figure 7). All these sub-watersheds experienced historically severe soil erosion rates, and deterioration of natural vegetation. Information on land use/land cover, slope and soil types was utilized to help in recommending proper soil and water conservation structures. The suggested soil conservation measures take into account the existing and use/land cover, i.e., rainfed cultivation, irrigated farming, soils, rangeland, and slope categories (0 - 25) connecting to each sub-basin. However, soil and water conservation measures were recommended in accordance with the priority ascribed so as to minimize the adverse effect on environmental resources. Expected benefits of improving and expanding soil conservation measures over sub-watersheds ranked as high and moderate priority (where rainfed farming is dominant) are illustrated in the following: control of soil erosion over sub-watersheds experienced high soil erosion loss to preserve the soil from future erosion; there was a reduction in sediment yield production to control sedimentation in W.Wala dam; and reduced peak flows across these sub-watersheds and the entire W.Wala catchment. The utilization of remote sensing and techniques verify the efficiency of this technology in prioritization of W.Wala sub-basins, and prove the capacity of morphometric parameters in prioritization analysis for the highland watersheds within a GIS environment. Furthermore, the limitations of morphometric analysis method were examined, and the validity of prioritization was tested statistically using Discriminant Analysis. The results revealed that the three priority classes are significantly different from each other. Thus prioritization based on morphometric analysis method is accepted. The present results are intended to help decision-makers to recognize priority sub-basins which need immediate adaptation of appropriate conservation measures, and land management practices. In conclusion, the findings of present research indicate that prioritization based on morphometric analysis is proven to be statistically valid, consistent and reliable, and of high capacity using the GIS platform. Accordingly, the morphometric analysis method is highly recommended for prioritization research.