Next Article in Journal
Ocean Recreation and the Economic Contributions of Visitation in and around Gray’s Reef National Marine Sanctuary
Previous Article in Journal
Effects of Temperature and Water Types on the Decay of Coronavirus: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping of Areas Vulnerable to Flash Floods by Means of Morphometric Analysis with Weighting Criteria Applied

1
Doctoral Program in Water Resources and Energy for Agriculture, Universidad de Concepcion, Av. Vicente Mendez 595, Chillan 3812120, Chile
2
CRHIAM Water Research Center, Universidad de Concepción, Victoria 1295, Concepción 4070386, Chile
3
Water Resources Department, College of Agricultural Engineering, Universidad Nacional Agraria La Molina, Av. La Molina s/n, Lima 15024, Peru
4
Water Resources Department, College of Agriculture Engineering, Universidad de Concepción, Av. Vicente Mendez 595, Chillan 3812120, Chile
5
Facultad de Ingeniería, Departamentol Ingeniería Civil, Universidad de Concepción, Concepción 4070386, Chile
*
Author to whom correspondence should be addressed.
Water 2023, 15(6), 1053; https://doi.org/10.3390/w15061053
Submission received: 27 January 2023 / Revised: 7 March 2023 / Accepted: 7 March 2023 / Published: 10 March 2023
(This article belongs to the Section Hydrology)

Abstract

:
Flash floods, produced by heavy seasonal rainfall and characterized by high speeds and destructive power, are among the most devastating natural phenomena and are capable of causing great destruction in very little time. In the absence of hydrological data, morphometric characterization can provide important information on preventive measures against flash floods. A priority categorization of hydrographic units in the Cañete River basin was carried out using morphometric analysis together with a weighted sum analysis (WSA) based on a statistical correlation matrix. The delineation of the drainage network was performed based on Digital Elevation Model (DEM) data from the Shuttle Radar Topography Mission (SRTM). The Cañete River basin was subdivided into 11 sub-basins, and 15 morphometric parameters were selected. The priority category (very high, high, and moderate) of each sub-basin was assigned according to the value of the composite factor obtained through WSA. The results of this analysis showed that 26.08% of the total area is under a very high flash flood risk (sub-basins 3, 9, and 11), 38.46% is under a high flash flood risk (sub-basins 5, 7, 8, and 10), and 35.45% is under a moderate flash flood risk. This study concludes that flash floods predominate in sub-basin 3 and that downstream areas present characteristics of river flooding (sub-basins 9 and 11).

1. Introduction

Flash floods are among the most devastating natural phenomena and are capable of causing great destruction in very little time [1]. These natural disasters are most often caused by extreme rainfall events (climate change), which occur during a short time in a relatively small area, resulting in excessive discharge, landslides, and mudflows [2]. Climate change is one of the greatest threats to the entire world, affecting Earth’s natural balance and ecosystems and disrupting all life [3]. This phenomenon is often measured by changes in primary variables such as temperature and precipitation [4]. The prediction and control of flash floods have been very difficult due to the highly dynamic nature of the climate and their sudden appearance; under these circumstances, a greater emphasis on prevention and protection studies in areas vulnerable to natural disasters is necessary [2].
In the absence of hydrological data, morphometric analysis can provide important information on the hydrological characteristics of a basin [5,6]. Quantitative morphometric assessment of a basin is a very effective method for interpreting various aspects of its drainage network and evaluating its hydrological behavior [7,8,9]. Studies on basin sub-unit prioritization that have been conducted based on morphometric analysis have aided planners and decision-makers in the development of management plans [10,11,12,13,14,15]. Thus, to achieve good drainage basin management, it is necessary to study their morphometry [5,7].
Various studies have evaluated susceptibility to flash floods in drainage basins based on an analysis of morphometric parameters, and quantitative assessment of the morphometric characteristics of a basin has allowed its hydrological response behavior to be defined [1,9,16,17,18,19,20]. Remote sensing and Geographic Information Systems (GIS) are ideal tools for morphometric analysis based on the treatment and quantification of topographic data [7,9,10,21,22,23,24,25]. Recent publications on basin prioritization have used various methodologies such as principal component analysis (PCA) [26,27,28,29], weighted sum analysis (WSA), also known as multicriteria analysis [9,16,18,22,23,24,28], simple additive weighting (SAW), and the technique for order of performance by similarity to ideal solution (TOPSIS) [30]. All these methodologies, together with statistical techniques, have allowed ranked priority scores to be assigned to each sub-basin according to its relationship with flash flood risk [20]. However, there are differences among the various applied methodologies in the quantity and types of parameters used in their analyses.
This investigation aims to identify zones vulnerable to flash floods based on the characterization of morphometric parameters (linear aspects, area, and relief), use weighted sum analysis (WSA) based on a statistical correlation matrix, and set a priority category for each hydrographic unit. In addition, in the absence of hydrological data, the alternative is to use topographic data as an input to identify flash flood risk zones by analyzing morphometric characteristics. The methodology adopted in this study is to use GIS tools, programming, and statistical analysis in mapping flash flood risk zones. The results show morphometric characteristics of 11 fourth-order sub-basins with very high, high, and moderate susceptibility to flash floods. The article is organized as follows: Section 2 presents information on the location of the study area, the morphometric parameters considered, the dataset, and the methodology. Section 3 presents and discusses the results. Finally, the conclusions of this study are described in Section 4.

2. Materials and Methods

2.1. Study Area

The Cañete River basin is part of the Peruvian Pacific slope, located between geographic coordinates 11°58′ and 13°10′ south and 76°25′ and 75°30′ west, in the Cañete and Yauyos provinces, in the Department of Lima, Perú. The river basin covers an area of 6192 km2. It has a maximum altitude of 5800 masl in the central Andes, and the river flows into the Pacific Ocean at 0 masl (Figure 1). The 79.5% of the basin located above 2500 masl is classified as wet. The Cañete River originates at Ticllacocha Lake and is located in the foothills of the Ticlla and Pichahuarco ranges at an altitude of 4429 masl. After flowing 235.67 km, with an average slope of 1.85%, it empties into the Pacific Ocean. Along its path, the Cañete River receives the contributions of numerous tributaries, including, on its right bank, the Miraflores, Yauyos, Huantuya (Carania), and Aucampi rivers and, on its left bank, the Alis, Laraos, Huantán, Tupe, Cacra, and Huangascar rivers, primarily [31,32].

2.2. Morphometric Parameters

These indicators are necessary to understand and estimate the hydrological and morphological characteristics of a basin; they also allow an interpretation of linear aspects, area, and relief of a drainage basin [1,14,15,20,33,34,35]. Remote sensing and Geographic Information Systems (GIS) are effective tools for delineating and understanding the morphometry of any drainage basin [7,9,35,36].
The morphometric parameters used in the study are presented in Table 1, Table 2 and Table 3.

2.3. Data Used

In this study, Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) data were used to delineate the basin, extract the drainage/stream network, and subdivide the catchment area into sub-basins. The SRTM V3 (SRTM Plus) product is provided by the National Aeronautics and Space Administration (NASA) and has a resolution of 1 arc-second (approximately 30 m) [44].
The DEMs were downloaded from Google Earth Engine (GEE), a cloud-based platform that facilitates access to a catalog of petabytes of publicly available geospatial data, including satellite and aerial images, environmental, meteorological, climate, and topographic variables, land cover, etc. [45,46,47]. GEE is accessed using an application programming interface (API) accessible by the internet and an associated web-based integrated development environment (IDE) that allows the quick creation of prototypes and visualization of results [45,47,48,49]. In this study GEE was accessed using the Kaggle web platform, which offers a customizable Jupyter Notebooks interface, using Python programming language [50]. A rectangular region that covers the entire Cañete River basin was downloaded. Data Availability shows the codes that were used.

2.4. Methodology

The methodology adopted in this study is described by the following steps: (1) Extraction of drainage network and sub-basin delineation. (2) Morphometric analysis using the QGIS tool. (3) Assignment of preliminary ranking of sub-basin priority. (4) Weighted sum analysis and final ranking, y. (5) Mapping of areas vulnerable to flash floods. Figure 2 presents the methodological diagram of this research.

2.4.1. Extraction of Drainage Network and Sub-Basin Delineation

The drainage network extraction and basin delineation were carried out following DEM processing methods (cell filling, flow direction, flow accumulation, stream order definition, stream segmentation, and basin delineation) using the QSWAT complement in QGIS 3.16 [51,52]. To define the streams in the study area, a threshold value corresponding to 1% (192,777 cells) of the DEM cell count was used. The Cañete River basin was subdivided into 11 sub-basins, designated as SC-1 to SC-11, and the entire drainage network was classified as a fourth-order catchment basin using the Strahler method [37].

2.4.2. Morphometric Analysis

A morphometric analysis is necessary to ascertain the sub-basin prioritization. This process is advantageous as the derived basin variables are in the form of ratios or dimensionless numbers, providing an effective comparison independent of scale [1,14,14,20,24,34,35]. Morphometric analysis offers a complete representation of the drainage network, geometry, and topography of the basin, allowing an interpretation of linear aspects, area, and relief of the basin, respectively [9].
Each of these parameters can directly or indirectly influence the surface runoff of the drainage basin [9], facilitating the understanding of how they influence flash flood risk. Table 1, Table 2 and Table 3 show a list of morphometric parameters that describe the dimensionless and multidimensional characteristics of the basin, which are grouped into 3 aspects. First, linear morphometric characterization: stream order (u), stream length (Lᵤ), and basin length (Lb), shown in Table 1 [37,38,39]. Second, morphometric characterization of area: basin area (A), basin perimeter (P), stream frequency, also known as flow frequency (Fs), drainage density (Dd), form factor (Ff), circularity ratio (Cr), drainage texture (Tr), elongation ratio (Er), and shape factor (Sf), shown in Table 2 [38,41,42,43]. Finally, basin relief (R), relief ratio ( R r ), and average slope (As) are morphometric relief parameters, shown in Table 3 [41].

2.4.3. Preliminary Ranking of Sub-Basin Priority

The preliminary ranking (PR) approach was used, in which morphometric parameters were divided into two groups according to their ability to (directly or indirectly) influence the adversity of the morphometric parameter conditions or degree of susceptibility to flash floods [16,18,19,22,24]. The first group of parameters (capable of direct influence) consisted of: basin relief (R), relief ratio ( R r ), drainage density (Dd), stream frequency (Fs), circularity ratio (Cr), drainage texture (Tr), and average slope (As). The second group of parameters (capable of indirect influence) consisted of: elongation ratio (Er), form factor (Ff), and shape factor (Sf).
The preliminary rankings of the parameters in the first group were assigned in such a way that for each parameter the sub-basin with the highest value was classified as 1, that with the next highest value was classified as 2, and so on for the remaining sub-basins. The opposite was done for the parameters in the second group, assigning rankings in such a way that for each parameter the sub-basin with the lowest value was placed in position 1, that with the next lowest value was placed in position 2, and so on for the remaining sub-basins [16,18,22].

2.4.4. Weighted Sum Analysis and Final Ranking

Weighted sum analysis (WSA) is a well-known method that provides consistency in addressing complicated problems to compare land surface processes in related entities such as drainage basins [26]. WAS, also known as multi-criteria decision making, is widely used to select the best alternatives among multiple options. Based on preliminary rankings, the correlation matrix and correlation coefficients were calculated [16,18]. The composite parameters ( W S A c p ) were calculated using the following equation (Equation (1)) [9,16,22,26]:
W S A c p = P R p 1 × W p 1 + P R p 2 × W p 2 + + P R p n × W p n
where W S A c p = composite parameter used for weighted sum analysis; P R = preliminary priority ranking of each morphometric parameter (p1, p2, … pn); and W = weights of the morphometric parameters obtained by means of the correlation matrix, which was calculated using the following equation (Equation (2)):
P a r a m e t e r   w e i g h t s   W = C o r r e l a t i o n   c o e f f i c i e n t   s u m C o r r e l a t i o n   t o t a l
Therefore, sub-basin priority was assigned by taking the W S A c p + values (capable of direct influence) corresponding to the first group and subtracting the W S A c p   values (capable of direct influence) corresponding to the second group, following the equation ( Priority   = W S A c p + W S A c p ). Subsequently, final ranking values were assigned, one for the lowest priority value, two for the following value, and so on. Finally, to categorize the priority type, the maximum category (very high) corresponds to the lowest priority levels [9,18,22]. Finally, a model was formulated to assess the final priority using the W S A c p value of each parameter, as shown below:
P r i o r i z a c i o n = R + R r + D d + F s + C r + T r + A s E r + F f + S f

3. Results and Discussion

The identification of areas in drainage basins vulnerable to flash floods using morphometric parameters is considered one of the most effective methods to characterize various geohydrological properties of a basin [14,24]. Therefore, this study has used different morphometric parameters that govern the hydrological response of a basin to prioritize the sub-basins of the Cañete River in terms of their vulnerability to flash floods.

3.1. Morphometric Analysis of the Basin

The morphometric analysis was carried out for 11 sub-basins, and the drainage network of the entire study area was considered a fourth-order basin (Figure 3). The morphometric analysis covered 15 parameters for each sub-basin, which was necessary to determine its dimensions, shape and area, and the characteristics of the drainage network. The results reveal that SB-9 is the smallest, with an area of 31.49 km2, while SB-11 is the largest, with an area of 1175.47 km2 (Table 4). More than half of the sub-basins have areas greater than 400 km2, with only SB-9 having an area smaller than 50 km2. Most of the sub-basins are large, with areas greater than 300 km2. The area of a sub-basin directly affects its susceptibility to flash floods [16,19].
Basin perimeter (P) values for the study area vary from 30.01 km (SB-9) to 281.88 km (SB-11). Basin length (Lb) is understood as the maximum length and is calculated from the farthest point on the basin border to the confluence point; the values vary from 9.31 km (SB-9) to 72.75 km (SB-11). The values are presented in Table 4.
Figure 3 and Table 5 show that there are 97 first-order streams, 41 second-order streams, 18 third-order streams, and 37 fourth-order streams. Table 5 also shows that the greatest number of streams was found in SB-11 (49), and the smallest number of streams was found in SB-9 (1).
Table 6 shows the lengths of the first- to fourth-order streams of the 11 sub-basins, determined using the GIS technique [51], which indicates that the total length of all streams varies from a minimum of 7.68 km (SB-9) to a maximum of 247.19 km (SB-11), and that the total length of all streams in the basin is 1113.75 km. The drainage length values show a significant interrelationship with drainage flow discharge and the erosion phase of the river basin [18]. In addition, it is seen that the five sub-basins with fourth-order streams (SB-3, SB-5, SB-7, SB-9, and SB-11) are possibly susceptible to flash floods [16,19]. Drainage length is a parameter that is directly related to the hydrological reaction of a watershed and its importance is approximately equal to half the reciprocal of the drainage density [19].
Drainage density (Dd) and stream frequency (Fs) were calculated for all sub-basins (Table 7). In the analysis of the sub-basin and drainage network map (Figure 3), the Cañete River basin was found to be a fourth-order basin, with a dendritic drainage pattern. Dd and Fs are significant parameters that contribute to the hydrological responses in each sub-basin [19]. The Fs and Dd values in each sub-basin are directly proportional; therefore, a greater stream number corresponds to a greater stream length. The high value of Fs in sub-basins 3, 7, and 11 indicates that they produce more runoff in comparison to the other sub-basins. Meanwhile, high Dd values are observed in sub-basins 9 and 11, indicating a well-developed network, which is conducive to high runoff concentrations that give rise to flash floods. The low value of Dd in sub-basins 4 and 6 indicates that they have highly permeable subsoil material, with dense plant cover and a low relief [53]. Therefore, high values of Dd and Fs are likely to increase susceptibility to high surface runoff and flash floods [19].
The circularity ratio (Cr) is influenced by stream length and frequency, geological structures, land use and cover, climate, relief, and basin slope [35]. For a perfectly circular basin, the value of the circularity ratio is 1 [43]. In this study, the Cr values of the sub-basins range from 0.184 to 0.439 (Table 7); all the values are less than 1, which indicates that the area is characterized by an elevated relief and that the drainage system is structurally controlled.
Infiltration capacity is the only important factor that controls the texture ratio (Tr) [9]. The texture ratio in the sub-basins range from 0.00 to 0.085 (Table 7). The highest relief values are attributed to the pronounced slopes of the sub-basins, which makes them more prone to flood risk [16,24]. Meanwhile, the lowest Tr values indicate that the basin is flat, with slope values close to zero degrees [35].
The elongation ratio (Er) is the measure of basin dimensions or basin shape [41]. According to [37], a sub-basin with an Er value above 0.9 is classified as circular, those with Er values between 0.9 and 0.8 are classified as oval, those with Er values between 0.7 and 0.8 are classified as less elongated, those with Er values between 0.5 and 0.7 are classified as elongated, and those with Er values below 0.5 are classified as more elongated. In this study, the Er values vary between 0.532 and 0.680, indicating that the sub-basins are classified as elongated.
The form factor (Ff) is a dimensionless ratio of the area of a basin to the square of its length [42] and can be effectively related to flood occurrence, erosion intensity, and sediment transport capacity in a basin [1]. The lower the form factor value, the longer the basin. Basins with a high form factor have high, short-duration maximum flows, while those with a low form factor have lower, longer-duration maximum flows [16,35]. In this study the sub-basins have Ff values ranging from 0.222 to 0.363 (Table 7), indicating that they have elongated shapes and suggesting lower, longer-duration maximum flows.
The shape factor (Sf) is a dimensionless ratio of the square of the length of a basin to its area [42]. Its values indicate the opposite of those of Ff, with a maximum Sf value corresponding to a minimum Ff value. The Sf values for each sub-basin vary between 2.752 (SB-9) and 4.502 (SB-11) (Table 7).
Basin relief (R) is defined as the change between its highest- and lowest-elevation points [9,18]. Table S1 (Supplementary Materials) details the minimum and maximum height values and other statistics for each sub-basin. R shows the potential energy of a drainage basin, which significantly influences the channel gradient and aspect and landform evolution; therefore, it directly affects surface runoff, flood patterns, and sediment transport [9,18,24]. The highest and lowest R values are found in SB-11 (4420.00 m) and SB-2 (2061.00 m) (Table 8). Maximum R values are indicative of the potential energy of a given sub-basin to move water and sediment along the slope.
The relief ratio (Rr) is estimated as the ratio of basin relief to basin length. According to [41], there is a correlation between the hydrological characteristics of a basin and the relief ratio. Therefore, it is presented as an indicator of the intensity of the erosion process in the basin [9,16]. High R values are characteristic of mountainous regions. The R values for all the sub-basins (Table 8) are between 0.045 (SB-1) and 0.258 (SB-9). Average slope refers to the amount of inclination of the physical feature or the topographic form with respect to the horizontal. Table S2 (Supplementary Materials) details the slope raster statistics in degrees for each sub-basin. Slope analysis is very important in morphometric studies. Slope elements are, in turn, controlled by climate-morphogenic processes in areas with rocks of varying resistance [1,54]. The average slopes of the sub-basins vary from 19.485° (SB-6) to 34.042° (SB-9) (Table 8). The pronounced slopes also favor a faster movement of surface runoff.

3.2. Assignment of Preliminary Sub-Basin Priority Rankings

The seven morphometric parameters (R, Rr, Dd, Fs, Cr, Tr, and As) are directly proportional to soil degradation and water factors. The rankings were assigned from greatest to lowest priority, i.e., rank 1 for the sub-basin with the maximum parameter value and rank 11 for the sub-basin with the minimum parameter value. For example, parameter R (Table 8) with the maximum value of 4420.0 m (SB-11) was assigned the highest priority (rank 1), the next descending value was assigned rank 2, and this went up to the minimum parameter value of 2061.0 m (SB-2), which was assigned the lowest priority (rank 11). The assignment of the rankings for the six parameters Rr, Dd, Fs, Cr, Tr, and As was assigned in a similar way as explained above, the results are shown in Table 9.
The three remaining parameters (Er, Ff, and Sf) have an inverse relationship with soil degradation and water factors. Rankings were assigned from lowest to highest priority, i.e., rank 1 for the sub-basin with the lowest parameter value and rank 11 for the sub-basin with the highest parameter value. For example, parameter Er (Table 7) with a minimum value of 0.532 (SB-11) was assigned the highest priority (rank one), and the next higher value was assigned rank two; the maximum parameter value of 0.680 (SB-9) was assigned the lowest priority (rank 11). The assignment of the rankings for the next two parameters (Ff and Sf) were similarly assigned, the results are shown in Table 9.
A correlation matrix of the 10 morphometric parameters capable of directly or indirectly influencing the vulnerability to flash floods is presented in Table 10. It was estimated based on the preliminary priority rankings (Table 9), and the correlation coefficients were obtained in a process carried out using Python code (Data Availability shows the codes used). The statistical correlation matrix shows that basin relief has a positive correlation with relief ratio, drainage density, stream frequency, texture ratio, average slope, elongation ratio, and form factor and an inverse correlation with circularity ratio and shape factor. The relief ratio has a positive correlation with basin relief, drainage density, circularity ratio, average slope, and shape factor and an inverse correlation with stream frequency, texture ratio, elongation ratio, and form factor. The drainage density, stream frequency, circularity ratio, texture ratio, average slope, elongation ratio, form factor, and shape factor correlations were also calculated and are shown in Table 10.
The final weights for each parameter were calculated by dividing the correlation coefficient sum of each parameter by the overall correlation total (Equation (2)). The final weight of R (weight = 0.171) was obtained though division (1.991/11.636). The final weights of the other morphometric parameters were obtained in a similar manner; the results are shown at the bottom of Table 10. The value 11.636 was obtained by summing the sum row (penultimate row Table 10).

3.3. Final Ranking Using Weighted Sum Analysis

The W S A c p values of the morphometric parameters (Equation (1)) were calculated according to their importance using the weighted sum model. This process used the preliminary priority values (Table 9) and final weights of each morphometric parameter (Table 10). For example, the W S A c p   value of R (SB-1) was obtained by multiplying 8 × 0.171 = 1.37, and similarly for the following W S A c p values of R of the other sub-basins; W S A c p   of Rr was multiplied 11 × 0.009 = 0.09. Therefore, for each parameter there is a corresponding weighting (w). The W S A c p values are shown in Table 11.
Based on the results of the composite weighted sum W S A c p of the different parameters (Table 11), a model was formulated to evaluate the final priority ranking (equation 3). The parameters were categorized into two groups: first group W S A c p   + , parameters that have the ability to directly influence flash floods; second group, W S A c p   , parameters that have the ability to indirectly influence flash floods. The priorities ranking (prioritization) is obtained by subtracting the two groups W S A c p   + W S A c p   ; the sub-basin with the lowest composite value receives the highest priority (one), that with the next lowest composite value receives the second rank, and so on for the classification of each hydrological unit.
The final classification of priorities was carried out in such a way that the lowest composite factor value received priority rank 1, the next lowest value received priority rank 2, and so on for the 11 sub-basins. As observed in Table 12, the highest priority rank (1) was assigned to SB-9, followed by SB-3, SB-11, SB-5, SB-7, SB-10, SB-8, SB-2, SB-6, SB-4, and SB-1. Figure 4 shows the final map of priority classifications of the 11 studied sub-basins.
From the results shown in Table 12, there are composite parameter values W S A c p that have the ability to directly or indirectly influence the degree of susceptibility to flash floods. Based on the composite factor value, the 11 sub-basins of the Cañete River basin were classified into three priority categories: very high, high, and medium [9]. In Table 12 it can be seen that three sub-basins (SB-9, SB-3, and SB-11) are in the “very high” category, four sub-basins (SB-5, SB-7, SB-10, and SB-8) are in the “high” category, and four sub-basins (SB-2, SB-6, SB-4, and SB-1) are in the “medium” category. The final priority category map of the 11 sub-basins (Figure 4) shows that the “very high” category accounts for 26.08% of the total area, the “high” category for 38.46%, and the “medium” category for 35.45%.
These results show that 64.55% of the sub-basins (seven sub-basins) are in zones with a very high or high propensity to soil erosion, indicating that there are potential areas to carry out soil protection measures for efficient basin management and development. In addition, with the final priority classification of the 11 sub-basins (Figure 4), and considering the continuity of the flow wave or flood that starts in SB-3 and moves downstream and is produced by a maximum rainfall event over the central region of the sub-basin, SB-5 and SB-7 could also be considered very high categories, as they follow the continuity of SB-3. Therefore, these results will be of great help to policy-makers, planners, and managers to address vulnerable areas through specific action plans for flood risk reduction.

4. Conclusions

Through the morphometric characterization carried out by interpreting the linear, area, and relief aspects of a hydrological basin, it was possible to identify areas vulnerable to flash floods in the occurrence of extreme events, as a methodological technique in basins with an absence of hydrological data and to know the hydrological behavior in a basin. In addition, the important role of GIS tools and statistical approaches in developing research was shown.
The results show that SB-3, SB-9, and SB-11 are susceptible to floods and soil loss. Flash floods predominate in the upstream sub-basin (SB-3), while the downstream sub-basin (SB-9 and SB-11) present characteristics of river floods. Both are destructive in the study area, affecting the human population and their residential units, farmland with standing crops, and infrastructure (highways, bridges, sewers, and water supply systems).
In addition, the sub-basins located in the middle part of the basin (SB-5, SB-7, SB-8, and SB-10) are categorized as high risk and the remaining sub-basins (SB-1, SB-2, SB-4, and SB-6) as medium risk. Thus, the Cañete River basin has sub-basins that present a very high, high, and moderate susceptibility to flash floods. Finally, it can be affirmed that the hydrological response and specifically the risk of flash floods and extreme events in a river basin depends on its morphometric characteristics. Therefore, this methodology is presented as an alternative to decision-makers in the application of suitable drainage basin management techniques in terms of soil and water conservation measures, allowing them to safeguard the study area and mitigate its degradation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15061053/s1; Table S1: DEM statistics for each sub-basin; Table S2: raster statistics of slope in grade for each sub-basin. Python codes: CorrelationAnalysis and DownloadDEM_CodeKaggle.

Author Contributions

Conceptualization, M.P.-M. and J.L.A.; methodology, M.P.-M.; software, M.P.-M.; validation, M.P.-M. and J.L.A.; formal analysis, M.P.-M. and D.R.-R.; investigation, M.P.-M., J.L.A., A.S., O.L., E.C.-V. and D.R.-R.; writing—original draft preparation, M.P.-M. and D.R.-R.; writing—review and editing, M.P.-M., J.L.A., A.S., O.L., E.C.-V. and D.R.-R.; visualization, M.P.-M., E.C.-V. and D.R.-R.; supervision, J.L.A. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the CRHIAM Water Research Center: ANID/FONDAP/15130015.

Data Availability Statement

The data presented in this study are openly available at Supplementary Materials codes, time series, algorithms, and others.

Acknowledgments

CRHIAM Water Research Center: Project ANID/FONDAP/15130015, Josselyn Portuguez Contreras, and Mariam De la Cruz.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Bisht, S.; Chaudhry, S.; Sharma, S.; Soni, S. Assessment of flash flood vulnerability zonation through Geospatial technique in high altitude Himalayan watershed, Himachal Pradesh India. Remote Sens. Appl. Soc. Environ. 2018, 12, 35–47. [Google Scholar] [CrossRef]
  2. Cahyono, C.; Adidarma, W.K. Influence analysis of peak rate factor in the flood events’ calibration process using HEC-HMS. Model. Earth Syst. Environ. 2019, 5, 1705–1722. [Google Scholar] [CrossRef]
  3. Teng, F.; Huang, W.; Ginis, I. Hydrological modeling of storm runoff and snowmelt in Taunton River Basin by applications of HEC-HMS and PRMS models. Nat. Hazards 2018, 91, 179–199. [Google Scholar] [CrossRef]
  4. Wang, N.; Lombardo, L.; Gariano, S.L.; Cheng, W.; Liu, C.; Xiong, J.; Wang, R. Using satellite rainfall products to assess the triggering conditions for hydro-morphological processes in different geomorphological settings in China. Int. J. Appl. Earth Obs. Geoinf. 2021, 102, 102350. [Google Scholar] [CrossRef]
  5. Prakash, K.; Rawat, D.; Singh, S.; Chaubey, K.; Kanhaiya, S.; Mohanty, T. Morphometric analysis using SRTM and GIS in synergy with depiction: A case study of the Karmanasa River basin, North central India. Appl. Water Sci. 2019, 9, 13. [Google Scholar] [CrossRef] [Green Version]
  6. Perucca, L.P.; Esper Angilieri, Y. Morphometric characterization of del Molle Basin applied to the evaluation of flash floods hazard, Iglesia Department, San Juan, Argentina. Quat. Int. 2011, 233, 81–86. [Google Scholar] [CrossRef]
  7. Rai, P.K.; Singh, P.; Mishra, V.N.; Singh, A.; Sajan, V.; Shahi, A.P. Geospatial approach for quantitative drainage morphometric analysis of Varuna river basin, India. J. Landsc. Ecol. 2019, 12, 1–25. [Google Scholar] [CrossRef] [Green Version]
  8. Rai, P.K.; Chandel, R.S.; Mishra, V.N.; Singh, P. Hydrological inferences through morphometric analysis of lower Kosi river basin of India for water resource management based on remote sensing data. Appl. Water Sci. 2018, 8, 15. [Google Scholar] [CrossRef] [Green Version]
  9. Aher, P.; Adinarayana, J.; Gorantiwar, S. Quantification of morphometric characterization and prioritization for management planning in semi-arid tropics of India: A remote sensing and GIS approach. J. Hydrol. 2014, 511, 850–860. [Google Scholar] [CrossRef]
  10. Shivhare, V.; Gupta, C.; Mallick, J.; Singh, C.K. Geospatial modelling for sub-watershed prioritization in Western Himalayan Basin using morphometric parameters. Nat. Hazards 2022, 110, 545–561. [Google Scholar] [CrossRef]
  11. Abdeta, G.C.; Tesemma, A.B.; Tura, A.L.; Atlabachew, G.H. Morphometric analysis for prioritizing sub-watersheds and management planning and practices in Gidabo Basin, Southern Rift Valley of Ethiopia. Appl. Water Sci. 2020, 10, 158. [Google Scholar] [CrossRef]
  12. Waiyasusri, K.; Chotpantarat, S. Watershed Prioritization of Kaeng Lawa Sub-Watershed, Khon Kaen Province Using the Morphometric and Land-Use Analysis: A Case Study of Heavy Flooding Caused by Tropical Storm Podul. Water 2020, 12, 1570. [Google Scholar] [CrossRef]
  13. Ahirwar, R.; Malik, M.S.; Shukla, J.P. Prioritization of Sub-Watersheds for Soil and Water Conservation in Parts of Narmada River through Morphometric Analysis Using Remote Sensing and GIS. J. Geol. Soc. India 2019, 94, 515–524. [Google Scholar] [CrossRef]
  14. Anees, M.T.; Abdullah, K.; Nawawi, M.; Rahman, N.N.N.A.; Ismail, A.Z.; Syakir, M.; Abdul Kadir, V. Prioritization of Flood Vulnerability Zones Using Remote Sensing and GIS for Hydrological Modelling. Irrig. Drain. 2019, 68, 176–190. [Google Scholar] [CrossRef]
  15. Chauhan, P.; Chauniyal, D.D.; Singh, N.; Tiwari, R.K. Quantitative geo-morphometric and land cover-based micro-watershed prioritization in the Tons river basin of the lesser Himalaya. Environ. Earth Sci. 2016, 75, 498. [Google Scholar] [CrossRef]
  16. Singh, G.; Pandey, A. Morphometric Characterization and Flash Flood Zonation of a Mountainous Catchment Using Weighted Sum Approach. In Geospatial Technologies for Land and Water Resources Management; Pandey, A., Chowdary, V.M., Behera, M.D., Singh, V.P., Eds.; Springer International Publishing: Berlin/Heidelberg, Germany, 2022; pp. 409–428. [Google Scholar] [CrossRef]
  17. El-Fakharany, M.A.; Hegazy, M.N.; Mansour, N.M.; Abdo, A.M. Flash flood hazard assessment and prioritization of sub-watersheds in Heliopolis basin, East Cairo, Egypt. Arab. J. Geosci. 2021, 14, 1693. [Google Scholar] [CrossRef]
  18. Jothimani, M.; Abebe, V.; Dawit, Z. Mapping of soil erosion-prone sub-watersheds through drainage morphometric analysis and weighted sum approach: A case study of the Kulfo River basin, Rift valley, Arba Minch, Southern Ethiopia. Model. Earth Syst. Environ. 2020, 6, 2377–2389. [Google Scholar] [CrossRef]
  19. Mahmood, S.; Rahman, A. Flash flood susceptibility modeling using geo-morphometric and hydrological approaches in Panjkora Basin, Eastern Hindu Kush, Pakistan. Environ. Earth Sci. 2019, 78, 43. [Google Scholar] [CrossRef]
  20. Mahmood, S.; Rahman, A. Flash flood susceptibility modelling using geomorphometric approach in the Ushairy Basin, eastern Hindu Kush. J. Earth Syst. Sci. 2019, 128, 97. [Google Scholar] [CrossRef] [Green Version]
  21. Khan, I.; Bali, R.; Agarwal, K.K.; Kumar, D.; Singh, S.K. Morphometric Analysis of Parvati Basin, NW Himalaya: A Remote Sensing and GIS Based Approach. J. Geol. Soc. India 2021, 97, 165–172. [Google Scholar] [CrossRef]
  22. Malik, A.; Kumar, A.; Kandpal, H. Morphometric analysis and prioritization of sub-watersheds in a hilly watershed using weighted sum approach. Arab. J. Geosci. 2019, 12, 118. [Google Scholar] [CrossRef]
  23. Sakthivel, R.; Jawahar Raj, N.; Sivasankar, V.; Akhila, P.; Omine, K. Geo-spatial technique-based approach on drainage morphometric analysis at Kalrayan Hills, Tamil Nadu, India. Appl. Water Sci. 2019, 9, 24. [Google Scholar] [CrossRef] [Green Version]
  24. Prasad, R.N.; Pani, P. Geo-hydrological analysis and sub watershed prioritization for flash flood risk using weighted sum model and snyder’s synthetic unit hydrograph. Model. Earth Syst. Environ. 2017, 3, 1491–1502. [Google Scholar] [CrossRef]
  25. Prakash, K.; Mohanty, T.; Singh, S.; Chaubey, K.; Prakash, P. Drainage morphometry of the Dhasan river basin, Bundelkhand craton, central India using remote sensing and GIS techniques. J. Geomat. 2016, 10, 122–132. [Google Scholar]
  26. Kumar, A.; Singh, S.; Pramanik, M.; Chaudhary, S.; Maurya, A.K.; Kumar, M. Watershed prioritization for soil erosion mapping in the Lesser Himalayan Indian basin using PCA and WSA methods in conjunction with morphometric parameters and GIS-based approach. Environ. Dev. Sustain. 2022, 24, 3723–3761. [Google Scholar] [CrossRef]
  27. Rahman, M.M.; Zaman, M.N.; Biswas, P.K. Optimization of significant morphometric parameters and sub-watershed prioritization using PCA and PCA-WSM for soil conservation: A case study in dharla River watershed, Bangladesh. Model. Earth Syst. Environ. 2022, 8, 2661–2674. [Google Scholar] [CrossRef]
  28. Malik, A.; Kumar, A.; Kushwaha, D.P.; Kisi, O.; Salih, S.Q.; Al-Ansari, N.; Yaseen, Z.M. The Implementation of a Hybrid Model for Hilly Sub-Watershed Prioritization Using Morphometric Variables: Case Study in India. Water 2019, 11, 1505–1519. [Google Scholar] [CrossRef] [Green Version]
  29. Meshram, S.G.; Sharma, S.K. Prioritization of watershed through morphometric parameters: A PCA-based approach. Appl. Water Sci. 2017, 7, 1505–1519. [Google Scholar] [CrossRef] [Green Version]
  30. Meshram, S.G.; Alvandi, E.; Meshram, C.; Kahya, E.; Fadhil Al-Quraishi, A.M. Application of SAW and TOPSIS in Prioritizing Watersheds. Water Resour. Manag. 2020, 34, 715–732. [Google Scholar] [CrossRef]
  31. CARE-Perú. Modelización Hidrológica de la Cuenca Cañete y Evaluación del Impacto del Cambio Climático; CARE-Perú: Lima, Peru, 2018. [Google Scholar]
  32. Portuguez-Maurtua, M. Aplicación de la geoestadística a modelos hidrológicos en la cuenca del río Cañete. Master’s Thesis, Universidad Nacional Agraria La Molina, Lima, Peru, 2017. [Google Scholar]
  33. Arulbalaji, P.; Padmalal, D. Sub-watershed Prioritization Based on Drainage Morphometric Analysis: A Case Study of Cauvery River Basin in South India. J. Geol. Soc. India 2020, 95, 25–35. [Google Scholar] [CrossRef]
  34. Meraj, G.; Romshoo, S.A.; Yousuf, A.R.; Altaf, S.; Altaf, F. Assessing the influence of watershed characteristics on the flood vulnerability of Jhelum basin in Kashmir Himalaya. Nat. Hazards 2015, 77, 153–175. [Google Scholar] [CrossRef]
  35. Gajbhiye, S.; Mishra, S.K.; Pandey, A. Prioritizing erosion-prone area through morphometric analysis: An RS and GIS perspective. Appl. Water Sci. 2014, 4, 51–61. [Google Scholar] [CrossRef] [Green Version]
  36. Odiji, C.A.; Aderoju, O.M.; Eta, J.B.; Shehu, I.; Mai-Bukar, A.; Onuoha, H. Morphometric analysis and prioritization of upper benue river watershed. Appl. Water Sci. 2021, 11, 41. [Google Scholar] [CrossRef]
  37. Strahler, A. Quantitative Geomorphology of Drainage Basin and Channel Networks. In Handbook of Applied Hydrology; McGraw-Hill: New York, NY, USA, 1964. [Google Scholar]
  38. Horton, R.E. Erosional development of streams and their drainage basins; hydrophysical approach to quantitative morphology. Geol. Soc. Am. Bull. 1945, 56, 275–370. [Google Scholar] [CrossRef] [Green Version]
  39. Nooka Ratnam, K.; Srivastava, Y.; Venkateswara Rao, V.; Amminedu, E.; Murthy, K. Check dam positioning by prioritization of micro-watersheds using SYI model and morphometric analysis remote sensing and GIS perspective. J. Indian Soc. Remote Sens. 2005, 33, 25–38. [Google Scholar] [CrossRef]
  40. Sreedevi, P.D.; Sreekanth, P.D.; Khan, H.H.; Ahmed, S. Drainage morphometry and its influence on hydrology in an semi arid region: Using SRTM data and GIS. Environ. Earth Sci. 2013, 70, 839–848. [Google Scholar] [CrossRef]
  41. Schumm, S.A. Evolution of drainage systems and slopes in badlands at Perth Amboy, New Jersey. Geol. Soc. Am. Bull. 1956, 67, 597–646. [Google Scholar] [CrossRef]
  42. Horton, R.E. Drainage-basin characteristics. Trans. Am. Geophys. Union 1932, 13, 350–361. [Google Scholar] [CrossRef]
  43. Miller, V. A Quantitative Geomorphic Study of Drainage Basin Characteristics in the Clinch Mountain Area Virginia and Tennessee. J. Geol. 1953, 65, 112–113. [Google Scholar] [CrossRef]
  44. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L.; et al. The Shuttle Radar Topography Mission. Rev. Geophys. 2007, 45, 5326–5350. [Google Scholar] [CrossRef] [Green Version]
  45. Amani, M.; Ghorbanian, A.; Ahmadi, S.A.; Kakooei, M.; Moghimi, A.; Mirmazloumi, S.M.; Moghaddam, S.H.A.; Mahdavi, S.; Ghahremanloo, M.; Parsian, S.; et al. Google Earth Engine Cloud Computing Platform for Remote Sensing Big Data Applications: A Comprehensive Review. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 5326–5350. [Google Scholar] [CrossRef]
  46. Mutanga, O.; Kumar, L. Google Earth Engine Applications. Remote Sens. 2019, 11, 18–27. [Google Scholar] [CrossRef] [Green Version]
  47. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  48. Zhao, Q.; Yu, L.; Li, X.; Peng, D.; Zhang, Y.; Gong, P. Progress and Trends in the Application of Google Earth and Google Earth Engine. Remote Sens. 2021, 13, 3778. [Google Scholar] [CrossRef]
  49. Tamiminia, H.; Salehi, B.; Mahdianpari, M.; Quackenbush, L.; Adeli, S.; Brisco, B. Google Earth Engine for geo-big data applications: A meta-analysis and systematic review. ISPRS J. Photogramm. Remote Sens. 2020, 164, 152–170. [Google Scholar] [CrossRef]
  50. Konrad Banachewicz, L.M. The Kaggle Book; Packt: Birmingham, UK, 2022; ISBN 9781801817479. [Google Scholar]
  51. Ismail, M.; Singh, H.; Farooq, I.; Yousuf, N. Quantitative morphometric analysis of Veshav and Rembi Ara watersheds, India, using quantum GIS. Appl. Geomat. 2022, 14, 119–134. [Google Scholar] [CrossRef]
  52. Basnet, K.; Paudel, R.C.; Sherchan, B. Analysis of watersheds in Gandaki province, Nepal using QGIS. Tech. J. 2019, 1, 16–28. [Google Scholar] [CrossRef]
  53. Alencar da Silva Alves, K.M.; Parodi Dávila, M.C.; Zimmermann García, E.D.; Rodrigues de Lira, D.; De Araujo Monteiro, K. Caracterización morfométrica de la cuenca del Salado Bajo, Región de Atacama, Chile. Investig. Geográficas 2021, 62, 90–105. [Google Scholar] [CrossRef]
  54. Gayen, S.; Bhunia, G.S.; Shit, P.K. Morphometric Analysis of Kangshabati-Darkeswar Interfluves Area in West Bengal, India using ASTER DEM and GIS Techniques. 2013. Available online: http://111.93.204.14:8080/xmlui/handle/123456789/582 (accessed on 30 November 2022).
Figure 1. Map showing the location of the Cañete River drainage system.
Figure 1. Map showing the location of the Cañete River drainage system.
Water 15 01053 g001
Figure 2. Methodological diagram for mapping of areas vulnerable to flash floods.
Figure 2. Methodological diagram for mapping of areas vulnerable to flash floods.
Water 15 01053 g002
Figure 3. Sub-basins and drainage network of the Cañete River.
Figure 3. Sub-basins and drainage network of the Cañete River.
Water 15 01053 g003
Figure 4. Map of zones vulnerable to flash floods in the Cañete basin.
Figure 4. Map of zones vulnerable to flash floods in the Cañete basin.
Water 15 01053 g004
Table 1. Linear morphometric parameters and formulas with references.
Table 1. Linear morphometric parameters and formulas with references.
ItemMorphometric ParameterUnitFormula and DefinitionReference
1Stream order (u)DimensionlessHierarchical range[37]
2Stream length (Lᵤ)KmStream length[38]
3Basin length (Lb)KmLb = 1.312 × A0.568
Lb = basin length (km) A = basin area (km2).
[39,40]
Table 2. Morphometric parameters related to area and formulas with references.
Table 2. Morphometric parameters related to area and formulas with references.
ItemMorphometric ParameterUnitFormula and DefinitionReference
1Basin area (A)Km2Estimated in GIS
2Basin perimeter (P)KmEstimated in GIS[41]
3Stream frequency or flow frequency (Fs) ( Km 2 )Fs = ΣNᵤ/A
Nᵤ = total number of stream segments of order “u” and A = basin area (km2)
[42]
4Drainage density (Dd) ( Km 1 )Dd = Σ L/A
L = total stream length; A = basin area
[42]
5Form factor (Ff)DimensionlessFf = A/ L b 2
A = basin area
L b = basin length
[42]
6Circularity ratio (Cr)DimensionlessCr = 4 π A / P 2
A = basin area (km2)
P = basin perimeter (km)
[43]
7Texture ratio (Tr)DimensionlessTr = N l / P
N l = total number of first order streams
P = basin perimeter
[38]
8Elongation ratio (Er)DimensionlessEr = 2 A / π / L b
A = basin area
π = 3.14
L b = basin length
[41]
9Shape factor (Sf)DimensionlessSf   = L b 2 /A[42]
Table 3. Morphometric parameters related to relief and formulas with references.
Table 3. Morphometric parameters related to relief and formulas with references.
ItemMorphometric ParameterUnitFormula and DefinitionReference
1Basin relief (R)MetersR = H − h
R = basin relief, H = maximum elevation in meters
H = minimum elevation in meters
[41]
2 Relief   ratio   ( R r )Dimensionless R r = H/ L b
R r = relief ratio, H = basin relief, L b = basin length
[41]
3Average slope (As)DegreesEstimated in GIS
Table 4. Area, perimeter, and length of each sub-basin.
Table 4. Area, perimeter, and length of each sub-basin.
Sub-BasinPerimeter (km)Area (km2)Length (km)
SB-1197.09943.3064.20
SB-2158.45448.5842.09
SB-3137.14364.6837.42
SB-4145.79419.9640.54
SB-5162.71586.7649.03
SB-6118.83324.5635.02
SB-7183.65596.7149.50
SB-8204.97616.6150.43
SB-930.0131.499.31
SB-10167.51517.6445.66
SB-11281.881175.4772.75
Table 5. Number of streams in each sub-basin by stream order.
Table 5. Number of streams in each sub-basin by stream order.
Stream Order (u)SB-1SB-2SB-3SB-4SB-5SB-6SB-7SB-8SB-9SB-10SB-11
1st order158657410100824
2nd order56241345038
3rd order91000004040
4th order005070701017
Total2915139157211911549
Table 6. Stream length in each sub-basin by stream order.
Table 6. Stream length in each sub-basin by stream order.
Stream Order (u)SB-1
(km)
SB-2
(km)
SB-3
(km)
SB-4
(km)
SB-5
(km)
SB-6
(km)
SB-7
(km)
SB-8
(km)
SB-9
(km)
SB-10
(km)
SB-11
(km)
1st order88.0135.0432.8436.1066.9329.1141.0754.750.0059.13128.42
2nd order31.8029.3911.9230.527.6622.3129.7917.160.0014.1139.91
3rd order54.5510.430.000.000.000.000.0041.790.0028.510.00
4th order0.000.0019.350.0036.670.0029.950.007.680.0078.86
Total174.3674.8664.1166.62111.2651.43100.80113.707.68101.74247.19
Table 7. Morphometric parameters that represent area and form aspects of the sub-basins.
Table 7. Morphometric parameters that represent area and form aspects of the sub-basins.
Sub-BasinDrainage Density
(Dd)
Stream Frequency
(Fs)
Circularity Ratio
(Cr)
Texture Ratio
(Tr)
Elongation Ratio
(Er)
Form Factor
(Ff)
Shape Factor
(Sf)
SB-10.1850.0310.3050.0760.5400.2294.369
SB-20.1670.0330.2250.0500.5680.2533.949
SB-30.1760.0360.2440.0440.5760.2603.840
SB-40.1590.0210.2480.0340.5700.2553.914
SB-50.1900.0260.2790.0430.5580.2444.096
SB-60.1580.0220.2890.0340.5800.2653.779
SB-70.1690.0350.2220.0540.5570.2444.106
SB-80.1840.0310.1840.0490.5560.2424.124
SB-90.2440.0320.4390.0000.6800.3632.752
SB-100.1970.0290.2320.0480.5620.2484.027
SB-110.2100.0420.1860.0850.5320.2224.502
Table 8. Morphometric parameters that represent sub-basin relief aspects.
Table 8. Morphometric parameters that represent sub-basin relief aspects.
Sub-BasinBasin Relief (R)
(m)
Relief Ratio (Rr) Average Slope (As)
(°)
SB-12906.000.04522.590
SB-22061.000.04920.816
SB-33287.000.08827.833
SB-42844.000.07022.800
SB-54127.000.08427.959
SB-63309.000.09419.485
SB-74063.000.08227.023
SB-84156.000.08222.011
SB-92405.000.25834.042
SB-103963.000.08721.225
SB-114420.000.06123.243
Table 9. Preliminary priority rankings.
Table 9. Preliminary priority rankings.
Sub-BasinRRrDdFsCrTrAsErFfSf
SB-1811572272210
SB-21110948410775
SB-37372673993
SB-4981011596884
SB-53549482557
SB-66211103101110102
SB-74783934448
SB-826661158339
SB-910115111111111
SB-105438769666
SB-11192110151111
Table 10. Correlation matrix of the morphometric parameters.
Table 10. Correlation matrix of the morphometric parameters.
RRrDdFsCrTrAsErFfSf
R1.0000.0090.2820.173−0.5640.3730.0820.6360.636−0.636
Rr0.0091.0000.136−0.1550.345−0.7820.255−0.709−0.7090.709
Dd0.2820.1361.0000.3270.0640.1450.5180.2820.282−0.282
Fs0.173−0.1550.3271.000−0.4820.5640.3180.2820.282−0.282
Cr−0.5640.3450.064−0.4821.000−0.6000.227−0.545−0.5450.545
Tr0.373−0.7820.1450.564−0.6001.000−0.1730.8730.873−0.873
As0.0820.2550.5180.3180.227−0.1731.000−0.064−0.0640.064
Er0.636−0.7090.2820.282−0.5450.873−0.0641.0001.000−1.000
Ff0.636−0.7090.2820.282−0.5450.873−0.0641.0001.000−1.000
Sf−0.6360.709−0.282−0.2820.545−0.8730.064−1.000−1.0001.000
Sum1.9910.1002.7552.027−0.5551.4002.1641.7551.755−1.755
Weight (w)0.1710.0090.2370.174−0.0480.1200.1860.1510.151−0.151
Table 11. W S A c p values in each sub-basin.
Table 11. W S A c p values in each sub-basin.
Sub-BasinRRrDdFsCrTrAsErFfSf
SB-11.370.091.181.22−0.100.241.300.300.30−1.51
SB-21.880.092.130.70−0.380.481.861.061.06−0.75
SB-31.200.031.660.35−0.290.840.561.361.36−0.45
SB-41.540.072.371.92−0.241.081.121.211.21−0.60
SB-50.510.040.951.57−0.190.960.370.750.75−1.06
SB-61.030.022.601.74−0.141.202.051.511.51−0.30
SB-70.680.061.890.52−0.430.360.740.600.60−1.21
SB-80.340.051.421.05−0.520.601.490.450.45−1.36
SB-91.710.010.240.87−0.051.320.191.661.66−0.15
SB-100.860.030.711.39−0.330.721.670.900.90−0.90
SB-110.170.080.470.17−0.480.120.930.150.15−1.66
Table 12. Final ranking and sub-basin priority areas.
Table 12. Final ranking and sub-basin priority areas.
Sub-Basin W S A c p   ( + ) W S A c p   ( ) PriorityPriority LevelPriority TypeArea (%)
SB-15.31−0.906.2211Medium15.65
SB-26.751.365.408Medium7.44
SB-34.342.262.082Very high6.05
SB-47.851.816.0410Medium6.97
SB-54.210.453.764High9.74
SB-68.502.715.789Medium5.39
SB-73.840.003.845High9.90
SB-84.42−0.454.887High10.23
SB-94.293.171.121Very high0.52
SB-105.060.904.156High8.59
SB-111.47−1.362.833Very high19.51
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Portuguez-Maurtua, M.; Arumi, J.L.; Stehr, A.; Lagos, O.; Chávarri-Velarde, E.; Rivera-Ruiz, D. Mapping of Areas Vulnerable to Flash Floods by Means of Morphometric Analysis with Weighting Criteria Applied. Water 2023, 15, 1053. https://doi.org/10.3390/w15061053

AMA Style

Portuguez-Maurtua M, Arumi JL, Stehr A, Lagos O, Chávarri-Velarde E, Rivera-Ruiz D. Mapping of Areas Vulnerable to Flash Floods by Means of Morphometric Analysis with Weighting Criteria Applied. Water. 2023; 15(6):1053. https://doi.org/10.3390/w15061053

Chicago/Turabian Style

Portuguez-Maurtua, Marcelo, Jose Luis Arumi, Alejandra Stehr, Octavio Lagos, Eduardo Chávarri-Velarde, and Daniela Rivera-Ruiz. 2023. "Mapping of Areas Vulnerable to Flash Floods by Means of Morphometric Analysis with Weighting Criteria Applied" Water 15, no. 6: 1053. https://doi.org/10.3390/w15061053

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

Article Metrics

Back to TopTop