Prioritizing sub-watersheds for soil erosion using geospatial techniques based on morphometric and hypsometric analysis: a case study of the Indian Wyra River basin

The hydrological availability and scarcity of water can be affected by geomorphological processes occurring within a watershed. Hence, it is crucial to perform a quantitative evaluation of the watershed’s geometry to determine the impact of such processes on its hydrology. Geographic information systems (GIS) and remote sensing (RS) techniques have become increasingly significant because they enable decision-makers and strategists to make accurate and efficient decisions. To prioritize sub-watersheds within the Wyra watershed, this research employs two methods: morphometric analysis and hypsometric analysis. The watershed was divided into eleven sub-watersheds (SWs). The prioritization of sub-watersheds in the Wyra watershed involved assessing several morphometric parameters, such as relief, linear, and areal features, for each sub-watershed. Furthermore, the importance of the sub-watersheds was determined by computing hypsometric integral (HI) values using the elevation–relief ratio method. The final prioritization of sub-watersheds based on morphometric analysis was determined through the integration of principal component analysis (PCA) and weighted sum approach (WSA). SW2 and SW9 have had higher priorities using morphometric analysis, whereas SW6, SW7, and SW10 have obtained higher priorities using hypsometric analysis. SW4 is the most common SW that shares the same priority. The most vulnerable sub-watersheds are those with the highest priority, and therefore, programmes for soil and water conservation should pay more attention to them. The conclusions of the study may prove useful to various stakeholders involved in initiatives related to watershed development and management.


Introduction
Land and water are two of the most valuable and necessary resources for life and various development endeavours (Nookaratnam et al. 2005;Mukta et al. 2022;Kudnar and Rajasekhar 2020;Rendana et al. 2023). Population growth has accelerated over time, leading in a shortage of both land and water assets. Rapid industrialization is also a time need, necessitating infrastructure; this, in turn, generates a feedback mechanism, putting more strain on precious land and water assets. It is appropriate to use the watershed technique to investigate different processes occurring at the surface of the ground, which is an area of the ground where the primary discharge is transferred to a single exit (Pande and Moharir 2017;Pande et al. 2020). Watersheds, or hydrological units, are regarded to be more efficient and appropriate for conducting necessary surveys and investigations, as well as planning and implementing various improvement initiatives including water and soil conservation, and assuring their long-term viability (Singha et al. 2022;Poongodi and Venkateswaran 2018;Shekar and Mathew 2022a;Krishnan et al. 2017). Therefore, watershed management should be given special attention to address water-related issues Gautam et al. 2023).
A watershed is an area of land where all the water that falls within its boundaries drains or flows downhill into a particular body of water, such as a river, lake, or ocean. Land degradation can have a significant impact on watersheds by altering the natural processes of water flow and nutrient cycling. When land is degraded, such as through deforestation, overgrazing, or soil erosion, it can lead to the loss of vegetation cover and soil fertility, resulting in decreased water infiltration and increased runoff. This can cause erosion, sedimentation, and the degradation of water quality in downstream areas (Rekha et al. 2011;Bhattacharya et al. 2020). Floods are another significant impact of land degradation, particularly in areas with high levels of erosion. Erosion can reduce the capacity of soil to absorb water, leading to increased runoff and higher flood risks . Land-use practices such as conservation tillage, crop rotation, and cover crops, reforestation, and best management practices can mitigate the impacts of land use on soil erosion (Obeidat et al. 2019). Watershed management plans can be developed to guide land-use decisions and prioritize conservation efforts (Choudhari et al. 2018;Pande et al. 2023). The characteristics of a watershed have an effect on hydrological cycle within it, which can be examined using morphometric analysis (Pande et al. 2021a, b;Awawdeh et al. 2015).
Morphology is the study of the earth's surface using mathematics to describe its topographic reliefs Obi et al. 2002;Clarke 1966;Agarwal 1998). A detailed representation of the structure of a watershed and its stream channel needs the estimation of parameters of the channel system (Strahler 1964). Distinct academics from around the world have conducted morphometric analyses of various river basins on various continents. According to several morphometric researches, watershed morphologies reflect varied geological and geomorphological processes over time (Horton 1945;Strahler 1957).
GIS is a powerful tool used for computerized mapping and geographical analysis, which has advanced capabilities (Awawdeh et al. 2014). One such capability is the use of GIS-based evaluations with satellite image data, such as the Shuttle Radar Topographic Mission (SRTM), that allows for quick and accurate analysis of hydrological systems (Grohmann 2004). GIS techniques have been used to assess different morphometric features of stream watersheds and watersheds because they offer a comfortable workspace and a useful tool for manipulating and analysing satellite information, especially for future information collection and recognition for a better awareness. The development of RS and GIS has enabled more accurate and affordable morphometric study of natural drains (Shekar and Mathew 2022c;Grohmann et al. 2007;Abdo 2020;Smith and Sandwell 2003;Shelar et al. 2022). Prioritization of watersheds based on morphometric features has also been done, as well as assistance in detecting soil erosion areas (Mishra et al. 2011;Sharma and Mahajan 2020;Wakode et al. 2011). Many researchers studied watershed prioritization based on morphometric parameters Mathew 2023a, 2023b;Kushwaha et al. 2022;Javed and Khanday 2011;Magalhaes et al. 2022;Abdo et al. 2023;Redvan and Mustafa 2021;Gupta et al. 2020;Javed et al. 2009;Esin and Akgul 2021;Mathew et al. 2022;Bogale 2021;Singh et al. 2021;Rais and Javed 2014;Sreedevi et al. 2005;Mathew and Shekar 2023;Jasmin and Mallikarjuna 2013;Lopez-Perez and Fernandez-Reynoso 2021;Sutradhar and Mondal 2023;Moharir et al. 2021).
The distribution of ground surface cross-sectional areas in respect to elevations is the subject of the hypsometric analysis (Strahler 1952). At different stages of erosion, it is utilized to define erosional landforms (Schumm 1956). Langbein (1947) was the first to establish the concept of hypsometry, which aided in the generation of parameters such as the hypsometric curve and hypsometric integral. HC is a relationship between the amount of soil mass in a watershed and the amount of erosion in relation to the watershed's remaining mass (Hurtrez et al. 1999;Ritter et al. 2002). The structure of HC for various drainage watersheds under similar hydrologic conditions can be compared to understand watershed in the past, soil displacement. As a result, the shape of HC describes the temporal variations in the original watershed's slope. According to the shape of the hypsometric curves, Strahler (1952) classified watersheds as young, peneplain, mature. The hypsometric integral can also determine the erosion cycle (Strahler 1952). HI is equivalent to the elevation-relief ratio (E) established by Pike and Wilson when the hypsometric curve is integrated (1971). The numbers from the hypsometric integral cycle erosion of soil show that the old catchment is completely stabilized, in equilibrium, and that it is at risk of soil erosion, while it is in equilibrium. The geologic stages of watershed development are characterized by HI, a geomorphological characteristic. It is significant in determining the erosion state of a watershed (Sharma et al. 2018;Shekar and Mathew 2022b).
In reality, the importance of parameters may vary across different sub-watersheds depending on their specific characteristics. Therefore, the study used two methodsmorphometric and hypsometric analysis-to prioritize sub-watersheds for soil erosion management. The study used several morphometric parameters, to assess the vulnerability of each sub-watershed to soil erosion. On the other hand, hypsometric analysis is based on the elevation data of the sub-watershed and examines the relationship between the area and the elevation range. This analysis helps in understanding the morphological characteristics of the sub-watershed and its hydrological behaviour. In the study, hypsometric analysis was used to determine the degree of erosion susceptibility of the sub-watersheds. Using both methods helps to identify sub-watersheds that have different characteristics and prioritize them based on their specific vulnerabilities to soil erosion. Moreover, these two methods provide a more comprehensive approach to prioritize sub-watersheds for soil erosion management. The objective of the present study is to prioritize sub-watersheds by conducting morphometric and hypsometric analyses of each sub-watershed for soil erosion. The prioritization of sub-watersheds based on morphometric analysis was carried out through the integration of principal component analysis (PCA) and weighted sum approach (WSA). Furthermore, the study employs morphometric and hypsometric analyses to determine the sub-watersheds that have the priority in common.

Study area
The Wyra watershed includes Telangana state and Andhra Pradesh. The Wyra River is shown in Fig. 1. The Wyra watershed is situated between the latitudes of 16° 40′ 00′′ and 17° 35′ 00′′ north and the longitudes of 80°05′00′′ and 80°55′00′′ east. The outlet of the Wyra watershed is 16° 43′ 48′′ latitudes and 80°19′37′′ longitudes, respectively. It has an overall area of 3403 Km 2 . The Wyra watershed region experiences a semi-arid climate, characterized by hot summers and mild winters. The region receives most of its rainfall during the monsoon season, which lasts from June to September. Based on Sentinel-2 imagery from 2021, Fig. 1 Wyra River basin's geographical location the land-use/land-cover (LULC) data can be classified into the following categories and respective percentages: water (2.38%), trees (11.51%), flooded vegetation (0.01%), crops (75.51%), built area (4.61%), bare ground (0.15%), and rangeland (5.83%) (Fig. 2). This information was obtained from the following source: https:// livin gatlas. arcgis. com/ landc over/ (Karra et al. 2021). Land use in the basin is dominated by agriculture, with the majority of the land under cultivation. The main crops grown in the basin include paddy, cotton, and pulses. According to the World Geologic Maps of the United States Geological Survey (USGS), the study area has two types of significant rocks. Quaternary sediments and undivided Precambrian are the geological age of two rocks (https:// certm apper. cr. usgs. gov/ data/ apps/ worldmaps/). The Wyra watershed is 33 to 792 m above sea level, according to the digital elevation model (DEM) from SRTM (https:// earth explo rer. usgs. gov/).

Morphometric analysis
The current study's methodology includes employing ArcGIS 10.4.1 software to perform automatic extraction procedures for analysing the features of the Wyra watershed. Processing over DEM was calculated to determine the morphometric analysis, as shown in Fig. 3. Sub-watersheds are categorized (SW 1 to SW 11). The relationship between morphometric parameters and soil erosion can be direct or inverse, depending on the specific parameter (Pande et al. 2021a, b;Obeidat et al. 2023). For example, linear parameters such as drainage density, bifurcation ratio, stream length, stream length ratio, stream number, stream order, drainage intensity, mean bifurcation ratio, length of overland flow, mean stream length ratio, stream frequency, constant of channel maintenance, drainage texture, infiltration number, and rho coefficient are directly correlated with soil erosion, as they increase the potential for soil detachment and transport. Similarly, relief parameters such as relief, relative relief, maximum elevation, ruggedness number, and minimum elevation can also have a direct correlation with soil erosion, as they influence the amount of water runoff and sediment deposition. On the other hand, shape parameters such as form factor, area of watershed, circulatory ratio, watershed length, compactness coefficient, elongation ratio, perimeter of watershed, lemniscate ratio, and shape index are typically inversely correlated with soil erosion, as they reflect the compactness and irregularity of the landscape. Table 1 provides a list of the many empirical approaches that were employed to identify these characteristics. Table 2 shows the estimated and reported Wyra River basin linear parameters (SW 1 to SW 11). In order of importance, SW with the high preliminary rank (PR) of relief and linear values was ranked first. In order of importance, SW with the low preliminary rank of shape values was ranked first. In this study, principal component analysis (PCA) was utilized to identify the significant parameters for morphometric analysis. As a multivariate statistical technique, PCA is employed to reduce the dimensionality of parameters. By converting the original data, PCA generates two or more principal components. The Kaiser criterion and varimax rotation of factor loading were used to select principal components with eigenvalues greater than 1 (Kaiser 1958). In order to enhance the correlation for defining the most significant parameters, a factor loading rotation was executed. Next, a weighted sum approach (WSA) was applied to the most significant parameters obtained from PCA. The final priority ranking and categorization were determined based on the compound values, which were calculated by multiplying the ranks from morphometric analysis with their corresponding weights obtained through cross-correlation analysis of these parameters. The resulting compound factor was then used for the final prioritization of sub-watersheds. The determination of the weighted value of significant parameters (Wsp) was achieved by means of cross-correlation analysis, as expressed by Aher et al. (2014).
where CV = compound value; PRsp = preliminary ranking of significant parameter; and Wsp = weight of significant parameter. For all sub-watersheds, the priority rank was Average of bifurcation ratio of all orders Strahler (1957) Average of the steam length ratio of all orders assigned based on the lowest CV being given priority rank 1, the second lowest being assigned priority rank 2, and so forth. Following this step, the sub-watersheds were grouped into three categories based on their CV values.

Hypsometric analysis
Hypsometric analysis is a method used to study the topographic relief of a landscape. It involves the analysis of the hypsometric curve, which is a plot of the cumulative area of a region at different elevations. In other words, it shows the proportion of the landscape at different elevations. Hypsometric analysis with a DEM is a powerful tool for understanding the topography and landscape evolution of a region. Hypsometric analysis is widely used in earth science research, particularly in the fields of geomorphology and hydrology (Pande et al. 2021a, b). To conduct hypsometric analysis with a DEM, the elevations of the region are first extracted from the DEM and sorted into elevation intervals. The area of each interval is then calculated, and the cumulative area and elevation for each interval are plotted on a graph to create a hypsometric curve. SRTM-DEM and ArcGIS 10.4.1 were used to create a hypsometric curve for the Wyra watershed. Using attribute feature classes that take these variables into account, HC for the watersheds under study was plotted. The elevation-relief ratio approach is used to compute the HI values in this study, as shown in Table 1. After obtaining the HI values, divide them into 3 equal intervals to assign ranking, as the SWs were divided into 3 groups. The highest priority is given to the maximum interval values, the next interval values are given a medium priority, and the minimum interval values are given a low priority.

Results and discussion
Morphological parameters are measurements that describe the shape and form of a landscape. In the context of the Wyra watershed, these parameters were computed using digital elevation models and GIS tools to analyse the topography of the area. Table 3 shows the computed morphological parameters for the Wyra watershed, which are categorized into three aspects: linear, relief, and areal.

Morphometric analysis
The Wyra drainage watershed's morphometric study offers a quantitative description of the watershed geometry, which makes it easier to understand the geomorphological characteristics of the watershed and its responses to various hydrological processes (Chatterjee and Tantuley 2006). For analysis and discussion, the morphometry of the Wyra watershed is divided into three categories.

Linear aspects
Stream order (U) The assignment of stream orders based on a hierarchical classification of streams is the first step in drainage watershed study. The approach suggested by Strahler (1964) was used to rate streams in the current study. The watershed's drainage system is dendritic to sub- Stream length (L u ) Horton's (1945) law was used to compute the length of the stream. One of a region's key hydrological features is a stream's length because it offers information on the characteristics of surface runoff. SW3 (568 km) has a long L u , while SW2 (191 km) has a short L u in this study. Bifurcation ratio (R b ) It is the proportion of the number of streams in one order to those in the order above it. SW10 (19.86) has a high R b in this study, while SW7 (14.33) has a low R b .

Stream length ratio (R l )
It is the one order's segment divided by the mean stream length of the lower order segment. In this study, SW9 (3.76) has a high R l, whereas SW2 (1.76) has a low R l .  Horton (1945) identified two fundamental principles linking the number of distinct orders in a stream catchment and the length of a stream. Figure 5 demonstrates a strong correlation between stream order and stream number, with coefficients of determination in the range of 0.999 for SW2 and 0.966 for SW9. Similarly, Fig. 6 shows a significant correlation between stream order and stream length, with coefficients of determination in the range of 0.999 for SW2 and 0.774 for SW9.
Stream frequency (F s ) The total number of orders' streams in a certain area. SW9 (1.08) has a high F s in this study, while SW6 (0.37) has a low F s .

Drainage density (D d )
It is the proportion of the total length of all streams in a catchment to its entire area. It is a critical characteristic that is directly related to runoff speed, which is followed by precipitation. In this study, SW3 (1.17) has a high D d, while SW6 (0.86) has a low D d .
Drainage texture (D t ) One of the essential terms in geomorphology is drainage texture, which refers to the relative spacing of drainage lines. SW3 (2.59) has a high drainage texture in this study, while SW6 (0.74) has a low drainage texture.
Length of the overland flow (L o ) It is equal to half of D d (Horton 1945). SW6 (0.58) has a long overland flow, while SW3 (0.43) has a short overland flow in this study. (Faniran 1968). SW9 (0.93) has high drainage intensity in this study, while SW6 (4.50) has low drainage intensity.

Drainage intensity (D i ) It is the proportion of F s to D d
Rho coefficient ( ) It is the ratio of stream length to R b . SW11 has a high rho coefficient in this study, while SW9 has a low rho coefficient.
Infiltration number (I f ) Faniran (1968) defines it as the result of the interaction between D d and F s . SW9 (1.25) has a high I f in this study, while SW6 (0.32) has a low I f .

Constant of channel maintenance (C cm )
It was initially proposed by Schumm in 1956, defined it as the inverse of D d . SW6 (1.16) has a high C cm in this study, while SW3 (0.85) has a low C cm .

Relief aspects
Relief (B h ) The difference in elevation between the watershed's maximum and minimum elevations is known as a relief. In this study, SW3 (0.72) has a high relief, while SW10 (0.12) has a low relief.

Relative relief (R hp )
From the highest point on the catchment's border to the stream's mouth, the greatest amount of watershed relief was attained. In this study, SW9 (0.50) has a high relative relief, while SW10 (0.09) has a low relative relief.

Ruggedness ratio (R n )
The product of the watershed relief and D d is the ruggedness number. In this study, SW3 (0.85) has a high ruggedness ratio, while SW10 (0.59) has a low ruggedness ratio.

Area of watershed (A)
The amount of runoff a catchment generates is directly influenced by the watershed region. A total area of 3403 km 2 is covered by the watershed. SW10 (497.30 Km 2 ) has a large watershed area, while SW2 (174.19 Km 2 ) has a small area in this study, as shown in Fig. 7.

Perimeter of a watershed (P)
The watershed perimeter is the outer limit of the catchment that defines its area. SW4 (198.80 km) has a large watershed perimeter, while SW2 (172.07 km) has a small watershed perimeter in this study, as shown in Fig. 8.

Watershed length (L b )
Among the fundamental dimensions of the main drainage channel, the catchment length is the most crucial. SW10 (44.63 km) has a large watershed length, while SW2 (24.59 km) has a small watershed length in this study, as shown in Fig. 9.
Circulatory ratio (R c ) It is the ratio of the area of the watershed to the surface area of a circle whose circumference equals the area of the catchment. In this study, SW2 (0.27) has a high R c , while SW4 (0.08) has a low R c .
Elongation ratio (R e ) According to Schumm (1965), it is the ratio of the diameter of a circle with the same area of the watershed to the maximum watershed length. In this current research, SW2 (0.61) has a high R e , while SW10 (0.56) has a low R e .
Form factor (F f ) According to Horton (1932), it is the ratio of catchment area to the square of catchment length. In this current research, SW10 (0.29) has a high F f , while SW2 (0.25) has a low F f . Chorely et al. (1957) uses lemniscate's value to calculate the catchment's gradient. In this study, SW10 (1.0) has a high K, while SW2 (0.87) has a low K.

Lemniscate ratio (K)
Compactness coefficient (C c ) Horton (1945) defined the compactness coefficient as the ratio of the catchment's perimeter to that of a comparable circular region. In this study, SW4 (3.44) has a high C c , while SW2 (1.94) has a low C c .

Shape index (S b )
The reciprocal of F f is the shape index. Horton was the one who first recommended it (1932). In this study, SW10 (4.01) has a high shape index, while SW2 (3.47) has a low shape index, as shown in Fig. 10.

Sub-watersheds prioritization of morphometric analysis based on PCA-WSA
The nineteen morphometric features representing linear, shape, and relief aspects were utilized to select sub-watersheds for soil conservation. Since soil erosion is directly correlated with the relief and linear features, rank one is given a higher value. However, because soil erosion is indirectly Fig. 9 A sub-watershed's watershed length correlated with the areal features, rank one is given a lower value (Nookaratnam et al. 2005). The higher value gets a 1 preliminary ranking for relief and linear features, and so forth. A preliminary ranking of 1 was assigned to the shape feature with the lowest value, and so forth ( Table 4).
The purpose of conducting PCA on all parameters is to assess their correlations, identify principal components, and reduce the parameter dimensionality to highlight the most significant ones. Table 5 displays the correlation matrix of all parameters. To analyse the correlation among geomorphic parameters, a correlation matrix is generated using SPSS 14.0 software. After analysing the correlation matrix of the 19 geomorphic parameters in Wyra watershed, it is evident that strong correlations (with a correlation coefficient exceeding 0.9) are present between D i and F s , I f and F s , L o and D d , C cm and D d , I f and L o , C cm and L o , I f and D i , R n and B h , C c and R c , F f and R e , K and R e , S b and R e , K and F f , S b and F f , S b and K. The good correlation (correlation coefficient between 0.75 and 0.9) is between and R lm , D d and F s , L o and F s , C cm and F s , I f and D d , R hp and B h , R n and R hp . Moderately correlated parameters (correlation coefficient more than 0.6) are P and R bm , D t and F s , D t and D d , D i and D d , B h and D d , R hp and D d , R n and D d , L o and D t , D i and D t , I f and D t , C cm and D t , D i and L o , R hp and L o , R n and L o , C cm and D i , R hp and C cm , R n and C cm . Correlation between parameters indicates that there may be shared information across multiple parameters. However, it is often difficult to group these parameters into meaningful components and assign physical interpretations. To address this, one practical approach is to use principal component analysis (PCA) on the correlation matrix to reduce the parameter dimension (Meshram and Sharma 2017). Therefore, in the next step, PCA has been applied to the correlation matrix.
To obtain the first factor loading matrix, the principal component analysis method was utilized, followed by orthogonal transformation to obtain the rotated loading matrix. Using the correlation matrix of 19 geomorphic parameters, the first unrotated factor loading matrix was derived. Table 6 shows that the first five components, with eigenvalues greater than 1, account for approximately 94.41% of the total variance in the Wyra watershed. Table 7 indicates that the first component exhibits a strong correlation (above 0.90) with D d , L o , I f , and C cm . It also shows a good correlation with F s and R hp , and a moderate correlation with D t , D i , B h , and R n . The second principal component exhibited a strong correlation with R e , F f , K, and S b , and a moderate correlation with R lm . The third principal component showed a strong correlation with R bm and a moderate correlation with R c and C c . The fourth and fifth principal components exhibited no correlation with any of the ranges.

Table 5
The correlation matrix among linear, areal, and relief parameters Parameters Based on the first factor loading matrix, it can be observed that is not correlated with any of the components. While some parameters show a high correlation with certain components, others exhibit a moderate correlation, and some parameters do not correlate with any component at all. As a result, it is difficult to determine the most significant parameters for each principal component. To establish better correlations and identify significant parameters, it is necessary to rotate the first factor loading matrix. The rotated factor loading matrix is presented in Table 8, which reveals that the first, second, third, fourth, and fifth principal components exhibited strong correlations with F s , K, B h , C c , and , respectively. These parameters are also considered significant and are utilized in WSA and sub-watershed prioritization.
The compound value (CV), which was determined by incorporating the preliminary rank (PR) and the weight of significant parameters (F s , K, B h , C c , and ), was utilized for the ultimate sub-watershed prioritization. The weightage of these critical parameters was determined by cross-correlation analysis among the five parameters listed in Table 9. By applying the weighted sum of essential parameters, the compound value (CV) was derived using Eq. (1). Subsequently, the sub-watershed with the lowest CV (3.82), i.e. SW2, was accorded the highest priority. Conversely, the sub-watershed with the highest CV (8.88), i.e. SW6, was deemed to be of the lowest priority.
The subdivision of the compound value into three categories for soil erosion is based on the level of erosion risk associated with each category. The three categories are known as the low, moderate, and high categories. The basis for this subdivision is the fact that different levels of soil erosion risk require different levels of intervention to prevent or mitigate erosion. If the compound value is in the low-risk category, it indicates that the erosion risk is minimal, and erosion control measures may not be necessary. In this case, the focus may be on maintaining current land-use practices and monitoring erosion levels to ensure they do not increase. If the compound value is in the moderate-risk category, it indicates that erosion control measures may be necessary to prevent erosion from causing significant damage. In this case, erosion control practices, such as vegetation management,    Figure 11 shows the Wyra watershed's final priority map of sub-watersheds based on the morphometric study. The degree of erosion in a sub-watershed is directly proportional to its priority level, with higher priority indicating a greater extent of erosion. Such sub-watersheds are therefore considered as potential areas where soil conservation measures should be enforced. The analysis of morphometric parameters revealed that sub-watersheds SW2 and SW9 exhibit a particularly high susceptibility to soil erosion in the study area. The underlying parameters that make SW2 and SW9 a high-priority sub-watershed could vary based on the specific morphometric parameters used in the analysis. For instance, if the sub-watershed has a high stream frequency, it could also increase the risk of flooding during heavy rainfall events. Additionally, if sub-watershed has a high drainage density, which indicates a large number of streams and channels in the sub-watershed, it may be more susceptible to erosion and sedimentation. On the other hand, if the sub-watershed has a complex drainage pattern or a high relief, it may indicate a greater potential for soil erosion and landslides. The high-risk sub-watersheds must receive immediate attention, followed by the medium-risk sub-watersheds, and so on, until sufficient time and resources are available for the remaining sub-watersheds.

Relation among morphometric features and hypsometric curve
The relief ratio and watershed volume represented by the hypsometric curve (HC), according to Vivoni et al. (2008), are useful in evaluating runoff and other hydrological processes. The hypsometric curves reveal not only the watershed's  Page 21 of 26 160 Fig. 12 The hypsometric integral values of each sub-watershed erosion status but also the tectonic, climatic, and lithological elements that influence it (Sarp et al. 2011). The drainage network and watershed geometry have a big impact on hypsometry. The aspect ratio decreases, the stream system becomes more branched, and the bifurcation ratio increases. The toe height will then be of increasing elevation at the downstream part of the watershed for a low aspect ratio (Roy 2002).

Prioritization of sub-watersheds based on hypsometric analysis
The hypsometric integral (HI) was obtained by using elevation-relief ratio (E) approach. HI in this study ranges from 0.1285 to 0.4018, as shown in Fig. 12. The ranking values for all eleven SWs were based on hypsometric integral values. High (0.311 to 0.402), medium (0.311 to 0.220), and low (0.220 to 0.128) were used to categorize the subwatersheds (Farhan et al. 2016). As stated in Table 11, highpriority sub-watersheds are SW6, SW7, and SW10, there are no sub-watersheds of a medium priority, and low-priority sub-watersheds are SW1, SW4, SW2, SW9, SW5, SW11, and SW8. The hypsometric analysis's final priority map is displayed in Fig. 13.

Common sub-watersheds
The most common sub-watersheds were determined based on morphometric and hypsometric analysis. The most common watersheds are SW4 (low priority), as shown in Table 12. SW4 is identified as a low-priority sub-watershed in both the morphometric analysis-based PCA-WSA and hypsometric analysis-based prioritization, and it suggests that this sub-watershed is less susceptible to land degradation compared to other sub-watersheds in the Wyra River basin. Therefore, the implementation of management practices in SW4 may not be an immediate priority. However, it is still important to monitor the condition of SW4 to ensure that it remains stable and does not become more vulnerable to land degradation in the future.

Conclusion
The most crucial aspect of organizing and implementing watershed improvement and management programs is prioritizing the watershed. To facilitate efficient planning for watershed management, this study employed geospatial techniques to conduct morphometric and hypsometric analyses of the sub-watersheds of the Wyra River basin. The prioritization of sub-watersheds based on morphometric analysis using PCA-WSA reveals that SW2 and SW9 are considered high-priority areas. However, when hypsometric analysis is used for prioritization, the high-priority areas are SW6, SW7, and SW10. The sub-watershed with the most common occurrence in both prioritization methods is SW4. To mitigate soil erosion, conservation measures such as artificial recharge structures (e.g. check dams and percolation tanks) can be implemented in the area. These measures have the potential to reduce erosion and can be effective tools for ensuring the long-term sustainability of the watershed. Moreover, the insights obtained from this study can be valuable for decision-makers in the Wyra watershed to implement effective management practices aimed at mitigating and preventing land degradation.
In future studies, it is recommended to consider incorporating social and economic factors could contribute to a more comprehensive understanding of the sub-watersheds and their priority for management interventions. Additionally, evaluating the effectiveness of management practices implemented based on subwatershed prioritization could help refine and improve the process for future watershed management initiatives.

Data availability
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Conflict of interest
On behalf of all authors, the corresponding author states that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Ethical approval Not applicable.

Consent for publication Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.