Evaluation of the impact of the intensive exploitation of groundwater and the mega-drought based on the hydrochemical and isotopic composition of the waters of the Chacabuco-Polpaico basin in central Chile.

A hydrogeochemical and isotopic study has been carried out to understand the hydrogeological functioning of a small alluvial aquifer in central Chile in a context of mega-drought and intensive exploitation of its waters. Additionally, two mine tailings dams from porphyry copper mining are situated in the area. The prolonged mega-drought, which has lasted for over thirteen years, has resulted in a significant decrease in rainfall recharge and a drop of up to 50 m in piezometric levels, although no serious groundwater contamination problems have yet been detected, except for a rise in nitrate contents (ranging between 23 and 45 mg/L NO3) attributed to return irrigation. Groundwaters are calcium-bicarbonate and calcium-sodium-bicarbonate in composition. The values of δ18O and δ2H of the alluvial aquifer indicate fractionation by evaporation that would be explained by the recirculation of water that occurs in the agricultural areas of the basin, where the excess irrigation water that go back to the aquifer presents fractionation by evaporation. The δ34S and δ18O of dissolved sulfate point to pyrite oxidation, which could be related to the pyrite present in the copper porphyry and recognized in the Andes Cordillera. The 87Sr/86Sr isotopic values of the alluvial aquifer waters are close to the isotopic fingerprint of the volcanic rocks of the Abanico Formation. However, the water from the wells located further downstream in the basin and close to the tailing dams show δ34S and δ18O of dissolved sulfate and 87Sr/86Sr consistent with Miocene intrusive mineralogies of the copper porphyry type. The groundwater chemistry does not show water seepage from the tailings dam. Therefore, a minor contribution of minerals related to the intrusive rocks is proposed, which would originate from the movement of fine particles by the wind from the dams to the valley floor. The 14C activities indicate that groundwater is recent.

Processes explaining the origin and evolution of groundwater composition in the Andean Precordillera and Altiplano of the Tarapacá Region of northern Chile.

In the arid area of northern Chile, groundwater resources in the Andean formations are essential for native populations, ecological services, mining, and other human activities. Validated conceptual hydrogeological models are required for current and future water and land management. This work aims to explain the processes controlling the origin and distribution of recharge and groundwater composition in the Andean Precordillera and Altiplano of the Tarapacá Region of northern Chile, using major solutes in spring, river, and well water, and the stable and radioactive isotopes of water oxygen, hydrogen, and dissolved inorganic carbon. The waters are mainly of the Na-Ca-SO4 type. Processes controlling the chemical evolution of waters are atmospheric dust contribution, evapo-concentration, and enhanced volcanic rock weathering, as well as halite dissolution in some locations. The isotopic composition of Precordillera eastern flank water samples follows an evaporation line, while those in the western flank, in the Altiplano, follow a line that is parallel to the local meteoric line, suggesting unsaturated zone evaporation processes of infiltrated rainfall. δ13CDIC contents (-2 to -27‰) indicate mixing processes, volcanic CO2 in the Altiplano, and calcite dissolution in some sectors. In the western depression, the only recharge is due to water infiltration in creek channels. In the highland areas, 5-25% of precipitation produces recharge. The estimated groundwater renewal time in the Precordillera was 3-14 kyr. The piezometric elevation in the Precordillera due to low-permeability intrusive rocks and local recharge prevents the east-west groundwater transfer from the Altiplano to the western depression and explains why the volcanic CO2 in the Altiplano basins is not observed on the western flank. These results provide new insights for the evolution of water quality in volcanic aquifers in arid environments and provide considerations for estimating groundwater residence times using radiocarbon in areas influenced by volcanic CO2.

Recharge and residence times of groundwater in hyper arid areas: The confined aquifer of Calama, Loa River Basin, Atacama Desert, Chile.

Groundwater recharge in hyper arid areas often depends on surface water infiltration and diffuse recharge of highly evaporated precipitation only contribute under favorable conditions. This happens in the Calama basin two-aquifer system, in the Central Andean area of northern Chile. A conceptual model of the groundwater system and its relationship with the Loa River is defined. We focus on the confined aquifer of the Calama basin, combining hydrodynamic, hydrogeochemical and isotopic methods. Radiocarbon (14C) activity data of dissolved inorganic carbon (DIC), in conjunction with chemical data, are applied to evaluate groundwater residence time within the confined aquifer. The Loa River recharges the Calama basin aquifers in its northeastern part, with water that has chemical and isotopic characteristics inherited from the arid environment and volcanic rocks in its upper basin. In the central and northeastern part of the confined aquifer, minor variations in chloride concentration suggest that the deep aquifer is well confined. The δ18O and δ2H values in groundwater of the confined aquifer show an increasing isotopic fractionation from the recharge area (around -10‰ δ18O) to those in the discharge area (between -8.5‰ and -8‰) in the southwestern part of the aquifer. The 14C activity continuously decreases down flow from the recharge by the Loa River. Adjusted DIC radiocarbon ages indicate a groundwater travel time between 1500 and 4000 years in the confined aquifer of Calama. Despite the limitations and uncertainties of radiocarbon in DIC to estimate groundwater transit times for the confined aquifer and considering complementary chemical and isotopic constraints, the DIC 14C provides acceptable values. The approach may be applicable in other confined aquifers in hyper-arid climates in which the formation of aquifer systems linked to river damming by geological action took place. This information is needed for sound management of the scarce groundwater resources.

Groundwater recharge and hydrodynamics of complex volcanic aquifers with a shallow saline lake: Laguna Tuyajto, Andean Cordillera of northern Chile.

Estimating groundwater recharge in arid regions with seasonal snow cover, as in the Andean Altiplano of northern Chile, is important for sustainable development policies and the effective management of scarce resources in a high water demanding zone, as fragile ecosystems depends on a stable water contribution. This research aims to evaluate and quantify rainfall and snowfall contribution to aquifer recharge while assessing the factors that control the hydrodynamics in such areas, based in the knowledge of the better documented Tuyajto Lake in the Tuyajto catchment/basin. The modeling framework involves an energy balance of the snow cover, a soil water balance and a groundwater flow and chloride transport model. The basin average annual recharge is about 23% of average precipitation. Snowmelt contribution to recharge is important at altitudes above 4700 m a.s.l. during September, while rainfall is more important in February and March, during short intense precipitation events. The hydraulic conductivity of ignimbrites and other volcanic formations are the most important hydrogeological parameters controlling lake level and spring flow rates, while albedo and snowpack surface roughness length on the energy balance causes the greatest variation of lake level. Evaporation is the process controlling the variability of the lake level, as aquifer contribution remains relatively constant and springs flow variability is not enough to cause the observed variations, except during November. Possible buried salts deposits on the eastern edges of Pampa Colorada and Tuyajto Lake, together with volcanic HCl contribution, justify the high measured groundwater chloride concentrations. A recharge 2-3 higher than the current one is necessary to justify a lake level 40 m above its modern value during the Last Glacial period, giving insight into past hydrological changes in the basin due to climate variability. The knowledge gained can be applied to other high altitude volcanic basins with seasonal snow cover.

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.

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.

Groundwater discharge in high-mountain watersheds: A valuable resource for downstream semi-arid zones. The case of the Bérchules River in Sierra Nevada (Southern Spain).

Aquifers in permeable formations developed in high-mountain watersheds slow down the transfer of snowmelt to rivers, modifying rivers' flow pattern. To gain insight into the processes that control the hydrologic response of such systems the role played by groundwater in an alpine basin located at the southeastern part of the Iberian Peninsula is investigated. As data in these environments is generally scarce and its variability is high, simple lumped parameter hydrological models that consider the groundwater component and snow accumulation and melting are needed. Instead of using existing models that use many parameters, the Témez lumped hydrological model of common use in Spain and Ibero-American countries is selected and modified to consider snow to get a simplified tool to separate hydrograph components. The result is the TDD model (Témez-Degree Day) which is applied in a high mountain watershed with seasonal snow cover in Southern Spain to help in quantifying groundwater recharge and determining the groundwater contribution to the outflow. Average groundwater recharge is about 23% of the precipitation, and groundwater contribution to total outflow ranges between 70 and 97%. Direct surface runoff is 1% of precipitation. These values depend on the existence of snow. Results are consistent with those obtained with chloride atmospheric deposition mass balances by other authors. They highlight the important role of groundwater in high mountain areas, which is enhanced by seasonal snow cover. Results compare well with other areas. This effect is often neglected in water planning, but can be easily taken into account just by extending the water balance tool in use, or any other, following the procedure that has being developed.

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.

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.