Next Article in Journal
Structural and Functional Variations of the Macrobenthic Community of the Adige Basin along the River Continuum
Next Article in Special Issue
Framework to Study the Effects of Climate Change on Vulnerability of Ecosystems and Societies: Case Study of Nitrates in Drinking Water in Southern Finland
Previous Article in Journal
A Case Study of the Effects of Management Interventions on the Phosphorus Dynamics at a Coastal, Eutrophic, Caribbean lagoon (Laguna Cartagena, Puerto Rico)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nutrient Load Mitigation with Wintertime Cover as Estimated by the INCA Model

1
Biodiversity Centre, Finnish Environment Institute, 00790 Helsinki, Finland
2
Natural Resources Institute Finland, Natural Resources, Water Quality Impacts, 31600 Jokioinen, Finland
3
Department of Aquatic Sciences and Assessment, Swedish University of Agricultural Sciences, P.O. Box 7050, 750-07 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
Water 2021, 13(4), 450; https://doi.org/10.3390/w13040450
Submission received: 23 December 2020 / Revised: 28 January 2021 / Accepted: 3 February 2021 / Published: 9 February 2021
(This article belongs to the Special Issue Current Trends in Catchment Biogeochemical and Hydrological Modelling)

Abstract

:
Increased nutrient loading causes deterioration of receiving surface waters in areas of intensive agriculture. While nitrate and particulate phosphorus load can be efficiently controlled by reducing tillage frequency and increasing vegetation cover, many field studies have shown simultaneously increased loading of bioavailable phosphorus. In the latest phase of the Rural Programme of EU agri-environmental measures, the highest potential to reduce the nutrient loading to receiving waters were the maximum limits for fertilization of arable crops and retaining plant cover on fields with, e.g., no-till methods and uncultivated nature management fields. Due to the latter two measures, the area of vegetation cover has increased since 1995, suggesting clear effects on nutrient loading in the catchment scale as well. We modeled the effectiveness of agri-environmental measures to reduce phosphorus and nitrogen loads to waters and additionally tested the performance of the dynamic, process-based INCA-P (Integrated Nutrients in Catchments—Phosphorus) model to simulate P dynamics in an agricultural catchment. We concluded that INCA-P was able to simulate both fast (immediate) and slow (non-immediate) processes that influence P loading from catchments. Based on our model simulations, it was also evident that no-till methods had increased bioavailable P load to receiving waters, even though total P and total N loading were reduced.

1. Introduction

Increased nutrient loading of both phosphorus (P) and nitrogen (N) cause deterioration of receiving surface waters in areas of intensive agriculture. In Finland, agriculture is most intensive and is the dominant land cover on clayey soils situated in the southwestern part of the country from where runoff enters the Baltic Sea. In this region, a major proportion (60–90%) of total phosphorus (TP) load is transported as particulate phosphorus (PP) with eroded soil particles [1], and cultivated fields contribute 66–100% of the total suspended sediments delivered from a small agricultural catchment [2]. Moreover, past and present P excess fertilization is reflected in high soil test p values (STP) of the fields in the area, which is again connected with high dissolved reactive P (DRP) concentrations in runoff from fields [3,4].
Increased P loads, the major nutrient controlling eutrophication in many aquatic systems (e.g., [5]), have been observed to cause eutrophication of surface waters. The Baltic Sea may be seasonally or spatially N limited [6]. Further, P is not always the limiting nutrient in freshwaters. In Finland, some smaller lakes are observed to be partly N limited [7], calling for measures to reduce the loading of both P and N from agriculture.
Nitrogen (N) has partly different sources and pathways than P. Most commonly, it is transported as nitrate (NO3), which is highly leachable and can also percolate into ground water. Nitrate can reach receiving waters because of agricultural activities, such as N fertilization and tillage, enhancing the decay of crop residues and soil organic matter (SOM). Under anaerobic conditions, NO3 may be denitrified to N2.
For cutting nutrient load from agriculture in Finland, the main tool since 1995 has been the Rural Programme of EU, whereby more than 90% of the farmers have joined and have been subsidized for using agri-environmental measures (AEM). In the latest phase of AEM (2014–2020), measures with the highest potential for reducing the nutrient loading on waters were upper limits for P and N fertilization of arable crops, winter surface cover of fields, and green set-aside as uncultivated nature management fields [8,9].
The most common measures to mitigate PP load to receiving waters by preventing erosion and suspended sediment transport include autumn ploughing or no-till that leaves an undisturbed soil surface over winter. For example, after changing over from intensive autumn tillage (ploughing) to permanent vegetation cover on clayey soils, erosion and PP loss decreased by 70–65% and 70–54%, respectively [10,11]. Moreover, N loads often decrease due to less intensive tillage [12].
While PP load can be clearly controlled by soil surface cover over winter, many field studies have shown simultaneously increased loads of DRP [11,13,14,15,16]. One explanation is that the shallower tillage depth results in both stratification of P in the upper soil layers, (e.g., [17,18,19]) and, in Finnish conditions, higher risk of surface runoff [11,20]. The estimated P balance (P fertilization − P uptake in yield) can be 3 kg ha−1 higher in no-till fields than in cultivated fields, due to lower crop yield, which together with plant residues that are left on the soil surface gradually increase soil P status in the top layer and thus DRP load potential from the fields [13]. When tests have shown only around 5% (0–13%) of PP bioavailable in river waters but DRP totally bioavailable [21], their magnitudes and proportions are important to the actual ecological consequences of the different measures in place in an agricultural catchment.
In North America, Lake Erie is an example of severe water quality problems due to harmful algal blooms that have been connected to controversial effects of surface cover. Even though agri-environmental measures (no-till, buffer strips, etc.) have been able to decrease considerably suspended sediment and PP load from catchments to the lake, DRP loads have been increasing. Algal blooms are directly linked with increasing riverine input of DRP due to agricultural management and stratification of P at the soil surface [17,22,23,24].
To estimate the efficiency of agricultural water protection measures, e.g., the success of AEM; water quantity and quality are frequently monitored in small catchments (e.g., [25]). Especially high frequency observations together with statistical tools can improve understanding of the processes that influence water quality [26]; however, that dense data is not always available at most locations, and/or time series are often too short.
Mathematical modeling can be useful to consider the combined effect of different pressures and water quality models are developed to determine the source, transformation, transport, and delivery of nutrients and suspended sediments in catchments and into water bodies. These models can be used to predict water quality in areas where monitoring is not feasible, or to help estimate the changes in water quality due to management actions. Compared to empirical models, process-based models can be more accurate to simulate conditions outside the current observations (e.g., [27]). Catchment-scale models, like INCA (Integrated Nutrients in Catchments) [28] and SWAT (Soil and Water Assessment Tool) [29], are especially valuable in assessing both water quality and quantity.
As catchments are heterogeneous systems and hydro-biogeochemical processes are non-linear, the use of dynamic, process-based models has also raised several questions concerning parameterization and calibration, such as equifinality and over-parameterization (e.g., [30]), and whether grid-based observations are valid in describing heterogeneous catchments [31,32,33,34]. Different sensitivity and uncertainty analyses are primarily concerned with the question of how model outputs are affected by variability in model parameters and input values and provide useful information when these components are not completely known. Sensitivity analysis [35] identifies which parameters predominately control model behavior, whereas uncertainty analysis relates prediction errors to uncertainty in the model structure and parameters. Good modeling practice requires that the modelers provide an evaluation of the confidence in the model, possibly assessing the uncertainties associated with the modeling process and with the outcome of the model itself (e.g., [36,37]).
Our main aim was to evaluate the effectiveness of AEM introduced by the Rural Programme to reduce P and N loads to receiving water courses using a chained set of catchment-scale models. Our motivation to concentrate on fields and agri-environmental measures is based on our previous study, where changes in general land use, animal density, point sources, or climate did not explain changes in water quality. Our second aim was to use sensitivity and uncertainty analysis methods to evaluate the performance of the dynamic, process-based INCA-P model for simulating P dynamics in an agricultural catchment against a large set of observational data at different levels.

2. Materials and Methods

2.1. Study Area

The Yläneenjoki catchment (233 km2) has been one of the follow-up study areas of the effects of AEM (earlier FAEP) [8]. The catchment is located in the upper reaches of the Eurajoki River basin in southwestern Finland, as shown in Figure 1. The Eurajoki River discharges into the Bothnian Sea that is mostly P but intermittently N limited [6]. The mean annual precipitation in the area is 700 mm and the mean annual temperature is +4 °C (data from Finnish Meteorological Institute). With a mean discharge of 1.94 m3 s−1 in 1991–2010 [38], the Yläneenjoki River flows to the Pyhäjärvi Lake, which is an important lake for water supply, commercial fishing, and recreation [39]. Eutrophication of the lake is ongoing. The Yläneenjoki River itself is estimated to achieve good ecological status by 2027 according to the river basin plan of the Water Framework Directive [40].
Main soil types in the catchment are clay (Eutric Cambisol 2, Vertic Cambisol) in the cultivated area and till and rock (Lithic Leptosol 1, Haplic Podzol 1, 2) in the forested area with some organic soil types (Fibric/Terric Histosol 1) [41]. Arable fields cover 28% of the catchment area, and the main crops are spring cereals with some sugar beet (Beta vulgaris L.) cultivation. The plant production is accompanied by animal husbandry, mainly pig and poultry, with animal density of 0.53 AU ha−1 in 2005 [42]. In the catchment, most of the farmers are committed to AEM included in the Rural Programme [42].

2.2. Data Available

2.2.1. Catchment-Scale Data

Meteorological data, including daily measurements of precipitation and temperature, were obtained from the Finnish Meteorological Institute (FMI), Jokioinen Observatory station (WGS84 60.81, 23.5).
Maps of soil types, land cover, and the digital elevation model all had a grid size of 25 × 25 m. Soil data were obtained from [41] and land cover data were obtained from the CORINE 2006 GIS (geographical information system). Land cover data were supplemented by more detailed field parcel data of the Natural Resources Institute Finland.
The farmers were interviewed about their annual cultivation methods since 1995 up to 2005 [42,43,44]. Data between 2008 and 2016 originated from annual statistics of Natural Resources Institute Finland. Along with AEM, there has been an increase in the share of no-till methods and decrease in autumn ploughing as the primary tillage method of the spring cereal fields, as shown in Figure 2, while the share of reduced tillage (e.g., stubble cultivation) has remained under 20% of the field area throughout the study period, and the main changes have happened between no-till and ploughing.
Measurements of river discharge (daily values) and suspended sediment, TP, total dissolved phosphorus (TDP), and DRP concentrations (grab samples: on average 33 per year in the Yläneenjoki River at Vanhakartano station) were available from national monitoring programs [45]. The main observation site (WGS84 60.87, 22.40) is located close to the river outlet, as shown in Figure 1. In the upper reaches there was one automatic measurement station (WGS84 60.88, 22.58), which collected daily data of discharge and turbidity over five periods between 15 March 2006 and 9 November 2007.
In the analysis of TP, polyphosphates and organophosphorus compounds were first converted to orthophosphate by acid peroxodisulphate digestion under pressure at 120 °C. Orthophosphate ions reacted in an acidic solution containing molybdate and antimony ions to form a phosphomolybdate complex. Reduction of the complex with ascorbic acid formed a colored molybdenum blue complex, the absorption of which was measured at 880 nm by a spectrophotometer. TDP was similarly analyzed after first filtering the sample through a Nuclepore polycarbonate membrane (0.4 µm). Phosphate P was analyzed after filtering the sample and used as a measure of DRP. Particulate phosphorus (PP) was calculated as difference of TP and TDP. From TP, around 69% was PP, 19% DRP, and 12% dissolved unreactive phosphorus (estimated as the difference between TDP and DRP). Turbidity was used as an estimate of suspended sediment concentration at the upper measurement location due to their good correlation (y = 0.86 × x + 2.38, r2 = 0.94). TP concentration correlated with turbidity (R2 = 0.68) and TDP with DRP (R2 = 0.78). TN determination was initiated by digestion with peroxodisulphate, followed by reduction of NO3 with a Cd amalgam and determination of NO2 by the azo color method.

2.2.2. Field Scale Data

Field scale data with runoff from 2002–2007 (green set-aside) and 2008–2017 (crop cultivation) were obtained from the Kotkanoja experimental field of Natural Resources Institute situated on a clay soil (Protovertic Luvisol, Cleyic, Cutanic) in Jokioinen. In the study years, spring cereals were cultivated in the field with either continuous no-till (two 0.5 ha plots) or autumn ploughing (two 0.5 ha plots) as primary tillage and winter cover options. The field had a gentle slope of 2% (1–4%) and a well-functioning sub-surface drainage system (the averaged share of sub-surface flow of total measured runoff, sub-surface plus surface runoff, was 63 and 75% in no-till and ploughing, respectively). Quantities of both the surface runoff and the sub-surface drainage flow were recorded throughout the years by tipping buckets and the flow-weighted concentrations of TP, DRP, and TN were measured using similar methods as described above with the exception of using a filter pore size of 0.2 µm for DRP analyses (partly published for the above time interval, [11]).

2.3. Models and Model Set-Ups

2.3.1. RUSLE

We used a modification of the Universal Soil Loss Equation (USLE) [46] to calculate the relative erosion sensitivity for each of the fields in the catchment based on land use, soil type, and topography. In the modified version RUSLE [47], spatial differences in climate are not considered, and the range of the factors has been adjusted for Finnish conditions. The relative erosion sensitivity (K) is calculated as a product of four constant factors, determined by soil type (E), slope (R), land use (M), and distance to water (V) as
K = E × R × M × V

2.3.2. PERSiST

We used a flexible rainfall-runoff modeling toolkit, PERSiST (the Precipitation, Evapotranspiration, and Runoff Simulator for Solute Transport), with the INCA family of models. PERSiST is designed for simulating present-day hydrology, projecting possible future effects of climate or land use change on runoff and catchment water storage [48]. PERSiST has limited data requirements and is calibrated using observed time series of precipitation, air temperature and runoff, at one or more points in a river network.

2.3.3. INCA-N and INCA-P

INCA-N [27,49] and INCA-P [50] are dynamic mass-balance models, which attempt to track the temporal variations in the hydrological flow paths and nutrient transformations and stores, in both the land and in-stream components of a river system. INCA models provide daily and annual land-use specific nutrient loads for all stores and transformation processes, and include concentrations in the soil water, ground water and direct runoff. Decay of the soil organic N pool is described by first-order kinetics, and the decay rate is limited by temperature and moisture. The change in mass of P adsorbed to or desorbed from the soil is calculated using a formulation of the Freundlich isotherm [51], where the change in mass adsorbed is proportional to the difference between TDP concentration in soil water and the equilibrium TDP concentration at which no adsorption or desorption occurs (EPC0,soil). It is a dynamic variable that depends on the mass of P in the labile soil store. The other soil store is the inactive P store. This allows for the simulation of longer term soil P dynamics and fits with monitoring data showing EPC0,soil increasing as soil P enrichment increases (e.g., [52]).

2.3.4. PEST

The Model-Independent Parameter Estimation system (PEST; [53]) was used to make a parameter sensitivity and uncertainty analysis. PEST is a non-linear parameter estimation tool for automatic model calibration and predictive analysis. It aims to minimize the weighted sum of squared differences between the model simulation and an observed set of data. The optimization problem is iteratively solved by linearizing the relationship between model output and parameters using a Taylor series expansion. The parameter vector is updated at each step using the Gauss–Marquardt–Levenberg algorithm [54,55]. The derivatives of the model outputs with respect to its parameters provide a measure of parameter sensitivity. Uncertainty analysis is conducted by keeping the objective function with an aim to minimize or maximize a certain set of observations.

2.4. Model Set-Up

We used the branching version 1.4.11 of the INCA models in this study. To reduce the number of parameters and thus the risk for over parameterization, and to enable multiple model runs, the calibration was simplified as much as possible, while still aiming to keep the characters of the system.
The INCA-P model was calibrated against daily observations of discharge and monthly to biweekly observations of suspended sediments, and TDP and TP concentrations observed close to outlet. There were also available short series of automated discharge and turbidity measurements in the upper sub-catchment. The calibration period was 2003–2008 and the validation period 2009–2011. The INCA-N model, in turn, has previously been calibrated against NO3-N and NH4-N concentrations in the same river [56].
The Yläneenjoki catchment was divided into two parts according to the erosion sensitivity map based on the modified USLE. In the less sensitive part of the study area, soil types are typically not prone to erosion and fields are not directly connected to the river. In the highly sensitive area, soil types are prone to erosion and fields are located on steep slopes and are typically also connected to surface water courses [57]. These two areas have different parameter values for erosion and transport capacity processes. In the Yläneenjoki catchment, 40% of the field area was allocated in the less sensitive area. In addition, the upper reach where the automatic measurement station was located was split to form a separate sub-catchment. Moreover, the river was divided into 3 sections according to channel width.
The average fertilizer use, crop yield, P uptake, and initial soil-test P (STP) were derived from farmer interviews (e.g., [42]). Land use was divided into three classes. The most erosion sensitive class included fields where spring cereals were cultivated on clay soils. Grass fields are less sensitive to erosion due to a longer growing period and denser vegetation cover. Sugar beet was cultivated on coarser soil types than spring cereals, and the fields do not have permanent or long-lasting vegetation cover. In addition, P balances that affect STP tend to be high for the sugar beet fields. P loads from forests were calibrated against observed loads [58] and P equilibrium [59]. River water and sediment equilibrium is based on the literature [60], but no observations of P exchange between sediment and water in the river were available.
In the PEST analysis conducted against suspended sediment and TP concentrations, three groups of parameters were used. Furthermore, the parameter ranges, initial parameter values, and parameter increments must be provided for PEST. Parameters with known values, e.g., catchment area and areal percentages of different crops, belonged to the first group. These values were not allowed to change. The second group included relatively well-defined parameters for which the range of values was set narrow and to have known limits. The third group includes parameters for which the values were not well-defined, and the range was set to either +/−25% or +/−50%.
Kotkanoja research field measurements were used to calibrate sub models and equations in the INCA models. Nutrient concentrations in single land use classes (crop–soil type–tillage method combination) were calibrated against concentrations measured in the research field. Equilibrium coefficients were calibrated against those obtained for different soil profile layers in the research field. In addition, the calibrated equilibrium coefficients were compared to measured equilibrium coefficients in corresponding laboratory trials [61].
We used different goodness-of-fit indices. Pearson’s coefficient of determination (R2) is an index of the degree of linear relationship between observed and simulated data, so that value 0 indicates no linear relationship. The Nash–Sutcliffe efficiency (NSE) is commonly used in hydrology. It is a normalized statistic that determines the relative magnitude of the residual variance compared to the measured data variance [62]. Values lower than 0.0 indicate that the mean observed value is a better predictor than the simulated value. RMSE (Root Mean Square Error) is one of the commonly used error index statistics, and it is commonly accepted that the lower the RMSE, the better the model performance. Percent bias (PBIAS) measures the average tendency of the simulated data to be larger or smaller than observations. The optimal value of PBIAS is 0.0, with low-magnitude values indicating accurate model simulation.

2.5. Scenarios

We created three scenarios that were all run for a 22-year-long period using the same initial conditions and climatic data but different agricultural land uses:
  • Autumn ploughing as the primary tillage for all spring cereal fields in the Yläneenjoki catchment;
  • Primary tillage as actually performed in the catchment in 1994–1995;
  • Primary tillage as actually performed in the catchment in 2008–2010.
Scenario 1 was set as the baseline, and the two other tillage scenarios were compared to the baseline.

3. Results

3.1. Calibration and Uncertainty Bounds

Split-sample calibration and validation periods cover years 1995–2007, when exact cultivation practices in the area were available. The starting point for sensitivity and uncertainty analysis was in turn the manually calibrated model for the period 2003–2007, when automatic data was available at the upper reaches. Manual calibration was further finetuned by automatic PEST calibration. The validation period was 1995–2002. Goodness-of-fit values R2, RMSE, percentage bias (PBIAS), and Nash–Sutcliffe efficiency (NSE) for discharge are presented in Table 1.
For phosphorus validation, measures with R2 ≤ 0.40 were considered as not satisfactory, whereas for nitrogen 0.30 < R2 ≤ 0.60 were satisfactory and 0.60 < R2 ≤ 0.70 were good. Values for PBIAS ≥ |30| were considered not satisfactory, |20| ≤ PBIAS < |30| satisfactory, |15| ≤ PBIAS < |20| good, and PBIAS < |15| very good in their accuracy according to published guidelines [63,64].
Discharge calibration was good based on NSE. Nitrate calibration was satisfactory, and TP and DRP calibration were not satisfactory. During the validation period, goodness-of-fit values of discharge and nutrients decreased considerably along the years. On the other hand, PBIAS of discharge was <10 during both calibration and validation periods, which resulted in “very good” calibration. In addition, NO3 simulations were very good, and TP and DRP simulations satisfactory. Simulated TP and DRP concentrations reached the same level and seasonality as the observations, as shown in Figure 3 and Figure 4. There was greater variation in the observed than in the simulated concentrations, but that variation fell between 95% uncertainty bounds, as shown in Figure 5.

3.2. Plot and Field Measurements

In the INCA-P soil compartment, the labile P pool corresponds to the easily soluble P pool from which soluble P leaching occurs. Soil-test P values can thus be used as indicators of the size of that pool. The simulated influence of soil P balance (the difference between the input of P in fertilizers and output in crop uptake) on the labile P pool was compared with the results of an empirical model [61], which related the P balance surplus (or deficit) in a farm to the edge-of-field loads of algal-available P. Based on long-term fertilizer trials, the model first estimated the change in soil-test P of the topsoil with the aid of the soil-surface P balance. Soil-test P was then used to approximate the concentration of DRP in surface runoff and drainage flow, as adjusted for different P application types. The INCA-P model produced slightly faster change in soil labile P, as indicated by STP, than the empirical model. Even though in the Yläneenjoki system the simulated DRP concentration was lower, the slope of the change was the same, as shown in Figure 6.
Simulated soil water TP, DRP, and NO3 loads corresponded to the calculated loads based on measured drainage water and surface runoff concentrations, as shown in Figure 6, Figure 7 and Figure 8. Simulated loads in the different tillage options of the INCA-P soil TP and DRP was similar to observed surface runoff concentrations from the Kotkanoja field. Simulated nitrate concentrations reflected sub-surface drain concentrations.
Simulated equilibrium isotherms for P were close to what [61] measured in agricultural soils (Figure 9), as shown in Figure 7, and in forest soils, the isotherm values were close to what [59] measured in Podzol soils in Myrtillus type forests. Further, equilibrium concentrations of P in river bottom sediments and in water were in line with measurements of [60]. Cross stands for mean and dots for outliers.

3.3. Parameter Sensitivity

The most sensitive parameters were those describing flow velocity in the river (flow), and catchment and soil sensitivity to erosion, as shown in Figure 10. The parameters describing fertilization, vegetation growth, or P processes in the soil were not among the 20 most sensitive parameters.

3.4. Scenario Results

According to the scenarios, change in tillage methods from autumn ploughing towards no-till methods decrease nitrate and suspended sediment loading from the Yläneenjoki catchment but increase DRP loading compared to the baseline with autumn ploughing only, as shown in Table 2. Similar change is seen in land use specific loads, as shown in Table 3, which are in the range of literature values for the corresponding cultivation methods [12].

4. Discussion

4.1. Model Development

We used a set of goodness-of-fit measures to evaluate the success of the model performance, as they measure different features in the fit between simulated and observed values. Pearson’s coefficient of determination (R2) has been widely used for model evaluation; however, it is oversensitive to high extreme values (outliers) and insensitive to additive and proportional differences between model predictions and measured data [65]. NSE is very commonly used in evaluating the fit of a hydrograph, which provides extensive information on reported values. PBIAS has the ability to clearly indicate poor model performance [66]. Even though R2 values in our study were low, BIAS gave relatively good fit between simulated and observed values. That was confirmed visually in a figure, where observed values lay between simulated uncertainty bounds.
Catchment models offer useful methods for exploring problems in space and time, but the quantity and quality of data are not sufficient to parameterize the models uniquely and the modeling addresses a wide range of uncertainties [34]. In addition, several empirical studies have shown that many models and many parameter combinations give equally good fits to data, indicating that it is impossible to find an optimal model or an optimal parameter set in hydrological modeling—a problem termed equifinality [67]. Thus, the use of experimental knowledge of nutrient processes, soft data, is recommended during calibration as it helps to constrain the model based on the available knowledge of the catchment processes [68]. In this study, we used data from experimental fields to calibrate P processes and their outcomes to observed levels.
Various methods from expert assessment to sensitivity analysis can be applied to evaluate uncertainty of deterministic model outputs, concluding that the best method for uncertainty evaluation is determined by the definition of a model, and the amount of available information [69]. A sensitivity analysis framework was recommended as a good starting point to consider uncertainty of spatially distributed environmental models in general [70]. In our study, the model outcome was most sensitive to parameters that defined river flow and catchment erosion processes. These are the processes that give the immediate response to simulated discharge and TP concentrations. The values of these parameters are easy to either measure or calibrate against the measurements. The model was not sensitive to parameters defining DRP processes, reflecting that these processes do not have immediate effect on P concentrations. Thus, change in environmental conditions/changes that influence these processes are difficult to catch in short-term simulations, or see in calibration.
Most of the P is transported to receiving waters attached to erosion material, and thus most of the P transport models are concentrated on describing these processes accurately. When modeling the fate of DRP, equilibrium isotherms both describe the dynamic, interactive P transport between the soil labile P pool and soil water, and can be defined in the laboratory for different land cover types. As equilibrium equations add complexity to modeling process, most of the models use more simple descriptions, either by using empirical equations that are based on soil physical properties, or by limiting the active land cover types (e.g., [30,71,72]). The interpretation of the results should thus be proportionate to how well the model describes the relevant processes.
The Yläneenjoki catchment is an erosion-sensitive area due to the high elevation differences and erosion-prone clay soil types [73]. As such, it is a typical example of an area where TP and suspended sediment load correlate. Even though field percentage is not exceptionally high, around 70% of TP load comes from fields [74]. The share of forests is presumably reflected in water quality, so that TDP concentration is higher than DRP concentration (assumedly P in organic matter). There have not been major changes or trends in anthropogenic pressures on forests, so we assume that changes in water quality are related to agriculture. In the previous study [57], general changes in land use, number of domestic animals, or climate did not explain changes in water quality.
We concluded that INCA-P was able to describe both fast, immediate processes but also slow, non-immediate processes that influence P loading from catchments. We recommend paying special attention when calibrating slow processes controlling P pools in the soil, because they are easily missed when focusing on fast, immediate processes in suspended sediment and PP loading.

4.2. Nutrient Load Mitigation by Winter Cover

Agri-environmental measures (AEM) of the Rural Programme have shown success in improving nutrient balances by introducing maximum fertilization levels, and in increasing the area of reduced tillage and wintertime vegetation cover [8]. However, while previous studies have shown that lower fertilization has the potential to reduce nutrient loads (e.g., [75,76]), the observed decrease in riverine nutrient loads has been modest [74]. In controlling the phosphorus loads, AEM have so far succeeded only in catchments with erosion-sensitive soils where the decrease in suspended sediment and PP loads exceeded the increase in DRP load [56].
According to this study in the Yläneenjoki catchment, the percentage increase in simulated algal available DRP load was higher than the percentage decrease in TP load when the land use in 2008–2010 statistics was compared to the baseline (autumn ploughing in spring cereal fields). Thus, according to [20], this change in tillage was enough to increase bioavailable P load to the Pyhäjärvi Lake up to 10%. No-till methods were relatively popular already in 1994, before launching AEM, and that change was enough to increase bioavailable P load by 0–2%. Statistics from different years are not directly comparable as they originate from different sources and have different resolutions. In any case, we can conclude, based on our model simulations, that no-till methods probably have a contribution to increased bioavailable P load to the Pyhäjärvi Lake.
This study reflects the contradictory nature of vegetation cover and especially no-till methods. When planning catchment-scale water protection methods, the tradeoff is between N, P, and suspended sediment loads, at least when there is no such favorable system change in, e.g., soil structure due to the introduced methods that would clearly influence the result in the long term, like the reduced proportion of surface runoff or the lower need of P fertilizer input for good crop growth. In our research field data [11], the immediate increasing effect induced by green set-aside (that is very similar to the management of the uncultivated nature management field of AEM) on DRP loss was clearly less, while the decreasing effect on TP loss was even more than for no-till.
If the focus is on N-limited parts of the Baltic Sea, no-till methods are one possible measure to reduce N loading. If the focus is on P-limited lakes or sea areas, other methods to decrease DRP loads should also be considered. Even though nutrient balances in Finland have decreased [8], input of P in mineral fertilizers and manure is still above the need of plants, if soil reserves are accounted for [77], and P fertilization needs to be further optimized to speed up the decline of STP values on fields where they are unnecessary high for good crop growth. By this means, it is possible to achieve a permanent reduction in the P loading potential of fields. Moreover, some researchers (e.g., [17]) have even recommended that no-till fields should be occasionally tilled or ploughed conventionally to cut down the stratification of P in the top layers. Obviously, the controversial effects of wintertime cover on nutrient loading to surface waters call for advanced decision tools for site-specific and optimized soil management planning.

5. Conclusions

An equilibrium isotherm describes the dynamic, interactive P transport between the soil labile P pool and soil water in the INCA-P model. We calibrated successfully these equations against a long-term field trial. The model outcome was still the most sensitive to parameters that defined river flow and catchment erosion processes. These are the processes that give the immediate response to simulated discharge and TP concentrations. We concluded that INCA-P was able to describe both fast, immediate processes but also slow, non-immediate processes that influence P loading from catchments.
According to this study in the Yläneenjoki catchment, the percentage increase in simulated algal available DRP load was higher than the percentage decrease in TP load, when the land use in 2008–2010 statistics was compared to the baseline (autumn ploughing in spring cereal fields). This change in tillage was enough to increase bioavailable P load to the Pyhäjärvi Lake up to 10%. We can conclude based on our model simulations, that no-till methods probably have a contribution to increased bioavailable P load to the Pyhäjärvi Lake.

Author Contributions

Conceptualization, K.R. and E.T.; methodology, K.R., E.T., and M.F.; data, R.L. and J.E.C.B.; formal analysis, K.R., J.E.C.B., and R.L.; writing, K.R., E.T., and M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Agriculture and forestry (Project MYTTEHO), Nordic Council of Ministers (Machu Picchu, Biowater), and EU FP7 (REFRESH).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available discharge and river water quality datasets were analyzed in this study. This data can be found here: https://www.syke.fi/avointieto (accessed on 5 September 2018). Statistical data is available at https://www.luke.fi/avoin-tieto/tilastopalvelu/ (accessed on 5 September 2018). The field scale data presented in this study are available on request from the corresponding author. The data are not publicly available due to content of personal information.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pietiläinen, O.-P.; Rekolainen, S. Dissolved reactive and total phosphorus load from agricultural and forested basins to surface waters in Finland. Aqua Fenn. 1991, 21, 127–136. [Google Scholar]
  2. Pietiläinen, O.-P.; Ekholm, P. Origin of eroded material in a small agricultural drainage basin in southwestern Finland. Aqua Fenn. 1992, 22, 105–110. [Google Scholar]
  3. Uusitalo, R.; Aura, E. A rainfall simulation study on the relationships between soil test P versus dissolved and potentially bioavailable particulate phosphorus forms in runoff. Agric. Food Sci. 2005, 14, 335–345. [Google Scholar] [CrossRef] [Green Version]
  4. Sharpley, A.N. Soil phosphorus dynamics: Agronomic and environmental impacts. Ecol. Eng. 1995, 5, 261–279. [Google Scholar] [CrossRef]
  5. Sharpley, A.N.; Rekolainen, S.; Tunney, H. Phosphorus in Agriculture and Its Environmental Implications, in Phosphorus Loss from Soil to Water; Tunney, H., Carton, O.T., Brookes, P.C., Johnston, A.E., Eds.; CAB International Press: Cambridge, UK, 1997; pp. 1–54. [Google Scholar]
  6. Tamminen, T.; Andersen, T. Seasonal phytoplankton nutrient limitation patterns as revealed by bioassays over Baltic Sea gradients of salinity and eutrophication. Mar. Ecol. Prog. Ser. 2007, 340, 121–138. [Google Scholar] [CrossRef] [Green Version]
  7. Pietiläinen, O.-P.; Räike, A. Typpi ja Fosfori Suomen Sisävesien Minimiravinteena. In Suomen Ympäristö, Ympäristönsuojelu; Finnish Environment Institute: Helsinki, Finland, 1999; Volume 313, 64p. [Google Scholar]
  8. Aakkula, J.; Kuussaari, M.; Rankinen, K.; Ekholm, P.; Heliölä, J.; Hyvönen, T.; Kitti, T.; Salo, T. Follow-Up Study on the Impacts of Agri-Environmental Measures in Finland. In Evaluation of Agri-Environmental Policies Selected Methodological Issues and Case Studies; OECD Publishing: Paris, France, 2012; pp. 111–127. [Google Scholar]
  9. Hyvönen, T.; Heliölä, J.; Koikkalainen, K.; Kuussaari, M.; Lemola, R.; Miettinen, A.; Rankinen, K.; Regina, K.; Turtola, E. Maatalouden Ympäristötoimenpiteiden Ympäristö- ja Kustannustehokkuus (MYTTEHO): Loppuraportti. In Luonnonvara- ja Biotalouden Tutkimus 12/2020; Luonnonvarakeskus: Helsinki, Finland, 2020. [Google Scholar]
  10. Puustinen, M.; Koskiaho, J.; Peltonen, K. Influence of cultivation methods on suspended solids and phosphorus concen-trations in surface runoff on clayey sloped fields in boreal climate. Agric. Ecosyst. Environ. 2005, 105, 565–579. [Google Scholar] [CrossRef]
  11. Uusitalo, R.; Lemola, R.; Turtola, E. Surface and Subsurface Phosphorus Discharge from a Clay Soil in a Nine-Year Study Comparing No-Till and Plowing. J. Environ. Qual. 2018, 47, 1478–1486. [Google Scholar] [CrossRef] [Green Version]
  12. Puustinen, M.; Turtola, E.; Kukkonen, M.; Koskiaho, J.; Linjama, J.; Niinioja, R.; Tattari, S. VIHMA—A tool for allocation of measures to control erosion and nutrient loading from Finnish agricultural catchments. Agric. Ecosyst. Environ. 2010, 138, 306–317. [Google Scholar] [CrossRef]
  13. Muukkonen, P.; Hartikainen, H.; Lahti, K.; Särkelä, A.; Puustinen, M.; Alakukku, L. Influence of no-tillage on the distribution and lability of phosphorus in Finnish clay soils. Agric. Ecosyst. Environ. 2007, 120, 299–306. [Google Scholar] [CrossRef]
  14. Ulen, B.; Aronsson, H.; Bechmann, M.; Krogstad, T.; Oygarden, L.; Stenberg, M. Soil tillage methods to control phosphorus loss and potential side-effects: A Scandinavian review. Soil Use Manag. 2010, 26, 94–107. [Google Scholar] [CrossRef]
  15. Uusi-Kämppä, J.; Jauhiainen, L. Long-term monitoring of buffer zone efficiency under different cultivation techniques in boreal conditions. Agric. Ecosyst. Environ. 2010, 137, 75–85. [Google Scholar] [CrossRef]
  16. Uusitalo, R.; Turtola, E.; Lemola, R. Phosphorus losses from a subdrained clayey soil as affected by cultivation practices. Agric. Food Sci. 2007, 16, 352–365. [Google Scholar] [CrossRef] [Green Version]
  17. Baker, D.B.; Johnson, L.T.; Confesor, R.B.; Crumrine, J.P. Vertical Stratification of Soil Phosphorus as a Concern for Dissolved Phosphorus Runoff in the Lake Erie Basin. J. Environ. Qual. 2017, 16, 1287–1295. [Google Scholar] [CrossRef] [Green Version]
  18. Thompson, C.A.; Whitney, D.A. Effects of 30 years of cropping and tillage systems on surface soil test changes. Commun. Soil Sci. Plant Anal. 2000, 31, 241–257. [Google Scholar] [CrossRef]
  19. Cook, R.L.; Trlica, A. Tillage and Fertilizer Effects on Crop Yield and Soil Properties over 45 Years in Southern Illinois. Agron. J. 2016, 108, 415–426. [Google Scholar] [CrossRef] [Green Version]
  20. Turtola, E.; Alakukku, L.; Uusitalo, R.; Kaseva, A. Surface runoff, subsurface drainflow and soil erosion as affected by tillage in a clayey Finnish soil. Agric. Food Sci. 2007, 16, 332–351. [Google Scholar] [CrossRef] [Green Version]
  21. Ekholm, P. Bioavailability of phosphorus in agriculturally loaded rivers in southern Finland. Hydrobiologia 1994, 287, 179–194. [Google Scholar] [CrossRef]
  22. Daloğlu, I.; Cho, K.H.; Scavia, D. Evaluating Causes of Trends in Long-Term Dissolved Reactive Phosphorus Loads to Lake Erie. Environ. Sci. Technol. 2012, 46, 10650–10666. [Google Scholar] [CrossRef]
  23. Jarvie, H.P.; Johnson, L.T.; Sharpley, A.N.; Smith, D.R.; Baker, D.B.; Bruulsema, T.W.; Confesor, R. Increased Soluble Phosphorus Loads to Lake Erie: Unintended Consequences of Conservation Practices? J. Environ. Qual. 2017, 46, 123–132. [Google Scholar] [CrossRef] [Green Version]
  24. Sekaluvu, L.; Zhang, L.; Gitau, M. Evaluation of constraints to water quality improvements in the Western Lake Erie Basin. J. Environ. Manag. 2018, 205, 85–98. [Google Scholar] [CrossRef]
  25. Bol, R.; Gruau, G.; Mellander, P.-E.; Dupas, R.; Bechmann, M.; Skarbøvik, E.; Bieroza, M.; Djodjic, F.; Glendell, M.; Jordan, P.; et al. Challenges of Reducing Phosphorus Based Water Eutrophication in the Agricultural Landscapes of Northwest Europe. Front. Mar. Sci. 2018, 5, 5. [Google Scholar] [CrossRef] [Green Version]
  26. Rode, M.; Wade, A.; Cohen, M.J.; Hensley, R.T.; Bowes, M.J.; Kirchner, J.W.; Arhonditsis, G.B.; Jordan, P.; Kronvang, B.; Halliday, S.J.; et al. Sensors in the Stream: The High-Frequency Wave of the Present. Environ. Sci. Technol. 2016, 50, 10297–10307. [Google Scholar] [CrossRef] [Green Version]
  27. Leavesley, G.H. Modeling the effects of climate change on water resources—A review. Clim. Chang. 1994, 28, 159–177. [Google Scholar] [CrossRef]
  28. Whitehead, P.G.; Wilson, E.J.; Butterfield, D. A semi-distributed Integrated Nitrogen model for multiple source assessment in Catchments (INCA): Part I-model structure and process equations. Sci. Total Environ. 1998, 210–211, 547–558. [Google Scholar] [CrossRef]
  29. Arnold, J.G.; Srinivasin, R.; Muttiah, R.S.; Williams, J.R. Large Area Hydrologic Modeling and Assessment: Part I. Model Development. JAWRA J. Am. Water Resour. Assoc. 1998, 34, 73–89. [Google Scholar] [CrossRef]
  30. Beven, K. Towards an alternative blueprint for a physically based digitally simulated hydrologic response modelling system. Hydrol. Process. 2002, 16, 189–206. [Google Scholar] [CrossRef] [Green Version]
  31. Beven, K. Changing ideas in hydrology—The case of physically-based models. J. Hydrol. 1989, 105, 157–172. [Google Scholar] [CrossRef]
  32. Beven, K. How far can we go in distributed hydrological modelling? Hydrol. Earth Syst. Sci. 2001, 5, 1–12. [Google Scholar] [CrossRef]
  33. Blöschl, G.; Grayson, R. Advances in distributed hydrological modelling-towards a new paradigm. In Water Resources and Environment Research ICWRER 2002, Proceedings of the 3rd International Conference on Water Resources and Environment Research, Dresden, Germany, 22–25 July 2002; Dresden University of Technology: Dresden, Germany, 2002. [Google Scholar]
  34. Guillaume, J.H.A.; Jakeman, J.D.; Marsili-Libelli, S.; Asher, M.; Brunner, P.; Croke, B.F.W.; Hill, M.C.; Jakeman, A.J.; Keesman, K.J.; Razavi, S.; et al. Introductory overview of identifiability analysis: A guide to evaluating whether you have the right type of data for your modeling purpose. Environ. Model. Softw. 2019, 119, 418–432. [Google Scholar] [CrossRef]
  35. Hamby, D.M. A review of techniques for parameter sensitivity analysis of environmental models. Environ. Monit. Assess. 1994, 32, 135–154. [Google Scholar] [CrossRef]
  36. Refsgaard, J.C.; Henriksen, H.J. Modelling guidelines—Terminology and guiding principles. Adv. Water Resour. 2004, 27, 71–82. [Google Scholar] [CrossRef]
  37. Refsgaard, J.C.; Henriksen, H.J.; Harrar, W.G.; Scholten, H.; Kassahun, A. Quality assurance in model based water man-agement—review of existing practice and outline of new approaches. Environ. Model. Softw. 2005, 20, 1201–1215. [Google Scholar] [CrossRef]
  38. Korhonen, J.; Haavanlammi, E. Hydrologinen Vuosikirja 2006–2012 (Hydrological Yearbook 2006–2010). In Suomen Ympäristö 8/2012; Luonnonvarat: Helsinki, Finland, 2012; 234p, ISBN 978-952-11-3988-8. [Google Scholar]
  39. Ventelä, A.-M.; Tarvainen, M.; Helminen, H.; Sarvala, J. Long-term management of Pyhäjärvi (southwest Finland): Euthrophication, restoration-recovery? J. Lake Reserv. Manag. 2007, 23, 428–438. [Google Scholar] [CrossRef] [Green Version]
  40. Water Framework Directive (WFD). Directive 2000/60/EC of the European Parliament and the Council of 23 October 2000 Establishing a Framework for Community Action in the Field of Water Policy. Off. J. Eur. Community 2000, L327, 1–73.
  41. Lilja, H.; Uusitalo, R.; Yli-Halla, M.; Nevalainen, R.; Väänänen, T.; Tamminen, P. Suomen Maannostietokanta: Maannoskartta 1:250,000 ja Maaperän Ominaisuuksia. In MTT:n Selvityksiä 114; MTT: Jokioinen, Finland, 2006; 70p. [Google Scholar]
  42. Mattila, P.; Rankinen, K.; Grönroos, J.; Siimes, K.; Karhu, E.; Laitinen, P.; Granlund, K.; Ekholm, P.; Antikainen, R. Viljelytoimenpiteet ja Vesistökuormitus Ympäristötukitiloilla Vuosina 2003–2005 (Changes in Cultivation Practices and Nutrient Loading to the Waters Due to the Agri-Environmental Support Scheme in 2003–2005, in Finnish with English Abstract). In Suomen Ympäristö 40; Suomen Ympäristökeskus: Helsinki, Finland, 2007; 101p. [Google Scholar]
  43. Pyykkönen, S.; Grönroos, J.; Rankinen, K.; Laitinen, P.; Karhu, E.; Granlund, K. Ympäristötuen Mukaiset Viljelytoimenpiteet ja Niiden Vaikutukset Vesistökuormitukseen Vuosina 2000–2002; Suomen Ympäristökeskus: Helsinki, Finland, 2004; 119p. [Google Scholar]
  44. Grönroos, J.; Rekolainen, S.; Palva, R.; Granlund, K.; Bärlund, I.; Nikander, A.; Laine, Y. Maatalouden Ympäristötuki: Toimenpiteiden Toteutuminen ja Vaikutukset vuosina 1995–1997; Suomen Ympäristö, Ympäristönsuojelu: Helsinki, Finland, 1998; Volume 239, 77p. [Google Scholar]
  45. Niemi, J.; Heinonen, P.; Mitikka, S.; Vuoristo, H.; Pietiläinen, O.-P.; Puupponen, M.; Rönkä, E. The Finnish Eurowaternet with Information about Finnish Water Resources and Monitoring Strategies; The Finnish Environment, Environmental Protection 445; Finnish Environment Institute: Helsinki, Finland, 2001; 62p. [Google Scholar]
  46. Renard, K.G.; Foster, G.R.; Yoder, D.C.; McCool, D.K. RUSLE revisited: Status, questions, answers, and the future. J. Soil Water Conserv. 1994, 49, 213–220. [Google Scholar]
  47. Räsänen, A. Kiintoaineen ja Kasviravinteiden Vesistökuormituksen Riskialuekartoitus Aurajoen Valuma-Alueella. In Bio- ja Ympäristötieteiden Laitos; Jyväskylän Yliopisto: Jyväskylä, Helsinki, Finland, 2010. [Google Scholar]
  48. Futter, M.N.; Erlandsson, M.A.; Butterfield, D.; Whitehead, P.G.; Oni, S.K.; Wade, A.J. PERSiST: A flexible rainfall-runoff modelling toolkit for use with the INCA family of models. Hydrol. Earth Syst. Sci. 2014, 18, 855–873. [Google Scholar] [CrossRef] [Green Version]
  49. Wade, A.; Durand, P.; Beaujoan, V.; Wessels, W.; Raat, K.; Whitehead, P.G.; Butterfield, D.; Rankinen, K.; Lepistö, A. Towards a generic nitrogen model of European ecosystems: New model structure and equations. Hydrol. Earth Syst. Sci. 2002, 6, 559–582. [Google Scholar] [CrossRef]
  50. Jackson-Blake, L.A.; Wade, A.J.; Futter, M.N.; Butterfield, D.; Couture, R.-M.; Cox, B.A.; Crossman, J.; Ekholm, P.; Halliday, S.J.; Jin, L.; et al. The INtegrated CAtchment model of phosphorus dynamics (INCA-P): Description and demonstration of new model structure and equations. Environ. Model. Softw. 2016, 83, 1–31. [Google Scholar] [CrossRef] [Green Version]
  51. House, W.A.; Denison, F.H. Factors influencing the measurement of equilibrium phosphate concentrations in river sediments. Water Res. 2000, 34, 1187–1200. [Google Scholar] [CrossRef]
  52. Lin, C.; Banin, A. Effect of Long-Term Effluent Recharge on Phosphate Sorption by Soils in a Wastewater Reclamation Plant. Water Air Soil Pollut. 2005, 164, 257–273. [Google Scholar] [CrossRef]
  53. Doherty, J.; Johnston, J.M. Methodologies for calibration and predictive analysis of a watershed model. JAWRA J. Am. Water Resour. Assoc. 2003, 39, 251–265. [Google Scholar] [CrossRef]
  54. Marquardt, D.W. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  55. Levenberg, K. A method for the solution of certain non-linear problems in least squares. Q. Appl. Math. 1944, 2, 164–168. [Google Scholar] [CrossRef] [Green Version]
  56. Rankinen, K.; Granlund, K.; Etheridge, R.; Seuri, P. Valuation of nitrogen retention as an ecosystem service on a catchment scale. Hydrol. Res. 2014, 45, 411–424. [Google Scholar] [CrossRef]
  57. Rankinen, K.; Gao, G.; Granlund, K.; Grönroos, J.; Vesikko, L. Comparison of impacts of human activities and climate change on water quantity and quality in Finnish agricultural catchments. Landsc. Ecol. 2015, 30, 415–428. [Google Scholar] [CrossRef] [Green Version]
  58. Finér, L.; Lepistö, A.; Karlsson, K.; Räike, A.; Tattari, S.; Huttunen, M.; Härkönen, L.; Joensuu, S.; Kortelainen, P.; Mattsson, T.; et al. Metsistä ja Soilta Tuleva Vesistökuormitus 2020. In Valtioneuvoston Selvitys- ja Tutkimustoiminnan Julkaisusarja; Valtioneuvoston Kanslia: Helsinki, Finland, 2020. [Google Scholar]
  59. Väänänen, R. Phosphorus Retention in Forest Soils and the Functioning of Buffer Zones Used in Forestry. Ph.D. Thesis, University of Helsinki, Helsinki, Finland, 2008. Dissertationes Forestales 60. [Google Scholar]
  60. Koski-Vähälä, J. Role of Resuspension and Silicate in Internal Phosphorus Loading. Ph.D. Thesis, University of Helsinki, Helsinki, Finland, 2001. [Google Scholar]
  61. Peltovuori, T. Phosphorus in Agricultural Soils of Finland—Characterization of Reserves and Retention in Mineral Soil Profiles. Ph.D. Thesis, University of Helsinki, Helsinki, Finland, 2006. [Google Scholar]
  62. Nash, J.E.; Sutcliffe, J.V. River Flow Forecasting through Conceptional Models. 1. A Discussion of Principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  63. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmet, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  64. Moriasi, D.N.; Gitau, M.W.; Pai, N.; Daggupati, P. Hydrologic and water quality models: Performance measures and evaluation criteria. Trans. Am. Soc. Agric. Biol. Eng. 2015, 58, 1763–1785. [Google Scholar]
  65. Le Gates, D.R.; McCabe, G.J., Jr. Evaluating the use of “goodness-of-fit” Measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 1999, 35, 233–241. [Google Scholar] [CrossRef]
  66. Gupta, H.V.; Sorooshian, S.; Yapo, P.O. Status of Automatic Calibration for Hydrologic Models: Comparison with Multi-level Expert Calibration. J. Hydrol. Eng. 1999, 4, 135–143. [Google Scholar] [CrossRef]
  67. Beven, K. A manifesto for the equifinality thesis. J. Hydrol. 2006, 320, 18–36. [Google Scholar] [CrossRef] [Green Version]
  68. Rankinen, K.; Karvonen, T.; Butterfield, D. An application of the GLUE methodology for estimating the parameters of the INCA-N model. Sci. Total. Environ. 2006, 365, 123–139. [Google Scholar] [CrossRef]
  69. Uusitalo, L.; Lehikoinen, A.; Helle, I.; Myrberg, K. An overview of methods to evaluate uncertainty of deterministic models in decision support. Environ. Model. Softw. 2015, 63, 24–31. [Google Scholar] [CrossRef] [Green Version]
  70. Koo, H.; Chen, M.; Jakeman, A.J.; Zhang, F. A global sensitivity analysis approach for identifying critical sources of uncertainty in non-identifiable, spatially distributed environmental models: A holistic analysis applied to SWAT for input datasets and model parameters. Environ. Model. Softw. 2020, 127, 104676. [Google Scholar] [CrossRef]
  71. Jackson-Blake, L.; Sample, J.E.; Wade, A.J.; Helliwell, R.C.; Skeffington, R.A. Are our dynamic water quality models too complex? A comparison of a new parsimonious phosphorus model, SimplyP, and INCA-P. Water Resour. Res. 2017, 53, 5382–5399. [Google Scholar] [CrossRef] [Green Version]
  72. Larsson, M.H.; Persson, K.; Ulén, B.; Lindsjö, A.; Jarvis, N.J. A dual porosity model to quantify phosphorus losses from macroporous soils. Ecol. Model. 2007, 205, 123–134. [Google Scholar] [CrossRef]
  73. Aura, E.; Saarela, K.; Räty, M. Savimaiden Eroosio. In MTT:n Selvityksiä 118; MTT: Jokioinen, Finland, 2006; 32p, Erratum in Macroporous Soils. Ecol. Model. 2007, 205, 123–134. [Google Scholar]
  74. Rankinen, K.; Keinänen, H.; Bernal, J.E.C. Influence of climate and land use changes on nutrient fluxes from Finnish rivers to the Baltic Sea. Agric. Ecosyst. Environ. 2016, 216, 100–115. [Google Scholar] [CrossRef]
  75. Valkama, E.; Uusitalo, R.; Ylivainio, K.; Virkajärvi, P.; Turtola, E. Phosphorus fertilization: A meta-analysis of 80 years of research in Finland. Agric. Ecosyst. Environ. 2007, 130, 75–85. [Google Scholar] [CrossRef]
  76. Valkama, E.; Uusitalo, R.; Ylivainio, K.; Virkajärvi, P.; Turtola, E. Nitrogen fertilization of grass leys: Yield production and risk of N leaching. Agric. Ecosyst. Environ. 2016, 230, 341–352. [Google Scholar] [CrossRef]
  77. Ylivainio, K.; Sarvi, M.; Lemola, R.; Uusitalo, R.; Turtola, E. Regional P Stocks in Soil and in Animal Manure as Compared to P Requirement of Plants in Finland. Baltic Forum for Innovative Technologies for Sustainable Manure Management. WP4 Standardisation of Manure Types with Focus on Phosphorus. In Natural Resources and Bioeconomy Studies 62/2015; Natural Resources Institute Finland: Helsinki, Finland, 2015; 32p. [Google Scholar]
Figure 1. Location of the Yläneenjoki catchment and its sub-catchments.
Figure 1. Location of the Yläneenjoki catchment and its sub-catchments.
Water 13 00450 g001
Figure 2. Development of the share of ploughing and no-till/reduced tillage on spring cereal fields of the Yläneenjoki catchment.
Figure 2. Development of the share of ploughing and no-till/reduced tillage on spring cereal fields of the Yläneenjoki catchment.
Water 13 00450 g002
Figure 3. Monthly averages of dissolved reactive P (DRP) concentrations in river water: (a) observed, (b) simulated.
Figure 3. Monthly averages of dissolved reactive P (DRP) concentrations in river water: (a) observed, (b) simulated.
Water 13 00450 g003
Figure 4. Monthly averages of total phosphorus (TP) concentrations in river water: (a) observed, (b) simulated.
Figure 4. Monthly averages of total phosphorus (TP) concentrations in river water: (a) observed, (b) simulated.
Water 13 00450 g004
Figure 5. Observed total phosphorus concentrations, simulated and 95% uncertainty bounds of simulated concentrations by PEST (Model-Independent Parameter Estimation) at Vanhakartano station.
Figure 5. Observed total phosphorus concentrations, simulated and 95% uncertainty bounds of simulated concentrations by PEST (Model-Independent Parameter Estimation) at Vanhakartano station.
Water 13 00450 g005
Figure 6. The effect of P balance on soil-test P (STP, based on the acid ammonium acetate method, pH 4.65, used in Finland, [62]) values (a) and DRP concentration (b) in river water.
Figure 6. The effect of P balance on soil-test P (STP, based on the acid ammonium acetate method, pH 4.65, used in Finland, [62]) values (a) and DRP concentration (b) in river water.
Water 13 00450 g006
Figure 7. Measured and simulated DRP monthly loads from (a) autumn ploughing, (b) no-till, and (c) green set-aside. Lines stand for minimum, first quartile (Q1), median, third quartile (Q3), and maximum values; cross stands for mean; and dots stand for outliers.
Figure 7. Measured and simulated DRP monthly loads from (a) autumn ploughing, (b) no-till, and (c) green set-aside. Lines stand for minimum, first quartile (Q1), median, third quartile (Q3), and maximum values; cross stands for mean; and dots stand for outliers.
Water 13 00450 g007
Figure 8. Measured and simulated NO3-N monthly load from (a) autumn ploughing, (b) no-till, and (c) green set-aside. Lines stand for minimum, first quartile (Q1), median, third quartile (Q3), and maximum values; cross stands for mean; and dots stand for outliers.
Figure 8. Measured and simulated NO3-N monthly load from (a) autumn ploughing, (b) no-till, and (c) green set-aside. Lines stand for minimum, first quartile (Q1), median, third quartile (Q3), and maximum values; cross stands for mean; and dots stand for outliers.
Water 13 00450 g008
Figure 9. Equilibrium isotherms of P in soil as modeled and measured in the Kotkanoja research field at different soil layers.
Figure 9. Equilibrium isotherms of P in soil as modeled and measured in the Kotkanoja research field at different soil layers.
Water 13 00450 g009
Figure 10. Relative sensitivity of the 20 most sensitive model parameters as labeled with the main component in the model which the parameter affects. Parameter groups are blue: flow velocity, gray: catchment erosion, red: land cover erosion, black: constant.
Figure 10. Relative sensitivity of the 20 most sensitive model parameters as labeled with the main component in the model which the parameter affects. Parameter groups are blue: flow velocity, gray: catchment erosion, red: land cover erosion, black: constant.
Water 13 00450 g010
Table 1. Goodness-of-fit values of PERSiST (Precipitation, Evapotranspiration, and Runoff Simulator for Solute Transport) and INCA-P model application.
Table 1. Goodness-of-fit values of PERSiST (Precipitation, Evapotranspiration, and Runoff Simulator for Solute Transport) and INCA-P model application.
Station CalibrationValidation
R2NSERMSEPBIASR2NSERMSEPBIAS
Automatic stationQ0.6230.45678.76951.33---
NO3-N0.262-95.8840.087---
VanhakartanoQ0.7240.71490.0659.680.3740.373128.3541.21
TP0.167-71.02725.420.1-79.1334.68
DRP0.1-94.18919.780.088-121.11621.71
NO3-N0.351-71.1260.300.168-201.4370.81
Table 2. Simulated change in nutrient and suspended sediment load from Yläneenjoki catchment in different land use periods compared to the baseline (with autumn ploughing only as the primary tillage method in spring cereal cultivation).
Table 2. Simulated change in nutrient and suspended sediment load from Yläneenjoki catchment in different land use periods compared to the baseline (with autumn ploughing only as the primary tillage method in spring cereal cultivation).
Land Use PeriodTPDRPSSNO3-N
[%][%][%][%]
1994–1995−4%2%−6%−4%
2008–2010−5%10%−11%−5%
Table 3. Simulated 22-year mean export coefficient values of different cultivation methods in the Yläneenjoki catchment.
Table 3. Simulated 22-year mean export coefficient values of different cultivation methods in the Yläneenjoki catchment.
Land Use ClassPPDRPSSNO3N
[kg/ha][kg/ha][kg/ha][kg/ha]
Spring cereals, autumn ploughing0.80.4489.118.3
Spring cereals, autumn stubble cultivation0.60.5418.015.2
Spring cereals, no-till0.50.7302.212.4
Grass0.40.8248.95.1
Set-aside grass0.30.6207.54.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rankinen, K.; Turtola, E.; Lemola, R.; Futter, M.; Cano Bernal, J.E. Nutrient Load Mitigation with Wintertime Cover as Estimated by the INCA Model. Water 2021, 13, 450. https://doi.org/10.3390/w13040450

AMA Style

Rankinen K, Turtola E, Lemola R, Futter M, Cano Bernal JE. Nutrient Load Mitigation with Wintertime Cover as Estimated by the INCA Model. Water. 2021; 13(4):450. https://doi.org/10.3390/w13040450

Chicago/Turabian Style

Rankinen, Katri, Eila Turtola, Riitta Lemola, Martyn Futter, and José Enrique Cano Bernal. 2021. "Nutrient Load Mitigation with Wintertime Cover as Estimated by the INCA Model" Water 13, no. 4: 450. https://doi.org/10.3390/w13040450

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