Morphometric Evaluation and Its Incidence in the Mass Movements Present in the Chicamocha Canyon, Colombia

: The dynamic behavior of the basins evaluated by their morphometric parameters establishes a relationship with the endogenous and exogenous factors of the earth that control the modeling of the rocky massif by weathering and erosion processes. The characterization of these relationships can deﬁne the degree of affectation of the surfaces and the agents that control them as a categorical element in the deﬁnition of scenarios within the planning of physical and natural territory. This evaluation considers parameters contained within the characterization of the relief, shape and texture of the drainage and the mobility of the ﬂow in the basin. As a result of this evaluation, three factors were obtained that control the processes in the basins: a tectonic structural factor followed by climate control and, ﬁnally, seismic activity that plays an important role in the mechanical weathering of surfaces. The correlation of these factors showed that the Umpal á , Guaca, R í o Negro, Cantabara and La Cureña basins display a high degree of physical and mechanical weathering of the rocky massif and that the Manco, Talarcuta, Las Pavas, Felisco and El Abra basins are more likely to generate torrential ﬂows. The joint evaluation of these parameters provides a useful tool to understand the dynamic behavior of basins and their impact on anthropogenic setting.


Introduction
In recent decades, hydrometeorological conditions have generated a considerable increase in the manifestations of processes of instability on surfaces (mass movements).These conditions are possibly being caused by a significant climate change that is represented by an increase in average, maximum and minimum air temperatures that provide rainfall scenarios with greater intensity and recurrence, especially in mountainous areas where they can easily trigger mass movement events due to major changes in the ground surface [1][2][3][4][5].
At the global level, the effects generated by mass movements have caused drastic damage to infrastructure of different kinds, as well as victims and human losses.A study conducted by the World Health Organization between 1998 and 2017 showed that mass movements affected 4.8 million people and caused more than 18,000 deaths [6,7].In addition to the above, there is a substantial increase in the manifestations of events related to climate activity, especially rainfall, as well as anthropogenic elements that are affected by seismic activity [8].These events have generated considerable loss of human life and infrastructure, with an estimated damage reported by Swiss Re for 2011 of USD 370 billion in losses of anthropogenic elements and the natural environment [9].Within these climatic conditions, floods are one of the most influential factors worldwide, since, according to a study that examined the period between 1985 and 2009, floods represent 40% of all natural disasters [10].
In a more regional context, climatic conditions have triggered events such as floods and mass movements that have affected the infrastructure and generated loss of human life with incidents reported for the years 1998, 2010, 2011, 2013 and 2020 by the Inventory System of Disasters-DesInventar [11].Within the context of characterization of an environment that is disturbed by climatic and/or internal conditions of the earth, it is necessary to estimate the degree of affectation that may occur on anthropic and natural environments [12,13].
To establish these degrees of affectation, it is necessary to start from models defined by parametric and non-parametric functions on a fundamental unit in a regional analysis, and which in this case are based on the study of the hydrographic basin.This is because the basin describes this dynamic behavior in the same system containing patterns and/or behaviors typical of the surrounding atmosphere.Conditions such as floods, which occur with greater recurrence, together with the degradation of the rock massif, which may be susceptible to mass movements, respond to events of the torrential flows, landslides and falls that affect these environments.One of the conditions within the analysis of hydrographic basins is their morphometric characterization, which can include several factors that define the modeling of their surfaces.
Factors such as tectonic and structural condition and climatic behavior are closely related to the control of the shape, size, slope, texture of the surfaces of the hydrographic basins, as well as condition the processes of weathering and erosion of these surfaces [14,15].The morphometric analysis from a quantitative perspective determines the properties or punctual, linear, area and volumetric parameters of the relief [16,17].This evaluation allows to explain the geological and geomorphological history of each river basin from the quantitative perspective providing a description of the physiography and/or the behavior of the drainage system [18][19][20].The evaluation of the basin employs defined parameters for the characterization of the drainage system, the analysis of the mobility of the flow, the texture of the channel and characterization of the relief.The included parameters are the Size of the basin (A), Basin length (Lb), ratio of the Drainage bifurcation (Rb), Form factor (Rf), Drainage density (Dd), Drainage texture (T), Torrential coefficient (Ct), Basin concentration time (Tc), Asymmetry factor (Af), Hypsometric curve, Compactness coefficient (Kc), Massive coefficient (Cm), Orographic coefficient (Co), Sinuosity of mountain front (Smf), Ruggedness number (Rn), Stream length-gradient Index (SL) and Reason for valley width and valley altitude (Vf) [16,[21][22][23][24][25].
With the evaluation of each one of these parameters added to the analysis of the factors that condition the dynamic control of each hydrographic basin, studies focused on the dynamic behavior of surfaces and flows can be carried out.This type of analysis serves as a useful tool when evaluating the functioning of the hydrological system of a region and thus contributing to the management of natural resources [26][27][28], which, in turn, provides a basis for integration into territorial policies, complying with current regulations on issues related to territorial planning and disaster risk management, and providing an input from the quantitative and qualitative description of the basins, hydrographic for their morphometric parameters.
Under this approach, the present study focuses on an evaluation of surface dynamics that create unstable conditions on the slopes (mass movements) and qualitatively determines the ways in which these processes can affect anthropogenic elements.It is clear that the dynamics of the basins in relation to slope instability events are closely related to the degradation of the rock massif and the concentration of rainfall in each hydrographic basin.Therefore, complementary to the research approach, these morphometric parameters that are related to tectonic, structural, climatic, and seismic activity elements are quantified as precursors of weathering and erosion processes that trigger mass movement of different types such as, landslides, flow, and fall.

Regional Setting
The study area is located north of Bogotá, Colombia in a sector within the Chicamocha Canyon region, Santander.The total study area extends over 1357.62 km 2 and is subdivided by 16 hydrographic basins (Figure 1).

Regional Setting
The study area is located north of Bogotá, Colombia in a sector within the Chicamocha Canyon region, Santander.The total study area extends over 1357.62 km 2 and is subdivided by 16 hydrographic basins (Figure 1).The study area has an altitude of between 542 m to 4448 m with an average of 1743 m, with a slope of between 0° to 78°, an average of 28.24° and a domain in the range of between 22° and 27°, with surfaces that are mostly oriented to the west (247.5° to 292.5°) (Figure 2).The drainage pattern is defined as being mainly subparallel, parallel to subdendritic.The region is characterized by having a climate strongly controlled by the topography and the Intertropical Convergence Zone (ITCZ), where temperature and rainfall (annual average values) can range from 8° to 16° C and 2000 mm to 2500 mm in the upper regions of the basins, up to 20°-26° and 400-900 mm in the lower parts, respectively [29].The correlation of precipitation data with evapotranspiration values (500 to 1600 mm/year) [30], according to the Holdridge classification, and areas of dry to very dry tropical xerophytic sub-desert forests are described, characterized by cacti and shrubs with very small leaves to avoid evapotranspiration; they are located in the lower parts of the basins.In the high mountain areas (upper part), there are areas of humid forests [31].
These climatic conditions generate important and rooted changes in the surfaces (lower region of the basins), leading to the creation of canyons and slopes of moderate to high slope associated with dry and warm climates; on the other hand, in the upper parts, there are wetter areas associated with topographic expressions of plateaus that develop on the Santandereano Massif [32].Due to these climatic and topographic contrasts, the surfaces of these regions develop superficial, shallow, skeletal, and very eroded soils, where there is a predominance of rhosolic entisol-type soils [33].In addition, Rangel-Ch and Cadena [33] assert that depending on the topographic aspect, two large groups of soils can be differentiated: (1) Present on slopes with gentle to steep slopes, with poorly evolved soils, low content of organic matter, rate of erosion equal to or greater than soil The study area has an altitude of between 542 m to 4448 m with an average of 1743 m, with a slope of between 0 • to 78 • , an average of 28.24 • and a domain in the range of between 22 • and 27 • , with surfaces that are mostly oriented to the west (247.5 • to 292.5 • ) (Figure 2).The drainage pattern is defined as being mainly subparallel, parallel to subdendritic.The region is characterized by having a climate strongly controlled by the topography and the Intertropical Convergence Zone (ITCZ), where temperature and rainfall (annual average values) can range from 8 • to 16 • C and 2000 mm to 2500 mm in the upper regions of the basins, up to 20-26 • and 400-900 mm in the lower parts, respectively [29].The correlation of precipitation data with evapotranspiration values (500 to 1600 mm/year) [30], according to the Holdridge classification, and areas of dry to very dry tropical xerophytic sub-desert forests are described, characterized by cacti and shrubs with very small leaves to avoid evapotranspiration; they are located in the lower parts of the basins.In the high mountain areas (upper part), there are areas of humid forests [31].
These climatic conditions generate important and rooted changes in the surfaces (lower region of the basins), leading to the creation of canyons and slopes of moderate to high slope associated with dry and warm climates; on the other hand, in the upper parts, there are wetter areas associated with topographic expressions of plateaus that develop on the Santandereano Massif [32].Due to these climatic and topographic contrasts, the surfaces of these regions develop superficial, shallow, skeletal, and very eroded soils, where there is a predominance of rhosolic entisol-type soils [33].In addition, Rangel-Ch and Cadena [33] assert that depending on the topographic aspect, two large groups of soils can be differentiated: (1) Present on slopes with gentle to steep slopes, with poorly evolved soils, low content of organic matter, rate of erosion equal to or greater than soil formation and where leaching processes are frequent.(2) At the edges of tributaries in sites with soils of medium content of organic matter, poor in soluble phosphorus and where the main agricultural work of the region has been developed.
formation and where leaching processes are frequent.(2) At the edges of tributaries in sites with soils of medium content of organic matter, poor in soluble phosphorus and where the main agricultural work of the region has been developed.
In the Cretaceous period (K), rocks are generated with different compositions such as clear quartz sandstones, with intercalations of gray to brown siltstones, white to gray sandstones; intercalations of black shale, locally sandy; black clays, dark gray limestones; quartzite white sandstones, quartzite white sandstones; laminated black clays, fossiliferous limestones; dark gray shale; clay limestone and chert; greenish grey shale, limonitic nodules [35][36][37].In the Quaternary period, materials are deposited that are the product of the processes of erosion and weathering of the surrounding massifs; they are distinguished by their transport mechanism that can be of fluvial or hydrogravitational type.Among these deposits are terraces and cones of dejection (Qal) with a lithology of rounded to subrounded clasts, of size between coarse sands to pebble, in a clayey sand matrix, as well as deposits of fluvial transport of size between fine sand to gravel, rounded to subrounded; colluvial, slopes, and collapse (Qc), with rock fragments of size between pebble and boulder, angular to subangular [34][35][36][37].
In the Cretaceous period (K), rocks are generated with different compositions such as clear quartz sandstones, with intercalations of gray to brown siltstones, white to gray sandstones; intercalations of black shale, locally sandy; black clays, dark gray limestones; quartzite white sandstones, quartzite white sandstones; laminated black clays, fossiliferous limestones; dark gray shale; clay limestone and chert; greenish grey shale, limonitic nodules [35][36][37].In the Quaternary period, materials are deposited that are the product of the processes of erosion and weathering of the surrounding massifs; they are distinguished by their transport mechanism that can be of fluvial or hydrogravitational type.Among these deposits are terraces and cones of dejection (Qal) with a lithology of rounded to subrounded clasts, of size between coarse sands to pebble, in a clayey sand matrix, as well as deposits of fluvial transport of size between fine sand to gravel, rounded to subrounded; colluvial, slopes, and collapse (Qc), with rock fragments of size between pebble and boulder, angular to subangular [34][35][36][37].The structural arrangement of the region is framed mainly by the layout of the Bucaramanga Fault which is a regional type structure with a combined component in the course (sinistral) and in the dip (reverse) patterns [39,40].On the other hand, there are Manco River Faults, Umpalá River and Perchiquez River, with a northeast direction that extends from the Bucaramanga fault, on igneous and metamorphic rocks.These faults, when they cut the Bucaramanga fault, produce right lateral displacements [36].Baraya Fault presents a northerly direction; with a low-angle dip (locally), it contactspre-Cretaceous rocks with Cretaceous.Finally, there is Morro de las Peñas Fault that presents similar kinematics to Baraya Fault, but the latter is a little further north than the previous one [36].
For the region, a geomorphological control of a structural type is presented where structurally controlled loin units, flexural hooks and controlled basins are developed.These units, such as the loin, generate large block slide as a back-tilting on the left surface of Chicamocha River in response to the decompensation of the balance profile of the slope, causing its staggered appearance [41].In addition, by the interaction of the faults present, The structural arrangement of the region is framed mainly by the layout of the Bucaramanga Fault which is a regional type structure with a combined component in the course (sinistral) and in the dip (reverse) patterns [39,40].On the other hand, there are Manco River Faults, Umpalá River and Perchiquez River, with a northeast direction that extends from the Bucaramanga fault, on igneous and metamorphic rocks.These faults, when they cut the Bucaramanga fault, produce right lateral displacements [36].Baraya Fault presents a northerly direction; with a low-angle dip (locally), it contactspre-Cretaceous rocks with Cretaceous.Finally, there is Morro de las Peñas Fault that presents similar kinematics to Baraya Fault, but the latter is a little further north than the previous one [36].
For the region, a geomorphological control of a structural type is presented where structurally controlled loin units, flexural hooks and controlled basins are developed.These units, such as the loin, generate large block slide as a back-tilting on the left surface of Chicamocha River in response to the decompensation of the balance profile of the slope, causing its staggered appearance [41].In addition, by the interaction of the faults present, fragile deformations are generated on the rocks surrounding the trace of the faults, which give rise to multiple fractures of a systematic and penetrative nature that originate a centimeter of tabular cataclastic foliation.With this relationship added to the climatic conditions, a large amount of particulate material develops, giving rise to deposits of types of dejection cones, alluvial fans, and terraces [41].

Materials and Methods
To address the morphometric characterization of the hydrographic basins from a quantitative analysis, the elements defined within the physiographic characteristics of the environment are evaluated based on a parametric analysis of the relief, the drainage system, and the texture of the channel.The evaluation of each one of these elements describes the degree of affectation exerted by endogenous and exogenous factors of the earth on them and which, in turn, control the processes of weathering and erosion of the surfaces.From this analysis, the dynamic behavior of each hydrographic basin and its response in qualitative terms of the degree of susceptibility to generating mass movements of the type of landslide, flow, and fall are materialized.For the evaluation of the morphometric conditions present in each of the hydrographic basins, equations arranged internationally for this type of characterizations were used (Table 1).The base cartography was obtained from the free data of the Agustín Codazzi Geographic Institute (IGAC) at a scale of 1:25,000 [42]; in turn, the Digital Elevation Model (DEM) was implemented with a resolution of 12.5 m [43], which was obtained from the portal ASF-ALOS PALSAR.These inputs within the morphometric characterization of the basins are the basis for scale of analysis (1:25,000) and the degree of resolution (12.5 m) of the present study.From the base cartography, the configuration of the drainage network and the altitude levels was obtained, and the limit of the hydrographic basins was defined, as well as other morphometric parameters that will be described below.The concatenation of these two inputs served as the basis for the calculation of each morphometric parameter of the hydrological basins.The morphometric calculations performed were developed on the ArcGIS V10.8 program.The Form factor (Rf) expresses the relationship between the area of the basin (km 2 ) and the length of the basin (km) [22].The form factor described by Pérez [55] was used for the evaluation of this parameter.The Drainage density (Dd) parameter indicates the relationship between the total length of the irregular and regular watercourses of the basin (km) and the total area of the basin (km 2 ) [22].This parameter provides a real measure of drainage efficiency when examining the capacity to dislodge a volume of water.The Drainage Texture (T) parameter shows a relationship between Drainage density (Dd) by Stream frequency (Fs) [23].This parameter is conditioned by factors such as climate (temperature and precipitation), vegetation, type of rock and soil, infiltration capacity, relief, and its state of development.
The Basin concentration time (Tc) evaluates the hydraulic behavior of each basin against each morphometric and morphological aspect.The concentration time is the travel time of a drop of rainwater that drains superficially from the farthest place in the basin to the point of exit [56].Different formulas were used to calculate this parameter, which are described below (a-i): (a) SCS-Ranser method: T c = 0.97K 0.385 , K = L 3 c H , where T c is the concentration time (hours), H is the difference between the highest and lowest elevation of the basin (feet) and L c is the length of the main channel (km) [57].
(b) California Culvert Practice method: T c = 0.87 , where T c is the concentration time (hours), H is the difference between the highest and the lowest level of the basin (m) and L c is the length of the main channel (km) [58].
(c) Kirpich's method: T c = 0.066 , where T c is the concentration time (hours), L is the length of the main channel to the divide (km) and S 0 is the average slope of the main channel (m/m) [59].
, where T c is the concentration time (hours), L c is the length of the main channel (km) and S 0 is the average slope of the main channel (in percentage) [60].
(e) Giandotti's method: T c = 4 , where T c is the concentration time (hours), L c is the length of the main channel (km), S 0 is the average slope of the main channel (m/m) and A is the area of the basin (km 2 ) [61].
(f) V.T Chow's method: T c = 0.8773 , where T c is the concentration time (hours), L c is the length of the main channel (km), CM is the highest bound of the main channel (m.a.s.l.) and Cm c is the lowest bound of the main channel (m.a.s.l.) [56]., where T c is the concentration time (hours), L c is the length of the main channel (km) and S 0 is the average slope of the main channel (in percentage) [63].
(i) Passini's method: T c = 0.108 , where T c is the concentration time (hours), L c is the length of the main channel (km), S 0 is the average slope of the main channel (m/m) and A is the area of the basin (km 2 ) [64].
The parameter Asymmetry factor (Af) describes tectonic tilting conditions perpendicular to the main drainage direction.It is defined by the basin area to the right of the main drainage downstream (Ar) and the total basin area (A) [48].Depending on its result, a tilt is defined to the left (>50) or to the right (<50).The Compactness coefficient (Kc) compares the shape of the basin with that of a circumference whose inscribed circle has the same area of the basin under study.It is defined with the perimeter of the basin and the area of the basin.The degree of approximation of this index to the unit will indicate the tendency to concentrate strong volumes of runoff water, being more accentuated the closer it is to the index, that is, the greater the concentration of water [21].The Massive coefficient (Cm) calculates the condition of the basin with respect to whether it is mountainous or, on the contrary, a flat basin [49].
The Orographic coefficient (Co) parameter expresses the degradation potential of the basin and considers the average height of the basin which directly influences the potential flow of water and the area whose inclination exerts direct action on surface runoff [50].The result expresses a little rugged relief for <6 values, as well as extensive basin, low slope, and lower erosion, and for >6 values, it expresses a rugged relief with a high degree of erosion [50].The Elongation relation (Re) is described as the quotient between the diameter of a circle with the same basin surface area and the maximum length of the basin, defined as the largest dimension of the basin, along a straight line from the outlet to the extreme boundary of the watershed, parallel to the main river [24].The Circularity ratio (Rc) shows the relationship between the basin area and the circle area with the same circumference as the basin perimeter.This relationship is influenced by current length and frequency, geological structure, land use and/or land cover, climate, and slope [51].
The Sinuosity of mountain front (Smf) parameter is defined as the ratio or relationship between the total length of the mountain front (Lmf) and the straight-line length measured between the ends of the mountain front (Ls) [25].The Index of the gradient of the longitudinal profile of the riverbed (SL) is used to detect the anomalous gradients within the path of the currents in a hydrographic basin from the morphometric point of view [65].The drainage network presents regular geometric properties which establish erosive processes that help provide data on the geological conditions of the evaluated area as well as its geomorphological evolution [66].The parameter of Reason for valley width and valley height (Vf), according to Bull and McFadden, (1977); Bull, (1978); cited in Keller and Pinter [48], allows to evaluate the distribution and degree of tectonic activity specifically on the mountain range-piedmont junction.This parameter differentiates between broad-based valleys with relatively high values that indicate degradation processes and low rates of uplift, and valleys in the form of "V" with relatively low values that suggest prolonged incision and active lifting.The equation defines the Reason for valley width and valley height (Vf), where the average height of the slope on the left side of the valley is Eld, the average height of the left side is Erd, the measure of elevation of the channel is Esc, and the width of the valley is Vfw.

Results
From the morphometric analysis carried out on each hydrographic basin following the proposed methodology and the equations described above, each morphometric aspect was determined as a regional evaluation.The following results are obtained from this analysis.
The Form factor parameter (Rf) described by Horton [22] determined that the hydrographic basins of Río Negro, Río Guaca, river and Seca, Felisco and El Abra stream have a slightly elongated geometry (northeast region) providing conditions of moderate to low occurrence to sudden floods of the channels.On the other hand, the basin of the Cantabara stream (widened geometry) has a greater tendency to generate a sudden flood of the channel (Table 2).
The Bifurcatio ratio (Rb) parameter describes the change ratio of each drainage order [16].This relationship shows that the Calichana stream-Seca stream (Rb, orders 1 and 2), Las Pavas stream-Perchiquez stream (Rb, orders 2 and 3), Chorrerana stream-Perchiquez stream (Rb, order 3 and 4) and Manco River-Umpalá River (Rb, order 4 and 5) present the highest Rb coefficients due to an orographic factor that is controlled by the geological conditions associated with the stroke of the Bucaramanga fault line, Perchiquez and Umpalá faults that affect rocks especially of Precambrian (PD) and Jurassic (J) age causing their steep terrains that generate very erodible soils (Table 2).
Performing a comprehensive analysis of the parameters calculated in Table 2 for the description of the drainage network, it is determined that the basins of the Guaca, Umpalá and Manco rivers together with the Cantabara, La Cureña and Chorrerana streams present scenarios conducive to the creation of torrential flows, as well as weathering processes and erosion of their surfaces, giving rise to a high concentration of stone material especially of rocks of Precambrian (PD) and Jurassic (J) origin.
Within the parameters for the characterization of the relief of the hydrographic basins, the Asymmetry factor (Af) describes a tendency of the basins to tilt to the right with a marked asymmetry.The basins with the highest Af coefficients are Umpalá, Quebrada Seca and Quebrada Talarcuta.This tilting factor is controlled by tectonic and structural conditions, especially by the stroke of the Bucaramanga fault, that present kinematics of sinistral course with a reverse component (Figure 4 and Table 3).The Compactness coefficient (Kc) signifies that the Cantabara, Felisco streams and the Río Negro River have a round-oval to oblong-oval basin shapes with indices closer to 1, indicating the tendency to concentrate strong volumes of water by runoff from their surfaces.The Massive coefficient (Cm) describes a tendency of the basins towards a mountainous to very mountainous relief, the streams El Abra, Seca and Las Pavas being the surfaces that present the highest coefficients in this parameter (Table 3).
The parameter of Orographic coefficient (Co) presents a similar behavior for all basins, with little rugged reliefs, extensive basins, low slope, and less erosion (Table 3).The elongation ratio (Re) describes a tendency towards elongated to very elongated basins, with the Manco River and the Talarcuta and Perchiquez streams having the lowest indices (Table 3).In turn, this parameter indicates a delay in the concentration of runoff due to the length of the channel, which generates a high probability of torrential flow when the channel reaches its maximum concentration level [67].The Circularity ratio (Rc) also defines the ability of the basin to evacuate runoff water.This parameter is based on the area of the basin.The coefficients obtained show a tendency towards intermediate to slightly elongated basins.The basins with lowest coefficients are the Manco River and the Perchiquez and Talarcuta streams, which will have a delay in runoff concentration due to the length of the channel (Table 3).For the parameter of Sinuosity of mountain front (Smf) it is observed that a large part of the basins is in an active state, except for the basins of the La Cureña, Seca, Cantabara and Calichana streams (little activity).The basins of the Chorrerana, Perchiquez and Talarcuta streams show very low coefficients that indicate greater activity of the mountain front in weathering processes and erosion of the rocky massif (Figure 5 and Table 3).
The hypsometric curve is one of the most important parameters in the morphometric assessment of the basins due to its relationship with the state of erosion of the channel.To carry out the analysis of all the basins, a normalization of the data was carried out, superimposing the trends of all the basins.Within the construction of the curves and their interpretation, the mature state (red line) was considered as the basis for the graphic output (Figure 6).
With the construction and evaluation of the hypsometric curves, three conditions are described: (1) basins in juvenile state, where there is great potential for erosion (bottom undermining of the bed), (2) basins in mature state, where there is balance (bottom and lateral undermining), and (3) basins in senile state, where lateral deposition and undermining conditions occur [45].With these conditions, it is observed that the basins present For the parameter of Sinuosity of mountain front (Smf) it is observed that a large part of the basins is in an active state, except for the basins of the La Cureña, Seca, Cantabara and Calichana streams (little activity).The basins of the Chorrerana, Perchiquez and Talarcuta streams show very low coefficients that indicate greater activity of the mountain front in weathering processes and erosion of the rocky massif (Figure 5 and Table 3).
The hypsometric curve is one of the most important parameters in the morphometric assessment of the basins due to its relationship with the state of erosion of the channel.To carry out the analysis of all the basins, a normalization of the data was carried out, superimposing the trends of all the basins.Within the construction of the curves and their interpretation, the mature state (red line) was considered as the basis for the graphic output (Figure 6).
With the construction and evaluation of the hypsometric curves, three conditions are described: (1) basins in juvenile state, where there is great potential for erosion (bottom undermining of the bed), (2) basins in mature state, where there is balance (bottom and lateral undermining), and (3) basins in senile state, where lateral deposition and undermining conditions occur [45].With these conditions, it is observed that the basins present in     The Index parameter of the gradient of the longitudinal profile of the river channel (SL) or Hack Index parameter quantifies the geomorphological conditions within the causes (peaks in erosion dynamics) that have been affected by tectonic activity, differential erosion of rocks and slope instability processes.Its implementation in this type of studies helps to detect knickzones which are the turning points on the longitudinal profile of the current (typical shape: concave upwards) [54].As a result of this analysis, active processes of the slopes (landslides) that can interact directly with the riverbed can be described, as well as the presence of geological structures and affectation of the shape of the channels [54] (Figure 7).The Index parameter of the gradient of the longitudinal profile of the river channel (SL) or Hack Index parameter quantifies the geomorphological conditions within the causes (peaks in erosion dynamics) that have been affected by tectonic activity, differential erosion of rocks and slope instability processes.Its implementation in this type of studies helps to detect knickzones which are the turning points on the longitudinal profile of the current (typical shape: concave upwards) [54].As a result of this analysis, active processes of the slopes (landslides) that can interact directly with the riverbed can be described, as well as the presence of geological structures and affectation of the shape of the channels [54] (Figure 7).Considering the above, it is determined that 1.
The Manco River has the largest number of knickzones towards the upper part of the basin.It is related to a slope change and altitudinal difference.Due to high concentration of knickzones in this part of the basin, it is more susceptible to generate processes of erosion of the channel that can trigger events of the type of torrential flows.

2.
Umpalá River presents on its tributaries on the left bank (downstream) a high concentration of knickzones.This may be related to the stroke of Umpalá Fault and the tilting present in the basin.

3.
Guaca River presents a condition like that of Umpalá River, with the exception that the knickzones are concentrated in the upper parts of each tributary.4.
The Rio Negro basin presents knickzones in several segments of the river, upper part of the basin, tributaries of the left surface (downstream) and the lower zone of the basin.It can be mainly related to the stroke of the Bucaramanga fault, in addition to multiple structures and lithologies present in the area.

5.
The Perchiquez stream basin has the same condition as that of the Umpalá, and the Talarcuta stream basin has a higher concentration of knickzones in the region of change in course of drainage.This change in direction and concentration of knickzones is associated with the stroke and kinematics of the Bucaramanga Fault.

6.
The other basins present on the right bank of Chicamocha River (downstream), such as the Seca, Chorrerana, Abra and Pavas streams, present a greater number of knickzones in the middle part of the basins.This is associated with a structural behavior (mainly the Bucaramanga fault), altitudinal difference and lithological change.7.
The basins of the Calichana, Cantabara, Felisco, Cureña, Lauchas and Basteña streams (left bank of the Chicamocha River) have the highest concentration of knickzones in the middle part of each basin.These knickzones are associated with the structural behavior of the Bucaramanga fault, together with an altitudinal difference.
The conditions evaluated by this morphometric parameter show areas of sectorized concentration described above, which indicates the regions highly susceptible to erosion processes (lateral and vertical) that are controlled by structural aspects, exerted mainly by the Bucaramanga fault stroke, as well as lithological aspects, since almost 50% of the rocks are of igneous and metamorphic origin, causing rigidity of the rocky massif.Added to these conditions are the physical characteristics of the environment described in Tables 2 and 3, as well as the tilting of the basins as currently reflected on the shape of the landscape.In some sectors, torrential flows can be generated by the combination of two or more conditions described.
Within the evaluation of the characterization of the relief of the hydrographic basins, the parameter of Drainage density (Dd), Horton [22] defines the efficiency of drainage on the basis of whether a basin can dislodge a volume of water in the shortest possible time.Within this parameter it is described that the basins of the Río Negro, Río Guaca, Río Umpalá, the Chorrerana, La Cureña and Perchiquez streams present a moderate to low drainage density, which indicates a greater time required to dislodge the volume of water contained in the hydrographic basin (Table 4).The Drainage texture (T) defines the level of runoff that a basin can have [23].From this parameter, it was determined that 70% of the basins (11 basins) present a very fine drainage texture, indicating basins with high levels of surface runoff, impermeable rock massif and low permeability soils.The Río Negro, Río Guaca and Quebrada Chorrerana basins present an intermediate texture, characteristic for areas with medium runoff levels, with an intermediate texture and permeability of soils (Table 4).For parameter of Reason for valley width and valley height (Vf), it was observed that the basins of the Manco, Río Negro and Umpalá, and the Calichana, Cantabara, La Cureña, Perchiquez and Talarcuta streams in some sectors of the channel present strong processes of incision and lifting with active fronts.The other basins including those previously named display degradation conditions, with expansion of the bottom and inactive front.The Ruggedness number (Rn) parameter indicates the structural complexity of the terrain associated with relief and drainage density [52].The result of Rn indicates that the Manco, Talarcuta Umpalá, and Guaca basins have a high structural complexity and are very susceptible to high erosion of their surfaces.On the contrary, the Chorrerana, Seca, Cantabara, and La Cureña basins have the lowest coefficients.Finally, the Melton ruggedness number (MRn) shows the index related to the accumulation and resistance of the simple flow of a basin, which, in turn, describes the susceptibility to generate different types of flow [53].From the result of MRn, it is interpreted that the basins La Basteña, Río Negro, Las Lauchas, Manco, Umpalá, Cantabara, and Guaca are prone to generate flows with low material content.The Seca, Talarcuta, Calichana, Felisco, La Cureña, and Perchiquez basins are prone to generate flows with moderate material content and the El Abra, Chorrerana, and Las Pavas basins are prone to generate flows with high material content.
With the results obtained from Tables 3 and 4, an integral analysis of these twelve parameters described for the characterization of the relief was carried out.From this analysis, it was determined that the Manco, Umpalá, Río Negro River basins and the Talarcuta and Perchiquez streams present the most unique conditions within each evaluation carried out.From the calculation of each parameter, a control is observed by processes of tectonic, structural, lithological and topographic type that provide their own and accentuated conditions in sectors within each of these hydrographic basins.The combination of each of these elements magnifies the processes of erosion of the rock mass and of the channels, giving rise to processes of lateral and vertical erosion of the channel that, together with the speed in the accumulation of the flow, result in a favorable scenario for torrential flow type events.
Within the parameters for the analysis of the flow mobility of the hydrographic basins, the Torrential coefficient (Ct) was evaluated.From this parameter, it can be seen that the basins of the Las Pavas, Felisco, El Abra and Seca streams have the highest torrential indices.On the other hand, the basins of the Guaca, Río Negro rivers and the Chorrerana and La Cureña streams present the lowest indices.Finally, the Basin concentration time (Tc) of each basin was estimated.From the results obtained, it can be seen that the basins of the Guaca (154.46 min), Manco (104.52 min), Umpalá (103.54 min) and Perchiquez streams (75.17 min) present the highest concentration times (the values present here correspond to the average of the results for each method).The basins of the Seca (15.80 min), El Abra (16.65 min), Las Pavas (25.16 min) and Calichana (26.75 min) streams present the lowest concentration times.At the same time, it is observed that the Manco River basin has an irregular dispersion of the data.This is associated with the geometry of the basin and its structural control (Table 4 and Figure 8a).For Figure 8b, there is a correlation between the concentration times of the basins and the torrential coefficient.From this correlation, it is observed that the basins of the Seca, Felisco, El Abra and Las Pavas streams show low concentration times but the highest rates of torrentiality.This is because they are smaller basins that are in an equilibrium phase (juvenile-mature) where the development of their surfaces generates both vertical and horizontal erosion processes; in turn, a structural control is included that is exerted on these basins, causing significant changes in their surfaces since they are located on the stroke of the Bucaramanga fault.

Discussion
With the results obtained in each of the parameters raised, very relevant aspects were detected when performing an evaluation of the geodynamic behavior present for the study region from the point of view of weathering and erosion processes.Within the analyses, it was observed that the morphometric conditions (geometry and surface contrasts) would be strongly controlled by a structural context as described in Tables 3 and 4. The action of regional and local faults such as Bucaramanga, Manco, Umpalá, Perchiquez, Baraya and Morro de las Peña faults create important changes in the surfaces, such as tilting, strong ground incision processes, control in the weathering pattern-erosion, changes in the geomechanical behavior of the rocks and sediment transport.These conditions are represented by the calculations performed with the parameters of the characterization of the relief (Table 3), where basins are described with accentuated processes in the mountain fronts, deepening of the channel due to a phase of dynamic balance of the channels, elongated geometries, tilted and with flexures in the direction of the channels.All these elements together offer a support in the morphometric construction of the landscape that derives in actions at a regional or local level, where the results of the present evaluation serve as a basis input in the ordering of the hydrographic basins.In turn, the descriptions provided by each morphometric parameter define sectors with a greater or lesser degree of affectation by processes that degrade the rocky massif, which provides useful tools in the planning and organization of both natural and physical resources of a territory.
To further understand the development of these basins, regional behavior must be integrated from the tectonic and lithological point of view.The region is located in the tectonic context of the Santander Massif, which is a structure within the eastern mountain range of Colombia.The Santander Massif presents conditions of cortical deformation due to lifting processes, transpressive and elapsed events and compressive deformation with folding in sedimentary rocks [39,41].These aspects create zones with different degrees of deformation within the massif, which give rise to metamorphic rocks of the mylonite augen, mylonite with non-deformed blocks, mylonite with banding and igneous and sedimentary rock with catalysis types [41].Due to these tectonic and structural conditions, previous studies in the region have identified sections within the Bucaramanga fault that present different levels of stresses and deformation (Villamizar, (2014) and González and Jiménez, (2014) cited by Velandia and Bermúdez, [39]; García-Delgado, et al., [68]).The study region is part of the sector called the Cepita section, which presents a high structural complexity due to the stroke of the Bucaramanga fault that exhumes part of the fragileductile crust zone.This description was constructed on the basis of outcrops with identified structures, in addition to the presence of pseudotachylites that served as indicators of this process (Villamizar, (2014) and González and Jiménez, (2014) cited by Velandia and

Discussion
With the results obtained in each of the parameters raised, very relevant aspects were detected when performing an evaluation of the geodynamic behavior present for the study region from the point of view of weathering and erosion processes.Within the analyses, it was observed that the morphometric conditions (geometry and surface contrasts) would be strongly controlled by a structural context as described in Tables 3 and 4. The action of regional and local faults such as Bucaramanga, Manco, Umpalá, Perchiquez, Baraya and Morro de las Peña faults create important changes in the surfaces, such as tilting, strong ground incision processes, control in the weathering pattern-erosion, changes in the geomechanical behavior of the rocks and sediment transport.These conditions are represented by the calculations performed with the parameters of the characterization of the relief (Table 3), where basins are described with accentuated processes in the mountain fronts, deepening of the channel due to a phase of dynamic balance of the channels, elongated geometries, tilted and with flexures in the direction of the channels.All these elements together offer a support in the morphometric construction of the landscape that derives in actions at a regional or local level, where the results of the present evaluation serve as a basis input in the ordering of the hydrographic basins.In turn, the descriptions provided by each morphometric parameter define sectors with a greater or lesser degree of affectation by processes that degrade the rocky massif, which provides useful tools in the planning and organization of both natural and physical resources of a territory.
To further understand the development of these basins, regional behavior must be integrated from the tectonic and lithological point of view.The region is located in the tectonic context of the Santander Massif, which is a structure within the eastern mountain range of Colombia.The Santander Massif presents conditions of cortical deformation due to lifting processes, transpressive and elapsed events and compressive deformation with folding in sedimentary rocks [39,41].These aspects create zones with different degrees of deformation within the massif, which give rise to metamorphic rocks of the mylonite augen, mylonite with non-deformed blocks, mylonite with banding and igneous and sedimentary rock with catalysis types [41].Due to these tectonic and structural conditions, previous studies in the region have identified sections within the Bucaramanga fault that present different levels of stresses and deformation (Villamizar, (2014) and González and Jiménez, (2014) cited by Velandia and Bermúdez, [39]; García-Delgado, et al., [68]).The study region is part of the sector called the Cepita section, which presents a high structural complexity due to the stroke of the Bucaramanga fault that exhumes part of the fragileductile crust zone.This description was constructed on the basis of outcrops with identified structures, in addition to the presence of pseudotachylites that served as indicators of this process (Villamizar, (2014) and González and Jiménez, (2014) cited by Velandia and Bermúdez, [39]).In turn, geomorphological ranges of recent tectonic activity were detected, such as triangular and trapezoidal facets, sealing ridges, pressure peaks, drainage deviation, linear valleys, and drainage control [41].
With this tectonic and structural perspective present in the study region, contrasts are established with the results of the calculated parameters.Parameters such as the Asymmetry factor (Fa), Hypsometric curve, Massive coefficient (Cm), Orographic coefficient (Co), Sinuosity of mountain front (Smf) and Reason for valley width and valley height (Vf) are closely related to recent tectonic activity.This relationship shows a tendency of tilting of the basins in the direction of stroke and kinematics of the Bucaramanga fault, as expected, except for the basins of Manco River and the Calichana, Las Lauchas and El Abra streams that may be due to a local structural control, a compositional variation of the lithology present or an increase in the rates of weathering and erosion on the surfaces due to zonal climate control.In addition, a relationship is observed between the relief (mountainous to very mountainous) that is a little rugged, with extensive basins, strong processes of incision, lifting, degradation of the slopes, undermining (vertical and horizontal), expansion of the bottom and active fronts that can obey this structural control.In turn, the exhumation rates were recorded by Velandia [69] and Amaya et al. [70] of approximately 0.3 km/Ma from approximately 25 Ma ago.As a parallel process to the exhumation of the rock massif, the weathering of the rocks is generated by pressure relief, where igneous and metamorphic rocks present a decompression by lithostatic load, which produces a fracturing by expansion or dilation of the rocks along the surfaces approximately parallel to the ground [71].
With the approximation performed on the basis of the exhumation rate of the rocky massif, an estimate of the lifting generated in the last 25 Ma is calculated.As a result, it is determined that these surfaces with a constant rate and dynamics for the current period would present a relief of approximately 7500 m.These data contrasted with the current topographic information of the massif for the study area, presenting an altitudinal difference of 3594 m (52.08%).This difference is associated with the degradation of the massif by the action of weathering and erosion processes in the hydrographic basins that are mostly in a juvenile state (hypsometric curve) (Figure 9a).With this result, a statistical modeling (approximation) of the rate of degradation of the massif per million years is carried out, based on the percentage of altitudinal difference and the current volume of the hydrographic basins (Figure 9b).These combined effects create scenarios with zones of fragile behavior that develop multiple fracturing planes (joints) resulting in a degradation of the rock mass by geomechanical aspects of the rocks.Generally, this behavior is associated with lithological units of igneous and metamorphic type, which display mechanical degradation.This is a factor considered in this study for the weathering and degradation of the rock massif.These processes generate a high production of fragmented stone material that is eroded by hydrogravitational, pluvial and fluvial agents.As a consequence of these processes, planes of instability are developed that give rise to mass movements mainly of the rock flow and Based on the idea that the production of sediments is not proportional to the exhumation rate of the rock massif (non-linear relationship), the behavior of the graph must be expressed exponentially.For a time T 0 (25 Ma ago) production is low, and for a time T 1 (present) it is greater with an increase in degradation per a million years.This is due to the dynamic behavior of the basins and their relationship to the stabilization of the channels due to the lifting.This behavior of degradation of the massif can be explained in a better way based on the equation of the exponential geometric function of population growth proposed by Robert Malthus and adapted for the present research.Since this function shows the conditions of a population that grows annually, it can be extrapolated as a function of degradation of the massif for every million years.The equation is D( t ) = D 0 × (1 + A) t , where D( t ) is the degradation as a function of time, D 0 is the value of the degradation at a time T 0 , A this is the percentage of altitudinal difference and t is the time in millions of years.The output of this function is displayed in Figure 9.
With the relationship generated in Figure 9 between the lifting of the massif and the degradation curve, a trend of degraded volume of each basin can be obtained.The equation in exponential function described as y = 3281.2e−0.419x estimates the volume of degraded material at any point in a timeline, with potential to estimate the volume in the future if the conditions of the site are not altered; in this case, if the exhumation rate is not altered by terrestrial dynamics.
These combined effects create scenarios with zones of fragile behavior that develop multiple fracturing planes (joints) resulting in a degradation of the rock mass by geomechanical aspects of the rocks.Generally, this behavior is associated with lithological units of igneous and metamorphic type, which display mechanical degradation.This is a factor considered in this study for the weathering and degradation of the rock massif.These processes generate a high production of fragmented stone material that is eroded by hydrogravitational, pluvial and fluvial agents.As a consequence of these processes, planes of instability are developed that give rise to mass movements mainly of the rock flow and fall types of some sectors; these movements give rise to geoforms such as talus cones, dejection cones, channels and alluvial-torrential fans (Figure 10).These combined effects create scenarios with zones of fragile behavior that develop multiple fracturing planes (joints) resulting in a degradation of the rock mass by geomechanical aspects of the rocks.Generally, this behavior is associated with lithological units of igneous and metamorphic type, which display mechanical degradation.This is a factor considered in this study for the weathering and degradation of the rock massif.These processes generate a high production of fragmented stone material that is eroded by hydrogravitational, pluvial and fluvial agents.As a consequence of these processes, planes of instability are developed that give rise to mass movements mainly of the rock flow and fall types of some sectors; these movements give rise to geoforms such as talus cones, dejection cones, channels and alluvial-torrential fans (Figure 10).On the other hand, the climatic conditions present in the study region also provide factors that control the rate of degradation of the rocky massif that is consistent with the evolution of the basins and their modeling.For the region, it is described that processes On the other hand, the climatic conditions present in the study region also provide factors that control the rate of degradation of the rocky massif that is consistent with the evolution of the basins and their modeling.For the region, it is described that processes such as insolation and alternation between wetting and desiccation create constant and progressive states of weathering and erosion.The climate of the Chicamocha Canyon is very exclusive, especially for this region where annual rainfall of 721 mm is recorded with an annual temperature of 25.4 • C. Considering the classification by Thornthwaite [72], this sector of the Chicamocha Canyon basin is described as having a semi-arid climate, without any water gain, and megathermal (DdA'a').There is no storage of useful water during eleven months of the year (November to September) and precipitation is equal to actual evapotranspiration, except in the month of October [73].This condition is accentuated by the surrounding topographic control that influences the behavior of the rains in quantity and distribution due to the dynamics of the wind that is affected by the relief.These orographic rains are subject to the mountainous reliefs that serve as a shield for the dry air currents that tend to rise, generating dry and warm climates in the lower parts, while in the upper parts they are humid and rainy [74].
With the climatic conditions described above, weathering processes of the accentuated insolation type are created due to a systematic repetition of the phenomenon of heating of the surfaces, causing a stress that leads to the rupture of the rocks (Figure 11).As a result of this process, a cataphyllar exfoliation or concentric desquamation is present together with a series of fissures that, when converged, create planes where the rocks are eroded (rock fragments).The different surface temperatures affect the minerals present in the rocks magnifying the thermal expansion, which creates a condition of repetitive stresses that could cause rupture by fatigue failure on adjacent minerals, especially those with a granitic composition (igneous rocks) [75].Additionally, the repetitive alternation of wet and dry states in rocks with high clay contents creates disintegration processes due to increased tension stress.This behavior should be expected for the basins located in the eastern region of the study area (basins of Guaca, Río Negro rivers and the Chorrerana stream) since they present lithologies with more clayey compositions.There is no storage of useful water during eleven months of the year (November to September) and precipitation is equal to actual evapotranspiration, except in the month of October [73].This condition is accentuated by the surrounding topographic control that influences the behavior of the rains in quantity and distribution due to the dynamics of the wind that is affected by the relief.These orographic rains are subject to the mountainous reliefs that serve as a shield for the dry air currents that tend to rise, generating dry and warm climates in the lower parts, while in the upper parts they are humid and rainy [74].With the climatic conditions described above, weathering processes of the accentuated insolation type are created due to a systematic repetition of the phenomenon of heating of the surfaces, causing a stress that leads to the rupture of the rocks (Figure 11).As a result of this process, a cataphyllar exfoliation or concentric desquamation is present together with a series of fissures that, when converged, create planes where the rocks are eroded (rock fragments).The different surface temperatures affect the minerals present in the rocks magnifying the thermal expansion, which creates a condition of repetitive stresses that could cause rupture by fatigue failure on adjacent minerals, especially those with a granitic composition (igneous rocks) [75].Additionally, the repetitive alternation of wet and dry states in rocks with high clay contents creates disintegration processes due to increased tension stress.This behavior should be expected for the basins located in the eastern region of the study area (basins of Guaca, Río Negro rivers and the Chorrerana stream) since they present lithologies with more clayey compositions.Of the morphometric parameters calculated, the Form factor (Rf), Drainage density (Dd), Drainage texture (T), Compactness coefficient (Kc), Ruggedness number (Rn), Melton Ruggedness number (MRn), Torrential coefficient (Ct) and Basin concentration time (Tc) are closely related to climatic factors and basin geometry.Two conditions can be extracted from this characterization.The first is associated with the tributaries of the basins present in the areas proximal to the bed of the Chicamocha River, where channeled mass movements of dry flow are evident.These movements are generated by gravitational processes and high-angle slopes, where the eroded material is accumulated at the tipping points of the slopes that subsequently, by sedimentary load, is set in motion, thus creating channeled dry flows of fragmented material.The second is associated with the upper parts of the basins, where precipitation levels are higher, leading to the generation of critical flows that, when detected with the stone material arranged in the lower parts of the basins, create torrential flows.Depending on the geometry of the basin, the speed, the concentration of the flow and the material transported, these basins present moderate to high conditions in the generation of torrential flows that bring with them an affectation on anthropic elements (houses, roads, bridges, and agricultural areas).
Finally, a factor that would be affecting the control of weathering and erosion of the rocky massif to a lesser extent is seismic activity.The study area is part of one of the regions with the highest seismic activity in the world.According to Richter [76], the sectors with the highest number of seismic events worldwide are defined as "seismic nests".In turn, this high seismic activity must have persistence over time [77].The seismic nest of Bucaramanga, as this sector is called, presents a high productivity of seismic events.Based on the catalog of the National Seismological Network of Colombia, between the years 1993 and 2021, which has a total of 153,469 events and historical data collected between 1826 and 1992 with a total of 11 events, an estimate of the seismic productivity is made by applying the equation of Gutenberg and Richter [78].An important fact to consider is that no data processing was carried out in the catalog, such as parametric and non-parametric controls, data standardization, elimination of anomalous points, selection by seismic source and selection of a single magnitude class.This is because it is sought to determine the seismic productivity present in the region as well as the maximum magnitude it reaches and its probability of occurrence (Figure 12).channeled dry flows of fragmented material.The second is associated with the upper parts of the basins, where precipitation levels are higher, leading to the generation of critical flows that, when detected with the stone material arranged in the lower parts of the basins, create torrential flows.Depending on the geometry of the basin, the speed, the concentration of the flow and the material transported, these basins present moderate to high conditions in the generation of torrential flows that bring with them an affectation on anthropic elements (houses, roads, bridges, and agricultural areas).Finally, a factor that would be affecting the control of weathering and erosion of the rocky massif to a lesser extent is seismic activity.The study area is part of one of the regions with the highest seismic activity in the world.According to Richter [76], the sectors with the highest number of seismic events worldwide are defined as "seismic nests".In turn, this high seismic activity must have persistence over time [77].The seismic nest of Bucaramanga, as this sector is called, presents a high productivity of seismic events.Based on the catalog of the National Seismological Network of Colombia, between the years 1993 and 2021, which has a total of 153,469 events and historical data collected between 1826 and 1992 with a total of 11 events, an estimate of the seismic productivity is made by applying the equation of Gutenberg and Richter [78].An important fact to consider is that no data processing was carried out in the catalog, such as parametric and non-parametric controls, data standardization, elimination of anomalous points, selection by seismic source and selection of a single magnitude class.This is because it is sought to determine the seismic productivity present in the region as well as the maximum magnitude it reaches and its probability of occurrence (Figure 12).From the estimation made of seismic productivity, it can be inferred that eighteen events of magnitude 2 or one event of magnitude 3 can be generated in one day.With From the estimation made of seismic productivity, it can be inferred that eighteen events of magnitude 2 or one event of magnitude 3 can be generated in one day.With these data, a surface scenario can be extrapolated where the rock mass would be in constant excitement due to the vibrations generated in the ground.These vibrations contribute to the processes of mechanical weathering of the rocks creating microcracks that expand by the repeated action of telluric movements.The formation of these microcracks is due to a mechanical induction of the rock massif, which generates a tensional stress that produces cracks caused by local stress that exceeds local resistance [79].This would be a non-inclusive factor within the dynamic activity of the watersheds, but to some extent would activate the states of weathering and erosion of the rocks due to the constant process of vibration of the rock massif.

Conclusions
The present study identified three factors that influence basin modelling, sediment generation and ground instability that give rise to mass movements.The first factor is the most influential and obeys a structural control exercised by the stroke and kinematics of the Bucaramanga fault, where it modifies the geometry, arrangement, and dynamics of the basins, creating planes of weakness that give rise to processes of weathering and accentuated erosion, as well as the generation of slope instability events (rock falls, rock flows and rockslides and/or detritus).The second factor is related to climate control, where the region presents special conditions of temperature, precipitation, and evapotranspiration, giving rise to a dynamic control of the basins, as well as the degradation of the rock massif because of insolation, wetting and drying of the rocks.The combination of the rainy states in the upper parts of the basins, together with the high productivity of sediments in the lower parts, create favorable scenarios for the generation of torrential flows that greatly affect elements of anthropogenic nature.The third factor is focused on the recurring seismic activity in the region.This activity generates rock dislocation processes due to tensional stress, increasing the productivity of rock fragments or contributing to the other two factors in accentuating and magnifying the weathering and mechanical erosion processes present in the region.The combination of the three factors establishes the line of evolution of the landscape present in the region, as well as controls the modeling and the rates of weathering and erosion of the basins, which, in turn, develop the dynamics of the channels that are prone to generating torrential flows.
Within the unified morphometric evaluation of the hydrographic basins, it was observed that the Manco River basin and the Talarcuta, Las Pavas, Felisco and El Abra streams present the geometries, surfaces, and dynamics with the greatest possibilities in the generation of torrential flows.Proof of this is the torrential flow generated by heavy rains over the Manco River basin and its tributaries on 26 February 2020, causing multiple losses of human lives and considerable damage to infrastructure [80].On the other hand, the basins that present less development to this type of movement are the Umpalá, Guaca, and Río Negro River basins and the Cantabara and La Cureña streams; however, these basins in their lower parts present a high degree of weathering of the rocky massif, giving rise to the formation of hydrogravitational deposits that are arranged on the margins of the channels and provide material for the evolution of torrential flows.
According to the evaluation of the morphometric parameters carried out for all the basins, it is estimated that when heavy or constant rains in the upper parts of the basins are combined with the material arranged by the processes of weathering and erosion that generate deposits in the lower areas of the basins, movements of torrential flow type occur.To some extent, within the evaluation, it is estimated that the development of these movements would be considerably controlled by the amount of rainfall that may occur in the upper parts of the basins.This material is suggested as a reference to consider within the analysis of territorial planning and the use of natural resources.

Figure 1 .
Figure 1.Location of the study area with delimitation of the hydrographic basins.

Figure 1 .
Figure 1.Location of the study area with delimitation of the hydrographic basins.
c is the concentration time (hours), S 0 is the average slope of the main channel (m/m) and A is the area of the basin (km 2 )[62].(h)Ventura-Heron method: T c = 0.3
the northern region (right bank, downstream of Chicamocha River) mostly have Condition 2, except Río Negro and Guaca rivers (Condition 1) and Quebrada Seca (Condition 3).The basins present in the southern region (left bank, downstream of Chicamocha River) such as Felisco, Basteña and Calichana streams have Condition 1, Cantabara and Lauchas streams present a transition phase between Conditions 1 and 2, and Cureña stream is in Condition 2. This shows that a large part of the basins is in the process of deepening to reach a channel balance.Sustainability 2023, 15, x FOR PEER REVIEW 12 of 27 in the northern region (right bank, downstream of Chicamocha River) mostly have Condition 2, except Río Negro and Guaca rivers (Condition 1) and Quebrada Seca (Condition 3).The basins present in the southern region (left bank, downstream of Chicamocha River) such as Felisco, Basteña and Calichana streams have Condition 1, Cantabara and Lauchas streams present a transition phase between Conditions 1 and 2, and Cureña stream is in Condition 2. This shows that a large part of the basins is in the process of deepening to reach a channel balance.

Figure 6 .
Figure 6.(a) Hypsometric curves for the northern region.(b) Hypsometric curves for the southern region.The red line indicates (mature state) the basis for the analysis of hypsometric behavior.

Figure 6 .
Figure 6.(a) Hypsometric curves for the northern region.(b) Hypsometric curves for the southern region.The red line indicates (mature state) the basis for the analysis of hypsometric behavior.

Figure 7 .
Figure 7. Index of the gradient of the longitudinal profile of the riverbed (SL) for each of the hydrographic basins.Figure 7. Index of the gradient of the longitudinal profile of the riverbed (SL) for each of the hydrographic basins.

Figure 7 .
Figure 7. Index of the gradient of the longitudinal profile of the riverbed (SL) for each of the hydrographic basins.Figure 7. Index of the gradient of the longitudinal profile of the riverbed (SL) for each of the hydrographic basins.

Figure 8 .
Figure 8.(a) Statistical estimation of Basin concentration time (Tc) for each of the equations described.Time in minutes.(b) Correlation between concentration times and torrential coefficient.

Figure 8 .
Figure 8.(a) Statistical estimation of Basin concentration time (Tc) for each of the equations described.Time in minutes.(b) Correlation between concentration times and torrential coefficient.

Figure 9 .
Figure 9. Correlation of exhumation and degradation values of the rock massif.Contrast between the estimate of the solid and the current expression with the degradation curve and current volume of degraded basin area.

Figure 9 .
Figure 9. Correlation of exhumation and degradation values of the rock massif.Contrast between the estimate of the solid and the current expression with the degradation curve and current volume of degraded basin area.

Figure 9 .
Figure 9. Correlation of exhumation and degradation values of the rock massif.Contrast between the estimate of the solid and the current expression with the degradation curve and current volume of degraded basin area.

Figure 10 .
Figure 10.Panoramic image showing the mechanical degradation on the surfaces of several hydrographic basins.Slope movements associated with rock flows can be seen.(a) Manco River basin (lower part).(b) Umpalá River basin (lower part).(c) El Abra stream basin.(d) Seca stream basin (lower part).(e) Talarcuta stream basin (lower part).(f) Chicamocha River heading northwest.

Figure 10 .
Figure 10.Panoramic image showing the mechanical degradation on the surfaces of several hydrographic basins.Slope movements associated with rock flows can be seen.(a) Manco River basin (lower part).(b) Umpalá River basin (lower part).(c) El Abra stream basin.(d) Seca stream basin (lower part).(e) Talarcuta stream basin (lower part).(f) Chicamocha River heading northwest.

Sustainability 2023 ,
15,  x FOR PEER REVIEW 22 of 27 such as insolation and alternation between wetting and desiccation create constant and progressive states of weathering and erosion.The climate of the Chicamocha Canyon is very exclusive, especially for this region where annual rainfall of 721 mm is recorded with an annual temperature of 25.4 °C.Considering the classification byThornthwaite [72], this sector of the Chicamocha Canyon basin is described as having a semi-arid climate, without any water gain, and megathermal (DdA'a').

Figure 11 .
Figure 11.Weathering by insolation (thermoclasty).(a) Concentric scaly.(b) Rupture due to fatigue failure.Of the morphometric parameters calculated, the Form factor (Rf), Drainage density (Dd), Drainage texture (T), Compactness coefficient (Kc), Ruggedness number (Rn), Melton Ruggedness number (MRn), Torrential coefficient (Ct) and Basin concentration time (Tc) are closely related to climatic factors and basin geometry.Two conditions can be extracted from this characterization.The first is associated with the tributaries of the basins present in the areas proximal to the bed of the Chicamocha River, where channeled mass movements of dry flow are evident.These movements are generated by gravitational processes and high-angle slopes, where the eroded material is accumulated at the tipping points of the slopes that subsequently, by sedimentary load, is set in motion, thus creating

Figure 12 .
Figure 12.Productivity of the seismic nest of Bucaramanga (evaluation with data from the 1826-1993 "historical catalogue" and the 1993-2021 "new registrations").(a) Linear fit of the Gutenberg-Richter equation by least squares.(b) Number of events expected from the solution of the Gutenberg-Richter equation.(c) Maximum expected magnitude.(d) Poisson distribution that estimates the probability of an event of a certain magnitude occurring.

Figure 12 .
Figure 12.Productivity of the seismic nest of Bucaramanga (evaluation with data from the 1826-1993 "historical catalogue" and the 1993-2021 "new registrations").(a) Linear fit of the Gutenberg-Richter equation by least squares.(b) Number of events expected from the solution of the Gutenberg-Richter equation.(c) Maximum expected magnitude.(d) Poisson distribution that estimates the probability of an event of a certain magnitude occurring.

Author
Contributions: Conceptualization, data collection, data analysis, writing draft and final manuscript, J.A.V.O.; concept, writing draft manuscript, A.M.M.-G.; writing of final manuscript, commentary, and revision, J.A.V.O. and A.M.M.-G.All authors have read and agreed to the published version of the manuscript.

Table 1 .
Methodology adopted for the calculation of morphometric parameters.

Table 2 .
Morphometric parameters calculated for all basins.Drainage network.
The green box indicates the quantitative valuation in a low range for this parameter.The red box indicates the quantitative valuation in a high range for this parameter.

Table 3 .
Morphometric parameters calculated for all basins.Relief characterization.
The green box indicates the quantitative valuation in a low range for this parameter.The red box indicates the quantitative valuation in a high range for this parameter.

Table 4 .
Morphometric parameters calculated for all basins.Drainage texture and flow mobility analysis.
The green box indicates the quantitative valuation in a low range for this parameter.The red box indicates the quantitative valuation in a high range for this parameter.