Groundwater quality trend and trend reversal assessment in the European Water Framework Directive context: an example with nitrates in Italy.

Groundwater resources are of utmost importance in sustaining water related ecosystems, including humans. The long-lasting impacts from anthropogenic activities require early actions, owing to the natural time lag in groundwater formation and renewal. The European Union (EU) policy, within the implementation of the Water Framework Directive (WFD), requires Member States to identify and reverse any significant and sustained upward trend in the concentration of pollutants, defining specific protection measures to be included in the River Basin Management Plans (RBMP). In Italy, official guidelines for trend and trend reversal assessment have been published recently. Statistical methods, such as the Mann-Kendall test for trend analysis and the Sen's method for estimating concentration scenarios, should be applied at the fixed terms stated by the WFD implementation cycles to identify upward trends, while the Pettitt test is proposed for the identification of trend reversal. In this paper, we present an application of a slightly modified version of the Italian Guidelines to a groundwater body in Northern Italy featuring nitrate pollution and discuss its advantages and limitations. In addition to Pettitt test, for the trend reversal analysis, we apply the Mann-Kendall test in two sections and compare the results. We conclude that this method seems more reliable than Pettitt test to identify a reversal point in quality time series. The overall procedure can be easily applied to any groundwater body defined at risk across Europe, for the assessment of the upward trends of pollutants and their reversal, even with little chemical monitoring data. Although focused on the EU legislative framework, this procedure may be relevant for a wider context, allowing to individuate upward trend as early warning for contamination processes in an integrated water resources management context.

Generation and dispersion of chemical and biological gradients in a large-deep multi-basin lake (Lake Como, north Italy): The joint effect of external drivers and internal wave motions.

The formation of chemical and biological heterogeneity in lakes can be favored by physical and morphometric constrains. This study describes the results of four whole-lake field campaigns carried out in Lake Como (north Italy) during thermal stratification. The aim was to analyze the distribution of chemical-biological variables in a multi-basin lake as a result of internal and external physical drivers and constrains. Lake Como has a y-like shape encompassing three main sub-basins: northern, south-eastern, and south-western. Field data underlined: the presence of chemical-biological gradients between the south-western basin and the rest of the lake and the propagation of a fresher water-plume (formed by the two main northern inflows) into the northern basin and then into the south-western closed arm. The use of a three-dimensional hydrodynamic model showed a periodic movement of this plume that tends to enter the south-western basin first and then to return toward north, moving forward and back through the junction of the three arms. The entrance into the eastern basin, instead, occurs only secondarily. Wavelet transform analysis revealed common periodicity between the plume movement and the action of different external and internal lake-stressors, including: the discharge of the main inflow (period centered at 1, 3.3, and 6.5 day), the wind intensity (0.5 and 1 day) and the two main basin-scale internal wave motions: (3.3 day and 7.1 day). The periodic movement of the fresher water-plume enhances the water exchange and reduces the chemical-biological gradients between the western closed basin and the main lake, playing a crucial role in distributing both inorganic and organic materials at the lake basin-scale.

Pore-scale simulations of concentration tails in heterogeneous porous media.

The retention of contaminants in the finest and less-conductive regions of natural aquifer is known to strongly affect the decontamination of polluted aquifers. In fact, contaminant transfer from low to high mobility regions at the back end of a contaminant plume (i.e. back diffusion) is responsible for the long-term release of contaminants during remediation operation. In this paper, we perform pore-scale calculations for the transport of contaminant through heterogeneous porous media composed of low and high mobility regions with two objectives: (i) study the effect of permeability contrast and solute transport conditions on the exchange of solutes between mobile and immobile regions and (ii) estimate the mass of contaminants sequestered in low mobility regions based on concentration breakthrough curves.

Long-term climate sensitivity of an integrated water supply system: The role of irrigation.

The assessment of the impact of long-term climate variability on water supply systems depends not only on possible variations of the resources availability, but also on the variation of the demand. In this framework, a robust estimation of direct (climate induced) and indirect (anthropogenically induced) effects of climate change is mandatory to design mitigation measures, especially in those regions of the planet where the groundwater equilibrium is strongly perturbed by exploitations for irrigation purposes. The main goal of this contribution is to propose a comprehensive model that integrates distributed crop water requirements with surface and groundwater mass balance, able to consider management rules of the water supply system. The proposed overall model, implemented, calibrated and validated for the case study of the Fortore water supply system (Apulia region, South Italy), permits to simulate the conjunctive use of the water from a surface artificial reservoir and from groundwater. The relative contributions of groundwater recharges and withdrawals to the aquifer stress have been evaluated under different climate perturbations, with emphasis on irrigation practices. Results point out that irrigated agriculture primarily affects groundwater discharge, indicating that ecosystem services connected to river base flow are particularly exposed to climate variation in irrigated areas. Moreover, findings show that the recharge both to surface and to groundwater is mainly affected by drier climate conditions, while hotter conditions have a major impact on the water demand. The non-linearity arising from combined drier and hotter conditions may exacerbate the aquifer stress by exposing it to massive sea-water intrusion.

Total phosphorus reference condition for subalpine lakes: a comparison among traditional methods and a new process-based watershed approach.

Different methods for estimating the total phosphorus (TP) reference conditions of lakes have rarely been compared. This work tests the uncertainty and accuracy of the most frequently used approaches (Morpho-edaphic index -MEI-, export coefficient, diatoms and pigment-inferred TP models) for 35 subalpine lakes. Furthermore, we propose a new process-based watershed approach that was tested on a subalpine environment and consists of combining a space for time substitution with a space for space substitution. The possible presence of uncontaminated or less contaminated environments inside or next to the watershed can be exploited by training a hydrological transport watershed model according to the uncontaminated conditions and then applying the calibration to the entire watershed, which reconstructs a natural or semi-natural TP load scenario. We found that the root mean square error (RMSE) for the MEI is 4 µg L(-1). However, its application is limited for lakes that present with an alkalinity ≤1 meq L(-1). For lakes with a higher alkalinity, we observed a loss of predictive capability that results from the lower solubility of phosphorus under conditions of high calcium content. The export coefficient model was applied with a mean export coefficient and presents similar prediction capabilities as the MEI. The chlorophyll-inferred TP model shows a higher uncertainty (RMSE = 8 µg L(-1)); however, it produced fewer underestimations and overestimations. With regards to the diatom-inferred TP model, we are only able to evaluate an uncertainty of 5 µg L(-1) at the European level. Finally, the proposed process-based watershed approach adequately predicted the reference condition of the selected lake and had an uncertainty lower than the other methods (2 µg L(-1)). We conclude by revealing the potential and limitations of this approach in the field of ecological lake modelling more and more attracted by TP pristine load inputs in studies on the effects of climate change and eutrophication of lakes.

Further characterization of the time transfer capabilities of precise point positioning (PPP): the Sliding Batch Procedure.

In recent years, many national timing laboratories have installed geodetic Global Positioning System receivers together with their traditional GPS/GLONASS Common View receivers and Two Way Satellite Time and Frequency Transfer equipment. Many of these geodetic receivers operate continuously within the International GNSS Service (IGS), and their data are regularly processed by IGS Analysis Centers. From its global network of over 350 stations and its Analysis Centers, the IGS generates precise combined GPS ephemeredes and station and satellite clock time series referred to the IGS Time Scale. A processing method called Precise Point Positioning (PPP) is in use in the geodetic community allowing precise recovery of GPS antenna position, clock phase, and atmospheric delays by taking advantage of these IGS precise products. Previous assessments, carried out at Istituto Nazionale di Ricerca Metrologica (INRiM; formerly IEN) with a PPP implementation developed at Natural Resources Canada (NRCan), showed PPP clock solutions have better stability over short/medium term than GPS CV and GPS P3 methods and significantly reduce the day-boundary discontinuities when used in multi-day continuous processing, allowing time-limited, campaign-style time-transfer experiments. This paper reports on follow-on work performed at INRiM and NRCan to further characterize and develop the PPP method for time transfer applications, using data from some of the National Metrology Institutes. We develop a processing procedure that takes advantage of the improved stability of the phase-connected multi-day PPP solutions while allowing the generation of continuous clock time series, more applicable to continuous operation/monitoring of timing equipment.