Morphometric Prioritization, Fluvial Classiﬁcation, and Hydrogeomorphological Quality in High Andean Livestock Micro-Watersheds in Northern Peru

: Anthropic activity a ﬀ ects the hydrogeomorphological quality of ﬂuvial systems. River and valley classiﬁcations are fundamental preliminary steps in determining their ecological status, and their prioritization is essential for the proper planning and management of soil and water resources. Given the importance of the High Andean livestock micro-watershed (HAL-MWs) ecosystems in Peru, an integrated methodological framework is presented for morphometric prioritization that uses a Principal Component Analysis (PCA) and Weighted Sum Approach (WSA), geomorphological ﬂuvial classiﬁcations (channel, slope, and valley), and hydrogeomorphological evaluations using the Hydrogeomorphological Index (IHG). Of six HAL-MWs studied in Leimebamba and Molinopampa (Amazonas region), the PCWSA hybrid model identiﬁed the San Antonio HAL-MW as a top priority, needing the rapid adoption of appropriate conservation practices. Thirty-nine types of river course were identiﬁed, by combining 13 types of valley and 11 types of riverbed. The total assessment of the IHG indicated that 7.6% (21.8 km), 14.5% (41.6 km), 27.9% (80.0 km), and 50.0% (143.2 km) of the basin lengths have “Poor”, “Moderate”, “Good”, and “Very good” quality rankings, respectively. The increase in the artiﬁcial use of river channels and ﬂood plains is closely linked to the decrease in hydrogeomorphological quality.


Introduction
In 2011, 11% of Earth's surface and 70% of water extracted from aquifers, rivers, and lakes were purposed for agriculture [1]. By 2050, a world population of 9.1 billion [2] and a 70% increase in food production, compared to 2009, is estimated, which will have a direct impact on land and water resource availability [3]. Thirty-one percent of Earth's 35 million km 3 fresh water resources, on which aquatic and terrestrial ecosystems as well as mankind and diversity depend, are concentrated in Latin America and the Caribbean [4]. In Peru, the average availability of renewable water sources was 65.726 m 3 /inhabitant/year in 2016, ranking 17th out of 180 countries [5]. However, this availability does not correspond to the spatial distribution of the population, is non-uniform in time (mainly because of variable precipitation) and has decreased due to population growth [5,6]. These circumstances have resulted in scarcity and water stress, thus generating social and productive inequalities [5,7].

Study Area
Our study area corresponding to Northern Peru's Amazonas region has an approximate surface area of 39.25 km 2 , much of it covered by unexplored tropical forests. It is geographically located between parallels 3 • 0 15 , 7 • 2 0 south latitude and meridians 77 • 0 15 and 78 • 42 15 west longitude, with an altitudinal gradient between 120 m.a.s.l. in the north, where humid lowland tropical forests predominate, and 4900 m.a.s.l. to the south where humid highland Andean tropical forests, cloud forests, and deforested grasslands predominate. Agriculture is the main economic activity in Amazonas, occupying 24.9% of the territory and being responsible for 51.22% of the regional Gross Domestic Product (GDP) [66]. There are four areas dedicated to livestock in the region [67]: (1) Pomacochas-Jumbilla, (2) Molinopampa-Mendoza, (3) Leimebamba, and (4) Chiriaco. The first three zones are located in areas of cold temperate climate, where dairy cattle predominate; while the last zone, which has a warm and humid climate, is dedicated to raising Zebu cattle. In Leimebamba and Molinopampa, situated in the province of Chachapoyas (Figure 1), open-field cattle raising (extensive cattle farming) is executed alongside the Andean crop farming of potatoes, corn, and beans [68]. These areas include HAL-MW, belonging to the Utcubamba River level 5 Hydrographic Unit (HU N5) [37], located on slopes and mountain tops with altitudes exceeding 2000 m.a.s.l. They cover large areas of Andean grasslands and scrublands, used as natural pastures managed by anthropic burning [67], which in some cases is complemented by the planting of small pasture areas and forages near stables used for young bovine management. In that regard, Chachapoyas is noteworthy for being an exceptionally suitable area for dairy and beef cattle farming, to such an extent that 42.22% of the economically active population is dedicated to the aforementioned activities [68]. However, as a result of poor agricultural practices, unsustainable logging, urban expansion, and the construction of road infrastructure, these high Andean ecosystems are being degraded [23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38].

Base Map and Satellite Framework
To construct the base map and satellite framework, we utilized the HU N5 of the Utcubamba River, contained in the Peruvian hydrographic watershed vector layer, obtained from the National Water Authority's (ANA) Geo-hydro portal [69] . Populated centers and the hydrography from the digitized 13h, 13i, and 14h quadrangles in the National Geographic Institute (IGN) topographical map series (scale of 1:100,000) were downloaded from the Ministry of Education's web portal [70]. Road and bridge infrastructure data were obtained from the Transport and Communication Ministry's website [71]. The Digital Elevation Model (DEM), generated by the Phased Array Type Lband Synthetic Aperture Radar (PALSAR) of Advanced Land Observing Satellite (ALOS) [72] from the Japan Aerospace Exploration Agency (JAXA), was also utilized. The data were downloaded from the National Aeronautics and Space Administration's (NASA) ASD Data Search Vertex web portal [73], with a 12.5 meter spatial resolution. To generate the Coverage and Land Use (LC / LU) maps, we used two images with a spatial resolution of 10 meters acquired on July 23 rd , 2017 from the Sentinel 2A satellite, Path 17, and Row MRP and MRN. These were acquired from the European Space Agency's (ESA) Copernicus Services Data Hub platform, through QGIS's Semi-Automatic Classification Plugin (SCP) [74].

Micro-Watershed Delimitation
Delimitation of the HAL-MWs was done using the DEM and coded from the Utcubamba River HU N5. This process was based on the Pfafstetter method [75,76] while using PgHydro Tools, a QGIS (version 2.18.10) plugin [77], to activate the PgHydro Extension functions for PostgreSQL/PostGIS. The linear water network layer was imported from Google Earth Pro (version 7.3.2.5576) and SAS Planet (version 190707) interfaces, and subsequently updated and complemented with manual mapping [78,79]. This procedure was critical in obtaining the detailed geomorphology of the channel at the micro-watershed level, because the base layer brought the smoothed and generalized rivers (scale 1:100000).

Base Map and Satellite Framework
To construct the base map and satellite framework, we utilized the HU N5 of the Utcubamba River, contained in the Peruvian hydrographic watershed vector layer, obtained from the National Water Authority's (ANA) Geo-hydro portal [69]. Populated centers and the hydrography from the digitized 13h, 13i, and 14h quadrangles in the National Geographic Institute (IGN) topographical map series (scale of 1:100,000) were downloaded from the Ministry of Education's web portal [70]. Road and bridge infrastructure data were obtained from the Transport and Communication Ministry's website [71]. The Digital Elevation Model (DEM), generated by the Phased Array Type L-band Synthetic Aperture Radar (PALSAR) of Advanced Land Observing Satellite (ALOS) [72] from the Japan Aerospace Exploration Agency (JAXA), was also utilized. The data were downloaded from the National Aeronautics and Space Administration's (NASA) ASD Data Search Vertex web portal [73], with a 12.5 meter spatial resolution. To generate the Coverage and Land Use (LC/LU) maps, we used two images with a spatial resolution of 10 meters acquired on July 23rd, 2017 from the Sentinel 2A satellite, Path 17, and Row MRP and MRN. These were acquired from the European Space Agency's (ESA) Copernicus Services Data Hub platform, through QGIS's Semi-Automatic Classification Plugin (SCP) [74].

Micro-Watershed Delimitation
Delimitation of the HAL-MWs was done using the DEM and coded from the Utcubamba River HU N5. This process was based on the Pfafstetter method [75,76] while using PgHydro Tools, a QGIS (version 2.18.10) plugin [77], to activate the PgHydro Extension functions for PostgreSQL/PostGIS. The linear water network layer was imported from Google Earth Pro (version 7.3.2.5576) and SAS Planet (version 190707) interfaces, and subsequently updated and complemented with manual mapping [78,79]. This procedure was critical in obtaining the detailed geomorphology of the channel at the micro-watershed level, because the base layer brought the smoothed and generalized rivers (scale 1:100000).

Micro-Watershed Priorization Model
HAL-MW prioritization was carried out based on linear, areal, and shape morphometric variables using the PCWSA hybrid model proposed by Malik et al. [65]. The linear variables that were measured were: maximum and minimum height (Hmax, Hmin), area (A), perimeter (P), basin length (Lb), Strahler order (u), length (L), slope of the main stream (Sl), stream length (Lu), stream length mean (Lsm), and the Bifurcation ratio (Rb), which depends on Lu and the total number of streams of order u (Nu) ( Table 1). The areal variables were: the mean slope of the basin (Sb), drainage density (Dd), stream frequency (Fs), texture ratio (Rt), and mean length of the overland flow (Lom) ( Table 1). The analyzed shape variables were: the form factor (Ff ), Circularity ratio (Rc), Compactness coefficient (Cc), and Elongation ratio (Re) ( Table 1). 1 ∆H and ΣLl are the equidistance and the total length of the contour lines that pass through the basin, respectively.
Preliminary Priority Ranks (PPR) were assigned to each HAL-MW, with the use of one linear variable (Rb), four areal variables (Dd, Fs, Rt, Lom), and four shape variables (Ff, Rc, Cc, Re) [58,64,65]. Higher values of linear and area variables indicate a greater potential for soil erosion (direct relationship), while morphometric shape variables have an inverse relationship. Therefore, the highest erosion potential of these variables was assigned rank 1 (highest priority), the next highest potential value was assigned rank 2, and so forth for all HAL-MWs [58,61]. The correlation matrix, the first Load Factor (FL), and the rotated FL of the nine morphometric variables were constructed using PCA. This allowed us to identify the most significant morphometric variables. The PCA was performed using the SPSS 22.0 software; the methodological background for this can be found in Malik et al. [65]. The WSA was later applied to significant morphometric variables, and the value of the Composite Factor (CF) was calculated for the final priority classification. CF is defined by the PPR of the significant morphometric variable and its weight (Wi) (Equation (1)) [58,64,65]. Wi is obtained by analyzing the cross-correlation matrix between the significant morphometric variables (one per each component) and is calculated as the quotient between the vertical sum of the correlations for each variable (ri) and the total sum of correlations of the matrix (rij) (Equation (2)): The final priority range for the six HAL-MWs was assigned based on the value of CF. The lowest value was assigned priority rank 1, the next lowest value was assigned priority rank 2, and so forth for all HAL-MWs [65].

Land Cover and Land Use (LC/LU) Classification
To generate LC/LU maps, we followed the methodological flowchart developed by Rojas et al. [24]. All spectral bands were atmospherically and automatically calibrated by applying the Dark Object Subtraction (DOS1) [87] correction in the QGIS' SCP [74], and then bands 2-8, 11, and 12 were combined to construct multispectral images. These were adapted to the existing geographical boundaries of the study area and georeferenced using a second order polynomial transformation based on 33 Earth Control Points. Pixels were resampled to a new location by interpolation, with a permissible Mean Square Error (MSE) < 0.15 [88].
Based on the CORINE Land Cover methodology adapted for Peru [89] and prior knowledge of the study area, five classes of LC/LU were identified: built Area (BA), Andean grassland/scrubs (AG/S), grasses and crops (PC), water bodies (WB), and forest (Fo). Multispectral images were classified using the Maximum Probability algorithm based on the spectral signature of 218 training areas mapped in the field. Then, with the purpose of minimizing position and classification errors [90], the images were visually interpreted taking into account morphological characteristics such as shape, size, tone and color, patterns, texture, geographical position, and the association of the different LC/LU types [91]. Only the polygons where classification errors occurred due to the spectral similarity of the classes were modified [24].
The thematic accuracy of the maps was evaluated with the construction of a Confusion Matrix based on 196 verification sites [88]. These were established through a systematic randomized non-aligned stratified sampling on the final classified map [92] and verified in the field and in Google Earth Pro and SAS Planet interface [24]. The Global Accuracy and the Kappa Index (k) [93] were calculated.

Geomorphological Classification of Fluvial Systems
The classification of a territory's fluvial system is an important step in determining its ecological state [55]. The classification of fluvial systems is based on three geomorphological aspects: micro-watershed geomorphology (river style), channel slope, and valley geomorphology (Table 2) [53,54]. The classification process was carried out in two stages allowing fluvial system characterization from a geomorphological perspective: zoning and typification [53]. In regards to zoning, a sectorization of the river system was made for each geomorphological parameter, while with the intersection of these results, the typification stage was carried out. The latter stage involved categorizing, internally and functionally, homogenous types of river channel depending upon geomorphological aspects.
The linear channel layer (.shp) was interpolated with the DEM to acquire the altimetric data of the channel and was divided into sections of 1 km, which is the ideal observation scale for the use of IHG [44]. In the table of attributes, basic descriptors of each section were calculated ( Figure 3a): length (L), altitude, and the east and north coordinates of the initial node (Xi, Yi, Zi) and end (Xf, Yf, Zf). Then, over the center of each section, cross sections were drawn with an offset of 750 m on both sides, with the help of the QGIS's RiverGIS Plugin [94]. These were interpolated with the DEM and the length (Lcs) between the maximum altitude on the left (Zl) and right (Zr) of each cross section was calculated, as well as Zl, Zr, and the central altitude (Zc) corresponding to the channel axis ( Figure 3b). length (Lcs) between the maximum altitude on the left (Zl) and right (Zr) of each cross section was calculated, as well as Zl, Zr ,and the central altitude (Zc) corresponding to the channel axis ( Figure  3b).  All sections were classified as single channels (N1) ( Table 2), due to the infrequency of other channel variants (626 m of 286.50 km). Sinuosity index (S), slope (P), and valley confinement (E) were estimated with the use of Equations (3), (4), and (5) respectively. To calculate valley bottom width (V) through the use of Equation (6), the linear cross-sectional layer (.shp) was imported to the Google Earth Pro interface, where the ruler tool was used to measure said horizontal distance in accordance with each section's valley morphology as seen in Figure 3b. These data were manually registered to the table of attributes.  (6), the linear cross-sectional layer (.shp) was imported to the Google Earth Pro interface, where the ruler tool was used to measure said horizontal distance in accordance with each section's valley morphology as seen in Figure 3b. These data were manually registered to the table of attributes.
Each morphological parameter was reclassified according to established zoning criteria ( Table 2). For the typification process, the fields S, P, E, and V were concatenated for each section. Moreover, adjacent sections with the same classification were grouped. Each group, or Functional Sector (FS), is an internally and functionally homogenous fluvial channel.

Hydrogeomorphological Quality Evaluation
The hydrogeomorphological quality of each FS of the HAL-MWs was evaluated with the use of the IHG Index [19,46,96]. In the laboratory, hydrological and infrastructure documentation, satellite images (recent and old to assess the change processes; Figure 4), and cartography (terrain topography, land use and road network) were acquired to distinguish pressures and impacts on the river system that may distance its functionality, continuity, naturalness, complexity, and dynamics from a reference state [44]. A hydrogeomorphological process and impact cartographic base was generated [47]. In the field, the channel margins of almost all evaluated kilometers were traveled, during July and August 2017, to apply the final IHG through observations of the current state of the river system. This stage allowed the confirmation of observations in the laboratory, resolving doubts, looking for symptoms of impacts, finding others not visible in images or maps, and combining the expected effects of the these. FSs with difficult access were evaluated only in the laboratory.
Each morphological parameter was reclassified according to established zoning criteria ( Table  2). For the typification process, the fields S, P, E, and V were concatenated for each section. Moreover, adjacent sections with the same classification were grouped. Each group, or Functional Sector (FS), is an internally and functionally homogenous fluvial channel.

Hydrogeomorphological Quality Evaluation
The hydrogeomorphological quality of each FS of the HAL-MWs was evaluated with the use of the IHG Index [19,46,96]. In the laboratory, hydrological and infrastructure documentation, satellite images (recent and old to assess the change processes; Figure 4), and cartography (terrain topography, land use and road network) were acquired to distinguish pressures and impacts on the river system that may distance its functionality, continuity, naturalness, complexity, and dynamics from a reference state [44]. A hydrogeomorphological process and impact cartographic base was generated [47]. In the field, the channel margins of almost all evaluated kilometers were traveled, during July and August 2017, to apply the final IHG through observations of the current state of the river system. This stage allowed the confirmation of observations in the laboratory, resolving doubts, looking for symptoms of impacts, finding others not visible in images or maps, and combining the expected effects of the these. FSs with difficult access were evaluated only in the laboratory. The IHG was used to assess three sections of a fluvial system: (1) Functional Quality (FQ), (2) Channel Quality (CQ), and (3) Riparian Quality (RQ), with three subsections each. A score of 10 points was assigned to each subsection. However, these points were subtracted when impacts and damage were observed in accordance with every subsection criterion in the IHG Index. Final hydromorphological quality is calculated according to each FS' final score in conformity with Table 3 [46].

Morphometry and Preliminary Priority Ranges (PPR) of the Micro-Watersheds
The HAL-MWs are located at 2198-4275 m.a.s.l in Leimebamba and 1954-3790 m.a.s.l. in Molinopampa ( Figure 5). According to Strahler [80], the maximum stream order is three (3). The IHG was used to assess three sections of a fluvial system: (1) Functional Quality (FQ), (2) Channel Quality (CQ), and (3) Riparian Quality (RQ), with three subsections each. A score of 10 points was assigned to each subsection. However, these points were subtracted when impacts and damage were observed in accordance with every subsection criterion in the IHG Index. Final hydromorphological quality is calculated according to each FS' final score in conformity with Table  3 [46].

Morphometry and Preliminary Priority Ranges (PPR) of the Micro-Watersheds
The HAL-MWs are located at 2198-4275 m.a.s.l in Leimebamba and 1954-3790 m.a.s.l. in Molinopampa ( Figure 5). According to Strahler [80], the maximum stream order is three (3). The morphometric characterizations of HAL-MWs, based on linear variables, area, and shape are reported in Table 4. The results show Rb values ranging from 2.750 (Atuen HAL-WS) to 5.00 (Timbambo HAL-MW). Higher values of Rb indicate greater soil erosion [84]. Dd indicates the closeness of spacing of channels [61], which varies between 0.387 (Cabildo HAL-MW) and 0.698 (Timbambo HAL-MW). Low Dd values occur in regions with dense vegetation, low relief, and highly resistant and permeable subsoils, while high Dd are found in regions with sparse vegetation, high relief, and weak and impermeable subsoils [80]. Melton [97] analyzed the direct relationship between Dd, Fs, and the runoff processes. The high Fs value of the HAL-MW in Pomacochas (0.288) indicates there is more runoff in comparison to other HAL-MWs. Rt is classified into five kinds of texture, ranging from very thick (<2), thick (2-4), moderate (4-6), and good (6)(7)(8), to very fine (>8) [98]. All evaluated HAL-MWs show Rt < 2. Lom is an important independent variable that affects the hydrological and physiographic developments of the drainage basin [81].  After morphometric analysis, PPRs were assigned to all six HAL-MWs (according to the concept of direct and inverse relationships), as indicated in Table 5.

Micro-Watershed Prioritization using PCWSA
The positive correlations between the linear morphometric, areal, and shape variables are shown in Table 6. A strong correlation (r > 0.9) is observed between Dd-Lom, Ff -Re, and Rc-Cc, and moderate correlations (r > 0.60) exist between Rb-Ff, Fs-Dd, Fs-Lom, Fs-Ff, Fs-Rc, Fs-Cc, Fs-Re, Rt-Rc, and Rt-Cc. Table 7 indicates that the top three components have values >1.5, and together, these represent about 93.949% of the total variance. However, at this stage, it was too difficult to classify the variables into components and add physical significance [65].  Therefore, the first FL (not rotated) and the rotated FL were constructed using principal component analysis (Table 8). Due to the fact that the third component (PC-3) of the first FL is moderately correlated with Dd and Lom, it is difficult to obtain a physically important component [65]. However, after analyzing the rotated FL matrix, the most important morphometric variables were Lom (PC-3), Ff (PC-2), and Rc (PC-1). This is in contrast to research done by Malik et al. [65], who, while using the same nine morphometric variables, found that the most important variables for nine sub basins in the Bino basin of India were Fs, Lom, and Ff. Cross-correlation between the three significant morphometric variables (Lom, Ff, and Rc) is shown in Table 9, while CF values and the final priority range are indicated in Table 10. From the six evaluated HAL-MWs, the established priority for the San Antonio HAL-MW is 1 and the Pomacochas HAL-MW has a priority of 6.   [100]. After a supervised classification, visual interpretation allowed us to correct errors generated by the spectral similarity between the "GC" and "AG/S" types of grass. This is a direct consequence of the frequent transitions and overlapping of the types of grass [24].

Land Cover/Land Use (LC/LU)
Natural land coverage ("Fo" plus "AG/S") constitutes 62.3% to 84.4% of the territory of each HAL-MW, as well as 75.2% and 69.0% of the territory evaluated in Leimebamba and Molinopampa, respectively (Table 11). The largest areas of "BA" (0.77 km 2 ) and "GC" (61.15 km 2 ) are found in the Ventilla HAL-MW. Agricultural and livestock activities are present throughout all HAL-MWs and present the greatest anthropic pressure on soil and water. Oliva et al. [31] evaluated four different productive systems in Molinopampa: forest (PS1), open field pasture (PS2), a silvopasture system with Pinus patula (PS3), and another with Alnus acuminata (PS4). Of these, PS2 recorded the highest soil compaction value (395 psi), apparent density (0.93 g/cm 3 ), and EC (0.36 µS/cm), as well as the lowest values of phosphorus (4.22 ppm), organic carbon (3.64%), organic matter (5.92%), and nitrogen (0.31 ppm). According to both this research and another investigation [30], pH levels tend to decrease in the high Andean Pinus patula plantations of Amazonas because they are closely linked to plantation age and vegetation density. In the San Antonio HAL-WS, Oliva et al. [29] studied seven stages of migratory agriculture and observed that this process generates significant changes in physiochemical soil characteristics at different depths (greater impact at 0-15 cm than at 15-30 cm) due to forest-cutting and burning practices.   The "AG/S" type is the most representative (53.2%) amongst the HAL-MWs in Leimebamba,   The "AG/S" type is the most representative (53.2%) amongst the HAL-MWs in Leimebamba, followed by "GC" (24.5%) and "Fo" (22.0 %). Ramírez [67], Salas et al. [25], and Rojas et al. [24] report that "AG/S 's" natural meadows in Amazonas are used as natural open field pastures and managed by periodic anthropic burning. Furthermore, Vasquez et al. [28] recorded 129 weed species (out of 148 herbaceous plants) in these natural grasslands. The average abundance and the number of weed species in PS2 and silvopastures (PS3, PS4, and others) were 41.32% and 22.07%, and 111 and 70, respectively. Mendoza et al. [27] found that between 1989 and 2016, 32.02 km 2 of forest were lost in Leimebamba, at a rate of 1.19 km 2 /year, attributed to agricultural and livestock pasture expansion ("GC"). In the case of the Molinopampa HAL-MW, the "Fo" type is the most representative with 37.2%. In the specific case of the San Antonio HAL-MW, García-Pérez et al. [23] characterized the local homogeneous palm forest (genus Ceroxylon) and indicated that the low diversity of species (C. peruvianum, C. quindiuense, C. vogelianum, and C. parvifrons) is due to interaction with activities such as agriculture and livestock. However, Sanín [101] explains that a similar density of adult C. quindiuense in deforested grasslands and in forests may be due to (i) adult palm trees being saved from logging and (ii) regeneration through underground meristems after pasture installation.
Bacteriological parameters were analyzed in the most dynamic river areas of "GC" and "BC" types of HAL-MW in Ventilla (Figure 4b) by Chávez et al. [33], and indicated that this area is considerably contaminated by the presence of cattle near riverbanks and city wastewater discharge. Studies of macroinvertebrates, and the physicochemical and microbiological properties of water in other HAL-MWs of Amazonas, such as Alto Imaza [32], El Chido and Allpachaca-Lindapa [34], Chinata, and Gocta [36] and Shocol [35], show that quality decreases as a consequence of anthropic pressure in "GC", and "BA". Lastly, Ibisate et al. [47] indicate that the loss of hydrogeomorphological quality is closely linked to the sociodemographic pressure caused by the proliferation of artificial uses in the channel and the flood plain, but also by changes in basin land uses.

Fluvial Typology and Functional Sectors (FS)
Twenty-three types of river typology were identified in the Leimebamba HAL-MWs as a result of the combination of nine types of valley and six types of riverbed characterized in the fluvial typification stage (Figure 7a; Table 12). Twenty-eight river typologies were identified in Molinopampa from 10 types of valley and nine types of riverbed (Figure 7b; Table 13). In sum, 39 river typologies were identified from a combination of 13 types of valley and 11 types of riverbed. This figure reflects a very high landscape and river diversity.
Twenty-three types of river typology were identified in the Leimebamba HAL-MWs as a result of the combination of nine types of valley and six types of riverbed characterized in the fluvial typification stage (Figure 7a; Table 12). Twenty-eight river typologies were identified in Molinopampa from 10 types of valley and nine types of riverbed (Figure 7b; Table 13). In sum, 39 river typologies were identified from a combination of 13 types of valley and 11 types of riverbed. This figure reflects a very high landscape and river diversity.  (Table 2). Tables 12 and 13 indicate the occupancy percentages of each type of riverbed and valley. The dominant river typology in Leimebamba is a high slope straight channel located in a very tight valley and narrow river bottom (S1P4-E2V2) with 12.74%, followed by a moderately high sloped winding channel located in a gently fitted valley with a medium width river bottom (S2P3-E4V3; 8.99%), and finally by a high sloped winding channel located in a moderately fitted valley with a narrow width river bottom (S2P4-E3V2; 8.74%). In conclusion, sinuous riverbends with moderately high slopes (S2P3; 32.54%) and fitted valleys with wide width river bottoms (E2V2; 23.55%) are predominant.  (Table 2).   Tables 12 and 13 indicate the occupancy percentages of each type of riverbed and valley. The dominant river typology in Leimebamba is a high slope straight channel located in a very tight valley and narrow river bottom (S1P4-E2V2) with 12.74%, followed by a moderately high sloped winding channel located in a gently fitted valley with a medium width river bottom (S2P3-E4V3; 8.99%), and finally by a high sloped winding channel located in a moderately fitted valley with a narrow width river bottom (S2P4-E3V2; 8.74%). In conclusion, sinuous riverbends with moderately high slopes (S2P3; 32.54%) and fitted valleys with wide width river bottoms (E2V2; 23.55%) are predominant.
In Molinopampa, the dominant typology, with 17.81%, is rivers with sinuous high slopes located in a very tight valley with a narrow river bottom (S2P4-E2V2), followed by winding channel rivers with moderately high slopes in moderately fitted valleys and medium-width river bottoms (S2P3-E3V3; 12.80%), and finally by straight channel rivers with moderately-high slopes in very tight valleys with narrow river bottoms (S1P3-E2V2; 8.74%). In general, sinuous rivers with high slopes (S2P4; 33.70%) located in very embedded valleys with narrow river bottoms are predominant (E2V2; 37.43%).
Ollero [52] stated that river geodiversity is one of the planet's richest natural heritage features and therefore classifying channels and valleys is fundamental for any river study. Even though the most recent geomorphological methodology [53][54][55][56][57] is based on lithological units at a 1:50,000 scale and topographic units are generated with a 5 m spatial resolution DEM (Lithotopo) [57], this approach was not applied for the present study. Regardless of the fact that topographic units (altitude, slope, and roughness) can be generated from the most detailed DEM available in Peru (ALOS PALSAR, 12.5 m), the National Geological Cuadrangle lithology is not yet available at a 1:50,000 scale for the entire Peruvian territory. In addition, this computing resource is not useful for detailed local scale (micro-basins and stretches <1 km) work.
In Leimebamba and Molinopampa, 53 FS and 65 FS were assembled within 48.96 km and 65.80 km of total channel length, respectively (Table 14). These had variable lengths, from 0.30 km (FS4 in Cabildo HAL-MW) to 8.02 km (FS25 in Ventilla HAL-MW). Ollero et al. [44] mention that the smaller the FS (greater detailed work), the more accurate the resulting evaluation. However, they also state that the level of detail is conditioned by the study objective, or even by the budget itself. For example, Barboza et al. in the Utcubamba basin [37] and in the Leiva MW [38] considered eight FS and 17 FS in 250 km and 56.48 km of the total channel, respectively. 3.5. Hydromorphological Quality Determination using IHG Figure 8 shows the hydrogeomorphological quality pattern of the channels in all six HAL-MWs. In general, a deterioration in quality is observed as the altitude descends, from the high channels (tributaries of order 1) to the medium and low channels, except in the Ventilla HAL-MW. This pattern was found by Barboza et al. [37] in the Utcubamba basin.
In Leimebamba, the total IHG assessment showed that 8.5%, 26.6%, 30.5%, and 34.4% of the lengths of the channels are assessed as being of "Poor", "Moderate", "Good", and "Very good" quality, respectively (Table 15). In Molinopampa, 7.2%, 8.4%, 26.6%, and 57.8% of the lengths of the channels have "Poor", "Moderate", "Good", and "Very good" quality assessments. None of the sections of the six assessed HAL-MWs had a "Very bad" quality. Hence, sections that were assessed as "Good" and "Very good" can be considered as reference sections for river restoration [47]. Figure 8 shows the hydrogeomorphological quality pattern of the channels in all six HAL-MWs. In general, a deterioration in quality is observed as the altitude descends, from the high channels (tributaries of order 1) to the medium and low channels, except in the Ventilla HAL-MW. This pattern was found by Barboza et al. [37] in the Utcubamba basin.    Very bad  ----------------Poor  ---------------  In the high channels of Leimebamba and Molinopampa, the riverbank and floodplain degradation is predominantly caused by pressure from livestock and agricultural activities, such as grazing, migratory agriculture, clearing, fires, etc., that alter soil structure [18,22], induce shrubland growth due to the disconnection with the phreatic zone [10], stimulate weed proliferation [28], and produce longitudinal discontinuities [37]. The most important immediate impacts registered are those derived from vehicular and pedestrian bridges, small weirs, channels, and longitudinal stone defenses. Dams, canals, irrigation systems, and other hydraulic works cause deterioration in the river current, which often has irreversible and sometimes unknown consequences [11,48]. One exception is the Timbambo HAL-MW, where aggregate extraction is non-existent due to its road inaccessibility. This particular activity (aggregate extraction), when done massively and indiscriminately, causes lateral river incision and tends to affect riverbank and riverbed natural sediment accumulation areas [8]. Furthermore, the extensive access roads to crop parcels and access trails to basins for cattle break transversal riverbank connectivity are noteworthy [44]. The gravel surface of the Leimebamba-Chuquibamba road that runs parallel to, and sometimes intersects, the Cabildo and Atuen HAL-MWs, alters their hydromorphological processes causing overflows, floods, and flood flows [44]. The same impacts on San Antonio and Ventilla HAL-MWs are caused by the asphalt-covered road that connects Chachapoyas and Mendoza. Irrespective of the negative impacts just described, road infrastructure favors livestock growth and small-scale agricultural migration [24]. The case of the Pomacochas HAL-MW differs as it is located near the final road section and the district capital of Leimebamba (Figure 4a). Along this urban section, the floodplain has disappeared, the margins have been completely channeled, there are abundant obstacles, and the entire floodplain has been raised or waterproofed. In the Timbambo HAL-MW, 5.8% and 94.2% of the lengths of the channels have "Good" and "Very good" quality assessments, and this therefore has the largest number of sections that can be considered as reference conditions for river restoration [47].

Riparian quality
Finally, it can be concluded that hydrogeomorphological river alterations originate because of socio-economical activities that consume water, sediments ("aggregates"), and territory (river space), and also due to community preferences for living next to rivers in risky situations that require infrastructural protection against floods and river dynamics [8].

Morphometric Prioritization, Fluvial Classification, and Hydrogeomorphological Quality
Watershed prioritization examines the intensity of each HAL-MW's erosion problem so that the range will be used to prioritize the treatment of each with soil and water conservation measures [65]. River and valley geomorphological classification is fundamental for any kind of fluvial study. Hence, it must be used in any land planning project as an essential reference instrument for the understanding, functioning, and enhancement of natural systems [53]. Hydrogeomorphological dynamics guarantee the protection of each and every one of the elements of the river system, and it is the key not only to its functioning but also to its ecological, landscape, and environmental value [19]. Various studies of morphometric prioritization utilizing multivariate statistics [58,[61][62][63][64][65], geomorphological river classifications [53][54][55][56][57], and evaluations of hydrogeomorphological quality using IHG have been carried out [37,38,44,45,[47][48][49]. From the aforementioned, two did not integrate fluvial classification prior to the use of the IHG [48,49]. However, the integration of the three different methodologies used in this study [19,53,65] has not yet been reported.
This study evaluated the feasibility of an integrated approach to morphometric prioritization, geomorphological river classification, and the hydrogeomorphological assessment of micro-watersheds. The proposed methodology allows us to identify and prioritize micro-watersheds susceptible to erosion and riverbed sections in need of conservation and/or restoration to inform and improve the decision-making process in the selected study area. However, future research could incorporate ecological, socioeconomic, and geospatial perspectives to further enhance this kind of modeling.

Conclusions
The PCWSA hybrid model reveals that the San Antonio HAL-MW and the Atuen HAL-MW are located in areas very susceptible to erosion. Therefore, they were assigned top priorities, and rapid soil adaptation and water conservation practices are recommended. Consequently, the priority order is as follows: Ventilla HAL-MW, Cabildo HAL-MW, Timbambo HAL-MW, and Pomacochas HAL-MW. Concerning geomorphological river classification, 39 types of river course were identified within the six HAL-MWs, as the result of combining 13 types of valley and 11 types of riverbed, thereby indicating very high landscape and river diversity. Total IHG assessment gave results of 7.6% (21.8 km), 14.5% (41.6 km), 27.9% (80.0 km), and 50.0% (143.2 km) of the total channel lengths being of a "Poor", "Moderate", "Good", and "Very good" quality, respectively. None of the sections of the six assessed HAL-MWs had a "Very bad" quality assessment, hence sections that were assessed as "Good" and "Very good" can be considered as reference sections for river restoration. The loss of hydrogeomorphological quality is closely linked to the sociodemographic pressures caused by the rise in anthropic modifications of the basin and floodplain.
Given the importance of the HAL-MW ecosystems in Amazonas and in Peru, an integrated methodological framework of morphometric prioritization, geomorphological river classification, and hydrogeomorphological evaluation of hydrographic micro-watersheds is presented. Ergo, due to its low complexity, this methodology can be replicated, with the use of necessary complements, in all Peruvian ecosystems and will ultimately contribute to adequate territory, soil, and water planning and management.