Abstract: The hydrogeological unit of Aguadulce (Campo de Dalías aquifers, SE Spain) has a complex geometry. This fact, together with a continuous rise in water demand due to intensive agriculture and tourism create problems for groundwater quantity and quality. In this paper classic geochemical tools managed by means of GIS software and geochemical simulations are combined to delineate, identify and locate the possible physicochemical processes acting in the Aguadulce groundwater. Two main aquifers can be distinguished: the carbonate or lower aquifer of Triassic age, and the calcodetritic or upper aquifer of Plio-Quaternary age. Groundwaters from the latter are more saline and, assuming all chlorinity originates from seawater intrusion, the seawater contribution to their composition would be up to 7%. Nevertheless the carbonate aquifer appears not to be homogeneous: it is compartmentalised into 4 zones where different processes explain the different groundwaters compositions. Zone 4 samples (E margin of the carbonate aquifer) resemble those of the Plio-Quaternary aquifer, where calcite precipitation, dolomite and gypsum dissolution and some cation exchange (water–rock interaction) together with seawater–freshwater mixing occur. In contrast, water–rock interaction predominates in zones 1 and 3 of the carbonate aquifer. Moreover, zone 2 samples, located between zones 1 and 3, are explained by water–rock interaction in addition to mixing with Plio-Quaternary aquifer waters. The combination of geochemical simulations with GIS and hydrogeochemical analyses has proven to be effective in identifying and locating the different physicochemical processes in the aquifer areas, thus improving understanding of hydrogeochemistry in complex aquifers.

Abstract: In this study a typical coastal karst aquifer, developed in lower Cretaceous limestones, on the western Mediterranean seashore (La Clape massif, southern France) was investigated. A combination of geochemical and isotopic approaches was used to investigate the origin of salinity in the aquifer. Water samples were collected between 2009 and 2011. Three groundwater groups (A, B and C) were identified based on the hydrogeological setting and on the Cl− concentrations. Average and maximum Cl− concentrations in the recharge waters were calculated (ClRef. and ClRef.Max) to be 0.51 and 2.85mmol/L, respectively). Group A includes spring waters with Cl− concentrations that are within the same order of magnitude as the ClRef concentration. Group B includes groundwater with Cl− concentrations that range between the ClRef and ClRef.Max concentrations. Group C includes brackish groundwater with Cl− concentrations that are significantly greater than the ClRef.Max concentration. Overall, the chemistry of the La Clape groundwater evolves from dominantly Ca–HCO3 to NaCl type. On binary diagrams of the major ions vs. Cl, most of the La Clape waters plot along mixing lines. The mixing end-members include spring waters and a saline component (current seawater or fossil saline water). Based on the Br/Clmolar ratio, the hypothesis of halite dissolution from Triassic evaporites is rejected to explain the origin of salinity in the brackish groundwater. Groundwaters display 87Sr/86Sr ratios intermediate between those of the limestone aquifer matrix and current Mediterranean seawater. On a Sr mixing diagram, most of the La Clape waters plot on a mixing line. The end-members include the La Clape spring waters and saline waters, which are similar to the deep geothermal waters that were identified at the nearby Balaruc site. The 36Cl/Cl ratios of a few groundwater samples from group C are in agreement with the mixing hypothesis of local recharge water with deep saline water at secular equilibrium within a carbonate matrix. Finally, PHREEQC modelling was run based on calcite dissolution in an open system prior to mixing with the Balaruc type saline waters. Modelled data are consistent with the observed data that were obtained from the group C groundwater. Based on several tracers (i.e. concentrations and isotopic compositions of Cl and Sr), calculated ratios of deep saline water in the mixture are coherent and range from 3% to 16% and 0% to 3% for groundwater of groups C and B, respectively. With regard to the La Clape karst aquifer, the extension of a lithospheric fault in the study area may favour the rise of deep saline water. Such rises occur at the nearby geothermal Balaruc site along another lithospheric fault. At the regional scale, several coastal karst aquifers are located along the Gulf of Lion and occur in Mezosoic limestones of similar ages. The 87Sr/86Sr ratios of these aquifers tend toward values of 0.708557, which suggests a general mixing process of shallow karst waters with deep saline fossil waters. The occurrence of these fossil saline waters may be related to the introduction of seawater during and after the Flandrian transgression, when the highly karstified massifs invaded by seawater, formed islands and peninsulas along the Mediterranean coast.