Search

Freshwater in coastal and island aquifers is a valuable resource whose availability is strongly conditioned by heterogeneity. More than 80 % of the Earth’s surface is of volcanic origin, but the effect of volcanic dykes on the geometry of the saline interface that separates freshwater from seawater is still underexplored. This paper analyzes the impact of volcanic dykes on the depth of the saline interface in coastal and island aquifers and, subsequently, on the availability of fresh groundwater. Hydrogeological and hydrochemical data from a gallery, perpendicularly crossing several tens of dykes, were integrated with numerical modeling on the volcanic island of El Hierro (Canary Islands, Spain). Measured hydraulic heads demonstrated that the presence of dykes increased the hydraulic gradient by more than an order of magnitude, with respect to an adjacent area not affected by dykes. Numerical assessment confirmed that the lower the hydraulic conductivity of the dykes, the greater the depth of the saline interface inland. This impact led to fresh groundwater reserves increasing inland, relative to a hypothetical case without dykes. Numerical simulations also demonstrated that dykes can prevent salinization of production wells in coastal and island aquifers, if they are correctly located. Locating production wells far enough inland in an area affected by dykes allowed a higher freshwater extraction rate than if dykes did not exist; near the coastline, the effect tended to be the opposite. These results will be key to improving the management of fresh groundwater resources in coastal volcanic aquifers, and especially on volcanic islands such as the Hawaiian Islands or the Macaronesian archipelagos.

The 2021 volcanic eruption in the Cumbre Vieja mountain range on La Palma Island (Canary Islands, Spain) raised concerns regarding the potential impact on groundwater resources. This study is the first part of a series of papers investigating those impacts, and focuses on the hydraulic impacts of the eruption, while subsequent papers will explore the geochemical consequences. Three boreholes equipped with sensors to measure hydraulic head, temperature, and electrical conductivity of groundwater were installed near the volcano. Monitoring started during the eruption and continued a year after it. Statistical analysis were performed to assess the relationship between the measured variables and real-time seismic-amplitude measurements (RSAM). In addition, the possibility of groundwater vaporization due to magma emergence was assessed with a groundwater flow numerical modelling of the island. Correlation coefficients were computed to assess the linear relationship between groundwater parameters and seismic signals, observing a statistically significant association, and suggesting near-instantaneous variations in parameters such as groundwater levels and EC. Different response patterns of groundwater levels were observed in recharging areas in highlands compared to discharge areas, showing an opposite correlation direction. Deduction of natural trends from the linear regression models of head and RSAM two months after the eruption revealed a more predictable impact on the groundwater system, as the hydrogeological system adjusts to the volcanic activity and its effects on the aquifer. The hydrogeologicalsimulation of the "magma pumping effect" suggested that groundwater extraction was possible, but the absence of an appropriate groundwater monitoring network made it impossible to determine the amount of water extracted from the aquifer. The uncertainty analysis showed values up to 2000 m3⋅day􀀀 1. These findings have important implications for understanding the negative impacts of volcanic eruptions on groundwater resources, highlighting the need for regular monitoring and assessment by hydrogeologists and water management professionals.

Managed aquifer recharge (MAR) and, in particular, recharge by direct injection into the aquifer through wells or boreholes allows for a series of very interesting solutions to solve various technical and environmental problems related to management of the urban water cycle. These problems include the overexploitation of water resources, marine intrusion or contamination of groundwater by irrigation returns. The island of Gran Canaria presents several of these problems; thus, the feasibility of implementing a recharge system has been studied to provide a solution to some of them, using a resource of great potential such as reclaimed water. A detailed characterization of groundwater quality in the study area was carried out, complemented by a field campaign with water sampling from the surrounding catchments, in situ analysis and subsequent laboratory analysis. Specifically, an MBR treatment with disinfection is proposed, where the final conclusions indicate that this is a technically and economically viable project, innovative in its application on islands, a priori with an acceptable productive recharge capacity, possibly scalable after the experimental phase and extrapolable to other locations with similar conditions. In addition, it presents a set of important environmental benefits with respect to conservation of and2 improvement of the state of the groundwater bodies in the studied area, as well contributing to knowledge of and research into water management in volcanic lands and islands.

Solar photovoltaic energy is increasingly positioned as an excellent renewable energy source to be implemented in energy communities. However, studies on photovoltaic potential are still required in order to enhance the clean energy transition in energy dependent systems, such as oceanic islands. This work is a study of photovoltaic potential of public building roofs in the World Heritage city of San Cristóbal de La Laguna in Tenerife (Canary Islands). The buildings selected for the study are sports centres, educational centres and civic centres, considering optimum and maximum production scenarios. The results showed the electricity generation capacity of each building, with 8811 and 861 MWh capacity for the maximum and optimum scenarios, respectively. In the optimal scenario, 56% of energy would be generated, of which 31% would be consumed on site, and the remaining 25% would be sent to the grid. The main difference is that in the case of optimal photovoltaic systems, economic profitability criteria take precedence, while in the case of maximum photovoltaic systems, the environmental benefits in terms of CO2 emissions into the atmosphere increase considerably. In the latter case, CO2 emissions are reduced by 771%, compared to 75% in the case of optimal photovoltaic systems.

Fully 3D and multidimensional geological models are key for energy, mining and water resources assessment and management, which has a great economic impact. 3D geological models provide critical information regarding geometric properties of the geological bodies of interest, such as ore deposits or reservoirs, allowing inference of volumes and availability to be exploited. However, defining a useful methodology to create a fully digital 3D geological model of volcanic islands is still very challenging due to the intrinsic complexity observed in the structure of volcanic systems. Two-dimensional geological maps, associated vertical geological cross-sectional diagrams and borehole and water-mining tunnels (galleries) logs were examined on El Hierro Island. A geological modelling procedure using the software GeoModeller obtained the first 3D geological model of El Hierro volcanic island, as a basis for volcanic formation definitions and geometric analysis. Surface, volume and center of mass estimations of the differentiated sub-units allowed a better understanding of the spatial distribution and evolution of the volcanic cycles of El Hierro Island over time. The 3D model provided geometry information correlatable with volcanic activity cycles and lava/pyroclast deposition mechanisms that ultimately defines the geometry of the volcanic units. This information provides new variables valuable for the investigation of the existing relationships between rift zones and giant lateral collapses. In addition, the obtained 3D geological model of a volcanic island will support the building of the first hydrogeological and geothermal model of El Hierro Island, allowing advancement of sustainability of isolated systems.

Carbon dioxide emissions linked to the transport sector are particularly relevant to islands. The Canary Islands have high level of tourism, with tourists who generally travel in rented vehicles on arrival in the archipelago. In addition, mobility of the local population in the islands is also always growing. Thus, transport is one of the sectors that emits the greatest amount of greenhouse gases. In this sense, the net carbon balance becomes a pioneering study in the framework of management of high-capacity road transport routes and will help plan new strategies, which pursue the neutrality and climate resilience of the road network. It will therefore be possible to design different future scenarios, according to traffic intensities detected, to mitigate their local effects by increasing green areas, which will contribute to greater absorption. This study calculated the carbon footprint linked to vehicles on the three high-capacity roads on the island of Tenerife, as well as the carbon footprint of the public lighting attached to these roads. The absorption of carbon by vegetation planted by the Road Service and the soil located on the margins of the road have also been calculated to explore the possibility of modifying the vegetation, using other types of local shrubs and trees that have a higher absorption factor. The results demonstrate that conifers have the highest absorption factor, the Canary Island pine being the best option in this case. However, it is not always possible to place trees on the roadside, and it is necessary to resort to smaller species. This study concludes that better mobility planning, an increase in electric vehicles, improved energy production from renewable sources, and the promotion and enhancement of absorption capacity of carbon sinks are key in the face of climate change. Integr Environ Assess Manag 2022;00:1-8. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental ToxicologyChemistry (SETAC).

The Canary Island archipelago holds 2.2 million inhabitants and received more than 14 million tourist visitors in 2019. Its volcanic origin and the presence of recent volcanic activity highlight its potential geothermal interest for heat production, therefore offering an attractive pathway for a renewable energy supply, not only for power generation and deep geothermal exploitation but also for low-enthalpy heating and cooling production. This work investigates nine touristic infrastructures in the Canary Islands which experienced a transition from conventional heat production systems to shallow geothermal energy (SGE) systems. The SGE systems established consisted in ultra-low temperature district heating and cooling (DHC) micro-networks involving, in most cases, small apartment blocks fed by shallow geothermal plants with an average cooling capacity of 602 kW and a heating capacity of 614 kW. The examination of this transition has exhibited an average saving of 374 MWh, 69,235 € and 280 tCO2 for each installation every year. This work provides evidence of the economical, energetical and environmental advantages of SGE technology in volcanic islands facing both and enormous heating and cooling demand due to tourism, in an energy-dependent system that also faces the challenge of the decarbonization of the heating production sector.

The Canary Islands have a water culture tied to the exploitation of their groundwater by means of wells and water galleries. However, the growth of tourism, the increase in the local population and the development of agriculture have led to the emergence of new ways of obtaining water, such as the desalination of seawater. The presence of these desalination plants covers the entire archipelago except for the island of La Palma, and sometimes they function as a complement to water needs, while in other cases they are the only source of drinking water available. To study the environmental impact of the production of drinking water through the exploitation of the aquifer and the desalination of seawater, the carbon footprint methodology was used following the guidelines of the GHG Protocol. The result has shown that seawater installations have the largest carbon footprint, mainly due to the high electricity consumption in the islands and the electricity mix of the archipelago which, as it does not rely entirely on renewable energy sources, increases CO2 emissions into the atmosphere due to the production of drinking water in the islands.

Global demands for energy efficient heating and cooling systems coupled with rising commitments toward net zero emissions is resulting in wide deployment of shallow geothermal systems, typically installed to a depth of 100 to 200m, and in the continued growth of the global ground source heat pump (GSHP) market. Ground coupled heat pump (GCHP) systems take up to 85% of the global GSHP market. With increasing deployment of GCHP systems in urban areas coping with limited regulations, there is growing potential and risk for these systems to impact the subsurface thermal regime and to interact with each other or with nearby heat‐sensitive subsurface infrastructures. In this paper, we present three numerical modelling case studies, from the UK and Canada, which examine GCHP systems’ response to perturbation of the wider hydrogeological and thermal regimes. The studies demonstrate how GCHP systems can be impacted by external influences and perturbations arising from subsurface activities that change the thermal and hydraulic regimes in the area surrounding these systems. Additional subsurface heat loads near existing schemes are found to have varied impacts on system efficiency with reduction ranging from

Aquifers constitute the main freshwater supply of oceanic islands. Maintaining groundwater quantity and quality is of critical concern as demographic and climatic changes place additional pressures on already fragile water resources systems. Islands with heterogeneous volcanic aquifers pose additional difficulties in assessing their water resources. This work proposes an approach for improving the hydrogeological characterization of heterogeneous volcanic aquifer systems by making use of recession coefficients from deep water gallery discharge. To demonstrate the usefulness of this approach, hydrographs and recession curves from groundwater discharge of 30 water galleries on La Palma (Canary Islands) were evaluated. This analysis allowed us to obtain the main hydrogeological parameters of a volcanic aquifer system, in terms of hydraulic diffusivity. A Maillet-Boussinesq model with an exponential decay law was adopted, according to field observations of drainage discharge. The alpha coefficients of recession values ranged between 10−3 and 4·10−4 day−1 and showed significant spatial correlation with insular geology. Additionally, hydraulic diffusivity values of island hydrogeological domains were obtained from recession coefficients using the Rorabaugh-Singh method. Weighted storage coefficients for volcanic materials were in the range of 3% to 7%, with an average transmissivity in the range of 15 to 150 m2·day−1. The methodology proposed has demonstrated its usefulness in coping with local uncertainty in hydraulic characterization of insular aquifers associated with volcanic heterogeneity. This is an improvement compared to standard pumping tests, thus providing hydraulic parameters prior to numerical analysis for water management planning. This research was partially supported by the European Union's Horizon 2020 research and innovation programme under grant agreement 101037424, project ARSINOE (Climate resilient regions through systemic solutions and innovations), and the Spanish Research Agency (project SAGE4CAN PID2020-114218RA-100).

The outermost regions of Europe are nine: Guadeloupe, French Guiana, Martinique, Mayotte, Reunion and Saint Martin (France), the Canary Islands (Spain), the Azores and Madeira (Portugal). These regions enrich the EU economically, culturally and geographically, hosting 80% of its biodiversity. However, due to their remoteness and other unique features, they pose challenges for their development. The particular case of the Canary Islands will be developed in the framework of the European H2020 project Arsinoe, where the hydrological cycle and agriculture in the Canary archipelago will be studied in El Hierro and La Palma. These two islands have been selected for the following reasons: i) El Hierro is a pioneer in presenting a self-sufficient energy model (La Gorona del Viento project) and is rich in groundwater, this being the most used water resource on the island; ii) La Palma has been selected because it is an island rich in groundwater (in fact, it is the only island in the Canary archipelago that does not have desalination plants) and where agriculture is very important (especially tropical crops such as banana, mango, avocado, etc.) and, due to the volcanic eruption that began in September 2021, the situation of the aquifer is uncertain, something that is worrying since La Palma depends on groundwater resources to guarantee the water demand of agriculture, local population and tourism. The effect of the eruption on the vulnerability of the aquifer of La Palma is still unknown, so it is desired to study in depth the effects on the aquifer in quantitative and qualitative terms therefore, ARSINOE will focus on the ecological transition and vulnerability of aquifers in volcanic islands and will put further efforts to the primary production including agriculture, forestry, fisheries and aquaculture, water management and clean energy infrastructure. ARSINOE will take into account the interdependence between water and agriculture. The agricultural sector is the largest water user in the Canary Islands, where wine, potatoes and tomatoes are the main exports. Therefore, greater sustainability within the water sector (through the water footprint and the carbon footprint) will positively affect the agricultural sector and, therefore, the water and energy situation of the archipelago. But aquifers of both islands are also at risk due to other circumstances, specifically those derived from climate change: greater saline intrusion (due to rising sea levels), losses in freshwater inputs (due to decreased rainfall), changes in physic-chemical conditions of all water bodies… All these effects will be studied on both islands, through the ARSINOE Project, and from a local point of view.

The use of shallow geothermal energy (SGE) resources in oceanic volcanic environments entails additional challenges when compared to continental sedimentary/plutonic settings. The efficiency of shallow geothermal heat exchangers heavily depends on the geology and hydrogeology of the terrain where are placed. Volcanic rocks in small oceanic islands (2) are the result of volcanism, erosion, and tectonic collapse. All these processes conform highly heterogeneous formations with complex hydrogeology whose thermal response to shallow geothermal systems requires a good understanding of heat transfer in such environments. The SAGE4CAN project will concentrate on SGE resource assessment taking into account heterogeneity characteristic of volcanic formations, both at local and insular scale. To this end, the Canary Islands are selected as representative volcanic oceanic islands, to define SGE implementation barriers including but not limited to (1) heterogeneities of thermal properties intrinsic to volcanic formations (volcanic dikes, red layers, landslides, etc.), (2) heat advection in the context of complex groundwater flow in the unsaturated (dominate in midlands and highlands) as well as in the saturated medium (coast), (3) enhanced geothermal gradients, (4) transient effects of urban and volcanic activity, (5) heating and cooling demand, (6) shallow geothermal energy installations design and optimization, as well as (7) energy transition strategies in energy-dependent islands. The SGE4CAN project will investigate novel approaches to overcome such boundary conditions of oceanic volcanic islands in the estimation of the renewability of the resources, developing novel procedures to conduct cost-efficient and open-access Thermal Response Tests (TRTs), investigate the performance of existent SGE systems, assessing environmental impacts associated with SGE use. The knowledge generated from this project will be used on its final stage to identify adequate strategies for the integration of SGE into heating and cooling policies and action plans, as well as to raise awareness about the technology so that it gets recognition.

One solution for reducing the current consumption of fossil fuels is a more frequent use of shallow geothermal energy. However, particularly regarding urban subsurface resources, increased use conflicts are predictable. Consequently, reasonable exploitation of subsurface resources requires an assessment of technologically achievable energy potentials with scientific based tools. We present application-oriented management tools which target on deriving shallow subsurface energy potentials. 3D groundwater flow and heat-transport models are used to capture groundwater flow and heat-transport dynamics on the city- and quarter-scale, 2D box models are used to quantify technically feasible extraction rates of well doublets for groundwater heat pump systems. For Basel (Switzerland), prospective large theoretical energy potentials can be derived for areas with high advective heat flux and high temperature gradients. Likewise, single city quarters are suitable for ‘active’ thermal use with well doublets, whereas thermal power potentials reach 1.2 MW. Regarding ‘passive’ installations of energy absorbers in subsurface structures located within the groundwater, energy potentials amount to 4 and up to 40 W m−2. The assessment results can be integrated into urban energy plans and support architects, city planners and potential users to acquire initial site-specific information on the technical feasibility of shallow geothermal energy systems.

The increase in the frequency of extreme precipitation events due to climate change, together with the continuous development of cities and surface sealing that hinder water infiltration into the subsoil, is accelerating the search for new facilities to manage stormwater. The Canary Islands (Spain) are taking advantage of the knowledge acquired in the construction of water mines to exploit a novel stormwater management facility, which we have defined as a dry gallery. Dry galleries are constituted by a vertical well connected to a horizontal gallery dug into highly permeable volcanic layers of the vadose zone, from where infiltration takes place. However, the lack of scientific knowledge about these facilities prevents them from being properly dimensioned and managed. In this work, we simulate for the first time the infiltration process and the wetting front propagation from dry galleries based on a 3D unsaturated flow model and provide some recommendations for the installation and sizing of these facilities. The fastest advance of the wetting front takes place during the earliest times of infiltration ( rates higher than 1000m∙d−1 and 2m3∙s−1. As time progresses, the propagation velocity and infiltration rate decrease as a consequence of the hydraulic gradient attenuation between the gallery and the aquifer. Therefore, stormwater infiltration is a highly transient process in which a sizing underestimation of 100%may be committed if unsaturated conditions or geological configuration are neglected.

Volcanic lakes in oceanic islands represent extremely important areas for biodiversity and offer exceptional conditions for nature-based tourism as one of the main pillars of economic growth in these regions. Nitrogen and phosphorus fertilizers are being used extensively at the Azores archipelago, similarly to other places in the world, to increase agricultural production and is causing severe pollution and eutrophication of surface freshwater reservoirs. This work concentrates on the evaluation of the efficiency of surface water diversion as a remediation measure to reduce nutrient loading and reverse eutrophication of Furnas crater lake on the island of São Miguel. Nutrient loading was monitored using an extensive water quality monitoring program in the main watershed of the eutrophicated lake that, together with watershed-scale mass balance methods and groundwater and solute transport models, allowed us to identify an average 98% efficiency in the reduction of nitrate loading. However, phosphorus total load in the discharged water to Furnas lake was observed to only be decreased by 33% due to the groundwater origin of phosphorus. Results from modelling suggested that nutrients were emitted from both point (nitrates) and diffuse (phosphorus) sources as surface runoff and as groundwater seepage, respectively. The results obtained recognized a partially successful surface water diversion, since groundwater path was not initially identified, thus highlighting the importance of groundwater flow regime in the design of such remediation measures. This work also provides a perspective on surface water diversion to revert eutrophication under a volcanic formation, where lakes can be naturally more nutrient rich.

Aviation safety is a priority that may be compromised by adverse weather conditions. This is the case for poor visibility in the vicinity of airports, which can pose a risk during takeoff and landing. For airports that are prone to fog because of their location, an accurate forecast of poor-visibility episodes is vital. However, the forecasting of low clouds is still a challenge in numerical weather prediction, especially when an airport is near complex terrain for which the use of non-hydrostatic mesoscale models is mandatory. All these factors are present at Tenerife Norte Airport, which is commonly affected by poor visibility from low clouds related to persistent trade winds and moist flows from the Atlantic Ocean. In this paper, several methods for estimating visibility based on mesoscale model outputs are tested. Use of the HARMONIE-AROME model is encouraged because of its excellent performance in the detection of poor-visibility episodes (False Alarm Ratio = 0.34–0.38; Frequency Of Misses = 0.22–0.38, depending on the model version and method used). In addition, the use of satellite application facilities is proposed for the nowcasting of low clouds affecting the airport area. Specifically, we used products that estimate cloud type, cloud top altitude, and integrated water vapor content in the boundary layer. Finally, an application is presented for the monitoring of weather conditions in real time to estimate poor-visibility risk.

© 2020 by the authors., The long-term sustainability of shallow geothermal systems in dense urbanized areas can be potentially compromised by the existence of thermal interfaces. Thermal interferences between systems have to be avoided to prevent the loss of system performance. Nevertheless, in this work we provide evidence of a positive feedback from thermal interferences in certain controlled situations. Two real groundwater heat pump systems were investigated using real exploitation data sets to estimate the thermal energy demand bias and, by extrapolation, to assess the nature of thermal interferences between the systems. To do that, thermal interferences were modelled by means of a calibrated and validated 3D city-scale numerical model reproducing groundwater flow and heat transport. Results obtained showed a 39% (522 MWh·yr−1) energy imbalance towards cooling for one of the systems, which generated a hot thermal plume towards the downgradient and second system investigated. The nested system in the hot thermal plume only used groundwater for heating, thus establishing a positive symbiotic relationship between them. Considering the energy balance of both systems together, a reduced 9% imbalance was found, hence ensuring the long-term sustainability and renewability of the shallow geothermal resource exploited. The nested geothermal systems described illustrate the possibilities of a new management strategy in shallow geothermal energy governance., Data analyses were undertaken within the EU H2020 GeoERA project MUSE–Managing Urban Shallow geothermal Energy. MUSE received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 731166 under the terms of GeoERA programme–ERA-NET Cofund Action

This document develops a methodology that evaluates the impact on the environment of the rivers produced by the creation of civil infrastructures. The methodology is based on the two-dimensional hydrodynamic calculation by using shallow water equations both in the conditions prior to the creation of the infrastructure, and in the new conditions after the infrastructure is created. Subsequently, several characteristics, such as water depth and velocity, among others, are compared between the initial and final conditions, and a two-dimensional zoning of the changes observed is obtained. The methodology herein presented is useful to verify the impact that the implantation of different infrastructures around the river currents could produce. In addition, it is also relevant for carrying out a study taking into account different infrastructure options related to river currents, as well as for selecting the most suitable one. By using the methodology presented, changes on the regime of the currents caused by the infrastructures can be deduced, including a qualitative and quantitative zoning of the changes, with a special emphasis on some characteristics, such as depth and velocity. The methodology is applied in a case study for the creation of a road bridge over the Jalon River in Spain.

Salt flats (salars) are endorheic hydrogeological systems associated with arid to hyperarid climates. The brines of salt flats account the 80 % of the world’s reserves of Li highly demanded by modern industry. About 40 % of the worldwide Li is extracted from the brine that fills the pores and cavities of the Salar de Atacama. However, the origin of the extreme Li-enrichment of these brines is still unknown.The thick accumulation of salts and brines in salt flats results from the groundwater discharge (phreatic evaporation) near the land surface for thousands to millions of years. The strong evaporation contributes the enrichment in major cations and anions as well as other rare elements (e.g. Li, B, Ba, Sr, Br, I and F) which are very attractive for mining exploitation. However, only evaporation cannot explain by itself the extreme concentrations of some of these elements and the strong decoupling between the most evaporated brines and the most Li-enriched brines in the Salar de Atacama. Several hypotheses have been proposed to explain the extreme Li-enrichment of the salt flat brines: (a) concentrated brines leaking down from salt flats located in the Andean Plateau, (b) leaching of hypothetical ancient salt flats buried among volcanic rocks, and (c) rising of hydrothermal brines from deep reservoirs through faults. However, none of them has been able probed neither validated by a numerical model till the date.The objective of this work is to discuss the feasibility of the different hypotheses proposed until now to explain the formation of the world's largest lithium reserve. To achieve this objective, two sets of numerical simulations of a 2D vertical cross-section of the entire Salar de Atacama basin are carried out to define (1) the origin and evolution of a salt flat and how climate cycles can affect the location of the most Li-concentrated brines by evaporation and (2) the establishment of the hydro-thermo-haline circulation of a mature salt flat basin.

Shallow geothermal systems are the most efficient and clean technology for the air-conditioning of buildings and constitutes an emergent renewable energy resource in the worldwide market. Undisturbed systems are capable of efficiently exchanging heat with the subsurface and transferring it to human infrastructures, providing the basis for the successful decarbonisation of heating and cooling demands of cities. Unmanaged intensive use of groundwater for thermal purposes as a shallow geothermal energy (SGE) resource in urban environments threatens the resources´ renewability and the systems´ performance, due to the thermal interferences created by a biased energy demand throughout the year. To ensure sustainability, scientifically-based criteria are required to prevent potential thermal interferences between geothermal systems. In this work, a management indicator (balanced sustainability index, BSI) applicable to groundwater heat pump systems is defined to assign a quantitative value of sustainability to each system, based on their intrinsic potential to produce thermal interference. The BSI indicator relies on the net heat balance transferred to the terrain throughout the year and the maximum seasonal thermal load associated. To define this indicator, 75 heating-cooling scenarios based in 23 real systems were established to cover all possible different operational conditions. The scenarios were simulated in a standard numerical model, adopted as a reference framework, and thermal impacts were evaluated. Two polynomial regression models were used for the interpolation of thermal impacts, thus allowing the direct calculation of the sustainability indicator developed as a function of heating-cooling ratios and maximum seasonal thermal loads. The BSI indicator could provide authorities and technicians with scientifically-based criteria to establish geothermal monitoring programs, which are critical to maintain the implementation rates and renewability of these systems in the cities.

The shallow subsurface comprising groundwater bodies as well as solid rock formations in the uppermost tens to hundreds of meters below surface offer significant resources for renewable heating, cooling and seasonal underground heat storage. Shallow geothermal energy (SGE) comprises the technologies to exchange heat between the subsurface and surface via closed loop or open loop heat exchangers. Although SGE just covered around 2% of the renewable heat production in the EU in 2018, its huge potential for low temperature heating and cooling supply is expected to lead to a significant market growth across Europe in the upcoming decade. Especially as SGE offers the unique possibility to supply heating, cooling and storing waste heat with one technology. SGE offers advantages especially in urban areas. It does not produce waste heat if applied for cooling, which can be considered as an important measure to mitigate urban heat islands. It consumes low amount of surface space for its installation and applying SGE is free of emissions in terms of aerosols or noise. Moreover, it can be combined with other renewables such as solar and waste heat or excess energy. In these cases, SGE acts as a seasonal heat storage.The increasing interest in SGE in urban areas, however, puts pressure on the subsurface, especially on shallow groundwater bodies. In that context, SGE systems may compete with each other as well as with water supply and subsurface installations. In many European countries, management approaches of SGE are either lacking or follow the first come first serve approach. Integrative management approaches follow an information and decision cycle, starting and ending at collecting, processing and providing geoscientific data on subsurface conditions to stakeholders, such as authorities, investors and city planners.GeoERA MUSE addresses integrative management approaches for the use of SGE by harmonizing concepts and testing them in 14 European cities facing different climatic, hydrogeologic and socio-economic boundary conditions. MUSE deals with mapping resources and limitations of SGE resources and displays them in modern web-based interfaces. Knowing resources and limitations referring to interference with other SGE systems or other shallow subsurface installations is the starting point for integrative management approaches, which include summation effects and abandon first come first serve. MUSE pilot areas follow the whole management cycle from creating subsurface data (e.g. subsurface temperatures, thermal rock properties), deriving resource models (amount of energy available for use), limitations of use (contaminated areas, problematic chemical composition of groundwater) and displaying the information gained at the EGDI web platform of EuroGeoSurveys. Furthermore, MUSE interacts with local stakeholders to transfer geoscientific data models into managing strategies. In that sense, the pilot areas act as role model for other urban regions in Europe. Additionally, MUSE creates joint concepts and standards to strengthen the role of Geological Survey Organisations towards a key player in managing an efficient and sustainable use of urban subsurface in general and SGE in urban areas in detail. MUSE has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 731166.

Different aspects of management policies for shallow geothermal systems are currently under development. Although this technology has been used for a long time, doubts and concerns have been raised in the last years due to the massive implementation of new systems. To assess possible environmental impacts and manage subsurface energy resources, collecting data from operating shallow geothermal systems is becoming mandatory in Europe. This study presents novel advances in the upscaling of operation datasets obtained from open-loop geothermal energy systems for an optimal integration in hydrogeological models. The proposed procedure allows efficient numerical simulations to be performed at an urban scale. Specifically, this work proposes a novel methodology to optimize the data treatment of highly transient real exploitation regimes by integrating energy transfer in the environment to reduce more than 90% registered raw datasets. The proposed methodology is then applied to and validated on five different real optimization scenarios in which upscaling transformation of the injection temperature series of 15-min sampling frequency has been considered. The error derived from each approach was evaluated and compared for validation purposes. The results obtained from the upscaling procedures have proven the usefulness and transferability of the proposed method for achieving daily time functions to efficiently reproduce the exploitation regimes of these systems with an acceptable error in a sustainable resource management framework.

The shallow subsurface in urban areas is increasingly used by shallow geothermal energy systems as a renewable energy resource and as a cheap cooling medium, e.g. for building air conditioning. In combination with further anthropogenic activities, this results in altered thermal regimes in the subsurface and the so-called subsurface urban heat island effect. Successful thermal management of urban groundwater resources requires understanding the relative contributions of the different thermal parameters and boundary conditions that result in the “present thermal state” of individual urban groundwater bodies. To evaluate the “present thermal state” of urban groundwater bodies, good quality data are required to characterize the hydraulic and thermal aquifer parameters. This process also involved adequate monitoring systems which provide consistent subsurface temperature measurements and are the basis for parameterizing numerical heat-transport models. This study is based on previous work already published for two urban groundwater bodies in Basel (CH) and Zaragoza (ES), where comprehensive monitoring networks (hydraulics and temperature) as well as calibrated high-resolution numerical flow- and heat-transport models have been analyzed. The “present thermal state” and how it developed according to the different hydraulic and thermal boundary conditions is compared to a “potential natural state” in order to assess the anthropogenic thermal changes that have already occurred in the urban groundwater bodies we investigated. This comparison allows us to describe the various processes concerning groundwater flow and thermal regimes for the different urban settings. Furthermore, the results facilitate defining goals for specific aquifer regions, including future aquifer use and urbanization, as well as evaluating the thermal use potential for these regions. As one example for a more sustainable thermal use of subsurface water resources, we introduce the thermal management concept of the “relaxation factor”, which is a first approach to overcome the present policy of “first come, first served”. Remediation measures to regenerate overheated urban aquifers are also introduced. The transferability of the applied methods to other urban areas is discussed. It is shown that an appropriate selection of locations for monitoring hydraulic and thermal boundary conditions make it possible to implement representative interpretations of groundwater flow and thermal regimes as well as to set up high-resolution numerical flow- and heat-transport models. Those models are the basis for the sustainable management of thermal resources.

As a result of the increasing use of shallow geothermal resources, hydraulic, thermal and chemical impacts affecting groundwater quality can be observed with ever increasing frequency (Possemiers et al., 2014). To overcome the uncertainty associated with chemical impacts, a city scale study on the effects of intensive geothermal resource use by groundwater heat pump systems on groundwater quality, with special emphasis on heavy metal contents was performed. Statistical analysis of geochemical data obtained from several field campaigns has allowed studying the spatiotemporal relationship between temperature anomalies in the aquifer and trace element composition of groundwater. The relationship between temperature and the concentrations of trace elements resulted in weak correlations, indicating that temperature changes are not the driving factor in enhancing heavy metal contaminations. Regression models established for these correlations showed a very low reactivity or response of heavy metal contents to temperature changes. The change rates of heavy metal contents with respect to temperature changes obtained indicate a low risk of exceeding quality threshold values by means of the exploitation regimes used, neither producing nor enhancing contamination significantly. However, modification of pH, redox potential, electrical conductivity, dissolved oxygen and alkalinity correlated with the concentrations of heavy metals. In this case, the change rates of heavy metal contents are higher, with a greater risk of exceeding threshold values.

A steady increase in the consumption of pharmaceuticals and personal-care products worldwide is increasing their occurrence in the biosphere. The current study describes the abundance of 42 selected emerging organic contaminants (EOCs), including human and veterinary antibiotics, UV-filters and analgesics in the groundwater of the urban aquifer of Zaragoza (Spain), which is affected by intensive exploitation of shallow geothermal resources. The presence of groundwater heat pump systems in the aquifer studied offered the opportunity to study the occurrence of EOCs in relation to groundwater temperature and other physicochemical effects derived from this technology. Analysis of the data obtained allowed us to identify statistically significant relationships between the presence of EOCs and temperature, as well as other physicochemical and geochemical properties of groundwater. The results obtained suggest that temperature is a minor factor controlling the degradation of the organic compounds analysed compared to the oxygen input from groundwater heat pump systems which is possibly increasing the aerobic redox conditions, thus preventing the degradation of organic pollutants. Intensive use of shallow geothermal resources therefore seems to contribute in the prevalence of such compounds in the aquifer close to geothermal systems.

An actual event that happened in the Roncal valley (Spain) is investigated and the results are compared between models with and without snowmelt. A distributed rainfall model is generated with the specific data recorded by the rain gauges of the catchment during the episode. To describe the process of water routing in the hydrological cycle of the basin, a model is used based on combinations of parallel linear reservoirs (PLR model), distribution by the basin, and tip-out into its drainage network configured using a digital terrain model (DTM). This PLR model allows simulation of the different actual reservoirs of the basin, including the snow and the contribution due to its melting which, in the model, depends on the temperature. The PLR model also allows for a water budget of the episode where, in addition to the effective rainfall contribution, the water that comes from the thaw is taken into account. The PLR model also allows determination of the amount of water that exists in the basin before and after the episode, data of great interest. When comparing the simulations with and without taking into account the thawing process, it is evident that the intervention of the snow reservoir has been decisive in causing a flood to occur.

The steady increase of geothermal systems using groundwater is compromising the renewability of the geothermal resources in shallow urban aquifers. To ensure sustainability, scientifically-based criteria are required to prevent potential thermal interferences between geothermal systems. In this work, a management indicator (balanced sustainability index, BSI) applicable to groundwater heat pump systems is defined to assign a quantitative value of sustainability to each system, based on their intrinsic potential to produce thermal interference. The BSI indicator relies on the net heat balance transferred to the terrain throughout the year and the maximum seasonal thermal load associated. To define this indicator, 75 heating-cooling scenarios based in 23 real systems were established to cover all possible different operational conditions. The scenarios were simulated in a standard numerical model, adopted as a reference framework, and thermal impacts were evaluated. Two polynomial regression models were used for the interpolation of thermal impacts, thus allowing the direct calculation of the sustainability indicator developed as a function of heating-cooling ratios and maximum seasonal thermal loads. The BSI indicator could provide authorities and technicians with scientifically-based criteria to establish geothermal monitoring programs, which are critical to maintain the implementation rates and renewability of these systems in the cities. © 2018 Elsevier Ltd, This work was performed by the Geological Survey of Spain (IGME) under the Cooperation Agreement framework between IGME and the Ebro Hydrographic Confederation (CHE) for the specific arrangement entitled “The application of a groundwater flow and heat transport numerical model for the simulation of management strategies of geothermal installations in the city of Zaragoza”. Alejandro García-Gil gratefully acknowledges MIKE by DHI for the sponsored FEFLOW license.