Next Article in Journal
Experimental Research on Deformation Characteristics of Using Silty Clay Modified by Oil Shale Ash and Fly Ash as the Subgrade Material after Freeze-Thaw Cycles
Next Article in Special Issue
Recognizing Crucial Aquatic Factors Influencing Greenhouse Gas Emissions in the Eutrophication Zone of Taihu Lake, China
Previous Article in Journal
The Concept of Financial Sustainability Measurement: A Case of Food Companies from Northern Europe
Previous Article in Special Issue
Water Quality Sustainability Evaluation under Uncertainty: A Multi-Scenario Analysis Based on Bayesian Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multivariate Geomorphometric Approach to Prioritize Erosion-Prone Watersheds

by
Jesús A. Prieto-Amparán
1,
Alfredo Pinedo-Alvarez
1,
Griselda Vázquez-Quintero
2,
María C. Valles-Aragón
2,
Argelia E. Rascón-Ramos
1,
Martin Martinez-Salvador
1 and
Federico Villarreal-Guerrero
1,*
1
Facultad de Zootecnia y Ecología, Universidad Autónoma de Chihuahua, Chihuahua 31453, Mexico
2
Facultad de Ciencias Agrotecnológicas, Universidad Autónoma de Chihuahua, Chihuahua 31350, Mexico
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(18), 5140; https://doi.org/10.3390/su11185140
Submission received: 15 August 2019 / Revised: 12 September 2019 / Accepted: 15 September 2019 / Published: 19 September 2019
(This article belongs to the Special Issue Sustainability in the Development of Water Systems Management)

Abstract

:
Soil erosion is considered one of the main degradation processes in ecosystems located in developing countries. In northern Mexico, one of the most important hydrological regions is the Conchos River Basin (CRB) due to its utilization as a runoff source. However, the CRB is subjected to significant erosion processes due to natural and anthropogenic causes. Thus, classifying the CRB’s watersheds based on their erosion susceptibility is of great importance. This study classified and then prioritized the 31 watersheds composing the CRB. For that, multivariate techniques such as principal component analysis (PCA), group analysis (GA), and the ranking methodology known as compound parameter (Cp) were used. After a correlation analysis, the values of 26 from 33 geomorphometric parameters estimated from each watershed served for the evaluation. The PCA defined linear-type parameters as the main source of variability among the watersheds. The GA and the Cp were effective for grouping the watersheds in five groups, and provided the information for the spatial analysis. The GA methodology best classified the watersheds based on the variance of their parameters. The group with the highest prioritization and erosion susceptibility included watersheds RH24Lf, RH24Lb, RH24Nc, and RH24Jb. These watersheds are potential candidates for the implementation of soil conservation practices.

1. Introduction

Soil erosion is considered one of the most important degradation processes in the world [1,2]. The soil resource is limited and its wide use is of utmost importance; it sustains biogeochemical processes and is the habitat for a great diversity of microorganisms [3]. Sustained soil development, conservation, and restoration is one of today’s main challenges for humankind.
Hydric erosive processes affect the fertile soil layer, which is a key factor in the primary production of ecosystems [4]. The production of goods and satisfiers for the population such as wood, food, fiber, fodder, water, and recreational areas, among others, in addition to industrial expansion and the need for infrastructure facilities, have increased land-use/land-cover changes, increasing the pressure over the soil [5]. This has caused experts to pay more attention to the growing trend of soil erosion and the importance of water and soil conservation for achieving sustainable development.
Integrated watershed management is an alternative for soil management [6,7,8]. Watersheds are one of the spatial units that are used for the planning and management of soil resources [9]. Management implies the characterization of the ecosystems inside the watershed and the understanding of the relationships between uplands, lowlands, land use/land cover, geomorphic processes, slope, and soil [10,11]. In watershed management, erosion control is one of the main components [3]. Thus, the hydrological planning and monitoring of a watershed is important for the development of environmental policies [12]. In this sense, the quantification of the watersheds’ characteristics is fundamental to understanding their dynamics and degradation levels. This knowledge serves to define and implement strategies to prevent soil erosive processes and promote the conservation and restoration of watersheds [13].
Morphometry is used in the analysis of the watershed configuration [14]. Such methodology was developed by Horton [15] and then modified by Strahler [16], and provides information on the behavior of the basin [17]. It is an important tool for identifying and prioritizing eroded watersheds [18].
Nevertheless, monitoring soil erosion in situ is costly and time consuming in large watersheds. Thus, the analysis of geomorphometry is often carried out based on geographic information systems (GIS) [19,20,21,22,23]. On a spatial scale, geomorphometric parameters, e.g., the Gravelius compactness coefficient [24] and elongation ratio [25], among others, are important to know the hydrological configuration of watersheds. The relationships among these parameters are useful for developing hydrological models, which allow prioritizing watersheds based on their condition, such as erosion susceptibility. To determine the aforementioned relationships, statistical methods, such as multivariate techniques, have been widely used worldwide [26,27]. For instance, Gavit et al. [28] characterized 13 geomorphometric parameters in 11 watersheds located in the Godavari river in Maharashta, India. Youssef et al. [29] estimated the erosion risk by using remote sensing technology, GIS, and geomorphometric parameters in 11 watersheds located in Sinai, Egypt. Makwana and Tiwari [30] used seven geomorphometric parameters to characterize 19 watersheds in the region of Gujarat, India. Sharma et al. [31] applied the multivariate technique of principal component analysis (PCA) to 13 geomorphometric parameters from eight watersheds located in the Madhya district, India. Meshram and Sharma [32] and Farhan et al. [33] used PCA to analyze the geomorphometric parameters of a group of watersheds located in the Shakkar Catchment River in India and Jordan.
A large number of shape, relief, and hydrological parameters are associated with watershed geomorphometry [34,35]. The statistical techniques of PCA and group analysis (GA) [36,37,38,39,40,41], as well as multivariate analysis of variance (MANOVA) [41] and the ranking methodology known as compound parameter (Cp) [42] have been widely used in recent years for the analysis of environmental data from watersheds. These techniques assist with analyzing the spatial variability of watersheds, their structure, as well as the relationships existing among them.
The most important basin in the state of Chihuahua, as a runoff source [43], is the Conchos River Basin (CRB). Yet, this basin has experienced water stress conditions during the past years. According to Ordoñez [44], approximately 8000 km2 (11.8%) of the basin high lands present strong erosion problems, which could impact waterflow and water quality. In these high lands, deforestation and land-use/land-cover changes had contributed to reduce the amount of infiltrated water, impacting on groundwater availability [43]. In the low lands of the basin, agriculture consumes 90% of the water harvested in the basin. Other consumers include the industrial and domestic sectors [43,45]. In addition, the international water trade between Mexico and the U.S. from 1944 commands Mexico to deliver annually from this basin a total of 432 × 106 m3 of water to the U.S. [46]. Therefore, specific knowledge about water management and the status of the basin’s soil erosion is required to implement strategies to solve water-related problems and to promote the sustainable development of the region.
The objective of this study was: (1) to describe the behavior of the 31 watersheds located along the CRB in the state of Chihuahua, Mexico, based on the values of their geomorphometric parameters; (2) to spatially classify the 31 watersheds into groups; and (3) to prioritize the groups of basins according to their erosion susceptibility. For the grouping, multivariate techniques and the compound parameter (Cp) were used and their efficacy compared.

2. Materials and Methods

2.1. Study Area

The CRB is located in the state of Chihuahua and Durango, Mexico, and is part of the 24th Rio Bravo-Conchos Hydrological Region [47] (Figure 1). The basin has an area of 67,800 km2 [48], with an altitudinal distribution ranging from 841 m to 2348 m [49]. It presents a diversity of climates ranging from temperate in the upper, semi-arid in the middle, and arid in the lower parts of the basin [50]. The physiography of the upper basin belongs to the mountainous zone made up of temperate forests dominated by species of pines and oaks. The middle part of the Altiplano or central valleys is made up of transition zones where oaks and bushes are present. Regarding the lower part, it belongs to the arid region and is made up of shrublands and grasslands [51]. The basin has a precipitation regime from June to September, with July and August being the wettest months, and fluctuating from 200 mm to 700 mm [48].

2.2. Data

Data of the CRB was obtained from the online GIS source of CONABIO [52]. Likewise, the data of the 31 watersheds composing the CRB (Figure 2) were obtained from the Watershed Water Flows Simulator [53]. The relief and hydrology type parameters were estimated by processing the necessary data from a Digital Terrain Model (DTM), with a resolution of 15 × 15 m, downloaded from INEGI [54]. The values of the basic parameters from the watersheds were obtained by using the Hydrology tool [55] of ArcMap© 10.3 (ESRI, Redlands, CA, USA; https://wwwesri.com/en-us/home). The values of the shape, relief, and linear type parameters were calculated from the equations listed in Table 1.
Basic parameters include the area (A), perimeter (P), watershed length (Lb2), stream order (u), main channel length (Lc), all channel lengths (Lu), and contour length (Li). The area and perimeter are the most important parameters of the watersheds to understand their hydrological design and reflect the volume of water that can be discharged in a rainfall event [61].
Shape parameters include the Gravelius compactness coefficient (Cc), elongation ratio (Re), shape factor (Rf), elongation index (Ia), unit shape factor (Ru), and circularity ratio (Rc). The Gravelius compactness coefficient is the relationship between the perimeter of the watershed and the perimeter of a circle of area equal to that of the watershed [24]. Low values of the Gravelius compactness coefficient, shape factor, and elongation index indicate a more elongated watershed, while high values correspond to a less elongated watershed. Watersheds with less elongated shapes have a shorter maximum flow duration, while elongated watersheds correspond to watersheds with longer flow duration [24]. The elongation ratio is the diameter of a circle with an area equal to that of the watershed divided by the length of the watershed [32]. Those watersheds with values close to or greater than the unit correspond to flat watersheds, while values that move away from the unit are watersheds with pronounced relief [25]. The shape factor is the relationship between the area and the length of the watershed [15]. The elongation index is also a relationship between the length and width of the watershed [15]. The unit shape factor is the relationship between the length and the area of the watershed. Values less than 2 indicate that they have weak flood discharge periods, while values greater than 2 indicate that their flood discharge is strong [30]. The circularity ratio is defined as the relationship between the area and the perimeter of the watershed. If the circularity of the watershed is low, then discharge will be slow, and the susceptibility to erosion will be low [61].
Relief parameters include the mean watershed slope (J), massivity coefficient (tgα), relief ratio (Rh), relative relief (Rr), and orographic coefficient (Co). The mean watershed slope indicates the degree of the terrain roughness. As the slope increases, the watershed will be prone to erosion. The massivity coefficient represents the relationship between the mean watershed height and its surface area, which is expressed as a percentage. Small values of the massivity coefficient correspond to very mountainous watersheds, while large values correspond to watersheds from moderately mountainous to flat. The relief ratio is directly related to the slope of the currents and the watershed surface, being an indicator of hydrological processes and erosion. The relief ratio, similar to the relative relief, has a direct correlation with the watershed erosion processes [25].
Linear parameters include drainage density (Dd), mean slope of the main channel (j), mean distance (Am), sinuosity of the main channel (Scp), Kirpich concentration time (Tc), average peak flow (Qp), texture ratio (T), rivers frequency (Fu), resistance number (Rn), general flow length (Lo), and drainage intensity (Di). High drainage density values indicate a high current density, and therefore a rapid response to precipitation events [62]. The mean slope of the main channel indicates the slope of the longest channel in the watershed. The high values of this parameter indicate that sediment flow and entrainment will quickly exit the watershed [15]. The sinuosity of the main channel indicates the velocity of flow movement in the channels. The lowest values of sinuosity correspond to channels where the flow travels at greater speed, whereas the channels with the highest values of sinuosity indicate the slow movement of the flow [57]. The Kirpich concentration time is the time when a drop of water falls at the furthest point until it reaches the exit point [58]. Average peak flow is defined as the mean maximum amount of water passing through a specific section [59]. The texture ratio corresponds to the relationship between the total number of streams and the watershed perimeter. Rivers frequency represents the total number of streams of all orders per unit area [15]. The resistance number is used to measure the flood potential of rivers. It has a direct relationship with erosion, where increasing its value would represent an increment in erosion susceptibility [15].

2.3. Watershed’s Description and Classification

Prior to the watersheds’ classification, their geomorphometric parameters were correlated [6,63]. Correlation indicates when part of the information contained in a set of geomorphometric parameters is also contained in the remaining ones [32]. That served to reduce the number of parameters included in the subsequent analysis.
To describe the variability of the geomorphometric parameters, principal component analysis (PCA) was performed. The technique of PCA reduces the data dimensionality, simplifies the dataset, and makes it easier to explain through a set of new principal components (PCs) [64,65]. The first principal component is the linear combination of the original geomorphometric parameters that contributes the most to the total variance; the second principal component, not correlated to the first, contributes the most to the residual variance, and so on until the total variance is analyzed. The PCs are orthogonal variables that could be obtained by multiplying the original variables, which are correlated, with coefficients similar to the ones described in Equation (1):
z i j = a i 1 x 1 j + a i 2 x 2 j + + a i m x m j
where z represents the coefficient of the component, a represents the weight of the component, x represents the measured value of the variable, i corresponds to the component number, j represents the sample number, and m represents the total number of variables.
The PCA was applied to the values of geomorphometric parameters to calculate the correlation matrix and to obtain the PCs [66]. Both the correlation analysis and the PCA were performed in the SAS© 9.1.3 software (The SAS Institute Inc., Cary, NC, USA, https://www.sas.com/en_us/home.html). Then, the PCs were mapped for interpretation.
For the classification, the first method used for the watersheds was a hierarchical group analysis, which was based on the Ward criterion [67]. The Ward criterion was applied to the GA by using the square Euclidean distance to explore the clustering of the 31 watersheds. The definition of the watershed groups was performed based on the coefficient of determination (R2) [68]. Finally, a multivariate analysis of variance (MANOVA) was used to analyze whether significant multivariate differences exist between the groups based on the values of their geomorphometric parameters [69].
The second classification method considered in this work was the compound parameter (Cp). Previous research has employed this approach for sustainable watershed planning and management [42]. Linear and shape parameters have been commonly used for this method, whereas the relief and basic parameters were additionally included in this study. Linear geomorphometric parameters have a direct relationship with erosion susceptibility, where high values are more likely to result where high erosion probabilities are present [61,70]. Thus, for watershed classification, the highest value of linear parameters was ranked as 1, the second highest was ranked as 2, and so on. Conversely, shape parameters have an inverse relationship with erosion [61,70], and low values are related to high susceptibility to erosion and vice versa. Then, the lowest value of the shape parameters was ranked as 1, the next lowest value was ranked as 2, and so on. Regarding the relief and basic parameters, the highest value was ranked as 1, the second highest value was ranked as 2, and so on. After this procedure was completed, the ranked values from each watershed were summed and then averaged to produce the Cp of each watershed. This average represents the collective impact of all the parameters, and is calculated according to Equation (2) [42]:
C P = 1 n i = 1 n R
where Cp is the compound parameter, R is the ranked value of each parameter, and n is the number of parameters.
Based on the Cp, the highest priority was assigned to the watersheds with the lowest Cp value, the second priority was assigned to the next higher Cp value, and so on. Then, the Cp was classified into five categories or groups of erosion susceptibility. The categories were defined as: Very High (Group 5), High (Group 4), Moderate (Group 3), Low (Group 2), and Very Low priority (Group 1), similar to classifications made in previous studies [71].

2.4. Comparison of the Classification Methods

To compare the grouping methods, the following procedure was followed.
First, an ANOVA was carried out for each geomorphometric parameter (independent variable), and separated for each grouping method. The source of variation or class was considered to be the group. Such analyses served to determine possible significant differences among the groups within each grouping method. After the ANOVA analyses were completed, the grouping method that achieved the highest number of significant p-values (α < 0.05) was considered the most effective for grouping the watersheds of the CRB.
A list of abbreviations can be found in Appendix B (Table A3).

3. Results

3.1. Characterization of Conchos River Basin Watersheds

The values of the basic, shape, relief, and linear-type parameters from the 31 watersheds used in this study are shown in Table A1. Watershed RH24Ia has the smallest area and perimeter with 78.62 km2 and 50.59 km, respectively. Meanwhile, the watershed with the largest area and perimeter is RH24Lb, with 5428.29 km2 and 640.77 km, respectively. The watershed with the highest stream order is RH24Kb, while watershed RH24Lb has the longest main channel. First-order streams do not have tributaries, and their flow depends on the secondary surface contributions that converge to them [61]. The watershed with the highest number of channels of order one is RH24Ia, whereas the watershed with the lowest number of channels is RH24Jb. The watershed lengths vary from 11.17 km to 133 km. Watershed RH24Kb presents the largest elongation ratio value, indicating that it is the flattest, while watershed RH24Ib has the lowest value for this parameter, indicating that it has the steepest slope [25]. The values of sinuosity of the main channel vary from 0.05 km to 4.1 km. The watershed with the shortest Kirpich concentration time is RH24Kb, while watershed RH24Lb has the longest. The lowest average peak flow value corresponds to watershed RH24Ia, whereas watershed RH24Lb presents the highest. The texture ratio values are between 4.98–75.1, which are considered as moderate to high values; the low values correspond to watershed RH24Na, while the highest value belongs to watershed RH24Lg, so the former is not susceptible to erosion, while the latter is. The watershed with the lowest resistance value is RH24Na, while the watershed with the highest value is RH24Lg.

3.2. Correlations and Principal Component Analyses

The data matrix of the 31 watersheds and the 33 parameters, which included the basic shape, linear, and relief type parameters, were analyzed through a correlation analysis. A set of parameters showing high correlations were identified. From each pair of highly correlated parameters, only one parameter was chosen, and the rest were eliminated to reduce the data dimensionality. After the reduction, the final number of geomorphometric parameters was 26, as listed in Table A2.
The PCA was performed on the 26 parameters and showed that the first five principal components were the most important for explaining the data variance (Figure 2). The most important parameter was selected according to its contribution to the principal component, as it is shown by the values of the eigenvectors (in bold) of the correlation matrix (Table A2). The first five principal components accounted for 88.44% of the total variance of the dataset. The linear parameters (hydrology) are the ones mainly explaining the CRB behavior (Table 2).

3.3. Watershed’s Classification Based on Group Analysis

Figure 3 shows the dendrogram, resulting from the group analysis (GA) of the 31 watersheds. Five groups were identified based on their basic, relief, shape, and hydrology-type parameters, and considering a value of R2 = 0.84. Group 1, with seven watersheds, has the largest amount of low values for the shape, relief, and linear-type parameters, as can be verified in Table 3. Group 2, with eight watersheds, presents the lowest average values of drainage density, sinuosity of the main channel, and general flow length. Group 3, also with eight watersheds, showed the highest values of elongation ratio, drainage density, mean slope of the main channel, and general flow length. Group 4, with four watersheds, has the highest values in maximum and minimum height, mean slope of the watershed, mean distance, unit shape factor, resistance number, and drainage intensity. In this group, the lowest values correspond to the elongation ratio. Group 5, also with four watersheds, presents the largest amount of high values for the basic, shape, relief, and linear-type parameters. The multivariate analysis of variance (MANOVA) showed significant differences among the groups of watersheds, showing a value of Wilks’ lambda equal to 0.0025, with a value of p < 0.0001.
The geospatial distribution of the groups is shown in Figure 4. Group 1 shows a homogeneous pattern in its distribution, which is concentrated in the central part of the study area. In contrast, Group 2 shows a dispersed distribution, mainly at the edges of the CRB. Group 3 is distributed in the northern, central, and southern parts, and is represented by small clusters of two or three watersheds. Group 4 corresponds to watersheds spatially dispersed over the basin. The watersheds of these groups are isolated from the other watersheds of the same group. Group 5 shows watersheds clustered in the southern part of the basin, except for the watershed RH24Jb, which is located in the northern region. The GA showed a clustered geospatial pattern for Groups 1, 3, and 5, who share characteristics in their parameters and space.

3.4. Watershed’s Classification Based on the Compound Parameter (Cp)

Considering the 26 geomorphometric parameters selected after the correlation analysis was performed, the value of the compound parameter (Cp) was calculated for the 31 watersheds of the CRB (Figure 5a). The watershed RH24Lg received the highest priority (1), followed by the watershed RH24Le (2). The watershed with the lowest priority (31) was watershed RH24Kg. A high priority is an indicator of a high degree of erosion susceptibility in the watershed.
The resulting Cp map (Figure 5a) was reclassified into the following five categories: Very High, High, Moderate, Low, and Very Low (Figure 5b). The spatial distribution of the groups was reclassified by natural breaks [71]. The watersheds classified as Very High show a homogeneous pattern in their distribution in the southwestern part of the watershed. Meanwhile, the watershed RH24Jb is isolated in the northwestern part. The High, Moderate, and Low classes show a dispersed distribution, with at least two of their watersheds clustered in space. The Very Low class shows a homogeneous distribution in space in the southeastern part of the study area, with only one dispersed watershed (RH24Hf).

3.5. Comparison of the Classification Methods

Regarding the GA classification method, the results from the ANOVA analyses performed on 14 geomorphometric parameters, out of 26, detected significant differences among the groups defined by the method. In the case of the Cp classification method, the ANOVA analysis of only two parameters detected significant differences among the groups defined by this classification method (Table 4).

4. Discussion

The prioritization of watersheds, based on susceptibility to erosion, has been carried out in different regions of the world [72,73], using different prioritization methods [74]. This study contributes to the lack of knowledge regarding the susceptibility to erosion in northern Mexico. This was assessed by implementing two methods for prioritization based on the analysis of a set of 33 parameters, which differ from other studies [5,7,75]. The inclusion of several geomorphometric parameters and their relationships within several connected watersheds enriched the study of their erosion susceptibility. In this sense, multivariate techniques have proved to be appropriate methods for establishing priorities, reducing the dimensionality of the dataset by losing the least amount of information [76].
This study integrated a multivariate analysis of several geomorphometric parameters that served to identify those watersheds, which may be prone to erosion. That was possible by evaluating the behavior of such geomorphometric parameters and representing them in a geospatial basis [77,78]. Their relationships provided significant information about the main sources of variability among the studied watersheds [74]. Previous research studies have reported that topography, geomorphology, and land use/land cover are the most important factors in the watershed susceptibility to erosion [79,80,81,82]. In this study, the factors with the greatest influence on the hydrological behavior of watersheds and their erosion susceptibility were the average peak flow and the all channel lengths, as it has also been found in previous studies [6,83].
The PCA is considered a statistical exploratory technique, whose results have helped explain the distribution of environmental attributes [69]. Results from the PCA were useful to identify the sources of variance, which were mainly represented by the dominant parameters influencing the data structure. Then, the basin’s hydrological configuration was explained by those geomorphometric parameters explaining the greatest portion of the variance among the watersheds. The PCA results from this study are consistent with the observations made by Meshram and Sharma [32] and Farhan et al. [33].
From the PCA analysis, PC1 and PC2 are mainly influenced by linear geomorphometric parameters. Some of the linear parameters with an influence on PC1 are the average peak flow. This is shown in Figure 3, where the lowest PC1 coefficients correspond to the watersheds with the lowest mean slope values of the small channels. Regarding PC2, drainage density is one of the linear parameters with an influence. Watersheds with low drainage density indicate the presence of permeable surface material, good vegetation cover, and low relief [84,85]. The map of PC2 (Figure 3) was highly influenced by drainage density, since the watersheds with low values of this parameter are located in the south–central part of the study area and grouped in such a map.
PC3 and PC4 are influenced by linear parameters such as mean distance and shape parameters such as elongation ratio, as well as topographic parameters such as minimum height and mean slope. These factors are associated with the main channel, relief, and slopes, among others. In Figure 3, the watersheds with the greatest heights and slopes correspond to the watersheds located in a mountainous zone, while the watersheds with the lowest elevations and slopes correspond to the arid and semi-arid zones of the state of Chihuahua [86]. Regarding PC5, it is mainly influenced by the elongation index (shape parameter) and the general flow length (linear parameter). The high values of the elongation index correspond to enlarged watersheds, which are related to high drainage densities. A watershed with a high drainage density implies a quick hydrological response to rainfall events, while non-enlarged watersheds correspond to fan-shaped watersheds, which are characterized by short channels [87].
Group analysis (GA) was one of the methodologies used in this study to group and then prioritize the watersheds. It was useful to relate watersheds that share the same characteristics based on their geomorphometric parameters. The groups delineated by the analysis have a unique combination in terms of their geomorphometric attributes [40]. The groups of watersheds follow a territorial pattern. Group 1 includes watersheds located in the zones with the least slopes, where the predominant economic activity is agriculture. Group 2 and 3 belong to watersheds located in transition zones because of their moderate slopes. Groups 4 and 5 are watersheds with rugged topography, with vegetation of shrublands and oak forests.
The compound parameter (Cp) was the second methodology employed to prioritize the watersheds. The value of Cp was calculated for each of the 31 watersheds composing the CRB (Figure 5a). Based on the value of Cp, watershed RH24Hf received the highest priority (1), followed by the watershed RH24Le (2). By comparing the results from GA and Cp, Group 4 was identically integrated by the two methodologies. This group is characterized by watersheds having the highest average values of maximum and minimum height, elongation ratio, elongation index, mean watershed slope, slope of the main channel and unit shape factor. The high values of these parameters correspond to watersheds with a high erosion susceptibility [88]. Conversely, Group 5 was formed by watersheds having the highest values of main channel length, watershed length, contour length, all channel lengths, Gravelius compactness coefficient, shape factor, massivity coefficient, mean distance, Kirpich concentration time, sinuosity of the main channel, average peak flow, texture ratio, river frequency and resistance number. This coincide with high values of Cp, which correspond to watersheds distributed in the southwestern zone of the study area and may have a low erosion susceptibility [61].
The two prioritization schemes used in this study gave similar results according to the spatial distribution of watershed groups. The prioritization of watersheds, obtained through GA and Cp, highlighted those watersheds with potential for the implementation of soil and water conservation practices. Based on the ANOVA analyses performed to statistically compare the GA and Cp methodologies, the former resulted in more effectively classifying the watersheds, since it permitted better differentiating the watershed groups.
Results from the GA show that erosion susceptibility is strongly related to linear parameters (surface hydrology) for southwestern watersheds, where steep slopes of both the watershed and the main channel influence soil erosion [89]. Watersheds RH24Lg, RH24Le, RH24Lf, RH24Mc, RH24Lb, RH24Nc, and RH24Ne have the steepest slopes, making them more prone to erosion [29].
One of the advantages of using the watershed as a territorial unit is the analysis of multiple geomorphometric parameters, which are related to the watershed’s hydrological configuration, topography, and shape. Most of the watershed surface attributes depend on local topographic conditions [5]. In this study, the Basin’s altitudinal gradient, a surface attribute, assists in exhibiting the contrasts among watersheds groups, while showing a homogeneous geographic distribution within them. The linear and shape-type parameters are important because of their influence on soil erosion.
The description and spatial grouping of the 31 watersheds through their 26 parameters using multivariate techniques proved to be useful to understand the main factors that control the variance in the CRB. Prioritization through the two types of grouping was also effective in detecting those watersheds susceptible to erosion. The proposed methodology for prioritizing watersheds on a geospatial basis is a feasible approach for identifying watersheds that are susceptible to erosion. However, prioritization with parameters that are based on shape, linear, and relief of the watersheds may not be sufficient. Thus, the incorporation of information regarding the activities on the territory of the CRB would help to improve the efficacy of the classification of watersheds based on their erosion susceptibility. Therefore, future research could include socioeconomic attributes that contribute to soil loss, such as agriculture [39]. Despite the limitations of this study, the contribution of this work represents an advance in the identification of the watersheds that are most susceptible to erosion in the CRB. This in turn contributes to land-use planning, which may help mitigate soil degradation processes.

5. Conclusions

The application of GA and Cp methodologies allowed integrating a large set of geomorphometric parameters, which served to classify watersheds according to their characteristics.
GA more effectively clustered the watersheds of the Conchos River Basin than Cp, since the groups formed by GA were more differentiated based on the analysis of the watersheds’ geomorphometric parameters. The results of GA show that watersheds RH24Lf, RH24Lb, RH24Nc, and RH24Jb might be subjected to strong erosion processes, and are potential candidates to be subjected to soil conservation practices.
The present study demonstrates the usefulness of integrating GIS and multivariate techniques to prioritize watersheds based on their erosion susceptibility. Such an integration approach showed the spatial relationships of the different geomorphometric parameters analyzed. Although the present study permitted a definition of watershed groups according to the values of their geomorphometric parameters and their relation with erosion susceptibility, the integration of additional variables in the analysis may provide a more insightful classification and thus a more reliable watershed prioritization. Such variables could include land use/land cover, soil type, lithology, geomorphology, and socioeconomic activities, among others.

Author Contributions

Conceptualization, A.P.-A. and F.V.-G.; Data curation, G.V.-Q.; Formal analysis, J.A.P.-A., A.P.-A., M.C.V.-A., A.E.R.-R., and M.M.-S.; Investigation, M.C.V.-A.; Methodology, M.M.-S.; Project administration, A.P.-A.; Software, J.A.P.-A. and G.V.-Q.; Supervision, F.V.-G.; Writing—original draft, J.A.P.-A.; Writing—review and editing, M.C.V.-A., A.E.R.-R., and F.V.-G.

Funding

This research received no external funding.

Acknowledgments

We deeply thank Consejo Nacional de Ciencia y Tecnología (CONACYT-Mexico) for supporting the doctoral graduate studies of the first author of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Values of 33 geomorphometric parameters used for the analysis of the 31 watersheds composing the Conchos River Basin.
Table A1. Values of 33 geomorphometric parameters used for the analysis of the 31 watersheds composing the Conchos River Basin.
WidAPLb2LcLiLuNuNo1HminHmaxHmedCcReRfIaJtgα
RH24Hf183.1971.9912.1812.66123.87505.4364933694011007803.001.2514.470.876.760.19
RH24Ia78.6250.6024.1911.1769.21189.8825211192010408003.220.417.041.598.800.09
RH24Ib1209.37259.75140.2034.152261.863915.1950392367147021408004.210.2835.410.9618.700.82
RH24Ja2344.14297.44132.1461.843197.226430.3777163854147021408003.470.4137.911.6313.641.59
RH24Jb4527.91373.79121.3056.885415.4214,899.5128,07413,371158022609003.130.6379.600.7111.962.87
RH24Jc1880.25290.9940.6082.931877.945304.8187844197161022409803.791.2122.673.669.991.17
RH24Ka2345.09283.55104.9945.813530.786993.14877144311620224010003.300.5251.190.8915.061.45
RH24Kb3447.12546.625.47106.052200.407188.1010,60151231645221010805.2512.1232.513.266.382.10
RH24Kc2928.21433.41141.0983.923015.006069.72915343571700232010804.520.4334.892.4110.301.72
RH24Kd1076.82190.5064.8837.642111.472426.44443121002020264014003.280.5728.611.3219.610.53
RH24Ke380.77120.2733.6022.50872.801023.8218128932000246015403.480.6616.921.3322.920.19
RH24Kf920.06232.6093.6564.151030.601673.31231111231690226011204.330.3714.344.4711.200.54
RH24Kg1167.54193.8843.2930.35717.781864.60233111201540180012803.200.8938.470.796.150.76
RH24La953.97169.4965.0731.48781.192467.01374718431530176013003.100.5430.301.048.190.62
RH24Lb5428.30640.78283.10133.1710,026.3213,139.5319,45597862070282013204.910.2940.763.2718.472.62
RH24Lc1760.18228.0293.5551.283884.502710.82287214282320282018203.070.5134.321.4922.070.76
RH24Ld2017.67324.6316.1075.214868.583104.04395219122400298018204.083.1526.832.8024.130.84
RH24Le2267.52286.04145.5173.495722.485511.9612,33861222200278016203.390.3730.852.3825.241.03
RH24Lf4296.09482.75158.2897.019577.979829.2222,33411,0212150286014404.160.4744.282.1922.292.00
RH24Lg2618.69322.00156.7185.607070.696394.0517,22184282420328015603.550.3730.592.8027.001.08
RH24Lh1457.89220.72111.3965.27946.942997.72551026921840236013203.260.3922.342.926.500.79
RH24Ma3778.37465.35215.32130.551552.455018.08735133741680214012204.270.3228.944.514.112.25
RH24Mb3808.53347.58126.6792.623022.418366.9115,13774431970256013803.180.5541.122.257.941.93
RH24Mc2854.62435.90117.9178.485141.297178.6812,80462642330302016404.600.5136.372.1618.011.23
RH24Md2485.07400.5920.76122.531646.853661.01573828222070288012604.532.7120.286.046.631.20
RH24Me804.95211.9295.4068.671032.981770.22275613552170282015204.210.3411.725.8612.830.37
RH24Na1430.41232.19104.8271.05574.841594.10216410231580204011203.460.4120.133.534.020.91
RH24Nb1115.90186.3772.0251.871298.612417.57336416661750226012403.150.5221.512.4111.640.64
RH24Nc4655.43494.37207.40110.907684.7910,752.0423,25311,6002020278012604.090.3741.982.6416.512.30
RH24Nd2761.21363.46166.11118.174289.875719.63825041462080280013603.900.3623.375.0615.541.33
RH24Ne2077.04315.23128.6093.973020.715385.38999649641930250013603.900.4022.104.2514.541.08
WidCoDdjaTcKScpQpTRuFuRnRhRrLoRcDi
RH24Hf183.192.762.630.903.120.96108.129.010.943.543034.9686.9015.281.380.441.28
RH24Ia78.622.420.992.737.692.1769.534.981.263.212511.6393.1020.551.210.391.33
RH24Ib1209.373.240.964.0330.214.11289.5819.400.984.176927.9762.678.241.620.231.29
RH24Ja2344.142.741.012.7328.212.14409.0725.941.283.295870.3734.617.191.370.331.20
RH24Jb4527.913.291.121.8025.412.13576.8375.110.856.207436.7439.736.051.650.411.88
RH24Jc1880.252.823.100.947.390.49364.6030.191.914.676319.7827.017.701.410.281.66
RH24Ka2345.092.981.182.1722.282.29409.1630.930.953.746679.7548.897.901.490.371.25
RH24Kb3447.122.0920.670.090.760.05500.2919.391.813.084608.4020.844.041.040.141.47
RH24Kc2928.212.070.882.6131.351.68459.4521.121.553.134808.9927.645.351.040.201.51
RH24Kd1076.822.251.911.9812.781.72272.5523.261.154.115948.7970.1513.861.130.371.83
RH24Ke380.772.692.741.726.701.49158.4115.071.154.766614.51109.3220.451.340.331.77
RH24Kf920.061.821.223.0920.171.46251.069.942.112.514110.2535.239.720.910.211.38
RH24Kg1167.541.601.201.2711.191.43284.3112.020.892.002874.6559.319.280.800.391.25
RH24La953.972.590.712.1118.792.07255.8522.111.023.934551.4755.9110.381.290.421.52
RH24Lb5428.302.420.533.8465.122.13634.1130.361.813.586825.9821.184.401.210.171.48
RH24Lc1760.181.541.072.2321.191.82352.2512.601.221.634343.0254.9912.370.770.431.06
RH24Ld2017.671.547.200.362.620.21378.2712.171.671.964584.5139.629.180.770.241.27
RH24Le2267.522.430.803.0633.331.98402.0443.131.545.446757.7137.839.721.220.352.24
RH24Lf4296.092.290.902.4133.981.63561.2246.261.485.206543.5229.485.921.140.232.27
RH24Lg2618.692.441.103.0631.201.83433.4253.481.676.588008.7738.3210.191.220.322.69
RH24Lh1457.892.060.932.9225.531.71319.2524.961.713.784852.6636.1610.691.030.381.84
RH24Ma3778.371.330.433.5057.311.65524.8415.802.121.952842.1516.394.600.660.221.46
RH24Mb3808.532.200.932.0528.211.37527.0243.551.503.975624.0327.647.371.100.401.81
RH24Mc2854.622.511.172.2124.451.50453.3829.371.474.497594.5738.486.931.260.191.78
RH24Md2485.071.477.800.423.090.17421.7314.322.462.314242.8223.507.190.740.191.57
RH24Me804.952.201.363.3619.591.39234.1413.002.423.426201.6441.0613.311.100.231.56
RH24Na1430.411.110.882.7724.951.48316.109.321.881.512273.4528.718.790.560.331.36
RH24Nb1115.902.171.422.1615.541.39277.6718.051.553.014896.2343.5712.131.080.401.39
RH24Nc4655.432.310.733.0445.231.87585.2647.041.634.996420.6125.075.621.150.242.16
RH24Nd2761.212.070.873.1635.741.41445.5822.702.252.995799.9823.707.701.040.261.44
RH24Ne2077.042.590.892.8229.091.37384.0431.712.064.816482.0526.607.931.300.261.86
Wid = Watershed identification, A = Area, P = Perimeter, Lb2 = Length of watershed, Lc = Main channel length, Li = Length of contour lines, Lu = Length of channels, Nu = Number of channels, No1 = Total number of channels of order 1 of the watershed, Hmin = Minimum height, Hmax = Maximum height, Hmed = Medium height, Cc = Gravelius compactness coefficient, Re = Elongation elation, Rf = Form factor, Ia = Elongation index, J = Medium slope of the watershed, tgα = Massive coefficient, Co = Orographic coefficient, Dd = Drainage density, j = Mean slope of main channel, a = Medium distance, TcK = Kirpich concentration time, Scp = Sinuosity of main channel, Qp = Average peak flow, T = Texture ratio, Ru = Unitary shape factor, Fu = River frequency, Rn = Resistance number, Rh = Relief ratio, Rr = Relative relief, Lo = Length of general flow, Rc = Circularity ratio, Di = Drainage intensity.
Table A2. Eigenvectors of the correlation matrix.
Table A2. Eigenvectors of the correlation matrix.
GPPC1PC2PC3PC4PC5
Lc0.26780.0223−0.2753−0.13300.0095
Lb20.2345−0.2721−0.0327−0.00590.0880
Li0.28720.03430.01490.1182−0.1631
Lu0.30140.03440.1560−0.1236−0.0398
Hmax0.2170−0.1034−0.05430.3656−0.1169
Hmin0.0664−0.0941−0.11830.4424−0.3506
Cg0.1308−0.24250.06770.00280.2681
Re−0.0340−0.19920.40390.00720.0603
Rf0.21800.13090.2002−0.1898−0.3142
Ia0.0510−0.3141−0.15330.10960.3442
J0.12690.1467−0.04570.4034−0.1466
tgα0.2759−0.07210.1407−0.2389−0.0850
Dd0.07200.31740.1766−0.02480.3345
j−0.0484−0.21180.40350.05930.0789
a0.12710.1028−0.4225−0.09730.1850
TcK0.24830.0091−0.2925−0.1728−0.0094
Scp0.06940.2821−0.2272−0.18220.0258
Qp0.3031−0.11040.0908−0.0989−0.1442
T0.26420.15390.14300.0413−0.0230
Ru0.0641−0.3193−0.15820.12610.3286
Fu0.17370.24780.12910.19610.2488
Rn0.23050.17240.08190.24240.1779
Rh−0.22070.25120.01330.1126−0.0453
Rr−0.24150.1661−0.09860.22730.0458
Lo0.07200.31740.1766−0.02480.3345
Di0.19110.07950.02410.29920.1179
GP = Geomorphometric parameter, PC1 = Principal component 1, PC2 = Principal component 2, PC3 = Principal component 3, PC4 = Principal component 4, PC5 = Principal component 5, Lc = Length of main channel, Lb2 = Length of watershed, Li = Length of contour lines, Lu = Length of channels, Hmin = Minimum height, Hmax = Maximum height, Cc = Gravelius compactness coefficient, Re = Elongation ratio, Rf = Form factor, Ia = Elongation index, J = Mean slope of watershed, tgα = Mass coefficient, Dd = Drainage density, j = mean slope of the main channel, a = Medium distance, TcK = Kirpich concentration time, Scp = Sinuosity of the main channel, Qp = Average peak flow, T = Texture ratio, Ru = Unit shape factor, Fu = River frequency, Rn = Resistance number, Rh = Relief ratio, Rr = Relative relief, Lo = General flow length, Di = Drainage intensity. Bold letters indicate the dominant coefficient.

Appendix B

Table A3. Abbreviations.
Table A3. Abbreviations.
AbbreviationMeaning
ANOVAAnalysis of variance
CpCompound parameter
CRBConchos River Basin
DTMDigital terrain model
GISGeographic information systems
GAGroup analysis
INEGI *National Institute of Statistics, Geography and Informatics
MANOVAMultivariate analysis of variance
CONABIO *National Commission for the Knowledge and Use of Biodiversity
PCAPrincipal component analysis
PCsPrincipal components
* Acronyms in Spanish.

References

  1. Alexakis, D.D.; Hadjimitsis, D.G.; Agapiou, A. Integrated use of remote sensing, GIS and precipitation data for the assessment of soil erosion rate in the catchment area of “Yialias” in Cyprus. Atmos. Res. 2013, 131, 108–124. [Google Scholar] [CrossRef]
  2. Ganasri, B.; Ramesh, H. Assessment of soil erosion by RUSLE model using remote sensing and GIS—A case study of Nethravathi Basin. Geosci. Front. 2016, 7, 953–961. [Google Scholar] [CrossRef]
  3. Gajbhiye, S.; Sharma, S.K.; Awasthi, M.K. Application of Principal Components Analysis for Interpretation and Grouping of Water Quality Parameters. Int. J. Hybrid Inf. Technol. 2015, 8, 89–96. [Google Scholar] [CrossRef]
  4. Gutiérrez, M.; Carreón-Hernández, E. Salinidad en el bajo Rio Conchos: Aportes y tendencias. Terra Latinoam. 2004, 22, 499–506. [Google Scholar]
  5. Biswas, S.; Sudhakar, S.; Desai, V.R. Prioritisation of subwatersheds based on morphometric analysis of drainage basin: A remote sensing and gis approach. J. Indian Soc. Remote Sens. 1999, 27, 155–166. [Google Scholar] [CrossRef]
  6. Adhami, M.; Sadeghi, S.H. Sub-watershed prioritization based on sediment yield using game theory. J. Hydrol. 2016, 541, 977–987. [Google Scholar] [CrossRef]
  7. Malik, A.; Kumar, A.; Kushwaha, D.P.; Kisi, O.; Salih, S.Q.; Al-Ansari, N.; Yaseen, Z.M. The Implementation of a Hybrid Model for Hilly Sub-Watershed Prioritization Using Morphometric Variables: Case Study in India. Water 2019, 11, 1138. [Google Scholar] [CrossRef]
  8. Robinne, F.-N.; Bladon, K.D.; Silins, U.; Emelko, M.B.; Flannigan, M.D.; Parisien, M.-A.; Wang, X.; Kienzle, S.W.; Dupont, D.P. A Regional-Scale Index for Assessing the Exposure of Drinking-Water Sources to Wildfires. Forests 2019, 10, 384. [Google Scholar] [CrossRef]
  9. Khosravi, K.; Pham, B.T.; Chapi, K.; Shirzadi, A.; Shahabi, H.; Revhaug, I.; Prakash, I.; Bui, D.T. A comparative assessment of decision trees algorithms for flash flood susceptibility modeling at Haraz watershed, northern Iran. Sci. Total Environ. 2018, 627, 744–755. [Google Scholar] [CrossRef]
  10. Chen, T.; Niu, R.Q.; Li, P.X.; Zhang, L.P.; Du, B. Regional soil erosion risk mapping using RUSLE, GIS, and remote sensing: A case study in Miyun Watershed, North China. Environ. Earth Sci. 2011, 63, 533–541. [Google Scholar] [CrossRef]
  11. Rahaman, S.A.; Ajeez, S.A.; Aruchamy, S.; Jegankumar, R. Prioritization of Sub Watershed Based on Morphometric Characteristics Using Fuzzy Analytical Hierarchy Process and Geographical Information System—A Study of Kallar Watershed, Tamil Nadu. Aquat. Procedia 2015, 4, 1322–1330. [Google Scholar] [CrossRef]
  12. Sharma, S.K.; Gajbhiye, S.; Nema, R.K.; Tignath, S. Assessing vulnerability to soil erosion of a watershed of tons River basin In Madhya Pradesh using Remote sensing and GIS. Int. J. Environ. Res. Dev. 2014, 4, 153–164. [Google Scholar]
  13. Keesstra, S.; Pereira, P.; Novara, A.; Brevik, E.C.; Azorin-Molina, C.; Parras-Alcántara, L.; Jordán, A.; Cerdà, A. Effects of soil management techniques on soil water erosion in apricot orchards. Sci. Total Environ. 2016, 551, 357–366. [Google Scholar] [CrossRef] [PubMed]
  14. Clarke, J.I. Morphometry from Maps. Essays in Geomorphology; Elsevier Publishing Co.: New York, NY, USA, 1966; pp. 235–274. [Google Scholar]
  15. Horton, R.E. Erosional development of streams and their drainage basins: A hydrophysical approach to quantitative morphology. Geol. Soc. Am. Bull. 1945, 56, 275–370. [Google Scholar] [CrossRef]
  16. Strahler, A.N. Quantitative geomorphology of drainage basins and channel networks. In Handbook of Applied Hydrology, 1st ed.; Chow, V.T., Ed.; McGraw-Hill: New York, NY, USA, 1964; Volume 1, pp. 4–39. [Google Scholar]
  17. Singh, P.; Gupta, A.; Singh, M. Hydrological inferences from watershed analysis for water resource management using remote sensing and GIS techniques. Egypt. J. Remote Sens. Space Sci. 2014, 17, 111–121. [Google Scholar] [CrossRef] [Green Version]
  18. Nautiyal, M.D. Morphometric analysis of a drainage basin using aerial photographs: A case study of Khairkuli basin District Deharadun. J. Indian Soc. Remote Sens. 1994, 22, 251–262. [Google Scholar] [CrossRef]
  19. Shrimali, S.; Aggarwal, S.; Samra, J. Prioritizing erosion-prone areas in hills using remote sensing and GIS—A case study of the Sukhna Lake catchment, Northern India. Int. J. Appl. Earth Obs. Geoinf. 2001, 3, 54–60. [Google Scholar] [CrossRef]
  20. Thakker, A.K.; Dhiman, S.P. Morphometric analysis and prioritization of miniwatersheds in Mohr watershed, Gujarat using Remote sensing and GIS techniques. J. Indian Soc. Remote Sens. 2007, 35, 313–321. [Google Scholar] [CrossRef]
  21. Sharma, S.K.; Rajput, G.S.; Tignath, S.; Pandey, R.P. Morphometric analysis of a watershed using GIS. J. Indian Water Res. Soc. 2010, 30, 33–39. [Google Scholar]
  22. Viramontes-Olivas, O.A.; Reyes-Gómez, V.M.; Escoboza-García, L.F.; Román-Calleros, J.A.; Pérez-Márquez, A.; Pinedo-Álvarez, C.; Sánchez-Fernández, P.B.; Miranda Ojeda, N. Hidrología de los suelos de las cuencas media y baja del Río Conchos, Chihuahua, México. Rev. Latinoam. Recur. Nat. 2008, 4, 31–45. [Google Scholar]
  23. Tilahun, S.A.; Guzman, C.D.; Zegeye, A.D.; Ayana, E.K.; Collick, A.S.; Yitaferu, B.; Steenhuis, T.S. Spatial and temporal patterns of soil erosion in the semi-humid Ethiopian highlands: A case study of Debre Mawi watershed. In Nile River Basin, 1st ed.; Assefa, M., Melesse, W., Abtew Shimelis, G., Eds.; Springer: Cham, Switzerland, 2014; pp. 149–163. [Google Scholar]
  24. Zavoianu, I. Morphometry of Drainage Bassins, 2nd ed.; Editura Academiei: Bucharest, Romania; Elsevier Science Publishers: Amsterdam, The Netherlands, 1985; pp. 9–235. [Google Scholar]
  25. Schumm, S.A. Evolution of drainage systems and slopes in badlands at perth amboy, new jersey. GSA Bull. 1956, 67, 597. [Google Scholar] [CrossRef]
  26. Saha, S.; Mukhopadhyay, S. Determination of Significant Geomorphic Parameters of Kunur River Basin through Principal Component Analysis. J. Eng. Comput. Appl. Sci. 2012, 4, 232–243. [Google Scholar]
  27. Sharma, S.K.; Tignath, S.; Gajbhiye, S.; Patil, R. Application of principal component analysis in grouping geomorphic parameters of Uttela watershed for hydrological modeling. Int. J. Remote Sens. Geosci. 2013, 2, 63–70. [Google Scholar]
  28. Gavit, B.K.; Purohit, R.C.; Singh, P.K.; Kothari, M.K.; Jain, H.K. Grouping of Significant Geomorphic Parameters using Multivariate Parameters using Multivariate. Res. J. Recent Sci. 2016, 5, 32–38. [Google Scholar]
  29. Youssef, A.M.; Pradhan, B.; Hassan, A.M. Flash flood risk estimation along the St. Katherine road, southern Sinai, Egypt using GIS based morphometry and satellite imagery. Environ. Earth Sci. 2011, 62, 611–623. [Google Scholar] [CrossRef]
  30. Makwana, J.; Tiwari, M.K. Prioritization of agricultural sub-watersheds in semiarid middle region of Gujarat using remote sensing and GIS. Environ. Earth Sci. 2016, 75, 137. [Google Scholar] [CrossRef]
  31. Sharma, S.K.; Gajbhiye, S.; Tignath, S. Application of principal component analysis in grouping geomorphic parameters of a watershed for hydrological modeling. Appl. Water Sci. 2014, 5, 89–96. [Google Scholar] [CrossRef] [Green Version]
  32. Meshram, S.G.; Sharma, S.K. Prioritization of watershed through morphometric parameters: A PCA-based approach. Appl. Water. Sci. 2017, 7, 1505–1519. [Google Scholar] [CrossRef]
  33. Farhan, Y.; Anbar, A.; Al-Shaikh, N.; Mousa, R. Prioritization of semi-arid agricultural watershed using morphometric and principal component analysis, remote sensing, and GIS techniques, the Zerqa River Watershed, Northern Jordan. Agric. Sci. 2017, 8, 113–148. [Google Scholar] [CrossRef]
  34. Bach, K.; Schawe, M.; Beck, S.; Gerold, G.; Gradstein, S.R.; Moraes, R. Vegetación, suelos y clima en los diferentes pisos altitudinales de un bosque montano de Yungas, Bolivia: Primeros resultados. Ecol. Boliv. 2003, 38, 3–14. [Google Scholar]
  35. Sanderson, M.; Santini, M.; Valentini, R.; Pope, E. Relationships between Forests and Weather. EC Directorate General of the Environment; Met Office: Devon, UK, 2012; Available online: http://ec.europa.eu/environment/forests/pdf/EU_Forests_annex1.pdf (accessed on 29 July 2019).
  36. Batayneh, A.; Zumlot, T. Multivariate statistical approach to geochemical methods in water quality factor identification; application to the shallow aquifer system of the Yarmouk basin of north Jordan. Res. J. Environ. Earth Sci. 2012, 4, 756–768. [Google Scholar]
  37. Oketola, A.A.; Adekolurejo, S.M.; Osibanjo, O. Water Quality Assessment of River Ogun Using Multivariate Statistical Techniques. J. Environ. Prot. 2013, 4, 466–479. [Google Scholar] [CrossRef] [Green Version]
  38. Castillo-Rodríguez, M.; López-Blanco, J.; Muñoz-Salinas, E. A geomorphologic GIS-multivariate analysis approach to delineate environmental units, a case study of La Malinche volcano (central México). Appl. Geogr. 2010, 30, 629–638. [Google Scholar] [CrossRef]
  39. Miranda, J.; Andrade, E.; López-Suárez, A.; Ledesma, R.; Cahill, T.; Wakabayashi, P. A receptor model for atmospheric aerosols from a southwestern site in Mexico City. Atmos. Environ. 1996, 30, 3471–3479. [Google Scholar] [CrossRef]
  40. Tritsch, I.; Sist, P.; Narvaes, I.D.S.; Mazzei, L.; Blanc, L.; Bourgoin, C.; Cornu, G.; Gond, V. Multiple Patterns of Forest Disturbance and Logging Shape Forest Landscapes in Paragominas, Brazil. Forests 2016, 7, 315. [Google Scholar] [CrossRef]
  41. Prieto-Amparán, J.A.; Santellano-Estrada, E.; Villarreal-Guerrero, F.; Martinez-Salvador, M.; Pinedo-Alvarez, A.; Vázquez-Quintero, G.; Valles-Aragón, M.C.; Manjarrez-Domínguez, C. Spatial Analysis of Temperate Forest Structure: A Geostatistical Approach to Natural Forest Potential. Forests 2019, 10, 168. [Google Scholar] [CrossRef]
  42. Altaf, S.; Meraj, G.; Ahmad Romshoo, S. Morphometry and land cover based multicriteria analysis for assessing the soil erosion susceptibility of the western Himalayan watershed. Environ. Monit. Assess. 2014, 186, 8391–8412. [Google Scholar] [CrossRef]
  43. Montero-Martínez, M.J.; Ibáñez-Hernández, Ó.F. La cuenca del río Conchos: Una mirada desde las ciencias ante el cambio climático, 1st ed.; Instituto Mexicano de la Tecnología del Agua: Morelos, Mexico, 2017; p. 10. [Google Scholar]
  44. Ordoñez, R. Modelación de la erosión hídrica en la Cuenca Alta del Río Conchos, Chihuahua. Master’s Thesis, Universidad Autónoma de Chihuahua, Chih, Mexico, June 2017. [Google Scholar]
  45. Mahlknecht, J.; Horst, A.; Aravena, R.; Hernández-Limón, G.; Hernández-Limón, G. Groundwater geochemistry of the Chihuahua City region in the Rio Conchos Basin (northern Mexico) and implications for water resources management. Hydrol. Process. 2008, 22, 4736–4751. [Google Scholar] [CrossRef]
  46. Sánchez, A. 1944 Water Treaty between Mexico and the United States: Present situation and future potential. Front. Norte 2006, 18, 125–144. [Google Scholar]
  47. CONAGUA (Comisión Nacional del Agua). Comisión de cuenca del río Conchos. In Diagnóstico Versión 1.0; Gerencia Estatal Chihuahua: Chihuahua, México, 2001. [Google Scholar]
  48. Núñez López, D.; Treviño Garza, E.J.; Reyes Gómez, V.M.; Muñoz Robles, C.A.; Aguirre Calderón, O.A.; Jiménez Pérez, J. Uso de modelos de regresión para interpolar espacialmente la precipitación media mensual en la cuenca del río Conchos. Rev. Mex. Cienc. Agrícolas 2014, 5, 201–213. [Google Scholar] [CrossRef]
  49. Viramontes-Olivas, O.A.; Escoboza-García, L.F.; Pinedo-Alvarez, C.; Pinedo-Alvarez, A.; Reyes Gómez, V.; Román Calleros, J.; Pérez-Márquez, A. Morfometría de la cuenca del río San Pedro, Conchos, Chihuahua. Tecnociencia Chihuah. 2008, 1, 21–31. [Google Scholar]
  50. Aboites-Aguilar, L. Hacia una historia del río Conchos; CIESAS, Juárez: Chihuahua, México, 2002. [Google Scholar]
  51. INEGI (Instituto Nacional de Estadística, Geografía e Informática). Guía para la Interpretación de la Cartografía uso del Suelo y Vegetación, Escala 1:250,000 Serie V; INEGI: Aguascalientes, México, 2014; p. 30. Available online: https://www.inegi.org.mx/contenidos/temas/mapas/usosuelo/metadatos/guia_interusosuelov.pdf (accessed on 31 July 2019).
  52. CONABIO (Comisión Nacional para el Conocimiento y Uso de la Biodiversidad). Portal de Geoinformación. Available online: http://www.conabio.gob.mx/informacion/gis/ (accessed on 29 July 2019).
  53. SIATL. Simulador de Flujos de Agua de Cuencas Hidrográficas. Available online: http://antares.inegi.org.mx/analisis/red_hidro/siatl/ (accessed on 29 July 2019).
  54. INEGI (Instituto Nacional de Estadística, Geografía e Informática). Sistema de descarga del Continuo de Elevaciones Mexicano. Available online: https://www.inegi.org.mx/app/geo2/elevacionesmex/ (accessed on 29 July 2019).
  55. Hamdy, O.; Zhao, S.; Salheen, M.A.; Eid, Y.Y. Identifying the Risk Areas and Urban Growth by ArcGIS-Tools. Geosciences 2016, 6, 47. [Google Scholar] [CrossRef]
  56. Miller, V.C. A Quantitative Geomorphic Study of Drainage Basin Characteristics in the Clinch Mountain Area, Virginia and Tennessee; Office of Naval Research (USA): Arlington, VA, USA, 1953; pp. 1–51. [Google Scholar]
  57. Mueller, J.E. An introduction to the hydraulic and topographic sinuosity indexes 1. Ann. Assoc. Am. Geogr. 1968, 58, 371–385. [Google Scholar] [CrossRef]
  58. Kirpich, Z.P. Time of concentration of small agricultural watersheds. Civ. Eng. 1940, 10, 362. [Google Scholar]
  59. Şen, Z. Modified hydrograph method for arid regions. Hydrol. Process. 2008, 22, 356–365. [Google Scholar] [CrossRef]
  60. Faniran, A. The Index of Drainage Intensity—A Provisional New Drainage Factor. Aust. J. Sci. 1968, 31, 328–330. [Google Scholar]
  61. Patel, D.P.; Gajjar, C.A.; Srivastava, P.K. Prioritization of Malesari mini-watersheds through morphometric analysis: A remote sensing and GIS perspective. Environ. Earth Sci. 2013, 69, 2643–2656. [Google Scholar] [CrossRef]
  62. Agarwal, K.K.; Bali, R.; Kumar, M.G.; Srivastava, P.; Singh, P.V. Active tectonics in and around Kimin-Ziro area, Lower Subansiri District, Arunachal Pradesh, NE India. Z. Geomorphol. 2009, 53, 109–120. [Google Scholar] [CrossRef]
  63. Dillon, W.R.; Goldstein, M. Multivariate Analysis. Methods and Applications; John Wiley & Sons Inc.: New York, NY, USA, 1984; pp. 1–608. [Google Scholar]
  64. Jolliffe, I.T. Principal Component Analysis, 2nd ed.; John Wiley & Sons Inc.: Charlottesville, VA, USA, 2002; pp. 150–166. [Google Scholar]
  65. Yidana, S.M.; Banoeng-Yakubo, B.; Akabzaa, T.M. Analysis of groundwater quality using multivariate an spatial analyses in the Keta basin, Ghana. J. Afr. Earth Sci. 2010, 58, 220–234. [Google Scholar] [CrossRef]
  66. Bro, R.; Smilde, A.K. Principal component analysis. Anal. Methods 2014, 6, 2812. [Google Scholar] [CrossRef]
  67. Ward, J.H. Hierarchical Grouping to Optimize an Objective Function. J. Am. Stat. Assoc. 1963, 58, 236–244. [Google Scholar] [CrossRef]
  68. Eder, B.K.; Davis, J.M.; Bloomfield, P. An Automated Classification Scheme Designed to Better Elucidate the Dependence of Ozone on Meteorology. J. Appl. Meteorol. 1994, 33, 1182–1199. [Google Scholar] [CrossRef] [Green Version]
  69. Johnson, R.A.; Wichern, D.W. Applied Multivariate Statistical Analysis, 6th ed.; Pearson Prentice-Hall Inc.: Upper Saddle River, NJ, USA, 2007; pp. 296–301. [Google Scholar]
  70. Nooka Ratnam, K. Check Dam positioning by prioritization ofmicro watershed using SYI model and morphometric analysis—Remote sensing and GIS perspective. J. Indian Soc. Remote Sens. 2005, 33, 25–38. [Google Scholar] [CrossRef]
  71. Fajardo, J.; Lessmann, J.; Bonaccorso, E.; Devenish, C.; Munoz, J. Combined Use of Systematic Conservation Planning, Species Distribution Modelling, and Connectivity Analysis Reveals Severe Conservation Gaps in a Megadiverse Country (Peru). PLoS ONE 2014, 9, e114367. [Google Scholar] [CrossRef] [PubMed]
  72. Paul, I.I.; Bayode, E.N. Watershed Characteristics and Their Implication for Hydrologic Response in the Upper Sokoto Basin, Nigeria. J. Geogr. Geol. 2012, 4, 147–155. [Google Scholar] [CrossRef]
  73. Darabi, H.; Shahedi, K.; Solaimani, K.; Miryaghoubzadeh, M. Prioritization of subwatersheds based on flooding conditions using hydrological model, multivariate analysis and remote sensing technique. Water Environ. J. 2014, 28, 382–392. [Google Scholar] [CrossRef]
  74. Rawat, K.S.; Tripathi, V.K.; Mishra, A.K. Sediment yield index mapping and prioritization of Madia subwatershed, Sagar District of Madhya Pradesh (India). Arab. J. Geosci. 2014, 7, 3131–3145. [Google Scholar] [CrossRef]
  75. Suji, V.R.; Sheeja, R.V.; Karuppasamy, S.; Suji, V.R.; Sheeja, R.V.; Karuppasamy, S. Prioritization using Morphometric Analysis and Land use/Land Cover Parameters for Vazhichal Watershed using Remote Sensing and GIS Techniques. Int. J. Innov. Res. Sci. Tech. 2015, 1, 2349–6010. [Google Scholar]
  76. Prieto-Amparán, J.A.; Rocha-Gutiérrez, B.A.; Ballinas-Casarrubias, M.D.L.; Valles-Aragón, M.C.; Peralta-Pérez, M.D.R.; Pinedo-Alvarez, A. Multivariate and Spatial Analysis of Physicochemical Parameters in an Irrigation District, Chihuahua, Mexico. Water 2018, 10, 1037. [Google Scholar] [CrossRef]
  77. Singh, P.K.; Kumar, V.; Purohit, R.C.; Kothari, M.; Dashora, P.K. Application of principal component analysis in grouping geomorphic parameters for hydrologic modeling. Water Resour. Manag. 2009, 23, 325. [Google Scholar] [CrossRef]
  78. Yunus, A.P.; Oguchi, T.; Hayakawa, Y.S. Morphometric Analysis of Drainage Basins in the Western Arabian Peninsula Using Multivariate Statistics. Int. J. Geosci. 2014, 5, 527–539. [Google Scholar] [CrossRef] [Green Version]
  79. Adediran, A.O.; Parcharidis, I.; Poscolieri, M.; Pavlopoulos, K. Computer-assisted discrimination of morphological units on north-central Crete (Greece) by applying multivariate statistics to local relief gradients. Geomorphology 2004, 58, 357–370. [Google Scholar] [CrossRef]
  80. Javed, A.; Khanday, M.Y.; Ahmed, R. Prioritization of sub-watersheds based on morphometric and land use analysis using remote sensing and GIS techniques. J. Indian Soc. Remot 2009, 37, 261–274. [Google Scholar] [CrossRef]
  81. Kompani-Zare, M.; Soufi, M.; Hamzehzarghani, H.; Dehghani, M. The effect of some watershed, soil characteristics and morphometric factors on the relationship between the gully volume and length in Fars Province, Iran. Catena 2011, 86, 150–159. [Google Scholar] [CrossRef]
  82. Welde, K. Identification and prioritization of subwatersheds for land and water management in Tekeze dam watershed, Northern Ethiopia. Int. Soil Water Conserv. Res. 2016, 4, 30–38. [Google Scholar] [CrossRef] [Green Version]
  83. Subyani, A.M.; Qari, M.H.; Matsah, M.I. Digital elevation model and multivariate statistical analysis of morphometric parameters of some wadis, western Saudi Arabia. Arab J. Geosci. 2012, 5, 147–157. [Google Scholar] [CrossRef]
  84. Luo, W. Quantifying groundwater-sapping landforms with a hypsometric technique. J. Geophys. Res. Space Phys. 2000, 105, 1685–1694. [Google Scholar] [CrossRef] [Green Version]
  85. Harlin, J.M.; Wijeyawickrema, C. Irrigation and groundwater depletion in caddo county, oklahoma. JAWRA J. Am. Water Resour. Assoc. 1985, 21, 15–22. [Google Scholar] [CrossRef]
  86. CONABIO (Comisión Nacional para el Conocimiento y Uso de la Biodiversidad). La biodiversidad en Chihuahua: Estudio de Estado, 1st ed.; Comisión Nacional para el Conocimiento y Uso de la Biodiversidad: Méx, México, 2014; pp. 23–25. [Google Scholar]
  87. Rai, P.K.; Mishra, V.N.; Mohan, K. A study of morphometric evaluation of the Son basin, India using geospatial approach. Remote Sens. Appl. Soc. Environ. 2017, 7, 9–20. [Google Scholar]
  88. Rai, P.K.; Chaubey, P.K.; Mohan, K.; Singh, P. Geoinformatics for assessing the inferences of quantitative drainage morphometry of the Narmada Basin in India. Appl. Geomat. 2017, 9, 167–189. [Google Scholar] [CrossRef]
  89. Tucker, G.E.; Bras, R.L. Hillslope processes, drainage density, and landscape morphology. Water Resour. Res. 1998, 34, 2751–2764. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Land-use/land-cover types of the Conchos River Basin (CRB), (b) Delimitation of the 31 watersheds of the CRB, (c) Location of the CRB in Mexico.
Figure 1. (a) Land-use/land-cover types of the Conchos River Basin (CRB), (b) Delimitation of the 31 watersheds of the CRB, (c) Location of the CRB in Mexico.
Sustainability 11 05140 g001
Figure 2. Classification of 31 watersheds based on the values of five principal components resulting from a principal component analysis. Conchos River Basin, Chihuahua, Mexico. (a) PC1, (b) PC2, (c) PC3, (d) PC4, (e) PC5. PC1 = Principal component 1, PC2 = Principal component 2, PC3 = Principal component 3, PC4 = Principal component 4, PC5 = Principal component 5.
Figure 2. Classification of 31 watersheds based on the values of five principal components resulting from a principal component analysis. Conchos River Basin, Chihuahua, Mexico. (a) PC1, (b) PC2, (c) PC3, (d) PC4, (e) PC5. PC1 = Principal component 1, PC2 = Principal component 2, PC3 = Principal component 3, PC4 = Principal component 4, PC5 = Principal component 5.
Sustainability 11 05140 g002
Figure 3. Dendrogram classifying 31 watersheds by group analysis. Conchos River Basin, Chihuahua, Mexico. The red line was drawn to define the groups.
Figure 3. Dendrogram classifying 31 watersheds by group analysis. Conchos River Basin, Chihuahua, Mexico. The red line was drawn to define the groups.
Sustainability 11 05140 g003
Figure 4. Geospatial distribution of watershed groups classified by group analysis. Very High (), High (), Moderate (), Low (), Very Low (). Conchos River Basin, Chihuahua, Mexico.
Figure 4. Geospatial distribution of watershed groups classified by group analysis. Very High (), High (), Moderate (), Low (), Very Low (). Conchos River Basin, Chihuahua, Mexico.
Sustainability 11 05140 g004
Figure 5. (a) Watersheds and their compound parameter (Cp). (b) Geospatial distribution of watershed groups by Cp reclassification. Very High 1 (), High 2 (), Moderate 3 (), Low 4 (), Very Low 5 (). Conchos River Basin, Chihuahua, Mexico.
Figure 5. (a) Watersheds and their compound parameter (Cp). (b) Geospatial distribution of watershed groups by Cp reclassification. Very High 1 (), High 2 (), Moderate 3 (), Low 4 (), Very Low 5 (). Conchos River Basin, Chihuahua, Mexico.
Sustainability 11 05140 g005
Table 1. Geomorphometric parameters.
Table 1. Geomorphometric parameters.
Geomorphometric ParameterEquationReferences
Basic Parameters 1
Area (A)A = Watershed surface area [km2]Horton [15]
Perimeter (P)P = Watershed perimeter [km]Horton [15]
Length (Lb2)Lb2 = Watershed length [km]
Stream order (u)u = Stream order [unitless]Strahler [16]
Main Channel Length (Lc)Lc = Main flow channel length [km]
All Channel Lengths (Lu)Lu = Length of all the flow channels in the watershed [km]Horton [15]
Contour Length (Li)Li = Contour lines’ length [km]
Number of Flow Channels (Nu)Nu = Number of flow channels [unitless]
Number of First-Order Flow Channels (No1)No1 = Number of total first-order flow channels in the watershed [unitless]
Maximum Height (Hmax)Hmax = Watershed maximum height [m]
Minimum Height (Hmin)Hmin = Watershed minimum height [m]
Medium Height (Hmed)Hmed = Watershed medium height [m]
Shape Parameters
Gravelius Compactness Coefficient (Cc)Cc = P/2√πAZavoianu [24]
Elongation Ratio (Re)Re = 1.1284 (√A/Lc)Schumm [25]
Shape Factor (Rf)Rf = A/Lb2Horton [15]
Elongation Index (Ia)Ia = Lb2/W
where: W = watershed width (Km)
Horton [15]
Unit Shape Factor (RU)RU = Lb2/A0.5Horton [15]
Circularity Ratio (Rc)Rc = 4πA/P2Miller [56]
Relief Parameters
Mean Watershed Slope (J)J = (ΣLi E/A) × 100
where: E = equidistance among contour lines (km)
Horton [15]
Massivity Coefficient (tgα)tgα = Hmed/A
Relief Relationship(Rh)Rh = Hmax/LbSchumm [25]
Relative Relief (Rr)Rr = Hmax/PSchumm [25]
Orographic Coefficient (Co)Co = Hmed × tgα
Linear parameter
Drainage Density (Dd)Dd = ΣLu/AHorton [15]
Mean Slope of the Main Channel (j)j = (Hmax − Hmin)/Lc × 100Horton [15]
Mean Distance (Am)Am = Lc/(√A)Horton [15]
Sinuosity of the Main Flow Channel (Scp)Scp = Lc/Lb2Mueller [57]
Kirpich Concentration Time (TcK)TcK = 0.066(Lb2/j)0.77Kirpich [58]
Average Peak Flow (Qp)Qp = 43A0.522Sen [59]
Texture Ratio (T)T = Nu/PHorton [15]
Rivers Frequency (Fu)Fu = Nu/AHorton [15]
Resistance Number (Rn) Rn = Hmax × DdSchumm [25]
General Flow Length (Lo)Lo = 1/2 × DdSchumm [25]
Drainage Intensity (Di)Di = Fu/DdFaniran [60]
1 The basic parameters were calculated with the GIS software.
Table 2. Principal components and geomorphometric parameters explaining the greatest portion of variance in the watersheds of the Conchos River Basin, Chihuahua, Mexico.
Table 2. Principal components and geomorphometric parameters explaining the greatest portion of variance in the watersheds of the Conchos River Basin, Chihuahua, Mexico.
PCEigenvalueVarianceAccumulated VarianceGeomorphometric Parameter (First)Geomorphometric Parameter (Second)
19.00100.34620.3462Average peak flow (Qp)Channel lengths (Lu)
26.08480.23400.5802Unit shape factor (Ru)Drainage density (Dd)
33.42220.13160.7119Mean distance (Am)Elongation ratio (Re)
42.86200.11010.8219Minimum height (Hmin)Mean slope (J)
51.62460.06250.8844Elongation index (Ia)General flow length (Lo)
PC = Principal component, First = Primary importance, Second = Secondary importance.
Table 3. Average values of the geomorphometric parameters by group.
Table 3. Average values of the geomorphometric parameters by group.
GidLb2LcLiLuHminHmaxCcReRfIaJtgαDd
146.1231.31371.431038.5014959953.220.7420.031.696.430.491.97
280.1665.561820.632660.91248914153.700.9123.283.1113.620.791.97
3102.4775.162765.796206.70229410603.951.9934.732.4112.071.482.59
4146.5688.945556.086201.08297015453.860.4030.303.1021.451.172.36
5192.5299.498176.1312,155.08268012304.070.4451.662.2017.312.452.58
GidjaTcKScpQpTRuFuRnRhRrLoDi
11.421.9211.741.51194.518.831.242.562673.6767.0113.480.991.30
22.442.1718.481.37313.2816.481.693.034835.2847.8111.260.981.51
33.702.1822.191.69417.9027.781.503.865915.1734.496.971.301.51
40.982.8731.181.68433.6137.171.734.877040.2634.588.631.182.04
50.822.7742.441.94589.3649.691.444.996806.7128.865.501.291.95
Gid = Group identification, 1 = Group 1 (Very low erosion susceptibility), 2 = Group 2 (Low erosion susceptibility), 3 = Group 3 (Moderate erosion susceptibility), 4 = Group 4 (High erosion susceptibility), 5 = Group 5 (Very high erosion susceptibility), Lc = Length of main channel, Lb2 = Length of watershed, Li = Length of contour lines, Lu = Length of channels, Hmin = Minimum height, Hmax = Maximum height, Cc = Gravelius compactness coefficient, Re = Elongation elation, Rf = Form factor, Ia = Elongation index, J = Mean slope of watershed, tgα = Mass coefficient, Dd = Drainage density, j = Mean slope of main channel, a = Medium distance, TcK = Kirpich concentration time, Scp = Sinuosity of main channel, Qp = Average peak flow, T = Texture ratio, Ru = Unitary shape factor, Fu = River frequency, Rn = Resistance number, Rh = Relief ratio, Rr = Relative relief, Lo = Length of general flow, Di = Drainage intensity.
Table 4. p-values resulting from the ANOVA analyses performed to detect differences among the groups defined by each classification method.
Table 4. p-values resulting from the ANOVA analyses performed to detect differences among the groups defined by each classification method.
Geomorphometric Parameters
LcLb2LiLuHmaxHminCcReRfIaJtgαDd
GA0.00910.0865<0.0001<0.00010.19540.33150.40110.26750.00050.42730.0750<0.00010.0224
Cp0.37110.35760.32070.14620.76020.10290.43680.49830.32140.83360.52360.12250.1563
jaTcKScpQpTRuFuRnRhRrLoDi
GA0.33160.77120.03560.5116<.0001<.00010.53270.00690.00040.12050.00860.02240.0265
Cp0.59020.64540.43770.21580.22870.14420.82770.01740.16260.49130.19790.15630.0302
GA = Group analysis, Cp = Compound parameter, Lc = Length of main channel, Lb2 = Length of watershed, Li = Length of contour lines, Lu = Length of channels, Hmin = Minimum height, Hmax = Maximum height, Cc = Gravelius compactness coefficient, Re = Elongation elation, Rf = Form factor, Ia = Elongation index, J = Mean slope of watershed, tgα = Mass coefficient, Dd = Drainage density, j = mean slope of main channel, a = medium distance, TcK = Kirpich concentration time, Scp = Sinuosity of main channel, Qp = Average peak flow, T = Texture ratio, Ru = Unitary shape factor, Fu = River frequency, Rn = Resistance number, Rh = Relief ratio, Rr = Relative relief, Lo = Length of general flow, Di = Drainage intensity.

Share and Cite

MDPI and ACS Style

Prieto-Amparán, J.A.; Pinedo-Alvarez, A.; Vázquez-Quintero, G.; Valles-Aragón, M.C.; Rascón-Ramos, A.E.; Martinez-Salvador, M.; Villarreal-Guerrero, F. A Multivariate Geomorphometric Approach to Prioritize Erosion-Prone Watersheds. Sustainability 2019, 11, 5140. https://doi.org/10.3390/su11185140

AMA Style

Prieto-Amparán JA, Pinedo-Alvarez A, Vázquez-Quintero G, Valles-Aragón MC, Rascón-Ramos AE, Martinez-Salvador M, Villarreal-Guerrero F. A Multivariate Geomorphometric Approach to Prioritize Erosion-Prone Watersheds. Sustainability. 2019; 11(18):5140. https://doi.org/10.3390/su11185140

Chicago/Turabian Style

Prieto-Amparán, Jesús A., Alfredo Pinedo-Alvarez, Griselda Vázquez-Quintero, María C. Valles-Aragón, Argelia E. Rascón-Ramos, Martin Martinez-Salvador, and Federico Villarreal-Guerrero. 2019. "A Multivariate Geomorphometric Approach to Prioritize Erosion-Prone Watersheds" Sustainability 11, no. 18: 5140. https://doi.org/10.3390/su11185140

APA Style

Prieto-Amparán, J. A., Pinedo-Alvarez, A., Vázquez-Quintero, G., Valles-Aragón, M. C., Rascón-Ramos, A. E., Martinez-Salvador, M., & Villarreal-Guerrero, F. (2019). A Multivariate Geomorphometric Approach to Prioritize Erosion-Prone Watersheds. Sustainability, 11(18), 5140. https://doi.org/10.3390/su11185140

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop