InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Engineering » Environmental Engineering » "River Basin Management", book edited by Daniel Bucur, ISBN 978-953-51-2605-8, Print ISBN 978-953-51-2604-1, Published: August 10, 2016 under CC BY 3.0 license. © The Author(s).

Chapter 8

Modelling Agri‐Environmental Measures for Minimizing Soil Erosion While Protecting Valuable Agricultural Land

By Matjaž Glavan, Polona Ojsteršek Zorčič and Marina Pintar
DOI: 10.5772/63380

Article top


The river Ledava basin study site.
Figure 1. The river Ledava basin study site.
A comparison of simulated and measured values for river flow (m3/s) in the river Ledava for (a) calibration period (2006–2010) and (b) validation period (2011–2013).
Figure 2. A comparison of simulated and measured values for river flow (m3/s) in the river Ledava for (a) calibration period (2006–2010) and (b) validation period (2011–2013).
A comparison of the simulated and the calculated average daily (a, b) and monthly (c, d) sediment load (ton/day) in the river Ledava (Pertoča) between June 2013 and March 2014, with (a, c) and without (b, d) one outlier value.
Figure 3. A comparison of the simulated and the calculated average daily (a, b) and monthly (c, d) sediment load (ton/day) in the river Ledava (Pertoča) between June 2013 and March 2014, with (a, c) and without (b, d) one outlier value.
Average annual sediment yield (ton/ha/year) on the HRU level (a) and the sub‐basin level (b).
Figure 4. Average annual sediment yield (ton/ha/year) on the HRU level (a) and the sub‐basin level (b).
Changes in the average annual amount of sediment transported with surface runoff from HRUs as modelled by scenarios for the river Ledava basin.
Figure 5. Changes in the average annual amount of sediment transported with surface runoff from HRUs as modelled by scenarios for the river Ledava basin.

Modelling Agri‐Environmental Measures for Minimizing Soil Erosion While Protecting Valuable Agricultural Land

Matjaž Glavan1, Polona Ojsteršek Zorčič2 and Marina Pintar1
Show details


Erosion processes in river basins and the consequent transport of sediment and sediment‐bound pollutants to reservoirs cause hydromorphological changes and eutrophication, as well as the loss of reservoir storage capacity. This chapter deals with the optimal selection and implementation of agri‐environmental measures in river basins to reduce sediment yield and load. The main aim of this was to contribute to more efficient river basin management by minimizing soil erosion, while protecting valuable agricultural land. This includes implementing measures at the most critical source areas, where they are most effective and necessary. The river Ledava basin was selected as the study area. It covers an area of 105 km2 in northeast Slovenia and southeast Austria. The results of monitoring the river Ledava discharge reveal that the average annual concentration of sediment in the water body exceeded the recommended value of 25 mg/l by 46.7%. Using the Soil and Water Assessment Tool (SWAT), we were able to determine critical source areas and simulate the effects of eight different agri‐environmental scenarios on sediment yield reduction. The results show that critical source areas comprise 12% of the river basin. Most of the scenarios reduced sediment load in the river Ledava where steeper slopes in the sub‐basin prevail and where high average annual sediment transport from hydrologic response units (HRUs) has been identified. The impact of the scenarios on the average annual sediment load (ton/year) in the river was lower than for the sediment yield (ton/ha) at the HRU level.

Keywords: erosion, sediment, agriculture, SWAT, agri‐environmental measures

1. Introduction

Soil erosion has a significant impact regarding the degradation of valuable agricultural land, where a combination of rainfall, soil type, slope properties, and land management can result in extensive loss of soil and associated nutrients [1, 2]. At the same time, erosion processes in the river basin and the transport of sediment and sediment‐bound pollutants cause hydromorphological changes and eutrophication of surface water [3, 4]. Although erosion is a natural process, the rate of soil loss is site specific and can significantly increase with inappropriate land management (e.g., removal of vegetation cover, overgrazing, fire, mineralization of organic matter, and compaction through mechanization) and climate change [5]. According to the Food and Agriculture Organization of the United Nations (FAO), erosion resulting from human activities is 100–700 times faster than the natural rate (0.1–1.0 ton/ha/year) [6]. Theoretically, soil erosion should be maintained at a rate that is equal to or below the natural rate at which new soil forms [5]. However, this balance is difficult to achieve because of differences in soil, slopes, land cover, and climate that are site dependent. Due to the Soil Strategy for England [7], estimated erosion rates in England and Wales range between 1 and 20 ton/ha/year for most agricultural fields. The mean soil loss rate in the European Union's erosion‐prone land (agricultural land, forests, and semi‐natural areas) has been found to be 2.46 ton/ha/year, resulting in a total soil loss of 970 Mt annually [8]. In Slovenia, about 4 ton/ha of soil are lost on average every year [9].

In order to mitigate the on‐site and off‐site effects of erosion processes, several control measures have been designed around the world. Measures adapted to specific farming systems to improve the environment are termed agri‐environmental measures (AEM), adopted as the instrument of the Common Agricultural Policy. The introduction of the Water Framework Directive (WFD, 2000/60/EC) has increased interest in soil conservation, mainly because of the impact of runoff and associated sediment and pollutants on water quality [1]. Policy interventions (Good Agricultural and Environmental Condition, GAEC) over the last decade in Europe have reduced the soil loss rate by 9.5% on average and by 20% for arable land [8]. Many other studies report the positive effect of measures (rotation, conservation tillage, cover crops, contour farming, vegetative filter strips, terracing, etc.) on soil loss reduction at the field scale [1, 1015] or the river basin scale [1621]. As the performance of measures is highly dependent on local circumstances (e.g., soil type, slope, crop and climate), combining different measures can lead to a greater reduction of soil loss and sediment yield [16, 19, 22].

AEM can, however, also produce negative socioeconomic effects [20], because to achieve the water quality target under the WFD, the extent of arable land should be reduced from 77.2 to 46%, pasture increased from 4 to 15% and afforestation from 10 to 21%. Similar negative economic effects were reported by [23].

To avoid negative economic effects and to achieve the maximum possible benefit, the allocation of selected measures should focus on areas where they are most effective and necessary. The optimal locations for AEM are areas that contribute disproportionately high sediment and nutrient loads, which are often referred to as critical source areas (CSAs). These areas have a particular type of soil, land use, land management, and slope and represent an overlap with areas prone to generating high volumes of runoff and erosion [24]. The use of river basin models, such as the Soil and Water Assessment Tool (SWAT) [25], can help in identifying CSAs and prioritizing areas for the cost‐effective implementation of AEM [2628].

Many studies worldwide have used SWAT to evaluate the impact of land‐use scenarios and mitigation measures on water quality and sediment yield reduction [16, 17, 19, 2932]. The assessment of AEM using SWAT can ensure the most cost‐effective allocation of measures, minimizing sediment yield, while protecting valuable agricultural land.

This chapter deals with the optimal selection and implementation of AEM in a river basin, in order to reduce sediment yields from land surfaces and sediment loads in the river. The main aim of this study was to contribute to more efficient river basin management by minimizing soil erosion while protecting valuable agricultural land. This includes the identification of CSA, where AEM are most effective and necessary.

2. Materials and methods

2.1. Study site

The river Ledava basin was selected for the study. It covers an area of 105 km2 in northeast Slovenia and southeast Austria (Figure 1). Some 33.7 km2 lies in Austria and 71.6 km2 in Slovenia. The highest elevation on the Austrian side is 598 m above sea level, in Gleichenberger Kogel, and on the Slovenian side, 418 m above sea level at the Sotinski breg. Over a third of the basin comprises arable land (37.8%), followed by forests (36.7%), meadows (12.1%), orchards (3.4%), and other land use (10%). The otherwise homogenous hilly landscape is steeper in the upper part of the river basin. Inclinations between 0–11, 11–24, 24–35, 35–50, and over 50% represent, respectively, 35.5, 37.0, 16.1, 7.9, and 3.5% of the study area. Poorly adhering tertiary and quaternary sediments, the topographical properties and the exchange of soils with low and high permeability increase the possibility of landslides and erosion processes. Landslides in this area are small in extent but numerous, with the majority (48.9%) occurring on inclinations at or above 9° or 15.4% [33, 34].

The basin area consists of higher and steeper terrain in the headwaters of the river Ledava in Austria, as a result of stiffer geological composition of phyllite carboniferous slates and basaltic tuff of volcanic origin, and homogeneous hills in Slovenia, dominated by tertiary marlstone and claystone. The bottom of the valley is covered with alluvial clay and sandy clay sediments. Almost a third (27.07%) of the area's soil consists of pseudogley district brown soils and district brown soils on Pliocene sediments, which are washed out and not very fertile. Because the soil in the river valley is relatively impermeable, the density of the river channels network (Figure 1) is considerably higher (1.77–2.35 km/km2) than the average in Slovenia (1.33 km/km2). The river Ledava (Limbach in German) is a left tributary of the Mura River and originates in the village Pichl, Steiermark, (in the Kapfenstein municipality) at an altitude of 430 m. The length of the river to the mouth of the reservoir is 17.4 km, 8 km of which runs through Austria. The ecomorphological situation [35] of the headwater part of the river's network channels is evaluated as sustainably managed. In its central section, the river is channelized along agricultural areas, with fortified banks managed as grassland and with occasional shrub cover.


Figure 1.

The river Ledava basin study site.

The river Ledava basin area has a continental climate, with average annual precipitation between 800 and 900 mm, falling mainly in the summer as showers and thunderstorms. The highest average temperatures (2003–2014) were measured in July (21°C) and the lowest in January (-0.1°C). The hydrological regime of the river is characterized by a rain–snow mixture, with the highest flows in early spring and late autumn. Low flows in the summer are the result of higher temperatures, higher evapotranspiration, and rainfall interception by vegetation. The maximum average monthly river flows (1993–2013) were measured between December and March (-0.47 to 0.42 m3/s), with the highest river flow recorded in August 2005 (29.3 m3/s). Due to the rapid hydrological response in the basin, the river Ledava flooded regularly. In order to protect the city of Murska Sobota and the cultivated land along the river, the Ledavsko jezero reservoir was built in 1977.

The results of simultaneous monitoring from June 2013 to May 2014 on the Ledava river and the Ledavsko reservoir reveal that the average annual concentration of sediment in both water bodies exceeded the recommended value (Fish Directive, 2006/44/EC) of 25 mg/l by 46.7% in the river and by 73.3% in the reservoir.

2.2. Model description

We selected the SWAT to evaluate the effects of different AEM on sediment yield and to evaluate CSAs in the river Ledava basin. SWAT is a physically based, spatially distributed basin scale model developed by the Agricultural Research Service (ARS) of the US Department of Agriculture (USDA) [36]. The main objective of the model development was to support river basin managers in evaluating the effects of alternative management decisions concerning water resources and nonpoint‐source pollution in large river basins over a longer time span [25]. In SWAT, a river basin is divided into sub‐basins based on topography and river networks. Sub‐basins are further divided into hydrologic response units (HRUs), based on unique land cover, soil, slope, and management practice combinations. Processes such as surface runoff and sediment yield are simulated for each HRU, and the contributions of each HRU are then aggregated for the sub‐basin by a weighted average. Water is then routed to the outlet of the river basin. For this study, we used the SWAT 2012 model version, ArcGIS 10.0 software and the ArSWAT 2012.10_0.15 interface.

Erosion and sediment yields were estimated for each HRU using the modified universal soil loss equation (MUSLE). Whereas USLE uses rainfall as an indicator of erosive energy, MUSLE uses the amount of runoff to simulate erosion and sediment yields. This increases the prediction accuracy of the model, eliminates the need for a delivery ratio, and provides an estimation of single storm sediment yield [37]. MUSLE is represented by Eq. (1):


where sed is the sediment yield on a given day (metric tons), Qsurf is the surface runoff volume (mm H2O/ha), qpeak is the peak runoff rate (m3/s), areahru is the area of the HRU (ha), KUSLE is the USLE soil erodibility factor [(0.013 ton/m2/h)/(m3‐ton/cm)], 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.

Each of these factors can be adjusted to represent the adoption of AEM such as slope terracing, contour farming, field crop strip sowing, and residue management. Overland flow routings in SWAT model are controlled by curve number (CN) and coefficient of Manning's roughness (n). AEM can be represented by modifications of CN and n values affecting decrease in surface runoff by increasing infiltration (e.g., contour farming, slope terracing, and field crop strip sowing) and flow rate decrease by runoff interception (e.g., crop residue management, and field crop strip sowing), respectively, [38]. The effectiveness of vegetated filter strips (VFS) varies for different forms of sediments and nutrients. A single VFS is for modelling of concentrated flow split into two segments. In segment one, 90% of the VFS area obtains the least flow and segment two with 10% of area obtaining the major runoff (25–75%). Three parameters can be modified in the SWAT model to adapt modelled VFS design. The first one is drainage of the spatial area‐to‐VFS area ratio (DAFSratio), the second is percentage of the field drained by the most heavily loaded 10% of the VFS (DFcon), and the third is percentage of the flow through the most heavily loaded 10% of the VFS that is fully channelized (CFfrac) [28]. AEM can be simulated for specific dates and by explicitly defining the appropriate management parameters for each HRU.

2.3. Database and data analysis

Recently available geographical information system (GIS) maps for topography, land use, and soils for Slovenia and Austria were used to represent the area (Table 1). The stream network was defined based on a predefined stream network obtained from the Environment Agency of the Republic of Slovenia (ARSO). The topographic information was used for automatic delineation of the river basin. Land use and soil maps were superimposed over the basin.

Data typeScaleSourceData description
Topography (DEM raster)Slovenia: 25 m
Austria: 1 m
The Surveying and Mapping Authority of the Republic of Slovenia, GIS-Steiermark, and GIS-BurgenlandElevation, overland and channel slopes, lengths
SoilsSlovenia: 1:25000
Austria: 1 km Raster
Ministry of Agriculture, Forestry and Food of the Republic of Slovenia; Biotechnical Faculty (University of Ljubljana) Austrian Research Centre for ForestsSpatial soil variability, soil types, and properties
Land useSlovenia: 1 m Raster (Graphical Units of Agricultural Land)
Austria: 1:5000
Slovenia: Ministry of Agriculture, Forestry and Food of the Republic of Slovenia
Austria: GIS Steiermark, GIS Burgenland, and Municipalities (St. Anna am Aigen, Kapfenstein, Neuhaus am Klausenbach)
Land use, Land cover classification, and spatial representation
Land management information/Chamber of Agriculture and Forestry of Slovenia-Agricultural advisory service in CankovaCrop rotations (harvesting, planting, management), fertilizer application (rates and time)
WeatherSlovenia 3 and
Austria 2 stations
Environment Agency of the Republic of Slovenia (ARSO), Zentralanstalt für Meteorologie und Geodynamik (ZAMG), AustriaDaily precipitation, Temperature (max., min.), relative humidity, wind, solar radiation from 2003-2014
River discharge1 PointEnvironment Agency of the Republic of SloveniaDaily flow data (m3/s) from 2003 to 2014
Waste water treatment plantsSlovenia: 3
Austria: 2
Environment Agency of the Republic of Slovenia; Land Steiermark-Amt der Steiermärkischen LandesregierungAverage daily discharge of orgP, sediment and orgN
Water qualityTwo monitoring stations (Pertoča
and Ledavsko jezero)
Bi-weekly monitoringTSS, NO3-, NO2-, PO42-TP, TN, temperature, dissolved O2 (2013-2014)

Table 1.

Model input data sources for the river Ledava sub-basin.

Data for 11 years (2003–2014) of daily observed precipitation, the temperature from five stations and the relative humidity, wind speed, and solar radiation data from two stations were obtained from the ARSO and Zentralanstalt für Meteorologie und Geodynamik (ZAMG) in Austria. In order to better understand the endpoint and transport of sediment and nutrients, data for typical management practices such as crop growth, fertilizer application, and tillage operations for different land use were gathered from the Chamber of Agriculture and Forestry of Slovenia—Agricultural Advisory Service, in Cankova.

Because the water quality data provided by ARSO were insufficient, bi‐weekly water quality monitoring was carried out for 1 year (2013–2014). Concentrations of total suspended sediments (mg/l) were measured bi‐weekly, from June 2013 to May 2014, at the sampling point of Pertoča (46°46′26,52" N 16°2′24,64" E), located on the river Ledava about 1.5 km upstream of the Ledavsko jezero reservoir. In total, 94 samples were taken and analyzed at the National Laboratory of Health, Environment and Food and Erico d.o.o.

2.4. Model setup and evaluation

For this study, the SWAT model was built up by dividing the river Ledava basin into 18 sub‐basins (Figure 1) and 5758 HRUs. A warm‐up period of 3 years (2003–2005) was used to initialize the model. For the calibration, data from 2006 to 2010 were used for flow calibration and from 2011 to 2013 for validation. For sediment calibration, data from 2013 to 2014 was used. Due to the small amount of measured sediment data over 1 year (June 2013–June 2014), validation for sediment was not performed.

The first step in the calibration and validation process in SWAT is to determine the most sensitive parameters for a given basin. Sensitivity analysis is helpful for identifying and ranking parameters that have a significant impact on specific model outputs, such as streamflow or sediment [39]. For the sensitivity analysis and calibration, special software called SWAT‐CUP is used, which offers semi‐automatic or combined manual/automatic calibration. Among the algorithms included in SWAT‐CUP, the most efficient to calibrate a project is Sequential Uncertainty Fitting (SUFI‐2), the use of which has been found to require fewer simulations to complete a calibration/uncertainty project [40] and is highly recommended for the calibration of SWAT models [25].

In this study, sensitivity analysis was performed using measured data for the river Ledava. The analysis was carried out for average daily flow and sediments. We used the default lower and upper bound parameter values for all parameters [37]. Table 2 illustrates the parameters that have the greatest impact on the model when they are changed. The sensitivity analyses demonstrate the great importance of the parameters that are associated with surface runoff and snow (SFTMP, SMTM, TIMP, SOL_AWC, SURLAG, and CN2) and channel routing (CH_BNK_BD, CH_COV, SPEXP, and SPCON).

The most sensitive parameters were used to perform model calibration. Calibration is a process to better parameterize a model to a given set of local conditions, thereby reducing prediction uncertainty. During calibration, parameters are varied within an acceptable range until a satisfactory correlation is achieved between measured and simulated data [18, 25]. In this study, calibration of the model was performed for flow and sediment on a daily and monthly basis.

ParametersDefinitionProcessRangeSensitivity rank
CN2.mgtSCS runoff curve number for moisture
condition 2
SURLAG.bsnSurface runoff lag timeRunoff0-122
SOL_AWC(...).solAvailable water capacity of the soil layerSoil0-0.33
TIMP.bsnSnow pack temperature lag factorSnow0.01-14
SFTMP.bsnSnowfall temperatureSnow-5 to 55
SMTMP.bsnSnow melt base temperatureSnow-5 to 56
SedimentSPCON.bsnCalculating maximum amount of
sediment reentrained
SPEXP.bsnCalculating sediment reentrained in
channel sediment routing
CH_COV1.rteChannel cover factorChannel-0.001 to 13
CH_BNK_BD.rteBulk density of channel bank sedimentChannel1.1-1.94

Table 2.

Parameters identified during sensitive analysis to be sensitive to changes.

The model performance for river flow and sediment base simulation was tested with graphical comparison and with the help of objective functions commonly used in hydrological studies, such as the Pearson coefficient of correlation (R2), Nash–Sutcliffe simulation efficiency (ENS), and Percent bias (PBIAS) [41]. Several authors [20, 36, 42] provide further discussions regarding the strengths and weaknesses of using R2, ENS, and PBIAS.

2.5. Scenarios for AEM

In order to achieve “good ecological potential” in accordance with the Water Framework Directive (2000/60/EC) and Fish Directive (2006/44/EC), management measures have to be undertaken. Choosing measures often present a problem, because their effectiveness will depend on the type and source of pollution, as well as on the size and characteristics of the drainage area. At the same time, measures should be cost‐effective, with minimum implementation and maintenance costs, with long‐term effectiveness and without requiring major interventions in the area.

Based on a literature review [43], a list of 13 erosion AEM was proposed (cover crops, crop rotation management, conservation tillage, strip cropping, residue management, grassland conversion, vegetative filter strips, grassed waterways, riparian buffer strips, terracing, sedimentation basins, bank erosion protection, and check dams). Due to the heterogeneous characteristics of the study area, climatic conditions, ways of implementing the measures, different methods of measurement, capabilities of the model, and statements in relevant literature, we considered a set of measures as guidelines for the selection of scenarios. For the final selection of measures, we took into consideration the following five criteria:

  1. Land use, crop production management, and technologies;

  2. Landscape characteristic (soil type, slope, and farmland area);

  3. AEM already implemented in the river basin;

  4. Model requirements;

  5. Data availability.

Scenario/parametersUsed valueDefault valueRange
S1-Vegetative filter strip on 0-11% slope
S2-Vegetative filter strip on 11-24% slope
S3-Vegetative filter strip on 0-24% slope
 VFSRATIO60, 40100-300
S4-Conservation tillage
 EFTMIX0.3Equipment based0-1
 DEPTIL 150 (field cultivator)150Equipment based0-750
 DEPTIL 220 (field cultivator)220Equipment based0-750
S5-Strip cropping or contour planting on 11-24% slope
 CN2-3 unitsHRU based0-100
 USLE_P0.801-1 to 1
S6-Terraces on 11-24% slope
 CN2-6 unitsHRU based0-100
S7-No cover crops in winter
 Winter cover crop has been removed from rotation
S8-Cover crops every year
 After harvest in fall winter cover crops (red clover) was applied. In spring leftover was plough into a ground as green manure

Table 3.

Parameters and values used in the scenario simulations for the river Ledava basin.

From the 13 measures, we excluded those for which effectiveness increases in line with an increasing area of implementation (grassed waterways and sedimentation pools). Such measures are unsuitable for the study area, due to its scattered land parcel structure. Because of the important proportion of agricultural land on slopes >11%, we chose terraces, although they are not typical for this landscape and are rarely used. According to data from 2014 on the implementation of AEM in the Rural Development Program (RDP) under the EU Common Agricultural Policy (CAP) obtained from the Agency of the Republic of Slovenia for Agricultural Markets and Rural Development, three measures have also been implemented in the area to improve soil fertility and reduce soil erosion. These are cover crops, crop rotation management, and winterkill cover crop. We also selected vegetative filter strips, because the literature review indicates that these are among the most commonly used measures to reduce the effects of diffuse pollution. Five different measures were finally selected, from which eight modelling scenarios were built (Table 3).

We simulated measures for areas that are the most sensitive to erosion, such as arable land and drained arable land on slopes from 0 to 11% and from 11 to 24%. The aim of the scenarios was to find the most effective measure or combination of measures to effectively reduce soil erosion, sediment transport into the river and sediment flow into the downstream lake. Given the fact that orchards represent only 2.46% and vineyards 0.82% of the basin area, we decided to simulate the measures only for the arable land. This comprises the most critical areas in terms of soil erosion and covers 38.02% of the basin area. The current average arable land management practice in the basin involves a 7‐year crop rotation.

3. Results and discussion

3.1. Calibration and validation

3.1.1. Hydrology calibration and validation

Objective functions show that the simulated total flow is within the acceptable range (Table 4; Figure 2). To achieve acceptable calibration and validation results, a list of model parameters was changed from default to final values (Table 5). Table 4 lists the calibration and validation values for the model performance for flow. Values >0.5 for the coefficient of determination (R2) are considered acceptable, while negative PBIAS values indicate a small overestimation of the simulated values. Nash–Sutcliffe efficiency (ENS) on monthly and daily time steps are in the acceptable range [41]; however, the ENS coefficient is very sensitive to values that stand out from the average [42]. SWAT simulated the streamflow trends very well, as simulated streamflow values do not exceed the measured streamflow data by more than 15% [41]. Comparing simulations run under different time steps shows that this element is important for understanding model performance [44]. Results show that on daily and yearly time steps, model simulations for streamflow are between good and very good in terms of correlating with measured values. The model simulations are weaker for monthly time steps. However, as our results will be presented as a yearly average, we can state that our model is calibrated sufficiently to be used for further scenario simulations. After the base model calibration was completed, the parameters remained fixed for further use in scenario modelling.

Objective functionCalibration (2005-2010)Validation (2011-2013)

Table 4.

Statistical values for the annual, monthly and daily time step calibration for river flows of the river Ledava (2006-2010).


Figure 2.

A comparison of simulated and measured values for river flow (m3/s) in the river Ledava for (a) calibration period (2006–2010) and (b) validation period (2011–2013).

River flow
 ALPHA_BFBaseflow alpha factor0-10.0480.5972
 CANMXMaximum canopy storage0-1000-0.2017
 CN2SCS runoff curve number for moisture condition 235-98Default-0.0753
 ESCOSoil evaporation compensation factor0-10.950.88
 GW_DELAYGroundwater delay0-5003190.35
 GW_REVAPGroundwater “revap” coefficient0.02-
 GWQMNThreshold depth of water in the shallow aquifer
required for return flow to occur
 SFTMPSnowfall temperature-5 to 511.66
 SMFMNMinimum melt rate for snow during the year
(occurs on winter solstice)
 SMFMXMaximum melt rate for snow during year
(occurs on summer solstice)
 SMTMPSnow melt base temperature-5 to 50.51.5933
 SNOCOVMXMinimum snow water content that corresponds
to 100% snow cover
 SOL_AWCAvailable water capacity of the soil layer0-1default0.1917
 SOL_KSaturated hydraulic conductivity (mm/h)0-2000default0.1883
 SURLAGSurface runoff lag time0.01-2440.2617
 TIMPSnow pack temperature lag factor0-110.1321
 ADJ_PKRPeak rate adjustment factor for sediment routing in
the sub-basin (tributary channels)
 PRF_BSNPeak rate adjustment factor for sediment routing
in the main channel
 SPCONLinear parameter for calculating the maximum amount
of sediment that can be reentrained during channel
sediment routing
 SPEXPExponent parameter for calculating sediment reentrained
in channel sediment routing
 CH_COV1Channel erodibility factor-0.05 to 3005.4
 CH_COV2Channel cover factor-0.001 to 3005.4
 CH_ERODMOChannel erodibility factor. Jan-Dec0-100.3

Table 5.

The parameters, range (min. and max. value), default and final values used in the model for the flow and sediment calibration and validation of the river Ledava basin.

Due to the smaller volume of data covered, the performance results for the model validation period are often lower than for the calibration period [25]. Nevertheless, the results of the validation in this study are in line with the calibration results. ENS values during calibration (0.57) and validation (0.50) indicate a good agreement between measured and simulated values, as demonstrated by the coefficient of determination (R2 > 0.5). Only the PBIAS results vary between calibration and validation values. For calibration of the model, simulated values were overestimated, whereas for validation, the values were slightly underestimated. The PBIAS for validation is in the range of good to very good related to the calibration results.

3.1.2. Sediment calibration

Objective functions for a daily and monthly time step sediment calibration (Table 6; Figure 3) show that the model is acceptable for predicting sediment. Statistics significantly improve when only one outlier value is not included in the calculation. Due to the limited amount of data for sediment, with only 1 year of measurements, validation was not carried out. Sediment load was calculated as the product of the measured concentrations of sediment (mg/l) and the average daily and monthly flows (m3/s). The coefficient of determination (R2) indicates a moderate to strong positive agreement for daily and monthly time steps. The Nash–Sutcliff coefficient (ENS) only exceeds 0.5 for monthly time steps when the outlier value (17. 2. 2014 amounted to 259.8 ton/day) is excluded. PBIAS describes the 18% deviation of the results for a daily time step, which is in the range for very good model performance [41]. The performance in simulation of trends with ENS values (ranging from 0.26 to 0.11 for monthly and daily time steps, respectively) reveals that the simulation was unsatisfactory. However, this is because the efficiency coefficient is sensitive to extreme values and yields sub‐optimal results when the dataset contains outliers [41]. The performance for R2 and ENS improves to a satisfactory level for daily time steps following the removal of outlier values. A better prediction of ENS could be expected with a larger number of data samples over longer period of time [45].

Objective functionSediment load (ton/day)
All valuesNo outlier valueAll valuesNo outlier value

Table 6.

Objective functions for daily and monthly values of sediment concentration (mg/l) and sediment load (ton/day) in the river Ledava at the Pertoča observation station (June 2013 to March 2014).

Figure 3 shows the impact of one outlier value on the model statistical performance results for sediment. This value is a result of heavy precipitation (an average daily rainfall of 17.14 mm from the previous day to midnight at the nearby meteorological station at Cankova). This was a local storm with a local effect on sediment load. The impact of the outlier value is also evident from the difference in the coefficients of ENS, R2 and PBIAS. In all cases, the results were substantially improved when the outlier value was removed from the dataset, which is understandable considering the effect of this value on the average.


Figure 3.

A comparison of the simulated and the calculated average daily (a, b) and monthly (c, d) sediment load (ton/day) in the river Ledava (Pertoča) between June 2013 and March 2014, with (a, c) and without (b, d) one outlier value.

3.2. Base scenario

Through analysis of the base scenario, we determined the CSAs, that is, where the source and transport areas that are connected to water bodies coincide. HRUs, where the annual average sediment yield exceeded 0.5 ton/ha, were considered to be CSAs (Figure 4). It turns out that CSAs represent only 12.2% of the entire river Ledava basin, whereas the most intense erosion occurs on arable land and encompasses 88.2% of all CSAs. These areas have unstable geology and soil types susceptible to soil erosion and for that reason are strongly subject to the action of weather forces [34]. The spatial variability of sediment yield is influenced by many factors, particularly land use, management practice, slope, and soil characteristics. The impact of land management and its graphical representations greatly depend on the accuracy of the spatial and attribute data used in the modelling process. The average annual sediment yield at the HRU level was divided into five classes (Figure 4).

The average annual sediment yield (0.28 ton/ha) in the river basin calculated with SWAT coincides with earlier calculations for this area based on land use and slope recalculations [34, 46]. The source of sediment in this river basin is spatially heterogeneous and on average amounts to between 0.3 and 0.9 ton/ha/year. This study shows that in certain HRUs, sediment yield can reach up to 4.10 ton/ha/year (Table 7; Figure 4). This value can be exceeded during periods of heavy rainfall, such as in 2009 when the highest value in the river reached 18.17 ton/ha. The most critical areas, with a sediment yield of 0.5 ton/ha, comprise an area of 1260 ha or 11.97%. They can be found on slopes between 11 and 24 and over 50%, and on all types of agricultural land. These areas contribute 60.9% of all the sediment transported from the HRUs into the river. The highest amount of sediment is transported from arable fields with drainage systems (0.60 ton/ha/year), other arable fields (0.40 ton/ha/year), and vineyards (0.37 ton/ha/year). Among the most erodible soil types are heavy gley and pseudogley soils. Sub‐basins (1, 2, 4 and 5) with a maximum transport of sediment at the HRU level are located in Austria, where the area is characterized by intensive agricultural land, heavy soils, and steep topography. Soils with a high proportion of clay and without vegetation cover and steep topography are more erodible [5, 9, 34, 30].


Figure 4.

Average annual sediment yield (ton/ha/year) on the HRU level (a) and the sub‐basin level (b).

Erosion class (ton/ha)Area (ha)Area from total (%)Average sediment yield (ton/ha)

Table 7.

Average annual values of sediment yield (ton/ha) at the HRU level presented with five classes (2006-2013).

3.3. Scenario evaluation

Evaluation of the measures used in the scenarios is dependent on achieving set goals and benchmarks. In the context of this research, the main objective was to reduce soil erosion and improve river water quality, which can be achieved by mitigating sediment yield and transport into the river Ledava.

New scenarios did not have any important impact on the predicted average annual flow in the river Ledava. The main reason lies in the configuration of the model parameters, because the modelling runoff curve CN was set at a uniform value for the entire river basin area. Among other things, the vegetation buffer zone in SWAT does not affect the surface runoff [37]. The small impact of measures used in the study on the river flow has also been confirmed in other studies [47, 48]. Whereas scenario S2 increased the simulated average flow rate, scenario S8 is the only one that decreased it. The review of average monthly values shows that scenarios S1, S2, S4, S5, and S6 reduce the river flow between May and October. This was influenced by increased evapotranspiration and intercepted precipitation by the vegetation canopy during the growing season [14]. Low flow rates between October and April result from the growth of cover crops on agricultural land (S8), which influences soil water retention. At the same time, cover crops can improve soil properties and can increase organic carbon content, soil stability, and the infiltration rate [49].

The sediment reduction efficiency of scenarios was estimated by comparing the simulation results of the base scenario with AEM scenarios (Tables 8 and 9). The efficiency (r) of the scenario was calculated according to [29] as follows:


where ybase represents the base scenario sediment transport (SYLD in tonnes/ha/year) and yscenario represents sediment transport (SYLD in tonnes/ha/year) according to the simulation results of the agri‐environmental scenarios. The calculation is presented as a percentage change.

Sediment transport from the HRU is reduced in all cases, except for scenario S7 (Table 8; Figure 5). It is apparent that the sediment yield (ton/ha/year) reduction at the level of the HRU is most affected by vegetative filter strips (S2: 43.4% and S3: 56.1%) and terraces (S6: 42.4%), followed by conservation tillage (S4: 20.3%), contour tillage (S5: 18.9%), and cover crops (S8: 11.9%). A combination of vegetative filter strips (S3) on slopes up to 11% (S1) and from 11 to 24% (S2) is much more effective than any other measure. A very similar impact on sediment yield reduction at HRU level can be observed between scenarios S2 and S6, although they greatly differ in terms of financial and labor input. All the scenarios were most effective in the sub‐basin dominated by agriculture land with slopes of up to 24%.

Most of the scenarios reduced sediment load in the river Ledava where steeper slopes in the sub‐basin prevail and where high average annual sediment transport from the HRU was identified (Table 9). The impact of the scenarios on the average annual sediment load (ton/year) in the river was lower than at the HRU level. In this regard, scenarios S3 and S6 were also more effective than other measures (S4, S5, and S8): A highly similar order in the impact of the scenarios on sediment load in the river was achieved, although with less variation from the baseline scenario. Scenario S3 was the most efficient at reducing sediment load, followed by S2 and S6. Unlike S6, vegetative filter strips were more efficient in the lowland basin, with a smaller transport of sediment from the HRU. The importance of cover crops is observed in the S7 scenario, where cover crops were removed from rotation. Overall, negative changes appeared in comparison with the baseline scenario impacting on increase in the sediment yield at HRU level and sediment load in the river Ledava at the outlet of sub‐basins.

Sub-basinAreaSediment yield
Scenarios percentage change in sediment yield (%)

Table 8.

The impact of the agri-environmental scenarios (S) on the average annual sediment yield in the sub-basin level expressed as a percentage (%) change to the baseline scenario (2006-2013) for the river Ledava basin.

[i] - The five largest changes in the scenarios are shown in bold.

In sub‐basins 1, 2, and 4, where erosion processes are most prominent, the simulated implementation of terraces on slopes between 11 and 24% decreased the amount of sediment in the river by 51, 51.4, and 54.1% (Table 8). By contrast, the effect of terraces on reducing the sediment transport from the HRUs (Table 8) was the largest in sub‐basin 15 (62%), where average annual sediment transport was 0.76 ton/ha. The average annual sediment transport from the HRUs in sub‐basins 11, 14, and 17 was relatively small (0.11–0.13 ton/ha), where the maximum amounts of sediment were transported from areas with slopes between 11 and 24% (from 0.40 to 0.49 ton/ha/year). On average, scenarios S5 and S8 had the smallest impact on sediment load reduction in the river (Table 9). Results show that the effect of inferior non‐construction measures, such as conservation tillage (S4) or contour or strip cropping (S5), achieve a smaller reduction effect than terraces (S6) or vegetation buffer zones (S1–3) [36]. Scenarios S7 and S8 show the importance of winter cover crops on soil erosion processes and sediment yield transport at the HRU level, as well as on simulated sediment load and concentration in the river Ledava.

Sub-basinScenario impact on sediment load at outflow
ton/yearPercentage (%) change from base

Table 9.

The impact of the scenarios on the average annual sediment load (ton/year) at the outlets from the sub-basins and at the inflow into sub-basins 18, expressed as a percentage (%) change from the baseline scenario (2006-2013).

[i] - The five largest changes in the scenario are shown in bold.

We conclude that the reduction in sediment transported at the sub‐basin level was an effect of the proper placement of modelled AEM in CSAs where soil erosion and the sediment transport are highly evident, and also an effect of a higher proportion of forest, which may reduce excessive runoff [50]. The differences between the effect of terraces on the sediment transport at HRU level and the sediment yield in the river at the outlet from a sub‐basin have also been reported in previous studies [20]. These situations, where AEM are more effective at the HRU level than at the sub‐basin outlet or whole‐basin outlet, are due to morphology of the terrain and soil properties [5153].


Figure 5.

Changes in the average annual amount of sediment transported with surface runoff from HRUs as modelled by scenarios for the river Ledava basin.

4. Conclusions

According to the agricultural, environmental, and water policies in member states of the EU‐28, AEM should be implemented according to need, the willingness of farmers to adopt the measures, existing management practice and the natural conditions in the relevant region. This study reveals the measures that are relevant to address the protection of agricultural soil in the selected river Ledava basin study area. The fact that the measures are more effective when they are introduced at a smaller scale (HRU level) and are necessary only in the CSAs, is important in terms of protecting valuable agricultural land in situ.

The main uncertainties in our research are connected with obtaining data and its comparability between two countries. The precision of all spatial layers varied greatly, as not one layer was at the same scale for the two countries. Wherever possible, we upscaled or downscaled the layers. To overcome this issue, countries in future should establish a common data protocol for transboundary river basins. One of the everlasting problems is state monitoring of water quality sampling. Monitoring is spatially rare, with wide time steps between sampling. There are frequent changes to the spatial position of the sampling point, river flow and water quality sampling points are often in different places, and in certain cases monitoring stops after few years. To overcome this, we carried out 1‐year sampling of sediment concentrations at several points along the river. However, this required our own financial and working power resources, which are limited. We recommend that the relevant state environmental agency responsible for surface water measurements should pay more attention to the location of monitoring points and should ensure long‐term, regular, and comparable time series. Calibration of the model results for soil erosion and sediment transport was done only with comparison to measured values in the river streamflow. Although we had three general maps representing erosion rates at a different spatial scale to previous studies [8, 9], actual soil erosion was never measured in the area. We propose a series of on‐site field measurements of soil erosion at different slope, land usage, and soil types to evaluate the results of this and other similar studies.

One of the important factors for the effectiveness of AEM is climate change. Climate change already dictates differing land use, crop types, management, and production technology. This is already affecting the processes of mass transfer from catchment areas into water bodies, and we can expect that these processes will escalate in the future. These effects can only be determined with additional modelling of climate change impact. For example, currently foreseen AEM have been found to be less effective under future climate change scenarios, however not to the extent of being unnecessary.

The criteria for individual parameters that define the AEM in the model are necessary in order to assess, evaluate, and select alternatives with regard to the objectives of the decision‐maker. Therefore, it is necessary to supplement the missing limit values for the water quality parameters, in particular those related to soil erosion, sediment yield, and sediment loads in surface water bodies.

The results show that CSAs occupy 12.2% of the river basin. Most of the scenarios predict reduced sediment in the river Ledava where steeper slopes in the sub‐basin prevail and where high average annual sediment transport from the HRUs is identified. The impact of the scenarios on the average annual sediment load (ton/year) in the river was lower than at the HRU level. An extremely similar order in the impact of the scenarios on total sediment load in the river was achieved, although with less change from the baseline scenario. Overall, negative changes and or reductions appeared according to the baseline scenario in the sediment yield at HRU level and sediment in the river Ledava at the outlet of the sub‐basin. This study enables us to better understand the processes of selecting environmentally effective measures and their contribution to the long‐term improvement of the surface water bodies’ ecological status in accordance with the Water Framework Directive. This also opens up the potential for selecting cost‐effective and socioeconomically friendly solutions addressing the rural development programs of EU member states under the Common Agricultural Policy.


This work was funded by the SPIRIT Slovenia, Public Agency of the Republic of Slovenia for the Promotion of Entrepreneurship, Innovation, Development, Investment and Tourism, and Savaprojekt Ltd.


1 - Posthumus H, Deeks LK, Fenn I, Rickson RJ. Soil Conservation in two English catchments: linking soil management with policies. Land Degrad Dev. 2011 Jan–Feb;22(1):97–110.
2 - Walling DE, Collins AL, Jones PA, Leeks GJL, Old G. Establishing fine‐grained sediment budgets for the Pang and Lambourn LOCAR catchments, UK. J Hydrol. 2006 Oct 30;330(1–2):126–41.
3 - Rickson RJ. Can control of soil erosion mitigate water pollution by sediments? Sci Total Environ. 2014 Jan 15;468:1187–97.
4 - [4] Annandale GW. Reservoir Sedimentation. Encyclopedia of Hydrological Sciences. Chichester, UK; John Wiley & Sons, Ltd; 2006.
5 - Morgan RPC. Soil erosion and conservation. Malden; Oxford; Carlton: Blackwell Publishing; 2005.
6 - Roose E. Land Husbandry: Components and Strategy. Rome: FAO; 1996.
7 - DEFRA. Safeguarding Our Soils: A Strategy for England. London: Department for Environment, Food & Rural Affairs; 2009. Available from:‐soil‐strategy‐090910.pdf
8 - Panagos P, Meusburger K, Van Liedekerke M, Alewell C, Hiederer R, Montanarella L. Assessing soil erosion in Europe based on data collected through a European network. Soil Sci Plant Nutr. 2014 Jan 2;60(1):15–29.
9 - Komac B, Zorn M, Fridl J. Soil erosion on agricultural land in Slovenia—measurements of rill erosion in the Besnica valley = Erozija prsti na kmetijskih zemljiscih v Sloveniji—meritve zlebicne erozije v dolini Besnice. Acta Geogr Slov. 2005(45):53–86.
10 - Wei W, Chen LD, Zhang HD, Yang L, Yu Y, Chen J. Effects of crop rotation and rainfall on water erosion on a gentle slope in the hilly loess area, China. Catena. 2014 Dec;123:205–14.
11 - Espejo‐Perez AJ, Rodriguez‐Lizana A, Ordonez R, Giraldez JV. Soil loss and runoff reduction in olive‐tree dry‐farming with cover crops. Soil Sci Soc Am J. 2013 Nov–Dec;77(6):2140–8.
12 - Stevens CJ, Quinton JN, Bailey AP, Deasy C, Silgram M, Jackson DR. The effects of minimal tillage, contour cultivation and in‐field vegetative barriers on soil erosion and phosphorus loss. Soil Tillage Res. 2009 Dec;106(1):145–51.
13 - Yang Q, Meng FR, Zhao ZY, Chow TL, Benoy G, Rees HW, et al. Assessing the impacts of flow diversion terraces on stream water and sediment yields at a watershed level using SWAT model. Agr Ecosyst Environ. 2009 Jul;132(1–2):23–31.
14 - Zhang JH, Frielinghaus M, Tian G, Lobb DA. Ridge and contour tillage effects on soil erosion from steep hillstopes‐in the Sichuan Basin, China. J Soil Water Conserv. 2004 Nov–Dec;59(6):277–84.
15 - Posthumus H, De Graaff J. Cost‐benefit analysis of bench terraces, a case study in Peru. Land Degrad Dev. 2005 Jan–Feb;16(1):1–11.
16 - Strauch M, Volk M. SWAT plant growth modification for improved modeling of perennial vegetation in the tropics. Ecol Model. [Article]. 2013 Nov;269:98–112.
17 - Kaini P, Artita K, Nicklow JW. Optimizing structural best management practices using SWAT and genetic algorithm to improve water quality goals. Water Resour Manag. 2012 May;26(7):1827–45.
18 - Glavan M, Pintar M. Modelling of surface water quality by catchment model SWAT. In: Kumarasamy M, editor. Studies on Water Management Issues. Rijeka, Croatia: InTech; 2012a.
19 - Lam QD, Schmalz B, Fohrer N. The impact of agricultural Best Management Practices on water quality in a North German lowland catchment. Environ Monit Assess. 2011 Dec;183(1–4):351–79.
20 - Tuppad P, Kannan N, Srinivasan R, Rossi CG, Arnold JG. Simulation of agricultural management alternatives for watershed protection. Water Resour Manag. 2010 Sep;24(12):3115–44.
21 - Volk M, Liersch S, Schmidt G. Towards the implementation of the European Water Framework Directive? Lessons learned from water quality simulations in an agricultural watershed. Land Use Policy. 2009 Jul;26(3):580–8.
22 - Bosch NS. The influence of impoundments on riverine nutrient transport: An evaluation using the Soil and Water Assessment Tool. J Hydrol. 2008 Jun 20;355(1–4):131–47.
23 - Posthumus H, Deeks LK, Rickson RJ, Quinton JN. Costs and benefits of erosion control measures in the UK. Soil Use Manag. 2015 Oct;31:16–33.
24 - Pionke HB, Gburek WJ, Sharpley AN. Critical source area controls on water quality in an agricultural watershed located in the Chesapeake Basin. Ecol Eng. 2000 Apr;14(4):325–335.
25 - Arnold JG, Moriasi DN, Gassman PW, Abbaspour KC, White MJ, Srinivasan R, et al. Swat: model use, calibration, and validation. Trans ASABE. 2012 Jul–Aug;55(4):1491–508.
26 - Panagopoulos Y, Makropoulos C, Mimikou M. Diffuse surface water pollution: driving factors for different geoclimatic regions. Water Resour Manag. 2011 Nov 1;25(14):3635–60.
27 - Ghebremichael LT, Veith TL, Watzin MC. Determination of critical source areas for phosphorus loss: Lake Champlain basin, Vermont. Trans ASABE. 2010 Sep–Oct;53(5):1595–604.
28 - White WJ, Morris LA, Pinho AR, Jackson CR, West LT. Sediment retention by forested filter strips in the Piedmont of Georgia. J Soil Water Conserv. 2007 Nov–Dec;62(6):453–63.
29 - Arabi M, Frankenberger JR, Enge BA, Arnold JG. Representation of agricultural conservation practices with SWAT. Hydrol Process. 2008 Jul 30;22(16):3042–55.
30 - Glavan M, White SM, Holman IP. Water quality targets and maintenance of valued landscape character—experience in the Axe catchment, UK. J Environ Manag. 2012 Jul 30;103:142–53.
31 - Ullrich A, Volk M. Application of the Soil and Water Assessment Tool (SWAT) to predict the impact of alternative management practices on water quality and quantity. Agric Water Manag. 2009;96(8):1207–17.
32 - Rousseau AN, Savary S, Hallema DW, Gumiere SJ, Foulon E. Modeling the effects of agricultural BMPs on sediments, nutrients, and water quality of the Beaurivage River watershed (Quebec, Canada). Can Water Resour J. 2013;38(2):99–120.
33 - Natek K. Endangerment and regenerative capacities of the surface in northeast Slovenia, defined on the basis of surface regeneration in Haloze after the storms of July 1989. Possibilities of regional and spatial development of Spodnje Podravje and Prlekija Region. Conference proceedings. Ljubljana, Slovenia: Association of the Geographical Societies of Slovenia, VB & S; 1996.
34 - Komac B, Zorn M. Landslides in the Pomurje region. Pomurje Region – Sustainable regional development along the Mura River. Conference proceedings. Ljubljana, Slovenia: Association of the Geographical Societies of Slovenia, Littera picta; 2009.
35 - Šajn‐Slak A, Smolar‐Žvanut N, Vrhovšek D, Krivograd‐Klemencčič A. = Ecohydrological research and conceptual solutions for ecoremediation in the upper carchment area of the Ledava river. International conference on ecoremediation. Maribor, Slovenia; Universty of Maribor, Faculty of Philosophy, International Centre for Ecoremediation; 2008.
36 - Gassman PW, Reyes MR, Green CH, Arnold JG. The soil and water assessment tool: historical development, applications, and future research directions. Trans ASABE. 2007 Jul–Aug;50(4):1211–50.
37 - Neitsch SL, Arnold JG, Kiniry JR, Williams JR. Soil and water assessment tool theoretical documentation—Version 2005. Texas, Temple: Agricultural Experiment Station, Blackland Research Center, Agricultural Research Service, Grassland, Soil and Water Research Laboratory; 2005.
38 - Xie H, Chen L, Shen ZY. Assessment of agricultural best management practices using models: current issues and future perspectives. Water. 2015 Mar;7(3):1088–108.
39 - van Griensven A, Meixner T, Grunwald S, Bishop T, Diluzio A, Srinivasan R. A global sensitivity analysis tool for the parameters of multi‐variable catchment models. J Hydrol. 2006 Jun 15;324(1–4):10–23.
40 - Abbaspour KC, Yang J, Maximov I, Siber R, Bogner K, Mieleitner J, et al. Modelling hydrology and water quality in the pre‐ailpine/alpine Thur watershed using SWAT. J Hydrol. 2007 Feb 15;333(2–4):413–30.
41 - Moriasi DN, Arnold JG, Van Liew MW, Bingner RL, Harmel RD, Veith TL. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans ASABE. 2007 May–Jun;50(3):885–900.
42 - Krause P, Boyle DP, Bäse F. Comparison of different efficiency criteria for hydrological model assessment. Adv Geosci. 2005;5:89–97.
43 - Ojstršek Zorčič P. Selection and Implementation of Ecoremediation Measures in the Reservoir Watershed. Ljubljana: University of Ljubljana; 2015.
44 - Glavan M, White S, Holman IP. Evaluation of river water quality simulations at a daily time step—experience with SWAT in the axe catchment, UK. Clean Soil Air Water. 2011 Jan;39(1):43–54.
45 - Van Liew MW, Arnold JG, Garbrecht JD. Hydrologic simulation on agricultural watersheds: choosing between two models. Trans ASAE. 2003 Nov–Dec;46(6):1539–51.
46 - Zorn M, Mikos M. Soil erosion measurements in a forest in Slovenian Istria. Journal of the Association of Surveyors of Slovenia. 2010 Sep; 68(7/8): 361–366, 379–392
47 - Bracmort KS, Arabi M, Frankenberger JR, Engel BA, Arnold JG. Modeling long‐term water quality impact of structural BMPs. Trans ASABE. 2006 Mar–Apr;49(2):367–74.
48 - Cho J, Lowrance RR, Bosch DD, Strickland TC, Her Y, Vellidis G. Effect of watershed subdivision and filter width on SWAT simulation of a coastal plain watershed1. J Am Water Resour Assoc. 2010 Jun;46(3):586–602.
49 - Dabney SM, Moore MT, Locke MA. Integrated management of in‐field, edge‐of‐field, and after‐field buffers. J Am Water Resour Assoc. 2006 Feb;42(1):15–24.
50 - Norris JE. Slope Stability and Erosion Control Ecotechnological Solutions. Dordrecht: Springer; 2008.
51 - Verstraeten G, Poesen J, Gillijns K, Govers G. The use of riparian vegetated filter strips to reduce river sediment loads: an overestimated control measure? Hydrol Process. 2006 Dec 30;20(20):4259–67.
52 - Chiang LC, Chaubey I, Hong NM, Lin YP, Huang T. Implementation of BMP strategies for adaptation to climate change and land use change in a pasture‐dominated watershed. Int J Environ Res Pub Health. 2012 Oct;9(10):3654–84.
53 - Bohanec M. Decisions and modeling. Ljubljana, Slovenia: DMFA ‐ založnistvo; 2012.