The Abelar pilot basin in Coruña (northwestern Spain) has been monitored for hydrological and hydrochemical data to assess the effects of eucalyptus plantation and manure applications on water resources, water quality, and nitrate contamination. Here, we report the machine learning analysis of hydrological and hydrochemical data from the Abelar basin. K-means cluster analysis (CA) is used to relate nitrate concentrations at the outlet of the basin with daily interflows and groundwater flows calculated with a hydrological balance. CA identifies three linearly separable clusters. Times series Gaussian process regression (TS-GPR) is employed to predict surface water nitrate concentration by incorporating hydrological variables as additional input parameters using a time series shifting. TS-GPR allows modelling nitrate concentrations based on shifted interflows and groundwater flows and chemical concentrations with R² = 0.82 and 0.80 for training and testing, respectively. Groundwater flow from five days prior to the current date, $Q_{g5}$, is the most important input parameter of the TS-GPR model. Interaction effects between the variables are found. TS-GPR validation with recent data provides results consistent with those of testing (R² = 0.85). Model inspection by permutation feature importance and partial dependence plots shows interactions between $Q_{g5}$ and Cl, and between Ca and Mg.

Groundwater flow and transport models are essential tools for assessing and quantifying the migration of organic contaminants at polluted sites. Uncertainties in the hydrodynamic and transport parameters of the aquifer have a significant effect on model predictions. Uncertainties can be quantified with advanced sensitivity methods such as Sobol’s High Dimensional Model Reduction (HDMR) and Variogram Analysis of Response Surfaces (VARS). Here we present the application of VARS and HDMR to assess the global sensitivities of the outputs of a transient groundwater flow model of the Gállego alluvial aquifer which is located downstream of the Sardas landfill in Huesca (Spain). The aquifer is subject to the tidal effects caused by the daily oscillations of the water level in the Sabiñánigo reservoir. Global sensitivities are analyzed for hydraulic heads, aquifer/reservoir fluxes, groundwater Darcy velocity, and hydraulic head calibration metrics. Input parameters include aquifer hydraulic conductivities and specific storage, aquitard vertical hydraulic conductivities, and boundary inflows and conductances. VARS, HDMR, and graphical methods agree to identify the most influential parameters, which for most of the outputs are the hydraulic conductivities of the zones closest to the landfill, the vertical hydraulic conductivity of the most permeable zones of the aquitard, and the boundary inflow coming from the landfill. The sensitivity of heads and aquifer/reservoir fluxes with respect to specific storage change with time. The aquifer/reservoir flux when the reservoir level is high shows interactions between specific storage and aquitard conductivity. VARS and HDMR parameter rankings are similar for the most influential parameters. However, there are discrepancies for the less relevant parameters. The efficiency of VARS was demonstrated by achieving stable results with a relatively small number of simulations.

The Mediterranean basin may be one of the few bioregions on earth where the steep environmental gradients and the biota are so tightly intermixed and interrelated within diverse ecological niches, ecosystems and landscapes. However, several natural or human-induced challenges threaten this unique ecosystem, such as climate change, exacerbating droughts, habitat and biodiversity losses. Nature-Based Solutions (NbS) represent a new concept, including old and new practices (i.e., green roofs), inspired by nature in harmony with humans to face contemporary societal and environmental challenges. EU-funded projects such as the DRYAD aim to bring upfront, upscale and highlight such practices in Mediterranean agroecosystems. This article presents a literature review on NbSs since 2000 with special emphasis on the Mediterranean basin.

Groundwater flow and transport models are essential tools for assessing and quantifying the mi-gration of organic contaminants at polluted sites. Uncertainties in the hydrodynamic and transport parameters of the aquifer have a significant effect on model predictions. Uncertainties can be quantified with advanced sensitivity methods such as Sobol’s High Dimensional Model Reduction (HDMR) and Variogram Analysis of Response Surfaces (VARS). Here we present the application of VARS and HDMR to assess the global sensitivities of the outputs of a transient groundwater flow model of the Gállego alluvial aquifer which is located downstream the Sardas landfill in Huesca (Spain). The aquifer is subject to the tidal effects caused by the daily oscillations of the water level in the Sabiñánigo reservoir. Global sensitivities are analysed for hydraulic heads, aquifer/reservoir fluxes, groundwater Darcy velocity and hydraulic head calibration metrics. Input parameters in-clude aquifer hydraulic conductivities and specific storage, aquitard vertical hydraulic conductiv-ities, and boundary inflows and conductances. VARS, HDMR and graphical methods agree to identify the most influential parameters which for most of the outputs are the hydraulic conduc-tivities of the zones closest to the landfill, the vertical hydraulic conductivity of the most permeable zones of the aquitard, and the boundary inflow coming from the landfill. The sensitivity of heads and aquifer/reservoir fluxes with respect to specific storage change with time. The aquifer/reservoir flux when the reservoir level is high shows interactions between specific storage and aquitard conductivity. VARS and HDMR parameter rankings are similar for the most influential parameters. However, there are discrepancies for the less relevant parameters. The efficiency of VARS was demonstrated by achieving stable results with a relatively small number of simulations.