Next Article in Journal
Trophic Niche Overlap between Invasive and Indigenous Fish in a Northwest Reservoir of China
Previous Article in Journal
Using a Self-Organizing Map to Explore Local Weather Features for Smart Urban Agriculture in Northern Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Delimitation of Small Headwater Catchments and Their Classification in Terms of Runoff Risks

Department of Landscape Water Conservation, Faculty of Civil Engineering, Czech Technical University in Prague, Thákurova 7, Prague 6, 166 29 Prague, Czech Republic
Water 2021, 13(23), 3458; https://doi.org/10.3390/w13233458
Submission received: 26 October 2021 / Revised: 22 November 2021 / Accepted: 30 November 2021 / Published: 6 December 2021
(This article belongs to the Section Hydrology)

Abstract

:
The hydrological similarity of catchments forms a basis for generalizing their hydrological response. This similarity of the hydrological response enables catchments to be classified from numerous perspectives, e.g., hydrological extremes or ecological aspects of catchments. A specific group is formed by so-called “first-order catchments”. This article describes the derivation process of small headwater catchments up to 5 km2 in size on the territory of the Czech Republic. The delimitation is based on the digital terrain model, the stream network, and the water reservoirs. The catchments derived in this way cover 80% of the country. Five mutually independent and sufficiently representative parameters were selected with Principal Components Analysis (PCA), and were used for the cluster analysis performed on two to eight clusters. Clustering Validity Indices (CVI) was used to determine the optimal number of clusters. Subsequently, each generated cluster was assessed for the potential risk of the occurrence of direct runoff, in five classes, on a scale from a moderate degree of risk to a high degree of risk. Six clusters were generated, which is the optimal number in terms of the CVI and their hydrological properties. In this case, 17% of the Czech Republic territory is assessed as lying within a high-risk area, 39% as lying within a medium-risk area, and 24% as lying within a below-average risk area in terms of the occurrence of direct runoff.

1. Introduction

The mutual hydrological similarity of catchments, derived from the similarity of their response to a precipitation event, forms a basis for generalizing their hydrological links, and enables findings to be transferred between catchments. This approach enables catchments to be classified both in terms of their potential impacts on the environment and in terms of their vulnerability to hydrological extremes. Numerous authors have studied the similarities of catchments, their characteristics, and responses from various perspectives. For example, Wagener et al. [1] are motivated in their classification by long-term processes that affect the responses of catchments. A similar principle for the defining, classifying, and sharing the attributes of catchments has been adopted within the CAMELS data set [2,3]. The data sets are created for individual catchments to describe six main classes of attributes at catchment scale: topography, climate, streamflow, land cover, soil, and geology.
The “First-Order Catchments” (FOC), which are also referred to as headwater catchments or as upper catchments by some authors, form a specific group of catchments [4]. FOC are the backbone of the hydrographic network, and are the primary areas for capturing or mitigating flood-related damage. These catchments are often the most environmentally sensitive and the most rapidly evolving parts of many landscapes. They tend to be on the front line of environmental change, and pose the greatest challenges for those involved in land management, policy, and planning [5]. FOC also provide essential ecosystem goods and services for downstream areas [6].
The important parameters affecting the hydrological response from FOC are the properties of the soil and the land cover, morphological characteristics, and precipitation. A large set of characteristics describes morphology characterization, including the average slope, hillslope lengths, and the mean topographic index [7]. The greatest risks in terms of discharges and factors that influence outflow from these catchments are caused by heavy rains and as a consequence of erosion. The precipitation-runoff relationships in the FOC are affected by the speed of the processes. Response to causal precipitation comes relatively quickly, and runoff is affected by intense sub-day precipitation, often of the convectional type. The most frequently-used tool for describing important attributes of precipitation are Intensity-Duration-Frequency (IDF) curves [8], which describe the relationships between rainfall intensity, rainfall duration, and frequency of repetition. The spatial distribution of the respective statistical attributes was studied on a worldwide scale in [9]. Kašpar [10] dealt with sub-daily precipitation, focusing on return period of precipitation using small river basins modelling. The hydrological response is also significantly affected by hyetograph on the local scale. Characteristic shapes of six hour precipitation were defined from long-term measuring with rain gauge and precipitation radar by Muller [11]. On the other hand, the advantage in very small river basins is greater homogeneity of precipitation events, and it is not necessary to use area reduction coefficients [12] to determine the causal precipitation.
Relatively few long-term observations have been made in these small catchments [13], or, for example, the GEOMON network [14]. Damage caused by flash floods and erosion, which are recorded by rescue services or are recorded by authorities can be used for localization, for example, monitoring of soil erosion [15] in the Czech Republic. Regulations on small streams and, potentially, in a catchment area, are more often proposed on basis of a hydrological model of design precipitation than real episodes. A simple method often used by design planners is the SCS-CN method [15], which has been continuously developed and tested [16,17]. By contrast, more complex, physically-based models, e.g., WEPP [18], SMODERP [19], TOPMODEL [20], and others, are used to describe the response of a catchment in a more precise manner.
The main potential risk of flood in the FOC is the occurrence of rapid runoff. Flood risk management based on a catchment scale approach is widely adopted [21,22] with principle being that run-off can be managed most effectively with a combination of measures in the area of the catchment and downstream flood protections [23]. Potential risks associated with the rapid occurrence of runoff in terms of the FOC characteristics that are more sloping have less infiltration capacity with less permeable land use and higher precipitation. The risk can also be assessed in terms of the potential threat to the infrastructure by flood or sediment flux [24].
In addition to the spatial delimitation of small headwater catchments, this article also introduces a set of parameters for classifying them in terms of the potential hydrological response. This article also presents a classification of catchments, and identifies the potential risks of the occurrence of rapid runoff.
The classification of streams is governed by various criteria. Classification based on the profiles of surface water bodies can be mentioned as an example [25]. Hierarchical categorization is also sometimes performed the other way round. First, the catchments of a major waterway are defined, and then they are further subdivided according to the importance of individual sub catchments.
This article is devoted to the derivation of small headwater catchments up to a catchment size of 5 km2 in the Czech Republic territory. For catchments of this size, the use of hydrological models for determining the hydrological characteristics is assumed [26]. Catchments less than 5 km2 in size can be considered so small that they can be viewed as homogeneous areas in terms of precipitation. The catchment classification adopted in the Czech Republic, which applies a hierarchical categorization at four levels, can serve as an example. The first level is represented by the major watercourses (the Elbe, Danube, and Odra Rivers), and these are subdivided at lower levels. Catchments differing in size from one km2 to tens of km2 are grouped in the most detailed fourth category. The classification of four level catchments is presented in [27]. The distribution of precipitation in the Czech Republic has been studied by [10].

2. Materials and Methods

Small headwater catchments (SHC) are the subcategory of FOC, and defined as catchments with a surface area of less than 5 km2 [8] which, at the same time, do not have any tributaries according to the definition of a “first-order catchment” [4]. Small headwater catchments on the territory of the Czech Republic (see Figure 1) were identified according to the methodology described below. Parameters affecting the hydrological response, particularly with respect to the potential risks of rapid runoff due to extreme precipitations, were subsequently assigned to the catchment areas defined in this way. Short-term intensive precipitation events are crucial for potential flash flood risk factors in small catchments.

2.1. Definition of a Catchment

According to [26], SHCs are not only catchments with an area of at least 5 km2, but also include all smaller catchments. This means, for example, that two catchments each with an area of 3 km2 after the confluence do exceed 5 km2, but are still regarded as two separate catchments falling within the SHC category. A total of six categories differing in size have been identified for defining an SHC; see Table 1. The classification of catchments into categories by their size enables comparisons to be made.
Areas smaller than Category 005 can be considered as elementary runoff areas, and they are not assessed as areas of separate catchments. In total, three data sources were used to delimit SHCs: (a) the digital terrain model (DTM); (b) the streams shapes; and (c) water reservoirs. The main input for the definition of SHC was ZABAGED® [28], with a resolution of 5 × 5 m. This DTM was generated on the basis of aerial LiDAR scanning. Due to anthropogenic interventions and changes in the landscape, the axes of streams have changed, and do not correspond with the natural flow paths generated on the digital terrain model itself. By contrast, the axes of streams that form part of ZABAGED® are based on measurements at more detailed scales, and reflect the present-day state of the stream network. In terms of flow direction, these paths are considered to be more accurate and more useful for generating SHCs than the DTM-based flow directions. While derived from the catchments, these stream axes must be included in the solution. The axes of the streams were projected onto DTM. The value of the terrain model pixels through which the axis of the stream passes was reduced by 100 m of altitude. The axes of the streams were therefore included in the generation of the flow direction, which was generated on the basis of DTM. In the next step, the potential places without runoff were discarded from the terrain model modified in this way, and subsequently an accumulation layer was created on the basis of the flow direction. The one-way flow direction tool was used to direct the flow (D8).
For each catchment category (see Table 1), the accumulation layer was reclassified so that the accumulation area values outside the group boundaries would have a NoData value, and the accumulation area corresponding to the given category would then be equal to 1.
The reclassified accumulation rasters generated for each size category were converted into vector lines. The lines were also intersected with the polygons of the water reservoirs included in the DIBAVOD (digital database of water management data). If the end part of the line was located in the water reservoir, the line was shortened to the point of entry of the flow line into the water reservoir. The point was identified by the intersection of the stream line with the polygon representing the respective water reservoir. The endpoints were determined for the modified lines in individual categories. The points identified in this way represent the end profiles of the SHC categories, to which the catchment boundaries were subsequently generated based on the model, taking into account the stream axes.

2.2. Characteristics of Small Catchments

The hydrological response of a small headwater catchment results from its morphological characteristics, soil properties, land use and the precipitation falling over the catchment. It can be assumed that the similar catchments will also have a similar hydrological response. The parameters for the classification of the SHC categories were therefore derived in terms of their potential hydrological response. The morphological characteristics were identified based on the DTM and the streams. These are, in particular, characteristics related to the elevation, the slopes, and the length of the flow paths. In addition, there are several shape coefficients.
The mean catchment width b:
b = A F l m a x
The catchment shape factor alfa:
α = A F l m a x 2
The shape coefficient gra proposed by Gravelius [29]:
g r a = P 2 A π
All three shape coefficients describe the shape of a catchment. The Gravelius coefficient represents the difference in the shape of the catchment from a circle. The mean catchment width factor identifies the extent to which the shape of a catchment diverges from a square, and the catchment shape factor alfa identifies the extent to which the shape of a catchment diverges from the power of maximal flow length path.
The stream network density (SND)parameter is a standard descriptive parameter. This parameter is calculated as the ratio of the total length of the streams to the catchment area:
S N D = L A = t o t a l   s t r e a m _ l e n g t h a r e a
The lag time (Tlag) is a parameter which is affected by the morphology of the catchment, and affects the runoff pattern. Tlag is used to describe the unit hydrograph according to the SCS-CN method [15]. The Tlag value is calculated from the equation [30]
T l a g = F l m a x * 0.8 ( S + 1 ) 0.7 1900 Y
In Equations (1)–(5):
Tlag: lag time (hours)
F l m a x : length of the longest flow path (km)
F l m a x * : length of the longest flow path length (feet)
L : total stream length (km)
Y: catchment mean slope [-]
S*: maximum potential retention height marked (inch)
The direct runoff potential of a catchment is described by the mean value CN. CN integrates information on land use (LU) and soil infiltration properties. In the Czech Republic, CN values were derived within the Strategy for Protection against Negative Impacts of Floods and Erosion Phenomena by Nature-Friendly Measures in the Czech Republic [31].
The final group of parameters consists of precipitation data. As short-term precipitation is the dominant source of runoff in small headwater catchments, six-hour design totals derived from rainfall radars with spatial resolution of 1 × 1 km were applied [10,32]. These data are available at rain.fsv.cvut.cz, and the mean value for each SHC was calculated and used for this study.
All monitored parameters are listed in Table 2. These values describe the mean value, the variance, or maximum value, according to the type of parameter.
In total, there were 28 parameters, which were subsequently tested in all catchment size categories for mutual dependence using regression analysis. The objective was to obtain a set of independent parameters, and to classify the catchments into groups according to their similarity, using cluster analysis.
The delimitation of small catchments, and the assignment and calculation of their characteristics from DTM and the CSC-CN of the catchment were processed using the ESRI environment tool (ArcGIS and ArcGIS Pro).
Descriptive statistics, regression analysis, and principal components analysis (PCA) [33,34] were used for to exclude correlating parameters, and to narrow down the selection of the number of variables for clustering analyse. Cluster analysis was performed using the K-mean method. Clustering validity indices (CVI) [35] were used to determine the optimal number of clusters. Statistical and clustering analyses were performed in the R environment.

3. Results

3.1. SHC Delimitation

Basic data on SHC derived according to the methodology described above are presented in Table 3. As the SHC categories are always derived separately, the resulting catchments overlap between the categories; a smaller resulting catchment may be part of a larger catchment in the higher categories. Therefore, in addition to the categories described above, a group of catchments was generated, in which only the largest catchments were represented. The nested catchments were eliminated. In this way, catchments smaller than 5 km2 situated in the monitored area in the Czech Republic have been preserved. The resulting group of catchments generated in this way is from now on referred to as the “Set of the Largest Catchments” (SoLC), and is also listed in Table 3. For clarity, and to give an idea of the representation of each size category in the resulting SoLC group, the table contains data on the number of elements in the respective category that form part of the SoLC.

3.2. Choice of Parameters

The parameters for each catchment in all size categories were derived according to Table 2. For the purposes of cluster analysis, representative and independent parameters were sought in the first step. Dependent parameters must be discarded. The search for the degree of agreement between the monitored parameters was performed both for individual categories (including SoLC) and for all the catchments together. In terms of the groups of dependent parameters, the individual categories do not differ from each other. Therefore, the relationships between the monitored parameters are similar for all size categories. The visualized parameter agreement is shown in Figure 2.
Figure 2 shows the parameters that are mutually correlated in the marked groups A–D. Principal Component Analysis (PCA) was used for selection of representative parameters. PCA was used for (i) all parameters in first round and (ii) by the groups (A–D) in second.
(i) It was found that the first two principal components of PCA of all data do not explain even 50% of the variability. The first five components: 75% of the variability; the remaining six components: 80% of the variability (Figure 3).
The contribution of individual parameters, especially in the case of the first two substitute variables, does not exceed 9% and 13%, respectively (see Figure 4). PCA analysis did not provide a clear choice of explanatory variables (see Figure 4).
(ii) The first two dimensions of the PCA analysis explain for groups A–D: 80% of the already embodied 98 variability. The effect of the first two components is given in more detail using biplots in Appendix C. The PCA analysis did not show any significant factors that should not be overlooked. The final selection of representative parameters was chosen by combining the significance of the element in the PCA analysis, with regard to the usability of the selected parameters in hydrological modelling.
The six parameters were selected. Group C is represented by CN_mean. Slope_mean can also characterise group C. Moreover, CN_mean represents the retention capacity. Precipitation data have a separate group—group D. The parameter Tag from group A was chosen in the same way. This parameter interpreted the input to hydrological models more than the slightly better explanatory variables Fl_len_max or Fl_len_STD. A close correlation (R2 > 90) was found between all precipitation data. Rainfall with a return period of 20 years (P20) was chosen with regard to the possible hydrological response. With shorter repetition times, thus precipitation totals, the hydrological model may not lead to runoff, a repetition time of 100 years is too extreme, and the time series are relatively short to derive this repetition. Three parameters from group B were selected. Parameters SND and Alpha, around which subgroups are formed, together correlate with the parameter Fl_len_noStream_STD. All three of these parameters are explanatory variables of the third PCA component. Fl_len_noStream_STD explains (20%), SND (12%), and Alpha (7%) of this component. SND and Alpha do not correlate with each other, and together have the same explanatory weight as Fl_len_noStream_STD. SND and Alpha are direct value of the parameter and not of the Standard Deviation Parameter, and are more relevant to assign a risk level to these parameters.
Five parameters that can be considered independent and sufficiently representative (groups A–D) are as follows:
  • Precipitation P20—there is a significant correlation between individual values of six-hour precipitation events due to the derivation of these data;
  • The CN mean, which represents a group of several other parameters. The CN value shows agreement both with the slope and with the altitude;
  • Lag time (Tlag)—this parameter characterizes the shape of the runoff hydrograph, and therefore the peak flow volume;
  • Stream Network Density (SND)—this parameter represents the density of the streams and amount of the stream network;
  • Shape coefficient (Alpha)—This parameter incorporates the characteristics of the flow path length and of the catchment shape.

3.3. Distribution of Parameters

In order to classify the catchments into groups in terms of their potential response, it is necessary to compare the distribution of the classification parameters among the individual catchment categories. If the selected classification parameters were to show a different distribution for each group of catchments, it would imply that different size categories manifested a different type of hydrological response to a precipitation event. The objective is, therefore, to compare the differences between the respective size categories. The distribution of the five parameters selected for the individual size categories is presented in Figure 5. Visualization in the form of histograms has been chosen, showing the distribution of values and box plots, displaying the median value, the interquartile range, and outliers. The number of classes is the same for all histograms. The numbers of elements in the displayed histograms are normalized so that the distribution of the catchment representation in individual classes can be compared.
As the results in Figure 6 clearly show, the differences between the catchment size categories are not significant. There is an obvious difference in the SND specification for the 0.5 and 1 km catchment size, where catchments without a hydrographic network occur to a greater extent. In the other size categories, catchments without a hydrographic network occur significantly less often. The second clear difference is in the Tlag value, which also partially includes the catchment size; the concentration times are generally higher in the larger categories.
A cluster analysis for each category separately would not provide information significantly different from SoLC in terms of runoff potential. Cluster analysis was therefore performed for SoLC. Each catchment size categories are represented at least 20% in SoLC.
Cluster analysis (K-mean) was performed by five representative parameters for two to eight clusters, with 25 initial training points. The formation of groups of catchments is shown in Figure 7. Groups are marked with letters. If a group is generated only by separation from a previously formed group, a numeric designation is added.
The groups generated during gradual clustering are characterized below. Geographical clustering is shown in Figure 8. The centroids of each cluster for the built clusters are listed in Appendix A.
  • Two Clusters—When the first two clusters are formed, Group A is generated, which is characterized by a somewhat higher CN number with lower precipitation volumes. Group B is characterized by higher precipitation volumes and a higher CN value (Figure 8a);
  • Three Clusters—Group A is divided mainly in terms of the shape characteristics of the catchment, the density of the stream network, and the lag time (Tlag) (Figure 8b);
  • Four Clusters—Group B1, which is characterized by lower precipitation volumes while maintaining a lower CN value, and, by contrast, Group B2 with higher precipitation totals and, at the same time, a higher CN value, are separated from Group B (Figure 8c);
  • Five Clusters—Group A1 is predominantly divided on the basis of the lag time. The resulting group A12 is characterized by a significant Tlag time, while the initial characteristics of Group A1 are fairly preserved in Group A11. Group A11 and group A12 defined in this way are preserved even after the catchments have been subdivided into more clusters (Figure 8d);
  • Six Clusters—A completely new Group D is generated, which is characterized by a relatively high SND value and, at the same time, relatively low precipitation totals, while maintaining a relatively high CN value. Group D, generated in this way, is preserved even after the catchments have been subdivided into more clusters (Figure 8e);
  • Seven Clusters—Group B2, which is characterized by relatively high precipitation volumes, has been noticeably divided. Together with some of the catchments of Group A2, it generates a new Group C, which is characterized by higher precipitation totals and, at the same time, higher CN values. Some of the catchments of initial Group B2 and some of the catchments of Group B1 generate Group B3, which retains parameters similar to those of initial Group B2. However, the number of catchments in initial Group B2 is so small that the group is renamed B3 (Figure 8f);
  • Eight Clusters—The newly generated Group C is regrouped into Groups C1 and C2. The newly arising Group C1 is also supplemented by some of the catchments of Group A2 which, similar to the initial Group C, are characterized by higher precipitation volumes and the CN value. Consequently, Group C 1 differs from Group C2 by the difference in the SND and Alpha parameters (Figure 8g).
The gradually generated groups of catchments are characterized by the mean values of the five selected parameters.
The Davies-Bouldin index was used to determine the optimal number of clusters using Clustering Validity Indices (CVI) (Figure 9). The optimal number of clusters from this is five. A suitable number of clusters is six, from which completely new class D arises (described in more detail in discussion).
The generated clusters of the catchments are further classified according to the potential risks of the occurrence of rapid runoff. In terms of the effect of the parameters on the risks associated with the rapid occurrence of runoff, the following holds for individual parameters:
  • SND—the higher the value, the denser the network of permanent streams, the more likely the potential runoff is to be concentrated in these paths, where runoff is expected. A higher value means a lower level of risk;
  • Tlag—the greater the lag time, the lower the anticipated peak flows;
  • Alpha—the more complex the shape of the catchment is, the longer the flow paths and thus the lower the culmination.
  • CN mean—the lower the CN mean value, the greater the retention capacity in the catchment, and the lower the risk of a potential threat;
  • 6 h rain—the more intensive the rain is, the higher the risk of a potential runoff response.
The mean value, which is considered as a medium level of risk, was calculated for individual parameters in the SoLC category. The degree of risk was identified for the individual parameters relative to this mean value of the respective parameter. For each parameter value corresponding to the centroid of the individual clusters, the quotient was identified using this mean value, which yielded the degree of risk of each parameter in the cluster. Combinations of five parameters where a negative assessment prevails were then considered risky, and vice versa. The overall degree of risk was divided into five categories, from low risk to high risk, as described below.
  • Low risk—the combination of the parameters of a potential runoff response implies a small risk. in terms of rapid direct runoff affecting the catchment. These areas appear to be unproblematic in terms of a rapid response, and there is no need to implement any measures;
  • Decreased risk—the combined parameters of a potential runoff response imply a rather small risk in terms of rapid runoff affecting the catchment. These areas appear to be unproblematic in terms of a rapid response, and there is low need to implement any measures;
  • Medium risk—the combined parameters of a potential runoff response are average, and an average degree of risk is assumed in terms of rapid runoff affecting the catchment;
  • Increased risk—the combined parameters of a potential runoff response imply a rather higher degree of risk in terms of rapid runoff affecting the catchment;
  • High risk—the combined parameters of a potential runoff response imply a great risk in terms of rapid runoff affecting the catchment. In these areas, a more detailed survey and more detailed monitoring of potential negative impacts of rapid runoff need to be carried out.
The parameter values for identifying the degree of risk are listed in Table 4.
The classification of the catchment groups formed by cluster analysis according to the degree of risk is presented in Appendix A. Groups with two–eight clusters are included. The classification of the catchment groups produced by cluster analysis, according to the degree of risk of the developing clusters, is presented in Appendix B. The geographical representation of the degree of risk of six clusters is shown in Figure 10. Group A2 forms a high-risk group, Groups A11, B1, and B2 form a medium risk group, and Groups A12 and D form a group with a “decreased level of risk”.

4. Discussion

Catchment classifications are more frequently performed in experimental hydrology. In the expanding CAMELS database, the catchments are also classified. In the case of CAMELS, there is a longer list of parameters [2]. Unlike the selected parameters presented here, the list of parameters is extended to include the hydrological data of long-term balances and parameters that have a more significant effect on long-term runoff and on other components of the balance. In most cases, larger catchments are investigated. Long-term observation campaigns are conducted considerably less frequently in small catchments than in larger ones.
Input data with different spatial resolutions were used to generate the boundaries of the catchments and their properties. The catchment boundaries were delimited on the basis of a terrain model with a resolution of 5 × 5 m, which is sufficiently detailed even for the delimitation of small catchments in size category 005. Based on the terrain model, other morphological characteristics were derived with the same resolution. If the D8 method were used with lower resolution, the generation of the catchment boundaries could be affected, mainly in the smallest category.
Within the investigated catchment size categories from 0.5 to 5 km2, a significantly higher number of catchments are grouped in the first category (005). The parameters entering the cluster analysis do not differ significantly in terms of the distribution of values between the categories. The smaller catchments are also form part of larger categories, and together they form SoLC (Set of Largest Catchments), where at least 20% of the number of catchments in each category is represented. The total area of the upper catchments included in SoLC is 63 thousand km2, which amounts to 80% of the area of the Czech Republic (78 thousand km2).
A total of 28 parameters were considered to affect hydrological response. Some of the parameters create groups of mutually corelated parameters. PCA provided only a partial answer to the selection of the most suitable parameters from whole dataset. That the selection of representative elements is based in combination of PCA, pair correlations of parameters, and usability of the parameter in hydrological models. Selection of the parameters is discussed in the results, which is most logical. Slope parameters, stream slope and the altitude are considered to be in relation to CN. This correlation corresponds to the land use in mountainous areas, which are mostly more sloping, and are mostly forested. Another group of mutually correlated parameters are the areal parameters of the catchments (size, area, flow path length, runoff accumulation), which that are correlated with Tlag. SND and Alpha correlate with standard deviation of sheet flow length. However, this pair (SND and Alpha) is more often used in hydrological models. Precipitation events forming a separate group of parameters are mutually strongly correlated. However, they do not correlate significantly with the other parameters.
In SHC, short-term rains, in particular, are crucial in terms of response and, potentially, in terms of the risk of increased flows. A specific flood event and a potential threat occur when the current conditions in the catchment are combined with the course of the precipitation. Short-term torrential rains are events particularly difficult to forecast. However, two differently classified catchments that are in the same initial condition, and are subjected to the same precipitation will show a different response to similar precipitation events.
The subsequent cluster analysis of the catchments in terms of their hydrological response risk shows that, based on the selected parameters. SHC (small headwater catchments) can be subdivided and can split into two groups. In these groups, Categories A2 and B2 are gradually separated, and their parameters cause them to fall into the group with a risk of an increased runoff threat due to torrential rains. CVI and Davies-Bouldin [35] were used to determine the optimal number of clusters. Other indices were tested (Calinski Harabasz and Ratkowsky Lance) [35] without clear conclusion about the most appropriate number of clusters. The optimal number of clusters from Davies-Bouldin is five. Groups C and D arising from A and B previously generated sub-groups in six and more clusters. Group D that arises in six clusters has the main characteristic of high SND value—catchment with high stream density (Appendix A), and should be preserved in the final grouping. Group C is later separated from two groups (A2 and B2) high risk catchments. The formation of two groups, C1 and C2, from Group C and partially from Group A2, with a total of eight groups, no longer provides new information in terms of potential threats.
In relation to the hydrological response type and the possible level of risk, the view of the final number of clusters is determined not only by the CVI result, but also by values of centroids. Statistical processing did not show any important parameter that would not be used in the analysis. A comparison with the recorded events [36] can be used in feature work to verify the risk classification of SHC. Classification of catchments according to parameters has a practical impact for potential prioritization in terms of the implementation of protective measures.

5. Conclusions

The Czech Republic is an example of a country where headwater catchments form a significant part of the territory. In the case of SoLC, 80% of the Czech Republic is covered in this category. SHC is the space for the primary accumulation of rainwater. At the same time, SHC tend to be more affected by rapid runoff, which, in turn, reduces the availability of water in the ‘headwater catchment’ area. The SHC classification presented here specifies the degree of threat, and reveals the probable hydrological response of each cluster of catchments. Regarding the first five clusters, the primary subdivision into two groups, A and B that arise early as the first two clusters. Within Group A and Group B, two subgroups are gradually formed, which are moderately high-risk groups. A sufficiently explanatory classification of SHC uses seven clusters, where a very low-risk group (D) is generated from the elements of Groups A and B; and Group C, with a high degree of risk, is separated from Groups A2 and B2. With six clusters covering the area of the Czech Republic, this approach places 17% of the territory in the high-risk category, 39% in the medium-risk category, and 24% in the below-average risk category. A total of 20% of the territory of the Czech Republic is not assessed, as it does not fall into the SHC. Within the Czech Republic, the agriculturally exploited areas of south Moravia, the Bohemian-Moravian Highlands, and north-west Bohemia can be considered as regions with a more significant degree of risk. The medium risk areas are mainly mountainous regions, with the exception of the Krušné Mountains, which are a low-risk area. The lower risk areas include the foothills, with the exception of the Podorlicko region, the Elbe River lowlands, the Brdy Mountains and West Bohemia.
The classification of small headwater catchments in terms of the threat by potential torrential rains is one of the perspectives. Another use of the spatial delimitation of these catchments could be for a subsequent classification, for example, in terms of the availability of water for irrigation or in terms of the application of other adaptation measures due to anticipated changes in climate.

Funding

This research was supported by the Czech Ministry of Agriculture (grant No. QK1910029).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data is available in the web application rain.fsv.cvut.cz (accessed on 26 October 2021).

Conflicts of Interest

The author declares no conflict of interest and 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.

Appendix A

Table A1. A table presenting the development of the parameter values in the cluster centroids during the progressive formation of two to eight clusters. The table also shows the number of elements in the respective in each group (Count) of the total number of SHCs in SoLC, the in SoLC. The value is also expressed as a percentage.
Table A1. A table presenting the development of the parameter values in the cluster centroids during the progressive formation of two to eight clusters. The table also shows the number of elements in the respective in each group (Count) of the total number of SHCs in SoLC, the in SoLC. The value is also expressed as a percentage.
Number of Clusters Group Count % of SoLC SND Tlag Alpha CN Mean 6 h Rain
2A22,85856.90.623.044.0775.2247.59
2B17,29443.11.593.563.2860.4753.87
3A1823920.51.154.865.5072.6247.86
3A217,78344.30.542.503.4175.0548.24
3B14,13035.21.593.293.1058.9054.31
4A1648416.11.015.045.7273.2547.17
4A214,68936.60.452.533.4175.7147.01
4B1971224.21.303.483.1253.9149.80
4B2926723.11.702.953.4870.6358.21
5A11532613.30.913.096.1874.8047.47
5A1238269.51.077.074.1069.0247.48
5A213,42733.40.452.513.2675.6547.16
5B1893322.21.323.203.1153.4850.02
5B2864021.51.712.923.4470.5158.44
6A11457811.40.663.126.3374.6647.66
6A1235708.91.017.194.1169.0347.52
6A212,03130.00.352.523.2675.8146.99
6B1771819.21.163.273.0652.2950.57
6B2752618.71.322.953.3870.2159.83
6D472911.82.492.833.7870.3847.72
7A1138299.50.703.216.5675.2648.00
7A1233338.31.037.334.1369.2248.12
7A2932723.20.402.663.3976.3644.77
7B1650116.20.943.123.4054.1946.95
7B3466811.61.643.442.7256.1160.04
7C842421.00.832.573.3875.0656.06
7D407010.12.672.833.9470.6349.15
8A1131277.80.723.236.8475.0748.07
8A1228177.01.007.674.1668.6647.96
8A2751518.70.492.893.8276.3543.39
8B1568314.20.883.163.5053.9846.36
8B3468411.71.553.392.6054.6558.89
8C1704117.50.362.292.8375.4552.08
8C2541213.51.403.094.0174.2358.42
8D38739.62.662.833.7169.0547.98

Appendix B

Table A2. Table of the development of the classification of catchment groups produced by cluster analysis according to the degree of risk during the progressive formation of two to eight clusters.
Table A2. Table of the development of the classification of catchment groups produced by cluster analysis according to the degree of risk during the progressive formation of two to eight clusters.
Number of Clusters Degree of Risk of Individual Parameters Mean Risk
SND Tlag Alfa CN Mean 6 h Rain
2A1.681.070.921.090.951.14increased risk
B0.650.921.140.881.070.93decreased risk
3A10.900.670.681.050.950.85low risk
A21.931.301.091.090.961.27high risk
B0.650.991.200.861.080.96medium risk
4A11.020.650.651.060.940.86decreased risk
A22.311.291.091.100.931.35high risk
B10.790.941.190.780.990.94decreased risk
B20.611.111.071.031.160.99medium risk
5A111.141.060.601.090.940.97medium risk
A120.970.460.911.000.940.86decreased risk
A22.311.301.151.100.941.36high risk
B10.791.021.200.780.990.96medium risk
B20.601.121.081.021.161.00medium risk
6A111.561.050.591.080.951.05medium risk
A121.030.450.911.000.940.87decreased risk
A22.921.301.141.100.931.48high risk
B10.891.001.220.761.010.98medium risk
B20.791.111.101.021.191.04medium risk
D0.411.150.991.020.950.91decreased risk
7A111.481.020.571.090.951.02medium risk
A121.000.450.901.010.960.86decreased risk
A22.591.231.101.110.891.38high risk
B11.101.041.100.790.930.99medium risk
B30.630.951.370.811.190.99medium risk
C1.241.271.101.091.111.16high risk
D0.391.150.951.030.980.90decreased risk
8A111.431.010.551.090.961.01medium risk
A121.030.430.901.000.950.86decreased risk
A22.111.130.981.110.861.24high risk
B11.181.031.070.780.921.00medium risk
B30.670.961.440.791.171.01medium risk
C12.831.431.321.101.041.54high risk
C20.741.060.931.081.160.99medium risk
D0.391.161.011.000.950.90decreased risk

Appendix C

Figure A1. Biplots of the first two dimensions of the PCA analysis for Groups AD.
Figure A1. Biplots of the first two dimensions of the PCA analysis for Groups AD.
Water 13 03458 g0a1

References

  1. Wagener, T.; Sivapalan, M.; Troch, P.; Woods, R. Catchment Classification and Hydrologic Similarity. Geogr. Compass 2007, 1, 901–931. [Google Scholar] [CrossRef]
  2. Addor, N.; Newman, A.J.; Mizukami, N.; Clark, M.P. The CAMELS data set: Catchment attributes and meteorology for large-sample studies. Hydrol. Earth Syst. Sci. 2017, 21, 5293–5313. [Google Scholar] [CrossRef] [Green Version]
  3. Alvarez-Garreton, C.; Mendoza, P.A.; Pablo Boisier, J.; Addor, N.; Galleguillos, M.; Zambrano-Bigiarini, M.; Lara, A.; Puelma, C.; Cortes, G.; Garreaud, R.; et al. The CAMELS-CL dataset: Catchment attributes and meteorology for large sample studies-Chile dataset. Hydrol. Earth Syst. Sci. 2018, 22, 5817–5846. [Google Scholar] [CrossRef] [Green Version]
  4. Strahler, A.N. Quantitative analysis of watershed geomorphology. Eos Trans. Am. Geophys. Union 1957, 38, 913–920. [Google Scholar] [CrossRef] [Green Version]
  5. Křeček, J.; Haigh, M.J. Reviewing the Contexts of Headwater Control. Environ. Reconstr. Headwater Areas 2000, 9–24. [Google Scholar] [CrossRef]
  6. Gumbo, A.D.; Kapangaziwiri, E.; Chikoore, H.; Pienaar, H. Assessing water resources availability in headwater sub-catchments of Pungwe River Basin in a changing climate. J. Hydrol. Reg. Stud. 2021, 35, 100827. [Google Scholar] [CrossRef]
  7. Beven, K.; Heathwaite, L.; Haygarth, P.; Walling, D.; Brazier, R.; Withers, P. On the concept of delivery of sediment and nutrients to stream channels. Hydrol. Process. 2005, 19, 551–556. [Google Scholar] [CrossRef]
  8. Chen, C.L. Rainfall intensity-duration-frequency formulas. J. Hydraul. Eng. 1983, 109, 1603–1621. [Google Scholar] [CrossRef]
  9. Courty, L.G.; Wilby, R.L.; Hillier, J.K.; Slater, L.J. Intensity-duration-frequency curves at the global scale. Environ. Res. Lett. 2019, 14. [Google Scholar] [CrossRef]
  10. Kašpar, M.; Bližňák, V.; Hulec, F.; Müller, M. High-resolution spatial analysis of the variability in the subdaily rainfall time structure. Atmos. Res. 2021, 248, 105202. [Google Scholar] [CrossRef]
  11. Muller, M.; Kašpar, M.; Bližňnák, V.; Kavka, P. Precipitation intensity during heavy rains in various altitudes. In Rainfall in Urban and Natural Systems, Proceedings of the 10th International Workshop on Precipitation in Urban Areas (UrbanRain15), Zurich, Switzerland, 1–5 December 2015; ETH—Institute of Environmental Engineering: Zürichc, Switzerland, 2015. [Google Scholar] [CrossRef]
  12. Asquith, W.H.; Famiglietti, J.S. Precipitation areal-reduction factor estimation using an annual-maxima centered approach. J. Hydrol. 2000, 230, 55–69. [Google Scholar] [CrossRef] [Green Version]
  13. Zumr, D.; Dostál, T.; Devátý, J. Identification of prevailing storm runoff generation mechanisms in an intensively cultivated catchment COST CA18135-Fire in the Earth System: Science & Society View project Payments for Ecosystem Services (Forests for Water) View project. Artic. J. Hydrol. Hydromech. 2015. [Google Scholar] [CrossRef] [Green Version]
  14. Oulehle, F.; Fischer, M.; Hruška, J.; Chuman, T.; Krám, P.; Navrátil, T.; Tesař, M.; Trnka, M. The GEOMON network of Czech catchments provides long-term insights into altered forest biogeochemistry: From acid atmospheric deposition to climate change. Hydrol. Process. 2021, 35. [Google Scholar] [CrossRef]
  15. Mishra, S.K.; Singh, V.P. Soil Conservation Service Curve Number (SCS CN) Methodology; Springer: Dordrecht, The Netherlands, 2003; Volume 42. [Google Scholar]
  16. Caletka, M.; Michalková, M.Š.; Karásek, P.; Fučík, P. Improvement of SCS-CN Initial Abstraction Coefficient in the Czech Republic: A Study of Five Catchments. Water 2020, 12, 1964. [Google Scholar] [CrossRef]
  17. Soulis, K.X. Soil Conservation Service Curve Number (SCS-CN) Method: Current Applications, Remaining Challenges, and Future Perspectives. Water 2021, 13, 192. [Google Scholar] [CrossRef]
  18. Laflen, J.M.; Elliot, W.J.; Flanagan, D.C.; Meyer, C.R.; Nearing, M.A. WEPP-predicting water erosion using a process-based model. J. Soil Water Conserv. 1997, 52, 96. [Google Scholar]
  19. Kavka, P.; Weyskrabová, L.; Zajíček, J. Comparison of 1D and 2D modelling with soil erosion model SMODERP. In Proceedings of the European Geoscience Union General Assembly, Vienna, Austria, 7–12 April 2013; pp. 7–12. [Google Scholar]
  20. Metcalfe, P.; Beven, K.; Hankin, B.; Lamb, R. A new method, with application, for analysis of the impacts on flood risk of widely distributed enhanced hillslope storage. Hydrol. Earth Syst. Sci. 2018, 22, 2589–2605. [Google Scholar] [CrossRef] [Green Version]
  21. Dadson, S.J.; Hall, J.W.; Murgatroyd, A.; Acreman, M.; Bates, P.; Beven, K.; Heathwaite, L.; Holden, J.; Holman, I.P.; Lane, S.N.; et al. A restatement of the natural science evidence concerning catchment-based ‘natural’ flood management in the UK. Proc. R. Soc. A Math. Phys. Eng. Sci. 2017, 473, 20160706. [Google Scholar] [CrossRef] [Green Version]
  22. Werritty, A. Sustainable flood management: Oxymoron or new paradigm? Area 2006, 38, 16–23. [Google Scholar] [CrossRef]
  23. Lane, S.N. “Natural flood management,” Wiley Interdiscip. Rev. Water 2017, 4, e1211. [Google Scholar] [CrossRef] [Green Version]
  24. Jáchymová, B.; Krása, J.; Dostál, T. and Bauer, M. Can Lumped Characteristics of a Contributing Area Provide Risk Definition of Sediment Flux? Water 2020, 12, 1787. [Google Scholar] [CrossRef]
  25. European Parliament, Directive 2000/60/EC of the European Parliament and of the council of 23 October 2000 establishing a framework for community action in the field of water policy. Off. J. Eur. Commun. L 2000, 327, 172. Available online: https://ec.europa.eu/environment/water/water-framework/index_en.html (accessed on 23 October 2000).
  26. Sedlák, K. TNV—Stream Management. 2010. Available online: http://eagri.cz/public/web/file/104408/TNV_75_2102.pdf (accessed on 15 November 2021).
  27. Strouhal, L.; Kavka, P.; Weyskrabová, L.; Pavel, M. Classification of small watersheds in Czech Republic for hydrological modelling. Water Manage. Tech. Econ. Inf. J. 2018, 60, 16–21. [Google Scholar] [CrossRef]
  28. Šilhavý, J.; Čada, V. New automatic accuracy evaluation of altimetry data: DTM 5G compared with ZABAGED® altimetry. In Lecture Notes in Geoinformation and Cartography; Springer: Berlin, Germany, 2015; Volume 211, pp. 225–236. [Google Scholar]
  29. Gravelius, H. Grundriß der Gesamten Gewässerkunde, Band 1: Flußkunde. Compend. Hydrol. 1914, I, 265–278. [Google Scholar]
  30. Chow, V.T.; Mays, L.W.; York, N.; San, S.L.; Bogota, F.A.; Lisbon, C.; Madrid, L.; City, M.; New, M.M.; San, D.; et al. Applied Hydrology; McGraw Hill Book Company: New York, NY, USA, 1988; ISBN 0-07-010810-2. [Google Scholar]
  31. Drbal, K. Strategie ochrany před negativními dopady povodní a erozními jevy přírodě blízkými opatřeními v České republice. Vodohospodářské Tech. Inf. 2016, 58, 4. [Google Scholar] [CrossRef]
  32. Bližňák, V.; Kašpar, M.; Müller, M. Radar-based summer precipitation climatology of the Czech Republic. Int. J. Climatol. 2018, 38, 677–691. [Google Scholar] [CrossRef]
  33. Pearson, K. LIII. On lines and planes of closest fit to systems of points in space. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1901, 2, 559–572. [Google Scholar] [CrossRef] [Green Version]
  34. Singh, P.K.; Kumar, V.; Purohit, R.C.; Kothari, M.; Dashora, P.K. Application of Principal Component Analysis in Grouping Geomorphic Parameters for Hydrologic Modeling. Water Resour Manag. 2009, 23, 325–339. [Google Scholar] [CrossRef]
  35. Dimitriadou, E.; Dolničar, S.; Weingessel, A. An examination of indexes for determining the number of clusters in binary data sets. Psychometrika 2002, 67, 137–159. [Google Scholar] [CrossRef] [Green Version]
  36. Žížala, D.; Kapička, J.; Novotný, I. Monitoring Soil Erosion of Agricultural Land in Czech Rebublic and Data Assessment of Erosion Events from Spatial Database. In Proceedings of the Proceedings from International Conference Soil—The non-Renewable Environmental Resource, Brno, Czech Republic, 7–9 September 2016; Mendel University: Brno, Czech Republic, 2016; pp. 354–370. [Google Scholar]
Figure 1. Study area of the Czech Republic with elevation model (CUZK—Zabaged®).
Figure 1. Study area of the Czech Republic with elevation model (CUZK—Zabaged®).
Water 13 03458 g001
Figure 2. Visualization of the correlation coefficients between individual parameters for all catchments regardless of size category. A negative correlation is shown in red, and a positive correlation is shown in blue. The stronger the correlation between two parameters, the darker and larger the symbol. Parameters in correlation are grouped. The main groups (AD) are separated by a black line. A description of the parameters is provided in Table 2.
Figure 2. Visualization of the correlation coefficients between individual parameters for all catchments regardless of size category. A negative correlation is shown in red, and a positive correlation is shown in blue. The stronger the correlation between two parameters, the darker and larger the symbol. Parameters in correlation are grouped. The main groups (AD) are separated by a black line. A description of the parameters is provided in Table 2.
Water 13 03458 g002
Figure 3. Percentage of explained variances with PCA, individual proportion (black) and cumulative (blue).
Figure 3. Percentage of explained variances with PCA, individual proportion (black) and cumulative (blue).
Water 13 03458 g003
Figure 4. Contribution of variables for first four dimensions of PCA for all parameters. (a) to first, (b) second, (c) third, (d) fourth substitute variables.
Figure 4. Contribution of variables for first four dimensions of PCA for all parameters. (a) to first, (b) second, (c) third, (d) fourth substitute variables.
Water 13 03458 g004
Figure 5. Correlation of five selected parameters. The correlation coefficients are numerically shown below the diagonal, and the correlations are visually presented above the diagonal, as in Figure 2. The selected parameters can be regarded as independent variables.
Figure 5. Correlation of five selected parameters. The correlation coefficients are numerically shown below the diagonal, and the correlations are visually presented above the diagonal, as in Figure 2. The selected parameters can be regarded as independent variables.
Water 13 03458 g005
Figure 6. The left column displays the box plots showing the median value, the interquartile range, and the 25 or 75 quartiles of the values. The right column displays the frequency histograms in each class. The values in the histograms are normalized to achieve a comparable number of values in each catchment size category.
Figure 6. The left column displays the box plots showing the median value, the interquartile range, and the 25 or 75 quartiles of the values. The right column displays the frequency histograms in each class. The values in the histograms are normalized to achieve a comparable number of values in each catchment size category.
Water 13 03458 g006
Figure 7. Sankey diagram of the development of clusters of SoLC. The widths of the bands are proportional to number of elements. The diagram also shows transfer of individual catchments and regrouping together with the growing number of clusters. The basic subdivision is already evident when two main groups (A, B) are formed. The groups generated up to five clusters, which are subdivisions of A or B. Group D, which is a combination of previously formed groups, is generated in the six clusters step. Group C arises when the number of clusters reaches seven and eight as combinations of Groups A2 and B2.
Figure 7. Sankey diagram of the development of clusters of SoLC. The widths of the bands are proportional to number of elements. The diagram also shows transfer of individual catchments and regrouping together with the growing number of clusters. The basic subdivision is already evident when two main groups (A, B) are formed. The groups generated up to five clusters, which are subdivisions of A or B. Group D, which is a combination of previously formed groups, is generated in the six clusters step. Group C arises when the number of clusters reaches seven and eight as combinations of Groups A2 and B2.
Water 13 03458 g007
Figure 8. Geographical representation of the development of the groups during the formation of two to eight (ag) clusters.
Figure 8. Geographical representation of the development of the groups during the formation of two to eight (ag) clusters.
Water 13 03458 g008
Figure 9. CVI (Davies-Bouldin).
Figure 9. CVI (Davies-Bouldin).
Water 13 03458 g009
Figure 10. Degree of risk of the Small Headwater Catchments for case of six clusters—optimal count.
Figure 10. Degree of risk of the Small Headwater Catchments for case of six clusters—optimal count.
Water 13 03458 g010
Table 1. SHC categories. Each category is characterized by area.
Table 1. SHC categories. Each category is characterized by area.
CategoryFrom km2Up to km2
0050.30.7
0100.71.3
0201.72.3
0302.73.3
0403.54.5
0504.55.5
Table 2. List of parameters that enter the SHC classification. Names of each parameter included type of parameter; where relevant, it also contains a description of the type of statistics. The statistics (standard deviation (STD), mean, maximum) are calculated in GIS. Parameters 1–12 are derived from DTM, or in combination with the stream axes, where the slope characteristics (15–18) were identified in the original DTM layer before the stream receded. Parameters 20–24 were obtained from Equations (1)–(4).
Table 2. List of parameters that enter the SHC classification. Names of each parameter included type of parameter; where relevant, it also contains a description of the type of statistics. The statistics (standard deviation (STD), mean, maximum) are calculated in GIS. Parameters 1–12 are derived from DTM, or in combination with the stream axes, where the slope characteristics (15–18) were identified in the original DTM layer before the stream receded. Parameters 20–24 were obtained from Equations (1)–(4).
No.GroupName and DescriptionNo.GroupName and Description
1Basic parametersPerimeter15Slope parametersSlope_mean—average slope
2Area16Slope_STD—standard deviation
3Elevation mean—average catchment elevation17Slope_stream_mean–average slope of the streams
4Elevation STD—deviation in elevation describing the flatness of the catchment18Slope_stream_STD—deviation in the slope of the streams
5Flow accumulation parametersFl_acc_mean—average flow accumulation19Morphological parametersMedium width—mean catchment width
6Fl_acc_STD—deviation in flow accumulation20Shape coefficient Alpha
7Flow length parametersFl_len_max—maximum flow path length21Gravelius coefficient—shape coefficient
8Fl_len_mean—average flow path length22Tlag—lag time
9Fl_len_STD—standard deviation in the length of the flow path23Stream parametersTotal stream length—total length of the streams
10Fl_len_noStream_max—maximum length of the surface runoff flow path24SND—stream network density
11Fl_len_noStream_mean—average length of the surface runoff flow path25Mean six-hour design precipitationP2—rainfall with returnPeriod (RP) 2 yr
12Fl_len_noStream_STD—standard deviation in the surface runoff flow path26P10—rainfall with RP 10 yr
13SCS-CN parametersCN_mean—average CN number for the whole catchment27P20—rainfall with RP 20 yr
14CN_STD—standard deviation28P100—rainfall RP period 100 yr
Table 3. Number of catchments and the total area of catchments in each category. The first two columns describe the number of elements in the respective category, and the last two columns indicates the number of elements in SoLC.
Table 3. Number of catchments and the total area of catchments in each category. The first two columns describe the number of elements in the respective category, and the last two columns indicates the number of elements in SoLC.
CategoryNumber of ElementsTotal Area (km2)Representation of Elements in the SoLC
Number%km2%
00572,62137,63216,89423772712
01031,28733,04610,9073511,03818
02011,56024,179393834805113
030653020,289218733665511
040543122,610227142908614
050395720,479395710020,47832
SoLC40,15463,031
Table 4. Values of individual parameters used to express the degree of risk in relation to the mean values of the parameters.
Table 4. Values of individual parameters used to express the degree of risk in relation to the mean values of the parameters.
Risk CoefficientLow RiskDecreased RiskMedium RiskIncreased RiskHigh Risk
<0.85<0.95<1.05<1.15>1.15
SND1.191.091.030.980.88
Tlag3.753.433.263.102.77
Alpha4.293.923.733.553.17
CN mean79.2072.3168.8765.4258.54
6 h rain42.7547.7850.3052.8157.84
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kavka, P. Spatial Delimitation of Small Headwater Catchments and Their Classification in Terms of Runoff Risks. Water 2021, 13, 3458. https://doi.org/10.3390/w13233458

AMA Style

Kavka P. Spatial Delimitation of Small Headwater Catchments and Their Classification in Terms of Runoff Risks. Water. 2021; 13(23):3458. https://doi.org/10.3390/w13233458

Chicago/Turabian Style

Kavka, Petr. 2021. "Spatial Delimitation of Small Headwater Catchments and Their Classification in Terms of Runoff Risks" Water 13, no. 23: 3458. https://doi.org/10.3390/w13233458

APA Style

Kavka, P. (2021). Spatial Delimitation of Small Headwater Catchments and Their Classification in Terms of Runoff Risks. Water, 13(23), 3458. https://doi.org/10.3390/w13233458

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