Using Stable Isotopes to Assess Groundwater Recharge and Solute Transport in a Density-Driven Flow-Dominated Lake-Aquifer System

Saline lakes are mostly located in endorheic basins in arid and semi-arid regions, where the excess of evaporation over precipitation promotes the accumulation of salts on the surface. As the salinity of these lakes increases, their mass balance changes, and biogeochemical processes may be intensified. In that sense, Petrola Lake (SE Spain) is a terminal lake located in an endorheic basin with elevated anthropic pressure, mainly derived from agricultural inputs and wastewater discharge. The goal of this study was to evaluate the interaction between groundwater and saline water from Petrola Lake to improve our knowledge of groundwater recharge processes by density-driven flow (DDF) in terminal lakes. A combination of hydrochemical (chloride concentration) and stable isotope (delta¹⁸OH2O and delta²HH2O) data were used. In order to test the conceptual model, a simple numerical experiment was performed using a one-dimensional column that represents the relationship between the lake and the aquifer incorporating the variable density coupling control in solute migration. The isotopic composition of 190 groundwater and surface water samples collected between September 2008 and July 2015 provides a regression line (delta²HH2O = 5.0 center dot delta¹⁸OH2O - 14.3 parts per thousand, R² = 0.95) consistent with dominant evaporation processes in the lake. The DDF towards the underlying aquifer showed a strong influence on the mixing processes between the groundwater and surface water. Nevertheless, groundwater chemistry at different depths beneath the lake remains almost constant over time, suggesting an equilibrium between DDF and regional groundwater flow (RGF). Modelling isotope changes allowed inferring the temporal pattern of saline water recharge, coinciding with the summer season when water loss through evaporation is most significant. Consequently, the transport of solutes suitable for chemical reactions is then feasible to deeper zones of the aquifer.