Next Article in Journal
Glacial Change and Its Hydrological Response in Three Inland River Basins in the Qilian Mountains, Western China
Previous Article in Journal
3D GIS Platform for Flood Wargame: A Case Study of New Taipei City, Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Capability of Natural Flood Management Approaches in Groundwater-Dominated Chalk Streams

1
School of Geography and Environmental Science, The University of Southampton, Southampton SO17 1BJ, UK
2
Romsey Office, The Environment Agency, Romsey SO51 8DU, UK
*
Author to whom correspondence should be addressed.
Water 2021, 13(16), 2212; https://doi.org/10.3390/w13162212
Submission received: 22 June 2021 / Revised: 20 July 2021 / Accepted: 23 July 2021 / Published: 13 August 2021
(This article belongs to the Section Hydrology)

Abstract

:
This study aims to address the gap in the Natural Flood Management (NFM) evidence base concerning its implementation potential in groundwater-dominated catchments. We generated a typology of 198 chalk catchments using redundancy analysis and hierarchical clustering. Three catchment typologies were identified: (1) large catchments, (2) headwater catchments with permeable soils, and (3) catchments with impermeable soils and surfaces (urban and suburban land uses). The literature suggests that natural flood management application is most effective for catchments <20 km2, reducing the likelihood of significant flood mitigation in large catchments. The relatively lower proportion of surface runoff and higher recharge in permeable catchments diminishes natural flood management’s likely efficacy. Impermeable catchments are most suited to natural flood management due to a wide variety of flow pathways, making the full suite of natural flood management interventions applicable. Detailed groundwater flood maps and hydrological models are required to identify catchments where NFM can be used in a targeted manner to de-synchronise sub-catchment flood waves or to intercept runoff generated via groundwater emergence. Whilst our analysis suggests that most chalk groundwater-dominated catchments in this sample are unlikely to benefit from significant flood reductions due to natural flood management, the positive impact on ecosystem services and biodiversity makes it an attractive proposition.

1. Introduction

Natural Flood Management (NFM) represents a paradigm of management strategies that aim to improve a catchment’s resilience to prolonged and/or heavy rainfall by restoring, enhancing, or altering a catchment’s natural hydrological and morphological characteristics [1]. Strategies aim to increase interception and infiltration by reducing rapid runoff generation, increasing catchment water storage, and slowing overland channel flows [2,3]. Flood peak reductions of up to 30–40% have been attributed to NFM methods, such as storage pond networks, tree shelter belts, channel realignment, leaky dams, and winter cover crops in rural and urban settings [2,4,5,6,7,8]. Moreover, NFM can offer co-benefits, such as improvements to water quality, reductions in soil degradation, and enhanced biodiversity [1,9,10]. Consequently, the demand for NFM implementation has increased due to its perception as a relatively low cost, low maintenance flood mitigation solution that protects and maintains hydrological and biological function of the rivers in which it is implemented [1,11,12]. Despite this growing demand for NFM interventions, the NFM evidence base consistently cites groundwater-dominated river systems, like chalk streams, as a gap in the knowledge due to hydrological differences, meaning the current evidence may not be directly applicable [11,13,14].
A small number of chalk streams can be found in Northern France and Russia, but the 224 chalk streams found in the UK account for the majority of this river type found globally [15]. Catchments underlain with chalk bedrock are predominantly groundwater dominated [16]. This means that large quantities of water are stored and transmitted through the chalk bedrock before being released into river channels as baseflow, generating a stable annual river regime with a lack of spate conditions and relatively small variations between high and low flows [17]. It can take months for changes in rainfall inputs to manifest as changes in the river regime [18], giving chalk streams a strong seasonal flow oscillation in conjunction with annual rainfall patterns in the UK. Discharge rises slowly over the winter, peaking in late April after heavier winter rainfall, and gradually recedes over the summer, with lowest flows in autumn after generally drier summers. Groundwater flooding occurs when rainfall recharge causes the groundwater level to rise, in turn increasing groundwater inputs into river systems and causing groundwater emergence in topographic low points, winterbournes, and activating springs [19,20].
NFM functions by enhancing a catchments’ natural ability to absorb shocks from storms by storing water in the catchment and then slowly releasing it [21]. A successful NFM scheme reduces discharge at at-risk locations by extending a flood duration and reducing the peak discharge (and water level) at a point at any given time [13,22]. This is achieved by increasing interception and infiltration, slowing overland channel flows, and manipulating channel and catchment surface roughness [2,3]. NFM is a process-based approach, meaning that each NFM scheme is designed by matching flow pathways, landscape features, and sources of flood waters that contribute to peak flows to specific NFM interventions that tackle those issues. Because of this, NFM interventions can be categorised according to the processes that they manipulate: (1) the reduction of rapid runoff generation, (2) increasing catchment water storage, and (3) strategies to reduce the conveyance of water downstream [11,13,14]. This is an important distinction because it allows specific sources of flood waters to be linked to specific solutions in the process of NFM design. For example, many arable fields suffer from sediment loss due to excess rapid overland flow. In this case, NFM interventions would focus on reducing runoff inputs to the river channel by intercepting and storing water. Sediment can be trapped and stored and surface roughness increased (slowing the flow of surface runoff) by using winter cover crops or across-slope tillage. The choice of method would be dictated by the specific soil properties of the site as well as management preferences. As such, NFM schemes are often tailor-made to each application scenario, and it is important that the features and water-transfer processes of the catchment where implementation is proposed are first established.
In contrast, the vast majority of the NFM evidence base is founded on research conducted in surface water-dominated catchments where floods are caused by the convergence of multiple surface-runoff inputs to the river channel [11,13,23]. As a result, the application of NFM in groundwater-dominated catchments has been highlighted as a one of the key knowledge gaps in the NFM evidence base, in acknowledgement that the flow pathways and processes present in groundwater-dominated catchments are significantly different to those found in catchments with fluvial floods. This can be demonstrated by the difference in Base Flow Index (BFI) values for the River Lambourne (a typical chalk stream) and the Belford Burn catchment in Northumberland, which has been used as an NFM experimentation catchment for many years. BFI measures the proportion of total channel flow contributed by groundwater sources. Belford Burn has a BFI of 0.313 [24], and the River Lambourne has a BFI of 0.98 [25]. Therefore, groundwater catchments typically transfer a large proportion of water below the ground surface. Of the 12 main NFM techniques summarised by Lane (2017) [13], eight of them rely on managing surface water. However, if most of the water transfer throughout the catchment occurs below the ground surface, these measures may have limited effect (Figure 1). Because NFM interventions are tailored to fit the specific sources and flood water, it follows that NFM schemes in groundwater catchments will focus on different combinations of interventions than those found in surface-runoff-dominated catchments (Figure 1) depending on the key properties of chalk catchments. It is therefore important to identify the morphological and hydrological features that affect groundwater recharge and the production of overland flows (increasing the probability of channel bank exceedance) to guide future NFM strategies in catchments dominated by groundwater processes.
This study is intended as a screening process for NFM in groundwater catchments by grouping catchments according to hydrological properties and matching them to NFM interventions that specifically tackle these flow pathways. To do this, we quantify the relationships between hydrological variability and key morphological characteristics that are amenable to NFM for 198 catchments with chalk bedrock in the Southeast of England. We use these results to classify the catchments and infer flow pathways and to make suggestions for the most appropriate NFM strategies for these river basins. Because NFM schemes must be designed on a case-by-case basis for the greatest effectiveness, this is intended as a broad-scale screening process to narrow down options and not as a comprehensive guide for choosing NFM interventions in all chalk stream catchments.

2. Materials and Methods

2.1. Study Area and Catchment Selection

Chalk groundwater-dominated catchments were first identified using topographic and geological datasets. Topographic catchment boundaries were provided by the Centre for Ecology and Hydrology National River Flow Archive (NRFA) [26]. A bedrock map (BGS 625k) from the British Geological Survey [27] was used to identify catchments that were underlain by chalk bedrock. To reduce uncertainty due to differing groundwater transfer processes [28] and because of the global rarity of chalk stream river systems, we focused on chalk groundwater-dominated river systems, excluding limestones or Permo-Triassic sandstones. We selected catchments with at least ≥70% of the catchment within the chalk bedrock and with gauging stations within 5 km downstream of the chalk bedrock (determined via the buffer tool in ArcGIS and visual inspection). Using these criteria, a total of 198 catchments were available for analysis, located predominantly in the Southeast of England (Figure 2).

2.2. Data Analysis

To relate hydrological variability to catchment morphological characteristics, we compiled and analysed the covariation among four key hydrological variables and twenty-one variables quantifying the physical catchment properties. Details and data sources are found in Table 1.

2.2.1. Hydrological Variables

All hydrological variables were derived from the Centre for Ecology and Hydrology UK National River Flow Archive (NRFA) [26]. Average discharge (Qmean) and maximum recorded discharge (Qmax) are among the most basic hydrological metrics used to characterise the hydrological regime. The base flow index (BFI) [25] measures the proportion of streamflow from groundwater contributions. This is relevant to NFM implementation because high BFI values are related to stable flows and limited surface-runoff generation and contributions to streamflow. Conversely, the Richard–Baker flashiness index (RBFI) [29] measures the sensitivity of a streamflow to rainfall inputs. Therefore, high RBFI values indicate significant sources of runoff resulting in spate conditions. As mentioned previously, understanding the processes by which water is transferred throughout the catchments is essential to NFM design. Chalk streams are often characterised and defined by features such as long times to peak, high base flow domination, and a lack of spate conditions [15]. However, large ranges in BFI and RBFI values within the chalk stream sample demonstrate that there is major spatial variation in chalk stream properties. These metrics will help characterise chalk catchments according to flow variability, average flows, the short-term response of flow to sharp bursts of rainfall, groundwater contributions to the flow regime, and (indirectly) the rate of recharge and runoff production. This helps inform which NFM interventions will be most suitable based on whether water transfer occurs predominantly above or below the ground surface. Details of index calculations are provided in Table 1.

2.2.2. Physical Catchment Properties

We selected 21 physical variables that directly influence hydrology and that have explained variations in catchment hydrology in other studies [29,30,31]. This included the percentage cover of 9 land uses (arable land, broad leaved and coniferous woodland, grassland, heathland, urban, suburban, marshland, and inland rock) from the 2015 Land Cover Map (LCM 2015) [32]. Heathland, grassland, and marsh land covers are made up of sub-categories of each classification, respectively (i.e., all separate grassland types are combined into a new, homogenous grassland classification). Combinations of sub-categories to generate the three new classes were done as specified in Appendix A of the LCM2015 documentation to generate more optimum data distributions. The hydrological impact of land cover is well documented in the literature (Table 1) and is relevant to NFM because land cover directly impacts flow pathways and, consequently, the choices of NFM interventions that may be suitable. Additionally, land use can inform the amount of space available for NFM application to inform what is possible. Large quantities of storage ponds, controlled flood plain zones, and afforestation may not be suitable NFM interventions in regions of arable land due to the need to preserve economically productive space, whereas other interventions, like winter cover crops, hedge edges, and no-till farming, may be more suitable. Percentage cover of 6 soil types (bypass flow common, bypass flow uncommon, bypass flow very uncommon, bypass flow variable, slowly permeable, and impermeable) were classified according to the hydraulic conductivity in the Hydrology of Soil Types (HOST) [24]. This is because hydraulic conductivity can directly dictate the efficiency of aquifer recharge and surface-runoff generation. Soil-saturated hydraulic conductivity has been identified as one of the key features that dictates chalk stream river regime and flood response [31]. It is therefore important to understand the key hydrological processes in the catchment that will in turn be useful to inform the choice of NFM interventions. Topographic catchment shape and drainage density were included because they influence the rate of propagation of water through the catchment and river channel network [33,34] (p. 304). Bedrock transmissivity and a proxy for aquifer abstraction rates were included because they provide information about how water propagates through the sub surface of the catchment, which in turn influences streamflow response [35,36]. Transmissivity values of the chalk bedrock vary between 230 m2/day to 2600 m2/day [36]. Catchments where transmissivity is low have reduced rates of recharge, subsurface water transfer, and reduced groundwater contributions to channel flows relative to regions with high-transmissivity bedrock [36]. The hydrological influence of each physical catchment variable is included in Table 1 as well as details of their calculation.

2.3. Statistical Analysis

We used redundancy analysis (RDA) [47] to characterise and explain the variation in hydrological properties in relation to physical properties. RDA is an asymmetric ordination method whereby the variation in one set of (explanatory) variables is used to directly explain the variation in another set of (response) variables. Essentially a combination of multiple regression and principal component analysis, RDA generates a matrix of the fitted values of all response variables, which are then subjected to a principal component analysis. Linear combinations of explanatory variables (physical catchment properties) that best explain the variation in the response variables (hydrological properties) are sought by the model in successive order. The main advantages of using RDA over principal component analysis in this context are that variance in the response variables is attributed to the explanatory variables, and presence or absence of a relationship between specific x and y variables can be tested [47]. It can therefore be used to directly link hydrological traits of chalk streams to the presence or absence of specific physical features, enabling more targeted NFM scheme design.
All analyses were undertaken in the R environment (version 3.6.1) [48]. After compiling the variables as outlined in Table 1, all variables were transformed using a Box Cox transformation and were centred and scaled to reduce the influence of outliers and place variables on a common scale for the RDA model [49]. To identify the explanatory variables that best describe the co-variation of the hydrological variables, a preliminary redundancy analysis (RDA) model was initially run using the rda function of the vegan package in R [50]. This global model, using all 21 catchment variables, was subjected to model validation. To achieve a parsimonious model containing important physical explanatory variables, the global model was subjected to a forward stepwise procedure using the forward.sel function of the Packfor R package [51] (p. 48). This process selects the model with the combination of variables with the highest R2 and p value [52] (p. 178). The remaining physical catchment variables after the stepwise procedure are those that directly correlate with and explain the variation in hydrological regime in chalk streams whilst maintaining the highest explanatory power (R2). The resulting forward-selected RDA model was then subjected to a validation process, including ANOVA tests of the significance of the relationships identified in the parsimonious model and the number of significant axes in the model using 1000 permutations. The covariance of the selected physical catchment variables in the model was ascertained using the variance inflation factor (VIF). Catchments were then plotted as a points on a biplot according to their RDA coordinates.

Catchment Classification

Cluster analysis was used to group catchments with similar river regimes according to catchment RDA coordinates. Catchments with similar combinations of physical catchment property variables are located near each other on an RDA biplot (have similar RDA coordinates). Before subjecting the data to cluster analysis, the data’s tendency for clustering was established using the Hopkins Statistic with the get_clust_tendency function of the factoextra package in R [53]. Four hierarchical clustering methods were compared for their suitability to the dataset according to two distance measures (cophenetic correlation and the Gower statistic). Cophenetic correlation measures the degree of agreement between the original unmodeled pairwise distances and the pairwise distances in the dendrogram. A high positive correlation indicates that pairwise distances have been preserved [54]. The Gower statistic is the sum of squared differences between the dissimilarity matrix and the cophenetic distance and is calculated using dendrogram hierarchical partition [55,56]. The clustering algorithms compared were single, complete, and average linkage agglomerative clustering and Ward’s minimum variance clustering [55]. The optimum number of clusters was established via a Mantel correlation. Here, the original distance matrix is compared to binary matrices computed from the dendrogram being cut at multiple different levels for different numbers of cluster allocations [55]. The optimum number of clusters is where the Mantel correlation is highest.
The uncertainty in the allocated clusters was estimated using silhouette widths. This measures the degree of membership of an object to its allocated cluster by comparing the average distance of an object to all other objects in the same cluster, to the average distance between it, and all the objects in the next closest cluster [57]. Accordingly, high silhouette width values indicate that that catchment has a high degree of membership in that group. Catchments with negative silhouette width values can be assumed to have been misclassified.
Once the clustering allocations were validated using the pairwise distance measures, the catchment groupings were mapped and used to describe the variation within chalk stream hydrology and flow pathways, and the potential for the application of NFM for each group was assessed according to its specific physical and hydrological qualities.

3. Results

The stepwise redundancy analysis revealed that the combinations of the following catchment properties significantly explained the co-variation of hydrological variables (p ≤ 0.001): area, uncommon bypass flow, impermeable soils, slowly permeable soils, station elevation, form factor, and urban land use (Table 2). The adjusted R2 values, representing the proportion of hydrological variance, which is explained by the physical catchment properties, were 0.682 for the global model and 0.675 for the reduced model. Only the first two axes were used for analysis (axes 1: 77.4%, axes 2: 35.9%), as the third axis and beyond were statistically insignificant. Variance inflation factor (VIF) values (Table 3) revealed that the remaining 9 variables have very little collinearity (VIF values < 5). Despite the fact that the variables of uncommon bypass flow and slowly permeable soils were mildly collinear, they were both retained because removal of either caused a substantial drop in the explanatory power of the model. The 15 catchment variables removed via the forward selection process did little to reduce the Adjusted R2 value, suggesting that the discarded variables mostly generated noise.

Clustering and Cluster Validation

Catchments were grouped into three clusters using the average linkage hierarchical clustering algorithm. The data’s tendency for clustering was established prior to this via the Hopkins statistic (0.752; significant = >0.5 and <1). Average linkage hierarchical clustering was selected out of the four algorithms because it returned the highest cophenetic correlation and smallest value for the Gower distance (Table 4). Despite four clusters being identified as the optimum number of clusters for this dataset according to the Mantel correlation (Figure 3a), three clusters were used because the fourth was not associated with any vectors on the biplot, making it difficult to interpret. Silhouette widths were used to quantify and map the uncertainty in catchment group allocations (Figure 4c). An average silhouette width of 0.50 indicates that catchments were generally classified correctly, with all misclassifications occurring in group 1 (Figure 3b). Silhouette widths below 0.31 were considered uncertain due to being close to the decision boundary, and catchments with silhouette widths below 0 were misclassified [52] (p. 70). A drop in average silhouette width from 0.53 with four clusters to 0.5 with three clusters indicates that group cohesion is slightly reduced by this decision. To mitigate against this, the nine misclassified catchments (Figure 3b) were reclassified and displayed as their nearest neighbour alternative classification in Figure 4b.

4. Discussion

The explanatory power of the parsimonious RDA model is within the acceptable range for RDA models that characterise hydrological variation [58,59,60]. The model has been used to group catchments with similar physical catchment properties and river regime properties using their coordinates from the RDA plot. Thus, the cluster classifications were used to infer catchment typology and hydrological conditions, including rates of aquifer recharge and rapid runoff generation. The cluster classifications are used as a screening process to identify which NFM interventions can theoretically be applied or ruled out of each catchment typology based on the dominant physical catchment properties identified by that group as shown on the RDA plot (Figure 4).
It is acknowledged that this approach relies heavily on the partitioning of catchments into groups. Hence, care was taken to ensure that uncertainty in cluster allocation was reduced by checking the data’s tendency to cluster. The number of clusters and the clustering algorithm that made the most statistical sense were selected. However, uncertainty in cluster allocation cannot be removed entirely, and 13.6% of the sample had uncertain group allocations (Figure 4c). It is understood that these catchments share traits that would allow them to be comfortably grouped in to two of the three groups and are therefore understood to be intermediate catchments that do not easily fit in to a single classification. In these cases, the combined traits of the two catchment classifications should be taken into account when selecting NFM options for these catchments (alternative groupings for highly uncertain catchments are provided in Appendix A). Catchments that have been reclassified can be confidently allocated to their group.

4.1. Group 1: Large Catchments

Group 1 catchments have high Qmean and Qmax values that are explained by their large size (Figure 4). These are the largest catchments within the sample, with topographic catchment areas ranging from 108 km2 to 1459 km2. Qmax has a strong negative correlation to station elevation, demonstrating that catchments with higher discharges tend to be those closer to sea level and further downstream, resulting in large water accumulation. Qmean, Qmax, and catchment area have weak correlations (shown by orthogonality) with all other variables in the model (Figure 4a), rendering it difficult to comment further.
Previous research suggests that NFM becomes less effective at reducing flooding in catchments greater than 20 km2 [11,61]. Furthermore, evidence shows that increasing the area impacted by NFM measures does not always increase the gains in flow attenuation [13,41]. The lack of substantial evidence for NFM implementation at a larger catchment scale has become a barrier to the general uptake of NFM [62]; however, it must be acknowledged that very few catchment-scale NFM schemes have been implemented [13,41]. The NFM evidence base at this scale is mostly provided by risk-based predictive models [63], which introduces computation limitations to catchment-scale NFM research and their subsequent findings [64]. Nonetheless, it is generally accepted that NFM is most effective for (and possibly limited to) catchments < 20 km2.
Metcalfe et al. (2018) [64] argued that a large proportion of the benefits of NFM come from its ability to de-synchronise sub-catchment flood waves, which theoretically works at every catchment scale [65]. De-synchronisation is achieved by implementing NFM schemes at a small scale in carefully chosen sub-catchments that are designed to attenuate individual sub-catchment flood waves, reducing the likelihood of flood wave synchronisation. Whilst it is recognised that the flood peak reductions from this process are often small, it can be enough to prevent bank overtopping in events up to 1:100 year scale [66]. Additionally, this method is only suitable in situations where the configuration of the sub-catchments currently causes synchronised flood waves. Dixon et al. (2016) [67] illustrate that, in some cases, NFM placement can synchronise flood waves that were previously staggered, resulting in an overall increase in flood risk. Where desynchronisation of flood waves is the goal, experimental modelling is required in the design process and prior to implementation of an NFM scheme to mitigate such risks [64,68].
As none of the catchments found in group 1 are <20 km2, the evidence is lacking to support the effective use of NFM within these river basins. However, this does not prevent these larger catchments from being broken down into multiple sub-catchments or for NFM opportunities to be identified and applied locally for point-source flooding. An example of this could be to implement bunds on a sloped field to prevent field runoff from flooding an adjacent road. In specific situations, it may be suitable to investigate the use of de-synchronising flood waves to reduce the incidence of bank overtopping. Such areas would be downstream river sections where multiple sub-catchments contribute to flood waves, and relatively small changes in river levels alter the risk of bank overtopping [64,66]. Extensive research and considerable time are required to design and implement effective NFM strategies for de-synchronisation [68], so it is only recommended where combined flood waves are known to be a problem. Low drainage density of permeable chalk catchments, however, limits the opportunities for de-synchronisation.

4.2. Group 2: Permeable Catchments

These catchments are the stereotypical chalk streams: small catchments that are dominated by groundwater recharge and influxes of groundwater that dominate the river regime via baseflow. Group 2 catchments have river regime variability closely related to high BFI values and are associated with a large proportion of uncommon bypass flow soil types and higher station elevations (interpreted here as a higher proportion of headwater catchments). Uncommon bypass flow is a classification of soils from the HOST soil classification that describes thin (aquifer within 2m), permeable, and unconsolidated soils with micro and macro pores [24]. These soils have a vertical saturated hydraulic conductivity of >10 cm day−1 [69], allowing highly efficient recharge via vertical drainage into the chalk aquifer. This is linked to reduced runoff generation [70], limiting the potential for NFM implementation because a large proportion of the water transfer throughout the catchment occurs below the ground surface. Due to the lack of significant runoff, NFM treatments that aim to reduce runoff generation or store surface water within the catchment are highly unlikely to have a significant effect on streamflow. Examples of such interventions include winter cover crops, changes in tillage practices, lengthening drainage pathways, planting across slopes, online and offline storage ponds, wetlands, and controlled flood zones [13]. What remains are NFM techniques that focus on in-channel interventions to reduce downstream conveyance [13], such as channel realignment, sustainable urban drainage systems (SuDS), de-culverting covered river channels, and increasing in-channel, riparian, and marginal vegetation [6,71]. Whilst in-channel river restoration schemes such as these can significantly reduce flood peaks in small catchments [8,72], it must be emphasised that restoration schemes and in-channel interventions deliver the best results as part of a suite of other NFM measures [73]. Therefore, whilst beneficial, in-channel measures alone are unlikely to deliver the optimum impact of NFM interventions.
Despite a lack of surface-runoff generation in highly groundwater-dominated catchments, such as those in group 2, surface water can occur due to groundwater emergence. In this case, the water table rises to intersect with the ground surface, forming static pools of water in topographic low points called turloughs [74] and intermittently flowing river channels called winterbournes [75] and can activate springs in weak points and fractures in the chalk [20,76]. These phenomena generally occur after large quantities of prolonged rainfall. Previous groundwater-emergence floods have been recorded after rainfall events that double and triple the long-term averages [74,77,78,79]. Under such conditions, NFM has been shown to be far less effective because engineered and natural stores of water, such as soil storage, storage ponds, log jams, and groundwater stores, become full and are overwhelmed [11,13]. The presence of turloughs, winterbournes, and springs demonstrate that water transfer throughout the catchment is not uniform across space and is dictated by the location, size, and activation of hydrogeological features such as fractures [19,80,81,82]. As a result, NFM schemes in these catchments will require in-depth, local knowledge to design them, with sensitivity to the local hydrogeological features. We suggest that NFM interventions in this group of catchments be focused in areas where groundwater emergence occurs as winterbournes and springs to intercept or store flood flows, particularly where the resulting surface water causes disruption or damage to property. For example, in-channel interventions, such as log jams and online storage ponds, could be installed in in the path of known winterbournes and springs to intercept flows when the water table is high. The potential benefits of this kind of NFM installation have not currently been tested.
For these catchments, it is suggested that NFM will have diminished effectiveness due to a lack of significant surface runoff. This limits suitable NFM interventions to mostly instream channel modifications to reduce rapid conveyance, which will be less effective compared to a full suite of NFM interventions [13,73]. In many cases, it may be more economical to map and maintain local knowledge of groundwater emergence for flood-risk mapping. Previously, major flood damages were caused by urbanisation of forgotten, dormant winterbournes and springs, which then activate under heavy rainfall [78]. These mapping efforts would also be instrumental in designing and locating potential NFM schemes near hydrogeological features, such as springs and winterbournes. Therefore, effective flood mitigation measures are associated with tracking water table heights, flood warnings, flood mapping, and making a concerted effort to understand spatial changes in hydrogeology. Long-term groundwater emergence and flood-risk maps are required for appropriate flood planning, for identifying locations for NFM schemes, and for reducing flood risk by restricting building and development on areas at risk of groundwater emergence.

4.3. Group 3: Less Permeable Catchments

Group 3 catchments are associated with high RBFI values (Figure 4), which is best explained by higher incidence of impermeable and slowly permeable soils, the presence of urban land use, and higher values of form factor (indicating more circular catchment shapes). These are the chalk catchments with the largest proportions of surface runoff. Higher RBFI values describe greater flow variability and steeper and higher magnitude discharge peaks. The relationship between soil permeability and rapid runoff generation is well known and documented, where impermeable soils and surfaces generate greater quantities of surface runoff and quick-flow catchment pathways [24,31,70,83,84]. Work previously performed on chalk stream catchments by Ascott et al. (2017) [31] concluded that impermeable superficial deposits, such as those found in group 3 catchments, slow the vertical conveyance of water into the aquifers, reducing recharge and groundwater dominance in the river regime and flood response. Inversely, reduced rate of recharge and absorption of rainfall sub-surface will increase surface-runoff generation.
As demonstrated by Lane (2017) [13], a large proportion of NFM interventions work by manipulating and intercepting surface-runoff pathways [61]. By virtue of a greater quantity of surface runoff, the full suite of NFM strategies are viable in group 3 catchments, including reduction of rapid runoff generation through soil management and increased catchment roughness, increasing catchment water storage using storage ponds, and reducing the conveyance of water downstream with in-channel and river restoration strategies. This allows many different potential combinations of NFM intervention for optimising results [73], meaning that this group of chalk catchments is the most suitable for the application of NFM in the study sample.

4.4. Applications of NFM in Chalk Catchments

According to the typology of catchments generated via redundancy analysis, three chalk catchments in the UK (Yeading Brook West at North Hillingdon, Catchwater at Withernwick, and the River Dour at Crabble Mill; Appendix A) inherently have the physical and hydrological features best suited to the current range of NFM measures. The study sample accounts for 198 of the estimated 224 chalk streams in the UK. In the study sample, only 25 chalk catchments are <20 km2, and of these, 22 are classified as permeable catchments, limiting them to mostly in-channel NFM interventions and highly targeted NFM schemes downstream of hydrogeological features. Overall, this suggests that implementing NFM in chalk groundwater-dominated catchments is likely to have sub-optimal results compared to other catchment types [3,12,85].
The findings of this analysis do not necessarily negate the use of NFM methods in chalk catchments. Where appropriate, de-synchronising sub-catchment flood waves can be implemented for larger catchments classed as impermeable by applying NFM at the local scale. NFM should also be considered on the local scale in areas where groundwater emergence as springs and winterbournes cause disruptions or for any other known source of surface water or runoff. Additionally, the environmental benefits of NFM on water quality, aesthetic improvements, reductions in soil erosion, and biodiversity are uncontested [62]. There are proven benefits for river managers and catchment partnerships in chalk streams regions to implement NFM for river restoration, water quality, and biodiversity improvements. The application of NFM is a cross-disciplinary, collaborative process, so the benefits of increased communication and the formation of catchment partnerships have been argued to be another of the co-benefits [61]. Wingfield et al. (2019) [62] argued that building the evidence base for NFM to irrefutability could take decades and that by hesitating to implement NFM due to a lack of current evidence, potential benefits are lost. Minor flood benefits gained through river restoration will likely have minimal effect for large-scale storms but may reduce the incidence of small-scale nuisance floods [66,67]. Therefore, river restoration is advantageous to flood reduction but should not be considered the main objective unless supported by local, detailed analysis. Additionally, even in groups 1 and 2 catchments, there are likely small-scale point-source locations of pluvial flooding, like flooding of fields or roads, and groundwater emergence that could be combatted with small-scale NFM schemes. However, implementation of such schemes will require better reporting and mapping of local sources of flooding.

5. Conclusions

The results of a redundancy analysis model were used to generate a typology of chalk catchments, resulting in three groupings according to broadly similar river regimes and physical catchment properties. Using these classifications as a screening tool, the likely effectiveness of applying NFM in each of these groups was discussed. The first class (group 1) is grouped by virtue of their larger size, meaning they are likely unsuited to NFM due to NFM being most effective for catchments <20 km2. It is acknowledged that these catchments can be broken down into smaller sub-catchments, possibly to de-synchronise sub-catchment flood waves. De-synchronisation is only recommended in conjunction with hydrological modelling prior to NFM design. It should also be considered that flood reductions from de-synchronising flood waves are most effective for smaller nuisance floods rather than larger flood events. Permeable catchments (group 2) are associated with smaller headwater catchments and high-permeability soils. NFM interventions are less suited due to a low proportion of surface-runoff processes. Large quantities of surface water can, however, be generated due to groundwater emergence at hydrogeological features, such as winterbournes or springs, which activate when the water table is high. Effective flood planning in these catchments is more likely to come in the form of in-channel NFM interventions, hydrogeological mapping, building, and planning restrictions and the development of groundwater emergence early warning systems. We suggest that NFM schemes in these catchments be small-scale and highly targeted to deal with runoff from activated hydrogeological features that would otherwise cause small-scale disruptions (i.e., flooding roads). Catchments classed as impermeable (group 3) are related to the presence of impermeable and slowly permeable soils as well as other less-permeable surfaces, such as urban land use. Due to this, runoff is generated, meaning that this category of catchments is the most suited to NFM in the chalk stream sample.
Overall, this study suggests that implementing NFM in chalk groundwater-dominated catchments is likely to have sub-optimal results compared to other catchment types. However, NFM implementation may be justifiable purely on the merit of the multiple environmental benefits, such as improved water quality, aesthetic improvements, reduced soil erosion, increased biodiversity, and collaboration across multiple river-management bodies. This paper provides a first order triage of the potential for NFM runoff-management methods in chalk catchments. Further work in this field will need to focus on hydrological models that represent the permeability of the soils and the influence of groundwater on stream flows, detailed groundwater emergence, and flood-risk mapping as well as NFM-implementation studies specifically on chalk streams.

Author Contributions

I.B., D.S., J.S. and J.L. conceived and designed the study. I.B. collated the data and undertook the analysis, with statistical modelling input from R.S. I.B. drafted the manuscript, which was then edited by R.S., J.S., J.L., T.S. and D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the UK Environment Agency.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data and the R code script generated for and used for the analysis in this study can be found on Figshare. 10.6084/m9.figshare.14200229.

Acknowledgments

We thank the Environment Agency and the School of Geography and Environmental Science, University of Southampton, for support to Imogen Barnsley. Thanks to Katharine Smith for her editorial assistance.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Catchments in Their Groups after Redundancy Analysis

Group 1—Large Catchments
No.Station NumberRiverLocationCatchment Area (km2)Alternative Group Where Classification Is Uncertain
126001West BeckWansford Bridge195.56
226002HullHempholme Lock397.23
326004Gypsey RaceBridlington257.15
426005Gypsey RaceBoynton248.132
526009West BeckSnakeholme Lock195.61
627087DerwentLow Marishes475.93
733003CamBossington807.32
833004LarkIsleham464.16
933006WisseyNorthwold Total259.36
1033007NarMarham147.39
1133013SapistonRectory Bridge196.16
1233014LarkTemple274.04
1333016CamJesus Lock769.19
1433019ThetMelford Bridge311.38
1533021RheeBurnt Mill308.05
1633022IvelBlunham539.63
1733023Lea BrookBeck Bridge131.893
1833024CamDernford199.59
1933028FlitShefford119.42
2033034Little OuseAbbey Heath707.76
2133044ThetBridgham274.99
2233057OuzelLeighton Buzzrad122.40
2334003BureIngworth161.27
2434004WensumCostessey Mill559.70
2534006WaveneyNeedham Mill376.08
2634011WensumFakenham162.932
2734014WensumSwantom Morely Total377.46
2834019BureHorstead Mill327.91
2936001StourStratford St Mary838.01
3036006StourLangham571.35
3136015StourLamarsh481.283
3238001LeeFeildes Weir1045.10
3338018Upper LeeWater Hall157.733
3438031LeeRye Bridge758.52
3539003WandleSouth Wimbledon153.73
3639010ColneDenham725.90
3739013ColneBerrygrove349.253
3839016KennetTheale1037.87
3939019LambournShaw235.43
4039023WyeHedsor134.182
4139027PangPangbourne175.68
4239030GadeCroxley Green182.35
4339031LambournWelford158.87
4439043KennetKnighton299.172
4539078Wey (North)Farnham192.602
4639103KennetNewbury534.13
4739104MoleEsher471.36
4839115PangBucklebury108.94
4940003MedwayTeston1261.311
5040008Great StourWye226.403
5140011Great StourHorton341.27
5240012DarentHawley187.31
5340016CrayCrayford123.492
5440018DarentLullington122.242
5541004OuseBarcombe Mills400.61
5641009RotherHardham360.77
5742010ItchenHighbridge & Allbrook Total339.91
5842012AntonFullerton186.162
5942016ItchenEaston234.172
6042024TestChilbolton Total478.43
6143003AvonEast Mills Total1459.44
6243004BourneLaverstock165.21
6343005AvonAmesbury326.47
6443006NadderWilton215.63
6543007StourThroop1062.09
6643008WylyeSouth Newton448.17
6743018AllenWalford Mill170.822
6843024WylyeStockton Park170.892
6944001FromeEast Stoke Total414.59
7044002PiddleBaggs Mill183.80
7144004FromeDorchester Total205.67
Group 2—Catchments with More Permeable Soils
No.Station NumberRiverLocationCatchment Area (km2)Alternative Group Where Classification Is Uncertain
126003Foston BeckFoston Mill59.59
226006Elmswell BeckLittle Driffield133.181
326008Mires BeckNorth Cave40.96
426013Driffield TroutDriffield53.34
526016Gypsey RaceKirby Grindalythe16.98
627073Brompton BeckSnainton Ings8.06
729001Waithe BeckBrigsley108.75
829003LudLouth56.423
933025BabinglyWest Newton Mill44.69
1033032HeachamHeacham56.16
1133033HizArlesey112.96
1233040RheeAshwell2.00
1333049Stanford WaterBuckenham Tofts46.45
1433052Swaffham LodeSwaffham Bulbeck33.12
1533054BabingleyCastle Rising48.54
1633056Quy WaterLode92.56
1733061ShepFowlmere One1.03
1833062Guilden BrookFowlmere Two3.40
1933064Whaddon BrookWhaddon14.53
2033065HizHitchin12.03
2133068Cheney WaterGatley End0.11
2234012BurnBurnham Overy83.87
2334018StiffkeyWarham86.133
2438017MimramWhitwell38.40
2539015WhitewaterLodge Farm46.99
2639029Tilling BourneShalford58.78
2739032LambournEast Shefford145.061
2839033Winterbourne StreamBagnor45.31
2939036Law BrookAlbury16.07
3039037KennetMarlborough136.431
3139039WyeHigh Wycombe67.73
3239061Letcombe BrookLetcombe Bassett3.99
3339065Ewelme BrookEwelme11.98
3439077OgMarlborough Poulton Farm63.95
3539091MisbourneQuarrendon Mill65.87
3639101AldbourneRamsbury53.09
3739102MisbourneDenham Lodge93.24
3839107HogsmillEwell8.44
3939112Letcombe BrookArabellas Lake3.10
4039113Manor Farm BrookLetcombe Regis1.38
4139114PangFrilsham90.06
4239118WeyAlton44.50
4339119WeyKings Pond (Alton)46.15
4439120Caker StreamAlton83.94
4539146Mill BrookBlewbury2.01
4639147Wendover SpringsWendover9.49
4740013DarentOtford98.19
4840014WinghamDurlock30.72
4940033DourCrabble Mill44.93
5041015EmsWestbourne57.89
5141023LavantGraylingwell86.20
5241033Costers BrookCocking2.74
5341034EmsWalderton42.46
5441037Winterbourne StreamLewes17.48
5542005Wallop BrookBroughton53.46
5642006MeonMislingford75.84
5742007AlreDrove Lane Alresford57.45
5842008Cheriton StreamSewards Bridge74.33
5942009Candover StreamBorough Bridge72.06
6042015DeverWeston Colley50.15
6142025Lavant StreamLeigh Park55.96
6242026Wallop BrookBossington61.13
6342027DeverBransbury122.361
6443010AllenLoverley Farm94.84
6543011EbbleBodenham105.55
6643012WylyeNorton Bavant114.011
6743014East AvonUpavon85.82
6844006Sydling WaterSydling St Nicholas12.05
6944008South WinterbourneWinterbourne Steepleton20.183
7044009WeyBroadwey8.00
71101003Lukely BrookCarisbrooke Mill14.86
Group 3—Catchments with Less Permeable Soils
No.Station NumberRiverLocationCatchment Area (km2)Alternative Group Where Classification Is Uncertain
126007CatchwaterWithernwick10.84
229002Great EauClaythorpe Mill80.42
333011Little OuseCounty Bridge Euston129.34
433027RheeWimpole128.49
533029StringsideWhitebridge95.41
633030Clipstone BrookClipstone40.35
733045WittleQuidenham27.45
833046ThetRedbridge14.43
933050SnailFordham57.88
1033051CamChesterford140.02
1133053GrantaStapleford113.98
1233055GrantaBabraham101.97
1333066GrantaLinton61.61
1433070LarkFornham St Martin111.07
1534001YareColney228.81
1634002TasShotesham153.19
1734005TudCostessey Park72.11
1836002GlemGlemsford85.61
1936004Chad BrookLong Melford50.33
2036005BrettHadleigh155.85
2136007Belchamp BrookBardfield Bridge58.16
2236008StourWestmill222.82
2336009BrettCockfield25.56
2436010Bumpstead BrookBroad Green27.58
2536011Stour BrookSturmer34.24
2636012StourKedington76.65
2736013BrettHigham191.78
2837012ColnePoolstreet64.49
2937016PantCopford63.80
3038002AshMardock78.00
3138003MimramPanshanger Park130.21
3238004RibWadesmill136.79
3338005AshEasneye84.65
3438006RibHerts Training School148.67
3538011MimramFulling Mill99.25
3638012Stevenage BrookBragbury Brook35.11
3738013Upper LeeLuton Hoo70.31
3838028Stanstead BrookGypsy Lane26.372
3938029QuinGriggs Bridge50.41
4038030BeaneHartham173.85
4139004WandleBeddington Park117.352
4239005Beverly BrookWimbledon Common39.49
4339012HogsmillKingston Upon Thames72.91
4439014VerHansteads134.55
4539028DunHungerford100.09
4639055Yeading Brook WestNorth Hillingdon1.31
4739088ChessRickmansworth96.91
4839089GadeBury Mill44.73
4939125VerRedbourn62.58
5039127MisbourneLittle Missenden47.262
5140015White DrainFairbrook Farm31.40
5240027Sarre PenCalcott19.59
5341003CuckmereSherman Bridge130.33
5441028Chess StreamChess Bridge25.01
5542001WallingtonNorth Fareham111.68
5642011HambleFrogmill55.25

References

  1. Forbes, H.; Ball, K.; McLay, F. Natural Flood Management Handbook; Scottish Environment Protection Agency: Stirling, Scotland, UK, 2015. [Google Scholar] [CrossRef]
  2. Quinn, P.; O’donnell, G.; Nicholson, A.; Wilkinson, M.; Owen, G.; Jonczyk, J.; Barber, N.; Hardwick, M.; Davies, G. Potential Use of Runoff Attenuation Features in Small Rural Catchments for Flood Mitigation; Newcastle University: Newcastle upon Tyne, UK, 2013. [Google Scholar]
  3. SEPA. SEPA Position Statement “The Role of SEPA in Natural Flood Management”; SEPA: Stirling, Scotland, UK, 2012. [Google Scholar]
  4. Blanc, J.; Arthur, S.; Wright, G. Natural Flood Management (NFM) Knowledge System: Part 3—The Effect of Land Drainage on Flood Risk and Farming Practice; CREW—Centre of Expertise for Waters: Aberdeen, Scotland, UK, 2012. [Google Scholar]
  5. Jackson, B.; Wheater, H.; McIntyre, N.; Chell, J.; Francis, O.; Frogbrook, Z.; Marshall, M.; Reynolds, B.; Solloway, I. The Impact of Upland Land Management on Flooding: Insights from a Multiscale Experimental and Modelling Programme. J. Flood Risk Manag. 2008, 1, 71–80. [Google Scholar] [CrossRef]
  6. Schafer, K. Experiments with Catch Crops for Erosions Control in Maize Production. Wirtschaftseigene Futter 1986, 32, 60–71. [Google Scholar]
  7. Schilling, K.E.; Gassman, P.W.; Kling, C.L.; Campbell, T.; Jha, M.K.; Wolter, C.F.; Arnold, J.G. The Potential for Agricultural Land Use Change to Reduce Flood Risk in a Large Watershed. Hydrol. Process. 2014, 28, 3314–3325. [Google Scholar] [CrossRef]
  8. Skinner, K.; Haycock, N. Sinderland Brook—1.8 Km of River and Floodplain Restoration Integrated into a New Housing Development. River Restor. News Newsl. River Restor. Cent. 2005, 20, 6–7. [Google Scholar]
  9. Hankin, B.; Arnott, S.; Whiteman, M.; Burgess-Gamble, L.; Rose, S. Working with Natural Processes—Using Flood Risk Evidence to Make the Case for NFM. 2017. Available online: https://assets.publishing.service.gov.uk/media/6036c7f9d3bf7f0aaa2a45fa/Working_with_natural_processes_using_the_evidence_base.pdf (accessed on 20 October 2017).
  10. Janes, V.J.; Grabowski, R.C.; Mant, J.; Allen, D.; Morse, J.L.; Haynes, H. The Impacts of Natural Flood Management Approaches on In-Channel Sediment Quality. River Res. Appl. 2016, 33, 89–101. [Google Scholar] [CrossRef] [Green Version]
  11. Dadson, S.J.; Hall, J.W.; Murgatroyd, A.; Acreman, M.; Bates, P.; Beven, K.; Heathwaite, A.L.; Holden, J.; Holman, I.; 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]
  12. Nicholson, A.R.; Wilkinson, M.E.; O’Donnell, G.M.; Quinn, P.F. Runoff Attenuation Features: A Sustainable Flood Mitigation Strategy in the Belford Catchment, UK. Area 2012, 44, 463–469. [Google Scholar] [CrossRef]
  13. Lane, S.N. Natural Flood Management. Wiley Interdiscip. Rev. Water 2017, 4, 1211. [Google Scholar] [CrossRef] [Green Version]
  14. Pitt, M. The Pitt Review: Learning Lessons from the 2007 Floods. 2008. Available online: http://webarchive.nationalarchives.gov.uk/20100812084907/http://archive.cabinetoffice.gov.uk/pittreview/_/media/assets/www.cabinetoffice.gov.uk/flooding_review/pitt_review_full pdf.pdf (accessed on 24 July 2018).
  15. O’Neill, R.; Hughes, K. The State of England’s Chalk Streams. 2014. Available online: http://assets.wwf.org.uk/downloads/wwf_chalkstreamreport_final_lr.pdf (accessed on 1 June 2020).
  16. Sear, D.A.; Armitage, P.D.; Dawson, F.H. Groundwater Dominated Rivers. Hydrol. Process. 1999, 13, 255–276. [Google Scholar] [CrossRef]
  17. Berrie, A.D. The Chalk-Stream Environment. Hydrobiology 1992, 248, 3–9. [Google Scholar] [CrossRef]
  18. Wood, P.; Hannah, D.; Agnew, M.; Petts, G. Scales of Hydroecological Variability Within a Groundwater-Dominated Stream. Regul. Rivers Res. Manag. 2001, 17, 347–367. [Google Scholar] [CrossRef]
  19. Lamontagne, S.; Taylor, A.R.; Cook, P.G.; Crosbie, R.S.; Brownbill, R.; Williams, R.M.; Brunner, P. Field Assessment of Surface Water-Groundwater Connectivity in a Semi-Arid River Basin (Murray-Darling, Australia). Hydrol. Process. 2012, 28, 1561–1572. [Google Scholar] [CrossRef] [Green Version]
  20. Naughton, O.; McCormack, T.; Gill, L.; Johnston, P. Groundwater Flood Hazards and Mechanisms in Lowland Karst Terrains. Geol. Soc. London Spec. Publ. 2018, 466, 397–410. [Google Scholar] [CrossRef]
  21. Barber, N.J.; Quinn, P.F. Mitigating Diffuse Water Pollution from Agriculture Using Soft-Engineered Runoff Attenuation Features. Area 2012, 44, 454–462. [Google Scholar] [CrossRef]
  22. Mclean, L.; Beevers, L.; Pender, G.; Haynes, H.; Beevers, L.G.; Pender, H.; Haynes, M.; Wilkinson, J. Variables for measuring multiple benefits and ecosystem services. In Proceedings of the Infrastructure and Environment Scotland 1st Postgraduate Conference, Heriot-Watt University, Edinburgh, UK, 3 June 2013; Available online: https://www.researchgate.net/Publication/271433811 (accessed on 14 May 2019).
  23. The Environment Agency. Working with Natural Processes-Using the Evidence Base to Make the Case for Natural Flood Management SC150005. 2017. Available online: https://assets.publishing.service.gov.uk/media/6036c7f9d3bf7f0aaa2a45fa/Working_with_natural_processes_using_the_evidence_base.pdf (accessed on 16 July 2018).
  24. Boorman, D.B.; Hollist, J.M.; Lilly, A. Hydrology of Soil Types: A Hydrologically-Based Classification of the Soils of the United Kingdom. 1995. Available online: http://nora.nerc.ac.uk/id/eprint/7369/1/IH_126.pdf (accessed on 23 October 2019).
  25. Gustard, A.; Bullock, A.; Dixon, J.M. Low Flow Estimation in the United Kingdom. Report; UK Institute of Hydrology: Wallingford, UK, 1992; Number 108. [Google Scholar]
  26. Fry, M.J.; Swain, O. Hydrological Data Management Systems Within a National River Flow Archive. In Role of Hydrology in Managing Consequences of a Changing Global Environment. Proceedings of the BHS Third International Symposium; Kirby, C., Ed.; British Hydrological Society: London, UK, 2010; pp. 808–815. [Google Scholar]
  27. Smith, A. A New Edition of the Bedrock Geology Map of the United Kingdom. J. Maps 2009, 5, 232–252. [Google Scholar] [CrossRef] [Green Version]
  28. Downing, R.A.; Edmunds, W.M.; Gale, I.N. Regional Groundwater Flow in Sedimentary Basins in the U.K. Geol. Soc. London Spec. Publ. 1987, 34, 105–125. [Google Scholar] [CrossRef]
  29. Baker, D.B.; Richards, R.P.; Loftus, T.T.; Kramer, J.W. A New Flashiness INDEX: Characteristics and Applications to Midwestern Rivers and Streams. JAWRA J. Am. Water Resour. Assoc. 2004, 40, 503–522. [Google Scholar] [CrossRef]
  30. Al-Saady, Y.I.; Al-Suhail, Q.A.; Al-Tawash, B.S.; Othman, A.A. Drainage Network Extraction and Morphometric Analysis Using Remote Sensing and GIS Mapping Techniques (Lesser Zab River Basin, Iraq and Iran). Environ. Earth Sci. 2016, 75, 1243. [Google Scholar] [CrossRef]
  31. Ascott, M.J.; Marchant, B.P.; Macdonald, D.; McKenzie, A.A.; Bloomfield, J. Improved Understanding of Spatio-Temporal Controls on Regional Scale Groundwater Flooding Using Hydrograph Analysis and Impulse Response Functions. Hydrol. Process. 2017, 31, 4586–4599. [Google Scholar] [CrossRef] [Green Version]
  32. NERC; CEH. Land Cover Map 2015: Dataset Documentation. 2017. Available online: www.ceh.ac.Uk (accessed on 1 April 2019).
  33. Ogden, F.L.; Pradhan, N.R.; Downer, C.W.; Zahner, J.A. Relative Importance of Impervious Area, Drainage Density, Width Function, and Subsurface Storm Drainage on Flood Runoff from an Urbanized Catchment. Water Resour. Res. 2011, 47. [Google Scholar] [CrossRef] [Green Version]
  34. Strahler, A.N. Quantitative Geomorphology. Quantitative geomorphology. In Geomorphology. Encyclopedia of Earth Science; Springer: Berlin/Heidelberg, Germany, 1968. [Google Scholar]
  35. Allen, D.J.; Brewerton, L.J.; Coleby, L.M.; Gibbs, B.R.; Lewis, M.A.; Macdonald, A.M.; Wagstaff, S.J.; Williams, A.T. The Physical Properties of Major Aquifers in England and Wales; British Geological Survey Technical Report WD/97/034; British Geological Survey: Keyworth, UK, 1997; p. 333. [Google Scholar]
  36. MacDonald, A.; Allen, D.J. Aquifer Properties of the Chalk of England. Q. J. Eng. Geol. Hydrogeol. 2001, 34, 371–384. [Google Scholar] [CrossRef] [Green Version]
  37. Holko, L.; Parajka, J.; Kostka, Z.; Škoda, P.; Blöschl, G. Flashiness of Mountain Streams in Slovakia and Austria. J. Hydrol. 2011, 405, 392–401. [Google Scholar] [CrossRef]
  38. Sabo, J.L.; Finlay, J.; Kennedy, T.; Post, D. The Role of Discharge Variation in Scaling of Drainage Area and Food Chain Length in Rivers. Science 2010, 330, 965–967. [Google Scholar] [CrossRef] [PubMed]
  39. Sherman, L.K. The Relation of Hydrographs of Runoff to Size and Character of Drainage-Basins. Trans. Am. Geophys. Union 1932, 13, 332–339. [Google Scholar] [CrossRef]
  40. 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. [Google Scholar] [CrossRef] [Green Version]
  41. O’Connell, P.E.; Ewen, J.; O’Donnell, G.; Quinn, P. Is There a Link Between Agricultural Land-Use Management and Flooding? Hydrol. Earth Syst. Sci. 2007, 11, 96–107. [Google Scholar] [CrossRef] [Green Version]
  42. Nirupama, N.; Simonovic, S.P. Increase of Flood Risk Due to Urbanisation: A Canadian Example. Nat. Hazards 2007, 40, 25–41. [Google Scholar] [CrossRef]
  43. Alaoui, A.; Caduff, U.; Gerke, H.H.; Weingartner, R. Preferential Flow Effects on Infiltration and Runoff in Grassland and Forest Soils. Vadose Zone J. 2011, 10, 367–377. [Google Scholar] [CrossRef]
  44. Granged, A.J.; Zavala, L.M.M.; Jordán, A.; Bárcenas-Moreno, G. Post-Fire Evolution of Soil Properties and Vegetation Cover in a Mediterranean Heathland After Experimental Burning: A 3-Year Study. Geoderma 2011, 164, 85–94. [Google Scholar] [CrossRef]
  45. Miller, J.; Kim, H.; Kjeldsen, T.R.; Packman, J.; Grebby, S.; Dearden, R. Assessing the Impact of Urbanization on Storm Runoff in a Peri-Urban Catchment Using Historical Change in Impervious Cover. J. Hydrol. 2014, 515, 59–70. [Google Scholar] [CrossRef] [Green Version]
  46. Wong, L.; Hashim, R.; Ali, F. A Review on Hydraulic Conductivity and Compressibility of Peat. J. Appl. Sci. 2009, 9, 3207–3218. [Google Scholar] [CrossRef]
  47. Wollenberg, A.L.V.D. Redundancy Analysis an Alternative for Canonical Correlation Analysis. Psychometrika 1977, 42, 207–219. [Google Scholar] [CrossRef]
  48. R Core Team. A Language and Environment for Statistical Computing. 2019. Available online: https://www.r-project.Org/ (accessed on 6 January 2021).
  49. Box, G.E.P.; Cox, D.R. An Analysis of Transformations. J. R. Stat. Soc. Ser. B Statistical Methodol. 1964, 26, 211–243. [Google Scholar] [CrossRef]
  50. Oksanen, J.; Blanchet, G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. 2019. Available online: https://www.mcglinnlab.org/publication/2019-01-01_oksanen_vegan_2019/ (accessed on 18 April 2020).
  51. Dray, S.; Legendre, P.; Blanchet, G. Packfor: Forward Selection with Permutation. 2016. (accessed on 7 January 2021).
  52. Legendre, P.; Legendre, L. Numerical Ecology; Legendre, P., Legendre, L., Eds.; Elsevier: Amsterdam, The Netherlands, 2012. [Google Scholar]
  53. Lawson, R.G.; Jurs, P.C. New Index for Clustering Tendency and Its Application to Chemical Problems. J. Chem. Inf. Comput. Sci. 1990, 30, 36–41. [Google Scholar] [CrossRef]
  54. Saraçli, S.; Doğan, N.; Doğan, I. Comparison of Hierarchical Cluster Analysis Methods by Cophenetic Correlation. J. Inequalities Appl. 2013, 2013, 203. [Google Scholar] [CrossRef]
  55. Borcard, D.; Gillet, F.; Legendre, P. Numerical Ecology with R. In Numerical Ecology with R; Springer Science and Business Media LLC: Berlin, Germany, 2011. [Google Scholar] [CrossRef]
  56. Gower, J.C. A General Coefficient of Similarity and Some of Its Properties. Biometrics 1971, 27, 857. [Google Scholar] [CrossRef]
  57. Rousseeuw, P.J. Silhouettes: A Graphical Aid to the Interpretation and Validation of Cluster Analysis. J. Comput. Appl. Math. 1987, 20, 53–65. [Google Scholar] [CrossRef] [Green Version]
  58. Gerisch, M. Habitat Disturbance and Hydrological Parameters Determine the Body Size and Reproduction Strategy of Alluvial Ground Beetles. ZooKeys 2011, 100, 353–370. [Google Scholar] [CrossRef] [Green Version]
  59. Mazouz, R.; Assani, A.A.; Rodríguez, M.A. Application of Redundancy Analysis to Hydroclimatology: A Case Study of Spring Heavy Floods in Southern Québec (Canada). J. Hydrol. 2013, 496, 187–194. [Google Scholar] [CrossRef]
  60. Pham, S.V.; Leavitt, P.R.; McGowan, S.; Neto, P.P. Spatial Variability of Climate and Land-Use Effects on Lakes of the Northern Great Plains. Limnol. Oceanogr. 2008, 53, 728–742. [Google Scholar] [CrossRef] [Green Version]
  61. Wilkinson, M.E.; Addy, S.; Quinn, P.F.; Stutter, M. Natural Flood Management: Small-Scale Progress and Larger-Scale Challenges. Scott. Geogr. J. 2019, 135, 23–32. [Google Scholar] [CrossRef]
  62. Wingfield, T.; Macdonald, N.; Peters, K.; Spees, J.; Potter, K. Natural Flood Management: Beyond the Evidence Debate. Area 2019, 51, 743–751. [Google Scholar] [CrossRef]
  63. Woodward, M.; Gouldby, B.P.; Kapelan, Z.; Khu, S.; Townend, I.; Park, H. The Use of Real Options in Optimum Flood Risk Management Decision Making. 2009. Available online: https://eprints.hrwallingford.Com/767/ (accessed on 31 October 2020).
  64. Metcalfe, P.; Beven, K.; Hankin, B.; Lamb, R. Simplified Representation of Runoff Attenuation Features Within Analysis of the Hydrological Performance of a Natural Flood Management Scheme. Hydrol. Earth Syst. Sci. 2018, 22, 2589–2605. [Google Scholar] [CrossRef] [Green Version]
  65. Pattison, I.; Lane, S.N.; Hardy, R.J.; Reaney, S.M. The Role of Tributary Relative Timing and Sequencing in Controlling Large Floods. Water Resour. Res. 2014, 50, 5444–5458. [Google Scholar] [CrossRef] [Green Version]
  66. Ferguson, C.R.; Fenner, R.A. The Potential for Natural Flood Management to Maintain Free Discharge at Urban Drainage Outfalls. J. Flood Risk Manag. 2020, 13, e12617. [Google Scholar] [CrossRef]
  67. Dixon, S.J.; Sear, D.A.; Odoni, N.A.; Sykes, T.; Lane, S.N. The Effects of River Restoration on Catchment Scale Flood Risk and Flood Hydrology. Earth Surf. Process. Landf. 2016, 41, 997–1008. [Google Scholar] [CrossRef]
  68. Hankin, B.; Lamb, R.; Craigen, I.; Page, T.; Metcalfe, P.; Chappell, N. A Whole Catchment Approach to Improve Flood Resilience in the Eden: Winning Entry to the Defra Floods Competition 2016; JBA Consulting: Skipton, UK, 2017. [Google Scholar]
  69. Bell, F.G. Engineering Properties of Soils and Rocks, 3rd ed.; Butterworths: London, UK, 1985. [Google Scholar]
  70. Misstear, B.D.R.; Brown, L.; Daly, D. A Methodology for Making Initial Estimates of Groundwater Recharge from Groundwater Vulnerability Mapping. Hydrogeol. J. 2008, 17, 275–285. [Google Scholar] [CrossRef]
  71. Caraher, R. An Examination of De-Culverting as a Mechanism of Flood Management; Dundalk Institute of Technology: Dundalk, Ireland, 2015. [Google Scholar]
  72. Acreman, M.C.; Riddington, R.; Booker, D.J. Hydrological Impacts of Floodplain Restoration: A Case Study of the River Cherwell, UK. Hydrol. Earth Syst. Sci. 2003, 7, 75–85. [Google Scholar] [CrossRef] [Green Version]
  73. DEFRA. Making Space for Water Developing a New Government Strategy for Flood and Coastal Erosion Risk Management in England a Consultation Exercise. 2004. Available online: www.defra.gov.Uk (accessed on 17 May 2019).
  74. Naughton, O.; Johnston, P.; McCormack, T.; Gill, L. Groundwater Flood Risk Mapping and Management: Examples from a Lowland Karst Catchment in Ireland. J. Flood Risk Manag. 2015, 10, 53–64. [Google Scholar] [CrossRef] [Green Version]
  75. Smith, P.A.; Dosser, J.; Tero, C.; Kite, N. A method to identify chalk rivers and assess their nature-conservation value. Water Environ. J. 2003, 17, 140–144. [Google Scholar] [CrossRef]
  76. Atkinson, T.C.; Smith, D.I. Rapid Groundwater Flow in Fissures in the Chalk: An Example from South Hampshire. Q. J. Eng. Geol. Hydrogeol. 1974, 7, 197–205. [Google Scholar] [CrossRef]
  77. Macdonald, D.; Dixon, A.; Newell, A.; Hallaways, A. Groundwater Flooding Within an Urbanised Flood Plain. J. Flood Risk Manag. 2011, 5, 68–80. [Google Scholar] [CrossRef] [Green Version]
  78. Marsh, T.J.; Dale, M. The UK Floods of 2000? 2001: A Hydrometeorological Appraisal. Water Environ. J. 2002, 16, 180–188. [Google Scholar] [CrossRef]
  79. Pinault, J.-L.; Amraoui, N.; Golaz, C. Groundwater-Induced Flooding in Macropore-Dominated Hydrological System in the Context of Climate Changes. Water Resour. Res. 2005, 41. [Google Scholar] [CrossRef] [Green Version]
  80. Arnott, S. An Assessment of Methods for Catchment-Scale Identification of Goundwater-Suface Water Interractions in Lowland Permeable Catchments; University of Exeter: Exeter, UK, 2008. [Google Scholar]
  81. Bradford, R.B. Controls on the Discharge of Chalk Streams of the Berkshire Downs, UK. Sci. Total Environ. 2002, 282–283, 65–80. [Google Scholar] [CrossRef]
  82. Molenat, J.; Gascuel-Odoux, C.; Ruiz, L.; Gruau, G. Role of Water Table Dynamics on Stream Nitrate Export and Concentration in Agricultural Headwater Catchment (France). J. Hydrol. 2008, 348, 363–378. [Google Scholar] [CrossRef]
  83. Griffiths, J.A.; Binley, A.; Crook, N.; Nutter, J.; Young, A.; Fletcher, S. Streamflow Generation in the Pang and Lambourn Catchments, Berkshire, UK. J. Hydrol. 2006, 330, 71–83. [Google Scholar] [CrossRef]
  84. Holman, I.P.; Hollis, J.M.; Bramley, M.E.; Thompson, T.R.E. The Contribution of Soil Structural Degradation to Catchment Flooding: A Preliminary Investigation of the 2000 Floods in England and Wales. Hydrol. Earth Syst. Sci. 2003, 7, 755–766. [Google Scholar] [CrossRef] [Green Version]
  85. Wilkinson, M.; Quinn, P.; Benson, I.; Welton, P.; Kerr, P. Runoff Management: Mitigation Measures for Disconnecting Flow Pathways in the Belford Burn Catchment to Reduce Flood Risk; British Hydrological Society: London, UK, 2010. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The processes of flooding in impermeable upland catchments (a) and in permeable chalk catchments (b) along with common NFM interventions to demonstrate how groundwater dominance in flood production may impact the efficacy of NFM.
Figure 1. The processes of flooding in impermeable upland catchments (a) and in permeable chalk catchments (b) along with common NFM interventions to demonstrate how groundwater dominance in flood production may impact the efficacy of NFM.
Water 13 02212 g001
Figure 2. Location of the 198 catchments included in this study underlain by chalk bedrock. Catchment boundaries are from the National River Flow Archive, and the bedrock map is from the British Geological Survey.
Figure 2. Location of the 198 catchments included in this study underlain by chalk bedrock. Catchment boundaries are from the National River Flow Archive, and the bedrock map is from the British Geological Survey.
Water 13 02212 g002
Figure 3. (a) The optimum number of average linkage clusters according to the Mantel statistic. The bars show the correlations between original distance matrix and binary matrices computed from the dendrogram cut at various levels. (b) Silhouette plot of the final three groups from the average linkage agglomerative clustering.
Figure 3. (a) The optimum number of average linkage clusters according to the Mantel statistic. The bars show the correlations between original distance matrix and binary matrices computed from the dendrogram cut at various levels. (b) Silhouette plot of the final three groups from the average linkage agglomerative clustering.
Water 13 02212 g003
Figure 4. (a) Biplot showing catchment locations along two RDA axes. Dark orange vectors represent hydrological variables, and blue vectors are the catchment properties, with points for catchments placed relative to the vectors based on their individual values for catchment properties. (b) The catchment classifications mapped spatially. Misclassified catchments are presented in their reclassified grouping and outlined in bold. (c) The uncertainty of catchment classification according to silhouette widths. Misclassified catchments (<0) were reclassified into their nearest neighbour alternative classification.
Figure 4. (a) Biplot showing catchment locations along two RDA axes. Dark orange vectors represent hydrological variables, and blue vectors are the catchment properties, with points for catchments placed relative to the vectors based on their individual values for catchment properties. (b) The catchment classifications mapped spatially. Misclassified catchments are presented in their reclassified grouping and outlined in bold. (c) The uncertainty of catchment classification according to silhouette widths. Misclassified catchments (<0) were reclassified into their nearest neighbour alternative classification.
Water 13 02212 g004
Table 1. Equations and methodology of variables compiled for the response and explanatory databases for the redundancy analysis.
Table 1. Equations and methodology of variables compiled for the response and explanatory databases for the redundancy analysis.
Variable Description and Influence on Flow
Mean Discharge (QMean) Q n
n = days in record
The average quantity of water in the river channel. Gives an indication of the discharge under normal conditions from the NRFA [26].
Maximum Recorded Discharge (QMax) Q Max Maximum discharge capacity of catchments at the gauging site during high flows from the NRFA [26].
Richard–Baker Flashiness Index (RBFI) RBFI = i = 1 n | q i q i 1 | i 1 n q i
q = dimensionless measure of discharge
i = time
n = number of discharge measurements
Measures the absolute daily fluctuations in streamflow, divided by the sum of all stream flow for the time series length [37]. Values range between 0 and 1. Values near 0 represent stable flow and those close to 1 represent highly changeable flows and spate conditions [29].
Base Flow Index BFI measures the proportion of river runoff derived from stored sources [25,26].
Catchment Area (Ca)-Watershed boundary shapefiles from the NRFA [26] facilitate calculations for other catchment properties. Smaller catchments are associated with steeper hydrograph rising limbs due to reduced catchment complexity [38,39].
Form Factor (Rf) Rf = C a B l 2 Measures the geometric shape of the catchment. Circular catchments (Rf = 0) have steep hydrograph peaks [34] (p. 304).
Drainage Density (Dd) Dd = S t L C a High drainage densities are linked to drainage efficiency and high peak flows [33].
Channel Slope m/mCalculated from a raster layer of SRTM 30m Digital Elevation Model [40] in ArcGIS using the zonal statistics function.
% Land Cover Type L a n d u s e a r e a C a × 100 Land-use effects on hydrology [32].
Arable Land and HorticultureReduced peak flows [41].
Broadleaf WoodlandReduced peak flows [42].
Coniferous WoodlandReduced peak flows due to large leaf surface area [42].
GrasslandNo overall influence [43].
HeathlandNo overall influence—the effects of vegetation are counteracted by shallow soils and low storage capacity [44].
UrbanIncreased peak flows—impervious surfaces and drainage systems [45].
SuburbanIncreased peak flows—impervious surfaces and drainage systems [45].
Marshland
Inland RockIncreased peak flows—imperviousness.
% Soil Cover Type H O S T   S o i l T y p e C a × 100 The Hydrology Of Soil Types (HOST) from CEH [24].
Bypass Flow CommonPermeable
(Vertical saturated hydraulic conductivity > 10 cm/day−1)
Bypass Flow UncommonPermeable
(Vertical saturated hydraulic conductivity > 10 cm/day−1)
Bypass Flow Very UncommonPermeable
(Vertical saturated hydraulic conductivity > 10 cm/day−1)
Bypass Flow VariableSemi-permeable
(Vertical saturated hydraulic conductivity 0.1–10 cm/day−1)
Slowly Permeable SoilsSemi-Permeable (Vertical saturated hydraulic conductivity 0.1–10 cm/day−1)
ImpermeableImpermeable (<0.1 cm/day−1 Vertical saturated hydraulic conductivity)
PeatSemi-permeable
(Vertical saturated hydraulic conductivity 0.1–10 cm/day−1) [46].
Transmissivity (m2/day) Transmissivity = g a 3 12 v
g = acceleration due to gravity
a = area
v = kinematic viscosity of the fluid
The rate at which water passes through the chalk bedrock [35,36].
Abstraction ScoreSubjectively assigned an arbitrary abstraction score based on the flow regime description on the gauging station info in the NRFA. Scores: flow added; natural flow; minor, moderate and major reduction of flows due to abstraction.Quantifies influence of abstraction on river regime [26].
Table 2. Explanatory variables chosen during the forward selection process along with their comparative explanatory power (adjusted R2) and significance (F statistic and p Values). The explanatory variables represent the most parsimonious model and explain the most variation in the hydrological response dataset.
Table 2. Explanatory variables chosen during the forward selection process along with their comparative explanatory power (adjusted R2) and significance (F statistic and p Values). The explanatory variables represent the most parsimonious model and explain the most variation in the hydrological response dataset.
VariableOrderR2Cumulative R2Adjusted Cumulative R2F Statisticp Value
Area20.3570.3570.353108.60.001
Uncommon Bypass Flow180.2310.5880.584109.30.001
Impermeable Soils220.0340.6220.58417.50.001
Slowly Permeable Soils210.0370.6590.65220.80.001
Station Elevation10.020.6690.6605.80.001
Form Factor40.0100.6790.6696.10.001
Urban140.0070.6860.6754.430.001
Table 3. VIF values of the explanatory variables kept in the parsimonious RDA model. VIF values exceeding 10 indicate collinearity between variables.
Table 3. VIF values of the explanatory variables kept in the parsimonious RDA model. VIF values exceeding 10 indicate collinearity between variables.
AreaUncommon Bypass FlowImpermeable SoilsSlowly Permeable Soils
VIF1.1434.6251.6914.482
Station ElevationForm FactorUrban
VIF1.2481.0291.057
Table 4. Gower distance and cophenetic correlation results for all four hierarchical clustering methods for comparison. A high cophenetic correlation indicates a strong relationship between the distance matrix and the cophenetic matrix describing the dendrogram split where each point is clustered differently. A lower Gower distance indicates an acceptable clustering algorithm.
Table 4. Gower distance and cophenetic correlation results for all four hierarchical clustering methods for comparison. A high cophenetic correlation indicates a strong relationship between the distance matrix and the cophenetic matrix describing the dendrogram split where each point is clustered differently. A lower Gower distance indicates an acceptable clustering algorithm.
Single LinkageComplete LinkageAverage LinkageWard’s Clustering
Gower distance30,878.36062,837.2406538.116
Cophenetic correlation0.7570.7480.802
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barnsley, I.; Spake, R.; Sheffield, J.; Leyland, J.; Sykes, T.; Sear, D. Exploring the Capability of Natural Flood Management Approaches in Groundwater-Dominated Chalk Streams. Water 2021, 13, 2212. https://doi.org/10.3390/w13162212

AMA Style

Barnsley I, Spake R, Sheffield J, Leyland J, Sykes T, Sear D. Exploring the Capability of Natural Flood Management Approaches in Groundwater-Dominated Chalk Streams. Water. 2021; 13(16):2212. https://doi.org/10.3390/w13162212

Chicago/Turabian Style

Barnsley, Imogen, Rebecca Spake, Justin Sheffield, Julian Leyland, Tim Sykes, and David Sear. 2021. "Exploring the Capability of Natural Flood Management Approaches in Groundwater-Dominated Chalk Streams" Water 13, no. 16: 2212. https://doi.org/10.3390/w13162212

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