The effects of irrigation on groundwater quality and quantity in a human-modified hydro-system: The Oglio River basin, Po Plain, northern Italy

Impacts of large-scale irrigation and climate change on groundwater quality and the hydrological cycle: A case study of the Alqueva irrigation scheme and the Gabros de Beja aquifer system

This research aims to analyze the impacts of the large-scale Alqueva Irrigation System (AIS) on the water cycle in selected sub-basins and the underlying Gabros de Beja aquifer system (GBAS) in Southern Portugal. The Alqueva reservoir and irrigation project is one of the largest strategic water reservoirs in Western Europe and the AIS's primary source. The closure of the dam in 2002 resulted in significant changes to the region's land use and agricultural practices, shifting from predominantly rainfed dry cereals to intensively irrigated olive and almond orchards. Therefore, this study used SWAT+ to simulate water flows from 1934 to 2021 and examined the evolution of groundwater quality and its correlation with irrigation, using data from about 50 wells from 2002 to 2021. Kriging spatial interpolation, Mann-Kendal and Sen's trend tests and the correlation technique were used. The findings revealed several noteworthy trends. First, there was a significant historical decrease in precipitation, which can be attributed to climate change. The simulation indicated a decrease in runoff and recharge, along with an increase in actual evapotranspiration due to irrigation. Furthermore, the hydrogeostatiscal analysis showed that the aquifer experienced significant salinization after the AIS implementation. In contrast, a preponderant decreasing trend in nitrate concentration was observed, which may be attributed to (a) a decrease in fertilizer use, (b) dilution effects and (c) alteration in nitrates pathways due to changes in crop types. Finally, the correlation analysis suggested that nitrates and chlorides were highly correlated to actual evapotranspiration and precipitation evolution, which may be associated to irrigation. In conclusion, the large-scale irrigation implementation and climate change significantly altered the water cycle of the study region. Overall, these findings addressed existing knowledge gaps and provided valuable insights that can be extrapolated to draw conclusions and generalize climate change and irrigation's effects on fluvial ecosystems.

Hydrochemical variations and driving mechanisms in a large linked river-irrigation-lake system

A linked river-irrigation-lake system exhibits intricate and dynamic hydrochemical variations, closely related to changes in natural conditions and anthropogenic activities. However, little is known about the sources, migration and transformation of hydrochemical composition, and the driving mechanisms, in such systems. In this study, the hydrochemical characteristics and processes in the linked Yellow River-Hetao Irrigation District-Lake Ulansuhai system were studied, based on a comprehensive hydrochemical and stable isotope analysis of water samples collected during spring, summer, and autumn. The results showed that the water bodies in the system were weakly alkaline with a pH range of 8.05-8.49. The concentrations of hydrochemical ions showed an increasing trend in the water flow direction. Total dissolved solids (TDS) were less than 1000 mg/L (freshwater) in the Yellow River and the irrigation canals, and increased to more than 1800 mg/L (saltwater) in the drainage ditches and Lake Ulansuhai. The dominant hydrochemical types varied from SO₄•Cl-Ca•Mg and HCO₃-Ca•Mg types in the Yellow River and the irrigation canals to Cl-Na type in the drainage ditches and Lake Ulansuhai. The ion concentrations in the Yellow River, the irrigation canals, and the drainage ditches were highest during summer, while ion concentrations in Lake Ulansuhai were highest during spring. The hydrochemistry of the Yellow River and the irrigation canals was mainly affected by rock weathering, while evaporation was the principal controlling factor in the drainage ditches and Lake Ulansuhai. Water-rock interactions including the dissolution of evaporites and silicates, the precipitation of carbonates, and cation exchange were the main sources of hydrochemical compositions in this system. Anthropogenic inputs had a low impact on the hydrochemistry. Therefore, greater attention should be paid in future to hydrochemical variations, especially salt ions, in the management of linked river-irrigation-lake system water resources.

Redistribution effect of irrigation on shallow groundwater recharge source contributions in an arid agricultural region

Recharge sources such as precipitation, mountain front recharge, mountain block recharge and confined water are the sources usually considered in quantitative studies of groundwater recharge. Changes in recharge processes caused by irrigation practices need to be fully considered for the accurate budgeting and management of water resources. Here, we put forward a conceptual framework for evaluating the shallow groundwater recharge process in arid irrigated agricultural areas using hydrochemical and stable isotope techniques, combined with an assessment of hydrogeological conditions and quantitative models. In general, the recharge effect of atmospheric precipitation on shallow groundwater in arid areas is relatively small. The contributions made by recharge sources in the studied river irrigated area, from greater to smaller, were confined groundwater (46.98 %), river water (45.48 %) and precipitation (7.55 %). The original range in groundwater recharge levels caused by river leakage also appeared to have expanded in response to the establishment of canal irrigation networks. Lateral groundwater flow and confined groundwater were the main recharge sources of shallow groundwater in areas fed by well irrigation and well-/spring-water irrigation (not taking into account any groundwater irrigation leakage). However, had the recharge of shallow groundwater by groundwater irrigation leakage, which reached 19.8-41.1 %, not been counted as contributing to actual groundwater recharge, the recharge contributions made by lateral groundwater flow and confined groundwater to shallow groundwater would have been significantly overestimated. This is because the groundwater recharge process has been modified by the various irrigation measures employed in arid agricultural areas, leading to a redistribution effect in groundwater recharge source contributions. This study provides a new perspective and intuitive data support for the development and utilization of water resources in arid regions.

Multi-decadal analysis of water resources and agricultural change in a Mediterranean semiarid irrigated piedmont under water scarcity and human interaction

Mediterranean piedmonts are an important hydro-agricultural systems. They constitute the junction between the mountains, where the streamflow is generated, and the surrounding plains, where the water is used. In Morocco, these traditional systems extend largely along the High Atlas Mountains. Yet, changing conditions in the Mediterranean basin as well as recurrent droughts in recent decades remain poorly understood in terms of hydrological and agricultural impacts, particularly in traditional hydro-agro-systems. The combined effects of climate variability and ineffective management of water resource dynamics may lead to increased water scarcity in these regions. The present work aims to assess the effects of climate variability and associated agricultural changes on water resources in a traditional irrigated piedmont of the Moroccan High Atlas. To that end, a trend analysis, together with change points detection, was carried out on annual and monthly precipitation, and streamflow from 1965 to 2018. Then, the standardized precipitation index (SPI) was employed to identify meteorological droughts. Also, groundwater, and spring discharge data were analyzed and discussed from 1973 to 2021. SPI outcomes revealed three major droughts, in 1981-1988, 1999-2008, and 2013-2018. Although the precipitation data showed no significant trend, except for Tahannaout station, the average annual precipitation over the piedmont area decreased by 28%. Similarly, streamflow decreased significantly by almost 40% for all stations, as did the Abainou spring's discharge. Consistent with that, groundwater level has declined dramatically over the past decades in the downstream piedmont. These decreases in water cycle components were tightly aligned with droughts. Yet, irrigation diversions were maintained in both dry and wet periods. Paradoxically, this decrease in water resources was associated with an agricultural transition from seasonal crops (cereals) to perennial crops (olive trees). This conversion is likely to amplify the water shortage, leading to groundwater resources overexploitation to satisfy the growing agricultural demand.

A Water Balancing Act: Water Balances Highlight the Benefits of Community-Based Adaptive Management in Northern New Mexico, USA

Quantifying groundwater recharge from irrigation in water-scarce regions is critical for sustainable water management in an era of decreasing surface water deliveries and increasing reliance on groundwater pumping. Through a water balance approach, our study estimated deep percolation (DP) and characterized surface water and groundwater interactions of two flood-irrigated fields in northern New Mexico to evaluate the regional importance of irrigation-related recharge in the context of climate change. DP was estimated for each irrigation event from precipitation, irrigation input, runoff, change in soil water storage, and evapotranspiration data for both fields. Both fields exhibited positive, statistically significant relationships between DP and total water applied (TWA), where one field exhibited positive, statistically significant relationships between DP and groundwater level fluctuation (GWLF) and between GWLF and total water applied. In 2021, total DP on Field 1 was 739 mm, where 68% of irrigation water applied contributed to DP. Field 2′s total DP was 1249 mm, where 81% of irrigation water applied contributed to DP. Results from this study combined with long-term research indicate that the groundwater recharge and flexible management associated with traditional, community-based irrigation systems are the exact benefits needed for appropriate climate change adaptation.

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.

Linking local natural background levels in groundwater to their generating hydrogeochemical processes in Quaternary alluvial aquifers

Calculating natural background levels (NBLs) in groundwater is vital for supporting a sustainable use of groundwater resources. Although NBLs are often assessed through a unique concentration value per groundwater body, where hydrogeochemical features are highly variable, spatial heterogeneity needs to be accounted for, leading to the calculation of so-called "local" NBLs. Despite much research devoted to the identification of the best performing techniques for local NBLs spatialization, a deep understanding of the link between local NBL values and their generating hydrogeochemical processes is often lacking and so is addressed here for the redox-sensitive species As, NH₄, Fe and Mn in the groundwater bodies of Lombardy region, N Italy. Local NBLs were calculated by a tired approach involving the hybridization of preselection and probability plot methods. Since the spatial variability of the target species depends mainly on redox conditions, a redox zonation was performed using multivariate statistical analysis. A conceptual model was developed and improved combing factor and cluster analysis. Results showed that NBLs for arsenic were up to 291 μg/L, reached in groundwaters under methanogenesis, a condition related to the prolonged degradation of peat buried in aquifer sediments. Ammonium NBLs up to 6.62 mg/L were generated by the upwelling of fluids from deep sediments hosting petroleum systems; ammonium NBLs up to 4.48 mg/L were generated as the accumulation of by-products of peat degradation. Iron and manganese NBLs up to, respectively, 6.0 and 1.51 mg/L were generated by the oxidation of younger and less stable Mn and Fe oxides within river valleys, mostly the Po River valley. The evaluation of local NBLs, and their association to generating natural hydrogeochemical processes/conditions, achieves a step forward from the commonly used approach of a single NBL per groundwater body, improving decision-support tools for sustainable groundwater management and protection.

Spatial Analysis of the Chemical Regime of Groundwater in the Karatal Irrigation Massif in South-Eastern Kazakhstan

Understanding the groundwater chemical regime on irrigated lands is relevant for assessing hydrogeological and reclamation processes. Water chemistry, natural and climatic conditions, and drainage characteristics define the complex hydrochemical connections which are formed between irrigation waters, soil, and water of the unsaturated zone and groundwater. This research paper presents the results of studies of seasonal changes in the chemical composition of groundwater in the Karatal irrigated massif. A detailed analysis of the salt regime and chemical composition of groundwater have been performed based on a large number of samples for 2019. Besides, data related to the Karatal River and drainage waters quality inside the basin have also been analyzed. Results demonstrate that the salinization processes in the massif are determined by hydrogeological conditions and drainage availability. In the north-eastern part of the massif, a high risk of soil salinization has been observed as a consequence of all the salinization processes which are determined by its hydrogeological conditions and drainage availability. Data analysis over soil samples and irrigation and drainage waters show that the hydrochemical regime of groundwater on Karatal irrigation massif is influenced by the presence of saline soil-forming rocks, the groundwater depth, the volume, and the quality of irrigation waters, the characteristics of the natural drainage, and evaporation processes.

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.

Unsupervised Machine Learning and Data Mining Procedures Reveal Short Term, Climate Driven Patterns Linking Physico-Chemical Features and Zooplankton Diversity in Small Ponds

Machine Learning (ML) is an increasingly accessible discipline in computer science that develops dynamic algorithms capable of data-driven decisions and whose use in ecology is growing. Fuzzy sets are suitable descriptors of ecological communities as compared to other standard algorithms and allow the description of decisions that include elements of uncertainty and vagueness. However, fuzzy sets are scarcely applied in ecology. In this work, an unsupervised machine learning algorithm, fuzzy c-means and association rules mining were applied to assess the factors influencing the assemblage composition and distribution patterns of 12 zooplankton taxa in 24 shallow ponds in northern Italy. The fuzzy c-means algorithm was implemented to classify the ponds in terms of taxa they support, and to identify the influence of chemical and physical environmental features on the assemblage patterns. Data retrieved during 2014 and 2015 were compared, taking into account that 2014 late spring and summer air temperatures were much lower than historical records, whereas 2015 mean monthly air temperatures were much warmer than historical averages. In both years, fuzzy c-means show a strong clustering of ponds in two groups, contrasting sites characterized by different physico-chemical and biological features. Climatic anomalies, affecting the temperature regime, together with the main water supply to shallow ponds (e.g., surface runoff vs. groundwater) represent disturbance factors producing large interannual differences in the chemistry, biology and short-term dynamic of small aquatic ecosystems. Unsupervised machine learning algorithms and fuzzy sets may help in catching such apparently erratic differences.

Overlapping redox zones control arsenic pollution in Pleistocene multi-layer aquifers, the Po Plain (Italy)

Understanding the factors that control As concentrations in groundwater is vital for supplying safe groundwater in regions with As-polluted aquifers. Despite much research, mainly addressing Holocene aquifers hosting young (<100 yrs) groundwater, the source, transport, and fate of As in Pleistocene aquifers with fossil (>12,000 yrs) groundwaters are not yet fully understood and so are assessed here through an evaluation of the redox properties of the system in a type locality, the Po Plain (Italy). Analyses of redox-sensitive species and major ions on 22 groundwater samples from the Pleistocene arsenic-affected aquifer in the Po Plain shows that groundwater concentrations of As are controlled by the simultaneous operation of several terminal electron accepters. Organic matter, present as peat, is abundant in the aquifer, allowing groundwater to reach a quasi-steady-state of highly reducing conditions close to thermodynamic equilibrium. In this system, simultaneous reduction of Fe-oxide and sulfate results in low concentrations of As (median 7 μg/L) whereas As reaches higher concentrations (median of 82 μg/L) during simultaneous methanogenesis and Fe-reduction. The position of well-screens is an additional controlling factor on groundwater As: short screens that overlap confining aquitards generate higher As concentrations than long screens placed away from them. A conceptual model for groundwater As, applicable worldwide in other Pleistocene aquifers with reducible Fe-oxides and abundant organic matter is proposed: As may have two concentration peaks, the first after prolonged Fe-oxide reduction and until sulfate reduction takes place, the second during simultaneous Fe-reduction and methanogenesis.

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.

Evaluation of Seasonal Groundwater Quality Changes Associated with Groundwater Pumping and Level Fluctuations in an Agricultural Area, Korea

This study was conducted to evaluate seasonal groundwater quality due to groundwater pumping and hydrochemical characteristics with groundwater level fluctuations in an agricultural area in Korea. Groundwater levels were observed for about one year using automatic monitoring sensors, and groundwater uses were estimated based on the monitoring data. Groundwater use in the area is closely related to irrigation for rice farming, and rising groundwater levels occur during the pumping, which may be caused by the irrigation water of rice paddies. Hydrochemical analysis results for two separate times (17 July and 1 October 2019) show that the dissolved components in groundwater decreased overall due to dilution, especially at wells in the alluvial aquifer and shallow depth. More than 50% of the samples were classified as CaHCO3 water type, and changes in water type occurred depending on the well location. Water quality changes were small at most wells, but changes at some wells were evident. In addition, the groundwater quality was confirmed to have the effect of saltwater supplied during the 2018 drought by comparison with seawater. According to principal component analysis (PCA), the water quality from July to October was confirmed to have changed due to dilution, and the effect was strong at shallow wells. In the study areas where rice paddy farming is active in summer, irrigation water may be one of the important factors changing the groundwater quality. These results provide a qualitative and quantitative basis for groundwater quality change in agricultural areas, particularly rice paddies areas, along with groundwater level and usage.

Hydrochemical composition and potentially toxic elements in the Kyrgyzstan portion of the transboundary Chu-Talas river basin, Central Asia

Water chemistry and the assessment of health risks of potentially toxic elements have important research significance for water resource utilization and human health. However, not enough attention has been paid to the study of surface water environments in many parts of Central Asia. Sixty water samples were collected from the transboundary river basin of Chu-Talas during periods of high and low river flow, and the hydrochemical composition, including major ions and potentially toxic elements (Zn, Pb, Cu, Cr, and As), was used to determine the status of irrigation suitability and risks to human health. The results suggest that major ions in river water throughout the entire basin are mainly affected by water-rock interactions, resulting in the dissolution and weathering of carbonate and silicate rocks. The concentrations of major ions change to some extent with different hydrological periods; however, the hydrochemical type of calcium carbonate remains unchanged. Based on the water-quality assessment, river water in the basin is classified as excellent/good for irrigation. The relationship between potentially toxic elements (Zn, Pb, Cu, Cr, and As) and major ions is basically the same between periods of high and low river flow. There are significant differences between the sources of potentially toxic elements (Zn, Pb, Cu, and As) and major ions; however, Cr may share the same rock source as major ions. The risk assessment revealed low non-carcinogenic and carcinogenic risks for human health; however, the maximum carcinogenic risk for As exceeded the allowable value, which requires further consideration. These results provide a scientific basis for the management of agricultural irrigation uses and also infill existing gaps regarding the hydrochemical composition in the Chu-Talas river basin, Central Asia.

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.

Heavy Metals and Related Human Health Risk Assessment for River Waters in the Issyk-Kul Basin, Kyrgyzstan, Central Asia

The water resources of Central Asia play an important role in maintaining the fragile balance of ecosystems and the sustainable development of human society. However, the lack of research on the heavy metals in river waters has a far−reaching influence on public health and the sustainable development in Central Asia. In order to reveal the possible sources of the heavy metals and to assess the associated human health risks, thirty−eight water samples were collected from the rivers of the Issyk−Kul Basin during the period with low river flow (May) and the period with high river flow (July and August), and the hydrochemical compositions and major ions of heavy metals were analyzed. No changes in hydrochemical facies were observed between the two periods and the river water type was calcium bicarbonate. Carbonate dissolution and silicate weathering controlled the variation of cations and anions in river waters from the Issyk−Kul Basin. There were some differences in the sources of heavy metals in water bodies between the two periods. During the period with low river flow, heavy metals (Cr) were closely clustered with major ions, indicating that they were mainly affected by water–rock interactions. During the period with high river flow, all heavy metals studied in this paper had different sources of major ions, and the heavy metals maybe influenced by human activities. From the human health risk assessment, the hazard quotients for all samples were less than 1, reflecting that there was no noncarcinogenic risk in the river waters of the Issyk−Kul Basin during the two sampling periods. However, the water samples with carcinogenic risk of arsenic exceeding the threshold (10⁻⁴) accounted for 21.1% of the total, indicating that there were some certain carcinogenic hazards for human health via water drinking with direct oral ingestion. The results are of certain significance for the utilization and protection of water resources in the basin as well as the protection of public health.

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.

Groundwater and surface water quality characterization through positive matrix factorization combined with GIS approach

This study aims at testing the effectiveness of Positive Matrix Factorization in characterizing groundwater and surface water quality, in terms of identifying main hydrochemical features and processes (natural and anthropogenic) that govern them. This method is applied in a hydro-system featured by a strong interrelation between groundwater and surface water and highly impacted by agricultural activities. Therefore, a holistic approach considering groundwater together with the surface water bodies, consisting in lake, several rivers and springs, was used. Multivariate statistical analysis, in particular Factor Analysis, has been proved to be effective in elaborating and interpreting water quality data highlighting the information carried within them, but it presents some limitations: it does not consider data uncertainty and it groups variables which are correlated positively and negatively. Moreover, in some cases the resulting factors are not clearly interpretable, describing each one various overlapping features/processes. Here, Positive Matrix Factorization is applied to groundwater and surface water quality data, and the results are compared to those obtained through a Factor Analysis in terms of both factor profiles and their spatial distribution through a GIS approach. Results of isotopes analysis are used to validate PMF output and support interpretation. Positive Matrix Factorization allows to consider data uncertainty and the solution respects two positivity constraints, based on the concept of chemical mass balance, which leads to a more environmentally interpretable solution. Results show that Positive Matrix Factorization identifies five different factors reflecting main features and natural and anthropogenic processes affecting the study area: 1) surface water used for irrigation, 2) groundwater subjected to reducing processes at advanced stages, 3) groundwater subjected to reducing processes at early stages, 4) groundwater residence time and 5) the effects of the agricultural land use on both groundwater and surface water. Positive Matrix Factorization leads to a more detailed understanding of the studied system as compared to Factor Analysis which identifies only three factors with overlapping information. Based on the results of this study, Positive Matrix Factorization could be a useful technique to perform groundwater and surface water quality characterization and to reach a deeper understanding of the phenomena that govern water chemistry.