Controls of Organic Carbon and Nutrient Export from Unmanaged and Managed Boreal Forested Catchments

: Understanding the anthropogenic and natural factors that affect runoff water quality is essential for proper planning of water protection and forest management, particularly in the chang-ing climate. We measured water quality and runoff from 10 unmanaged and 20 managed forested headwater catchments (7–12,149 ha) located in Finland. We used linear mixed effect models to test whether the differences in total organic carbon (TOC), total nitrogen (TN) and total phosphorus (TP) export and concentrations observed can be explained by catchment characteristics, land use, forest management, soil fertility, tree volume and hydrometeorological variables. Results show that much of variation in TOC, TN and TP concentrations and export was explained by drainage, temperature sum, peatland percentage and the proportion of arable area in the catchment. These models explained 45–63% of variation in concentrations and exports. Mean annual TOC export in unmanaged catchments was 56.4 ± 9.6 kg ha − 1 a − 1 , while in managed it was 79.3 ± 3.3 kg ha − 1 a − 1 . Same values for TN export were 1.43 ± 0.2 kg ha − 1 a − 1 and 2.31 ± 0.2 kg ha − 1 a − 1 , while TP export was 0.053 ± 0.009 kg ha − 1 a − 1 and 0.095 ± 0.008 kg ha − 1 a − 1 for unmanaged and managed, respec-tively. Corresponding values for concentrations were: TOC 17.7 ± 2.1 mg L − 1 and 28.7 ± 1.6 mg L − 1 , for TN 420 ± 45 µg L − 1 and 825 ± 51 µg L − 1 and TP 15.3 ± 2.3 µg L − 1 and 35.6 ± 3.3 µg L − 1 . Overall concentrations and exports were signiﬁcantly higher in managed than in unmanaged catchments. Long term temperature sum had an increasing effect on all concentrations and exports, indicating that climate warming may set new challenges to controlling nutrient loads from catchment areas.


Introduction
Boreal catchments are mosaics of upland forests, peatlands and water bodies such as streams, ponds, lakes and rivers. Forest management practices impact carbon (C) and nutrient export to aquatic systems [1][2][3][4]. The effect of different forest management practices on C, nitrogen (N) and phosphorus (P) export to water courses have been studied extensively during the past decades. Paired catchment studies have shown that clearcutting increases N and P export to surface waters for several years [5][6][7][8][9]. Some studies have also reported elevated total organic carbon (TOC) and dissolved organic carbon (DOC) concentrations in runoff from clearcut areas [1,7,[9][10][11]. Moreover, site preparation connected to forest regeneration as well as forest fertilization can increase N and P export [12,13]. The paired catchment studies facilitate statistical detection of treatment effects if the response is greater than the intrinsic uncertainty in the experiment [8]. However, less intensive management, and small, long-term effects are easily masked behind the uncertainty of the experiment and it remains unclear how well the results can be generalized outside the study catchments. Assessing the nutrient export from managed boreal catchments is essential, as they not only affect the water quality inside Finland, but also the Baltic Sea. In the Baltic Sea basin 10 Mha of mires (from total of 19.5 Mha) have been drained for forestry [14].
In Finland, one third of forests are located on peatlands; 23% of total forest area has been drained, mostly during 1960s-1970s [14]. Around the Baltic Sea basin, in total, half of the peatland areas have been drained [13]. Unmanaged peatlands are no longer drained in Finland, but existing ditch networks are maintained. Some studies have reported increases in TN and/or TP concentrations in runoff following ditch network maintenance (DNM) [15,16], but most studies suggest that DNM reduces DOC export [17][18][19]. It has been a common assumption that after the first-time drainage water quality returns to the background level (i.e., the level prior to operation) over time [20]. However, a recent study [21] suggested that N and P concentrations in runoff from old drained peatlands may slowly increase as a function of time. Drainage causes also long-term changes in the quality of organic matter in runoff [22]. There is a clear knowledge gap between the controlled paired catchment studies and monitoring studies concerning response of water quality to long-term effects of drainage and low intensity harvesting such as thinning, small area clear-cutting and partial harvesting. It remains unclear how these effects accumulate over the long time in the catchment scale. This requires more research on C and nutrient export from catchments under normal forest management and, as a reference from unmanaged catchments.
Catchment scale water quality monitoring studies have shown that C, N and P export can vary markedly among both unmanaged [12,23] and managed forests [12,24], implying that part of the variation in water quality can be attributed to factors independent of forest management. Thus, all sources of DOC should be considered in the catchment level management [25]. Many studies have suggested that large part of the variation in DOC and TOC concentration in runoff can be explained by the peatland percentage in the catchment [26][27][28][29][30][31][32]. Other possible drivers for spatial and temporal variation in DOC leaching include temperature [23,33], runoff [23], and atmospheric sulfate and N deposition [34]. Moreover, factors related to net primary production, e.g., stand age [35], may be linked to DOC leaching [36]. Potential factors related to N and P export are also numerous, for instance, tree species, stand volume [23,37], stand age [38,39], site fertility [23], and catchment topography [40,41] have been identified.
Due to the large number of potential explanatory variables and limitations in data availability (either water quality monitoring data or data on catchment properties), comprehensive analyses on the controls of C and nutrient export in boreal forests have remained relatively rare. Furthermore, given that the climate in northern Europe is expected to change towards warmer and wetter, especially for winter and spring, it is increasingly important to understand the combined effects of forest management and the prevailing environmental conditions on water quality. Some climate change induced changes to seasonal variation of nutrient export have already been reported in river catchments [42]. Thus, with advancing climate change, it is essential to gain knowledge also on climatic factors affecting the water quality and nutrient export. This knowledge will be useful, as forestry faces new challenges in finding suitable management practices, even more so as current nutrient load assessments in Finland do not account for climatic drivers [24].
In the present study we analysed water quality data and runoff data from 10 unmanaged and 20 managed forest-dominated headwater catchments in Finland and studied relationships between water quality and export together with land use, tree stand, site and catchment characteristics, and climatic variables. The aims of the study were (1) to determine how much TOC, N and P concentrations and export differ between unmanaged and managed forest catchments and (2) to quantify to what extent climatic vari-ables, catchment characteristics and anthropogenic impacts explain variations in TOC, N and P concentrations and export. The data set covered a large geographical area (>900 km from south to north and >400 km from east to west) including both unmanaged and managed catchments with variable properties obtained from high resolution GIS-data (16 m × 16 m grid). The runoff and water quality measurements were carried out intensively for several years, thus providing a unique data set. The dataset belongs to a monitoring network for water quality in forested catchments and is available at https://metsainfo.luke.fi/fi/vesistokuormitukset, accessed on 2 August 2021.

Study Areas
The study was conducted in 10 unmanaged and 20 managed forested catchments located across Finland ( Figure 1). The catchments represent typical forests and forest management in Finland. In Finland rotation length from planting to harvest is 70-100 years. Forests are typically grown in even aged stands, which are thinned 2-3 times during the rotation before the final harvest. Most abundant species are Scots pine (Pinus sylvestris) and Norway spruce (Picea abies) and birch (Betula sp.). One fourth of the forests are located on drained peatlands where the open ditch drains are typically cleaned 1-2 times per rotation. The area of the catchments ranged from 7 to 12,150 ha. The peatland cover varied from 5 to 48% in unmanaged and 5 to 67% in managed catchments. Most of the managed catchments included drained peatlands. Coniferous species, Norway spruce (Picea abies Karsten) and Scots pine (Pinus sylvestris L.), were dominant in all except one catchment, where Silver birch (Betula pendula Roth.) was the most common tree species (Table S1). The total tree volume varied between 25 and 165 m 3 ha −1 in managed and between 2 and 174 m 3 ha −1 in unmanaged catchments (Table S1, Supplementary Materials). Annual mean air temperature in the catchments ranged from −1.4 to 7.2°C over the monitoring period (2014-2018). The daily air temperatures originate from Finnish Meteorological Institute (FMI) data that are interpolated to 10 × 10 km raster throughout Finland using existing meteorological stations. Mean annual runoff was 249-493 mm in unmanaged catchments and 186-503 mm in managed catchments (data available online at http://kartta.luke.fi/vesidata/, accessed on 2 August 2021). DNM operations were carried out in five catchments during the study period, and the total area of DNM was <5% of the catchment area. Forest regeneration was carried out in 14 catchments and fertilization in 2 catchments. The monitoring catchments were established in forest dominated areas, but small patches of arable areas were found in 12 catchments, where their areal share was mostly <3%. More detailed catchment specific factors, such as tree stand characteristics, land use and management practices are presented in Table S1.

Water Quality Monitoring
At all catchments discharge at the catchment outlet was monitored continuously using V-notch weirs and automatic water level loggers. Discharge monitoring started in all catchments between January 2014 and June 2014. For this study we included data until the end of 2018. On average, 20 water samples (data available: https://metsainfo.luke.fi/ fi/vesistokuormitukset, accessed on 2 August 2021) were collected from each catchment outlet annually. One sample a week was collected during spring high flows, one sample in two weeks intervals during autumn and once per month for rest of the year. More frequent sampling during high flows than during low flows was used to improve the estimates of annual export loads [43]. The water samples were distributed to regional accredited laboratories (Eurofins laboratoris, www.eurofins.fi/environmentandindustry/ environment/water-analysis/, accessed on 2 August 2021) where samples were analysed for pH, concentrations of unfiltered TOC, total nitrogen (TN), ammonium nitrogen (NH 4 -N), nitrate-nitrite nitrogen (NO 2-3 -N), total phosphorus (TP) and phosphate-phosphorus (PO 4 -P). Samples were filtered with Nuclepore polycarbonate filters with 0.4 µm pore size.
where L is annual load (kg ha −1 a −1 ), F is unit conversion factor for annual loads ( 3600 s divided by catchment area (ha)), K is length of hourly discharge timeseries, N is number of water samples, C -i instantaneous concentration at time of water sampling (g L −1 ) and Q p mean discharge (L s −1 ) between water samplings. Annual mean concentrations were computed by dividing L by total annual runoff.

Spatial Data on Catchment Characteristics
Weather data were acquired from Finnish Meteorological Institutes open data base. The mean long-term temperature sum and annual precipitation were calculated for the catchments from interpolated daily weather time series extending from 1980 to 2019. Topographic/terrain maps by the National Land Survey of Finland were used to determine percentages for arable area, and water bodies in the catchments. Drainage intensity (DI) was computed by dividing the total length of ditches (width < 2 m) by the catchment area. Forest stand characteristics within the catchments were obtained from multi-source National Forest Inventory (2013) 16 m × 16 m raster data [44]. The data contains information about site main class (mineral soil, peatlands) and site fertility class from 1 to 6. Each raster cell was characterized by stand volume, tree species, site main class (mineral, bog, mire, open peatland) and site fertility class (1, . . . , 6; fertility decreasing with increasing value). For example, proportion of rich peat was calculated as frequency of raster points belonging peatland and site fertility 1 or 2 divided by raster points in the catchment.

Statistical Analyses
Linear mixed effect models (LME) were used for comparisons of annual concentrations and export between unmanaged and managed areas. This was performed with lmer function from lme4-package [45] in R. Multiple comparisons were performed with multcomp package [46]. Measurement year was set as the random factor, while treatment (unmanaged or managed) was set as fixed factor. Concentrations and exports that did not meet the requirements of a parametric test were log transformed. In addition, we created LMEs by using catchment characteristics and weather data to identify the factors that best describe the annual TOC, TN and TP concentrations and export. In this analysis (Equation (2)), measurement year and catchment id were set as the random effect and the fixed effects were arable area, drainage intensity, peat and mineral soil fertility classes, long term temperature sum, annual precipitation, tree volume and management. We used stepwise selection beginning from Equation (2) and progressing based on model Akaike's Information Criterion (AIC) [47,48] and p-values of fixed effects. The best model was selected based on lowest AIC value, highest R2 and the significance of the fixed effects. For this, we used the drop1 function in R [48]. In all the accepted models, the overall p-value was < 0.05, and explanatory variables had a p-value of <0.1. To avoid multicollinearity, the Variance Inflation Factor (VIF) was calculated and only factors with VIF < 5 were accepted to the model. We chose the potential explanatory variables so that they are not mutually correlated, contain combined ecological and geological information (site fertility in mineral and peat soils, stand volume) and they provide information about management (arable land, drainage, stand volume). The fertility classes also combine ecological and geological information: coarse soils are mostly nutrient poor, while fine soils are nutrient rich. The chosen meteorological factors also carry information of geographical location. The original model was: where y is the dependent factor (TOC, TN or TP export load (kg ha −1 a −1 ) or concentration (mg L −1 or µg L −1 ), a is the model intercept, parameters b 1 -b 9 are regression coefficients of the fixed factors and u 1 and u 2 are the random factors (measurement year, catchment id), Arable is the percentage of arable land in catchment, Drainage is the drainage intensity (m ha −1 ), Rich p is ratio of rich peatlands (fertility classes 1-2,) of total catchment area. Poor m is the ratio of rich mineral soils (fertility classes 1-2 ) of total catchment area. The fertility classes follow the classification presented by Muukkonen and Mäkipää [49]. Temp is the mean long-term temperature sum (1980-2019), Prec is the annual precipitation (mm, 2014-2018). Managed is a dummy variable with value 1 for managed area and 0 for unmanaged area. Tvol is the mean tree volume (m 3 ha −1 ) and Peat p is the percentage of peatlands in the catchment. Catchment characteristics are presented in Table S1 and the correlations of explanatory variables are presented in Figure S1.

Comparison of Unmanaged and Managed Catchments
The comparison of TOC, TN and TP concentrations ( Figure 2) showed that the concentrations in runoff were consistently higher in managed than in unmanaged catchments (p < 0.001). The mean TOC concentration in runoff during the monitoring period for unmanaged catchments was 17.7 ± 2.1 mg L −1 and 28.7 ± 1.6 mg L −1 for managed catchments. The mean TN concentrations were 420 ± 45 µg L −1 and 825 ± 51 µg L −1 for unmanaged and managed catchments, while the corresponding values for TP were 15.3 ± 2.3 µg L −1 and 35.6 ± 3.3 µg L −1 . Furthermore, the export of TOC, TN and TP were higher for managed than unmanaged areas. The mean TOC export (Figure 3) was 56.4 ± 9.6 kg ha −1 a −1 for unmanaged and 79.3 ± 3.3 kg ha −1 a −1 for managed catchments. The difference was statistically significant (p < 0.001). The managed catchments also had significantly higher TN and TP export than unmanaged catchments (p < 0.001). The mean TN export was 1.43 ± 0.2 kg ha −1 a −1 for unmanaged and 2.31 ± 0.2 kg ha −1 a −1 for managed catchments. TP exports were 0.053 ± 0.009 kg ha −1 a −1 and 0.095 ± 0.008 kg ha −1 a −1 for unmanaged and managed catchments, respectively.

Controls of TOC Concentration and Export
The TOC concentration was most affected by drainage, ratio of rich peatlands, longterm temperature sum and peatland percentage. This model explained 61% of the variation in TOC concentrations (AIC = 993). All of the factors had a positive effect on TOC concentration (Table 1). Likewise the TOC export was best explained by ratio of rich peatlands, temperature sum and peatland percentage, all of which increased TOC export (Table 2). This model explained 48% of the variation in TOC export (AIC = 1330). The goodness of fit for the models are presented in Figure 4, while the residual plots are presented in Figure S2.

Controls of TN and TP Concentration and Export
TN concentrations were increased by area of arable lands, drainage, temperature sum and peatland percentage. The model explained 61% of the variation in TN concentrations (AIC = 1845). Arable area, temperature sum and peatland percentage also increased TN export. This model explained 45% of variation in TN export (AIC = 397).
TP concentrations were best predicted by the proportion of arable land, drainage and temperature sum. These together explained 60% of variation in TP concentrations (AIC = 1150). Aforementioned variables increased TP concentrations (Table 1). Arable land, ratio of rich peatlands, temperature sum and tree volume increased TP export, while the ratio of poor mineral soils had a decreasing effect ( Table 2). These explained 54% of variation (AIC = −364).

Carbon and Nutrient Export from Unmanaged and Managed Boreal Catchments
The managed catchments had higher TOC concentrations and export than the unmanaged ones. However, when statistical models were applied, management did not appear among the factors most affecting concentrations and exports. Although seemingly contradictory, this simply implies, that there are other factors (i.e., ratio of rich peatlands, peatland percentage and temperature sum) that have a more profound effect on catchment level export and concentration. For example, the background load of TC from peatlands may be so high, that the effect of management does not necessarily come forward [50] in the stepwise regression.
The observed mean TOC exports (mean 56.5 kg ha −1 a −1 ) in unmanaged areas were close to previously reported values (mean 62 kg ha −1 a −1 ) based on data collected from 21 to 28 pristine Finnish catchments [23,24,28]. Two of these studies had more catchments than the present study [23,28], but they had no observations from catchments located in the northernmost parts of Finland where the TOC export is typically low. This may explain the slightly lower TOC export in unmanaged catchments in our study compared to the previous studies. The third study had observations also from northern Finland and very long time series   [24], which could as well explain the differences. In managed catchments, the mean TOC export obtained in our study (79.4 kg ha −1 a −1 ) is within the range reported in other boreal catchments (23-148 kg ha −1 a −1 ) [4,51].
The statistical comparisons further indicate that TN and TP concentrations and export from managed boreal catchments are higher than those from unmanaged catchments (Figures 2 and 3). TN export from managed catchments (2.3 kg ha −1 a −1 ) found in our study is close to the mean value of 1.8 kg ha −1 a −1 reported previously [10]. Furthermore, for the unmanaged catchments, the range and mean of TN and TP export found in the present study ( Figure 3) agree with the previously published values [23,28].
Given the good geographical coverage of our catchments ( Figure 1) and the consistent water sampling strategy applied (see Section 2.2), we argue that the data are hitherto the best available for estimation of TOC, N and P export from the variety of managed and unmanaged forested catchments in Finland. However, we acknowledge that the time series presented in the study are still relatively short (5 years) and although the data is suitable for comparative purposes, more data needs to be collected and analysed to draw conclusions on the effects of forestry operations to C and nutrient export.

Controls on TOC Export
The LME showed that 61% of TOC concentration was explained by ratio of rich peatlands, long-term temperature sum, drainage and peatland percentage, while 48% of TOC export was explained by the rich peatlands, peatland percentage and temperature sum. The rich peatlands increased the TOC concentrations and export. Rich peatlands often are characterized by Carex peat, which releases more DOC than Sphagnum peat typical in poorer sites [52,53]. Fertile peatlands have been also found to be hotspots for TOC export in subarctic regions [54] and continental climates [32]. A positive correlation between peat coverage and TOC (or DOC) export is well known in many different environments [26][27][28][29]31,32], which was also apparent in our TOC models. Opposed to peatlands, mineral soils tend to retain DOC and therefore DOC concentrations in streams do not necessarily reflect the DOC release from organic matter [55].
Furthermore, TOC export and concentrations increased with increasing temperature sum. It has been found that TOC concentration in runoff increases with increasing temperature in undisturbed forests [23]. A positive correlation has also been found with TOC concentration in runoff and mean air temperature in peat extraction areas [56]. However, there are contrasting studies on the overall effect of temperature on decomposition and carbon balance in peat soils [57,58]. Nonetheless, warmer climate accelerates decomposition of soil organic matter [59,60], leading to higher release and potential transport of TOC. Supporting this view, previous studies have suggested that the increasing TOC trends have been caused, at least in part, by simultaneous increase in air and stream water temperature [33,61]. It is possible that the ongoing global warming will lead to increasing TOC exports in boreal catchments, as the climate in northern latitudes is predicted to experience the most drastic increases in temperature (IPCC in Climate Change, 2013). There are indications that the terrestrial boreal region is already reacting to ephemeral effects of global warming [62].
The TOC concentrations were increased by drainage intensity as also found by other studies [24]. Drainage is linked to hydrological conditions in the catchment and thus affect leaching [63], resulting in higher leaching of TOC. In less intensively drained areas TOC is not leached readily because organic matter decomposition is slowed down by higher water table [64]. While in more intensively drained areas water flow path and water residence time are shorter [65] thus increasing leaching, in less drained sites these are longer, allowing more efficient retention [66] and biodegradation of TOC during the transport [52]. However, some previous studies show no increase in organic carbon release due to drainage [15,51,67]. Given that DNM also may reduce DOC leaching [17][18][19] and that it is usually performed simultaneously to clear-cutting on drained areas (to compensate decreasing transpiration), DNM operations may counterbalance the potential increase in TOC release caused by clear-cutting [7]. For example, Åström et al. noted that after water table lowering the residence time of water in upper organic layers decreases, leading to decreased leaching of soluble matter [17]. At the same time humic substances may be retained in the mineral soil part of the ditches [17].

Controls of Nitrogen and Phosphorus Export
Arable land increased the concentrations and exports of TN and TP. Agriculture is the most significant source of N to surface waters, especially in the southern and western parts of Finland [68]. The higher N and P export and concentrations connected to arable lands are most commonly related to fertilization [10,13,69]. The results indicate that even a small proportion of arable land in catchment area, where management operations (fertilization, soil preparation etc.) are done yearly, can markedly increase TN and TP exports from the whole catchment [70]. Furthermore, temperature sum had an increasing effect on TN and TP concentrations and exports, likely because higher temperature increases biomass production, organic matter decomposition and the length of growing season. Of the biomass related factors in the model, only TP export was affected (decreased) by the tree volume, probably caused by P uptake by trees. Previously, it has been shown that forest harvesting increases labile P sources in soil [71]. Further, temperature sum is not only linked to the biomass production and growing season length, but also to the geographical location of the catchment.
Ratio of rich peatlands had an increasing effect on TP export, while TN concentrations were increased by the total peatland coverage. Fertile peatlands contain large nutrient stocks which are potentially subjected to leaching [72]. Peat accumulation on fens has been found to control the long-term P storage [73], whereas nutrient poor peatlands are prone to P leaching [74]. Iron and aluminium in peat can decrease P leaching [75], and poor peatlands have very low concentrations of Fe and Al, and therefore, may actually have higher P leaching than fertile sites [76]. Positive correlation between soil fertility and N leaching have been previously shown to occur among unmanaged boreal forests [23] and among drained peatland forests [15]. Other studies have suggested that peatland coverage per se is important for determining N and P losses from forested catchments [61,77,78]. In addition, increasing ratio of poor mineral soils decrease TP export. P stocks in poor mineral soils are small, and soils are typically coarse textured [79] where water dominantly moves downwards and P is effectively retained in the B-horizon of podzol profile [13]. Next, we suggest that, when appropriate data are available, different peatland types and their fertility could be also analysed separately to improve the assessment of TN and TP export from boreal forests. It also seems that peatland and mineral soil fertility class should not be overlooked in models and they should be considered together with the intensity of management operations.
Though statistical comparison showed that TN and TP concentrations and exports were higher in managed than unmanaged catchments, this was not observed in the mixed effect models, even though forest management practices have been found, e.g., to increase N leaching [80,81]. In boreal catchments the organic fractions of N and P dominate the export, and the export of the organic fractions is mostly governed by natural sources, such as peatlands [50]. Especially in peatland dominated catchments, the management covering a small share of the catchment at the time can easily be masked behind the natural background exports.
Drainage intensity increased both TN and TP concentrations, following previous findings [24]. The N cycle is strongly affected by management such as drainage and clearcutting that affect through decreased plant nutrient uptake, enhanced mineralization and leaching, and increased denitrification [82,83]. Increased N leaching may last 10-15 years after clear-cutting [9,84]. Likewise, it has been discovered that TN and TP concentrations are much higher among drained peat-dominated catchments in comparison with unmanaged (undrained) peatlands [21]. We suppose that the effect of drainage on N and P concentrations is related to increased nutrient mineralization and aerobic decomposition of the organic surface layers [64]. Moreover, the results of the present study are consistent with Nieminen et al. [21,85] in suggesting that peatland drainage is at least in part responsible for the increased TN and TP export from boreal forests.

Conclusions
This study covered overall 30 catchments across Finland, producing a geographically encompassing survey of water quality in both unmanaged and managed boreal forested catchments. Responsible forest management requires knowledge of how and why water quality varies between different catchments. On average TN and TP exports were about two times higher in managed than in unmanaged catchments. While the biogeochemical mechanisms underlying C and nutrient export cannot be resolved with the catchment level monitoring, we were able to identify the dominant natural and anthropogenic factors contributing to water quality. The results should be applicable to forest-dominated catchments in similar geologic and climatic conditions. In addition, the increasing effect of temperature sum on nutrient export and concentrations shows also that global warming may present us with new challenges in water protection in the future. Our models show that including soil and peatland fertility classes may provide additional information of the controls of nutrient export. We conclude that while it is difficult to isolate causal mechanisms and their interactions from the regression approach, the relatively simple and readily available catchment characteristics and climate variables identified in this study can be very helpful in assessing C and nutrient export from different boreal forested catchments. The obtained models provide information on the most important factors affecting the concentrations and exports in boreal catchment areas. Thus, these models may be utilized to identify risk areas and to estimate the dependence of export of catchment characteristics and location.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/w13172363/S1, Table S1: Catchment characteristics of unmanaged and managed boreal catchments used in the study. Figure S1: The correlations of fixed factors used in the models, Figure

Data Availability Statement:
The water quality monitoring data are updated on an annual basis and are freely available at http://kartta.luke.fi/vesidata/, accessed on 2 August 2021.

Acknowledgments:
The authors would like to thank Kristian Karlsson for support and advice on data analysis. The GIS data used in the study were produced and provided by the Natural Resources Institute Finland, Finnish Meteorological Institute, Geological Survey of Finland, and the National Land Survey of Finland.

Conflicts of Interest:
The authors declare no conflict of interest.