We analysed long-term (1992–2020) changes in fertiliser use, wastewater treatment, and river water nutrient status in Estonia (N-E Europe) in the context of changing socio-economic situations and legislation. We hypothesised that improved regulation of fertiliser usage and wastewater treatment are reflected as declining riverine nutrient concentrations, with the largest relative improvements occurring in catchments with initially high proportions of point source loading. We analysed nutrient dynamics in 16 rivers differing by catchment land use, population and livestock densities. Data on fertiliser use and wastewater treatment originated from the Statistics Estonia database, and riverine nutrient concentrations from the State Environmental Monitoring Database. We clustered the rivers by their catchment properties and analysed trends in their nutrient status. Point source nutrient loading reductions explained most of the decline in riverine nutrient concentrations, whereas application of mineral fertilisers has increased, hindering efforts to reach water quality and nutrient load targets set by the EU Water Framework Directive and the Baltic Sea Action Plan. Highest nitrogen concentrations and strongest increasing trends were found in rivers within the Nitrate Vulnerable Zone, indicating violation of the EU Nitrates Directive. To comply with these directives, resource managers must address non-point source nutrient loading from river watersheds.

  • Drop in point source loading explained the decline in riverine nitrogen (N) and phosphorus (P) since 1994.

  • Fertiliser and wastewater management measures failed short to meet the water quality and nutrient load targets set by the EU Water Framework Directive.

  • Highest N concentrations and strongest increasing trends were found in rivers within the nitrate vulnerable zone violating the EU Nitrates Directive.

Graphical Abstract

Graphical Abstract
Graphical Abstract

Anthropogenic nutrient loading in aquatic systems is a combination of discharges from point and non-point sources. Point sources, such as industries and wastewater treatment plants, are easier to identify and implement management approaches to reduce nutrient losses. Runoff from agricultural areas is the major cause of non-point source nutrient pollution worldwide and has long been considered as a major driver of eutrophication in downstream rivers and lakes (Barcala et al. 2020; Xia et al. 2020). In Europe, reduced loading from point sources has shifted agriculture into the key role in nutrient losses from rural areas to surface waters over recent decades (European Environment Agency 2005).

By 2050, the global N pollution level is projected to rise by 50% of that in 2010, with agriculture accounting for 60% of this increase; nearly half of N fertiliser applied to fields is lost via gaseous emissions or discharges to water bodies (Martínez-Dalmau et al. 2021). To alleviate eutrophication, it is often necessary to reduce the inputs of both N and P (Paerl et al. 2016).

Different dynamics in N and P concentrations have been observed in Europe. For example, P concentrations increased until 2004 in the Vltava River, Slapy reservoir (Czech Republic), then decreased, while N concentration increased throughout the period (1960–2019; Kopáček et al. 2021). N concentrations and losses from agricultural catchments decreased in Denmark and Sweden from 1988 to 2011, likely due to implementation of agro-environmental measures (Stålnacke et al. 2014). However, increasing N concentrations and catchment losses were detected in the Baltic countries, coinciding with intensification of agriculture in the Baltic States (Stålnacke et al. 2014). These results indicated that targeted strategies towards reducing N losses from agricultural land would improve surface water quality.

In the Baltic Sea catchment, a major part of the nutrient load originates from diffuse sources, mainly agriculture. Point sources now comprise 4–5% of nutrient loads to the Baltic Sea and has decreased substantially in recent decades (HELCOM 2018). Riverine discharges contribute 70.3% of total N (TN) and 94.8% of total P (TP) inputs to the Baltic Sea (HELCOM 2018).

Major socio-economic changes took place in Estonia in the early 1990s. The re-establishment of Estonia's independence closed a large part of the Russian markets to Estonian production, and agricultural production deceased. At the same time, lands expropriated in the 1940s were returned to their former owners and their descendants. As many of them lacked the skills and opportunities to cultivate the land, large areas of former agricultural land lied fallow, and few farmers could afford to buy fertilisers. These changes led to smaller meadows and pastures along riverbanks and lake shores being overgrown with reeds and brushwood. Estonia implemented the Water Act in 1994 and adopted European directives to replace the soviet-era laws (Supplementary material, Table 1). The spreading of fertilisers on frozen fields and snow was banned in 2001. After joining the EU, water consumption taxes increased, which reduced water extraction and the amount of wastewater (Supplementary material, Figure 3). With the support of EU funds, new water treatment plants were constructed, and existing plants were reconstructed. All of these measures were expected to improve water quality in the rivers.

Table 1

Sampling sites in rivers sorted by catchment size

River/sampling siteSite numberWFD typeCatchment type cluster numberCatchment area, km2Number of water samplesPopulation density, in h/km2Forest, %Wetland, %Arable, %Pasture and grassland, %
Emajõgi Kavastu 18 3C 8,539 348 21.9 42 11 19 23 
Väike-Emajõgi mouth 16 3C 1,270 256 24.1 46 13 32 
Võhandu mouth 3C 1,144 347 29.8 42 10 37 
Pirita mouth 34 2C 794 287 83.4 44 15 16 21 
Keila mouth 32 2C 682 345 34.1 38 12 24 22 
Reiu Lähkma 52 2B 548 177 12.2 62 17 10 
Kunda mouth 38 2C 528 344 15.5 49 13 17 20 
Vihterpalu mouth 40 2B 474 233 57 26 
Valgejõgi mouth 36 2C 453 345 27.3 47 21 13 16 
Avijõgi mouth 10 2C 366 344 8.2 58 14 13 14 
Pühajõgi mouth 47 2C 196 348 183.1 48 10 19 13 
Saarjõgi Kaansoo 53 2B 191 145 2.9 61 17 13 
Velise Valgu 43 2B 135 160 2.1 66 13 11 10 
Vodja Vodjaa 46 1C 52 151 10.4 22 19 35 24 
Preedi Varangua 24 1C 34.8 166 8.4 35 34 22 
Oostriku Oostrikua 25 1C 29.7 143 12.1 23 53 17 
River/sampling siteSite numberWFD typeCatchment type cluster numberCatchment area, km2Number of water samplesPopulation density, in h/km2Forest, %Wetland, %Arable, %Pasture and grassland, %
Emajõgi Kavastu 18 3C 8,539 348 21.9 42 11 19 23 
Väike-Emajõgi mouth 16 3C 1,270 256 24.1 46 13 32 
Võhandu mouth 3C 1,144 347 29.8 42 10 37 
Pirita mouth 34 2C 794 287 83.4 44 15 16 21 
Keila mouth 32 2C 682 345 34.1 38 12 24 22 
Reiu Lähkma 52 2B 548 177 12.2 62 17 10 
Kunda mouth 38 2C 528 344 15.5 49 13 17 20 
Vihterpalu mouth 40 2B 474 233 57 26 
Valgejõgi mouth 36 2C 453 345 27.3 47 21 13 16 
Avijõgi mouth 10 2C 366 344 8.2 58 14 13 14 
Pühajõgi mouth 47 2C 196 348 183.1 48 10 19 13 
Saarjõgi Kaansoo 53 2B 191 145 2.9 61 17 13 
Velise Valgu 43 2B 135 160 2.1 66 13 11 10 
Vodja Vodjaa 46 1C 52 151 10.4 22 19 35 24 
Preedi Varangua 24 1C 34.8 166 8.4 35 34 22 
Oostriku Oostrikua 25 1C 29.7 143 12.1 23 53 17 

Note: Site numbers according to the Environmental Monitoring Information System (KESE 2021). Coding of Water Framework Directive (WFD) type names: 1 – catchment area < 100 km2; 2 – catchment area 100–1,000 km2; 3 – catchment area > 1,000 km2; B – brown waters, chemical oxygen demand > 25 mg O/l; C – clear waters, chemical oxygen demand < 25 mg O/l.

aNitrate vulnerable zone (NVZ).

After the economic downturn in the 1990 and 2000s, agriculture recovered in Estonia due to EU subsidies (although smallest within the EU). These investments encouraged agricultural redevelopment, and fertilisers were again affordable for farmers. Concentration of agriculture by large landowners enabled new agronomic techniques (e.g., catch crops, granular fertilisers) and smart dosing to avoid over-fertilisation.

Here, we analyse trends in long-term changes in wastewater treatment, N and P fertiliser use, and river water nutrient status in the context of changing legislation and socio-economic status. We hypothesise that (1) more precise regulation of fertiliser use and improved wastewater treatment have led to declining TN and TP concentrations, and (2) the largest water quality improvements occurred in catchments with a high proportion of point source pollution.

Estonia (45,339 km2) has approximately 2,000 km2 of surface waters (Figure 1). Slightly more than half of the land is coniferous and mixed forests, with agricultural land (23%) and wetlands (22%) comprising much of the rest. In central Estonia, large limestone areas are characterised by karst phenomena, where the thickness of glacial drifts can be insignificant, and the groundwater is unprotected. The karst area largely overlaps with the nitrate vulnerable zone (NVZ), which includes areas designated as being at risk from agricultural nitrate pollution (Figure 1). The NVZ in Estonia was determined considering soil and ground conditions, ground and surface water vulnerability, and intensity of agriculture (Water Act 2019) and covers 3,250 km2, or 7.5% of the land area.
Figure 1

Location map of the study area and water quality monitoring stations in the rivers. Site numbers according to the Environmental Monitoring Information System (KESE 2021) (n = 16). Nitrate vulnerable zone surrounded with dashed line.

Figure 1

Location map of the study area and water quality monitoring stations in the rivers. Site numbers according to the Environmental Monitoring Information System (KESE 2021) (n = 16). Nitrate vulnerable zone surrounded with dashed line.

Close modal

Long-term (1981–2010) annual average air temperature is between 4.9 and 7.1 °C (Estonian Weather Service 2021). Snow cover lasts from 61 to 130 days, with high spatial variation (Tarand et al. 2013; Viru & Jaagus 2020). Estonia is sparsely populated, with 30.6 inhabitants per km2 (2021).

There are approximately 31,000 km of natural running waters, and their density is 0.72 km km−2. The total annual runoff is about 12 km3, accounting for about 40% of annual precipitation (Timm et al. 2019).

Rivers

TN and TP concentrations were analysed in 16 rivers, which were selected based on a sampling frequency of at least four times per year (Table 1; Figure 1). Six rivers flow directly into the Baltic Sea, three into Lake Peipsi, one into Lake Võrtsjärv, and the other six are tributaries of different rivers. Three sampling sites (Nos 24, 25, and 46) were situated in the NVZ. Details on catchment land use of these sites are provided in Iital et al. (2010). According to the EU Water Framework Directive (Council Directive 2000/60/EC 2000, WFD) typology in Estonia, based on catchment size and water colour (humic matter content measured as permanganate oxygen demand), the selected sites belonged to four different WFD types (Table 1).

TN and TP concentrations in river water were obtained from the Environmental Monitoring Information System (KESE 2021). Historical data on daily flow rates of rivers for the period 1992–2020 were obtained from the Estonian Weather Service (2021).

Data collection

Data on agricultural activities were obtained from the Statistics Estonia database (2021). Data on N and P in livestock manure (PM0646: Nitrogen and phosphorus in livestock manure) were available only for 6 years (2014–2019) and are not directly comparable with previous data from 1992 to 2011 (KK19: Nitrogen and phosphorus brought into the soil with organic fertilisers) because the methodology for measurement changed. No data were available for 2012 and 2013 (Supplementary material, Figure 1). Data for the numbers of livestock and poultry were available from 1992 to 2020 (PM09: Livestock and poultry). In accordance with the Ministry of Rural Affairs Regulation No. 71 (2014; Estimated nutrient values for different types of manure, methodology for calculating manure storage capacity and conversion factors for livestock units; Riigiteataja I: 16.07.2014, 8; Appendix 9, Coefficients for converting animals into livestock units based on animal faeces), the amount of manure produced by livestock and poultry was calculated. The number of livestock units (LSUs) for each year back to 1992 was also calculated (Supplementary material, Figure 2). Consistent with the available 6 years of data, one LSU produced an average of 120 kg N and 23 kg P per year. These numbers were multiplied by the LSU numbers for the years 1992–2013 to generate a homogenous time series for N and P applications via manure. Data on application of N and P mineral fertilisers, use of water, and N and P discharged with waste waters for 1992–2020 were obtained from the Statistics Estonia database (PM065: Use of mineral fertilisers for the production; KK47: Water use by field of water use; KK23: Waste water treatment; KK25: Pollution load to surface water bodies).

Data analysis

River sites differed by catchment characteristics (area, population density, percentages of land use categories; see Table 1), so clustering was used to typify the catchments. In contrast to the WFD typology, which is based on natural features (catchment size, water colour), our clustering approach considered long-term human impacts reflected in population density and land use categories in the catchment. For cluster analysis, catchment characteristics for all 16 river sites were standardised by subtracting the mean of each variable and dividing the residuals by the standard deviation. Using the k-means clustering procedure in STATISTICA 64 (Version 13, TIBCO Software Inc. 2017), the sites were grouped into three clusters. The number of clusters was optimised using the Elbow method (Kodinariya & Makwana 2013). K-mean clustering was selected because it shows the role of different factors in defining the clusters. Differences between clusters by individual factors were determined using the Fischer LSD test at 95% confidence.

Trends in river time series were analysed with the Mann–Kendall non-parametric test, and step changes in means were evaluated using the Distribution-Free CUSUM non-parametric test. A step change identifies an inflection point in a time series characterised by the largest difference between the means in two parts of the record. Both tests were run using trend/change detection software (TREND) provided in the e-Water Toolkit (https://ewater.org.au/products/ewater-toolkit/, accessed 20 December 2021).

Changes in agriculture and fertiliser application

Both the extension and intensity of Estonian agriculture declined after re-establishing its sovereignty in 1992 (Supplementary material, Table 2; Figure 2). The land area fertilised with mineral and organic fertilisers, as well as the use of fertilisers per hectare, decreased more than two-fold. The application of mineral fertilisers decreased initially, reaching a minimum in 1996, because farmers could not afford agrochemicals (Lehtsaar & Jullinen 1998). This situation promoted the development of organic farming (Zobena 1998), which is still popular in the Baltic countries and exemplifies how societal practices can be facilitated by economic situations. Altogether, the use of fertilisers decreased more than 3.5-fold. Numbers of livestock also decreased, and the amount of manure declined accordingly (Supplementary material, Figure 1), resulting in less N and P being applied to fields. This trend reversed in the early 2000s, and more than 20,000 t of N and 4,000 t of P are applied with organic fertilisers annually (Figure 2; Supplementary material, Table 2).
Table 2

Concentrations of TN, TP, and the TN/TP ratio in three river clusters (mean ± standard error)

ClusterTN (mg L−1)TP (μg L−1)TN/TP
2.18 ± 0.03 36 ± 1 83 ± 2 
2.01 ± 0.03 98 ± 3 33 ± 1 
3.50 ± 0.05 54 ± 2 199 ± 12 
ClusterTN (mg L−1)TP (μg L−1)TN/TP
2.18 ± 0.03 36 ± 1 83 ± 2 
2.01 ± 0.03 98 ± 3 33 ± 1 
3.50 ± 0.05 54 ± 2 199 ± 12 
Figure 2

Use of N and P in mineral (N-min and P-min) and organic fertilisers (N-org and P-org) in Estonia (1992–2020).

Figure 2

Use of N and P in mineral (N-min and P-min) and organic fertilisers (N-org and P-org) in Estonia (1992–2020).

Close modal

After Estonia joined the EU in 2004, the amount of arable land and utilised agricultural area, as well as the use of mineral fertilisers, increased with the aid of agricultural subsidies. There was a brief decrease in fertiliser usage after the global economic recession in 2008. By 2020, mineral N and P fertiliser applications had increased by 2.5 and 3.7 times, respectively, over 1996 values. In recent years, mineral N applications (around 40,000 t year−1) exceeded that of organic N. Mineral and organic P were applied almost equally at around 4,000 t year−1 in recent years, but mineral P fertiliser applications have doubled since 2009 (Figure 2). In 2019, the amount of LSU in Estonia was 180,000 (Supplementary material, Figure 2), the area fertilised with mineral and/or organic fertilisers was 6,198 km2, and 63,119 t of N and 8,269 t of P were applied on agricultural land (Supplementary material, Table 2).

Nutrient discharges with wastewater

From 1994 to 2020, TN and TP discharges decreased three- and six-fold, respectively, with a step change (p< 0.01) in 2004 (Figure 3) due to a rapid increase in wastewater treatment plant construction in the early 2000s (Supplementary material, Figure 4). Since 2013, N and P discharges have stabilised at 1,200–1,300 and 60 t year−1, respectively. A larger decrease in P compared to N in wastewater has led to a distinct increase in the N/P ratio (Figure 3). Although the total water usage, as well as industrial and domestic water use, decreased by half during the 1990s (Supplementary material, Figure 3), the volume of wastewater did not decrease substantially during the study period (Figure 3). Water usage in Estonia is primarily cooling water in thermal power plants, followed by industrial and domestic uses (Supplementary material, Figure 3).
Figure 3

Nutrient load to Estonian surface water bodies with discharged wastewater (WW); arrow indicates the step change from 2004 to 2005 (p < 0.01) occurring in both nutrient series.

Figure 3

Nutrient load to Estonian surface water bodies with discharged wastewater (WW); arrow indicates the step change from 2004 to 2005 (p < 0.01) occurring in both nutrient series.

Close modal

Natural vs. anthropogenic river types

The WFD-compliant river typology in Estonia in mainly based on two natural characteristics: catchment size and permanganate oxygen demand as a proxy for the content of humic substances or water colour (Pinnaveekogumite 2020). However, water quality and nutrient dynamics are largely controlled by anthropogenic impacts, such as long-term population density and land use in the catchment. For the study period, both parameters can be considered static and used for typifying the rivers. We entered these catchment parameters into cluster analysis to categorise the selected river sections by their degree of anthropogenic influence (although catchment size was also included). Three river clusters were identified, differing significantly by percentages of forest and non-cultivated agricultural land in their catchment areas (Figure 4). Additional variables, such as area, population density, percentages of wetlands and cultivated agricultural lands contributed to distinguishing clusters from each other although not all differences were significant.
Figure 4

Standardised mean values of river characteristics used for clustering. For each variable, different letters denote differences between groups determined using the Fischer LSD test at 95% confidence. Variables are in arbitrary order and connected with lines to better visualise the contrasting patterns.

Figure 4

Standardised mean values of river characteristics used for clustering. For each variable, different letters denote differences between groups determined using the Fischer LSD test at 95% confidence. Variables are in arbitrary order and connected with lines to better visualise the contrasting patterns.

Close modal

Cluster 1 consisted of seven WFD type 2 rivers characterised by low population density, highest percentages of wetlands and forests, and lowest percentages of agricultural lands in their catchments (Table 1). All four humic water rivers were included in Cluster 1. Cluster 2 consisted of five rivers, including all three WFD type 3 rivers. These rivers exhibited highest average population density and percentage of pasture and grasslands, and lowest wetland percentage, in their catchments, whereas proportions of forests and cultivated lands were average. Cluster 3 included the three WFD type 1 rivers of the NVZ, plus R. Keila. This cluster had the highest proportion of cultivated lands and the lowest proportion of forests in their catchments (Figure 4).

Changes in river water quality

Mean concentrations of TN and TP, and the TN/TP ratio, differed (p< 0.01) between the clusters (Table 2). The least impacted rivers (Cluster 1) had the lowest TP, rivers in densely populated watersheds (Cluster 2) had the lowest TN and TN/TP ratio, but the highest TP, whereas rivers in the most agricultural watersheds (Cluster 3), including the NVZ rivers, had the highest TN and TN/TP ratio.

These findings are in line with the different sources and pathways of these elements in the catchments. N and P are synanthropic elements, as their concentrations in surface waters tend to increase with population densities in the catchment (Chen et al. 2016). The capacity of agricultural soils to retain excess N is low, and leaching losses from common grain-production systems typically range from 10 to 30% of total N inputs (Meisinger & Delgado 2002). P forms complexes with iron, aluminium, calcium, and organic matter and is well-bound to a variety of soil types. Leaching of P to surface waters increases only after exceeding soil saturation levels (Djodjic et al. 2004). P loading may also be associated with discharges of untreated wastewaters (Edwards & Withers 2007) and, hence, can be considerably reduced by wastewater treatment.

TP concentrations declined in all river clusters from 1992 to 2020, especially in Cluster 2, where TN concentrations also declined (Figure 5). In the other two clusters, TN increased marginally (Cluster 1, p = 0.152; Cluster 3, p= 0.052), coinciding with increased mineral N fertiliser use. As a result, TN/TP ratio increased in most cases, and most strongly in Clusters 2 and 3. River discharges showed no trend in any of the clusters.
Figure 5

Mann–Kendall trends in river nutrient concentrations (TP and TN), nutrient ratio (N/P), and discharge (Q) from 1992 to 2020. Dashed red lines denote strength of the trends, with the NS range indicating no trend. Please refer to the online version of this paper to see this figure in colour: http://dx.doi.org/10.2166/nh.2022.057.

Figure 5

Mann–Kendall trends in river nutrient concentrations (TP and TN), nutrient ratio (N/P), and discharge (Q) from 1992 to 2020. Dashed red lines denote strength of the trends, with the NS range indicating no trend. Please refer to the online version of this paper to see this figure in colour: http://dx.doi.org/10.2166/nh.2022.057.

Close modal
Step changes of nutrient concentrations occurred in several rivers (Figure 6). Most of the downward step changes in both TN and TP occurred in Cluster 2 between 2003 and 2010, whereas upward step changes were frequently observed in Cluster 3. No significant step change was found for the discharge of any river (not shown).
Figure 6

Step changes (p < 0.05) in nutrient concentrations and their ratios in rivers from 1992 to 2020. Number on the bars denote cluster numbers of the individual rivers, and two numbers on one bar mean two rivers. Height of the bar denotes strength of the change.

Figure 6

Step changes (p < 0.05) in nutrient concentrations and their ratios in rivers from 1992 to 2020. Number on the bars denote cluster numbers of the individual rivers, and two numbers on one bar mean two rivers. Height of the bar denotes strength of the change.

Close modal

Observed changes in river water quality reflect sensible outcomes of measures applied in the watersheds. Some improvements in river water quality were evident, with major changes observed in rivers draining catchments with high population and livestock densities (Cluster 2). Livestock densities were indicated by the percentage of non-cultivated agricultural lands (i.e., pastures and grasslands). Both high population and livestock densities suggest a major role for point sources, such as runoff from municipal settlements and cattle farms (Kiedrzyńska et al. 2014), which are theoretically easier to mitigate (Bowes et al. 2008). In this cluster, TP declined, but TN step changes also occurred in the mid-2000s, coinciding with construction of wastewater treatment plants (Supplementary material, Figure 4).

The use of organic fertilisers did not change remarkably during the study period, while that of mineral N fertilisers has increased since the late 2000s (Figure 2). Simultaneously, TN concentrations increased in river water (Figure 6), especially in Cluster 3, with the highest proportion of cultivated lands and lowest proportion of forests in their catchments, including rivers in the NVZ (Figure 4).

Despite these efforts in fertiliser and wastewater management, the main objective of the WFD, to achieve good water quality status by 2015, was not met (Tsakiris 2015). HELCOM (2021) pointed out that nutrient loading, dominated by mineral fertilisers inputs, still exceeds the goals of the Baltic Sea Action Plan, and a further load reduction has been agreed upon (ICES Advice 2019). From 2000 to 2010, the use of mineral fertilisers around the Baltic Sea has increased by 68% for N and by 70% for P (Hong et al. 2017), despite efforts of conservationists and governments to reduce impacts on water resources. No decreasing trends have been observed, regardless of new agro-requirements, such as keeping field books, fertilisation plans, compliance with fertilisation restrictions in the NVZ, requirement for winter vegetation cover, etc.

In a study of 53 Estonian rivers over 1984–2006, Iital et al. (2010) found declining trends for TN in 18 rivers and for TP in 13 rivers. Increasing trends occurred in two rivers for TN and six rivers for TP. Annual N losses from agricultural land ranged from 10 to 40 kg ha−1 (Iital et al. 2014), but variability was greater for P, ranging from 0.09 to 7.5 kg ha−1 (Pengerud et al. 2015). Bechmann et al. (2014) studied trends in N balance and its impact on N concentrations in streams. In four of 14 long-term time series in streams, Bechmann et al. (2014) found positive relationships between annual N balance in the catchment and TN concentrations.

The above studies apply to mineral fertilisers, but the role of organic fertilisers has not been studied as intensively (see Supplementary material, background information). However, to reverse large-scale eutrophication of water bodies, both N and P from mineral and organic fertilisers and wastewaters need to be reduced (Haque 2021; Löw et al. 2021; Sabo et al. 2021). Nutrient losses could also be reduced by partial replacement of mineral fertilisers with effective manure application (Hong et al. 2017). In Estonia, the Water Act allows 175 kg N and 25 kg P to be applied per hectare; for manure, the resulting N/P ratio is 7. The actual N/P ratio is 5, so it is recommended to calculate manure applications according to P requirements to avoid P over-application. In practice, the calculations are often based on N. P over-application may be acceptable in some cases, but it should be strictly monitored to avoid P losses to surface waters. According to the European Nitrate Directive (Council Directive 91/676/EEC 1991), the maximum allowed amount of N in mineral fertilisers in the NVZ is 100 kg per hectare of the area under cultivation, and the use of organic fertilisers is restricted by allowing not more than 1.5 LSU per hectare of the area under cultivation. The threshold NO3 concentration in ground and surface water is 50 mg L−1.

The collapse of soviet-era agriculture in Estonia in the early 1990s contributed to declines in N and P concentrations in rivers. Adoption of EU legislation on water use and protection created a good basis for further improvement of river water quality. Success was achieved in reducing point source pollution from large farms and municipalities, supporting our hypothesis that regulatory measures would improve water quality. However, despite more strict governance on fertiliser use, the decline of nutrients in river water ceased after Estonia joined the EU in 2004, as the financial means to buy fertilisers improved. From 1996 to 2020, the use of mineral N and P fertilisers has increased in Estonia by factors of 2.5 and 3.7, respectively. A particular concern is increased N concentrations in rivers within the NVZ. Among three river clusters, changes were minimal in Cluster 1, representing rivers least impacted by humans, but median TP declined in this cluster. The largest changes were observed in rivers draining catchments with high population and livestock densities; i.e., with a high proportion of point source pollution, thus supporting our second hypothesis. The mitigation of point source pollution resulted in declines in both N and P concentrations, whereas step changes coincided with or followed construction and upgrades of new and existing wastewater treatment plants, respectively, in the early 2000s. We found highest N concentrations and strongest increasing trends in Cluster 3, comprising rivers draining the most agricultural catchments, including those within the NVZ. These results indicate that the rules set by the EU Nitrates Directive (Council Directive 91/676/EEC 1991) to protect surface and groundwater quality in karst areas have been insufficient to reduce nutrient concentrations in surface waters within these watersheds.

This study was funded by the Estonian Research Council grants PRG709 and PRG1167, by the European Regional Development Fund through the Estonian University of Life Sciences ASTRA project ‘Value-chain based bio-economy’. This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 951963. The Estonian Ministry of Environment and the Estonian Environment Agency supported data collection in the national monitoring programme.

T.N. and P.N. conceptualised the study; S.V. and P.N. studied methodology and did formal analysis; S.V. did investigation and data curation, and wrote original draft; T.N., P.N., M.M. and M.T. wrote, reviewed and edited the article; S.V., M.T. and P.N. visualised the study; T.N. was involved in project administration. All authors have read and agreed to the published version of the manuscript.

Data cannot be made publicly available; readers should contact the corresponding author for details.

The authors declare there is no conflict.

Barcala
V.
,
Rozemeijer
J.
,
Osté
L.
,
der Grift
B. V.
,
Gerner
L.
&
Behrends
T.
2020
Processes controlling the flux of legacy phosphorus to surface waters at the farm scale
.
Environmental Research Letters
16
,
015003
.
doi:10.1088/1748-9326/abcdd4
.
Bechmann
M.
,
Blicher-Mathiesen
G.
,
Kyllmar
K.
,
Iital
A.
,
Lagzdins
A.
&
Salo
T.
2014
Nitrogen application, balances and their effect on water quality in small catchments in the Nordic–Baltic countries
.
Agriculture, Ecosystems & Environment
198
,
104
113
.
doi:10.1016/j.agee.2014.04.004
.
Bowes
M. J.
,
Smith
J. T.
,
Jarvie
H. P.
&
Neal
C.
2008
Modelling of phosphorus inputs to rivers from diffuse and point sources
.
Science of The Total Environment
395
,
125
138
.
doi:10.1016/j.scitotenv.2008.01.054
.
Chen
Q.
,
Mei
K.
,
Dahlgren
R. A.
,
Wang
T.
,
Gong
J.
&
Zhang
M.
2016
Impacts of land use and population density on seasonal surface water quality using a modified geographically weighted regression
.
Science of the Total Environment
572
,
450
466
.
Council Directive 91/676/EEC of 12 December 1991
Concerning the protection of waters against pollution caused by nitrates from agricultural sources
.
Official Journal L
375
,
0001
0008
.
Council Directive 2000/60/EC of 22 December 2000 of the European Parliament and of the Council establishing a framework for Community action in the field of water policy
.
Official Journal L
327
,
1
73
.
Djodjic
F.
,
Börling
K.
&
Bergström
L.
2004
Phosphorus leaching in relation to soil type and soil phosphorus content
.
Journal of Environmental Quality
33
,
678
684
.
Edwards
A. C.
&
Withers
P. J. A.
2007
Linking phosphorus sources to impacts in different types of water body
.
Soil Use and Management
23
,
133
143
.
doi:10.1111/j.1475-2743.2007.00110.x
.
Estonian Weather Service
2021
www.ilmateenistus.ee (accessed 15 November 2021). Available from: https://www.ilmateenistus.ee/siseveed/ajaloolised-vaatlusandmed/vooluhulgad/ (accessed 26 September 2021)
.
European Environment Agency
2005
Source Apportionment of Nitrogen and Phosphorus Inputs Into the Aquatic Environment. EEA Report no. 7
.
European Environment Agency
,
Copenhagen
,
Denmark
, pp.
1
48
.
HELCOM
2018
HELCOM Thematic Assessment of Eutrophication 2011–2016. Baltic Sea Environment Proceedings No. 156
.
HELCOM
2021
Baltic Sea Action Plan—2021 Update. Adopted in Lübeck
.
Hong
B.
,
Swaney
D. P.
,
McCrackin
M.
,
Svanbäck
A.
,
Humborg
C.
,
Gustafsson
B.
,
Yershova
A.
&
Pakhomau
A.
2017
Advances in NANI and NAPI accounting for the Baltic drainage basin: spatial and temporal trends and relationships to watershed TN and TP fluxes
.
Biogeochemistry
133
,
245
261
.
doi:10.1007/s10533-017-0330-0
.
ICES
2019
Baltic Sea Ecoregion – Ecosystem Overview. Report ICES Advice: Ecosystem Overviews
.
Iital
A.
,
Pachel
K.
,
Loigu
E.
,
Pihlak
M.
&
Leisk
Ü
.
2010
Recent trends in nutrient concentrations in Estonian rivers as a response to large-scale changes in land-use intensity and life-styles
.
Journal of Environmental Monitoring
12
,
178
188
.
doi:10.1039/B912923E
.
Iital
A.
,
Klõga
M.
,
Pihlak
M.
,
Pachel
K.
,
Zahharov
A.
&
Loigu
E.
2014
Exploring the impact of agricultural land use on the nitrogen content and trends in streams in Estonia
.
Agriculture, Ecosystems & Environment
198
,
44
53
.
KESE
2021
Environmental Monitoring Information System
.
Available from: https://kese.envir.ee/ (accessed 26 May 2021)
.
Kiedrzyńska
E.
,
Kiedrzyński
M.
,
Urbaniak
M.
,
Magnuszewski
A.
,
Skłodowski
M.
,
Wyrwicka
A.
&
Zalewski
M.
2014
Point sources of nutrient pollution in the lowland river catchment in the context of the Baltic Sea eutrophication
.
Ecological Engineering
70
,
337
348
.
doi:10.1016/j.ecoleng.2014.06.010
.
Kodinariya
T.
&
Makwana
P.
2013
Review on determining of cluster in K-means clustering
.
International Journal of Advance Research in Computer Science and Management Studies
1
,
90
95
.
Kopáček
J.
,
Hejzlar
J.
,
Porcal
P.
&
Znachor
P.
2021
Biogeochemical causes of sixty-year trends and seasonal variations of river water properties in a large European basin
.
Biogeochemistry
154
,
81
98
.
doi:10.1007/s10533-021-00800-z
.
Lehtsaar
J.
&
Jullinen
V.
1998
Rural financial problems and the experience of their solution in Estonia
. In:
Ivanova, I. (ed.)
.
Latvijas tautsaimnieciba: stavoklis, problemas, risinajumi, Jelgava (Latvia), 8 Apr 1998
.
Latvian University of Agriculture
,
Jelgava
.
Löw
P.
,
Osterburg
B.
&
Klages
S.
2021
Comparison of regulatory approaches for determining application limits for nitrogen fertilizer use in Germany
.
Environmental Research Letters
16
,
055009
.
doi:10.1088/1748-9326/abf3de
.
Martínez-Dalmau
J.
,
Berbel
J.
&
Ordóñez-Fernández
R.
2021
Nitrogen fertilization. A review of the risks associated with the inefficiency of its use and policy responses
.
Sustainability
13
,
5625
.
doi:10.3390/su13105625
.
Meisinger
J. J.
&
Delgado
J. A.
2002
Principles for managing nitrogen leaching
.
Journal of Soil and Water Conservation
57
(
6
),
485
498
.
Paerl
H. W.
,
Scott
J. T.
,
McCarthy
M. J.
,
Newell
S. E.
,
Gardner
W. S.
,
Havens
K. E.
,
Hoffman
D. K.
,
Wilhelm
S. W.
,
Wurtsbaugh
W. A.
2016
It takes two to Tango: when and where dual nutrient (N & P) reductions are needed to protect lakes and downstream ecosystems
.
Environmental Science & Technology
50
,
10805
10813
.
doi:10.1021/acs.est.6b02575
.
Pengerud
A.
,
Stålnacke
P.
,
Bechmann
M.
,
Blicher-Mathiesen
G.
,
Iital
A.
,
Koskiaho
J.
,
Kyllmar
K.
,
Lagzdins
A.
&
Povilaitis
A.
2015
Temporal trends in phosphorus concentrations and losses from agricultural catchments in the Nordic and Baltic countries
.
Acta Agriculturae Scandinavica, Section B – Soil & Plant Science
65
,
173
185
.
doi:10.1080/09064710.2014.993690
.
Pinnaveekogumite nimekiri, pinnaveekogumite ja territoriaalmere seisundiklasside määramise kord, pinnaveekogumite ökoloogiliste seisundiklasside kvaliteedinäitajate väärtused ja pinnaveekogumiga hõlmamata veekogude kvaliteedinäitajate väärtused. Keskkonnaministri määrus nr. 19, 16.04.2020. Riigi Teataja I, 21.04.2020, 61, Available from: https://www.riigiteataja.ee/akt/121042020061.
Sabo
R. D.
,
Clark
C. M.
&
Compton
J. E.
2021
Considerations when using nutrient inventories to prioritize water quality improvement efforts across the US
.
Environmental Research Communications
3
,
045005
.
Stålnacke
P.
,
Aakerøy
P. A.
,
Blicher-Mathiesen
G.
,
Iital
A.
,
Jansons
V.
,
Koskiaho
J.
,
Kyllmar
K.
,
Lagzdins
A.
,
Pengerud
A.
&
Povilaitis
A.
2014
Temporal trends in nitrogen concentrations and losses from agricultural catchments in the Nordic and Baltic countries
.
Agriculture, Ecosystems & Environment
198
,
94
103
.
Statistics Estonia database
.
2021
Available from: https://andmed.stat.ee (accessed 7 May 2021)
.
Tarand
A.
,
Jaagus
J.
&
Kallis
A.
2013
Eesti klima minevikus ja Tänapäeval
.
Tartu Ülikooli Kirjastus, Tartu
, p.
631
.
TIBCO Software Inc
.
2017
Statistica (Data Analysis Software System), Version 13
.
Available from: http://statistica.io.
Timm
H.
,
Järvekülg
R.
,
Pall
P.
&
Vilbaste
S.
2019
Eesti jõed (Estonian Rivers)
.
Varrak
,
Tallinn
. p.
339
.
Tsakiris
G.
2015
The status of the European waters in 2015: a review
.
Environmental Processes
2
,
543
557
.
Viru
B.
&
Jaagus
J.
2020
Spatio-temporal variability and seasonal dynamics of snow cover regime in Estonia
.
Theoretical and Applied Climatology
139
,
759
771
.
Water Act
2019
Riigikogu
.
Available from: www.riigiteataja.ee/en/eli/527122019007/consolide (accessed 27 April 2022)
.
Xia
Y.
,
Zhang
M.
,
Tsang
D. C. W.
,
Geng
N.
,
Lu
D.
,
Zhu
L.
,
Igalavithana
A. D.
,
Dissanayake
P. D.
,
Rinklebe
J.
,
Yang
X.
&
Ok
Y. S.
2020
Recent advances in control technologies for non-point source pollution with nitrogen and phosphorous from agricultural runoff: current practices and future prospects
.
Applied Biological Chemistry
63
(
1
),
1
13
.
Zobena
A.
1998
Organic farming in the baltic countries social aspects of development
.
Biopolitics
6
. .
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