Long-term evolution of a carbonate reservoir submitted to fresh, saline and thermal waters interactions – Jurassic carbonates in the coastal area of the Gulf of Lion margin (southern France)

– Securing and managing underground water resources requires a good knowledge of the structure, texture and connections of the reservoir, in order to develop realistic and reliable hydrogeological models. On the coastline of the Gulf of Lion Margin (S. France), the Balaruc-les-Bains deep karst reservoir is subjected to interactions between fresh, marine and deep thermal waters, respectively. Water resource usage for drinking, spa resort, and ﬁ sh-farming raises important economic and social issues. These were addressed by an integrated research program, involving drilling of an exploratory borehole across the Jurassic carbonate reservoir. This contribution analyses the 750 m cores, in order to (i) characterise the architecture and evolution of the karst reservoir and (ii) investigate the paleo-ﬂ uids circulations, witnessed by calcite and dolomite mineralization in the fractures, karst cavities, and as cement of tectonic beccia. The structure of the reservoir is characterised by the superposition of several aquifers separated by marly intervals. At shallow level, the initial grainstone is incompletely dolomitized in metre-thick intervals, while limestone in the 210 – 340 m interval was completely dolomitized at an early stage. Dolomite has been subjected to penetrative extensional cataclastic deformation, while the preserved limestone is affected by normal faulting, resulting from NNE – SSW extension. Distinct types of karsts have been documented, from the top of the reservoir (paleo-lapiaz ﬁ lled with Burdigalian marine marls), down to 500 m depth (paleo-endokarst ﬁ lled with continental silts). The upper reservoir (75 – 150 m) is intensely karsti ﬁ ed, and includes 0.1 to 1 m-wide cavities, where present day water ﬂ uxes are documented. Analyses of calcite and dolomite crystallisation under natural light and cathodoluminescence indicate precipitation from distinct ﬂ uids: formation water in chemical equilibrium with the host rock, water rich in oxides and hydroxides, ascending hydrothermal ﬂ uid and corrosive water of meteoric origin. Alternate dolomitization and calcitization observed in the upper reservoir suggests alternate ﬂ ows of karstic freshwater and marine salt-water. Vertical, metre-long and centimetre wide open cracks are presently used for large water ﬂ ows; several generations of syntaxial calcite growth provide evidence for varying chemistry of the circulating ﬂ uids. Structural cross cutting relationships allowed us to establish a relative chronology of events, which can be correlated with the regional geodynamic evolution. The study reveals that the present-day reservoir architecture results from the superimposition of structures formed during the Early Cretaceous extension, Maastrichtian-Eocene Pyrenean shortening, and Oligocene rifting of the Gulf of Lion. The reservoir was also shaped by successive karsti ﬁ cation episodes and marine transgressions. Although the present-day hydrological system is controlled by, and reactivates structures inherited from a long-term evolution, it is characterised by frequent turn-overs of the water ﬂ ow, tuned by high-frequency external forcings such as sea-level changes driven by Pleistocene glacio-eustasy, or varying precipitation rates.


Introduction
Securing water resources is currently a major challenge, especially in the semi-arid to arid Mediterranean countries where population and economic development concentrate along the coastline. In southern France, excessive use of underground water resources may lead to seawater intrusion and salinisation of sedimentary coastal aquifers (Petelet-Giraud et al., 2016;de Montety et al., 2008). In the case of coastal karst aquifers, salinisation may also occur, by direct seawater inflow in the karst at surface (Drogue and Bidaux, 1986;Arfib and Gilli, 2010) or at depth through inherited karst network (Fleury et al., 2007;Gilli, 2020). Alternatively, or additionally, salinisation of coastal karst aquifers may result from mixing surface waters with fossil seawater, trapped in the reservoir (Khaska et al., 2013). It may also reflect interaction with deep saline thermal waters (Khaska et al., 2015). The interaction of a significant proportion (several percent) of deep brines with shallow karst freshwater is commonly described in karstic reservoirs (Bicalho et al., 2019). Also, upflow of fluids from deep reservoirs may bring warmer geothermal waters in the nearsurface karstic reservoirs (Cantarero et al., 2013;Petrini et al., 2013).
Mixing processes between chemically and thermally distinct fluids depends on the geological structure and texture of reservoirs. In tectonically complex and polyphase areas, it is pivotal to consider the 3D position of the segmented parts of the reservoir offset by faults, to understand and model hydrological systems (e.g. Husson et al., 2018).
The geological evolution of the reservoir, which involves processes such as sedimentation, subsidence, uplift, weathering and deformation, spans tens to hundreds of millions of years and controls the present-day hydrogeological setting. However, the hosted hydrological system, which depends on parameters such as sea-level change or amount of precipitation, is varying on much shorter time scales (10 4 to 10 5 years). Consequently, interactions between the long-term evolution of the reservoir structure and the short-term evolution of the hydrological system may induce periodic or temporary changes, and cumulative effects, while only the present-day hydrological system can be monitored.
The signature of past interactions between host-rock and fluids can be sought in the mineralization. Cross-cutting relationships of successive generations of veins provide a relative chronology (Cantarero et al., 2013). Developments of U-Pb dating of calcite filling veins (Roberts and Walker, 2016;Hansman et al., 2018) open new perspectives on absolute dating of tectonic events (Beaudoin et al., 2018;Parizot et al., 2021) and help reconstructing the evolution of paleohydrological systems (Pagel et al., 2018;Yang et al., 2022). This could help reconstructing the evolution of the paleohydrosystems.
Balaruc-les-Bains is the largest spa town in France (more than 50,000 visitors annually). It is located along the shoreline of the Gulf of Lion margin (NW Mediterranean), where Neogene post-rift marls onlap the inherited Mesozoic carbonate platform, along a complex topography that includes carbonate promontories, lagunas and low-land alluvial plains ( Fig. 1). Underground hydrology of the area (Aquilina et al., 2002;Pétré et al., 2020 and references therein) is characterised by the interaction of: (i) karstic freshwaters recycled from precipitations over the carbonate plateaus hinterland, (ii) marine and brackish waters percolating from the Mediterranean Sea and the Thau Laguna, respectively, and (iii) hydrothermal hot waters ascending from a deep reservoir. Water management must address the balance between: (i) supplying drinking water for a large population, which increases with summer tourism during the drought summer period, and (ii) preserving the quality of the hydrothermal water for the spa.
The overall objective of the Dem'Eaux Thau Project, which encompasses the present study, is to qualify and quantify the interactions of the different waters in the karst reservoir to improve water(s) management (Ladouche et al., 2019). An exploration well was drilled down to 750 m deep in the Jurassic carbonate reservoir in order to investigate the structure and the hydrological setting of the deep reservoir. The objective of the present contribution is to presents the preliminary results of the structural study of the 750 m long core sampled along the borehole. Core visual and microscopic analyses allowed to determine the structure of the reservoir and to understand the different post-depositional deformation, alterations and paleofluid circulations. Correlation of the successive events determined from the core analyses with the geodynamic evolution of the area helps to reconstructing the long-term geological evolution of the reservoir.

Geological and hydrogeological context
The studied karstic reservoir results from a polyphase tectonic evolution starting with the deposition of marine carbonates and marls on the NW Tethys margin, during the Jurassic and Valanginian (Debrand-Passard and Courbouleix, 1984). Bathonian to present formations are exposed in the Thau watershed (Gottis et al., 1967). The continuous Mesozoic marine sedimentation on the Tethys margin stopped during the middle part of the Cretaceous, as a result of uplift, erosion and weathering in continental environment, and led to the development of (i) a remarkable regional unconformity and (ii) bauxites formation (Combes, 1990). The latter was trapped in the karstified carbonates of Late Jurassic age, which makes a distinctive geodynamic marker across the study area (Lajoinie and Laville, 1979) (Fig. 1). This episode known as the "Durancian Uplift" (Gignoux, 1926) was recently dated in the study area to Albian . It is related to the late Early Cretaceous rifting of the Pyrenean domain (Lagabrielle et al., 2010;Duretz et al., 2020) and offshore Provence domain (Fournier et al., 2016), both located south of the study area. Symmetrically, the Durancian Uplift is also related to the kilometre-scale uplift of the hinterland, as rift shoulders (Barbarand et al., 2001;Marchand et al., 2020) in the Massif Central, north and northwest of the study area.
Later, the area underwent Pyrenean NS shortening, during Maastrichtian to Eocene times (Arthaud and Séguret, 1981;Arthaud and Laurent, 1995). The Pyrenean shortening is characterised by folding and thrusting of the Mesozoic carbonate cover, above Triassic detritals and evaporites decollement (Hemelsdaël et al., 2021b) (Fig. 1). Finally, the area recorded the Chattian-Aquitanian rifting of the Gulf of Lion margin, characterised by NW-SE extension producing half-grabens bounded by NE-trending normal faults, and filled by continental syntectonic breccia and marls (Benedicto et al., 1996;Benedicto et al., 1999). The study area does not expose synrift basins, but it is partly covered by post-rift Burdigalian marine sediments that onlap paleo-reliefs. The latter correspond to the hanging-walls of the previous Pyrenean thrusts (Hemelsdaël et al., 2021b). Post-rift sequence thickens offshore (to the SE) in the direction of increasing thermal subsidence following the opening of the Ligurian-Provencal oceanic back-arc basin (Séranne, 1999). The post-rift sequence is truncated by the Messinian erosional surface (MES), which developed in response to the desiccation-related drop of the Mediterranean sea-level (Clauzon, 1973). During the Messinian, the hydrological network adapted to the large scale base-level drop (1000 m) by: (i) incision of the detrital Cenozoic sediments (Clauzon, 1982) and (ii) deep karstic incision in the Mesozoic carbonates (Audra et al., 2004). The MES is observed onshore in the study area (Ambert et al., 1998) and it is filled by the Pliocene to Pleistocene progradational prism, made of marine clays, littoral to alluvial-plain marls and sandstones (Fuchey et al., 2001;Lofi et al., 2011). The Pliocene sequences present similar hydrologic properties than the underlying Burdigalian interval, i.e. both marl-dominated formations act as a seal for the underlying karstified Jurassic carbonates.
The hydrogeological setting is characterised by three distinct underground waters converging towards the studied Jurassic reservoir (Fig. 2) (Bérard, 1995;Ladouche et al., 2001;Aquilina et al., 2002;Aquilina et al., 2003;Pétré et al., 2020):  the carbonate reservoir can be intruded by saltwater wedge (Pinault et al., 2004), especially along the southeastern shoreline (Drogue and Bidaux, 1986). Remarkably, during specific time intervals, the Vise spring absorbs brackish waters, with downward fluxes that may exceed 300 l/s (Ladouche et al., 2021). This reversal, locally known as "inversac" (Gèze, 1987) has been documented seven times since the middle of the XXth century (Bérard, 1995) and has been continuously monitored since 2018 (flow, temperature, electrical conductivity, pressure, and chemistry; Lamotte, 2019). This new equipment allowed to accurately record the last "inversac" that started in November 2020 and stopped in March 2022. No periodicity or specific duration of the phenomenon has been evidenced. The documented reversals of flow occurred during periods of low freshwater piezometric level in the underground reservoir, whether it was natural (pluriannual drought periods) or artificial (e.g. pumping in neighbouring mines; Pétré et al., 2020); the third fluid consists in the saline thermal waters (30-50°C) that are exploited at Balaruc-les-Bains spa. Historically, thermal waters were naturally flowing from a spring, but it was later pumped from the top of the carbonate reservoir sealed by the Neogene marly sediments, through a series of production wells. Hydrogeochemistry suggests that thermal fluids are interacting with shallow meteoric water (Aquilina et al., 2002) and that the "pure" thermal end-member is a mixture of seawater and meteoric paleo-waters (Aquilina et al., 2002). Thermal waters are rising up from a deep reservoir, through a major reverse fault: the Thau Thrust (Figs. 1 and 2) (Pétré et al., 2020;Hemelsdaël et al., 2021b).

Boreholes
This contribution relies on a group of 4 vertical boreholes drilled from a single 10 Â 10 m platform in Balaruc-les-Bains. Wells #1 and #2 were bored by destructive drilling (45 m and 170 m depth, respectively) and were used for testing the hydrological behaviour of shallow reservoirs. Well #3 was drilled and entirely cored from the surface down to 300 m. Well #4, designed to explore the deepest part of the reservoirs, was bored at 3 m distance of well #3 by destructive drilling down to 290 m, then cored from 290 m down to 750 m depth. Considering the short distance between wells #3 an #4 and the 10 m overlap of the cores, which allowed accurate correlation, the two cored wells are integrated within one synthetic log, and interpreted as a single borehole from surface down to 750 m depth. Drilling was lubricated, and cuttings were extracted, with mud. Extra water outflows were recorded at the wellhead and interpreted as a result of the borehole crosscutting a water vein.

Cores
Cores diameters are 85 mm and 61 mm in the 0-200 m and the 200 m to bottom intervals, respectively. Cores are not oriented with respect to the North; the vertical cores display a very low dipping stratification plane, in agreement with the <10°dip observed in neighbouring outcrops. Examination of the cores took place on the site, where a high-resolution 1/100scale log was measured.
Sedimentary petrology and sedimentology were analysed (Fig. 3). However, the cores of the upper 400 m are characterised by an unsuspected amount of deformation, which often prevents classical diagnostic sedimentological observations. Two distinct deformation styles are identified: (i) penetrative brecciation displaying metre-wide strain gradient, (ii) discrete faulting; cores being naturally cut along the fault plane allow observation of the slickensides, which consistently indicate a normal sense of motion. Fractures, stylolites and veins are also analysed and recorded. Finally, macroscopic analyses of the cores allow record of alteration features and specifically karstification, with special attention to karst cavities infill.

Optical Borehole Imager
A series of logging tools were used in the borehole, however, this contribution focuses only on Optical Borehole Imager (OBI) (Gaillot et al., 2005;Gaillot et al., 2007) as a mean to orientate the structures. Indeed, optical images of the borehole acquired by OBI, which is associated with a magnetometer, displays planar structures as sinusoides. The amplitude of the sinusoid reflects the dip (knowing the borehole diameter), while the maxima and minima provide orientation with respect to the magnetic North. High quality optical imaging allows to correlate the sinusoids from OBI with the structures observed on cores (Fig. 4) where the sense of motion can be determined by direct observation of the slickensides. It was thus possible to analyse the kinematics of structures and discuss the orientation of the stresses that generated them. Well #3 was imaged from 101 m to 298 m and well #4 was imaged from 310 m to 750 m. Unfortunately, the 298-310 m interval could not be imaged.

Optical microscopy
Standard and polished thin-sections of core samples (100 mm thick) were examined under plane polarized light without impregnation. Classical optical microscopy with an Olympus BH2 microscope supported sedimentological and microstructural analyses.

Cathodoluminescence
Cathodoluminescence (CL) of calcite or dolomite crystals was used to investigate carbonate cements stratigraphy. Colors and intensities are commonly correlated with Mn2þ and Fe2þ content. It has been argued that high concentrations of Mn2þ activates luminescence, whereas concentrations of 1 to 1.5 wt% Fe2þ inhibits luminescence of calcite (Richter et al., 2003), although this usually need to be confirmed with additional geochemical data (Machel, 2000). In his

Alterations Discontinuities
Miocene marl filling karst at top of reservoir

Back Scattered Electron imaging
Back Scattered Electron (BSE) imaging was acquired at the MEA platform of the University of Montpellier, using a HitaschiS-2600N equipment. BSE imaging allowed analyses of deformation mechanisms and distinction of carbonate minerals. The latter was confirmed with qualitative elements maps generated by Energy Dispersive X-ray Spectrometry (EDS).

Structure of the reservoir 4.1 Lithostratigraphy
The lateral extent of the reservoir is controlled by integrating surface geology, seismic reflection, passive seismic imaging and gravimetry interpretations (Hemelsdaël et al., 2021a;Hemelsdaël et al., 2021b), and hydrogeological tracing and monitoring (Ladouche et al., 2019;Pétré et al., 2020). The 3 m apart boreholes #3 and #4 have been entirely cored, allowing the high resolution lithostratigraphic and structural analyses of one synthetic well-record presented here, from top to bottom (Fig. 3).
In the Balaruc area, the Jurassic karstic reservoir is covered by a 74 m thick Miocene sequence made of marine marl, interbedded with coarse to medium sandstone, and topped by several metres of Pliocene fluvial conglomerates. Although the Miocene conglomerate and sandstone contain an aquifer, the whole Mio-Pliocene interval acts as a non-permeable cover for the underlying limestone.
From 74 to 157 m, the borehole samples the Tithonian carbonates dipping about 10°to the NE. The top is affected by karstification displaying small-cale lapiaz filled with the overlying Miocene marls. These distinctive yellow marls are present in karstified fractures down to 90 m depth. The cores display a complex succession of light-grey oolitic grainstones interbedded with dolomitized intervals, which are characterised by a reddish, high-porosity, crumbly material, and presence of centimetre-size dissolution cavities. Several discontinuities in the core-record correspond to karstic cavities, correlated with water flows encountered during drilling. The major one (1.7 m hollow) is located at around 135 m depth. The base of this interval corresponds to an From 157 to 210 m, below the normal fault, the lithology consists of thick massive beds of grey-beige mudstone of Upper Kimmeridgian age, dated by correlation with surface exposure. The lithostratigraphic change on each side of the fault argues for several tens of metres offset across the fault. Below the mudstone, occurs a thicker interval of brown marly limestone, the basal part of which presents clayey interbeds up to 4 cm thick.
From 210 to 345 m, extends a circa 135 m thick, reddish, entirely dolomitized, interval. Up to 6 fine clay levels represent the marly interbeds of the original limestone. A distinctive feature is the penetrative brittle deformation that affects the whole interval, although with varying intensity (see below). At around 300 m depth, a major water flow has been recorded during drilling. Although no direct stratigraphic evidence could be found in the cores, its position below the upper Kimmeridgian suggests that this interval covers the Lower Kimmeridgian to the Upper Oxfordian.
From 345 to 440 m, are found alternating 5-10 m-thick beige to grey marly fossiliferous limestones and dark to black marls. The base is marked by a minor fault zone. This interval of Oxfordian age is classically considered as an lowpermeability unit (Dörfliger et al., 2008).
From 440 to 550 m, a massive beige mudstone displays 2 to 5 cm large karstic cavities, completely filled up with red siltstone. In the following, we name karst or karstic cavities, centimetre-to metre-size cavities that result from dissolution of lithified carbonate. This mudstone shows intense fracturing with calcite mineralization filling the fractures. It is locally dolomitized and it overlies marly limestones with thin clayey interbeds. This interval is of Upper Callovian age. When exposed in outcrop, the Lower Oxfordian corresponds to 20 m of distinctive brown chert limestones with a high silica content (Gottis et al., 1967). This formation is missing here, potentially due to the extensional fault observed at 440 m.
From 550 to 630 m, the cores record a thick interval of black marls with calcarenites at the base, related to the late Middle Callovian. The calcarenites comprise polygenic intraclasts derived from a carbonate platform and bioclasts. The proportion of carbonates increases upwards, and the transition with the overlying Upper Callovian limestone is progressive.
From 630 to the bottom at 747 m, core samples alternate between dark grey marly wackestones and black marls, in several metres thick packages. While drilling, limestone intervals have revealed the presence of several artesian water flows. Drilling and coring stopped within a >50 m thick black marls sequence. Comparison with neighbouring outcrops and correlation with the regional lithostratigraphic column suggest a Middle Callovian age.

Deformation
The carbonates of the reservoir display distinct modes of deformation, according to the dominant lithology.
Limestones are affected by bedding-plane stylolites. They are more frequent in the marly mudstone intervals (Fig. 5a, and 5b). They are associated with thin (typically 1 mm) and short (typically several centimetre long) sub-vertical calcite veins, or centimetre-thick calcite-cemented pull-apart zones connecting several stylolites (Rispoli, 1981). Uncommon vertical stylolites, showing horizontal teeth are also observed in the limestone. Brecciation in limestone results in sub-angular, centimetrescale, clasts of limestone and white sparitic calcite, floating within a white calcite cement (Fig. 5c). Reddish silt and clay may also be present between the clasts. Such breccia contains stylolites that affect both clasts and calcite cement. These overall bedding-parallel stylolites, are sometimes rotated into oblique shear-zones (Hausegger et al., 2010) (Fig. 5d). All observed features, i.e. stylolite teeth, orientation and rotation, or pull-apart zones infilled by calcite cement between oblique planes or stylolites, reveal a normal sense of shear displacement. Fault cores, typically several tens of centimetres wide, suggest a displacement of several tens of metres, according to the relationship between faultcores width and offset (Torabi and Berg, 2011). This type of deformation is particularly well expressed in the 157-210 m interval, in the footwall of the major fault, and, to a lesser extent, in the 440-550 m interval. In addition, such brecciation is observed locally along centimetre-scale fine shear-zones sub-parallel to bedding, in the 550-630 m interval (Middle Callovian) made of dark marly limestones and marls.
Another type of breccia is locally observed in limestones. It displays clasts with jigsaw geometry of the host-rock limestone and clasts of the earlier milky calcite veins, which are cemented by euhedral calcite. Such cement is observed penetrating into the host-rock, through multidirectional calcite small veins (Fig. 5e). Such features correspond to hydraulic brecciation of the limestone host-rock due to fluid pressure increase in veins or in the formation (Jébrak, 1997;Chauvet, 2019). Hydraulic fracturing is particularly well developed in the 450-500 m interval and also in several much thinner intervals ( 1 m) between 540 m and 570 m.
Finally, all limestone intervals are affected by extension fractures filled with different types of calcite (see below). They are dominantly sub-vertical and often cogenetic withbeddingparallel stylolites. Other vertical or steeply-dipping planar fractures are open or not completely filled. Although mostly observed in limestones, some centimetre-wide open fractures affect the dolomite (e.g. at 300 m).
Dolomite is affected by pervasive brittle deformation, which overprints all initial sedimentary textures. The main dolomite interval (210-345 m) is entirely damaged, and it displays gradients of deformation, ranging from weak fracturing, through crushed breccia, to ultra-cataclasite (classification of Sibson, 1977, Figs. 3 and 6). Low-strain pervasive fractures (0.01-0.1 mm wide) are randomly oriented and form centimetre-scale spacing networks that individualise small ( 1 cm) blocks of dolomite (Fig. 6a). Increasing amount of deformation results in dolomite micro-blocs individualization and rotation leading to a centimetre-scale puzzle-like texture (Fig. 6b). Next increase of deformation sees fracturing of the blocks accompanied by clasts translation and rotation (brecciation), within a matrix of reddish cataclased dolomite, representing up to 50% of the volume. In thin sections, fractures are seen filled with cataclastic dolomite material (Fig. 6c). Strain increase is expressed by increasing amount of fracturing and of cataclastic matrix (up to 75%), and clast size decrease (comminution, Fig. 6d). Maximum strain is expressed by cohesive, red, up to 10 cm wide, cataclastic to ultracataclastic shear zones, i.e. with no more than a few percent of host material preserved in the fault core. Thinsections confirm the intense cataclasis of the dolomite, with clasts size lower than 1 mm in the matrix. These clasts show rhombohedral shape and clast-size distribution displays gradient consistent with micro-scale shear-bands, suggesting that comminution is achieved by fracturing along the dolomite crystal lattice.

Kinematics
Rotation of the clasts into the cataclastic fault zones suggests a normal sense of motion. However, their low angle of dip (<40°) is less than the classical 60°dip of andersonian normal faults. The cores that are affected by such cataclastic zones generally split along the plane of maximum strain localization, thus revealing slickenlines with high values of rake (>70°). .
From the borehole optical imaging, we picked and compiled both orientation and dip of 42 cataclastic fault zones affecting the dolomite reservoir. On the basis of core examination, we extrapolate that all cataclastic fault zones record a normal sense of displacement and a 90°rake. Such a hypothesis is not consistent with stress-tensor inversion methods (e.g. Célérier et al., 2012), however, it allows us to determine the orientation of extension with the FaultKin8 TM software (Allmendinger et al., 2012) (Fig. 7a). Stereonet shows quite dispersed cataclastic fault zones orientations, however, tension axes display a maximum around a NE-SW trend and a pression axis close to vertical, consistent with NE-SW extensional tectonic event.
Similar method including a similar hypothesis, applied to the brecciated faults observed in the limestone intervals, provides a similar direction of extension (Fig. 7b).
The faults observed in both limestones and dolomites, therefore correspond to a NNE-SSW oriented extension, expressed by distinct deformation modes depending on the lithology. These observations thus point to a penetrative deformation related to an extensional event that affected several metres-thick sections distributed throughout the whole reservoir (see log on Fig. 3a).
Only one normal fault of significant offset is observed on cores and on optical imaging of both parallel boreholes (#3 and #4) that are located 3 m apart. The fault borehole images provide orientations of 068°30°N and 081°60°N, respectively. The significant mismatch in fault orientation in the two boreholes suggests an uneven fault surface, with an average ENE-WSW orientation and 45°dip to the north. In association with the 90°pitch of the slickensides observed on cores, this orientation suggests a NW-SE extension (Séranne, 1999).
Extension fractures are observed on the cores. They are cemented with calcite and sometimes display polyphase filling. Extension fractures are also identified in the borehole optical imaging, although the details of cementation (color, crystal habits) cannot be distinguished. They were picked and compiled with the WELLCAD TM software. Figure 8 represents the poles of 500 extension veins. The poles of veins are arranged in a girdle around the stereonet, which indicates high-angle to subvertical veins. There are also few poles in the central area of the stereonet, corresponding low-dip to subhorizontal structures. Cross-examination and correlation with the cores reveal that structure with lower angle of dip corresponds to either bedding parallel shear-planes coated with calcite, or, to reactivated and crystallised fractures, specifically in the dolomite interval. In order to simplify the analyses and interpretations, we have selected the extension veins dipping 60 to 90°, which represent the overall dominant population (n = 320). Figure 8b shows 3 main groups of poles of subvertical extension veins that emerge from the girdle: 1: N-S poles, corresponding to E-W extension veins (red on Fig. 8b); 2: E-W poles, corresponding to N-S extension veins (black on Fig. 8b); and 3: SE-NW poles (green on Fig. 8b). The E-W oriented extension veins are consistent with a N010°-N0190°extension; we therefore relate them with the extension of the same direction responsible for the formation of the normal faults. Observation on cores shows that the subvertical N-S extension veins overprint the faults zones, both brecciated fault cores in limestones and cataclastic faults cores in dolomites. No cross-cutting relationships could be observed with the major normal fault, however, the latter is consistent with the third cluster of tension gashes corresponding to NW-SE extension, and we suggest that they belong to the same extensional event.
5 Alteration of the carbonate reservoir

Karstification
The cores reveal three distinct types of karst, according to their filling and to their position within the reservoir.
The top of the Jurassic carbonates reservoir is karstified, with dissolution structures and extending 15 m downward (from À75 to À90 m). In addition to the metre-scale karst cavities that were detected during drilling operations (see above), the cores show evidence of Lapiaz, or karst cavities of several centimetres wide, are filled with the overlying marine sandy marls of Burdigalian age (Fig. 9a). Karstification was achieved in subaerial conditions, prior to the transgression, of which the marine sediments unconformably overlie and seal the palaeo-epikarst.
Deeper karst structures extend in the limestones of the upper reservoir (À120 to À210 m, Fig. 9b). Centimetre-wide dissolution cavities affect both deformed (faulted and brecciated) and pristine carbonate, as well as calcite veins. Such karsts are presently hydrologically active, as shown by water flows during drilling (e.g. around 137 m). Others karsts are plugged by polyphase filling, including pink laminated silts, Burdigalian marls and calcite cements.
The major dolomite reservoir (210-347 m) does not show evidence of karstification, although dolomite outcrops observed in the close surroundings of the borehole display metre-scale structures interpreted as resulting from dissolution.
Very deep karst structures are found in the À450 to À500 m limestone interval. This endokarst is characterised by 10 cm wide dissolution cavities, which are superimposed onto fractures partially filled with white coarse calcite spars (Fig. 9c). Cavities are clogged by bright red silty clays, which may contain fine clasts of oxi-hydroxydes derived from laterite crust. These observations indicate that (i) meteoric fluid-driven dissolution was guided by a previously acquired fracturation; (ii) the base level was lower than the observed interval, whether this resulted from sea-level drop and/or the carbonate massif was then at higher elevation; (iii) the karst network was affected by high-energy hydrodynamics able to carry surface sediments 500 m downward, into the fractured and karstified massif. Oxidising conditions also support a water-table level lower than the observed karst fill (Travé et al., 2021). Another type of karst infill is observed in that 450 to 500 m deep interval, in the upper Callovian formation. The lack of emersion structures (desiccation, brecciation, transgressive lag) in this marine interval indicates that such karsts are not early, syndepositional features. The karst filling consists of lithified laminated marl, of ochre to pink color, and which are older than the red continental sits karsts-infill (Fig. 9c). Finally, the carbonates of the upper reservoir displays connected, centimetre-scale cavities (Fig. 9e). Where dolomitized, the host-rock shows high porosity textures, whereas limestones of that interval shows cavities with corroded walls and no infill nor cement. While drilling, such intense dissolution intervals revealed water influx; the larger one was recorded at 132-134 m depth (Figs. 3 and 9e).

Dolomitization and calcitization
The cores reveal dolomitization affected two types of limestones. Firstly, the oolite grainstone of Tithonian age from the upper reservoir (Fig. 10a) shows multiple, complex and irregular dolomitization fronts (74 to 157 m). Such dolomitized intervals are found within the interval containing 0.1 to 1 mwide karst conduits that are associated with present-day water circulation, as detected during drilling. Thin-sections show the original oolitic porous texture preserved as ghosts within a texture of large crystals of dolomite (Fig. 10b). Calcite replacement by unimodal, small, planar euhedral crystals of dolomite, in several interdigitated zones within the shallowest part of the reservoir may suggests dolomitization driven by marine water-flow through the highly porous limestone (Vahrenkamp and Swart, 1994), but other kinds of dolomitizing fluids cannot be ruled out. Dedolomitization or dolomite calcitization (Coniglio, 2003;Radwan et al., 2021) is demonstrated by corrosion of dolomite crystals, subsequently coated by calcite growth and internal calcite zones within rhombs. Other rhomboedric crystals display alternate fringes of Ca and Mg distribution in SEM imaging, suggesting alternating calcite and dolomite (Fig. 10c-e). Such microscopic evidence of successive dolomitization/dedolomitization is observed across the whole interval of the upper reservoir (75-157 m), although it is more frequent in the top of the Jurassic reservoir, between 75 and 90 m. This oolitic upper reservoir is therefore highly porous and affected by a network of karst conduits, with evidence of past and present circulations of both dolomitizing fluid and freshwater (karstification).
Secondly, the lower reservoir between 210 and 335 m depth features a complete dolomitization of mudstone and clayey limestone of Lower Kimmeridgian and Upper Oxfordian age. Dolomite texture appears as unimodal, planar, euhedral small (10 to 50 mm) rhombs, in a matrix consisting of the initial micritic limestone (Fig. 10f). Dolomitization resulted in a compact, saccaroid rock, which is intensely fractured. Fractures are presently filled with fine crystallined calcite. No evidence of dedolomitization processes was observed in this lower reservoir. Drilling operations recorded evidence of present-day water circulation in one major, 5 cmwide, sub-vertical crack, at 298 m depth. This crack displays calcite crystallisation along the walls and remains open in the centre. The homogeneous dolomitization of the mudstone is constrained upward by clayey centimetre-thick interbeds that may have acted as seals. This suggests that dolomitization of the 130 m-thick interval, beneath a >150 m-thick cover, could result from a deep-seated process, linked to an ascending fluid flow 6 Signature of different fluids The mineralizations observed in veins and fractures are signatures of the different fluid flows within the carbonate reservoir. Analyses of the distribution of these mineralizations, as well as their geometric, macroscopic and microscopic features, observed under natural and cathodoluminescence light, provide clues about the relative chronology, the nature and origin of the fluids (e.g. Labaume et al., 2004;Travé et al., 2007;Laurent, 2015). Five types of crystallisation have been distinguished in the cores: calcite crystallised in the sub-vertical fractures associated with bedding-parallel stylolites presents the same redorange luminescence as the surrounding mudstone (Fig.  11a). This supports a closed hydrological system with short distance transfer by diffusion of calcite dissolved in the stylolites to the neighbouring fractures where it precipitated (Gratier et al., 1999;Toussaint et al., 2018); white scalenohedral calcite is observed filling geodes and centimetre-wide veins in the lower dolomitic reservoir (Fig. 12). Geodes usually appear at the top of brecciated zones. Dolomite puzzle-shaped clasts float in the calcite that fills the veins and secondary, smaller-scale calcite veins penetrate the host dolomite (Fig. 12c). This suggests hydraulic fracturing; (e.g. Jébrak, 1997). Scalenohedral calcite appears in cathodoluminescence with a dark-red color (Fig. 11b). The absence of varying luminescence zones indicates a constant chemistry of the mineralizing fluid throughout crystallisation, and/or a fast crystallisation process; orange calcite crystals are observed cementing the breccia developed in the damage zone of normal faults related to the NNE-SSW extension. In these intervals, orange calcite is also observed sealing the normal faults planes and filling the re-open sub-horizontal stylolites. It is also found coating the walls of one major crack (298 m). In plane polarized light, the orange calcite appears to incorporate significant concentrations of iron oxi-hydroxides (Fig.  11c). In CL, this calcite appears as mottled, with weakly luminescent areas. In rare occurrences, orange calcite can also be observed as a secondary cement, overprinting translucent calcite (see below), in centimetre-wide subvertical veins; translucent white calcite is observed filling up sub-vertical fractures and matrix porosity (Fig. 11d). In addition, this calcite often crystallises within dissolution structures, such as corrosion pits in the host limestone or in earlier calcite crystals. In plane polarized light, this calcite appears as large crystals filling completely the vein, whereas it is totally opaque in CL (Fig. 11d). Such black CL may indicate a precipitation of fresh-water (Labaume et al., 2004;Bussolotto et al., 2015), but additional geochemical analyses are needed to confirm this. Translucent calcite is found throughout the whole borehole, although it is prevalent in the upper reservoir (75-165 m); euhedral dogtooth calcite crystallised in subvertical, open, centimetre-wide cracks, presents successive growth generations, characterised by varying colors (white, light orange) sometimes separated by dark films. In natural light, this calcite appears as large clear crystals and the dark films correspond to patches of iron oxi-hydroxydes (Fig. 11e). In CL, the large crystals show a succession of thin zones that range from very luminescent to opaque, through dark red. The drastically different CL of successive calcite growth zones in the sub-vertical veins characterises temporal variation of fluid composition, and/or varying redox conditions in groundwaters at the time of growth (Meyers, 1991;Milodowski et al., 2018).
The cross-cutting relationships of the mineralized veins and fractures allows to establish a relative chronology of events and associated fluid-flows (Fig. 13).

Chronology of events
The observed structural relationships indicate a relative chronology of tectonic events while structural analyses allow to determine directions of extension and shortening, which can be correlated with regional geodynamics, in order to construct a calendar for the evolution of the reservoir (Fig. 14).
Following Late Jurassic carbonate deposition and early diagenesis (not addressed here), bed-parallel, compactionrelated, stylolites and related calcite-cemented veins are the earliest observed structures.
The whole reservoir, including compaction-related stylolites is affected by normal faults in limestones and cataclastic shear-bands in dolomites. Borehole Optical Imaging has allowed to orientate such structures corresponding to a NNE-SSW extension (Fig. 7), which represents the major and most distributed deformation across the reservoir. These extensional structures are consistent with the WNW-ESE-oriented extension veins and fractures (Fig. 8) and we interpret them as resulting from the same tectonic setting. Crystallization of the orange calcite cements the brecciated fault cores and clearly overprints bedding-parallel stylolites and associated calcite veins. We interpret the orange calcite cement as contemporaneous with these normal faults Such NNE-SSW direction of extension is not documented in the study area (Arthaud and Laurent, 1995). However, there are tectono-stratigraphic evidences of syn-bauxite deposition block-faulting in Languedoc (Hemelsdaël et al., 2021b). In addition, similar direction of extension has recently been recorded in the nearby Late Jurassic limestones that are karstified and filled by bauxite, dated to the Albian (Marchand,   cal.) in fault breccia comprising mudstone clasts, some of them including earlier calcite veins (cal.) and dark oxi-hydroxides elements; in CL, this calcite appears as non-homogeneous bright red areas including smaller darker or black zones; d: translucent calcite (trans.cal.) filling corrosion structures; e: dogtooth calcite (dog.cal.), the CL color bands indicate varying chemistry of the fluid, including bright orange (thermal), dark red and black (meteoric origin). CL photographs have been lightened by 10 to 30% to improve visibility but colors have not been modified. 2019; Marchand et al., 2020). Furthermore, direct U-Pb dating of calcite on slickensides of strike-slip and normal faults consistent with NE-SW extension, in several localities of Languedoc yield Aptian to Cenomanian ages (Parizot, 2021). We therefore relate the NNE-SSW extension to the late Early-Cretaceous pre-orogenic rifting of the Pyrenees and offshore Provence (e.g. Tavani et al., 2018 and references herein), simultaneously to the uplift of the south-European margin, locally known as the "Durancian Uplift" (e.g. Marchand et al., 2020 and references herein). In Languedoc, this uplift induced denudation of the Early Cretaceous carbonate cover by combined erosion and karstification, down to Malm sequences (Chanvry et al., 2020;Marchand et al., 2020). However, the studied cores do not provide evidence of karstification related to Durancian; in particular, the distinctive marker of this eventbauxitehas never been observed in the small karstic cavities of the cores. This is consistent with the palaeogeography derived from analyses of the neighbouring outcrops, where karsts filled with bauxite do not extend more than several metres below the bauxite deposits, which feature the palea-baselevel (Marchand, 2019). Cretaceous karstification has affected only the sedimentary interval above the Jurassic sequence sampled by the cores.
Extensional shear-bands and faults affect the dolomite of both the lower and upper reservoir. Furthermore, cataclasis and comminution of the dolomite material documented in the lower reservoir (Fig. 5d) indicates that dolomitization predates or is contemporaneous with the extensional tectonics corresponding to the "Durancian Uplift". Dolomite intervals in the upper reservoir display similar structural relationship and deformational evidences of a pre-tectonic dolomitization event. However, dolomite in the upper reservoir also reveals evidence of successive calcitization and dolomitization phases that are markedly absent in the lower dolomite reservoir.
The setting of the lower dolomite reservoir, covered by a marly limestone acting as a cover suggest deep dolomitization by an ascending fluid. This may correspond to a burial dolomite model, where basinal brines migrate from the deep and subsiding basin (Stacey et al., 2020), assisted by faults (Davies and Smith, 2006;Benjakul et al., 2020). The location of the study area on the western margin of the Tethys during Jurassic may account for the rise of thermal fluids (in the sense of (Machel and Lonnee, 2002) from the highly subsiding eastern part of the basin, during Jurassic times (Debrand-Passard and Courbouleix, 1984). Alternatively, dolomitization of the lower dolomite reservoir may be related to ascending basinal fluids related to the renewed subsidence that occurred south of the study area, during Valanginian-Barremian times (Séranne et al., 2002;Tavani et al., 2018;Marchand et al., 2020), prior to the extensional tectonics of the "Durancian Uplift". There is no unequivocal evidence that the pre-tectonic, first stage of dolomitization affecting the upper reservoir was contemporaneous. Timing of the additional calcitization/ dolomitization events overprinting the upper reservoir is poorly controlled.
The subvertical N-S-oriented veins that affect the whole reservoir (Fig. 8) are tentatively interpreted as a result of NS compression. The unoriented cores display some vertical stylolites and rare reverse micro-faults (unfortunately, not visible on borehole imaging, thus not oriented) that confirm the existence of a compressive episode. Such compressional structures are overprinted onto the Cretaceous extensional structures, as observed in the outcrops (Marchand, 2019;Hemelsdaël et al., 2021b). We correlate this deformation with the Pyrenean orogeny, which is the only shortening event that affected the Mesozoic carbonate platform of Languedoc, between the Maastrichtian and the Eocene (Arthaud and Laurent, 1995;Arthaud and Séguret, 1981;Hemelsdaël et al., 2021b). Folding and thrusting may have been accompanied by karstification of the uplifted and exposed Jurassic limestones. In the 450-500 m interval, cores show karstic cavities with laminated mudstones, later covered by laterite-related silts (Fig. 9d). The age of such deep karsts is questionable. Lack of Miocene-age tracers in the infill, such as the marine yellow marls observed in the shallower karsts, argues for a karstification phase older than Messinian. The laterite-derived infill corresponds to Late Cretaceous to Paleocene residual covers found in the area (Bruxelles, 2001;Séranne et al., 2002). The 500 m vertical amplitude of karstification suggests a combination of massif uplift and base level-drop, which could be correlated with the Maastrichtian-Early Paleocene initial phase of the Pyrenean shortening (Arthaud and Laurent, 1995;Hemelsdaël et al., 2021b) and the repeated base-level drops that have been documented during the Danian and the Selandian in the area (Combes et al., 2007;Husson et al., 2012;Husson, 2013).
Another episode of extensional tectonics is marked by a single major normal fault at 157 m depth, that offsets all previous tectonic structures. Its NE-strike and its relative chronology are consistent with the Oligocene-Early Miocene NW-SE rifting of the Gulf of Lion (Séranne, 1999;Hemelsdaël et al., 2021b).
Identification of Miocene marls in the karstified limestone as deep as 210 m provides clues for its evolution: karstification implies a base level drop of ≥210 m, associated with hydrodynamics able to transport downwards sediments deposited on the surface, and which postdates the Miocene. This is consistent with the Messinian base-level drop, which is known to have generated deep karst in the emerged carbonate platform of the area (Audra et al., 2004;Husson et al., 2018).
Finally, dissolution cavities in the upper reservoir overprints the Messinian-related karst (Fig. 9b). They are interpreted as the most recent karstification and the presently active water flows recorded during drilling, may suggest that the karstification process is still active.

Long-term evolution versus present-day water circulation
The structure of the carbonate reservoir therefore results from a long-term evolution, and the crystallisations attest of fluid circulations that have occurred from the Jurassic deposition and compaction to present-day. Interactions and/ or successive flows of marine, meteoric and ascending waters seem to have controlled the post-depositional long-term evolution of the reservoir; yet present-day fluid interactions involve the same type of waters. It is unlikely that the paleofluid circulations remained similar throughout the long-term evolution, until today. Indeed, the present-day fluid interactions are controlled by the large-scale structures of the reservoir, which have been acquired during successive stages. Specifically, present-day thermal waters are drained upward by the Thau Thrust, a Pyrenean (Eocene) structure (Pétré et al., 2020;Ranchoux et al., 2022;Ladouche et al., 2022); the thermal waters are sourced in a deeply buried limestone Jurassic reservoir Ladouche et al., 2022) that 3D geological modeling locates in the hanging-wall of the Nîmes fault, beneath more than 2 km of syn-and post-rift sequences, as a result of the Oligocene-Aquitanian rifting of the Gulf of Lion (Fig. 2) (Hemelsdaël et al., 2021b). Furthermore, the hydrosystem that produces the thermal waters at Balaruc-les-Bains cannot predate the youngest structural element of the reservoir, i.e. deposition of the Pliocene cover.

Dolomitization/dedolomitization
Our observations suggest that a first stage of dolomitization of the Malm limestone occurred during Early Cretaceous, probably in relation with the Durancian uplift. In the upper reservoir (75-157 m), preservation of texture in the oolitic carbonate suggests low-temperature dolomitization at shallow level, seawater being the source of Mg (Machel, 2004). In addition, the alternate dolomitization/dedolomitization processes observed only in this upper reservoir indicates superimposition of an alternation of meteoric, freshwater and marine water circulations. It could be hypothesised that it is related to the temporary flux reversals of the submarine Vise spring (Pétré et al., 2020). During normal (i.e. most frequent) setting, meteoric water is circulating in the upper reservoir, the carbonate host-rocks are karstified by CO 2 -rich waters, and freshwater is outpouring at the bottom of the laguna (Fig. 15). During reversal of the spring ("inversac"), brackish to marine water penetrates the reservoir through the Vise spring, and extends across a kilometre-wide area (Pétré et al., 2020). However, such reversals of the sub-marine spring are too shortterm to account for the dolomitization of the underlying carbonates, leached by the salty waters at temperatures less  than 25°C (Kaczmarek and Thornton, 2017) (Fig. 15). Alternatively, longer period oscillations, such as the Pleistocene glacio-eustasy sea-level changes, could be the driving mechanism for the oscillatory dolomitization and calcitization, as described by (Soreghan et al., 2000). During interglacial high-stand, as today, seal-level reaches the carbonate massifs outlining the coastline; temporal imbalance between karstic fresh water level and marine water pressure results in marinewater invasion of the reservoir and dolomitization. During glacial low-stands, which record up to 100 m sea-level drops, the shoreline regresses up to 70 km basinward (Rabineau et al., 2005). Consequently, the area of investigation is exclusively submitted to meteoric water, which reactivates karstification of the carbonates and allows calcitization, in the upper reservoir. Alternatively, the dolomitization-calcitization episodes affecting the upper reservoir may correlate with longer period events. Specifically, transgressions allow seawater to penetrate carbonates that were previously exposed and karstified (Machel, 2004;Moss and Tucker, 1996). The present-day large-scale structure of the study area was acquired after the Gulf of Lion rifting, in Oligocene-Aquitanian times. Since that time, marine transgressions have covered several times the exposed and previously karstified Jurassic limestones of the study area (Husson, 2013). During the Burdigalian post-rift transgression, sea-level rise was recorded by onlapping shore sediments as high as 80 m elevation, on the hills surrounding the study area (Combes et al., 2007). Later, following the Messinian desiccation event and deep karst formation, reflooding of the Mediterranean allowed deposition of littoral Pliocene sediments some 10 m higher than the present-day sealevel (Gottis et al., 1967). Therefore, the carbonates massifs of the study area and the Vise spring have been covered by marine water, which allowed salty water intake and interaction with the upper carbonate reservoir, leading to dolomitization (Fig. 15). Conversely, Messinian desiccation of the Mediterranean and glacial periods of the Pleistocene witnessed continentalization of the area, with formation or reactivation of karsts by circulation of meteoric waters and karstification of the exposed Jurassic limestone (Lofi et al., 2012).

Conclusion
The karst reservoir of Thau-Balaruc, located across the shoreline of the Gulf of Lion passive margin, concentrates three distinct water fluxes (karstic freshwater, marine water,    Fig. 15. a: conceptual model of alternating fresh and saline water circulation in the upper reservoir, responsible for successive episodes of dolomitization and calcitization of the carbonates. The lithological section in the centre corresponds to the studied borehole (same color code as for Fig. 3). Alternance of water circulation through the Vise spring in the Thau laguna and in the underlying reservoir could be tuned to the "inversac" phenomenon or (more likely), to longer period processes such as glacio-eustacy during Pleistocene, or regional tectonically-driven sea-level changes. See text for explanations; b: present-day setting of the Balaruc hydrosystem in relation with the architecture of the reservoir. The karstic waters circulate preferentially in the upper reservoir and through subvertical open crack in the dolomite reservoir; minor flows are recorded in the very deep reservoir. Upward flow of thermal water circulates across the dolomite reservoir through the same subvertical cracks and reaches surface level at natural thermal springs. and deep thermal water, each with important economic and social stakes. Hydrological modelling of this hydrosystem requires an improved knowledge of the reservoir, which has been addressed by drilling an exploratory 750 m deep borehole entirely cored. Preliminary analyses of the cores, specifically of the small-scale structures, of the karst fillings and of the calcite mineralizations has led to the following observations: the carbonate reservoir consists of Jurassic limestones covered by Neogene marls. Limestones have been subjected to two distinct types of dolomitization: (i) an early dolomitization of mudstones affected a lower reservoir, prior to a phase of pervasive cataclastic deformation and (ii) a Neogene to possibly Present dolomitization partially affected grainstones of the shallow reservoir; cross-cutting relationships of the structures and of the associated calcite mineralization allowed to establish a relative chronology, which can be correlated with the regional geodynamic evolution. Following the compaction-related subhorizontal stylolites and associated vertical calcite veins, the major deformation stage of the reservoir corresponds to a late Early-Cretaceous NNE-SSW extensional event related to the "Durancian uplift", overprinted by NS compression related to the Pyrenean event expressed by sub-vertical tension gashes. The basinscale geometry was finally acquired during the Oligocene-Aquitanian rifting of the Gulf of Lion; the study of calcite mineralization and chemical alterations of the reservoir shows similarities between the paleofluids that circulated through the reservoir and current fluids of the present hydrosystem. However, it is most unlikely that the latter is older than the Pliocene cover of the reservoir; dolomitization of the upper reservoir by Mg-rich seawater is favoured by alternate circulation of: (i) karstic freshwaters that create and develop a dense network of conduits as well as drastically increasing the contact surface for fluid-rock interaction and (ii) marine salty waters that penetrate deeply in the reservoir, using the karstic network; -450-500 m-deep paleo-karst filled with continental silts witness fluxes of meteoric corrosive water, driven by an even deeper outlet. Regional geodynamics suggests Paleocene and/or Messinian ages; the up-flowing geothermal water circulates across the lower dolomite reservoir through the Thau thrust, an inherited Pyrenean structure.
Our study indicates that the long-term geodynamic evolution has shaped the Jurassic karstic reservoir of the Balaruc peninsula, while the lithostratigraphic structures and deformation facies controlled the texture of this reservoir, and organised the connections within the current hydrogeological system. Although the present-day hydrosystem reactivates structures inherited from the long-term geodynamic evolution, it is subjected to turn-overs on a much faster rate, driven by high-frequency external forcings such as sea-level changes driven by Pleistocene glacio-eustasy, or climate changes controlling varying amounts of precipitation. Under the effects of the current climate change, the expected sea-level rise combined with aridification may raise concern for the future of the studied hydrosystem.
Further ongoing geochemical studies of the mineralizations will help constraining sources, residence time, paleotemperatures of paleo-fluids interacting with the host-rocks.