Peatland drainage can affect the natural state of hydrological conditions and nutrient loading but is rarely included in catchment-scale models. To understand the gap, we aimed to use the soil and water assessment tool (SWAT) to observe drained peatlands and their properties to predict nutrients and suspended solids (SS) in the peatland-dominated Simojoki catchment. We integrated drainage networks in SWAT to (i) identify drainage parameters; (ii) determine annual loads and mean concentrations of SS, organic phosphorus (Org-P), total phosphorus (TP), organic nitrogen (Org-N), and total nitrogen (TN); (iii) understand spatial variation of nutrients and SS; and (iv) investigate uncertainty ranges for the estimates. The calibrated model showed a 9.6 per cent bias between the simulated flow and the observations with low-to-medium loading variations for the water quality parameters. For Org-N and TN, the highest loading per year was at the downstream outlet, whereas for SS, Org-P, and TP, it was higher at the upstream outlet of the catchment. This approach of representing the drained peatland in SWAT indicated maximum spatial distributed load in the peat soil in the clear-cut area and can be beneficial in future hydrological modelling efforts in identifying the status of nutrients or SS.

  • Integration of peatland drainage in SWAT can predict nutrients and SS estimations.

  • Calibrated parameters resulted in a good–satisfactory status for nutrients and SS.

  • Uncertainty in estimates captured 35–81% of the observations.

  • Spatial and temporal distribution of loading shows the status of key locations.

  • Highest specific loading: Nutrients in clear-cut peat soil, SS in peat extraction.

Soil material that contains at least 30% content by mass of dead organic matter (Joosten & Clarke 2002) or organic matter with less than 20–35% mineral content can be considered peat (Turetsky et al. 2014). Peat deposition occurs when peat formation is higher than peat decay and is generally slow under waterlogged conditions (Holden et al. 2004). Favourable wet conditions lead to the accumulation of organic matter in specific areas, resulting in the formation of peatlands. Of the 3% peatland cover on Earth (Greenup et al. 2000), Finland has more than 30% peatland in its territory. In northern Finland, the largely drained peatlands usually cover more than 20,000 ha of the area and are over 30 cm in depth (Turunen 2008).

Peatlands are hotspots for biodiversity and ecosystem services (Kadlec & Wallace 2008). They are substantial carbon reservoirs in boreal and subarctic regions, constituting one-third of the global soil carbon pool (Yu et al. 2010). Peatland degradation occurs worldwide because of the drainage of peatlands (Ramchunder et al. 2012). In boreal and temperate zones, approximately 15 million hectares of peatland have been drained for forestry since the 1950s. Drainage activity peaked during the 1970s. In the past two decades, almost half of the total peatland area in Finland has been affected by drainage (Paavilainen & Päivänen 1995).

Peatland drainage modifies the water-holding capacity of peat and the amount of water that leaves it (Landry & Rochefort 2012). Subsurface flow dominates the runoff generation mechanism in drained peatlands, mainly because of the lower water table conditions compared with surface flow-dominated undisturbed peatlands (Holden et al. 2006). Previously drained areas may alter the natural state of hydrological conditions and nutrient loading (Paavilainen & Päivänen 1995; Nieminen et al. 2017b; Finér et al. 2021). Different forest management activities in peatlands can also affect nutrient leaching (Kreutzweiser et al. 2008).

Small catchment research on peatlands (Marttila et al. 2018) clearly shows the short-term (Finér et al. 2021) and long-term (Nieminen et al. 2017a) effects of forestry on water quality. Although long-term national monitoring programs are active at river outlets, little information on peatland drainage in larger catchments is available. Several models have been developed, ranging from the conceptual level to describing the processes of water movement in catchments with different soil types, land use, and management (Singh 2018) through the simulation of hydrology and water quality. A few modelling studies have focussed on peatland-dominated catchments to analyse the impacts of forestry (Lepistö & Syri 2001). However, the effect of the drainage network as the primary driver of change has not been accounted for directly. Often, monitoring large catchments fails to pinpoint the influence of specific land use operations. In recent years, there has been an increase in organic nitrogen exports from peatland-dominated catchments (Lepistö et al. 2021) in the Baltic Sea region. These require attention but are rarely included in catchment-scale simulations. This has led to an unclear understanding of how drainage of peatlands affects water quality at the catchment level.

This study was motivated by a lack of knowledge and the need for better modelling approaches to assess spatiotemporal changes in nutrient and suspended solids (SS) status in a peatland-dominated catchment. The soil and water assessment tool (SWAT) has been successfully implemented under various conditions (Molina-Navarro et al. 2017), including in peatlands (Melaku et al. 2020). However, to our knowledge, it has not yet been used in drained peatlands, especially in areas with intensive peatland forestry.

In this study, we aimed to observe the major land uses in peatlands and their specific properties, using the SWAT model to compare nutrient and SS predictions against observations in peatland-dominated catchments. The Simojoki catchment and one of its sub-catchments, Hosionkoski, in northern Finland, were used to assess the SWAT model performance in a peatland-dominated catchment. The specific objectives were to (i) identify the specific parameters of the SWAT model that can represent the effect of drainage at the catchment scale; (ii) determine annual loads and annual mean concentrations of SS, organic phosphorus (Org-P), total phosphorus (TP), organic nitrogen (Org-N), and total nitrogen (TN); (iii) understand how nutrients and SS vary across different geographical units across the catchment; and (iv) investigate the uncertainty ranges for the estimates.

Study area

The Simojoki catchment (3,160 km2) in northern Finland was used as the study area for modelling in SWAT (Figure 1). The unregulated salmon river runs 193 km from its origin in Lake Simojärvi to Bothnian Bay, dropping 176 m into the sea. The major human impact on catchments is from peatland drainage. Drainage operations began around 1950, with intensive drainage peaking during the 1970s. Peatland drainage was conducted throughout the catchment area to support forest production and peat extraction. An average of approximately 0.5% of the total catchment area is cut annually (Rankinen et al. 2006).
Figure 1

(a) The Simojoki catchment and one of its sub-catchments Hosionkoski with streams, flow, and water quality measurement stations, and weather stations; (b) existing drainage network within 1 m DEM based on a database compiled in a previous study (Bhattacharjee et al. 2021a); (c) land use; (d) soil; and (e) slope (%) maps of the catchment.

Figure 1

(a) The Simojoki catchment and one of its sub-catchments Hosionkoski with streams, flow, and water quality measurement stations, and weather stations; (b) existing drainage network within 1 m DEM based on a database compiled in a previous study (Bhattacharjee et al. 2021a); (c) land use; (d) soil; and (e) slope (%) maps of the catchment.

Close modal

Granites and gneisses are the most common bedrock types and are usually covered with till. There was sand over the bedrock near the catchment outlet. Forests on mineral soils (39%) and peatland forests (53%) dominate (Lepistö et al. 2014), whereas agriculture covers only approximately 3% of the catchment area, mainly near the catchment outlet (Perkkiö et al. 1995). The catchment contains 41.5% mixed soil types, whereas peat covers only 32.87%. A combination of a thin peat layer (THP) and mixed soil type (MST) covered 13.61% of the catchment area. In contrast, other combinations of multiple soil layers cover a very small portion of the catchment (<1%). We found that 54.26% of the catchment area had a slope of less than 2.5%, whereas 25.64% had a slope of 2.5–5% (Figure 1).

Model description

SWAT is one of the most globally used hydrological models for simulating the influence of hydrology and land use on nutrient concentration and loading. Catchment-scale modelling in SWAT includes the dynamics of quantitative and qualitative information and has the potential to simulate drained peatlands. However, the implementation of drained peatlands in a semi-distributed SWAT model is challenging. In addition, model conceptualization, parameter specifications, and management activities are aspects of uncertainty that can arise when applying SWAT (Gupta et al. 2005), particularly for peatland and organic soils.

The SWAT uses Manning's equation to define the flow rate and velocity. Water is routed through the river network using the variable storage routing method or the Muskingum River routing method. To compute sediment yield, SWAT uses a modified universal soil loss equation (MUSLE; Williams & Berndt 1977). The sediment yield of a specific catchment depends on surface runoff, peak flow rate, soil erosion, slope length, and crop management factors. Multiple processes are followed in SWAT to estimate the number of different forms of nitrogen (N). NO3-N in SWAT represents the product of the water volume and average concentration in runoff, lateral flow, and percolation. SWAT uses a loading function for Org-N that depends on the concentration of Org-N in the top soil layer, sediment yield, and enrichment ratio (Williams & Hann 1978). As phosphorus (P) is associated with the sediment phase, soluble P estimation depends on the labile P concentration and runoff volume (Knisel 1982).

Model set-up

Integration of drainage network and catchment delineation

The Simojoki catchment contains many ditches (drainage networks and individual drainage lines) in peatlands, particularly in the forest area (Figure 1(b)). For different land uses, the area drained with open ditches varied substantially throughout the catchment (Figure 1). To implement existing drained peatlands in the SWAT model, we:

  • - compiled a ditch database identified from aerial images of the catchment based on a method developed in an earlier study (Bhattacharjee et al. 2021a);

  • - collected a high-resolution digital elevation model (DEM) (1 m × 1 m) from the National Land Survey of Finland (NLSF [Paituli] 2019);

  • - the drainage (ditch) database contains the ditch length in a shapefile. Thus, to integrate the ditch database into the DEM, we burned the ditch database shapefile using the delineation option in ArcSWAT (Version-2012.10_7.24; SWAT2012.exe Revision 681);

  • - the burned DEM was masked with the existing catchment boundary obtained from NLSF, particularly when the land was flat;

  • - created a flow accumulation raster (1,000 ha) to identify potential stream outlets;

  • - streams and outlets were created based on the flow accumulation raster;

  • - verified streams and outlets using the existing database from NLSF [Paituli] (2019); and

  • - delineated the watershed based on the final outlet of Bothnian Bay, which resulted in 111 sub-catchments.

This process essentially sets up the necessary spatial and topographic data for the SWAT model, ensuring that the drainage networks are accurately represented. Integrating a ditch database into the DEM is a critical step for modelling the hydrology of peatlands in the catchment. The simulation period was also considered so that the model could represent the existing ditch condition in the catchment as closely as possible. The verification steps with existing databases add a quality control layer to the process. Therefore, the final streams were generated based on the hydrological network of the drainage lines of the burned DEM.

Hydrologic response unit (HRU) in SWAT

The hydrologic response unit (HRU) is the smallest spatial unit in SWAT and consists of a unique combination of land use, soil type, and the slope band of a sub-catchment based on user-defined thresholds. Detailed spatial and temporal land use analysis of the catchment was provided in an earlier study (Bhattacharjee et al. 2021b). Based on land use analysis from Landsat (20 m × 20 m) (Bhattacharjee et al. 2021b), six classes were considered from 2015 to represent the catchment (NLSF [Paituli] 2019; Figure 1(c)). Soil data (25 m × 25 m) were collected from the Geological Survey of Finland [GTK] (2023). Ten soil types dominated the catchment (Figure 1(d)) and based on earlier studies (Menberu et al. 2021), we matched them to soils already in the SWAT soil database to be integrated with the model. Five slope bands (Figure 1(e)) were generated using the DEM. These slope bands were the most representative elevation variations in the catchment.

We used threshold values of 15, 20, and 20% for the land use, soil type, and slope, respectively, to define the final HRUs. These threshold values indicate that any land use, soil type, or slope band representing a percentage below these threshold values in each sub-catchment was deleted to simplify the model, subsequently proportionally increasing the remaining land use, soil type, and slope band. As mentioned in an earlier study, commonly used thresholds might lead to substantial information loss in catchment landscapes of the catchment (Her et al. 2015). Thus, the thresholds were chosen carefully, resulting in 634 HRUs representing the drained peatland areas within each HRU and sub-catchment.

Weather, river flow, and water quality data

We used open-access gridded weather data (10 km × 10 km) provided by the Finnish Meteorological Institute (FMI) from the Paituli server (NLSF [Paituli] 2019) for all meteorological inputs (precipitation, temperature, solar radiation, and relative humidity) except wind speed. We used simulated wind speed data based on a modified exponential equation in SWAT (Neitsch et al. 2011). These meteorological stations cover the entire Simojoki catchment at 26 different locations. The long-term (1996–2015) average annual temperature in the area varied between −3.83 and 4.94 °C, and the long-term annual average precipitation ranged from 580 to 746 mm. The mean monthly temperature (1996–2015) varied between 20.5 °C (July) and −16.3 °C (January). Average monthly precipitation ranged from approximately 30 mm (January) to 83 mm (July).

We collected daily streamflow and concentration data of Org-N, TN, Org-P, TP, and SS from the Hertta database (Suomen ympäristökeskus [SYKE] 2020). These monitoring data were used as the observed time series in the SWAT model to account for temporal changes during calibration and validation in two different periods. Since 1996, the sampling frequency at both outlets (Figure 1(a)) of the Simojoki catchment has been 11–25 samples per year (Lepistö et al. 2014). We used two monitoring stations to calibrate the model (i) Hosionkoski sub-catchment in the central part of the catchment and (ii) Simojoki at the catchment outlet (average summary is shown in Table 1). To estimate the load from the measured streamflow and concentration data, we applied a regression method (Walker 1996) and followed the steps described in a previous study (Bhattacharjee et al. 2021b).

Table 1

A summary of the monitoring data from 1996 to 2015

Catchment and sub-catchmentFlow (m3/s)SS (mg/l)Total N (μg/l)Total P (μg/l)
Hosionkoski 25.7 9.53 626.5 48.5 
Simojoki 40.6 48.08 826 106 
Catchment and sub-catchmentFlow (m3/s)SS (mg/l)Total N (μg/l)Total P (μg/l)
Hosionkoski 25.7 9.53 626.5 48.5 
Simojoki 40.6 48.08 826 106 

SWAT modelling, calibration, and validation at Simojoki

We used the Penman–Monteith equation to estimate potential evapotranspiration and the variable storage routing method to define flow routing. First, we focussed on calibrating the SWAT model for flow and then for SS, Org-P, TP, Org-N, and TN sequentially. The full SWAT calibration period was from 1985 to 2005, with the first 10 years as the warm-up period.

We chose the SUFI2 algorithm in the SWAT-CUP (V-5.1.6.2) to perform semi-automatic calibration (Abbaspour 2015). First, we analysed the sensitivities of these parameters. As mentioned in earlier studies (Arnold et al. 2012), there are some calibration parameters in SWAT-CUP for which the model result can vary to a certain extent. Sensitivity analysis is an effective method for understanding the effects of the parameters in the model. Initially, we identified multiple parameters for use in the sensitivity analysis. We used parameter values from published literature (Menberu et al. 2021) to set up the initial parameter ranges of the physical and hydrological properties of peatlands. We used a Latin hypercube (LH) sampling scheme to calculate the sensitivity index of the parameters for global sensitivity analysis. The most sensitive parameters, based on t-stat and p-value for the calibration of flow, nitrogen, phosphorus, and SS, are shown in Supplementary Appendix Figure A1 where t-stat is the coefficient of a parameter divided by its standard error. If a coefficient is large compared with its standard error, then it is different from 0 and the parameter is sensitive. The p-value tests the null hypothesis (the coefficient is equal to zero, with no effect). A low p-value (<0.05) indicates rejection of the null hypothesis which represents the parameter associated with the response variable or affected by it (Abbaspour 2015).

We found some sensitive parameters (i.e., LAT_ORGP) that should be used only for calibration when the observed data are available to justify the selected ranges. As we had data from earlier studies that focussed on parameters related to peatland forestry (Huttunen et al. 2016; Räty et al. 2020), we used these parameters for calibration.

Based on the propagation of the parameter uncertainties, two prominent factors were considered when calibrating the model response variables, as recommended in a previous study (Abbaspour 2015):

  • At least 70% of the measured data (observations) were between the 2.5th and 97.5th prediction percentiles (95PPU) of the model results. The measurement factor by which the uncertainties are addressed is known as the p-factor.

  • The standard deviation of the measured data was greater than the difference between the 2.5th and 97.5th percentile predictions. The gap is known as the thickness of the 95PPU band or the r-factor.

Next, the ranges of parameters that fulfilled the above criteria were considered to check the reasonable agreement between the modelled and measured outputs. The Kling–Gupta efficiency (KGE) (Gupta et al. 2005) was selected as the objective function in the SWAT-CUP as recent studies have reported that KGE is suitable for evaluating the variability, peak, and mean of observations (Pechlivanidis et al. 2011).

Next, we started flow calibration based on the initial values of the sensitive parameters. We used 1,000 simulations in the SWAT-CUP for the first iteration. Because the initial model output did not match the observations well, we ran another iteration with the same number of simulations. After verifying the results of the second iteration, we run another 500 simulations. In each iteration, the parameter ranges were altered to obtain the most plausible model output based on the objective function. Finally, 23 parameters with reasonable ranges that affect the simulated flow (Table 3) were selected.

After calibrating the flow parameters, we started with the SS calibration and then moved to phosphorus and nitrogen. In addition to the flow parameters, we found nine parameters for calibrating the SS, nine for phosphorus, and three for nitrogen. We used three iterations with 500 simulations at each calibration stage to calibrate the nutrients and SS.

Finally, a split-sample validation was performed based on calibrated parameter values (Supplementary Appendix Table A1) in an independent dataset. The applicability of the calibrated model was tested from 2006 to 2015.

Model evaluation and analysis

In addition to KGE, the Nash–Sutcliffe efficiency (NSE) was also used to evaluate the model. Recent studies have argued that using NSE could lead to a poor representation of flows, whereas KGE does not have this issue (Gupta et al. 2005; Pechlivanidis et al. 2011). Thus, we used two metrics for model evaluation to add robustness to the calibrated model. NSE varied between −∞ and 1.0. A value of 1 represents a perfect fit between the simulated and observed values, whereas values less than zero indicate unacceptable performance (Moriasi et al. 2015). We also evaluated model performance using R2 and PBIAS. The R2 evaluates the dynamics between the observed and simulated values, whereas the PBIAS quantifies the difference between the observed and simulated values. Thus, a PBIAS value of zero indicates the best result. Based on the output of these statistical indices, we selected the best simulations to finalize our calibration for the studied catchment.

For the calibration period (1996–2005), the ranges of the fully calibrated parameters were considered to generate the uncertainty (95PPU) boundary, whereas, for the validation (2006–2015), we used a single parameter value for the model evaluation (Supplementary Appendix Table A1). We also added the results of another calibrated model without peatland drainage. We integrated the results of both models with and without peatland drainage only in Table 2 in the Results section because it is important to know how the modelling approaches represent the drained peatland.

Table 2

Performance metrics of the SWAT model (including peatland drainage) of daily flow, nutrients, and SS

Performance metricsHosionkoski
Simojoki
Hosionkoski
Simojoki
CalValCalValCalValCalVal
FlowSS
R2 0.73 (0.67) 0.67 (0.59) 0.72 (0.58) 0.64 (0.54) 0.64 (0.53) 0.33 (0.23) 0.12 (0.09) 0.40 (0.35) 
NSE 0.74 (0.6) 0.63 (0.57) 0.76 (0.53) 0.67 (0.53) 0.49 (0.52) 0.31 (0.21) 0.12 (0.09) 0.4 (0.23) 
PBIAS 9.9 (16.5) 12 (23.3) 9.6 (18.6) 9.7 (14) 52.1 (−17.2) 30.7 (34.6) 1.1 (−3.3) 18.6 (−35.4) 
KGE 0.84 (0.73) 0.74 (0.6) 0.81 (0.68) 0.77 (0.61) 0.35 (0.57) 0.33 (0.22) 0.04 (0.03) 0.33 (0.32) 
 Organic N Total N 
R2 0.82 (0.65) 0.63 (0.53) 0.57 (0.53) 0.68 (0.54) 0.82 (0.55) 0.64 (0.54) 0.55 (0.55) 0.66 (0.56) 
NSE 0.81 (0.65) 0.62 (0.52) 0.56 (0.53) 0.65 (0.53) 0.74 (0.62) 0.61 (0.51) 0.55 (0.53) 0.65 (0.63) 
PBIAS 0.6 (8.9) 9.4 (19.4) 11.2 (23.3) 16.1 (16.1) 20 (−24) 6.1 (−16) 3.2 (10.8) 9.8 (19.8) 
KGE 0.9 (0.71) 0.74 (0.64) 0.59 (0.52) 0.63 (0.61) 0.73 (0.63) 0.78 (0.68) 0.61 (0.51) 0.66 (0.62) 
 Organic P Total P 
R2 0.75 (0.65) 0.66 (0.56) 0.42 (0.34) 0.58 (0.47) 0.61 (0.4) 0.65 (0.45) 0.29 (0.23) 0.6 (0.3) 
NSE 0.15 (−75.6) 0.63 (−5.37) 0.38 (−90) 0.56 (−5.1) 0.41 (−20.5) 0.65 (0.02) 0.19 (0.27) 0.51 (−1.1) 
PBIAS 13.2 (−900) 12.2 (−280) 17.3 (−19) 5.8 (−305) 20.4 (−460) 7.3 (−41) 38.7 (−11.6) 25.6 (−50.6) 
KGE 0.35 (−26) 0.77 (2.52) 0.34 (−15) 0.56 (2.47) 0.62 (−12.8) 0.7 (0.36) 0.23 (0.21) 0.4 (0.04) 
Performance metricsHosionkoski
Simojoki
Hosionkoski
Simojoki
CalValCalValCalValCalVal
FlowSS
R2 0.73 (0.67) 0.67 (0.59) 0.72 (0.58) 0.64 (0.54) 0.64 (0.53) 0.33 (0.23) 0.12 (0.09) 0.40 (0.35) 
NSE 0.74 (0.6) 0.63 (0.57) 0.76 (0.53) 0.67 (0.53) 0.49 (0.52) 0.31 (0.21) 0.12 (0.09) 0.4 (0.23) 
PBIAS 9.9 (16.5) 12 (23.3) 9.6 (18.6) 9.7 (14) 52.1 (−17.2) 30.7 (34.6) 1.1 (−3.3) 18.6 (−35.4) 
KGE 0.84 (0.73) 0.74 (0.6) 0.81 (0.68) 0.77 (0.61) 0.35 (0.57) 0.33 (0.22) 0.04 (0.03) 0.33 (0.32) 
 Organic N Total N 
R2 0.82 (0.65) 0.63 (0.53) 0.57 (0.53) 0.68 (0.54) 0.82 (0.55) 0.64 (0.54) 0.55 (0.55) 0.66 (0.56) 
NSE 0.81 (0.65) 0.62 (0.52) 0.56 (0.53) 0.65 (0.53) 0.74 (0.62) 0.61 (0.51) 0.55 (0.53) 0.65 (0.63) 
PBIAS 0.6 (8.9) 9.4 (19.4) 11.2 (23.3) 16.1 (16.1) 20 (−24) 6.1 (−16) 3.2 (10.8) 9.8 (19.8) 
KGE 0.9 (0.71) 0.74 (0.64) 0.59 (0.52) 0.63 (0.61) 0.73 (0.63) 0.78 (0.68) 0.61 (0.51) 0.66 (0.62) 
 Organic P Total P 
R2 0.75 (0.65) 0.66 (0.56) 0.42 (0.34) 0.58 (0.47) 0.61 (0.4) 0.65 (0.45) 0.29 (0.23) 0.6 (0.3) 
NSE 0.15 (−75.6) 0.63 (−5.37) 0.38 (−90) 0.56 (−5.1) 0.41 (−20.5) 0.65 (0.02) 0.19 (0.27) 0.51 (−1.1) 
PBIAS 13.2 (−900) 12.2 (−280) 17.3 (−19) 5.8 (−305) 20.4 (−460) 7.3 (−41) 38.7 (−11.6) 25.6 (−50.6) 
KGE 0.35 (−26) 0.77 (2.52) 0.34 (−15) 0.56 (2.47) 0.62 (−12.8) 0.7 (0.36) 0.23 (0.21) 0.4 (0.04) 

Cal refers to calibration, and Val refers to validation. Values in parentheses represent the indices for the model without peatland drainage.

We further developed a Python script to analyse the calibration and validation results (Supplementary Appendix C) and extract loading information for any reach, sub-catchment, or HRU from the SWAT output database. The output results of the script can be used to map information regarding the nutrients and SS at different geographical boundaries.

Calibration and validation

Model performance

The performance metrics for flow, nutrients, and SS in the SWAT-CUP varied from very good to unsatisfactory (Moriasi et al. 2015) for the Hosionkoski sub-catchment and the entire Simojoki catchment. We also added indices from the SWAT-CUP model without peatland drainage (in parentheses, Table 2). For all indices, the model performance was better when peatland drainage was implemented within the model. For the flow and nitrogen, the performance metrics were satisfactory or higher for calibration and validation. However, the estimations for phosphorus and SS varied for different indices (Table 2). As the observed SS concentrations were low and mainly contained organic fractions, the SS budget was not closed, which affected the modelling accuracy of the SS and phosphorus.

The R2 and NSE values were higher for the model with peatland drainage, possibly indicating a better simulation of peak flows (Zhang et al. 2015). The KGE values were also a good match throughout the calibration and validation of flow and other nutrients, except for SS and total phosphorus. The PBIAS yielded by the SWAT-CUP also relied on the compensation of over- and underestimations rather than a satisfactory representation. The inclusion of model results with and without peatland drainage is visible in graphical illustrations (Figures 24).
Figure 2

SWAT simulated flow with the observed measurements at the Hosionkoski sub-catchment and the entire Simojoki catchment. The 95PPU boundary is only plotted for the calibration period 1996–2005.

Figure 2

SWAT simulated flow with the observed measurements at the Hosionkoski sub-catchment and the entire Simojoki catchment. The 95PPU boundary is only plotted for the calibration period 1996–2005.

Close modal
Figure 3

Simulated load of Org-N and TN at the Hosionkoski sub-catchment and the Simojoki catchment.

Figure 3

Simulated load of Org-N and TN at the Hosionkoski sub-catchment and the Simojoki catchment.

Close modal
Figure 4

Modelled Org-P, TP, and SS load at the Hosionkoski sub-catchment and the Simojoki catchment.

Figure 4

Modelled Org-P, TP, and SS load at the Hosionkoski sub-catchment and the Simojoki catchment.

Close modal

Hydrological mass balance

Based on the input and output components of the SWAT model, GW flow (shallow aquifer) had the highest share of water balance in the Simojoki catchment, followed by ET and the other components (Table 3).

Table 3

Water balance components after simulating the SWAT model

Water balance component (IN), mm/yearWater balance component (OUT), mm/year
Rainfall 374.24 Surface runoff 4.46 
Snowfall 259.26 Lateral flow 52.41 
  GW (shallow aquifer) flow 307.31 
  GW (deep aquifer) flow 1.45 
  Evapotranspiration (ET) 259 
 633.5  624.63 
Water balance component (IN), mm/yearWater balance component (OUT), mm/year
Rainfall 374.24 Surface runoff 4.46 
Snowfall 259.26 Lateral flow 52.41 
  GW (shallow aquifer) flow 307.31 
  GW (deep aquifer) flow 1.45 
  Evapotranspiration (ET) 259 
 633.5  624.63 

We also compared the simulated potential ET with the potential ET available from gridded data (1 km × 1 km) from the FMI (Pirinen et al. 2022). From 1996 to 2005 (April to September of each year), we compared the average monthly potential ET (Supplementary Appendix Figure E1) and found that the difference was always less than 10%.

Hydrology and water quality

During the calibration period (1996–2005), based on the model results with peatland drainage, the uncertainty boundary (95PPU) captured 74% of the observed flow at Simojoki station, whereas the percentage was 70% for the upstream Hosionkoski catchment. The percentage was almost the same (69%) during the validation period in the Simojoki catchment area. The baseflow was separated from the simulated flow. After separation, we found that the simulated baseflow ratio was within 20% of the observed baseflow ratio. This can be considered as a benchmark for satisfactory calibration, as suggested in an earlier study (Bracmort et al. 2006). Gridded meteorological inputs were used to achieve better spatial coverage. Thus, these data minimized the uncertainty issues that could affect the simulated flow. ArcSWAT considers a meteorological station near the centroid of each sub-catchment. For low-flow conditions, the modelling results followed the pattern of the observations, but not for high-flow conditions caused by snowmelt (Figure 2).

For Org-N and TN, the uncertainty band estimated by SWAT-CUP included more than 82% of the observations in both catchments for calibration. The calibrated Org-N differed from the observations by 0.6% in the Hosionkoski sub-catchment and 11.2% in the Simojoki catchment. For validation, the differences were 9.4 and 16.1% for Hosionkoski and Simojoki, respectively. The difference between the simulated and observed TN values varied from 20 to 3.2% from upstream to downstream (Figure 3).

The ratio between the average Org-N concentration of the observations and simulations was 1.04 for the Simojoki catchment, whereas it was 0.9 for the Hosionkoski catchment (Supplementary Appendix Figure B1). The average TN concentration ratio between the simulations and observations for both catchments is 1.1.

The difference in TP load between the simulations and observations in the Hosionkoski catchment was approximately 20%, whereas it was higher (38%) in the Simojoki catchment (Figure 4). We found that 60% of the observations fell within the 95% prediction uncertainty (95PPU) boundary generated from the calibration results.

The simulated average TP concentration was 21.8 μg/l in the Simojoki catchment, which was 4.6% higher than the observations. For the Hosionkoski catchment, the simulated TP values differed by 25% from observations (Supplementary Appendix Figure B2). The simulation and observation analysis showed a higher concentration in Hosionkoski than in the Simojoki catchment.

For SS, the simulated model captured only 37 and 38% of the observations from 1996 to 2005 within the uncertainty boundary (95PPU) in the Hosionkoski sub-catchment and the Simojoki catchment, respectively (Figure 4). Although the prediction of the average SS load differed from the observations by 1.1% in the Simojoki catchment, the difference was high (∼50%) in the Hosionkoski catchment, as imbalances might also exist in the nutrient pool of the SWAT at the sub-catchment level.

The average simulated concentration of SS was 7.6 mg/l, whereas it was 3.8 mg/l for the observed measurements in the Simojoki catchment (see Supplementary Appendix Figure B2). The simulated concentration in the Hosionkoski sub-catchment was 1.5 times higher than the observed concentrations. Although the simulated suspended sediment (SS) concentrations exhibited significant variations compared with the observations in both catchments, the ratio of the average SS concentration between the observed and simulated data remained nearly constant.

Loads and concentrations of nutrients and SS

Annual and monthly variation

The effect of melted snow was visible in the monthly variations in nutrients and SS, which showed the highest load in May. In February, the load was minimal for both catchments (Figure 5 and Supplementary Appendix Figure D1). For all nutrients and SS, the monthly variation in average concentration followed the same pattern (Figure 5) as the monthly average load.
Figure 5

Annual (left) and monthly average (right) nutrient and SS loads and concentrations at the Simojoki outlet.

Figure 5

Annual (left) and monthly average (right) nutrient and SS loads and concentrations at the Simojoki outlet.

Close modal

Spatial and temporal distribution

The variations in nutrient and SS loads between the headwaters and downstream reaches changed substantially from 1996 to 2015, as the peatland drainage activities were higher in the downstream Simojoki. Based on the suggested colour pattern of the ecological status of the surface water bodies by a previous water framework directive (analysis period 2012–2017; Suomen ympäristökeskus [SYKE] 2021), we mapped the load variation throughout the catchment (using the Python script in Supplementary Appendix C). The blue colour indicates the excellent status of nutrients and SS in the reaches, whereas the red colour indicates a poor status (Figure 6).
Figure 6

Spatial and temporal variation of Org-N, TN, Org-P, TP, and SS in 1996 and 2015. The left section shows the variations of load at the beginning (1996) and the end (2015) of the simulation. The variation of colour is chosen based on the suggested colour pattern of the ecological status of the surface water body by a previous Water Framework Directive (Suomen ympäristökeskus [SYKE] 2021).

Figure 6

Spatial and temporal variation of Org-N, TN, Org-P, TP, and SS in 1996 and 2015. The left section shows the variations of load at the beginning (1996) and the end (2015) of the simulation. The variation of colour is chosen based on the suggested colour pattern of the ecological status of the surface water body by a previous Water Framework Directive (Suomen ympäristökeskus [SYKE] 2021).

Close modal

The Org-N load was 80% higher in 2015 than that in 1996 at the Simojoki outlet, and this increase was even more pronounced for the TN load. However, the load ratio between the most downstream and upstream reaches was almost the same in 1996 and 2015 for both Org-N and TN. We found the highest load of 748 ton/year Org-N and 907 ton/year TN at the Simojoki outlet in the catchment.

Similar to Org-N and TN, Org-P and TP were higher in the middle and downstream parts of the Simojoki catchment. Between the most upstream and downstream reaches, the ratio of Org-P and TP loads was higher in 2015 than that in 1996. A total of 159 ton/year of Org-P and 169 ton/year of TP, the highest total load in a year, were found at the outlet of the Simojoki catchment.

The SS load at the furthest downstream reach in the Simojoki catchment was four times the SS load found at the furthest upstream reach in 2015 and almost five times the simulated load in 1996. In contrast to the nutrients, the highest total average SS load was found at a reach (reach 98) which was located above the most downstream outlet of the Simojoki catchment (Supplementary Appendix Figure F1).

Load variation at the sub-catchment level

The highest loads of Org-N and TN were found in the furthest downstream sub-catchment, whereas for phosphorus and SS, the case was different. For TP and TN, the average load was 0.02 and 0.1 ton/ha/year, respectively (Figure 7). There was a slight difference between the organic and total nutrient contents within each sub-catchment. The average SS was 1.92 ton/ha/year.
Figure 7

SS and nutrient variation in the sub-catchments and land uses from 1996 to 2015. The top panel shows the spatial distribution of loading in each sub-catchment. The bottom two panels represent loads of nutrients and SS from different land uses for the same period. The top-right axis with transparent colours shows the percentage area of each land use in the Simojoki catchment.

Figure 7

SS and nutrient variation in the sub-catchments and land uses from 1996 to 2015. The top panel shows the spatial distribution of loading in each sub-catchment. The bottom two panels represent loads of nutrients and SS from different land uses for the same period. The top-right axis with transparent colours shows the percentage area of each land use in the Simojoki catchment.

Close modal

The loading of nutrients in the peat soil in the clear-cut areas was the highest compared with all other land use categories. However, forests in the peat soil cover a large portion (31%) of the catchment. The loading of the mineral soil was also higher in the clear-cut areas than in the forest. From 1996 to 2015, the differences in loading between forests on peat soil and forests on mineral soil were 0.002 ton/ha/year for Org-N and TN, 0.0002 ton/ha/year for Org-P and TP, and 0.022 ton/ha/year for SS. Although peat extraction had the lowest area coverage within the Simojoki catchment, we observed the highest loading (1.03 ton/ha/year) of SS in that land use category.

The SWAT model has been used in various catchments worldwide and under northern conditions (Grizzetti et al. 2003). However, previous studies have not evaluated the SWAT performance in peatland-dominated catchments with intensive drainage activities, especially for various land uses in organic soils (forestry). The approach of burning the drainage network (collected from Bhattacharjee et al. (2021a)) into the DEM in a forestry catchment and integrating the process within the SWAT model resulted in good to satisfactory model performance for flow and nitrogen and slightly lower performance for phosphorus and SS. Specifically,

  • I. The sensitive parameters identified in this study can be used to simulate the nutrients and SS in similar catchments. The gap that remains in this study is answering the question of why certain ranges were selected for specific parameters based on the sensitivity analysis. We relied on previous studies (e.g., Huttunen et al. 2016; Räty et al. 2020; Menberu et al. 2021) to select parameter ranges related to peat properties. Thus, the uncertainties in the current SWAT model can be reduced if there is an opportunity to obtain physical measurements of these parameters to support selected ranges.

  • II. The nutrient and SS movement processes in multi-pools in SWAT depend heavily on the chosen parameters and their ranges (Supplementary Appendix Table A1), especially for forest land and organic soils. This process is more complex when the catchment has a drained peatland area of more than 30% (Rankinen et al. 2006). For example, owing to the organic fraction and the low concentration of SS and phosphorus, there are mismatches between the modelled output and the estimated loads from the monitored observations. Many parameters were dropped during the calibration process because of their insensitivity to peatland forestry. The uncertainty (95PPU) in estimates captured 35–81% of the observations that involved all plausible sources of uncertainties related to the model structure and the trade-off between parameter values and the objective function (Gupta et al. 2005).

  • III. The current SWAT model showed a PBIAS of 9.6% between the simulated flow and the observations. TIMP (snowpack temperature lag factor) was the most sensitive parameter for calibrating flow. The amount of snowmelt is affected by the lagging factor because it accounts for the snowpack density, snowpack depth, and other exposures (Neitsch et al. 2011). As the range of TIMP was between 0.64 and 0.84 in the current SWAT model, we found that the mean air temperature for each day during the studied period had a larger impact on the snowpack temperature on that specific day than the snowpack temperature for the preceding day.

  • IV. The simulation of SS loading is challenging as transported SS in peatland-dominated catchments typically consists mostly of organic particles (Marttila & Kløve 2015) with very fine particle sizes (Marttila et al. 2016). These physical processes, including the slow settling velocity and wash load nature, cannot be fully simulated using the SWAT. Compared with the observations (9,128 ton/year), the SWAT model generated approximately 11,292 ton/year SS load at the Simojoki outlet, whereas the simulated load was 1.4 times the observed load at the Hosionkoski outlet.

  • V. Within the soil profile, the soil surface contained the highest amount of nutrients and organic matter. As organic soils have a distinct process from mineral soils (Menberu et al. 2021), the movement of nitrogen from the fresh organic pool to the inorganic mineral pool in the HRUs was insignificant compared with the concentration of Org-N in the baseflow (parameters of nitrogen calibration in Supplementary Appendix Table A1). The differences between the simulated load and the observations for Org-N and TN varied by 11.2 and 3.2%, respectively. In a previous study (Huttunen et al. 2016), the authors found a 1,012 ton/year gross TN load between 1991 and 2011, whereas the current SWAT model resulted in a 907 ton/year TN load. A minimal change in the Org-N concentration in the baseflow substantially affected model calibration. Saturated hydraulic conductivity and bulk density in the peat soil also affected nitrogen processes in the catchment, as nitrogen was highly reactive in the peat soil.

Effects of drainage activities

Some inorganic processes might have occurred in the mineral soil layer of the catchment, but during 1996–2015, denitrification was the dominant process. When the temperature was 5 °C or higher on a specific day, nitrification and volatilization could also occur in the catchment (Neitsch et al. 2011). Most plant-essential nutrients are cations that are attracted to and sorbed onto negatively charged soil particles. As plants extract these cations from the soil solution and the peat soil particles are positively charged, the leached nitrogen in the current SWAT model was not noticeable compared with nitrogen in the surface runoff. Although the Simojoki catchment has lakes covering 5.7% of its total area, nutrient retention within the catchment is not sensitive to lake effects. The parameters responsible for calibrating Org-N also fit well with TN calibration. However, peat properties are site specific and may vary owing to the different peat decomposition levels in different peatlands.

In response to a concentration gradient, phosphorus usually moves through the migration of ions over small distances on the soil surface in SWAT (Neitsch et al. 2011). Solubility is limited to phosphorus in most environments compared with nitrogen, which is highly mobile (Chaubey et al. 2006). Peat soils contained more than 50% phosphorus, whereas mineral soils contained approximately 40%. Runoff in the drained peat soils was higher than that in the undisturbed soils, which also affected the phosphorus concentration. Similarly, the concentration of Org-P in the baseflow also influenced phosphorus calibration, as it was the most sensitive parameter. Because peat soils have high water-holding capacity, they usually have little surface runoff. After the snow melting period, TP was higher than in the other seasons, during which stagnant water in peat soils might affect phosphorus loading (Huttunen et al. 2016).

Previous studies have already mentioned (Joensuu et al. 2002; Nieminen et al. 2018) that drainage activities and ditch maintenance are the main reasons for the increased nutrient and SS loads and concentrations in peatland-dominated catchments. In a previous study (Lepistö et al. 2008), the authors noticed an increasing trend in Org-N and found that changes in hydrological dynamics and soil temperature during low flows might affect its concentration. Temporal estimates using the current SWAT model for the Simojoki catchment showed an increasing nitrogen trend.

Clear-cutting in peat soils has not been implemented on a large scale in the Simojoki catchment area. However, the loading is predicted to increase in the future (Lepistö et al. 2008). From the spatial distribution of loading in different land uses (Figure 7), the clear-cut peat soil contained the maximum loading of all nutrients. In a previous study based on an export coefficient database, the authors also found the maximum loading for the same land use (Launiainen et al. 2014; Bhattacharjee et al. 2021b). However, the loading estimates from the export coefficient differ from those of the current SWAT model. For example, the authors mentioned 34 kg/ha/year TN, 3.9 kg/ha/year TP, and 2,400 kg/ha/year SS, whereas we noticed 50 kg/ha/year TN, 6 kg/ha/year TP, and 1,200 kg/ha/year SS from the forest in the peat soil. A plausible explanation for the difference in loads could be the nutrient decay processes in large rivers, which were not covered by the export coefficients. The export values were monitored in different empirical studies that included seasonal loading information from different land use activities but did not include decay processes from land to sea in the surface water system. In addition, the SWAT model is needed in this current study to provide a more detailed, process-based understanding of nutrient export and its spatial and temporal variability, which the export coefficient model cannot fully capture. While both methods offer valuable insights, the SWAT model's ability to simulate hydrological processes, nutrient decay and transformation, and the effects of management practices under changing conditions make it essential for addressing the research questions in this study.

Future perspectives

Historical land use activities may show impacts for years or even decades. For example, the long-term prediction of Org-N in this study resulted in increasing trends from drained peatlands, which was in line with another recent study (Lepistö et al. 2021). Therefore, from a land manager's perspective, it is beneficial to know the spatial and temporal distribution of the drainage network, nutrients, and SS loading. Especially in the catchments containing a high percentage of peat soil, a proper management plan is essential.

Therefore, the modelling study can be used to understand the comprehensive status of a catchment, starting from knowing peatland drainage conditions to considering future projections. The study can also be used to monitor or achieve the good ecological status of surface water bodies based on the guidelines of the Water Framework Directive (WFD 2022–2027) (Suomen ympäristökeskus [SYKE] 2021).

For instance, the catchment manager can use the modelling approach to focus on the catchment-scale consequences of land use changes for nutrient loads and SS exports. However, including more detailed spatial information about drainage activities and the effects of hydrometeorological inputs would make the proposed approach more practical. Thus, the manager can apply the peatland drainage activity in the SWAT model in future hydrological modelling efforts by considering the impact of climate change.

The application of the SWAT tool in a drained peatland catchment demonstrated satisfactory calibration for flow and nitrogen, although there were differences in modelling SS and phosphorus. Despite these uncertainties, the selected parameters and ranges provide a good foundation for simulating physical processes in peatland-dominated catchments.

The observed increasing trend in Org-N during long-term predictions of drained peatlands is particularly important. Despite showing a PBIAS of 9.6% for simulated flow, the current SWAT model exhibited variations in loading between simulations and observations for SS, Org-N, TN, Org-P, and TP. Spatially, the clear-cut area's peat soil displayed the highest loading compared with other land uses, with the downstream outlet of the Simojoki catchment experiencing the maximum annual loading of nutrients, excluding SS and phosphorus.

While forestry practices are inherently dependent on local interpretations, our approach to integrating peatland drainage into SWAT holds promise for future hydrological modelling endeavours. The spatial estimation of nutrient distribution and SS loading has substantial environmental management implications, offering valuable insights for formulating catchment management strategies and implementing water protection measures.

This work was part of the Nordic Centre of Excellence BIOWATER, funded by Nordforsk under project number 82263.

J.B. and H.M. conceived the idea and developed the methodology. J.B. and H.M. performed formal analysis and wrote the first draft of the manuscript. E.M.N. reviewed and edited the manuscript. B.K. supervised and edited the manuscript.

The authors declare that they have the consent to participate.

The authors declare that they have the consent to publish.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Abbaspour
K. C.
(
2015
)
SWAT-CUP: SWAT Calibration and Uncertainty Programs – A User Manual
.
Swiss Federal Institute of Aquatic Science and Technology (EAWAG)
,
Dübendorf, Switzerland
.
Arnold
J. G.
,
Moriasi
D. N.
,
Gassman
P. W.
,
Abbaspour
K. C.
,
White
M. J.
,
Srinivasan
R.
,
Santhi
C.
,
Harmel
R. D.
,
Griensven
A. v.
,
Liew
M. W. V.
,
Kannan
N.
&
Jha
M. K.
(
2012
)
SWAT: Model use, calibration, and validation
,
Trans. ASABE
,
55
,
1491
1508
.
https://doi.org/10.13031/2013.42256
.
Bhattacharjee
J.
,
Marttila
H.
,
Haghighi
A. T.
,
Saarimaa
M.
,
Tolvanen
A.
,
Lepistö
A.
,
Futter
M. N.
&
Kløve
B.
(
2021a
)
Development of aerial photos and LIDAR data approaches to map spatial and temporal evolution of ditch networks in peat-dominated catchments
,
J. Irrig. Drain. Eng.
,
147
,
04021006
.
https://doi.org/10.1061/(ASCE)IR.1943-4774.0001547
.
Bhattacharjee
J.
,
Marttila
H.
,
Launiainen
S.
,
Lepistö
A.
&
Kløve
B.
(
2021b
)
Combined use of satellite image analysis, land-use statistics, and land-use-specific export coefficients to predict nutrients in drained peatland catchment
,
Sci. Total Environ.
,
779
,
146419
.
https://doi.org/10.1016/j.scitotenv.2021.146419
.
Bracmort
K. S.
,
Arabi
M.
,
Frankenberger
J. R.
,
Engel
B. A.
&
Arnold
J. G.
(
2006
)
Modeling long-term water quality impact of structural BMPs
,
Trans. ASABE
,
49
,
367
374
.
https://doi.org/10.13031/2013.20411
.
Chaubey
I.
,
Migliaccio
K. W.
,
Green
C. H.
,
Arnold
J. G.
&
Srinivasan
R.
(
2006
)
Phosphorus modeling in soil and water assessment tool (SWAT) model
. In:
Radcliffe, D. E. & Cabrera, M. L. (eds.)
,
Modeling Phosphorus in the Environment
,
CRC Press
,
New York
, pp.
163
187
.
Finér
L.
,
Lepistö
A.
,
Karlsson
K.
,
Räike
A.
,
Härkönen
L.
,
Huttunen
M.
,
Joensuu
S.
,
Kortelainen
P.
,
Mattsson
T.
,
Piirainen
S.
,
Sallantaus
T.
,
Sarkkola
S.
,
Tattari
S.
&
Ukonmaanaho
L.
(
2021
)
Drainage for forestry increases N, P and TOC export to boreal surface waters
,
Sci. Total Environ.
,
762
,
144098
.
https://doi.org/10.1016/j.scitotenv.2020.144098
.
Geological Survey of Finland [GTK]
(
2023
)
Soil Data
[WWW Document]. Available at: https://hakku.gtk.fi/en/locations/search (Accessed 17 January 2023)
.
Greenup
A. L.
,
Bradford
M. A.
,
Mcnamara
N. P.
,
Ineson
P.
&
Lee
J. A.
(
2000
)
The role of Eriophorum vaginatum in CH4 flux from an ombrotrophic peatland
,
Plant Soil
,
227
,
265
272
.
https://doi.org/10.1023/A:1026573727311
.
Grizzetti
B.
,
Bouraoui
F.
,
Granlund
K.
,
Rekolainen
S.
&
Bidoglio
G.
(
2003
)
Modelling diffuse emission and retention of nutrients in the Vantaanjoki watershed (Finland) using the SWAT model
,
Ecol. Modell.
,
169
,
25
38
.
https://doi.org/10.1016/S0304-3800(03)00198-4
.
Gupta
H. V.
,
Beven
K. J.
&
Wagener
T.
(
2005
)
Model calibration and uncertainty estimation
. In: Anderson, M. G. (ed.),
Encyclopedia of Hydrological Sciences
,
John Wiley & Sons
,
New York
, pp.
2015
2031
.
Her
Y.
,
Frankenberger
J.
,
Chaubey
I.
&
Srinivasan
R.
(
2015
)
Threshold effects in HRU definition of the soil and water assessment tool
,
Trans. ASABE
,
58
,
367
378
.
https://doi.org/10.13031/trans.58.10805
.
Holden
J.
,
Chapman
P. J.
&
Labadz
J. C.
(
2004
)
Artificial drainage of peatlands: Hydrological and hydrochemical process and wetland restoration
,
Prog. Phys. Geogr. Earth Environ.
,
28
,
95
123
.
https://doi.org/10.1191/0309133304pp403ra
.
Holden
J.
,
Evans
M. G.
,
Burt
T. P.
&
Horton
M.
(
2006
)
Impact of land drainage on peatland hydrology
,
J. Environ. Qual.
,
35
,
1764
.
https://doi.org/10.2134/jeq2005.0477
.
Huttunen
I.
,
Huttunen
M.
,
Piirainen
V.
,
Korppoo
M.
,
Lepistö
A.
,
Räike
A.
,
Tattari
S.
&
Vehviläinen
B.
(
2016
)
A national-scale nutrient loading model for Finnish watersheds – VEMALA
,
Environ. Model. Assess.
,
21
,
83
109
.
https://doi.org/10.1007/s10666-015-9470-6
.
Joensuu
S.
,
Ahti
E.
&
Vuollekoski
M.
(
2002
)
Effects of ditch network maintenance on the chemistry of run-off water from peatland forests
,
Scand. J. For. Res.
,
17
,
238
247
.
https://doi.org/10.1080/028275802753742909
.
Joosten
H.
&
Clarke
D.
(
2002
)
Wise Use of Mires and Peatlands – Background and Principles Including a Framework for Decision-Making
.
International Mire Conservation Group and International Peat Society
.
Kadlec
R. H.
&
Wallace
S.
(
2008
)
Treatment Wetlands
.
CRC Press
,
New York
.
Knisel
W. G.
(
1982
)
CREAMS: A Field-Scale Model for Chemicals, Runoff, and Erosion from Agricultural Management Systems
.
United States Department of Agriculture, Conservation Research Report 26, 640 pp. Available online: https://www.tucson.ars.ag.gov/unit/publications/PDFfiles/312.pdf
.
Kreutzweiser
D. P.
,
Hazlett
P. W.
&
Gunn
J. M.
(
2008
)
Logging impacts on the biogeochemistry of boreal forest soils and nutrient export to aquatic systems: A review
,
Environ. Rev.
,
16
,
157
179
.
https://doi.org/10.1139/A08-006
.
Landry
J.
&
Rochefort
L.
(
2012
)
The Drainage of Peatlands: Impacts and Rewetting Techniques
.
Peatland Ecology Research Group
,
Quebec
.
Launiainen
S.
,
Sarkkola
S.
,
Laurén
A.
,
Puustinen
M.
,
Tattari
S.
,
Mattsson
T.
,
Piirainen
S.
,
Heinonen
J.
,
Alakukku
L.
&
Finér
L.
(
2014
)
KUSTAA-työkalu valuma-alueen vesistö-kuormituksen laskentaan
.
Suomen Ympäristökeskuksen julkaisuja, 33/2014. 55 pp. ISBN: 978-952-11-4374-8. https://helda.helsinki.fi/bitstream/handle/10138/144108/SYKEra_33_2014.pdf?sequence=1.
Lepistö
A.
&
Syri
S.
(
2001
)
Modelling effects of changing deposition and forestry on nitrogen fluxes in a northern river basin
,
Water Air Soil Pollut.
,
130
,
1691
1696
.
https://doi.org/10.1023/A:1013902424916
.
Lepistö
A.
,
Kortelainen
P.
&
Mattsson
T.
(
2008
)
Increased organic C and N leaching in a northern boreal river basin in Finland
,
Global Biogeochem. Cycles
,
22
,
1
10
.
https://doi.org/10.1029/2007GB003175
.
Lepistö
A.
,
Futter
M. N.
&
Kortelainen
P.
(
2014
)
Almost 50 years of monitoring shows that climate, not forestry, controls long-term organic carbon fluxes in a large boreal watershed
,
Global Change Biol.
,
20
,
1225
1237
.
https://doi.org/10.1111/gcb.12491
.
Lepistö
A.
,
Räike
A.
,
Sallantaus
T.
&
Finér
L.
(
2021
)
Increases in organic carbon and nitrogen concentrations in boreal forested catchments – Changes driven by climate and deposition
,
Sci. Total Environ.
,
780
,
146627
.
https://doi.org/10.1016/J.SCITOTENV.2021.146627
.
Marttila
H.
,
Nieminen
M. L.
,
Ärje
J.
,
Meissner
K.
,
Tuukkanen
T.
,
Saukkoriipi
J.
&
Kløve
B.
(
2016
)
Dynamics and temporal changes in suspended sediment transport in northern Finland: Transport of very fine particulate matter
,
Boreal Environ. Res.
,
21
,
541
555
.
Marttila
H.
,
Karjalainen
S. M.
,
Kuoppala
M.
,
Nieminen
M. L.
,
Ronkanen
A. K.
,
Kløve
B.
&
Hellsten
S.
(
2018
)
Elevated nutrient concentrations in headwaters affected by drained peatland
,
Sci. Total Environ.
,
643
,
1304
1313
.
https://doi.org/10.1016/j.scitotenv.2018.06.278
.
Melaku
N. D.
,
Wang
J.
&
Meshesha
T. W.
(
2020
)
Improving hydrologic model to predict the effect of snowpack and soil temperature on carbon dioxide emission in the cold region peatlands
,
J. Hydrol.
,
587
,
124939
.
https://doi.org/10.1016/J.JHYDROL.2020.124939
.
Menberu
M. W.
,
Marttila
H.
,
Ronkanen
A.-K.
,
Haghighi
A. T.
&
Kløve
B.
(
2021
)
Hydraulic and physical properties of managed and intact peatlands: Application of the Van Genuchten-Mualem models to peat soils
,
Water Resour. Res.
,
57
,
e2020WR028624
.
https://doi.org/10.1029/2020WR028624
.
Molina-Navarro
E.
,
Andersen
H. E.
,
Nielsen
A.
,
Thodsen
H.
&
Trolle
D.
(
2017
)
Quantifying the combined effects of land use and climate changes on stream flow and nutrient loads: A modelling approach in the Odense Fjord catchment (Denmark)
,
Sci. Total Environ.
,
621
,
253
264
.
https://doi.org/10.1016/j.scitotenv.2017.11.251
.
Moriasi
D. N.
,
Gitau
M. W.
,
Pai
N.
&
Daggupati
P.
(
2015
)
Hydrologic and water quality models: Performance measures and evaluation criteria
,
Trans. ASABE
,
58
,
1763
1785
.
https://doi.org/10.13031/trans.58.10715
.
Neitsch
S. L.
,
Arnold
J. G.
,
Kiniry
J. R.
&
Williams
J. R.
(
2011
)
Soil & Water Assessment Tool Theoretical Documentation
.
Texas Water Resources Institute
,
College Station
.
Nieminen
M.
,
Sallantaus
T.
,
Ukonmaanaho
L.
,
Nieminen
T. M.
&
Sarkkola
S.
(
2017a
)
Nitrogen and phosphorus concentrations in discharge from drained peatland forests are increasing
,
Sci. Total Environ.
,
609
,
974
981
.
https://doi.org/10.1016/j.scitotenv.2017.07.210
.
Nieminen
M.
,
Piirainen
S.
,
Sikström
U.
,
Löfgren
S.
,
Marttila
H.
,
Sarkkola
S.
,
Laurén
A.
&
Finér
L.
(
2018
)
Ditch network maintenance in peat-dominated boreal forests: Review and analysis of water quality management options
,
Ambio
,
47
,
535
545
.
https://doi.org/10.1007/s13280-018-1047-6
.
NLSF [Paituli]
(
2019
)
Paituli – Spatial Data Download Service Portal
[WWW Document]. Available at: https://avaa.tdata.fi/web/avaa/-/paituli-paikkatietopalvelu (Accessed 29 March 2019)
.
Paavilainen
E.
&
Päivänen
J.
(
1995
)
Peatland Forestry, Ecological Studies
.
Berlin, Heidelberg
:
Springer Berlin Heidelberg
.
Pechlivanidis
I. G.
,
Jackson
B. M.
,
Mcintyre
N. R.
&
Wheater
H. S.
(
2011
)
Catchment scale hydrological modelling: A review of model types, calibration approaches and uncertainty analysis methods in the context of recent developments in technology and applications
,
Glob. Nest J.
,
13
,
193
214
.
https://doi.org/10.30955/gnj.000778
.
Perkkiö
S.
,
Huttula
E.
&
Nenonen
M.
(
1995
)
Simojoen vesistön vesiensuojelusuunnitelma [Plan for Water System Protection in the Simojoki River System]
.
Vesi- ja ympäristöhallitus
,
Helsinki
.
Pirinen
P.
,
Lehtonen
I.
,
Heikkinen
R. K.
,
Aapala
K.
&
Aalto
J.
(
2022
)
Daily gridded evapotranspiration data for Finland for 1981–2020
,
Ilmastokatsaus
,
4
,
2020
2022
.
https://doi.org/10.35614/issn-2341-6408-ik-2022-11-rl
.
Ramchunder
S. J.
,
Brown
L. E.
&
Holden
J.
(
2012
)
Catchment-scale peatland restoration benefits stream ecosystem biodiversity
,
J. Appl. Ecol.
,
49
,
182
191
.
https://doi.org/10.1111/j.1365-2664.2011.02075.x
.
Rankinen
K.
,
Kenttämies
K.
,
Lehtonen
H.
&
Nenonen
S.
(
2006
)
Nitrogen load predictions under land management scenarios for a boreal river basin in northern Finland
,
Boreal Environ. Res.
,
11
,
213
228
.
Räty
M.
,
Järvenranta
K.
,
Saarijärvi
E.
,
Koskiaho
J.
&
Virkajärvi
P.
(
2020
)
Losses of phosphorus, nitrogen, dissolved organic carbon and soil from a small agricultural and forested catchment in east-central Finland
,
Agric. Ecosyst. Environ.
,
302
,
107075
.
https://doi.org/10.1016/J.AGEE.2020.107075
.
Singh
V. P.
(
2018
)
Hydrologic modeling: Progress and future directions
,
Geosci. Lett.
,
5
,
1
18
.
https://doi.org/10.1186/s40562-018-0113-z
.
Suomen ympäristökeskus [SYKE]
(
2020
)
Open Environmental Information Systems
[WWW Document]. Available at: https://www.syke.fi/fi-FI/Avoin_tieto/Ymparistotietojarjestelmat (Accessed 20 February 2020)
.
Turetsky
M. R.
,
Benscoter
B.
,
Page
S.
,
Rein
G.
,
Van Der Werf
G. R.
&
Watts
A.
(
2014
)
Global vulnerability of peatlands to fire and carbon loss
,
Nat. Geosci.
,
8
,
11
14
.
https://doi.org/10.1038/ngeo2325
.
Turunen
J.
(
2008
)
Development of Finnish peatland area and carbon storage 1950–2000
,
Boreal Environ. Res.
,
13
,
319
334
.
Walker
W. W.
(
1996
)
Simplified Procedures for Eutrophication Assessment and Prediction: User Manual
.
US Army Engineer Waterways Experiment Station
,
Concord
.
Williams
J. R.
&
Berndt
H. D.
(
1977
)
Sediment yield prediction based on watershed hydrology
,
Trans. Am. Soc. Agric. Eng.
,
20
,
1100
1104
.
https://doi.org/10.13031/2013.35710
.
Williams
J. R.
&
Hann
R. W.
(
1978
)
Optimal Operation of Large Agricultural Watersheds with Water Quality Constraints
.
Texas Water Resources Institute
,
Temple
.
Yu
Z.
,
Loisel
J.
,
Brosseau
D. P.
,
Beilman
D. W.
&
Hunt
S. J.
(
2010
)
Global peatland dynamics since the Last Glacial Maximum
,
Geophys. Res. Lett.
,
37
,
43584
.
https://doi.org/10.1029/2010GL043584
.
Zhang
D.
,
Chen
X.
,
Yao
H.
&
Lin
B.
(
2015
)
Improved calibration scheme of SWAT by separating wet and dry seasons
,
Ecol. Modell.
,
301
,
54
61
.
https://doi.org/10.1016/J.ECOLMODEL.2015.01.018
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).

Supplementary data