Jorge Jódar

Isotopic and hydrogeochemical characterization of high-altitude karst aquifers in complex geological settings. The Ordesa and Monte Perdido National Park (Northern Spain) case study

The Ordesa and Monte Perdido National Park, located in the Southern Pyrenees, constitutes the highest karst system in Western Europe. No previous studies regarding its geochemical and isotopic groundwater characterization are available in this area. This work presents the results of field and sampling campaigns carried out between July 2007 and September 2013. The groundwater presents high calcium bicarbonate contents due to the occurrence of upper Cretaceous and lower Paleocene-Eocene carbonate materials in the studied area. Other relevant processes include dissolution of anhydrite and/or gypsum and incongruent dissolution of Mg-limestone and dolomite. The water stable isotopes (δ(18)O, δ(2)H) show that the oceanic fronts from the Atlantic Ocean are responsible for the high levels of precipitation. In autumn, winter, and spring, a deuterium excess is found in the recharge water, which could be related to local atmospheric transport of low-altitude snow sublimation vapour and its later condensation on the snow surface at higher altitude, where recharge is mostly produced. The recharge zones are mainly between 2500m and 3200ma.s.l. The tritium content of the water suggests short groundwater transit times. The isotopic composition of dissolved sulphate points to the existence of regional fluxes mixed with local discharge in some of the springs. This work highlights the major role played by the altitude difference between the recharge and discharge zones in controlling the chemistry and the vertical variability of the isotopic composition in high-altitude karst aquifers.

Groundwater flow in a closed basin with a saline shallow lake in a volcanic area: Laguna Tuyajto, northern Chilean Altiplano of the Andes

Laguna Tuyajto is a small, shallow saline water lake in the Andean Altiplano of northern Chile. In the eastern side it is fed by springs that discharge groundwater of the nearby volcanic aquifers. The area is arid: rainfall does not exceed 200mm/year in the rainiest parts. The stable isotopic content of spring water shows that the recharge is originated mainly from winter rain, snow melt, and to a lesser extent from some short and intense sporadic rainfall events. Most of the spring water outflowing in the northern side of Laguna Tuyajto is recharged in the Tuyajto volcano. Most of the spring water in the eastern side and groundwater are recharged at higher elevations, in the rims of the nearby endorheic basins of Pampa Colorada and Pampa Las Tecas to the East. The presence of tritium in some deep wells in Pampa Colorada and Pampa Las Tecas indicates recent recharge. Gas emission in recent volcanoes increase the sulfate content of atmospheric deposition and this is reflected in local groundwater. The chemical composition and concentration of spring waters are the result of meteoric water evapo-concentration, water-rock interaction, and mainly the dissolution of old and buried evaporitic deposits. Groundwater flow is mostly shallow due to a low permeability ignimbrite layer of regional extent, which also hinders brine spreading below and around the lake. High deep temperatures near the recent Tuyajto volcano explain the high dissolved silica contents and the δ(18)O shift to heavier values found in some of the spring waters. Laguna Tuyajto is a terminal lake where salts cumulate, mostly halite, but some brine transfer to the Salar de Aguas Calientes-3 cannot be excluded. The hydrogeological behavior of Laguna Tuyajto constitutes a model to understand the functioning of many other similar basins in other areas in the Andean Altiplano.

Groundwater origin and recharge in the hyperarid Cordillera de la Costa, Atacama Desert, northern Chile

The Cordillera de la Costa is located along the coastline of northern Chile, in the hyperarid Atacama Desert area. Chemical and isotopic analyses of several small coastal springs and groundwater reservoirs between 22.5 °S and 25.5 °S allow understanding groundwater origin, renewal time and the probable timing of recharge. The aquifers are mostly in old volcanic rocks and alluvial deposits. All spring waters are brackish, of the sodium chloride type due to intensive concentration of precipitation due aridity and for deep groundwater to additional water-rock interaction in slowly renewed groundwater and mixing with deep seated brines. The heavy δ18O and δ2H values in spring water are explained by recharge by the arrival of moist air masses from the Pacific Ocean and the originally lighter values in the deep wells can be associated to past recharge by air masses coming from the Atlantic Ocean. Current recharge is assumed almost nil but it was significant in past wetter-than-present periods, increasing groundwater reserves, which are not yet exhausted. To explain the observed chloride content and radiocarbon (14C) activity, a well-mixed (exponential) flow model has been considered for aquifer recharge. The average residence time of groundwater feeding the springs has been estimated between 1 and 2kyr, up to 5kyr and between 7 and 13kyr for deep well water, assuming that current recharge is much less than during the previous wetter period. The recharge period feeding the coastal springs could have been produced 1 to 5kyr BP, when the area was already inhabited, and recharge in the Michilla mine was produced during the 10 to 14.5kyr BP CAPE (Central Andean Pluvial Event) pluvial events of the central Andes. The approximate coincidence of turnover time with the past wet periods, as revealed by paleoclimate data, points to significant recharge during them.

Land surface temperature as an indicator of the unsaturated zone thickness: A remote sensing approach in the Atacama Desert

Land surface temperature (LST) seems to be related to the temperature of shallow aquifers and the unsaturated zone thickness (∆Zuz). That relationship is valid when the study area fulfils certain characteristics: a) there should be no downward moisture fluxes in an unsaturated zone, b) the soil composition in terms of both, the different horizon materials and their corresponding thermal and hydraulic properties, must be as homogeneous and isotropic as possible, c) flat and regular topography, and d) steady state groundwater temperature with a spatially homogeneous temperature distribution. A night time Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image and temperature field measurements are used to test the validity of the relationship between LST and ∆Zuz at the Pampa del Tamarugal, which is located in the Atacama Desert (Chile) and meets the above required conditions. The results indicate that there is a relation between the land surface temperature and the unsaturated zone thickness in the study area. Moreover, the field measurements of soil temperature indicate that shallow aquifers dampen both the daily and the seasonal amplitude of the temperature oscillation generated by the local climate conditions. Despite empirically observing the relationship between the LST and ∆Zuz in the study zone, such a relationship cannot be applied to directly estimate ∆Zuz using temperatures from nighttime thermal satellite images. To this end, it is necessary to consider the soil thermal properties, the soil surface roughness and the unseen water and moisture fluxes (e.g., capillarity and evaporation) that typically occur in the subsurface.

Vertical variation in the amplitude of the seasonal isotopic content of rainfall as a tool to jointly estimate the groundwater recharge zone and transit times in the Ordesa and Monte Perdido National Park aquifer system, north-eastern Spain

The time series of stable water isotope composition relative to meteorological stations and springs located in the high mountainous zone of the Ordesa and Monte Perdido National Park are analyzed in order to study how the seasonal isotopic content of precipitation propagates through the hydrogeological system in terms of the aquifer recharge zone elevation and transit time. The amplitude of the seasonal isotopic composition of precipitation and the mean isotopic content in rainfall vary along a vertical transect, with altitudinal slopes for δ18O of 0.9‰/km for seasonal amplitude and -2.2‰/km for isotopic content. The main recharge zone elevation for the sampled springs is between 1950 and 2600m·a.s.l. The water transit time for the sampled springs ranges from 1.1 to 4.5yr, with an average value of 1.85yr and a standard deviation of 0.8yr. The hydrological system tends to behave as a mixing reservoir.

Correlation of the seasonal isotopic amplitude of precipitation with annual evaporation and altitude in alpine regions.

The time series of stable water isotope composition relative to IAEA-GNIP meteorological stations located in alpine zones are analyzed in order to study how the amplitude of the seasonal isotopic composition of precipitation (Aδ) varies along a vertical transect. A clear relationship between Aδ and local evaporation is obtained, with slopes of -0.87 ‰/100mm/yr and -7.3 ‰/100mm/yr for Aδ(18)O and Aδ(2)H, respectively. When all sampling points of the vertical transect receive the same moisture sources, then a linear relationship between Aδ and elevation is obtained, with vertical gradients of 0.16 ‰/100mm/yr and 1.46 ‰/100mm/yr forAδ(18)O and Aδ(2)H, respectively.

Groundwater dynamics in a hydrologically-modified alpine watershed from an ancient managed recharge system (Sierra Nevada National Park, Southern Spain): Insights from hydrogeochemical and isotopic information

In many of the alpine watersheds of Sierra Nevada (Southern Spain) exists an ancient network of dug canals that collect, transport and facilitate the recharge the snowmelt in the underlying aquifer during the spring season. This practice, known as careos, in the lower part of the watersheds supply drinking water as spring discharge during the dry season. To study how this managed recharge technique modifies the natural response of these basins this work focuses on characterizing the hydrological behavior of one of the sites, the Berchules watershed. The mechanisms for mineralization of groundwater are based on geochemical processes such as evapo-concentration in the soil layer and silicate mineral weathering due to dissolved CO2 originated from both soil biogenic processes and the atmosphere. Groundwater presents a main hydrogeochemical calcium‑magnesium-bicarbonate type facies, which is associated to groundwater flowing through the upper weathered silicates and quickly drained through springs located in the uplands and in the intermediate altitude catchment zone. Additionally, in the lower part of the basin some springs discharge mineralized groundwater with a sodium-calcium-bicarbonate composition associated to regional groundwater flow. In natural conditions, this hydrogeological system behaves as a sloping aquifer, occurring recharge between 1400 and 2500 m a.s.l. The springs discharge groundwater with an isotopic content and temperature in coherence with the local rainfall isotopic and thermal atmospheric altitudinal lines. Nevertheless, once the careo recharge begins the affected springs reveal the fingerprint of the concentrated recharge system by blurring the fingerprint of both the isotopic and thermal altitudinal dependence in the springs discharge. This validates the previous conceptual model and supports average recharge values of 141 ± 140 mm/yr and total average water resources of 181 ± 111 mm/yr which include a 40% increase in the study period due to the effect of the acequias de careo.

Groundwater Isotopic Characterization in Ordesa and Monte Perdido National Park (Northern Spain)

The Ordesa and Monte Perdido National Park constitutes the largest calcareous mountain range of Western Europe, where the highest altitude karst of Europe is found. No previous studies regarding groundwater isotopic characterization in this area are known. This work presents the results of two preliminary campaigns carried out during July 2007 and April 2012. The water stable isotopes (δ18O, δ2H) show that the oceanic fronts from the Atlantic are responsible for the high levels of precipitation. In autumn, winter, and spring time, a deuterium excess is found in recharge water, which could be related to snow sublimation and its later condensation on the snow surface. The recharge zones are between 2,500 m and 3,200 m asl. The water tritium content points to short groundwater transit times.

Combining point and regular lattice data in geostatistical interpolation

Abstract This work studies how to include both point and areal measurements when estimating gaussian fields by kriging. To achieve this objective, three geostatistical approaches are considered for the areal distributed data: (a) regionalized measurements that are geographically referenced by their centroid as if they were point measurements, (b) regionalized measurements that are explicitly accounted by formally computing all the needed covariances (i.e. area-to-area, area-to-point and point-to-point covariances, respectively) and (c) regionalized measurements that are used as an external drift variable. Results indicate that the measurement error corresponding to the areal data plays a key role to decide when the spatial support of the areal measurements is relevant. For small measurement errors, it is necessary to explicitly consider the spatial support of the areal measurements to avoid large estimation variances. For large measurement errors, the difference between defining areal measurements by using their actual spatial support and defining areal measurements by referencing them by their centroids (i.e. gravity centre) is small. In this situation, it is possible to use the areal measurements as an external drift instead of merging both types of information (i.e. point and areal data) as measurements for kriging. In this case, the cross-validation analysis shows a larger coefficient of determination, similar average kriging variance and smaller mean square error than the obtained in the case of merging point and areal measurements for kriging.