Space and time variations of watershed N and P budgets and their relationships with reactive N and P loadings in a heavily impacted river basin (Po river, Northern Italy)

Climate change impacts on eutrophication in the Po River (Italy): Temperature-mediated reduction in nitrogen export but no effect on phosphorus

Rivers worldwide are under stress from eutrophication and nitrate pollution, but the ecological consequences overlap with climate change, and the resulting interactions may be unexpected and still unexplored. The Po River basin (northern Italy) is one of the most agriculturally productive and densely populated areas in Europe. It remains unclear whether the climate change impacts on the thermal and hydrological regimes are already affecting nutrient dynamics and transport to coastal areas. The present work addresses the long-term trends (1992-2020) of nitrogen and phosphorus export by investigating both the annual magnitude and the seasonal patterns and their relationship with water temperature and discharge trajectories. Despite the constant diffuse and point sources in the basin, a marked decrease (-20%) in nitrogen export, mostly as nitrate, was recorded in the last decade compared to the 1990s, while no significant downward trend was observed for phosphorus. The water temperature of the Po River has warmed, with the most pronounced signals in summer (+0.13°C/year) and autumn (+0.16°C/year), together with the strongest increase in the number of warm days (+70%-80%). An extended seasonal window of warm temperatures and the persistence of low flow periods are likely to create favorable conditions for permanent nitrate removal via denitrification, resulting in a lower delivery of reactive nitrogen to the sea. The present results show that climate change-driven warming may enhance nitrogen processing by increasing respiratory river metabolism, thereby reducing export from spring to early autumn, when the risk of eutrophication in coastal zones is higher.

Variability of hydrographic and biogeochemical properties in the North-western Adriatic coastal waters in relation to river discharge and climate changes

The dynamics of hydrographic and biogeochemical properties in a Northwestern coastal area of the Adriatic Sea were investigated. The time series data from continuous observation (2007-2022) allowed the investigation of annual trends and seasonal cycles along a coastal transect influenced by local river discharge. Various statistical models were used to investigate water temperature, salinity, chlorophyll a, dissolved organic, inorganic and particulate nutrients, precipitation and river discharge. It was found that the local river discharge regime played an essential role in interannual, and seasonal biogeochemical dynamics associated with global climate change in the Mediterranean region. A significant trend towards oligotrophic conditions was detected, as evidenced by the downward trend in the river mouth and on the sea of chlorophyll a (-0.2 μg L-1 in the sea), dissolved organic and inorganic nitrogen and phosphorus (i.e., -0.43 μM yr-1 of DON in the sea and -6.67 of DIN μM yr-1 in the river mouth or -0.07 μM yr-1 of DOP and -0.02 μM yr-1 of DIP in the river mouth) and silicate (-2.47 μM yr-1 in the river mouth) concentrations. Salinity showed a long-term increase in the sea (0.08 yr-1), corresponding to a significant decrease in water discharge from the local river (-0.27 m3 s-1 yr-1) and precipitation (-0.06 mm yr-1). The dissolved organic and inorganic nutrients highlighted a different seasonal accumulation under the river runoff regime. The nutrient enrichment was predominantly driven by river contribution. Data analysis showed that the coastal biogeochemical properties dynamics were mostly influenced by river discharge and precipitation regimes, which in turn are driven by climate change variability in the North-western Adriatic Sea.

Assessment of the eutrophication status at Mediterranean sub-basin scale, within the European Marine Strategy Framework Directive

The aim of this work is to define harmonized reference conditions and assessment thresholds for selected criteria elements of the Marine Strategy Framework Directive (MSFD) Descriptor 5 (Eutrophication) in the Eastern Mediterranean Sea and to test if a tool for integrated assessment of the status of marine systems can be used as a common methodological approach. In this frame, we tested two statistical approaches in order to set threshold values for four criteria of Descriptor 5: nutrients, chlorophyll a, transparency and dissolved oxygen in the bottom waters. It is noteworthy that this work revealed the need to apply common procedures in data treatment and assessment evaluation. This is the first attempt to set common methods for the assessment of eutrophication in the Eastern Mediterranean, which is essential in marine environments, especially those shared by several countries. To this end, we have applied common criteria and metrics and established thresholds "Good" and "Moderate" for nutrients, chlorophyll a, transparency and dissolved oxygen in the bottom waters for the different Water Types of the Adriatic and Aegean Seas (I, II, IIIW, IIIE), based on datasets provided by Italy, Slovenia, Croatia and Greece. The selected criteria elements were common for all countries, providing a unified approach to GES assessment of two case study areas: the Adriatic Sea and the Saronikos Gulf. Dissolved Inorganic Nitrogen (DIN) threshold values of 15.6, 6.85, 1.61 and 2.11 μmol L-1 were set for the Water Types I, II, IIIW and IIIE, respectively. We also tested if an aggregation tool for GES assessment, such as Nested Environmental status Assessment Tool (NEAT), could be used as a common methodological approach. The comparison of NEAT with TRIX showed good comparability. In this end, NEAT can be used as a useful and much needed assessment tool for assessing eutrophication status of the marine.

Temporal and spatial differences in nitrogen and phosphorus biogeochemistry and ecosystem functioning of a hypertrophic lagoon (Curonian Lagoon, SE Baltic Sea) revealed via Ecological Network Analysis

In coastal lagoons, eutrophication and hydrology are interacting factors that produce distortions in biogeochemical nitrogen (N) and phosphorus (P) cycles. Such distortions affect nutrient relative availability and produce cascade consequences on primary producer's community and ecosystem functioning. In this study, the seasonal functioning of a coastal lagoon was investigated with a multielement approach, via the construction and analysis of network models. Spring and summer networks, both for N and P flows, have been simultaneously compiled for the northern transitional and southern confined area of the hypertrophic Curonian Lagoon (SE Baltic Sea). Ecological Network Analysis was applied to address the combined effect of hydrology and seasonality on biogeochemical processes. Results suggest that the ecosystem is more active and presents higher N and P fluxes in summer compared to spring, regardless of the area. Furthermore, larger internal recycling characterizes the confined compared to the transitional area, regardless of the season. The two areas differed in the fate of available nutrients. The transitional area received large riverine inputs that were mainly transferred to the sea without the conversion into primary producers' biomass. The confined area had fewer inputs but proportionally larger conversion into phytoplankton biomass. In summer, particularly in the confined area, primary production was inefficiently consumed by herbivores. Most phytoplanktonic N and P, in the confined area more than in the transitional area, were conveyed to the detritus pathway where P, more than N, was recycled, contributing to the unbalance in N:P stoichiometry and favouring N-fixing cyanobacteria over other phytoplankton groups. The findings of this study provide a comprehensive understanding of N and P circulation patterns in lagoon areas characterized by different hydrology. They also support the importance of a stoichiometric approach to trace relative differences in N and P recycling and abundance, that promote blooms, drive algal communities and whole ecosystem functioning.

A diatom-based approach to refine nutrient concentrations compatible with the "good" status of Northern Italy rivers

The identification of ecologically sound thresholds represents an important step toward improving the ecological status of rivers through appropriate measures to contain nutrient loads. The aim of the present study was to estimate phosphorus and nitrogen concentrations compatible with the achievement of the "good" ecological status of rivers from data collected in the Po River District, the largest hydrographic system in Italy. For this purpose, relationships between the diatom index used in Italy for the national assessment of the stream ecological status, the ICMi (Intercalibration Common Metric index), and total phosphorus and nitrate concentrations were analyzed using monitoring data collected between 2009 and 2019. The Po River Basin encompasses five distinct river types, from Alpine to Mediterranean to Lowlands, characterized by different anthropogenic pressures and water quality. Through regression analysis between the ICMi and nutrient concentrations, we estimated ranges of the latter values corresponding to a "good" ecological status for each river type. The resulting thresholds are far more stringent than the limits set by the Italian legislation for water quality classification. This is particularly true for total phosphorus, whose threshold value should be roughly halved for all river types. For nitrates, the results are more differentiated according to river type: the estimated thresholds are much more stringent than those currently in use for siliceous Alpine and Mediterranean rivers. Moreover, the availability of such a large database allowed also to assess the influence of one nutrient over the other on the diatom community and to highlight some critical issues in the formulation of ICMi for Mediterranean rivers.

Groundwater is a hidden global keystone ecosystem

Groundwater is a vital ecosystem of the global water cycle, hosting unique biodiversity and providing essential services to societies. Despite being the largest unfrozen freshwater resource, in a period of depletion by extraction and pollution, groundwater environments have been repeatedly overlooked in global biodiversity conservation agendas. Disregarding the importance of groundwater as an ecosystem ignores its critical role in preserving surface biomes. To foster timely global conservation of groundwater, we propose elevating the concept of keystone species into the realm of ecosystems, claiming groundwater as a keystone ecosystem that influences the integrity of many dependent ecosystems. Our global analysis shows that over half of land surface areas (52.6%) has a medium-to-high interaction with groundwater, reaching up to 74.9% when deserts and high mountains are excluded. We postulate that the intrinsic transboundary features of groundwater are critical for shifting perspectives towards more holistic approaches in aquatic ecology and beyond. Furthermore, we propose eight key themes to develop a science-policy integrated groundwater conservation agenda. Given ecosystems above and below the ground intersect at many levels, considering groundwater as an essential component of planetary health is pivotal to reduce biodiversity loss and buffer against climate change.

Looking back to move forward: Restoring vegetated canals to meet missing Water Framework Directive goals in agricultural basins

Nitrate pollution and eutrophication remain pressing issues in Europe regarding the quality of aquatic ecosystems and the safety of drinking water. Achieving water quality goals under the Water Framework Directive (WFD) has proven to be particularly challenging in agricultural catchments, where high nitrate concentrations are the main reason for the failure of many water bodies to meet a good ecological status. Canals and ditches are common man-made features of irrigated and drained landscapes and, when vegetated, have recently been identified as denitrification hotspots. By combining experimental data and GIS-based upscaling estimation, the potential capacity of the canal network to reduce nitrate loads was quantified in several scenarios differing in the level of nitrate pollution and in the extent of the canal network length where conservative management practices are implemented. The analysis was carried out in the irrigated lowlands of the Po River basin, which is the largest hydrographic system in Italy and a global hotspot for nitrogen inputs and eutrophication. Scenario simulations showed that maintaining aquatic vegetation in at least 25 % of the canal network length, selecting sites with high nitrate availability (>2.4 mg N L-1), would promote a greater potential for permanent N removal. The increased denitrification capacity would meet the load reduction target required to achieve a WFD good ecological status in waters draining into the Adriatic Sea during the spring-summer months, when the eutrophication risk is higher. Promoting denitrification in the canal network by postponing the mowing of in-stream vegetation to the end of the growing season could be an effective mitigation strategy to improve water quality in agricultural basins and contribute to achieving the WFD goals.

Virus-pathogen interactions improve water quality along the Middle Route of the South-to-North Water Diversion Canal

Bacterial pathogens and viruses are the leading causes of global waterborne diseases. Here, we discovered an interesting natural paradigm of water "self-purification" through virus-pathogen interactions over a 1432 km continuum along the Middle Route of the South-to-North Water Diversion Canal (MR-SNWDC) in China, the largest water transfer project in the world. Due to the extremely low total phosphorus (TP) content (ND-0.02 mg/L) in the MR-SNWDC, the whole canal has experienced long-lasting phosphorus (P) limitation since its operation in 2015. Based on 4443 metagenome-assembled genomes (MAGs) and 40,261 nonredundant viral operational taxonomic units (vOTUs) derived from our recent monitoring campaign, we found that residential viruses experiencing extreme P constraints had to adopt special adaptive strategies by harboring smaller genomes to minimize nucleotide replication, DNA repair, and posttranslational modification costs. With the decreasing P supply downstream, bacterial pathogens showed repressed environmental fitness and growth potential, and a weakened capacity to maintain P acquisition, membrane formation, and ribonucleotide biosynthesis. Consequently, the unique viral predation effects under P limitation, characterized by enhanced viral lytic infections and an increased abundance of ribonucleotide reductase (RNR) genes linked to viral nuclear DNA replication cycles, led to unexpectedly lower health risks from waterborne bacterial pathogens in the downstream water-receiving areas. These findings highlighted the great potential of water self-purification associated with virus-pathogen dynamics for water-quality improvement and sustainable water resource management.

Decoding river pollution trends and their landscape determinants in an ecologically fragile karst basin using a machine learning model

Karst watersheds accommodate high landscape complexity and are influenced by both human-induced and natural activity, which affects the formation and process of runoff, sediment connectivity and contaminant transport and alters natural hydrological and nutrient cycling. However, physical monitoring stations are costly and labor-intensive, which has confined the assessment of water quality impairments on spatial scale. The geographical characteristics of catchments are potential influencing factors of water quality, often overlooked in previous studies of highly heterogeneous karst landscape. To solve this problem, we developed a machining learning method and applied Extreme Gradient Boosting (XGBoost) to predict the spatial distribution of water quality in the world's most ecologically fragile karst watershed. We used the Shapley Addition interpretation (SHAP) to explain the potential determinants. Before this process, we first used the water quality damage index (WQI-DET) to evaluate the water quality impairment status and determined that COD_Mn, TN and TP were causing river water quality impairments in the WRB. Second, we selected 46 watershed features based on the three key processes (sources-mobilization-transport) which affect the temporal and spatial variation of river pollutants to predict water quality in unmonitored reaches and decipher the potential determinants of river impairments. The predicting range of COD_Mn spanned from 1.39 mg/L to 17.40 mg/L. The predictions of TP and TN ranged from 0.02 to 1.31 mg/L and 0.25-5.72 mg/L, respectively. In general, the XGBoost model performs well in predicting the concentration of water quality in the WRB. SHAP explained that pollutant levels may be driven by three factors: anthropogenic sources (agricultural pollution inputs), fragile soils (low organic carbon content and high soil permeability to water flow), and pollutant transport mechanisms (TWI, carbonate rocks). Our study provides key data to support decision-making for water quality restoration projects in the WRB and information to help bridge the science:policy gap.

Spatial characteristics of nutrient budget on town scale in the Three Gorges Reservoir area, China

Accurately quantifying nutrient budget is an essential step toward sustainable nutrient management in large watersheds increasingly disturbed by human activity. A town-scale nutrient budget framework based on the Soil and Water Assessment Tool was developed for 2010-2012 in the Three Gorges Reservoir area in China (TGRA). Moran's I spatial correlation test and Geodetector spatial heterogeneity test were employed to systematically analyze the spatial characteristics of the resulting nutrient budget. The Moran's I value of total nitrogen (TN) and total phosphorus (TP) gradually increased from input to output in the range of 0.091-0.232 and 0.102-0.484, respectively. Towns with higher TN and TP inputs were largely concentrated in the main urban area of Chongqing because of its high population density. By contrast, towns with higher TN and TP outputs were concentrated in the head of the TGRA. The Moran's I values of the TN and TP retention coefficients (R) were 0.433 and 0.524, respectively, demonstrating clear spatial consistency. Towns with a "High-high" spatial consistency pattern and positive R value were concentrated in the tail and hinterland, while those with a "Low-low" spatial consistency pattern and negative coefficient value were located mainly in the head of the TGRA. This phenomenon was mostly caused by differences in regional elevation, the normalized difference vegetation index, and soil erosion factor. The interaction effect between any two of these three factors on nutrient retention (Geodetector q-value) was greater than 60%. Therefore, future nutrient management should be based on a full understanding of regional biophysical conditions, especially in large areas. These findings provide a new perspective on fine nutrient management.

Short-Term Effects of the EU Nitrate Directive Reintroduction: Reduced N Loads to River from an Alluvial Aquifer in Northern Italy

The Po Plain (northern Italy) is one of the largest aquifers in Europe, and 67% of the utilized agricultural land in this area is classified as a nitrate vulnerable zone (NVZ). However, it hosts intensive agriculture and livestock farming. In a stretch of the Mincio River (a tributary of the Po River), hydraulic heads and physico-chemical parameters of river and groundwater were monitored for a hydrologic year (2020–2021), to evaluate the effects of manure fertilization and flooding irrigation on surface- and groundwater chemistry. From 2020 the Nitrate Directive’s fertilization limit was reintroduced and a comparison has been performed comparing surface- and groundwater data from the 2019 fertilization period (before limit reintroduction) and 2020 (after). Results suggest that in 2021 the phreatic aquifer displayed elevated nitrate (NO3−) concentrations, exceeding 50 mg L−1, although average values were lower than those of 2019. Nitrate loads in the Mincio River reached 6670 kg NO3− d−1 and resulted from the overfertilization in the surrounding area and the quick transfer of nitrogen from groundwater to the river. As compared to 2019, the river loads decreased by 59%, suggesting that the introduction of fertilization limits can produce measurable, short-term responses in alluvial aquifers.

Assessment of Spatio-Temporal Variability of Faecal Pollution along Coastal Waters during and after Rainfall Events

More than 80% of wastewaters are discharged into rivers or seas, with a negative impact on water quality along the coast due to the presence of potential pathogens of faecal origin. Escherichia coli and enterococci are important indicators to assess, monitor, and predict microbial water quality in natural ecosystems. During rainfall events, the amount of wastewater delivered to rivers and coastal systems is increased dramatically. This study implements measures capable of monitoring the pathways of wastewater discharge to rivers and the transport of faecal bacteria to the coastal area during and following extreme rainfall events. Spatio-temporal variability of faecal microorganisms and their relationship with environmental variables and sewage outflow in an area located in the western Adriatic coast (Fano, Italy) was monitored. The daily monitoring during the rainy events was carried out for two summer seasons, for a total of five sampling periods. These results highlight that faecal microbial contaminations were related to rainy events with a high flow of wastewater, with recovery times for the microbiological indicators varying between 24 and 72 h and influenced by a dynamic dispersion. The positive correlation between ammonium and faecal bacteria at the Arzilla River and the consequences in seawater can provide a theoretical basis for controlling ammonium levels in rivers as a proxy to monitor the potential risk of bathing waters pathogen pollution.

The achievement of Water Framework Directive goals through the restoration of vegetation in agricultural canals

Decreasing nitrate concentrations is one of the most relevant Water Framework Directive (WFD) goals, which today is still unreached in several European countries. Vegetated canals have been recognized as effective filters to mitigate nitrate pollution, although rarely included in restoration programs aimed at improving water quality in agricultural watersheds. The Po di Volano basin (713 km², Northern Italy) is a deltaic territory crossed by an extensive network of agricultural canals (~1300 km). The effectiveness in buffering nitrate loads via denitrification was assessed for different levels of in-stream emergent vegetation maintenance by employing an upscale model based on extensive datasets of field measurements. The scenarios differed for the canal network length (5%, 20%, 40%, and 60%) where conservative management practices were adopted by postponing the mowing operations from the middle of summer, as nowadays, to the early autumn, i.e., the vegetative season end. The scenario simulations demonstrated that the capacity to mitigate diffuse nitrate pollution would increase up to four times, compared to the current condition (5% scenario), by postponing the vegetation mowing to the end of the vegetative season in 60% of the canal network length. By preserving the in-stream vegetation in 20% of the canal network, its denitrification capacity would equal the nitrate load reduction target required for achieving, from May to September, the good ecological status according to the WFD in waters delivered to the coastal areas. Changing the timing of vegetation mowing may create a large potential for permanent nitrate removal via denitrification in agricultural landscapes, thus protecting the coastal areas when the eutrophication risk is higher. Conservative management practices of in-stream vegetation might be promoted as an effective low-cost tool to be included in the WFD implementation strategies.

Governance and groundwater modelling: Hints to boost the implementation of the EU Nitrate Directive. The Lombardy Plain case, N Italy

The EU Nitrate Directive has been ruling for almost 30 years, nevertheless nitrate concentration in the Lombardy Plain did not decrease. Together with failures of management implementation, a possible cause for such field observations is that management actions were taken without adequately considering the actual hydrogeological dynamics. To consider this aspect, the paper presents a groundwater flow and transport numerical model of a specific area of the Lombardy Plain. The aim of this model is to demonstrate how modelling, as a management tool, can be useful in the governance process. The groundwater model, using well-known MODFLOW-MT3D codes, is based on existing hydrogeological information, while a nitrogen mass balance has been performed at municipal scale to determine the agricultural N surplus to the subsurface. The model adequately reproduces head levels and nitrate concentrations in observation wells for a 10-year simulation period, showing that 4.5% of the N annual input remains stored in the system. The model indicates the efficiency of rivers and springs to export N out from the system at an estimated rate of 77.5% of the annual N inputs. Back to governance, the model shows that management data at municipal level (e.g. irrigation rates, groundwater withdrawal, N net recharge) provide a satisfactory scale for successfully reproducing nitrate evolution. Hence those variables that can be object of debate during a governance process can be treated as input data to the numerical model. Therefore, backcasting exercises can be conducted to check whether the model outcome fits with the expected results of specific management actions. The model highlights how the N mass balance evolves, providing clues on which factors can be managed to reduce nitrate concentrations and meet the Directive's requirements. Numerical groundwater models, as an option to address water management issues, ultimately contribute to solve the information and capacity governance gaps.

Soil Denitrification, the Missing Piece in the Puzzle of Nitrogen Budget in Lowland Agricultural Basins

Denitrification is a key process buffering the environmental impacts of agricultural nitrate loads but, at present, remains the least understood and poorly quantified sink in nitrogen budgets at the watershed scale. The present work deals with a comprehensive and detailed analysis of nitrogen sources and sinks in the Burana–Volano–Navigabile basin, the southernmost portion of the Po River valley (Northern Italy), an intensively cultivated (> 85% of basin surface) low-lying landscape. Agricultural census data, extensive monitoring of surface–groundwater interactions, and laboratory experiments targeting N fluxes and pools were combined to provide reliable estimates of soil denitrification at the basin scale. In the agricultural soils of the basin, nitrogen inputs exceeded outputs by nearly 40% (~ 80 kg N ha−1 year−1), but this condition of potential N excess did not translate into widespread nitrate pollution. The general scarcity of inorganic nitrogen species in groundwater and soils indicated limited leakage and storage. Multiple pieces of evidence supported that soil denitrification was the process that needed to be introduced in the budget to explain the fate of the missing nitrogen. Denitrification was likely boosted in the soils of the studied basin, prone to waterlogged conditions and consequently oxygen-limited, owing to peculiar features such as fine texture, low hydraulic conductivity, and shallow water table. The present study highlighted the substantial contribution of soil denitrification to balancing nitrogen inputs and outputs in agricultural lowland basins, a paramount ecosystem function preventing eutrophication phenomena.

Human disturbance on phosphorus sources, processes and riverine export in a subtropical watershed

Phosphorus (P) is a key nutrient in freshwater systems, often acting as the limiting nutrient. The dominant sources of P in the Jiulong River watershed (S.E. China) are anthropogenic. Dissolved and particulate P species were measured in the West (WJR) and North (NJR) rivers during the wet and dry seasons of 2018 and at their river outlets during a storm (June 2019). Sources of P pollution were characterized from mainly single source subcatchments (dry season). The Agriculture source (WJR) had a total P of 114.7 ± 13.1 μg P L^-1, which was mainly dissolved inorganic P (DIP) from excess fertilizer washed from the fields. By contrast, the West Urban source (sewage effluent) was mainly particulate (POP) and dissolved organic P (DOP). The effect of reservoirs in the main NJR was to decrease total particulate P (TPP) and DIP and increase POP, due to increased sedimentation of particles and biological uptake. An increase in all P species was observed at the beginning of the storm, followed by a decrease on the rising hydrograph due to dilution. The final concentration of all P species was higher than baseflow, confirming that storms increase the P flux out of the watershed. P was initially washed off the fields during the storm, and during the falling hydrograph P increased due to interflow and other longer-term sources. The high DIN:DIP ratio confirmed the key importance of P inputs from human activities in substantially altering P sources and cycling, and hence the importance of science-based management to alleviate the eutrophication problem.

Influence of nitrogen-based fertilization on nitrates occurrence in groundwater of hilly vineyards

Nitrogen losses from intensive agricultural production may end up as high nitrate (NO₃^-) concentrations in groundwater, with a long-term impact on groundwater quality. The main objective of this study was to evaluate the impact of fertilization practices used for grape cultivation on groundwater quality of Tidone Valley, northwest of Italy, following an integrated socio-hydrogeological approach that consists on (i) the involvement of 175 farmers in the description of agricultural and fertilization practices, using a survey of ad hoc questionnaires, (ii) the evaluation of NO₃^- occurrence in groundwater and (iii) the identification of NO₃^- sources through isotopic and hydrochemical analysis. In this area, as for certain particular Apennines shallow aquifers, groundwater is of reduced interest due to its limited storage capacity and there are insufficient wells currently monitored by the local Environmental Agency (ARPAE) to evaluate the impact of agricultural fertilization on existing local aquifers. Farmers' questionnaires results highlighted an extensive use of inorganic nitrogen fertilization and a tendency of farmers to follow their own experience for fertilization. Chemical analyses revealed high variability of major and trace elements concentrations isotope data. NO₃^- concentrations were significantly higher in deeper wells with respects to shallow wells. Isotopic results indicated that groundwater NO₃^- origin is inorganic, in agreement with the land use and the declared viticultural practices. Comparing groundwater NO₃^- occurrence from the studied area with values of entire Emilia-Romagna Region, only 7.7% of groundwater samples showed values above the EQS. (50 mg NO₃^-/L) between Nov 2017 and Sept 2018, while in the entire region 11.5% of groundwater samples showed values above the EQS in the same period. Considering that the vineyards surface in the studied area represents almost 75% of the entire regional vineyard surface, the obtained results suggest a low to moderate impact of viticulture on NO₃^-concentration of regional groundwater.

Denitrification, Nitrogen Uptake, and Organic Matter Quality Undergo Different Seasonality in Sandy and Muddy Sediments of a Turbid Estuary

The interaction between microbial communities and benthic algae as nitrogen (N) regulators in poorly illuminated sediments is scarcely investigated in the literature. The role of sediments as sources or sinks of N was analyzed in spring and summer in sandy and muddy sediments in a turbid freshwater estuary, the Curonian Lagoon, Lithuania. Seasonality in this ecosystem is strongly marked by phytoplankton community succession with diatoms dominating in spring and cyanobacteria dominating in summer. Fluxes of dissolved gas and inorganic N and rates of denitrification of water column nitrate (D_w) and of nitrate produced by nitrification (D_n) and sedimentary features, including the macromolecular quality of organic matter (OM), were measured. Shallow/sandy sites had benthic diatoms, while at deep/muddy sites, settled pelagic microalgae were found. The OM in surface sediments was always higher at muddy than at sandy sites, and biochemical analyses revealed that at muddy sites the OM nutritional value changed seasonally. In spring, sandy sediments were net autotrophic and retained N, while muddy sediments were net heterotrophic and displayed higher rates of denitrification, mostly sustained by D_w. In summer, benthic oxygen demand increased dramatically, whereas denitrification, mostly sustained by D_n, decreased in muddy and remained unchanged in sandy sediments. The ratio between denitrification and oxygen demand was significantly lower in sandy compared with muddy sediments and in summer compared with spring. Muddy sediments displayed seasonally distinct biochemical composition with a larger fraction of lipids coinciding with cyanobacteria blooms and a seasonal switch from inorganic N sink to source. Sandy sediments had similar composition in both seasons and retained inorganic N also in summer. Nitrogen uptake by microphytobenthos at sandy sites always exceeded the amount loss via denitrification, and benthic diatoms appeared to inhibit denitrification, even in the dark and under conditions of elevated N availability. In spring, denitrification attenuated N delivery from the estuary to the coastal area by nearly 35%. In summer, denitrification was comparable (~100%) with the much lower N export from the watershed, but N loss was probably offset by large rates of N-fixation.

Groundwater-dependent ecosystems as transfer vectors of nitrogen from the aquifer to surface waters in agricultural basins: The fontanili of the Po Plain (Italy)

The most spread groundwater-dependent ecosystems in the River Po valley are semi-natural lowland springs called "fontanili". They provide specific habitats and support high biodiversity, but are often strongly impaired by agricultural pollution. In the present study we seasonally monitored the discharge and nitrogen concentration of 48 fontanili from the Adda and the Ticino river basins. We observed a wide spatial variability of both NO₃-N concentrations and flows. The annual NO₃-N loads ranged from <1 to 75 t y^-1 and < 1 to 29 t y^-1 in the Adda and Ticino basins respectively. In the springs characterized by variable discharge the N loads were exported mostly during the summer season when water table level was elevated mainly due to irrigation. Upscaling the mean NO₃-N load to each river catchment based on the total number of springs, we obtained an aerial export of 33.2 ± 6.0 and 12.5 ± 3.2 kg y^-1 ha^-1. Such loads accounted for the 30.4 and 21.5% of the N surplus estimated for the Adda and Ticino basins respectively. Random Forest analysis was performed to identify the most important environmental variables influencing the nitrate contamination in the spring waters. A total of 22 explanatory variables related to N sources, land uses, intrinsic hydrogeologic and soil proprieties, in "situ" and remotely sensed variables were considered. The percent of soil cultivated with maize in a 500 m radius buffer area surrounding the sampling site, the N from manure and the distance of each spring from the main river were the most effective factors in controlling the NO₃-N concentration in the fontanili water. The outcomes of this work open up to achievable management prospects for the protection and recovery of fontanili waters, and can be particularly useful for water managers in identifying areas and sites where restoration plans should be a priority.

The effects of hydrological extremes on denitrification, dissimilatory nitrate reduction to ammonium (DNRA) and mineralization in a coastal lagoon

Hydrological extremes of unusually high or low river discharge may deeply affect the biogeochemistry of coastal lagoons, but the effects are poorly explored. In this study, microbial nitrogen processes were analyzed through intact core incubations and ¹⁵N-isotope addition at three sites in the eutrophic Sacca di Goro lagoon (Northern Adriatic Sea) both under high discharge (spring) and after prolonged low discharge (late-summer) of the main freshwater inputs. Under high discharge/nitrate load, denitrification was the leading process and there was no internal recycling. The site located at the mouth of the main freshwater input and characterized by low salinity exhibited the highest denitrification rate (up to 1150 ± 81 μmol N m^-2 h^-1), mostly sustained by nitrification stimulated by burrowing macrofauna. In contrast, we recorded high internal recycling under low discharge, when denitrification dropped at all sites due to low nitrate concentrations, reduced bioturbation and nitrification. The highest recycling was measured at the sites close to the sea entrance and characterized by high salinity and particularly at the clams cultivated area (up to 1003 ± 70 μmol N m^-2 h^-1). At this site, internal recycling was sustained by ammonification of biodeposits, bivalve excretion and dissimilatory nitrate reduction to ammonium (DNRA), which represented 30% of nitrate reduction. Flash floods and high nitrate loads may overwhelm the denitrification capacity of the lagoon due to the reduced residence time and to the saturation of microbial enzymatic activity, resulting in high transport of nitrate to the sea. Prolonged dry periods favor large internal recycling, due to a combination of high temperatures, low oxygen solubility and low bioturbation, which may prolong the extent of algal blooms with negative effects on lagoon biogeochemical services. We conclude that hydrological extremes, which are expected to become more frequent under climate change scenarios, strongly alter N cycling in coastal sediments.

Climatic and Anthropogenic Impacts on Environmental Conditions and Phytoplankton Community in the Gulf of Trieste (Northern Adriatic Sea)

During the last century, human activities have exerted an increasing pressure on coastal ecosystems, primarily inducing their eutrophication, with a more recent partial mitigation of this phenomenon where improvements of environmental management practices were adopted. However, a reanalysis of the pressures on coastal zones and surrounding drainage basins is needed because of the alterations induced nowadays by the climate changes. A comparative analysis of long-term oceanographic and environmental data series (1986–2018) was performed, in order to highlight the effects of anthropogenic and climatic disturbances on the phytoplankton community in the Gulf of Trieste (GoT). After the 1980s, the decline in phytoplankton abundance was matched to increasing periods of low runoff, an overall deficit of the precipitation and to a decrease in phosphate availability in the coastal waters (−0.003 µmol L−1 yr−1), even in the presence of large riverine inputs of nitrogen and silicates. This trend of oligotrophication was reversed in the 2010s by the beginning of a new and unexpected phase of climatic instability, which also caused changes of the composition and seasonal cycle of the phytoplankton community. Beyond the management of nutrient loads, it was shown that climatic drivers such as seawater warming, precipitation and wind regime affect both nutrient balance and phytoplankton community in this coastal zone.

Reactive Silica Traces Manure Spreading in Alluvial Aquifers Affected by Nitrate Contamination: A Case Study in a High Plain of Northern Italy

In the northern sector of the Po River Plain (Italy), widespread intensive agriculture and animal farming are supported by large amounts of water from Alpine lakes and their emissaries. Flood irrigation and excess fertilization with manure affect both the hydrology and the chemical quality of surface and groundwater, resulting in diffuse nitrogen pollution. However, studies analyzing the mechanisms linking agricultural practices with vertical and horizontal nitrogen paths are scarce in this area. We investigated groundwater quality and quantity in an unconfined, coarse-grained alluvial aquifer adjacent to the Mincio River (a tributary of the Po River), where steep summer gradients of nitrate (NO3−) concentrations are reported. The effects of manure on solutes’ vertical transport during precipitation events in fertilized and in control soils were simulated under laboratory conditions. The results show high SiO2 and NO3− leaching in fertilized soils. Similarly, field data are characterized by high SiO2 and NO3− concentrations, with a comparable spatial distribution but a different temporal evolution, suggesting their common origin but different processes affecting their concentrations in the study area. Our results show that SiO2 can be used as a conservative tracer of manure spreading, as it does not undergo biogeochemical processes that significantly alter its concentrations. On the contrary, nitrate displays large short-term variations related to aquifer recharge (i.e., flood irrigation and precipitation). In fact, aquifer recharge may promote immediate solubilization and stimulate nitrification, resulting in high NO3− concentrations up to 95.9 mg/L, exceeding the Water Framework Directive (WFD) thresholds. When recharge ends, anoxic conditions likely establish in the saturated zone, favoring denitrification and resulting in a steep decrease in NO3− concentrations.

Seasonal and Interannual Trends of Oceanographic Parameters over 40 Years in the Northern Adriatic Sea in Relation to Nutrient Loadings Using the EMODnet Chemistry Data Portal

Long-term data series (1971–2015) of physical and biogeochemical parameters were analyzed in order to assess trends and variability of oceanographic conditions in the northern Adriatic Sea (NAS), a mid-latitude shallow continental shelf strongly impacted by river discharges, human activities and climate changes. Interpolation maps and statistical models were applied to investigate seasonal and spatial variability, as well as decadal trends of temperature, salinity, chlorophyll-a and nutrients. This analysis shows that sea surface temperature increased by +0.36% year−1 over four decades. Annual mean flow of the Po River markedly changed due to the occurrence of periods of persistent drought, whereas the frequency of flow rates higher than 3000 m3 s−1 decreased between 2006 and 2015. Moreover, we observed a long-term decrease in surface phosphate concentrations in Po River water (−1.34% year−1) and in seawater (in summer −2.56% year−1) coupled, however, to a significant increase in nitrate concentration in seawater (+3.80% year−1) in almost all seasons. These changes indicate that the nutrient concentrations in the NAS have been largely modulated, in the last forty years, by the evolution of environmental management practices and of the runoff. This implies that further alteration of the marine environment must be expected as a consequence of the climate changes.

Introducing Life Cycle Assessment in Costs and Benefits Analysis of Vegetation Management in Drainage Canals of Lowland Agricultural Landscapes

Nitrate pollution remains an unsolved issue worldwide, causing serious effects on water quality and eutrophication of freshwater and brackish water environments. Its economic costs are still underestimated. To reduce nitrogen excess, constructed wetlands are usually recognized as a solution but, in recent years, interest has been raised in the role of ditches and canals in nitrogen removal. In this study, we investigated the environmental and economical sustainability of nitrogen removal capacity, using as a model study a lowland agricultural sub-basin of the Po River (Northern Italy), where the role of aquatic vegetation and related microbial processes on the mitigation of nitrate pollution has been extensively studied. Based on the Life Cycle Assessment (LCA) approach and costs and benefits analysis (CBA), the effectiveness of two different scenarios of vegetation management, which differ for the timing of mowing, have been compared concerning the nitrogen removal via denitrification and other terms of environmental sustainability. The results highlighted that postponing the mowing to the end of the vegetative season would contribute to buffering up to 90% of the nitrogen load conveyed by the canal network during the irrigation period and would reduce by an order of magnitude the costs of eutrophication potential.

Identification of the Nitrogen Sources in the Eocene Aquifer Area (Palestine)

Groundwater is the main source of water in many countries all over the world. Prevention of the pollution of this source is essential for a sustainable utilization. Nitrate pollution of groundwater is a common problem due to the association between intensive agriculture to achieve food security and fertilization. For an efficient management of groundwater pollution from nitrate, the first step would be to quantify the different sources of nitrogen in the aquifer of concern. This paper aims at demonstrating a general approach based on Geographic Information Systems (GIS) to characterize the spatial distribution of the nitrogen amounts in the area of the Eocene aquifer (Palestine). The aquifer is heavily utilized for agricultural and domestic water supply. Fertilization in the study area is a widespread practice. As a result, the aquifer is undergoing a nitrate pollution problem. The methodology relies mainly on specifying all the sources of nitrogen in the aquifer area using GIS to account for spatiality. Thereafter, GIS attribute tables and Excel spreadsheets were utilized to quantify the magnitudes of nitrogen from the different sources. Maps of the corresponding on-ground nitrate, ammonium, organic nitrogen and total nitrogen were developed for the study area. The results indicate that the total on-ground annual nitrogen loading in the study area is about 3260 tons of which 38% is attributed to fertilizers (chemical and manure) where the dominant form of nitrogen is NH4 (58.3%). The average total on-ground nitrogen loading is 7028 kg-N/km2·year. The estimated annual nitrate leaching to the aquifer is 1968 kg-N/km2. The areas of high sources of nitrogen have long-term impacts on the degradation of the water quality of the aquifer. It is therefore essential to build up on the outcomes of this work and to develop a nitrate fate and transport model for the Eocene aquifer. This model will enable the stakeholders to arrive at the efficient alternatives to manage the nitrate contamination of the aquifer.

In Search for the Missing Nitrogen: Closing the Budget to Assess the Role of Denitrification in Agricultural Watersheds

Although representing a paramount mechanism against nitrogen excess in agricultural landscapes, soil denitrification is still a largely unknown term in nitrogen balances at the watershed scale. In the present work, a comprehensive investigation of nitrogen sources and sinks in agricultural soils and waters was performed with the aim of gaining insights into the relevance of soil denitrification in a highly farmed sub-basin of the Po River delta (Northern Italy). Agricultural statistics, water quality datasets, and results of laboratory experiments targeting nitrogen fluxes in soils were combined to set up a detailed nitrogen budget along the terrestrial–freshwater continuum. The soil nitrogen budget was not closed, with inputs exceeding outputs by 72 kg N·ha−1·year−1, highlighting a potential high risk of nitrate contamination. However, extensive monitoring showed a general scarcity of mineral nitrogen forms in both shallow aquifers and soils. The present study confirmed the importance of denitrification, representing ~37% of the total nitrogen inputs, as the leading process of nitrate removal in heavily fertilized fine-texture soils prone to waterlogged conditions.

Denitrification and benthic metabolism in lowland pit lakes: The role of trophic conditions

Over recent decades, a great number of pit lakes have been formed, as a result of sand and gravel quarrying in river floodplains that are often also heavily exploited for agriculture. These lakes can act as nutrient filters and regulate the nitrogen pollution resulting from agricultural fertiliser use. In this paper we report the main outcomes of a study of the major nitrogen pathways in five pit lakes of differing trophic status, located along a lowland stretch of the Po river (Northern Italy). Benthic nitrogen fluxes and denitrification rates were determined in the hypolimnion and denitrification and reactive nitrogen assimilation by microphytobenthos in the littoral zone. We tested the hypothesis that lake depth and trophic status can impair denitrification and/or reactive nitrogen assimilation, compromising the function of the lakes as nutrient filters. In the studied lakes, denitrification and reactive nitrogen assimilation by primary producer communities accounted for substantial nitrogen removal rates, which were among the highest reported in the literature. Benthic nitrogen fluxes and denitrification varied between and within lakes, with depth. The littoral zone and surface waters also supported primary production, favouring nitrogen assimilation and temporal retention in the primary producer biomass. In all lakes, denitrification rates decreased from littoral to hypolimnetic sites. Denitrification rates and net nitrogen assimilation also diminished from oligotrophic to eutrophic conditions. To some extent, in eutrophic lakes there was a transfer of primary production from the benthos to the water column and the benthic system became heterotrophic, reducing the capacity for net nitrogen removal. Overall these results highlight that floodplain pit lakes can provide ecosystem services formerly supplied by natural wetlands. An important factor for management is the development of extensive littoral and shallow water zones, which are critical for maximising the nitrogen removal.

Feces from Piscivorous and Herbivorous Birds Stimulate Differentially Phytoplankton Growth

Aquatic birds may impact shallow ecosystems via organic and nutrient enrichment with feces. Such input may alleviate nutrient limitation, unbalance their ecological stoichiometry, and stimulate primary production. Herbivorous and piscivorous birds may produce different effects on aquatic ecosystems due to different physiology, diet and feces elemental composition. We analyze the effects of droppings from swans (herbivorous) and cormorants (piscivorous) on phytoplankton growth via a laboratory experiment. These birds are well represented in the Curonian Lagoon, where they form large colonies. As this lagoon displays summer algal hyper-blooms, we hypothesize an active, direct role of birds via defecation on algal growth. Short-term incubations of phytoplankton under low and high feces addition produces different stimulation of algal growth, significantly higher with high inputs of cormorant feces. The latter produces a major effect on reactive phosphorus concentration that augments significantly, as compared to treatments with swan feces, and determines an unbalanced, N-limited stoichiometry along with the duration of the experiment. During the incubation period, the dominant algal groups switch from blue-green to green algae, but such switch is independent of the level of feces input and from their origin. Heterotrophic bacteria also are stimulated by feces addition, but their increase is transient.

Investigation of watershed nutrient export affected by extreme events and the corresponding sampling frequency

Although the real-time monitoring technique has been widely applied due to the improvement of sensors, development of traditional sampling methods is still worth of being discussed due to the economically feasibility. Currently, extreme events (e.g. extreme rainfall caused by climate change) play a relatively important role in nutrient export. However, impacts of extreme events on the optimization of sampling strategy is still not well addressed despite the uncertainty of different frequency sampling programs has been sufficiently discussed in previous studies. Therefore, the corresponding impact of extreme events impact on the optimization of sampling strategy was investigated by examining temporal (i.e., inter-annual and seasonal) variations of available data. Uncertainty of nutrient flux estimates under different sampling frequencies was explored by subsampling daily monitoring data. Results showed that uncertainty in flux estimates differed between nitrogen and phosphorus. The relative error (RE) in annual TN flux estimates ranged from -4.2% to 2.4% (once per three-day) to -21.4-31.1% (monthly sampling), while the RE in annual TP flux estimates varied from -14.1% to 8.2% (once per three-day) to -65.9%-163.4% (monthly sampling). Biweekly and weekly sampling routines are considered the optimal sampling program for total nitrogen (TN) and for total phosphorus (TP) when the extreme events impact were not been considered. The uncertainty of flux estimates with different sampling frequencies increased with the increasing extreme events. High level of uncertainty occurred in years with the most extreme events in 2012 (RE: 21.4-69.0% for TN, 33.3-96.6% for TP), while the lowest can be found in 2011 (RE: 0-20.7% for TN, 0-48.3% for TP) (with the fewest extreme events). In addition, uncertainty in TN and TP flux estimates was generally greater during summer season than during other seasons. These results highlighted the important role of extreme events in nutrient export. Approximately half of the annual TN and TP flux occurred in some extreme days that only accounted for less than 20% in the same year. The onset of these extremes of nutrient export was likely due to the stormflow with addition of external fertilizer and the direct discharge of surface ponding water from paddy fields during special periods of time. These results would be helpful for the optimization of sampling strategy.

Is Flood Irrigation a Potential Driver of River-Groundwater Interactions and Diffuse Nitrate Pollution in Agricultural Watersheds?

In the Po plain, northern Italy, rivers within agricultural basins display steep summer increases in nitrate (NO3−) concentrations. Flood irrigation in overfertilized, permeable soils may drive such diffuse pollution, facilitating interactions between NO3−-rich groundwater and surface waters. We discuss multiple, indirect evidence of this mechanism in the Adda, Oglio, and Mincio rivers. These rivers drain agricultural soils with elevated nitrogen (N) surpluses, averaging 139, 193, and 136 kg ha−1 in the Adda, Oglio, and Mincio watersheds, respectively. The three rivers cross a transitional area between highly permeable and impermeable soils, where summer NO3− concentrations may increase by one order of magnitude over short distances (8–20 km). Upstream of this transitional area, a major fraction of the river flow is diverted for flood irrigation, a traditional and widespread irrigation technique for permeable soils. We speculate that diverted water solubilizes soil N excess, recharges the aquifer, and transfers soil N surplus into groundwater, resulting in NO3− pollution. Groundwater–river interactions were estimated experimentally, via water and NO3− budgets in 0.3 to 1 m3 s−1 km−1 and in 1500 to 5400 kg NO3−–N day−1. The data suggest a pronounced east–west gradient of groundwater to river diffuse water inputs among the three adjacent basins, reflecting the soil permeability and the width of the river–groundwater interaction zone. Given the large stock of NO3− in groundwater, management interventions performed at the basin scale and aimed at decreasing N excess will not produce an immediate decrease in river NO3− pollution.

The future depends on what we do today - Projecting Europe's surface water quality into three different future scenarios

There are infinite possible future scenarios reflecting the impacts of anthropogenic multiple stress on our planet. These impacts include changes in climate and land cover, to which aquatic ecosystems are especially vulnerable. To assess plausible developments of the future state of European surface waters, we considered two climate scenarios and three storylines describing land use, management and anthropogenic development ('Consensus', 'Techno' and 'Fragmented', which in terms of environmental protection represent best-, intermediate- and worst-case, respectively). Three lake and four river basins were selected, representing a spectrum of European conditions through a range of different human impacts and climatic, geographical and biological characteristics. Using process-based and empirical models, freshwater total nitrogen, total phosphorus and chlorophyll-a concentrations were projected for 2030 and 2060. Under current conditions, the water bodies mostly fail good ecological status. In future predictions for the Techno and Fragmented World, concentrations further increased, while concentrations generally declined for the Consensus World. Furthermore, impacts were more severe for rivers than for lakes. Main pressures identified were nutrient inputs from agriculture, land use change, inadequately managed water abstractions and climate change effects. While the basins in the Continental and Atlantic regions were primarily affected by land use changes, in the Mediterranean/Anatolian the main driver was climate change. The Boreal basins showed combined impacts of land use and climate change and clearly reflected the climate-induced future trend of agricultural activities shifting northward. The storylines showed positive effects on ecological status by classical mitigation measures in the Consensus World (e.g. riparian shading), technical improvements in the Techno World (e.g. increasing wastewater treatment efficiency) and agricultural extensification in the Fragmented World. Results emphasize the need for implementing targeted measures to reduce anthropogenic impacts and the importance of having differing levels of ambition for improving the future status of water bodies depending on the societal future to be expected.

An ounce of prevention is worth a pound of cure: Managing macrophytes for nitrate mitigation in irrigated agricultural watersheds

Although ubiquitous elements of agricultural landscapes, the interest on ditches and canals as effective filters to buffer nitrate pollution has been raised only recently. The aim of the present study was to investigate the importance of in-ditch denitrification supported by emergent aquatic vegetation in the context of N budget in agricultural lands of a worldwide hotspot of nitrate contamination and eutrophication, i.e. the lowlands of the Po River basin (Northern Italy). The effectiveness of N abatement in the ditch network (>18,500 km) was evaluated by extrapolating up to the watershed reach-scale denitrification rates measured in a wide range of environmental conditions. Scenarios of variable extents of vegetation maintenance were simulated (25%, 50% and 90%), and compared to the current situation when the natural development occurs in only 5% of the ditch network length, subjected to mechanical mowing in summer. Along the typical range of nitrate availability in the Po River lowlands waterways (0.5-8 mg N L^-1), the current N removal performed by the ditch network was estimated in 3300-4900 t N yr^-1, accounting for at most 11% of the N excess from agriculture. The predicted nitrate mitigation potential would increase up to 4000-33,600 t N yr^-1 in case of vegetation maintenance in 90% of the total ditch length. Moreover, a further significant enhancement (57% on average) of this key ecosystem function would be achieved by postponing the mowing of vegetation at the end of the growing season. The simulated outcomes suggest that vegetated ditches may offer new agricultural landscape management opportunities for effectively decreasing nitrate loads in surface waters, with potential improved water quality at the watershed level and in the coastal zones. In conclusion, ditches and canals may act as metabolic regulators and providers of ecosystem services if conservative management practices of in-stream vegetation are properly implemented and coupled to hydraulic needs.

Flow Regime and Nutrient-Loading Trends from the Largest South European Watersheds: Implications for the Productivity of Mediterranean and Black Sea’s Coastal Areas

In the last century, large watersheds in Southern Europe have been impacted by a combination of anthropogenic and climatic pressures, which have rapidly evolved to change the ecological status of freshwater and coastal systems. A comparative analysis was performed for Ebro, Rhône, Po and Danube rivers, to investigate if they exhibited differential dynamics in hydrology and water quality that can be linked to specific human and natural forces acting at sub-continental scales. Flow regime series were analyzed from daily to multi-decadal scales, considering frequency distributions, trends (Mann–Kendall and Sen tests) and discontinuities (SRSD Method). River loads of suspended matter, nutrients and organic matter and the eutrophication potential of river nutrients were estimated to assess the impact of river loads on adjacent coastal areas. The decline of freshwater resources largely impacted the Ebro watershed on annual (−0.139 km3 yr−1) and seasonal (−0.4% yr−1) scales. In the other rivers, only spring–summer showed significant decreases of the runoff coupled to an exacerbated flow variability (0.1–0.3% yr−1), which suggested the presence of an enhanced regional climatic instability. Discontinuities in annual runoff series (every 20–30 years) indicated a similar long-term evolution of Rhône and Po rivers, differently from Ebro and Danube. Higher nutrient concentrations in the Ebro and Po (+50%) compared to Rhône and Danube and distinct stoichiometric nutrient ratios may exert specific impacts on the growth of plankton biomass in coastal areas. The overall decline of inorganic phosphorus in the Rhône and Po (since the 1980s) and the Ebro and Danube (since the 1990s) mitigated the eutrophication in coastal ecosystems inducing, however, a phase in which the role of organic phosphorus loads (Po > Danube > Rhône > Ebro) on coastal productivity could be more relevant. Overall, the study showed that the largest South European watersheds are differently impacted by anthropogenic and climatic forces and that this will influence their vulnerability to future changes of flow regime and water quality.