Next Article in Journal
Salt Tolerance and Desalination Abilities of Nine Common Green Microalgae Isolates
Next Article in Special Issue
Effects of Residue Cover on Infiltration Process of the Black Soil Under Rainfall Simulations
Previous Article in Journal
A Prediction–Correction Solver for Real-Time Simulation of Free-Surface Flows in River Networks
Previous Article in Special Issue
Evolution Pattern of Tailings Flow from Dam Failure and the Buffering Effect of Debris Blocking Dams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Impact of Man–Made Ponds on Soil Erosion and Sediment Transport in Limnological Basins

1
National Council for Scientific Research, Remote Sensing Center, Beirut 11-8281, Lebanon
2
Centre de Recherches en Sciences et Ingénierie, Lebanese University Faculty of Engineering II, Roumieh 1205, Lebanon
3
Centre d’Études et de Développement des Territoires et de l’Environnement, Université d’Orléans, 45100 Orléans, France
4
Faculty of Science and Technology, Free University of Bozen, 39100 Bozen, Italy
*
Author to whom correspondence should be addressed.
Water 2019, 11(12), 2526; https://doi.org/10.3390/w11122526
Submission received: 17 October 2019 / Revised: 26 November 2019 / Accepted: 27 November 2019 / Published: 29 November 2019
(This article belongs to the Special Issue Modeling of Soil Erosion and Sediment Transport)

Abstract

:
The impact of ponds on basins has recently started to receive its well-deserved scientific attention. In this study, pond-induced impacts on soil erosion and sediment transport were investigated at the scale of the French Claise basin. In order to determine erosion and sediment transport patterns of the Claise, the Coordination of Information on the Environment (CORINE) erosion and Soil and Water Assessment Tool (SWAT) models were used. The impact of ponds on the studied processes was revealed by means of land cover change scenarios, using ponded versus pondless inputs. Results show that under current conditions (pond presence), 12.48% of the basin corresponds to no-erosion risk zones (attributed to the dense pond network), while 65.66% corresponds to low-erosion risk, 21.68% to moderate-erosion risk, and only 0.18% to high-erosion risk zones. The SWAT model revealed that ponded sub-basins correspond to low sediment yields areas, in contrast to the pondless sub-basins, which yield appreciably higher erosion rates. Under the alternative pondless scenario, erosion risks shifted to 1.12%, 0.52%, 76.8%, and 21.56% for no, low, moderate, and high-erosion risks, respectively, while the sediment transport pattern completely shifted to higher sediment yield zones. This approach solidifies ponds as powerful human-induced modifications to hydro/sedimentary processes.

1. Introduction

Understanding the impact of land occupation (land use/cover) on basin processes, such as rainfall-runoff and soil erosion, is an integral part of land and water management-oriented decisions [1,2]. The processes of soil erosion and sediment transport take part in van Rijn’s (1993) [3] sedimentary cycle, and often are the main causes of soil loss in basins [4]. Although these processes are of natural origin [5], the interaction between climate, soil, topography, land use, and land cover significantly influences erosion rates and sediment loads [6,7]. Soil loss due to these processes is a frequent problem that hydrologists, land planners, and basin managers will need to contend with [8]. Accordingly, soil erosion quantification in erosion-prone areas, with the highest accuracy possible, provides a complete knowledge of soil loss hotspots and allows prioritized treatment measures for supporting land processes in the concerned basin [9]. Subsequently, the reliable assessment and representation of sediment yields—which depend on the cascading effect of soil erosion—allows an in-depth understanding of the soil erosion-sediment link at the basin scale [9].
For representing both processes, several models have been developed and used extensively to replace the conventional assessment methods, i.e., the Water Erosion Prediction Project (WEPP) [10,11], the Universal Soil Loss Equation (USLE) [12], the Modified Universal Soil Loss Equation (MUSLE) [13], the Revised Universal Soil Loss Equation (RUSLE) [14,15] and the Soil and Water Assessment Tool (SWAT) [16]. SWAT is one of the most widely used basin models. It has been applied extensively in modeling the impact of land occupation changes, under different scenarios and different contexts [17,18]. The widespread use of SWAT can be justified by its sensitivity and flexibility towards the land occupation input [19], its adaptability to different contexts—even to those with data scarcity [20]—its simple data requirements and ease of computation [21], as well as the straightforward calibration through its stand-alone SWATCUP interface [22].
Despite abundant research regarding the impact of land occupation on soil loss, few studies focus on the particular case of small water bodies and their effect as a land occupation class [23]. Small water bodies, like ponds and wetlands, are considered as the most amplified form of human-induced modifications to the hydro-sedimentological system of basins [24]. Ponds represent a total of over 90% of global standing water bodies, 30% of global standing waters by surface area [25], and form the most widespread aquatic habitat dominating the continental standing waters in Europe [26]. Despite their well-documented significance [27], abundant numbers, and increasing proliferation [28], ponds have not received considerable scientific attention compared to rivers and lakes [26]. It is worth mentioning that research regarding ponds in Europe has tripled in the last decade [29], where results showed that ponds contribute significantly to several basin related processes [30]. Examples of these processes are sediment interception [31], removal of pollutants for river protection [32], nutrient recycling [33], greenhouse gas emission [34], regulation of hydrological flows [35], biogeochemistry [30], and climate [36]. In addition to their environmental role, ponds have a well-known value for housing and sustaining biodiversity, supporting livelihoods, local economies, and taking part in the socio-cultural heritage of the settings in which they are located [37].
Under the hydro-sedimentological scope, particularly, ponds have shown to retain as much as 90% of sediments transported in basins [38]. Consequently, ponds have been heavily blamed for rupturing the ecological and sedimentary continuum of the basins to which they belong [39]. The disruptive effect of ponds is due to the increase of residence time of waters, resulting in a decline in the temporal variation of the main discharge [40]. Accordingly, the deceleration of overland flow allows suspended particles to settle under the effect of their weight, causing a reduction in the amount of sediments entrained by water, making ponds sediment sinks [40]. However, this effect strongly depends on their position in the basin, their depth, volume, slope [41], as well as the surrounding land occupation [42]. Winfield Fairchild and Velinsky (2006) [43] showed that ponds located upstream of rivers—considering their sediment retention capacity—are capable of creating a state of imbalance in the geochemical and hydro-sedimentary status of the underlying rivers. Consequently, the Directive Cadre sur l’Eau (DCE) [44] stresses the need to assess the impact of hydromorphological elements that are capable of influencing hydrologic pathways, river morphology, width, and continuity.
Beyond the contribution of isolated ponds, connected networks of ponds were found to contribute to basin processes at higher rates than lakes or even rivers [45,46]. Particularly in France, ponds are mainly concentrated in three regions: the Sologne region, Brenne (Central France), and Dombes (Eastern France). In response to DCE recommendations, this study aims to assess the impact of man–made ponds on soil erosion and sediment transport, at the scale of the Indre portion of the Claise basin. This part corresponds to the Brenne Natural Regional Park that houses 4500 waterbodies (ponds, marshes, and small water surfaces), 2179 of which are located in Claise, being part of an interconnected network. To evaluate erosion risks in the Claise, the Coordination of Information on the Environment CORINE (1992) Erosion Risk model will be used, since it presents a simplification of the reliable USLE model [1] and given that no-erosion field data, for the Claise, were available for the study. SWAT is employed to assess the impact of ponds on the Claise’s hydro-sedimentary regime. This choice is owed to SWAT’s ability to simulate the physical processes that occur in ponds, which in turn allows an accurate representation of pond containing basins [38]. This is achieved through the SWAT Pond (.pnd) input file that makes SWAT one of the few hydrological models having an input for ponds [47]. Related studies have assessed the behavior of ponds using SWAT [38,48] and highlighted the efficiency of SWAT on this part.
The impact of the Claise’s ponds on erosion and sediment transport will be assessed by testing alternative scenarios, where the land occupation input for both models will be simulated with and without ponds. By this approach, a quantification of the pond impact can be obtained. The presented work serves as a decision-oriented tool for basins similar to the Claise, where pond proliferation has been halted until a proper understanding of their effect is established. In addition, analysis of soil erosion risks and sediment transport is useful for conservation measures that aim towards prolonging the useful life of these small water bodies or for ceasing their proliferation.

2. Materials and Methods

2.1. Study Area

The Brenne portion of Claise basin (Figure 1) covers an area of 707 km2 and is one of the three entities housing the 4500 ponds of the Brenne Natural Regional Park, which are mostly grouped in the form of interconnected chains [36]. Within the Claise, 2179 of the Brenne waterbodies, along with the Five Bonds channel (Blizon), are the main contributors [49] to the 87.6 km long Claise River having an average discharge of 4.5 m3/s [50]. The basin is mostly dominated by a degraded oceanic climate with a mean temperature of 11 °C and average annual precipitation of 700 mm [51]. As to its topographical profile, the basin is considered to be a flat plain, with 99% of its surface falling into the 0%–5% slope class, while its altitude varies between 76 m and 181 m. The six poorly-permeable soil classes of the basin are mostly dominated by Luvisols [36,52]. At a combined state, the Claise’s challenging pedology, flat topographical setting, and quasi–impermeable lithology have resulted in the stagnation of incoming water leading to the formation of natural ponds [53]. According to Bennarrous (2009) [37], however, these ponds are not only a product of natural processes but also an anthropogenic adaptation to a poorly drained domain and a source of economic livelihood (aquaculture) in an environment of limited productivity. As a result of intensive pond proliferation throughout time, the basin has acquired a particular hydrographic network characterized by an abundance of different kinds of water bodies. Despite its richness, nonetheless, the hydrographic network of the basin is randomly organized and presents severe fragmentation [54]. In contrast to the evolving pond proliferation, the land occupation setting of the Claise has been relatively unchanged for the last 19 years, mainly displaying a dominance of an interlocked mosaic of forests, grasslands, and agricultural areas [54].

2.2. The SWAT Model: Basins and Impoundments

SWAT is one of the few hydrological models that have an input for ponds [47], which makes it an ideal tool for this study. SWAT is a semi-physical deterministic distributed and continuous hydrologic model that functions at a daily time step with options for sub-hourly routing, as well [55]. It is a quite complex model with significant input data demands. However, the basic components can be readily obtained and are relatively simple [56]. Essential inputs consist of weather data, digital elevation model (DEM), soil, and land occupation maps [57]. Based on these inputs, SWAT divides the basin into sub-basins, which further divides into smaller hydrological response units (HRUs) [58]. HRUs are defined as units with homogeneous land occupation, topography, and pedological properties [59]. These are generated in order to lump somewhat similar areas scattered through the basin into a single unit, simplifying the model’s run by avoiding unpractical simulations while accounting for the diversity of different factors in the basin [60]. Within SWAT, ponds are defined as waterbodies, being an integral part of a sub-basin’s hydrological network; they are capable of intercepting surface runoff [60], thereby modifying the hydro-sedimentary behavior of the basin. Since this work targets the erosion/sedimentary behavior of the SWAT model, only their related equations will be presented. Further details can be found in Neitsch et al. (2011) [47].
Soil erosion is calculated based on the Modified Universal Soil Loss Equation (MUSLE) using the following formula [55]:
sed = 11.8 × (Qsurf · qpeak · areahru)0.56 × KUSLE × CUSLE × PUSLE × LSUSLE × CFRG
where sed is the HRU sediment yield (t); Qsurf is the surface runoff volume (mm); qpeak is the peak runoff rate (m3/s); areahru is the area of the HRU (ha); KUSLE is the Universal Soil Loss Equation (USLE) soil erodibility factor; CUSLE is the USLE cover and management factor; PUSLE is the USLE support practice factor; LSUSLE is the USLE topographic factor; and CFRG is the coarse fragment factor.
It should be noted that MUSLE virtually calculates the sediment yields due to soil erosion that reach the main streams of the sub-basins in a unit of time. This is the reason why sediment delivery ratios are not required, contrarily to USLE and RUSLE [61].
The mass balance equation of sediment in ponds is described by the following formula [47]:
sedwb = sedwbi + sedflowinsedstlsedflowout
where sedwb is the amount of sediment at the end of the day; sedwbi is the amount of sediment at the day’s beginning; sedflowin is the amount of sediment provided from inflows; sedstl is the amount of settled sediments; and sedflowout is the amount of sediment transported as outflow. All components are expressed in metric tons.

2.3. The CORINE Erosion Model

The CORINE erosion model is a simplification of the USLE model. In the CORINE model, erosion risks are classified on a scale of 0–3, with 0 corresponding to the no-erosion class, 1 to low erosion risks, 2 to moderate erosion risks, and 3 to high-erosion risks [62]. For the estimation of erosion using the CORINE model, parameters such as soil erodibility, climate erosivity, topography (slope), and LU/LC (vegetation cover) are required [61]. Each parameter, in turn, consists of several sub-factors, and is classified according to the CORINE model into respective indices (Figure 2). Once established, the soil erodibility index is combined with climate erosivity and slope, to yield the potential soil erosion risk map. Subsequently, the potential soil erosion risk map indices (0–3) are crossed with those of the vegetation cover to yield the actual soil erosion map.

2.4. Input Data and Adaptation to the SWAT and CORINE Models

The input data used in this study, along with a short description, are summarized in Table 1.
ArcSWAT version 2012.10_3.19 was applied in an ArcGIS 10.3 (Environmental Systems Research Institute, Inc. (Esri), 380 New York Street, Redlands, California, USA) environment to perform simulations at a daily time step for the period 1970–2018. The first seven years of the simulation, from 1970 to 1976, were used for the model’s warm up. The calibration and validation phases of the model were carried out by means of comparison between simulated discharges and measured discharge data from the station L6202040–La Claise au Grand–Pressigny (Pont de Fer), obtained from Eau France Banque HYDRO–Ministère de l’Ecologie, de l’Energie, du Développement Durable et de l’Aménagement du Territoire (MEEDDAT)/Direction Générale de la Prévention des Risques (DGPR)/Service des Risques Naturels et Hydrauliques (SRNH). Due to the unavailability of sediment records for this study, only a hydrologic calibration was performed. However, surface runoff and stream flow—regardless of the approach used to model the sedimentary cycle—are the driving forces of runoff and streambed erosion, as well as of overland and stream sediment transport [65]. Hence, one might say that the calibration of the sediment part of the study might be implicit, but nonetheless robust. According to Hallouz et al. (2018) [66], a calibrated hydrologic model gives a certain degree of reliability to the sedimentary output. However, the availability of sediment data is often a constraint in relative studies, as measurements, either in the form of gross or net erosion or in the form of total sediment discharge, in streams, are often nonexistent. Even in cases where such data is available, it is in the form of sparse discrete measurements and rarely in the form of continuous sediment graphs, where they could be used for a robust calibration. In the case of the Claise, this would be a very challenging task, given the large number of ponds. For this reason, SWAT was calibrated according to the calibration scheme of Jalowska and Yuan (2019) [38], where a complete representation of an impoundment effect (here, ponds) is ensured. Nonetheless, the absence of sediment records poses some short of limitation that needs to be considered by the decision-makers of the basin.

2.4.1. DEM: Topographic Effect and Use in Both Models

In the CORINE model, DEM is used to extract the topographic parameter by means of slope computation. The slope angles of the study area were obtained using the “slope” tool of ArcGIS. Subsequently, the slope raster was classified, with respect to the CORINE’s model indices, into (1) very gentle to flat (<5%), (2) gentle (5%–15%), (3) steep (15%–30%), and (4) very steep (>30%).
In SWAT, the DEM is used for the extraction of the topographic parameters and for representing the basin’s physical parameters such as slopes, flow direction and accumulation, delineation of the hydrologic network, and basin partitioning [67]. After inputting the DEM into SWAT, the basin was initially divided into 25 sub-basins using SWAT’s watershed delineator. However, given the presence of a great number of waterbodies, and their impact on output accuracy, the number of sub-basins was increased to 35. The purpose of increasing the sub-basins’ number was to increase the spatial resolution for a more detailed representation of the basins’ processes [68]. Since the SWAT model allows only one pond per sub-basin [39] during the sub-basins’ delineation, the largest number of ponds was included in each in order to ensure maximal representation of pond-induced processes. Further, all ponds within each sub-basin were lumped and the outlets of the sub-basins were chosen to coincide with those of the ponds in order to maximize their representation and account for their effect. The SWAT delineated hydrologic network was, in turn, verified against a pre-defined stream network in order to ensure and improve its accuracy.

2.4.2. Pedology: Adaptation to the Different Requirements of the Models Used

The pedological composition of the Claise was determined from the Harmonized World Soil Database (HWSD) [52]. For the CORINE model, soil texture was determined using the USDA textural triangle after inputting the respective percentages of sand, silt, and clay for each soil group obtained from the HWSD. Texture was then classified with respect to the CORINE indices. The parameters of depth and stoniness were treated and classified similarly. After obtaining the texture, depth, and stoniness sub-factors, the three layers were input into the “raster calculator” in analogy to Figure 2.
For the SWAT model, the hydrologic soil group of each class was assigned following the United States Department of Agriculture [69] Soil Conservation Service (SCS) soil survey. Bulk densities were computed using the Soil Water Characteristics software following Saxton and Rawl (2006) [70] equations that were shown to have adequate accuracy for bulk density computation by Al Sayah et al. (2019) [71]. Likewise, the available water capacity and saturated hydraulic conductivities were also computed using the same software. Organic carbon content was determined by multiplying the organic matter content derived from the HWSD by 0.58, since organic carbon forms around 58% of the soil’s organic matter as a rule of thumb [72]. Sand, silt, and clay percentages, as well as the rock fragment contents were extracted from the HWSD. The USLE_K factor was computed following the formula presented in the SWAT documentation [55], while soil surface albedo was determined using the following formula [73]:
Soil albedo (0.3 − 2.8 μm) = 0.069·(color value) − 0.114
After building the soil database, all parameters were reclassified in SWAT using user-adapted look-up tables for integration into the SWAT database.

2.4.3. Weather Data: Forcing on-Field Weather Data to a Meteorological Model

Weather data was obtained from two sources, on-field weather station data and Météo France’s Système d’Analyse Fournissant des Renseignements Adaptés à la Nivologie (SAFRAN) model [74]. The SAFRAN model was used since it contains records for the period 1970–2018 while the weather stations, located next to a pond network, were used to account for the climatic regulating effect reported by Nedjai et al. (2018) [36]. SAFRAN was used in order to provide a large time span for the study. Moreover, the SAFRAN was validated by test-correlation with weather station data using Pearson’s correlation coefficient “r” [75]. Prior to correlation, harmonization of both datasets was performed since the weather stations record parameters at the hourly time step, while SAFRAN simulates at the daily time step. Results of the correlation are presented in Table 2.
For the CORINE model, the rainfall and temperature parameters are used for computation of the MFI and BGI indices following their respective formulas (Figure 2); these indices were then used to calculate erosivity following the workflow of Figure 2. In SWAT, the weather database was built using the WGN maker macro-tool [76] and input into the SWAT database.

2.4.4. Land Use and Land Cover: A Particularly Rich Natural Limnological Setting

The land occupation map of the Claise was obtained by on-screen digitizing of aerial photography at 0.5 m resolution. Such fine scale was used in order to ensure a detailed representation of the study area, particularly to account as accurately as possible for the scale of ponds.
For the CORINE model, the obtained land occupation map was reclassified into fully protected (forest, permanent pasture, and scrublands) and not-fully protected (cultivated or bare land) areas as demonstrated in Table 3.
In the case of SWAT, land use and land cover classes were reclassified into SWAT’s classes. Particular attention was given to the water (in SWAT terms WATR) classes. As reported by Jalowska and Yuan (2019) [38], the reason for this is that though SWAT allows the creation of HRUs with WATR, water bodies should be modeled either as reservoirs or ponds. Almendinger et al. (2014) and Jalowska and Yuan (2019) [38,48] have shown that an accurate representation of basin processes requires the integration of the SWAT model’s impoundments function. Furthermore, Wang et al. (2008) [77] highlighted the importance of considering impoundments, such as ponds, by testing scenarios of impoundment integration versus impoundment disregarding in a basin covered by only 3% of impoundments. Their results confirmed that simulations were considerably affected even with this small cover. In the same manner, Jalowska and Yuan (2019) [38] simulated different scenarios following integration or absence of impoundments like reservoirs rather than just a normal “water” land occupation class. They noted that disregarding impoundments leads to a series of uncertainties starting by an inaccurate SWAT performance, which in turn leads to inefficient calibration efforts, and overall inaccurate model performance.
After computing the slope, erosivity and soil erodibility indices, the potential soil erosion risk map was constructed. By combining the potential soil erosion risk map to the vegetation cover layer in the “raster calculator” tool and crossing each ones indices, the actual soil erosion risk map was obtained.

3. Results

3.1. CORINE Erosion Model Outputs for the Claise Basin

In the following sections, the erosion assessment for the Claise is presented and each component is detailed accordingly.

3.1.1. Soil Erodibility

Soil erodibility parameters are presented in Figure 3. The dominant textural class in the basin was found to be loam (68.4%) followed by loamy sand (28.75%), while the remainder percentages are clay (2.48%) and sand (0.34%). In terms of texture, since most of the study area is covered by loam with respect to [62], the study area predominantly falls into the highly erodible class. In terms of soil depth, 94.6% of the Claise fits to the slightly erodible class (1000 mm > 750 mm), while the remainder 5.4% rests within the moderately erodible class (300 mm, corresponding to the 250–750 mm class). As far as stoniness is concerned, only 6% of the Claise is under the fully protected cover.
Table 4 provides the numerical description of the soil’s erodibility components.
By inputting the texture, depth and stoniness parameters into the “raster calculator” tool of ArcGIS, the soil erodibility raster was generated and then reclassified using the “reclassify” tool (Figure 4).
As seen in Figure 4, the generated soil erodibility map revealed that about 96.83% (684.37 km2) of the study area is covered by highly erodible soils, while the remainder 3.17% (22.63 km2) correspond to moderate-risk areas. This fact is mainly due to the textural distribution of soil and the little stone cover of the Claise basin.

3.1.2. Erosivity in a Degraded Oceanic Climate Setting

Meteorological data of the study area is presented in Table 5. As seen, no dry month exists in the Claise and highest temperatures are recorded during the month of July making the basin subject to continuous precipitation.
MFI was found to be 80, corresponding to the low variability class. This signifies evenly distributed rainfall, thus reducing risks of climate-induced soil erosion [78]. The Bagnouls–Gaussen Aridity Index (BGI) on the other hand, was found to be “0”, signifying that the study area corresponds to the humid regime as a result of its oceanic influence. In terms of CORINE erosivity, the Modified Fournier Index (MFI) belongs to class 2 variability (low), while BGI corresponds to class 1 (humid). Accordingly, the erosivity index of the Claise watershed was found to be 2, corresponding to the low erosivity class.

3.1.3. Topography (Slope): A Reduced Effect in a Flat Setting

Using the “slope” tool in ArcGIS, slope angles were extracted from the DEM. Figure 5 displays the slope of the Claise basin and the adapted reclassified raster following CORINE classification.
From Figure 5, it can be observed that the Claise presents a dominantly flat topography with a maximum slope of 15°. In CORINE terms, 99.3% (702.051 km2) of the study area corresponds to “very gentle to flat” slope class (<5°) and 0.7% (4.94 km2) to the “gentle” slope category (5°–15°). The particular setting of low relief and small slope significantly plays a role in the reduction of erosion generated by runoff [79], despite the high soil erodibility risks.

3.1.4. Potential Soil Erosion Map

By overlaying soil erodibility (Figure 4), topography and erosivity, using the “raster calculator” tool, the potential soil erosion risk map was obtained (Figure 6), following the formula presented in Figure 2.
From Figure 6, low-erosion risks present 3% (21.22 km2), moderate risks present 96.5% (682.25 km2) and high risks form 0.5% (3.53 km2). These percentages, however, reflect only the potential soil erosion risks which according to CORINE (1992) [62] do not take into consideration the vegetative cover at this stage. The slope effect and equal repartition of precipitation, as reflected by the MFI values and the absence of arid periods (low BGI), are potentially responsible for moderating the soil erodibility map yielding the dominantly moderate potential soil erosion risks.

3.1.5. Vegetation Cover: The Presence/Absence Effect of Ponds

This factor is the focal point of this study since by slightly changing this parameter the outcome changes significantly. At this point, two vegetative cover scenarios are presented: the first presents the actual setting of the Claise accounting for the presence of ponds, while the second simulates a scenario where ponds are removed to assess the difference in erosion outcomes with and without their presence. This last step allowed to quantify the impact of ponds on erosion in the Claise basin. In the second scenario, the land occupation group “ponds” was changed to their surrounding class (grasslands).
The large area of ponds, which makes up around 11% of the Claise (Table 3), displays their potential role as modifiers basin processes. Figure 7 presents the two considered vegetation covers as inputs for crossing with the potential soil erosion risk map to yield two different actual soil erosion risk maps. These maps were then compared to evaluate the effect of the presence and absence of ponds and their role on erosion.
Class 0 refers to land cover categories that are not considered in CORINE; these categories are urban areas and water bodies, while classes 1 and 2 refer to the fully protected and not fully protected covers. From Figure 7b, 72.5% of the Claise basin corresponds to the fully protected class, while 27.5% of the study area is occupied by not fully protected cover. To assess the impact of pond presence, the second scenario of replacing ponds by their surrounding dominant cover was performed.

3.1.6. Actual Soil Erosion Maps Under Current and Alternative Scenarios

The two actual soil erosion risks maps were produced by multiplying the respective indices of the potential soil erosion risk map and the two vegetation cover scenarios using the “raster calculator” tool. Figure 8 reveals the outcome under both scenarios.
The Actual soil erosion map was cross-checked against: Institut National de la Recherche Agronomique (INRA) 2000 erosion maps and the combined GIS sol–INRA–SOeS 2011 maps [80]. These were produced following the Modèle d’Évaluation Spatiale de l’Aléa d’Érosion des Sols (MESALES). The comparison between the INRA maps, and the produced actual soil erosion risk map, is presented in Table 6. By this comparison, it is concluded that there is a good agreement between these maps. In addition, since INRA maps are not completely adequate to be considered at the basin scale [80], the established erosion maps are considered for the no-erosion zones, overcoming, this way, the challenge of coarse representation.
As can be seen from Table 6, a total of 134 validation points were chosen. These were divided into 86, 37, and 11 low-erosion, moderate-erosion, and high-erosion zones validation points, respectively. A large part of the moderate-erosion class was misinterpreted as the low-erosion class. This discrepancy is due to the fact that the no-erosion zones do not exist in the INRA maps, but are instead classified as low-erosion zones. Therefore, the error margin in the moderate-erosion class from Table 6 is justified by the finer scale representation of the produced maps, compared to the INRA maps. The overall accuracy was determined to be 75%, while the computed Cohen’s kappa coefficient [81] was found to be 0.7; this indicates a substantial agreement between the INRA maps, and the produced actual soil erosion risk map. The kappa coefficient was used since it tests inter-rater reliability; i.e., the coefficient represents the extent to which the generated data are correct representations of the measured data. In the case of this study, the generated data is the actual soil erosion risk map, while the measured data consists of the validation points obtained from the INRA maps.
From Figure 8, three main results can be drawn:
(1) The role of vegetation cover in changing erosion risks is solidified. This is particularly reflected by the setting of the Claise basin due to the agricultural and grass cover. By comparing the actual soil erosion risk map with the potential soil erosion risk map and statistics, a shift of erosion risk classes is observed. As mentioned in Section 3.1.4, considering the potential soil erosion risks, and ignoring the vegetation cover, low, moderate and high-risk areas take over 3%, 96.5% and 0.5% of the total basin area, accordingly. When the vegetation cover layer was taken into account the resulting actual soil erosion risk shifted to 65.66%, 21.68% and 0.18%, for low, moderate and high-risk areas. These observations solidify that vegetative cover is the most influential aspect for erosion assessment. In further detail regarding the vegetation cover layer, areas corresponding to agricultural classes are seen to have higher erosion risks than areas with different land cover types. This is in agreement with Verheijen et al. (2009) [82] observations that despite the considerable effect of soil type, topography and climatic conditions, the major influencer of soil erosion is the vegetative cover, especially cultivated areas.
(2) The remainder 12.48% of the actual soil erosion risk map is the no-erosion zone. As seen in Figure 7a, most of the no-erosion zone corresponds to the concentration area of ponds, while the remainder 1.48% represents the Claise River. The ponded area represents 88.23 km2 of the Claise under no risk of erosion, making these ponds a counter-erosion zone.
(3) At a graphical scale, a complete shift from low to moderate risks, in the greatest part of the basin, is observed. Table 7 presents the statistical difference between the actual soil erosion risks with current vegetation cover and those of the pondless scenario. From Table 7, the effective role of ponds as an erosion counter-measure is revealed.
Additionally, the impact of ponds on erosion at the scale of the basin is revealed. Not only did the no-erosion and low-erosion classes decrease by 11.36% and 65.14%, in the absence of ponds scenario, but also the moderate and high-erosion risks increased by 55.12% and 21.38%, respectively. These changes are due to several reasons:
(1) The most evident reason is that ponds effectively and directly nullify splash erosion in the areas they occupy.
(2) Their widespread, yet dense, positioning throughout the basin, counteracts runoff erosion in a twofold way: first, by intercepting eroded soils by overland flow, retaining this way, the transported material and preventing them from reaching the streams, and second, by slowing surface runoff and thus, abating its erosive force. Despite the fact that the low-slope topography does not particularly favor runoff erosion, let alone high velocities of overland flow, this obviously has some effect, especially in cases of intense rainfall events.
(3) Their dense aggregation in the basin attributes them the role of cascade check dams, containing sediments.
(4) Their large density and chain sequence where the retention effect is greatly amplified (factor of 2179 ponds) [38].
(5) The highly erodible setting of the basin resulting from a challenging pedology.

3.2. Sediment Transport in a Limnologically Rich Setting

Sediment transport in the Claise was simulated using the SWAT model. The model was calibrated using the SUFI-2 algorithm of SWATCUP following Jalowska and Yuan (2019) [38] proposed sequential calibration for settings characterized by the presence of water impoundments. As mentioned previously, due to the unavailability of sediment measurements for this study, only a hydrologic calibration was performed with the most related sensitive parameters (Table 8). Results of the hydrologic calibration yielded an R2 of 0.7 during calibration and 0.67 during validation.
As a result of an exhaustive hydrologic calibration, the average sediment yield of the Claise basin for the studied period under current and alternative conditions is presented in Figure 9.
As can be seen from Figure 9, low sediment yields correspond to ponded sub–basins in contrast to pondless sub-basins which present lower soil loss rates. However, this gradient varies according to the setting of the ponds in these sub–basins. Despite that sub-basins 1 and 9 contain a small number of ponds, these correspond to areas of high sediment yield as a result of their dominant agricultural cover and their location in the basin’s steepest areas. This shows that the presence of ponds alone is not sufficient to alter the sediment yields. In fact, the setting of ponds under the form of a collective dense network is the main modifier of sediment transport patterns. In order to further highlight the impact of ponds, alternative scenario testing was performed where a pondless land occupation setting was re-inputted to the calibrated SWAT model in order to determine subsequent sediment transport changes. Results of sediment yields, highlighting the impact of ponds, are presented in Table 9.
In the low sediment yield sub-basins, where the vast majority of ponds lies, occupying a large portion of the sub-basins they are located in, the land cover change from pond to grasslands is far too great and the effect of ponds presence/absence is quite obvious. Accordingly, the role of ponds as integral features of the landscape processes is solidified. Having determined their effect, an increased exposure for their integration into management plans is recommended, while an understanding of basin scale sediment dynamics offers insights regarding the role of ponds in water resources. Such implications could matter not only for the Claise and similar basins, but also for expectations regarding the landscape role of sediment retention basins, and for the general understanding of sediment transport in rain-dominated contexts. This, in turn, could be considered as a contribution to hydrologic modelling efforts, where small waterbodies, especially artificial ones, are often neglected, or even smoothed out, from digital elevation models.

4. Conclusions

An investigation of the pond-induced effects on soil erosion and sediment transport of limnologically rich basins was presented. Through this task, recommendations of the European framework for the Thematic Strategy on Soil Protection were addressed by revealing the different levels of soil loss, represented by providing an insight to the investigation of erosion-prone regions and sediment yield zones of different levels. Furthermore, recommendations of the DCE regarding the behavioral understanding of hydromorphological alternating factors (ponds as hydro-sedimentary elements), at the basin scale, were also considered. Despite the Claises’ weakly structured pedology, resulting in high erodibility, the Claise was found to have low erosion risks and sediment transport rates due to several reasons, like the evenly distributed rainfall, the relatively flat topography, and most importantly due to the presence of its dense pond network that acts as a natural measure against soil erosion. This was solidified by simulating a scenario where ponds were substituted by grasslands. After replacing ponds by a protective vegetative cover, all erosion risks and sediment yield classes of the Claise significantly varied: no and low-erosion risk zones decreased, while moderate and high-erosion risk zones increased. Regarding sediment transport, the replacement of ponds by grasslands led to a complete disappearance of low sediment yield zones and considerable increases of moderate and high sediment yield zones. Accordingly, the “safe” soil loss status of the Claise can be attributed to the low soil loss rates of the Claise basin to the presence of these ponds. Despite their protection against soil erosion, however, their presence in very large numbers might cause a distortion in the sediment balance of the underlying rivers. Such cases may lead to sediment starvation and force the river to engage in increasing streambed and bank erosion [83]. Therefore, the use of the presented approach may serve as an efficient tool towards the orientation of future decisions regarding the proliferation or cease of ponds depending on their effect.

Author Contributions

Conceptualization: R.N. and M.J.A.S.; data curation: M.J.A.S.; methodology: M.J.A.S. and R.N.; project administration: R.N., C.A., and M.K.; supervision: R.N., K.K., C.A., and M.K.; validation: M.J.A.S.; visualization: M.J.A.S. and K.K.; writing—original draft: M.J.A.S. and K.K.; writing—review and editing: M.J.A.S., K.K., R.N., C.A., and M.K.

Funding

This research is part of a PhD thesis funded by the National Council of Scientific Research—Lebanon (CNRS—L), Agence Universitaire de la Francophonie (AUF), Lebanon, and the Lebanese University. It is also part of the Dynétangs project funded by the French Centre–Val–de–Loire region.

Acknowledgments

The authors would like to express their gratitude for the funding agencies, the editor, and reviewers for leveraging the quality of this work and to the Brenne Natural Park for their help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gyamfi, C.; Ndambuki, J.M.; Salim, R.W. Simulation of sediment yield in a semi–arid river basin under changing land use: An integrated approach of hydrologic modelling and principal component analysis. Sustainability 2016, 8, 1133. [Google Scholar] [CrossRef]
  2. Dwarakish, G.S.; Ganasri, B.P. Impact of land use change on hydrological systems: A review of current modeling approaches. Cogent Geosci. 2015, 1, 1115691. [Google Scholar] [CrossRef]
  3. van Rijn, L.C. Principles of Sediment Transport in Rivers, Estuaries and Coastal Seas; Aqua Publications: Amsterdam, The Netherlands, 1993; p. 690. [Google Scholar]
  4. Esa, E.; Assen, M.; Legass, A. Implications of land use/cover dynamics on soil erosion potential of agricultural watershed, northwestern highlands of Ethiopia. Env. Syst. Res. 2018, 7, 21. [Google Scholar] [CrossRef]
  5. Aksoy, H.; Mahe, G.; Meddi, M. Modeling and Practice of Erosion and Sediment Transport under Change. Water 2019, 11, 1665. [Google Scholar] [CrossRef]
  6. Zhang, S.; Li, Z.; Lin, X.; Zhang, C. Assessment of Climate Change and Associated Vegetation Cover Change on Watershed–Scale Runoff and Sediment Yield. Water 2019, 11, 1373. [Google Scholar] [CrossRef]
  7. Bussi, G.; Dadson, S.J.; Prudhomme, C.; Whitehead, P.G. Modelling the future impacts of climate and land–use change on suspended sediment transport in the River Thames (UK). J. Hydrol. 2016, 542, 357–372. [Google Scholar] [CrossRef]
  8. Ffolliott, P.F.; Brooks, K.C.; Neary, D.G.; Tapia, P.R.; Chevesich, P.G. Soil Erosion and Sediment Production on Watershed Landscapes: Processes and Control, UNESCO Special Technical Publication No. 32; USDA: Washington, DC, USA, 2013; p. 73.
  9. Aga, A.O.; Chane, B.; Melesse, A.M. Soil erosion modelling and risk assessment in data Scarce Rift Valley Lake Regions, Ethiopia. Water 2018, 10, 1684. [Google Scholar] [CrossRef]
  10. Laflen, J.M.; Lane, L.J.; Foster, G.R. WEPP: A new generation of erosion prediction technology. J. Soil Water Conserv. 1991, 46, 34–38. [Google Scholar]
  11. Flanagan, D.C.; Gilley, J.E.; Franti, T.G. Water Erosion Prediction Project (WEPP) Development History, Model Capabilities and Future Enhancements. Am. Soc. Agric. Biol. Eng. 2007, 50, 1603–1612. [Google Scholar] [CrossRef]
  12. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses from Cropland East of the Rocky Mountains–Guide for Selection of Practices for Soil and Water Conservation, Agriculture Handbook No. 282; USDA: Washington, DC, USA, 1965.
  13. Williams, J.R. Sediment–Yield Prediction with Universal Soil Loss Equation Using Runoff Energy Factor. Present and Prospective Technology for Predicting Sediment Yields and Sources. In Proceedings of the Sediment Yield Workshop, Oxford, MS, USA, 17–18 January 1975; pp. 244–252. [Google Scholar]
  14. Renard, K.G.; Foster, G.R.; Weesies, G.A.; Mc Cool, D.K.; Yoder, D.C. Predicting Soil Erosion by Water: A Guide to Conservation Planning with the Revised Universal Soil Loss Equation (RUSLE); US Government Printing Office: Washington, DC, USA, 1997; Volume 703.
  15. Renard, K.G.; Yoder, D.C.; Lightle, D.T.; Dabney, S.M. . Universal Soil Loss Equation and Revised Universal Soil Loss Equation. In Handbook of Erosion Modelling; Morgan, R.P.C., Nearing, M., Eds.; Blackwell Publishing: Hoboken, NJ, USA, 2011; pp. 137–166. [Google Scholar]
  16. Arnold, J.G.; Srinivasan, R.; Muttiah, R.S.; Williams, J.R. Large area hydrologic modeling and assessment part I: Model development. J. Am. Water. Resour. Assoc. 1998, 34, 73–89. [Google Scholar] [CrossRef]
  17. Samaras, A.G.; Koutitas, C.G. Modeling the impact of climate change on sediment transport and morphology in coupled watershed–coast systems: A case study using an integrated approach. Int. J. Sediment. Res. 2015, 29, 304–315. [Google Scholar] [CrossRef]
  18. Kaffas, K.; Righetti, M.; Avesani, D.; Spiliotis, M.; Hrissanthou, V. Coupling CFSv2 with ArcSWAT for seasonal hydrological forecasting in a Mediterranean basin. In Proceedings of the 11th World Congress on Water Resources and Environment (EWRA 2019), Madrid, Spain, 25–29 June 2019; pp. 459–460. [Google Scholar]
  19. Napoli, M.; Massetti, L.; Orlandini, S. Hydrological response to land use and climate changes in a rural hilly basin in Italy. Catena 2017, 157, 1–11. [Google Scholar] [CrossRef]
  20. Alemaw, D.; Ayana, E.K.; Legesse, E.S.; Moges, M.M.; Tilahun, S.A.; Moges, M.A. Estimating reservoir sedimentation using bathymetric differencing and hydrologic modeling in data scarce Koga watershed, Upper Blue Nile, Ethiopia. J. Agric. Environ. Int. Dev. 2016, 110, 413–427. [Google Scholar] [CrossRef]
  21. Dhami, B.S.; Pandey, A. Comparative Review of Recently Developed Hydrologic Models. J. Indian Water Resour. Soc. 2013, 33, 34–42. [Google Scholar]
  22. Abbaspour, K.C.; Vaghefi, S.A.; Srinivasan, R. A guideline for successful calibration and uncertainty analysis for soil and water assessment: A review of papers from the 2016 international SWAT conference. Water 2017, 10, 6. [Google Scholar] [CrossRef]
  23. Biggs, J.; von Fumetti, S.; Kelly–Quinn, M. The importance of small waterbodies for biodiversity and ecosystem services: Implications for policy makers. Hydrobiologia 2016, 793, 3–39. [Google Scholar] [CrossRef]
  24. Mujere, N.; Eslamian, S. Climate Change Impacts on Hydrology and Water Resources. In Handbook of Engineering Hydrology Modeling, Climate Change and Variability; Eslamian, S., Ed.; CRC–Press: Boca Raton, FL, USA, 2014; pp. 114–125. [Google Scholar]
  25. Downing, J.A.; Prairie, Y.T.; Cole, J.J.; Duarte, C.M.; Tranvik, L.J.; Striegl, R.G.; McDowell, W.H.; Kortelainen, P.; Caraco, N.F.; Melack, J.M.; et al. The global abundance and size distribution of lakes, ponds, and impoundments. Limnol. Ocean. 2006, 51, 2388–2397. [Google Scholar] [CrossRef]
  26. Miracle, M.R.; Oertli, B.; Céréghino, R.; Hull, A. Preface: Conservation of European ponds–current knowledge and future needs. Limnetica 2010, 29, 1–8. [Google Scholar]
  27. Céreghino, R.; Biggs, J.; Oertli, B.; Declerck, S. The Ecology of European Ponds: Defining the Characteristics of a Neglected Freshwater Habitat, Hydrobiologia; Springer: Berlin/Heidelberg, Germany, 2008; Volume 597, pp. 1–6. [Google Scholar] [CrossRef]
  28. Williams, P.; Whitfield, M.; Biggs, J. How Can We Make New Ponds Biodiverse? A Case Study Monitored Over 7 Years, Hydrobiologia; Springer: Berlin/Heidelberg, Germany, 2008; Volume 597, pp. 137–148. [Google Scholar] [CrossRef]
  29. Céréghino, R.; Boix, D.; Cauchie, H.M.; Martens, K.; Oertli, B. The Ecological Role of Ponds In a Changing World, Hydrobiologia; Springer: Berlin/Heidelberg, Germany, 2014; Volume 723, pp. 1–6. [Google Scholar] [CrossRef]
  30. Sharpley, A.N.; Bergström, L.; Aronsson, H.; Bechmann, M.; Bolster, C.H.; Börling, K.; Djodjic, F.; Jarvie, H.P.; Schoumans, O.F.; Stamm, C.; et al. Future agriculture with minimized phosphorus losses to waters: Research needs and direction. Ambio 2015, 44, 163–179. [Google Scholar] [CrossRef]
  31. Oertli, B.; Frossard, P.A. Mares et Etangs: Ecologie, Conservation, Gestion, Valorisation, 1st ed.; Presses Polytechniques et Universitaires Romandes: Lausanne, Switzerland, 2013; p. 512. [Google Scholar]
  32. Hill, M.J.; Hassall, C.; Oertli, B.; Fahrig, L.; Robson, B.J.; Biggs, J.; Samways, M.J.; Usio, N.; Takamura, N.; Krishnaswamy, J.; et al. New policy directions for global pond conservation. Conserv. Lett. 2018, 11, 8. [Google Scholar] [CrossRef]
  33. Indermuehle, N.; Oertli, B.; Biggs, J.; Céréghino, R.; Grillas, P.; Hull, A.; Nicolet, P.; Scher, O. Pond conservation in Europe: The European Pond Conservation Network (EPCN). Int. Vereinigung Theor. Angew. Limnol. Verh. 2008, 30, 446–448. [Google Scholar] [CrossRef]
  34. Xiao, S.; Yang, H.; Liu, D.; Zhang, C.; Lei, D.; Wang, Y.; Peng, F.; Li, Y.; Wang, C.; Li, X.; et al. Gas transfer velocities of methane and carbon dioxide in a subtropical shallow pond. Tellus B Chem. Phys. Meteorol. 2014, 66, 23795. [Google Scholar] [CrossRef]
  35. Mujere, N.; Eslamian, S. Impact of Urbanization on Runoff Regime. In Handbook of Engineering Hydrology Modeling, Climate Change and Variability; Eslamian, S., Ed.; CRC–Press: Boca Raton, FL, USA, 2014; pp. 114–125. [Google Scholar]
  36. Nedjai, R.; Azaroual, A.; Chlif, K.; Bensaid, A.; Al Sayah, M.; Ysbaa, L. Impact of Ponds on Local Climate: A Remote Sensing and GIS Contribution Application to the Ponds of Brenne (France). J. Earth. Sci. Clim. Chang. 2018, 9, 10. [Google Scholar]
  37. Bennarrous, R. La Grande Brenne aux Périodes Préindustrielles (Indre) Contribution à L’histoire des Paysages, des Etangs et des Relations Sociétés/Milieux Dans une Zone Humide Continentale. Approches Historique, Archéologique et Paléo–Environnementale; Université de Paris I–Panthéon Sorbonne: Paris, France, 2009. [Google Scholar]
  38. Jalowska, A.M.; Yuan, Y. Evaluation of SWAT Impoundment Modeling Methods in Water and Sediment Simulations. J. Am. Water. Resour. Assoc. 2019, 55, 209–227. [Google Scholar] [CrossRef]
  39. Clifford, C.C.; Heffernan, J.B. Artificial Aquatic Ecosystems. Water 2018, 10, 1096. [Google Scholar] [CrossRef]
  40. Ebel, J.D.; Lowe, W.H. Constructed Ponds and Small Stream Habitats: Hypothesized Interactions and Methods to Minimize Impacts. J. Water. Resour. Prot. 2013, 5, 723–731. [Google Scholar] [CrossRef]
  41. Ismail, T.; Othman, M.A.; Fadzil, A.B.; Zainuddin, Z.M. Deposition of Sediments in Detention Pond. Malaysian. J. Civ. Eng. 2010, 22, 95–118. [Google Scholar]
  42. Declerck, S.; De Bie, T.; Ercken, D.; Hampel, H.; Schrijvers, S.; Van Wichelen, J.; Gillard, V.; Mandiki, R.; Losson, B.; Bauwens, D.; et al. Ecological characteristics of small farmland ponds: Associations with land use practices at multiple spatial scales. Biol. Conserv. 2006, 131, 523–532. [Google Scholar] [CrossRef]
  43. Fairchild, G.W.; Velinsky, D.J. Effects of Small Ponds on Stream Water Chemistry. Lake. Reserv. Manag. 2006, 22, 321–330. [Google Scholar] [CrossRef]
  44. Commission Européenne. La Directive–Cadre Européenne sur l’Eau; Commission Européenne: Brussels, Belgium, 2000. [Google Scholar]
  45. Biggs, J.; Williams, P.; Whitfield, P.; Nicolet, P.; Weatherby, A. 15 years of pond assessment in Britain: Results and lessons learned from the work of Pond Conservation. Aquat. Conserv. Mar. Freshw. Ecosyst. 2005, 15, 693–714. [Google Scholar] [CrossRef]
  46. Jeffries, M. Local–Scale Turnover of Pond Insects: Intrapond Habitat Quality and Inter–Pond Geometry are Both Important, Hydrobiologia; Springer: Berlin/Heidelberg, Germany, 2005; Volume 543, pp. 207–220. [Google Scholar]
  47. Neitsch, S.L.; Arnold, J.G.; Kiniry, J.R.; Williams, J.R. Soil & Water Assessment Tool Theoretical Documentation Version 2009. Texas Water Resources Institute Technical Report No.406: College Station, TX, USA, 2011. [Google Scholar]
  48. Almendinger, J.E.; Murphy, M.S.; Ulrich, J.S. Use of the Soil and Water Assessment Tool to Scale Sediment Delivery from Field to Watershed in an Agricultural Landscape with Topographic Depressions. J. Environ. Qual. 2014, 43, 9–17. [Google Scholar] [CrossRef] [PubMed]
  49. Guichané, R. L’aménagement hydraulique de la Claise tourangelle et de ses affluents du Moyen–Âge à nos jours/Mills on the Claise and its tributaries in Indre–et–Loire from the Middle Ages to modern times. Rev. Archeol. Centre France 1993, 32, 109–152. [Google Scholar] [CrossRef]
  50. SANDRE. Fiche Cours D’eau la Claise (L6––0200); SANDRE: Paris, France, 2012. [Google Scholar]
  51. Joly, D.; Brossard, T.; Cardot, H.; Cavailhes, J.; Hilal, M.; Wavresky, P. Les types de climats en France, une construction spatiale–Types of climates on continental France, a spatial construction. Cybergéo Eur. J. Geogr. 2010, 501, 1–23. [Google Scholar]
  52. Fischer, G.; Nachtergaele, F.; Prieler, S.; van Velthuizen, H.T.; Verelst, L.; Wiberg, D. Global Agro–Ecological Zones Assessment for Agriculture (GAEZ 2008); IIASA: Laxenburg, Austria; FAO: Rome, Italy, 2008. [Google Scholar]
  53. Dauphin, P.; Mansons, J.; Pellé, B.; Airault, V.; Trotignon, J.; Boyer, P.; Chatton, T.; Issa, N. Document D’objectifs des Sites Natura 2000 FR2410003 “Brenne” et FR2400534 “Grande Brenne.”. Available online: http://natura2000.mnhn.fr/uploads/doc/PRODBIOTOP/1847_DOCOB_BRENNE.pdf (accessed on 29 November 2019).
  54. Bouscasse, H.; Defrance, P.; Amand, B.; Grandmougi, B.; Strosser, P.; Beley, Y. Amélioration des connaissances sur les fonctions et usages des zones humides: Évaluation économique sur des sites tests le cas des étangs de la Grande Brenne. Available online: http://oaidoc.eau-loire-bretagne.fr/exl-doc/DOC00021327_s3.pdf (accessed on 29 November 2019).
  55. Arnold, J.G.; Kiniry, J.R.; Srinivasan, R.; Williams, J.R.; Haney, E.B.; Neitsch, S.L. Soil & Water Assessment Tool Input/Output Documentation Version 2012; Texas Water Resources Institute: College Station, TX, USA, 2012; p. 541. [Google Scholar]
  56. Kaffas, K.; Hrissanthou, V.; Sevastas, S. Modeling hydromorphological processes in a mountainous basin using a composite mathematical model and ArcSWAT. Catena 2018, 162, 108–129. [Google Scholar] [CrossRef]
  57. Gassman, P.W.; Reyes, M.R.; Green, C.H.; Arnold, J.G. The Soil and Water Assessment Tool: Historical Development, Applications, and Future Research Directions. Am. Soc. Agric. Biol. Eng. 2007, 50, 1211–1250. [Google Scholar] [CrossRef] [Green Version]
  58. Epelde, A.M.; Cerro, I.; Sánchez–Pérez, J.M.; Sauvage, S.; Srinivasan, R.; Antigüedad, I. Application of the SWAT model to assess the impact of changes in agricultural management practices on water quality. Hydrol. Sci. J. 2014, 60, 825–843. [Google Scholar] [CrossRef]
  59. Srinivasan, R.; Zhang, X.; Arnold, J. SWAT ungauged: Hydrological budget and crop yield predictions in the Upper Mississippi river basin. Trans. ASABE 2010, 53, 1533–1546. [Google Scholar] [CrossRef]
  60. Neitsch, S.L.; Arnold, J.G.; Kiniry, J.R.; Rinivasan, R.S.; Williams, J.R. Soil and Water Assessment Tool User’s Manual–Version 2000; Texas Water Resources Institute: College Station, TX, USA, 2002. [Google Scholar]
  61. Kaffas, K.; Hrissanthou, V. Annual sediment yield prediction by means of three soil erosion models at the basin scale. In Proceedings of the 10th World Congress of EWRA “Panta Rhei”, Athens, Greece, 5–9 July 2017; pp. 1346–1352. [Google Scholar]
  62. CORINE. CORINE Soil Erosion Risk and Important Land Resources-in the Southern Regions of the European Community; European Community-European Environment Agency: Copenhagen, Denmark, 1992; p. 541. [Google Scholar]
  63. Fournier, F. Climat et Erosion: La Relation Entre l’Erosion du sol par l’Eau et Les Precipitations Atmospheriques; Presses Universitaires De France: Paris, France, 1960; p. 201. [Google Scholar]
  64. Gaussen, H. Bioclimatic Map of Mediterranean Zone; UNESCO and FAO: Paris, France, 1963. [Google Scholar]
  65. Arnold, J.G.; Moriasi, D.N.; Gassman, P.W.; Abbaspour, K.C.; White, M.J.; Srinivasan, R.; Santhi, C.; Harmel, R.D.; van Griensven, A.; van Liew, M.W.; et al. SWAT: Model use, Calibration, and Validation. Am. Soc. Agric. Biol. Eng. 2012, 55, 1491–1508. [Google Scholar]
  66. Hallouz, F.; Meddi, M.; Mahé, G.; Alirahmani, S.; Keddar, A. Modeling of discharge and sediment transport through the SWAT model in the basin of Harraza (Northwest of Algeria). Water. Sci. 2018, 32, 79–88. [Google Scholar] [CrossRef] [Green Version]
  67. Buakhao, W.; Kangrang, A. DEM Resolution Impact on the Estimation of the Physical Characteristics of Watersheds by Using SWAT. Adv. Civ. Eng. 2016, 2016, 9. [Google Scholar] [CrossRef] [Green Version]
  68. Polanco, E.I.; Fleifle, A.; Ludwig, R.; Disse, M. Improving SWAT model performance in the upper Blue Nile Basin using meteorological data integration and subcatchment discretization. Hydrol. Earth. Syst Sci. 2017, 21, 4907–4926. [Google Scholar] [CrossRef] [Green Version]
  69. USDA. USDA Soil Conservation Service. National Engineering Handbook, Section 4: Hydrology; NEH Notice 4–102; USDA: Washington, DC, USA, 1972.
  70. Saxton, K.; Rawls, W.J. Soil Water Characteristic Estimates by Texture and Organic Matter for Hydrologic Solutions. Soil. Sci. Soc. Am. J. 2006, 70, 1569–1578. [Google Scholar] [CrossRef] [Green Version]
  71. Al Sayah, M.J.; Abdallah, C.; Khouri, M.; Nedjai, R.; Darwich, T. Application of the LDN concept for quantification of the impact of land use and land cover changes on Mediterranean watersheds–Al Awali basin–Lebanon as a case study. Catena 2019, 176, 264–278. [Google Scholar] [CrossRef]
  72. Périé, C.; Ouimet, R. Organic carbon, organic matter and bulk density relationships in boreal forest soils. Can. J. Soil. Sci. 2008, 88, 315–325. [Google Scholar] [CrossRef]
  73. Post, D.F.; Fimbres, A.; Matthias, A.D.; Sano, E.E.; Accioly, L.; Batchily, A.K.; Ferreira, L.G. Predicting Soil Albedo from Soil Color and Spectral Reflectance Data. Soil. Sci. Soc. Am. J. 2000, 64, 1027–1034. [Google Scholar] [CrossRef]
  74. Durand, Y.; Brun, E.; Mérindol, L.; Guyomarc’h, G.; Lesaffre, B.; Martin, E. A meteorological estimation of relevant parameters for snow models. Ann. Glaciol. 1993, 18, 65–71. [Google Scholar] [CrossRef] [Green Version]
  75. Pearson, K. Notes on regression and inheritance in the case of two parents. Proc. R. Soc. Lond. 1985, 58, 240–242. [Google Scholar]
  76. SWAT. SWAT Weather Data and Tools [Internet]. Available online: https://swat.tamu.edu/software/ (accessed on 8 March 2019).
  77. Wang, X.; Yang, W.; Melesse, A.M. Using Hydrologic Equivalent Wetland Concept within SWAT to Estimate Streamflow in Watersheds with Numerous Wetlands. Trans. Am. Soc. Agric. Eng. 2008, 51, 55–72. [Google Scholar] [CrossRef]
  78. Ahmed, S.I.; Rudra, R.P.; Gharabaghi, B.; Mackenzie, K.; Dickinson, W.T. Within–Storm Rainfall Distribution Effect on Soil Erosion Rate. ISRN Soil. Sci. 2012, 2012, 1–7. [Google Scholar] [CrossRef] [Green Version]
  79. Akale, A.T.; Dagnew, D.C.; Belete, M.A.; Tilahun, S.A.; Mekruia, W.; Steenhuis, T.S. Impact of Soil Depth and Topography on the Effectiveness of Conservation Practices on Discharge and Soil Loss in the Ethiopian Highlands. Land 2017, 6, 78. [Google Scholar] [CrossRef] [Green Version]
  80. Gis Sol; INRA; SOeS. Aléa D’érosion des Sols en 2000; Institut national de la recherche agronomique: Orléans, France, 2000. [Google Scholar]
  81. Cohen, J. A coefficient of agreement for nominal scales. Educ. Psychol. Meas. 1960, 20, 37–46. [Google Scholar] [CrossRef]
  82. Verheijen, F.G.A.; Jones, R.J.A.; Rickson, R.J.; Smith, C.J. Tolerable versus actual soil erosion rates in Europe. Earth Sci. Rev. 2009, 94, 23–38. [Google Scholar] [CrossRef] [Green Version]
  83. Kondolf, G.M. Hungry water: Effects of dams and gravel mining on river channels. Environ. Manag. 1997, 21, 533–551. [Google Scholar] [CrossRef]
Figure 1. Study area.
Figure 1. Study area.
Water 11 02526 g001
Figure 2. CORINE model framework, MFI: Modified Fournier Index [63]; Pi is the total precipitation in month i and P ¯ is the mean annual total precipitation; BGI: Bagnouls–Gaussen Index [64]; ti the mean temperature for the month i; Pi the total precipitation in month i and ki the proportion of the month i in which 2ti − Pi > 0 [62].
Figure 2. CORINE model framework, MFI: Modified Fournier Index [63]; Pi is the total precipitation in month i and P ¯ is the mean annual total precipitation; BGI: Bagnouls–Gaussen Index [64]; ti the mean temperature for the month i; Pi the total precipitation in month i and ki the proportion of the month i in which 2ti − Pi > 0 [62].
Water 11 02526 g002
Figure 3. Spatial distribution of soil erodibility parameters: (a) texture, (b) depth, (c) stoniness, and class descriptions of the Claise with their respective indices.
Figure 3. Spatial distribution of soil erodibility parameters: (a) texture, (b) depth, (c) stoniness, and class descriptions of the Claise with their respective indices.
Water 11 02526 g003
Figure 4. Soil erodibility map of the Claise classified according to the CORINE indices.
Figure 4. Soil erodibility map of the Claise classified according to the CORINE indices.
Water 11 02526 g004
Figure 5. (a) Slope angle of the Claise and (b) the corresponding CORINE indices description.
Figure 5. (a) Slope angle of the Claise and (b) the corresponding CORINE indices description.
Water 11 02526 g005
Figure 6. Potential soil erosion (PSE) risk map of the Claise basin.
Figure 6. Potential soil erosion (PSE) risk map of the Claise basin.
Water 11 02526 g006
Figure 7. (a) Current Claise vegetation cover and (b) alternative vegetation cover and corresponding CORINE classification; 1: fully protected, 2: not fully protected.
Figure 7. (a) Current Claise vegetation cover and (b) alternative vegetation cover and corresponding CORINE classification; 1: fully protected, 2: not fully protected.
Water 11 02526 g007
Figure 8. Actual soil erosion risk map of the Claise under (a) current vegetation cover and (b) the alternative pondless scenario.
Figure 8. Actual soil erosion risk map of the Claise under (a) current vegetation cover and (b) the alternative pondless scenario.
Water 11 02526 g008
Figure 9. Sediment yields of the Claise basin under (a) current conditions and (b) the pondless scenario.
Figure 9. Sediment yields of the Claise basin under (a) current conditions and (b) the pondless scenario.
Water 11 02526 g009
Table 1. CORINE erosion model and SWAT input data.
Table 1. CORINE erosion model and SWAT input data.
Input DataSourceDate Description
Topography—DEMInstitut Géographique National (IGN)—France201025 m resolution
Soil mapHarmonized World Soil Database (HWSD) [52]200830 arc-second raster database
Weather dataWeatherlink Pro 2 weather stations, source: R. Nedjai (Dynétangs)2016–2018Daily time step
Système d’Analyse Fournissant des Renseignements Adaptés à la Nivologie (SAFRAN) model (Durand et al., 1993)1970–2018
Land Use/Land Cover mapDigitized from aerial photography and cross-checked against ancillary maps2018Ortho-rectified aerial photography-Resolution: 0.5 m verified against ancillary maps: CORINE 2012 land use/land cover maps
The correspondence of the Claise basin to a natural park renders any human induced modifications on soils minor. Further, the Claise’s land occupation pattern is relatively stable. Therefore, climate is the only variable factor in the study area. Hence, the temporal difference of the utilized datasets does not cause temporal induced biases.
Table 2. Validation testing of the SAFRAN model against measured on field-data for revealing SAFRAN’s validity for use.
Table 2. Validation testing of the SAFRAN model against measured on field-data for revealing SAFRAN’s validity for use.
Weather Stations ParametersSAFRAN ParametersParameter LabelCorrelation, r
Temp OutT_QAverage temperature (°C)0.98
Hi TempTSUP_H_QMaximal temperature (°C)0.98
Low TempTINF_H_QMinimal temperature (°C)0.97
Out HumHU_QRelative humidity (%)0.97
Wind SpeedFF_QWind speed (m/s)0.83
RainPE_QEfficient rainfall (mm)0.64
Solar Rad.SSI_QIncoming solar radiation (J/cm2)0.98
Table 3. Numerical distribution of the land occupation setting of the Claise basin.
Table 3. Numerical distribution of the land occupation setting of the Claise basin.
Land Occupation ClassArea (km2)Percentage (%)CORINE Vegetation Cover Index
Clear broad—leaved forest2.210.311
Clear mixed forest0.900.131
Coniferous forest26.943.811
Dense broad—leaved forest163.5823.111
Dense mixed forest27.003.821
Field crops in medium to large terraces19.142.702
Fruit trees0.200.032
Grassland207.0429.261
Inland marshes4.010.571
Low density urban tissue3.240.462
Medium density urban tissue1.760.252
Mineral extraction site0.090.012
Non–irrigated field crops151.1721.362
Pond79.4711.230
River0.550.080
Scrubland2.800.401
Scrubland with some bigger dispersed trees15.362.171
Urban expansion site0.0102
Urban sprawl on clear wooded lands0.0102
Urban sprawl on field crops1.010.142
Urban sprawl on grassland1.200.172
Table 4. Numerical distribution of soil erodibility parameters and corresponding covered areas.
Table 4. Numerical distribution of soil erodibility parameters and corresponding covered areas.
ParameterCORINE ClassArea (km2)Percentage (%)
Soil Texture1: Slightly erodible (clay)17.532.48
2: Moderately erodible (loamy sand and sand)205.8829.12 (28.77 and 0.35)
3: Highly erodible (loam)483.5968.4
Total707100
Soil Depth1: Slightly erodible (>1000 mm)668.8294.6
2: Moderately erodible (250–750 mm)38.185.4
Total707100
Stoniness1: Fully protected (>10%)43.266.12
2: Not fully protected (<10%)663.7393.88
Total707100
Table 5. Average temperature and precipitation for the Claise basin (1970–2018).
Table 5. Average temperature and precipitation for the Claise basin (1970–2018).
MonthJanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
Average temp. (°C)4.75.27.810.414.317.920.119.716.111.87.74.8
Precipitation (mm)4166771106152617822045210818261356818482361
Table 6. Verification of the established erosion map.
Table 6. Verification of the established erosion map.
ClassesLowModerateHighValidation Points
Low7610086
Moderate1514837
High001111
Total912419134
Table 7. Actual soil erosion risk (ASE) for the Claise basin under current and simulated vegetation cover.
Table 7. Actual soil erosion risk (ASE) for the Claise basin under current and simulated vegetation cover.
Erosion RisksASE with Current Vegetation Cover—Area (km2)ASE with Simulated Vegetation Cover (Absence of Ponds)—Area (km2)ASE with Current Vegetation Cover—Percentage (%)ASE with Simulated Vegetation Cover (Absence of Ponds)—Percentage (%)
None88.237.9212.481.12
Low464.213.665.660.52
Moderate153.2754321.6876.8
High1.272152.40.1821.56
Table 8. Sensitivity analysis for calibration of the SWAT model.
Table 8. Sensitivity analysis for calibration of the SWAT model.
Parameter Namet-Statp-ValueFitted Value
R__WET_NVOL.pnd−37.520.000.823
V__GW_DELAY.gw−31.390.0052.325
R__PND_SED.pnd1.620.11318.100
V__GWQMN.gw1.590.111.461
R__SPCON.bsn−1.350.180.005
R__CN2.mgt1.290.200.564
R__PND_NSED.pnd1.230.222420.81
R__USLE_K.sol−0.970.330.209
V__ALPHA_BF.gw−0.900.370.529
R__NDTARG.pnd0.440.6649.890
R__PND_FR.pnd0.150.880.644
R__USLE_P.mgt−0.110.910.188
R__IGRO.mgt−0.040.970.080
Table 9. Shifts of sediment yields as a result of pond presence or absence.
Table 9. Shifts of sediment yields as a result of pond presence or absence.
Sediment YieldPond Presence (Area km2; Percentage %)Pond Absence (Area km2; Percentage %)
Low365.07; 63.310; 0
Moderate89.23; 15.47454.30; 78.79
High122.27; 21.20122.27; 21.20

Share and Cite

MDPI and ACS Style

Al Sayah, M.J.; Nedjai, R.; Kaffas, K.; Abdallah, C.; Khouri, M. Assessing the Impact of Man–Made Ponds on Soil Erosion and Sediment Transport in Limnological Basins. Water 2019, 11, 2526. https://doi.org/10.3390/w11122526

AMA Style

Al Sayah MJ, Nedjai R, Kaffas K, Abdallah C, Khouri M. Assessing the Impact of Man–Made Ponds on Soil Erosion and Sediment Transport in Limnological Basins. Water. 2019; 11(12):2526. https://doi.org/10.3390/w11122526

Chicago/Turabian Style

Al Sayah, Mario J., Rachid Nedjai, Konstantinos Kaffas, Chadi Abdallah, and Michel Khouri. 2019. "Assessing the Impact of Man–Made Ponds on Soil Erosion and Sediment Transport in Limnological Basins" Water 11, no. 12: 2526. https://doi.org/10.3390/w11122526

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