Geochemistry of the Bagnères-de-Bigorre thermal waters from the North Pyrenean Zone (France)

Thermal springs are poorly known in the sedimentary sites of the Pyrenees. In this paper, we describe the ‘Bagnères‐de‐Bigorre’ springs which occur in a remarkably active seismotectonic context. A chemical and isotopic study of 15 spring waters (both cold and thermal, ranging in temperature from 7.0 to 49.9°C), and continuous monitoring of a single spring allow us to characterise water–rock interactions, fluids paths and mixing processes. Three groups of waters are distinguished: (I) SO₄²⁺–Ca²⁺–Cl^– thermal waters (II) SO₄²⁺–Cl^––Ca²⁺ thermal waters and (III) HCO₃^––Ca²⁺ cold shallow waters. Their characteristics suggest interactions with Mesozoic evaporite and carbonate formations. O and D isotopes from thermal waters indicate a local meteoric origin of Atlantic signature and a recharge elevation of 800 to 1000 m, which corresponds to a single feeding area. Their δ¹³C values (?2.8 to ?9.6‰) are consistent with carbonate dissolution, slight fractionation and a surficial organic input leading to δ¹³C depletion. Sr isotopes (0.70751 to 0.70777), Na⁺/Cl^– and (Ca²⁺ + Mg²⁺)/SO₄^2– ratios as well as thermodynamic calculations show that the dissolution of anhydrite and halite‐bearing Triassic layers control the chemical composition of group‐I and ‐II waters. The contrasting trends of cation/Cl^– ratios and TDS of waters from groups I and II suggest the existence of two different circulation paths at depth as well as dilution with surficial waters similar to group III. Calculated mixing proportions show that three waters from group I are diluted from 17 to 66%, whereas all waters from group II are mixed. The aquifer temperature is estimated to be in the range 55–64°C using the retrograde and prograde solubilities of anhydrite and chalcedony, respectively. Accordingly, the mean depth of the reservoir is around 1.7 km, which allows us to constrain the depth of the Triassic layer.