Abstract
Deeply buried sandstone reservoirs are targeted in the Upper Rhine Graben (URG) for geothermal and hydrocarbon resources. These reservoirs are affected by a convective heat flow along fault zones and have a complex diagenetic and deformation history recorded in their paragenetic sequence. Here, the focus is made on siderite and barite cementation characterisation, which trace paleo geothermal circulations within the fracture network affecting the Buntsandstein Gp. sandstones. A double approach on geochemistry and fracture network features is used to characterise fluid-flow episodes in the rift basin and on its shoulders. Barite sulphur isotopic signature suggests a common source for all the locations. However, Rare Earth Elements distribution patterns, oxygen isotopic ratios, and fluid inclusion study suggest two distinct flow regimes for fluids associated with barite precipitation along the shoulders and at depth in the middle of the graben. The barite has a higher content in total REE and contains non-saline fluid inclusions on the graben shoulders, suggesting that fluid circulations within the border faults interact with sulphate rich layers and precipitate at temperature above 150 °C. In deep-seated samples from the central part of the basin, barite fluid inclusions show a wide range of salinities, suggesting a higher contribution of sedimentary brines and precipitation at lower temperatures (< 150 °C). According to their REE signature, these barite mineralisations are associated with siderite and apatite with a diagenetic source. A conceptual model for fluid circulation within the basin is built from this new dataset. Fast and deep down- and up-flows occur along the major border faults, locally leaching evaporitic horizons. A part of the infiltrated meteoric waters reaches the centre of the basin, where it then mixes with the brines in sedimentary. This new characterisation of fluid pathways in the targeted reservoir brings insights into geothermal circulation compartmentalisation at the basin scale.
Résumé
Les réservoirs gréseux profonds du Graben du Rhin supérieur (URG) sont des cibles pour l’exploration et l’exploitation des ressources géothermiques et d’hydrocarbures. Ces réservoirs, affectés par des transferts convectifs de chaleur le long des zones de faille, ont une histoire diagénétique et de déformation complexe, enregistrée par la séquence paragénétique des fractures. Dans cette étude, l’accent est mis sur la caractérisation des cimentations de sidérite et de barytine qui retracent des paléo-circulations géothermiques au sein du réseau de fractures affectant les grès du Buntsandstein Gp. Ces minéralisations sont étudiées avec une double approche géochimique et structurale (failles et réseaux de fractures associés), afin de caractériser les épisodes d’écoulement de fluides pour différents éléments structuraux dans la partie centrale du fossé et sur ses épaules Une source commune est suggérée pour le soufre par la signature isotopique de la barytine, et e qu’importe la provenance des échantillons. La répartition des terres rares, les rapports isotopiques de l’oxygène et l’étude des inclusions fluides suggèrent cependant deux régimes d’écoulement des fluides minéralisateurs de la barytine. Sur les épaules du fossé, la barytine a une teneur plus élevée en terres rares et contient des inclusions de fluides non salins. Ceci suggère que les fluides au niveau des failles de la bordure du Graben interagissent avec les couches riches en sulfate et précipitent à des températures élevées (> 150 °C). Dans les réservoirs gréseux profonds, situés dans la partie centrale du fossé, les inclusions fluides dans les barytines présentent une large gamme de salinités et des températures plus basses de piégeage (< 150 °C), ce qui suggère une contribution plus importante des saumures sédimentaires. Ces minéralisations de barytine sont associées à des sidérites et à de l’apatite d’origine diagénétique, selon leur signature en terres rares. Ces données sont utilisées pour construire un modèle conceptuel de circulation des fluides dans le Graben, selon lequel l’eau météorique s’infiltre rapidement le long des failles bordières, jusqu’à la base du bassin, avec un lessivage local des horizons évaporitiques, puis une précipitation à partir de ces fluides géothermiques qui remontent durant l’activité de la faille. Une partie de ces fluides météoriques atteint le centre du bassin, où ces fluides ascendants se mélangent aux saumures sédimentaires. Cette nouvelle caractérisation des remplissages de fractures au sein de ce réservoir sédimentaire apporte des perspectives sur la compréhension de la compartimentation des circulations géothermales à l’échelle du bassin.
Introduction
Fluids are involved in widespread deformation processes in many geological contexts such as rift basins, foreland basins and oceanic ridges. Their origins are multiple: basement derived, connate or meteoric, and they are commonly mixtures between these sources (Person and Garven, 1992; Pribnow and Schellschmidt, 2000; Lampe et al., 2001; Bouch et al., 2006; Staude et al., 2009; Pfaff et al., 2010; Bons et al., 2014). For rift systems, the enhanced thermal regime at the rift initiation stage provides a high heat flow caused by lithospheric thinning (Ranalli and Rybach, 2005; Cloetingh et al., 2010). This thermal regime can be recorded in mineralisation in fault and fractures planes and in the matrix of sedimentary formations (Gleeson et al., 2001; Wilkinson, 2003; McKinley et al., 2011; Olivarius et al., 2015; Kristensen et al., 2016). The tracking of the origin of fluids is crucial to understand the processes that affect the hydraulic properties of the sedimentary formations and their evolution in space and time. These aspects have to be integrated into characterisation studies for mineral and non-mineral (hydrocarbon, water) resource exploration and production (Bense et al., 2013; Griffiths et al., 2016; Vidal and Genter, 2018). Such processes are proposed as significant for the development of Mississippi Valley Type ore-deposits (Paradis et al., 2007; Pfaff et al., 2010; Boiron et al., 2011). For geothermal energy purposes, the characterisation of the fluid origin, temperature, and migration pathways, is essential to estimate resources. For petroleum systems analysis, reservoir compartmentalisation can affect field productivity due to sedimentary features and fault zones architecture, including vein occurrence.
The fluid origin of the mineralisation hosted in rift basins is complex because of the diversity of fluid sources and fluid transfer conditions. Thus, it requires a double approach on geochemical and structural characterisation from the fault system scale to the micro-fracture scale. This study aims to characterise the origin and conditions of fracture mineralisation within the Buntsandstein Gp. sandstones within the Upper Rhine Graben (URG). These Permo-Triassic sandstones are a target for hydrocarbon and geothermal energy exploration and production. These sandstones are the first reservoirs at the bottom of the pre-rift sedimentary deposits (Böcker et al., 2016). They cover up the basement in several already well-studied locations (e.g. Soultz, Rittershoffen) (Vidal et al., 2015; Griffiths et al., 2016; Kushnir et al., 2018; Vidal and Genter, 2018). These sandstones outcrop on both shoulders of the graben; in the Vosges and Pfalz areas in the West, and the Schwarzwald and Odenwald Massifs in the East. Several studies have already discussed paleo- and current fluid-flows schema through faults in the area (Staude et al., 2009; Bons et al., 2014; Walter et al., 2016, 2018, 2019). Active thermal sources on the graben’s shoulders also indicate that fluid flow and mineralisation are still ongoing (Baatartsogt et al., 2007; Loges et al., 2012). However, these studies proposed mineralisation schemes in the basement and the Schwarzwald area’s sedimentary cover, which have a different Cenozoic evolution than the areas sampled here. These mineralising events were not explored within deep-seated sandstone reservoirs, which though are targets for geothermal energy production. This study brings new insights into the fluid types responsible for the fracture network cementation in the area to derisk reservoir operations.
The data collected here relate to barite, carbonate, and apatite precipitation in the pre, syn and post-rift phases of the basin. The elements obtained make it possible to discuss the source of the precipitated products, the nature of the fluids that permitted the transport of this material and the temperature conditions of precipitation.
By integrating previous works and new geochemical analysis in a structural context, the nature of the source(s) and the understanding of fluid flow regimes in the Buntsandstein Gp. sandstones is improved. After defining the geological context, focusing on the current knowledge of the fluid flow episodes affecting the Buntsandstein Gp., the sampled structural settings and the mineralisations geochemical signature will be presented.
The nature and fluid precipitation temperature in the fracture and fault vicinity and the basin-scale variability of the mineralisations will be subsequently compared to previous studies outcomings to elaborate a renewed model of flow pathways in the URG.
Fluid flow in fault zones in the Upper Rhine Graben
The Upper Rhine Graben is a segment of the European Cenozoic Rift System (ECRIS). It is a typical example of intracontinental rifting. It is bordered northward by the Hunsrück-Taunus Massif, and southward by the Jura (Fig. 1). This NNE striking graben is 300 km long, with a width of 30 to 40 km (Ziegler, 1992; Sissingh, 1998; Derer et al., 2005). The Vosges Mountains border the basin on the West, and on the East, the Odenwald and Schwarzwald Massifs. The entire area has a complex geological history involving several tectonic activity phases, resulting in fluid flows of different nature and origin (Fig. 2).
Structural settings
Several sets of faults and fractures networks affect the area, with a first-order compartmentalisation with NNE-SSW oriented faults, and a secondary network of fault and fracture striking N000°E, N040°E, N120°E and N150-N170°E, (Illies, 1972; Schumacher, 2002; Lopes Cardozo and Behrmann, 2006) (Fig. 1). In addition, the sub-basins segmenting the URG are controlled by pre-existing Variscan lineaments (Bertrand et al., 2018) affecting the basement (Edel et al., 2007; Skrzypek, 2011).
During Permian, the geodynamical context switched from relaxation and collapse of the Variscan Belt to an extensive intra-plate endorheic basin, with an intense tectonic activity (Schumacher, 2002; Ziegler, 2005). A series of basins developed, compartmentalised by reactivated NW-SE and NE-SW striking faults inherited from Variscan orogeny. The Buntsandstein Gp. sandstones were deposited in a high to moderate sinuosity fluvial system under semi-arid to arid conditions (Bourquin et al., 2006, 2009). These sandstones are overburden by Middle and Upper Triassic carbonate evaporitic deposits, characterising marine transgression and carbonate-evaporitic deposits. The thermal subsidence regime started in Triassic and continued up to Late Jurassic. A phase of uplift and erosion affected the area from Late Jurassic-Early Cretaceous up to Early Eocene. This erosion event resulted in a SE dipping of Mesozoic formations. The uplift centre was located NW-NNW of the current URG. This uplift phase was associated with volcanic activity in the Rhenish Massif (Lutz et al., 2010, 2013). The rift initiation started in the Eocene in response to the Alpine N-S compression and led to the deposit of the detrital Eocene Basal Formation. During Oligocene, the subsidence began to increase regionally. Marine ingressions flooded the whole basin, from the North Sea and the Alpine Sea (Roussé, 2006). These ingressions resulted in the deposition of clayey organic-rich carbonates, sandstones, and marls even across the graben margins.
Due to a change of the regional stress regime to a transtensive sinistral shear, and the initiation of NE-SW striking mini-basins in the late Oligocene and Miocene, the depocenters shifted to the northern part of the URG. The southern segment of the URG underwent a gradual uplift at this period, beginning in the Burdigalian (17 Ma). This uplift probably reactivated NE to ENE striking faults (Rotstein et al., 2005).
The current geometry of the URG settled during the middle Miocene uplift. Simultaneously, the northern part of the URG (north of Soultz-sous-Forêts) underwent a third phase of subsidence from Pliocene to Quaternary (Sissingh, 1998; Schumacher, 2002; Cloetingh et al., 2005, 2006; Bourgeois et al., 2007).
During this burial and uplift history, the Buntsandstein Gp. sandstones were affected by a series of tectonic deformation events accompanied by fault-related fluid circulations (Clauer et al., 2008; Blaise et al., 2016; Walter et al., 2018).
Fluid flows phases and geothermal activity
In the URG, the heat flow is relatively high (up to 150 mW/m2) (Lucazeau and Vasseur, 1989; Baillieux et al., 2013; Harlé et al., 2019). Local thermal anomalies are marked by temperatures exceeding 150 °C at 2000 m depth (Pribnow and Schellschmidt, 2000; Baillieux et al., 2013; Guillou-Frottier et al., 2013). These thermal anomalies affect the Buntsandstein Gp. sandstone reservoirs and are of crucial importance for geothermal resource operations. Their restricted spatial locations suggest that the URG temperature distribution is not homogeneous (Vidal, 2017). The isotherms established from well temperature measurements show positive anomalies that are often concentric and centered on the fault alignment, particularly at Soultz and Rittershoffen. The processes proposed to explain these anomalies and their structures are fluid flows at the basement-sediment interface (Dubois et al., 1996; Bächler et al., 2003; Baillieux et al., 2013; Guillou-Frottier et al., 2013). The geochemistry of brines and formation water also gives an idea of the first order of fluid flow regimes at the basement-sediment interface (Sanjuan et al., 2016). Thus, it appears that current fluid flows are structured at the basin and local scale.
At the basin scale, along E-W sections, the meteoric water infiltrates from the reliefs of both shoulders of the rift.
In the direction N-S to NNW-SSE, the NE-SW oriented faults and the sedimentary layers, which act as cover; control the first-order fluid flows (Meixner et al., 2018). Salinity data collected on fluids from the geothermal sites between Soultz and Landau suggest that they correspond to primary brines from water supersaturated in evaporites mixed with meteoric waters with very low salinity (Buntsandstein deep aquifer) (Sanjuan et al., 2010, 2016; Walter et al., 2018). The brines resulting from the dissolution of halite following the various transgression-regression cycles from Triassic to Oligocene are also a part of the fluid mixture (Sanjuan et al., 2016; Walter et al., 2016, 2018). Given the variations in geochemical properties of these current brines, several reservoirs are candidates (Walter et al., 2018), each with a specific type of fluids: crystalline basement aquifer, Buntsandstein aquifers, Muschelkalk aquifers, clay sequences, marls and evaporite levels of the Keuper, Jurassic limestone aquifers. Depending on the studies, these aquifers are not described with the same details. Some authors group the Middle Triassic and Upper Triassic under the name Triassic aquifer. Moreover, locally, tertiary aquifer influence is marked by sulphate isotopic signature (Staude et al., 2011; Loges et al., 2012).
The different temperature gradients and the results of geochemical studies of the geothermal fluids suggest the scheme illustrated on Figure 1. Convection loops might be present at the basement-sedimentary cover interface with recharge points on the graben shoulders, and upward migrations of sedimentary brine to horst structures (e.g. Soultz horst). The upward migrations of sedimentary brine might be consistant with thermal anomalies spotted along the URG (Soultz, Rittershoffen, Landau). The same drains control geothermal and hydrocarbon fluids migrations. These fluids’ deep circulations would be spatially limited to an area between Cronenbourg and Eschau for the southern boundary and Riedstadt-Groß-Gerau for the northern border. As it is located in this area, the Roemerberg oil-field (see location Fig. 1) is under the influence of these deep hot brines (Böcker et al., 2016).
The main reservoir for the hot brine feeding the geothermal sites would be situated in the eastern part of the graben, where temperatures reach 225 °C, and more than 4 km depth for the Buntsandstein Gp. (Sanjuan et al., 2010, 2013; Dezayes et al., 2015). These hot brines would then migrate from the graben centre to the northwestern edge, infiltrate the granitic basement, and rise upwards to the basement-sediment interface at the Soultz and Landau geothermal sites (Sanjuan et al., 2016). At the block (field) scale, local upwelling of the isotherms indicates local circulation systems, along some of the NE-SW orientation faults (Sanjuan et al., 2016). Many of the geochemical data from this study suggest a contribution of Cenozoic fluids and recent meteoric waters in Mesozoic and granitic basement aquifers. The potential connection between these two scales remains unknown due to the significant uncertainties about these fluid flow pathways and timing. For current fluid flow, a convection system, still active today in the URG basin, is also assessed by previous works (Sanjuan et al., 2010; Vidal et al., 2015; Freymark et al., 2017; Vidal and Genter, 2018).
A question remains regarding the evolution of these fluid circulation systems before the rifting, during the rift opening, and at the present time.
For the URG area, six major stages were reported (Fig. 2 and Tab. 1) (Staude et al., 2009, 2011; Loges et al., 2012; Walter et al., 2016, 2017, 2019; Burisch et al., 2017b; Dezayes and Lerouge, 2019):
Carboniferous stage: According to fluid inclusions in hydrothermal mineralisation and stable isotope data, Variscan fluids are generally of low salinity. Homogenisation temperatures range from 150 to 350 °C. There is no evidence that the Variscan fluid system was open to the surface, so after correction from lithostatic pressure, trapping temperatures reached 250–350 °C. Such conditions have been documented for similar veins along the Variscan Belt (Cathelineau et al., 2004; Schwinn et al., 2006; Boiron et al., 2011).
Permian stage: Permian fractures are filled by quartz, Sb–Ag, minor barite and fluorite assemblage and present aqueous fluid inclusions of two types, with varying salinities (Baatartsogt et al., 2007). Most fluid inclusions in these veins show low salinities (< 10 wt.%). Homogenisation temperatures range from 150 to 350 °C. The second minor group of inclusions shows calculated salinities ranging between 23.6 and 27.2 wt.%. Homogenisation temperatures range from 90 to 150 °C.
Jurassic-Early Cretaceous stage: The URG area records several Late Jurassic high-temperature hydrothermal events all over Western Europe. These events are linked to far-field stress changes related to Tethys and Atlantic rifting (Guillocheau et al., 2000; Ziegler et al., 2004; Staude et al., 2009; Brockamp et al., 2011; Cathelineau et al., 2012). These events are a possible explanation for local illitisation of Rotliegend and Buntsandstein sandstones (illitisation dated from 157 to 137 Ma) in the southern Schwarzwald and the Vosges Massifs (Brockamp and Clauer, 2005; Clauer et al., 2008; Brockamp et al., 2011). Barite–quartz–sulfide assemblages, dated from 150 to 110 Ma, also record these hydrothermal events overlapping the illitisation period partly (Schwinn et al., 2006; Danišík et al., 2010). The commonly presented genetic model consists in a deep down-flow of basin brines within the Variscan basement from Triassic to Jurassic times, where fluids reached equilibrium with the basement (Stober and Bucher, 2004; Baatartsogt et al., 2007; Staude et al., 2011; Bons and Gomez-Rivas, 2013; Bons et al., 2014; Walter et al., 2018, 2019). During uplift phases in the Cretaceous, the faults’ reactivation would drain upwards the Ba and F enriched brines that would mix with sedimentary basin brine rich in sulphate. This is marked by fluorite and then barite precipitations in fractures. In this model, fluids up-flows from the basement are very saline and hot, while sedimentary and meteoric waters are colder and of lower salinity. For some authors, precipitations in fault zones occurred continuously from Upper Permian to Upper Jurassic (Baatartsogt et al., 2007).
Late Cretaceous-Paleogene stage: This stage is related to the Paleogene rifting along NE-SW to NNE-SSW oriented fault systems (Staude et al., 2011; Burisch et al., 2017b). Barite–quartz–carbonates and barite–quartz–fluorite assemblages correspond to the main paragenesis observed in the Schwarzwald. Most of these mineralisation are accompanied by Pb ores, and more sparsely As, Zn, Cu, Bi, and Ni ores. The fluid salinity varies from 0 to 20 wt.%, and homogenisation temperatures from 50 to 150 °C (Staude et al., 2009). Several Paleogene Mississippi Valley Type (MVT) deposits are referenced on the eastern border of the URG (Schwinn et al., 2006; Baatartsogt et al., 2007; Pfaff et al., 2010).
Miocene-Pliocene stage: Salinity of 13–18 wt.%, and homogenisation temperatures from 50 to 150 °C (Staude et al., 2009) characterise the barite–quartz veins of this stage. Carbonate phases were also identified as late-stage mineralisations in the Schwarzwald district (Staude et al., 2012; Burisch et al., 2017b).
The fracture infills studied here likely belong to one or several of the three latter stages, as the studied sandstones of the Buntsandstein Gp. are younger than the two first phases. These pathways which control present day fluid mixing, might have been active since the Miocene, as the URG kept the same structural organisation and stayed under similar regional stress (Schumacher, 2002; Reicherter et al., 2008).
Material and methods
Few datasets with a detailed structural background are available on the western shoulder and deep-seated of fracture mineralisations in the Buntsandstein Gp., despite their contribution to fracture permeability, especially regarding barite (Griffiths et al., 2016). A multi-disciplinary approach is developed to determine the source of fluids causing mineralisation within the fractures and filling the matrix. This approach was applied on samples from different localities, referenced on Figures 1 and 3, and Table 2.
Backscattered electron (BSE) and cathode-luminescence (CL) observations were performed on a hot cathode on Tescan VEGA3 SEM device with a current intensity of 15–17 nA and a voltage of 15 kV.
18O/16O, 34S/32S (from SO2) were determined using a Thermo Scientific MAT 253 stable isotope ratio mass spectrometer system (CRPG laboratory). For each sample, barite crystals were isolated, then powdered, and 2 to 6 tin capsules containing 1 mg of the powdered material were analysed. The number of measurements per sample and the error associated with the measure is indicated in Table 3. In-situ34S/32S were determined using Secondary Ion Mass Spectrometer (SIMS), Cameca IMS 1270 ion microprobe at the CRPG laboratory. The results are expressed with the conventional δ notation vs V-SMOW (O) and CDT (S).
Bulk and in-situ δ34S on barite mineralisations were measured on 24 samples across the Upper Rhine Graben to trace the fluid’s potential source. Bulk isotope analyses on δ18O on the same samples coupled with fluid-inclusion microthermometry were used to estimate the temperature range of precipitation and deduce the fluid signature.
Fluid inclusions (FI) microthermometry was performed on an Olympus BX50 optical microscope, equipped with a Linkam MDS600 stage, within a temperature measurement range from −90 to 250 °C. Homogenisation temperatures (Th) were measured, and then samples were cooled down to −90 °C. Ice melting temperature (Tmice) was measured during the reheating phase.
Temperature does not impact the REE distribution pattern. The composition in REE is not affected by further fluid-circulations after crystallisation, except for systems with an extremely high water/rock ratio (Azmy et al., 2011; Tostevin et al., 2016). Globally, REE are substituting for the metallic ion in the carbonate lattice, and remain stable during diagenesis (Tostevin et al., 2016). REE chemistry is thus an excellent tool to decipher the fluid origin.
REE chemistry was performed on 21 samples using LA-ICPMS, to characterise barite, carbonate, and apatite mineralisations. The device used is a GEOLAS Pro nanosecond excimer laser, with a 193 nm wavelength. He aliments the system, and induced plasma is analysed by an AGILENT 7500 mass spectrometer (GeoRessources laboratory, Université de Lorraine). NIST glass was used as standard, and quantification was calculated using internals standards, which combined EDS measured composition (SDD type EDS spectrometer) and WDS analysis (Oxford Wave WDS spectrometer) on a JEOL J7600F SEM (SCMEM, Université de Lorraine).
Results
Petrography and mineralogy
The analysed sandstones consist of fine to medium-grained sandstones belonging to the Buntsandstein Gp. Detrital grains are mostly composed by quartz (> 85%), feldspars (av. 12%) and lithic fragments (av. 3%). The fracture sampled apparent length varies from 5 cm to 3 m, and their opening ranges from 100 μm to 2 cm. No difference in the major elementary composition between matrix and fracture siderite was detected from EDS analyses. Thus, all of these carbonates are qualified as siderite.
For the sample of granite, only the barite infill of the fracture was analysed. Macroscopically, the granite is a porphyric granite with amphibole and biotite (Granite des Ballons). An N095°E oriented structure carries the mineralisation (Fig. 3C).
In the entire samples set, cathodoluminescence observations allowed to identify zonations neither in barite nor in siderite. For the latter, the high content in Fe can explain the non-luminescence.
Siderite mineralisationis localised within the matrix (sometimes with a nodular shape), in oil-bearing partially open fractures, and fully cemented fractures. Apatite crystals were observed in samples from deeply buried sandstones only. Siderites as fracture infills (Sd2, Sd3) were observed only in samples from the centre of the basin (Roemerberg, Soultz, Lipsheim, Fig. 1). Samples from the URG shoulders and from sandstones that are currently deeply seated show barite mineralisation. Barite is mainly present as fracture infill (Figs. 4B–4E). Locally pervasive barite mineralisation occurs in Buntsandstein Gp. sandstones. When pervasive, barite mineralisation occupies pore space and seals quartz and feldspars overgrowths.
The structural feature orientations are as follows (Fig. 3D): the deformation bands have an N140-N170°E orientation and were observed only in deep-seated sandstones. Deformation bands pre-date the fracture features geometrically (Fig. 4C).
These three fracture cementation phases are related to abnormal diagenetic evolution. The orientation of the fractures cemented by barites is restrained to the N000-N010°E, N140-N170°E intervals for deep-seated sandstones, and N040°E, N000-N010°E, N140-N170°E within rift shoulders sandstones. Fractures containing barite remain mostly only partially cemented, with an average opening from 100 μm to 25 mm. The orientations of fractures cemented by siderite are N000-N010°E, N080-N090°E and N140-N170°E, and siderite often clogs the fracture porosity. The average fracture opening ranges from 50 μm to 1 cm.
The host-rock diagenesis consists of a feldspar and quartz syntaxial cementation phase, a matrix siderite cementation (Sd1) (Fig. 5), followed by illitisation and precipitation of apatite crystals in the illite mesh (Fig. 6A). This matrix diagenesis pre-dates the brittle deformation and fracture cementation by two siderite phases and one barite phase (Sd2, Bar, Sd3). Analysed sandstones are also affected by deformation bands (Fig. 4C) which concentrate apatite mineralisation (Fig. 6B).
The first fracture infills are constituted by quartz overgrowth on quartz detrital grains at the fracture border (phase I), based on the observed geometrical relationships between mineral phases. The second event is siderite precipitation (Sd2) (phase II), in the fractures, but also penetrative in the matrix pore network (Figs. 5E and 5F). The orientation of the fractures affected by these mineralising events is N000-N010°E, N080-N090°E and N140-N170°E. This generation of siderite also presents authigenic apatite (Ap2) (phase II) crystals which appear to be cogenetic of the siderite precipitation. This first fracture cementation is followed by barite precipitation, with needle or platy shape (phase III), that may build bridges spanning along the fracture width (Figs. 4A, 4B, 4D and 4E). Barite mineralisation is limited to fractures oriented N000-N010°E, N140-N170°E. Barite needles also trap apatite crystals (Ap3) of small size (< 10 μm length) in the fractures. The third generation of siderite (Sd3) plugs the fractures (phase IV), in the orientations N000-N010°E and N140-N170°E and traps cogenetic apatite (Ap4) (phase IV).
Isotopic composition
The source of sulphur influences δ34S signature. δ18O signature is influenced by the origin of the fluid and by the precipitation temperature of the mineral (Boschetti et al., 2011).
The δ34S and δ18O isotopic compositions of barite have a broad spectrum from 8.6 to 14‰ for δ18O and from 8.4 to 18.5‰ for 34S (Fig. 7 and Tab. 3). In-situ δ34S is in the same range as bulk isotopic composition measured on powdered samples.
La Fonderie’s sample (LFX005), from a vein in the carboniferous granite, has a δ34S value of 8.4‰. Samples from Reichenberg have δ34S values ranging from 13.4 to 15.3‰. Samples from Marlenheim outcrops have δ34S values ranging from 14.9 to 18.5‰. For Roemerberg samples, δ34S is varying between 13 and 16.2‰. For Soultz samples, δ34S is varying between 13.9 and 15.4‰. For Lipsheim samples, δ34S is varying between 16.9 and 17.4‰.
La Fonderie’s sample (LFX005) has a δ18O value of 8.6‰. For Roemerberg samples, δ18O is varying between 10.4 and 11.3‰. For Soultz samples, δ18O is varying between 9.2 and 10.9‰. For Lipsheim samples, δ18O is varying between 9.7 and 10.5‰. Samples from Marlenheim outcrops have δ18O values ranging from 11.2 to 14‰. Samples from Reichenberg have δ18O values ranging from 10.8 to 13.8‰.
Average δ34S values are similar for samples coming from the shoulders and those coming from the central part of the basin while the average δ18O values for the shoulder samples are higher than those measured for samples of the central part of the basin (Fig. 8). Despite the large interval of δ34S for each location, no zonation in δ34S content could be identified on barite infills nor by in-situ analysis or cathodoluminescence observations.
REE content and anomalies
REE pattern for barites is presented in Figure 8A. The Eu positive anomaly is likely to be an artefact due to the interference between Eu and Ba2+ during LA-ICPMS protocol. Therefore, these Eu anomalies for barite will not be discussed. The total REE amount in barite (ΣREE) is 16.5 ppm on average and ranges from 10.9 to 127 ppm. ΣREE is higher for barite from the shoulders than from the central URG (Fig. 8A). Using the LREE/HREE ratio expressed by the LaCN/LuCN ratios, the distinction between samples coming from the shoulders and those sampled in the deep part of the basin is made. Central basin’s samples are organised along a line with a LaCN/LuCN ratio ranging from 1.10−1 to 5.101 (Fig. 8B), while the shoulder samples are distributed above this line with higher LuCN values.
Siderites have a mean ΣREE of 55.9 ppm, and ΣREE ranges from 0.1 to 408 ppm. Three generations of siderite were analysed, and generations 2 and 3 have two different microstructural contexts, e.g. matrix cement or fracture cement (Figs. 9A and 9B). Crystals belonging to Siderite 1, which only occur as matrix cement, have a typical sedimentary REE pattern, with a LaCN/LuCN ratio relatively constant, around 10 for a large variability of LuCN value. Despite the few minerals sampled, the second pattern, constituted by matrix siderite 2 containing authigenic apatite, has a bell shape, with moderate content in HREE and MREE. These siderites have a LaCN/LuCN ratio between 4.10−1 and 3.
Siderite 2 sampled in fractures, which contain apatite have a bell-shaped pattern, with high content in MREE and HREE, but total ΣREE value is higher than Siderite 2 sampled in the matrix. These siderites have a LaCN/LuCN ratio evolving following a power-law versus LuCN, with LaCN/LuCN between 1.10−1 and 8, for a LuCN varying from 1 to 2.101.
For Siderite 3, sampled in fully cemented fractures, the REE patterns also have a bell-shaped pattern, with depletion of LREE compared to fractures containing oil. La content can be varying from 6.10−6 to 8.10−1 compared to chondrite standard. For the fractures containing oil, LaCN/LuCN range from 2.10−2 to 1.10−1 (Fig. 9B).
Fully cemented fracture cements have a LaCN/LuCN ratio ranging from 2.10−6 to 1, with most points ranging from 6.10−4 to 1, and they have the widest dispersion (Fig. 10). The REE pattern presents a bell shape on the three samples where apatite crystals are big enough to be measured, with high content in MREE (Eu: 10−3 ppmsample/ppmCN) (Fig. 10A). The total ΣREE content is higher than 103 compared to chondrite. The content in La is variable for apatite sampled in compaction bands. For apatites sampled in fractures, the LaCN/LuCN ratio ranges from 1.10−1 to 5.10−1 (Fig. 10B).
Microthermometry
Microthermometry was performed on both siderite and barite samples. However, measurements on siderites were not successful because of leakage during the freezing step. Thus, only barite Th and Tmice distributions for analysed samples are used (Fig. 11).
The measured fluid inclusions size varies from 5 to 15 μm long. Inclusions generally have ellipsoidal shape, either singular, isolated primary inclusion or forming pseudo-secondary fluid inclusion plans. The inclusions are aqueous, as Raman spectroscopy detected no CO2. In sample C 335 (Roemerberg), the pseudo secondary inclusions form fluid inclusions plans with co-existing hydrocarbon-bearing and aqueous inclusions.
In URG shoulder samples (e.g. X015 and X009 from Marlenheim, and B002 from Reichenberg), Th range from 140 to 240 °C, with a majority of Th restrained from 165 to 200 °C. Tmice ranges from −3.3 to −0.6 °C, reflecting very low salinities (from 1.1 to 5.4 wt.%eq. NaCl). Soultz’s sample (EPS11379) has Th ranging from 110 up to 165 °C, and Tmice from −10.2 to −0.2 °C. The fluids trapped in these inclusions have a low to moderate salinity (0.4 to 14.1 wt.%eq. NaCl). Samples from Lipsheim’s well (LIP 1495.28 and LIP1497.44) have a large interval of Th, from 80 up to 190 °C, with a major family of inclusions with av. Th of 140 °C. Tmice ranges from −24 to −0.1 °C, reflecting the fluid’s varying salinity in the inclusions (from 1.9 to 25 wt.%eq. NaCl). Roemerberg (C2335) sample has a large interval of Th from 100 to 150 °C, and Tmice from −13.1 to −9.7 °C. These Tmice reflect a salinity comprised between 13.6 and 17 wt.%eq. NaCl, equivalent to a moderate salinity of the fluids.
Two populations of inclusions emerge from the overall dataset. The first population with high Th ranging from 180 to 240 °C, with very low saline to non-saline fluids, differs from the second population of milder Th ranging from 80 to 150 °C, and with varying salinity, from very low salinity to high salinity. For the second population, the colder the Th is, the higher is the salinity.
The gaseous phases of inclusions within barite were also analysed with Raman spectrometry. H2O mainly composes the gas phases trapped. Neither CO2 nor H2S were detected.
The samples from deeply seated fractured sandstones show mostly inclusions from the second family, whereas the samples from the URG shoulders have inclusions belonging to the first described family of inclusions. The sample from Roemerberg (C2335) also shows primary and pseudo-secondary oil trapping inclusions and inclusions plans, with Th from 100 to 120 °C, av. 110 °C. These oil containing inclusions belong to the same fluid inclusion plans as aqueous inclusions, which Th range from 115 to 135 °C.
Discussion
The data collected on barite, siderite, and apatite precipitation in the pre, syn and post-rift phases of the basin provides constraints on the source of the precipitated products, the nature of the fluids that allowed the transport of this material and the temperature of precipitation. These new elements, integrated with previously published data, serve to reassess the URG fluid circulation model.
Paragenetic sequence
The paragenetic sequence is available in Figure 12. The pre-rift sequence divides into three stages (all included in phase 0):
1. A first planar clay grain coating, within fine-grained sandstones of the Buntsandstein Gp. (Sizun, 1995; Haffen, 2012; Soyk, 2015; Bossennec et al., 2018). Early cementation (Sd1) forms nodules and calcrete horizons associated with the development of soils on sand deposits (Haffen, 2012; Soyk, 2015). An early syntaxial feldspar overgrowth around detrital feldspar grains takes place (Soyk, 2015).
2. Precipitation of syntaxial quartz overgrowths on the border of detrital quartz grains, which is concomitant with feldspar illitisation (Soyk, 2015), and is probably fed by solutes originating from chemical compaction of quartz grains, which takes place in surrounding layers (Blaise et al., 2016).
3. Second feldspar overgrowth and authigenic illite precipitation (Bossennec et al., 2018). These last stage illites are tracing hydrothermal fluid circulations during Jurassic (Blaise et al., 2016; Clauer et al., 2008; Schleicher et al., 2006a, 2006b). The fracture infill observations suggest a paragenesis within the sandstones that completes the paragenetic sequence previously published. Four phases mark this fracture paragenesis.
Phase (I), with the deformation bands, forms first and present, in some samples, a concentration of authigenic apatite minerals (Fig. 6). Quartz overgrowths on quartz detrital on fracture borders also take place. During phase (II), siderite (Sd2) precipitates in every fracture direction and is penetrative into the surrounding matrix, cogenetic with apatite. At phase (III), barite precipitates, cogenetic with apatite, mainly in N000-N010°E, N040°E and N140-170°E orientations. Then, the last siderite (post-barite) fracture cementation occurs (phase IV), trapping the last apatite generation (N000-N010°E, N040°E and N140-170°E orientations). The lack of siderite in outcrop samples may be the sign of carbonates dissolution, from which Fe–Mn oxides remnants are observed in several outcrops (Soyk, 2015).
Fracture diagenesis is limited to quartz and barite precipitation on the shoulders, whereas the paragenetic sequence is complete within deep-seated sandstones.
Origin of the fluids
The new chemical data gathers various information on the fluids’ origin. Barite mineralisations are present on the URG shoulders and within deep-seated reservoirs, but have different precipitation contexts (temperature, fluid origin).
Barite isotopic signature
The δ34S composition of the samples ranges from 10 to 18‰, with an average of 15.5‰ (Fig. 7). These values are like the ones published by previous studies on hydrothermal mineralisations of barite, along the northern fault system separating the URG from the Taunus Mountains (Loges et al., 2012) and on fractures and veins sampled in the Schwarzwald (Schwinn et al., 2006). Several sulfur sources are identified and comprise sulfides and sulfates, such as early diagenetic pyrite in the Middle Triassic limestone with isotopic values from −36.0 to−7.6‰, Sulfate from the Middle Triassic, of which δ34S values are 18.5–21.6‰ (Staude et al., 2011) or sulfate from the Upper Triassic with values ranging from 14.3 to 17.4‰ (Boschetti et al., 2011; Boschetti, 2013). Values from Oligocene URG sediments are 11.2–13.8‰ according to Staude et al. (2011), but the 20–20.6‰ interval is typical for Tertiary sedimentary rocks, according to Loges et al. (2012).
In detail, analysed barite could differentiate into two groups using the δ18O values. Those with values lower than 12‰ δ18O, the samples from Lipsheim, Soultz and Roemerberg came from the deep part of the basin. Those with values upper than 12‰ δ18O came from the shoulders.
The jointed δ18O and δ34S signature of the barite cement are in accordance to Upper Triassic evaporitic sulphates origin, even if a light contribution of the Cenozoic sediments or primary sulfide is possible. This hypothesis of a Triassic sedimentary source for brines is also suggested for the Schwarzwald mineralisations (Bons et al., 2014; Walter et al., 2016). These evaporitic horizons are most probably the sulphur source for barite within fractures of the Buntsandstein Gp. sandstones, for mineralisations located on the rift shoulders, and within deep-seated sandstones. These conclusions agree with the current brines sampled in wells within the URG that show similar δ34S (Sanjuan et al., 2016).
Barite, Siderite and Apatite REE composition
The diagrams, presented Figure 13, are divided into 5 fields, describing the following conditions: In the field IIa, apparent negative Ce anomaly, due to positive La anomaly; in the field IIb, a negative La anomaly produces an apparent positive anomaly in Ce; in the field IIIa, the positive Ce anomaly is real; in the field IIIb, the negative Ce anomaly is real; in the field IV, the positive La anomaly is hiding a positive Ce anomaly.
The Ce anomaly for barites from the graben’s shoulders and the Central Graben reflects a considerable variation in positive and negative Ce anomalies. Sampled barites belong to domains IIa, IIIb and IV (Fig. 13A). The La anomaly is less positive for shoulders samples than for basin’s ones. Some barite crystals have a marine signature (from the Ce negative anomaly), and some barites have a positive Ce anomaly masked by a positive La anomaly.
In the central part of the basin, barite mineralisations are associated with siderites and apatites. A negative Ce anomaly signs the contribution of seawater derived brine as mineralising fluid, and the hidden positive La anomaly is a characteristic of authigenic carbonate (Dolníček et al., 2014).
The first generation of siderite, located in the matrix, is included in the domain IIa and has a typical diagenetic signature, issued from marine brine (negative Ce anomaly). These marine brines are likely to be infiltrated during the Late Triassic, as it is the case on the Schwartzwalder’s shoulder (Walter et al., 2019). Siderites of the next phases (Sd2 and Sd3) have Ce and La anomalies in the domain IV and sign their diagenetic character. The interactions with silici-clastic material from the host-rock may explain the less negative Ce anomaly (Carpentier et al., 2014). However, some samples located in fully cemented fractures belong to the domain IIIb.
None of the siderite sampled present a strong Eu anomaly, suggesting that they are not influenced by hydrothermal activity. Usually, for carbonates, a positive Eu anomaly is a characteristic feature of high-temperature hydrothermal precipitation (Bau et al., 2003; Tostevin et al., 2016).
These siderite mineralisations (Sd2 and Sd3) can contain apatite crystals. Positive anomalies in Ce and La mostly characterise sampled apatite crystals, with (Ce/Ce*SN, Pr/Pr*SN) belonging to domain IV (Fig. 13B). The La anomaly of apatite also hides here a positive Ce anomaly. A majority of analysed apatites crystals have Eu anomaly around 1, with only a few samples with positive Eu anomaly. This anomaly is typical for diagenetic apatite as Eu does not fractionate under 200 °C (Joosu et al., 2016). With an MREE enriched bell-shaped pattern, La, Ce, and Eu anomaly sign a diagenetic-sedimentary origin of the brine (Göb et al., 2013).
The main limiting criteria to form apatite is the presence of P, F and Ca. These elements can be supplied by feldspar alteration from the basement, from compaction processes during diagenesis, and felspar grain dissolution (McConnell, 1973; Shields and Stille, 2001; Joosu et al., 2016). The source and precipitation mechanisms of P and REE can explain the bell-shaped and the observed enrichment in MREE. During the burial, REE and P can be released in the pore waters by the reductive dissolution of Fe–Mn–oxides and oxyhydroxides (Haley et al., 2004; Joosu et al., 2016). Diagenetic phosphates REE spectrums often exhibit a bell-shape, due to MREE enrichments (Joosu et al., 2016; Tostevin et al., 2016). These oxides and oxyhydroxides can also retain Ce (Caetano et al., 2009), and are frequently associated with smectite on grain coating in the early stages of diagenesis of the Buntsandstein Gp. sandstones (Bossennec et al., 2018; Soyk, 2015). The burial under reducing conditions of these coatings could be a potential source for the element forming diagenetic apatite in sandstones (Soyk, 2015). This aspect of phosphor source for apatite needs further investigation.
Overall, this signature suggests a strong role of basin fluids as a source for barite, siderite and apatite mineralisations. The REE pattern and anomalies in La, Ce, and Eu on barite, siderites and apatites sign the main role of sedimentary brines as the source of matter for the mineralisations. Fluids with diagenetic or basement origins interacted and induce variations of REE signatures.
The fluid-mixing pattern is observed only in samples from deep-seated reservoirs, located in the buried part of the URG, e.g. Roemerberg, Soultz or Lipsheim samples. This observation suggests that sedimentary brines, influenced by seawater trapped in sedimentary formations, are more involved in the barite precipitation process in the centre of the basin than along its shoulders. The rift shoulders circulations are mainly driven by hot, low to non-saline fluids, and the mixing rate is lower in faults bordering the shoulders.
Fluid typing
Jointed analysis of the IF data and δ18Obarite measured on barite mineralisation, is used to determine the fluid origin and the temperature conditions for precipitation (Figs. 14 and 15).
The IF study allows to distinguish two different poles of fluids. Inclusions from the sampled fractures of the URG border faults, such as in Reichenberg or Marlenheim, are characterised by high Tmice (close to 0 °C) and high Th from 160 to 240 °C. Samples from deeply seated fractures present two types of fluid inclusions, a) “cold” Th (from 80 to 110 °C) and low Tmice (−22 to −10 °C), reflecting saline and mild temperature brines, and b) intermediate Th (140–150 °C) and with high Tmice (−5 to 0 °C), reflecting mild temperature and low salinity fluids. FI of samples from deep seating sandstones reflect the trend observed for barite in veins sampled in the Buntsandstein Gp. sandstone, and in the underlying basement (Dubois et al., 1996), and from hydrothermal carbonate mineralisation events reported for the end-member of low-salinity and 80–110 °C fluids (Pfaff et al., 2010).
The FI from samples located on the rift shoulders (e.g. Marlenheim, Reichenberg) reflect the hot end-member of hydrothermal calcite sampled in the Schwarzwald (Pfaff et al., 2010), the quartz in veins from Soultz (Dubois et al., 1996; Cathelineau and Boiron, 2010) the Schwarzwald quartz (Baatartsogt et al., 2007) and the “cold” low salinity end-member of post Variscan barite sampled in the Schwarzwald (Schwinn et al., 2006).
The Tmice of the studied samples do not reach the low values of Schwarzwald hydrothermal calcites (Staude et al., 2012) nor the pattern of barites from the Thuringer Basin (Majzlan et al., 2016).
Coupled with the FI study, the oxygen isotopic composition of barite helps estimate mineralising fluids’ isotopic composition.
Calculated temperatures of equilibrium are compared to measured fluid inclusion homogenisation temperatures (Fig. 14).
The δ18O of the analysed barite crystals have a wide dispersion from 9.2 to 13.8‰. Barites sampled in the deep-seated sandstones have a lower δ18O than the barites sampled on the URG shoulders
The identified fluids are:
1. low-saline fluids, with high Th (180 to 240 °C), which are associated with δ18Obarite from 9 to 17‰, and most probably have δ18Ofluid ranging from −10 to −5‰. These δ18Ofluid values are typical of hot, hydrothermal fluids;
2. higher salinity fluids, with colder Th (from 80 to 140 °C), and δ18Obarite from 9 to 13‰.
3. The deduced δ18O fluid is ranging from −5 to +3‰, which is typical for a fluid mixture composed by a) sedimentary brines usually having δ18Ofluid varying from 0 to +2‰, and b) meteoric derived fluids with δ18Ofluid below −5‰.
4. The first group of hot fluids, with low salinity, is a hot geothermal fluid of meteoric origin. Such fluid presence suggests a pulse signature of meteoric downflow deep in the basement, followed by the flow of hot geothermal fluids reaching URG shoulders, signed by the high temperature (around 200 °C) and low salinity (−4.6 to −1 °C) fluid inclusions. This type of fluids is also recorded within Lower Triassic sandstones in hydrothermal precipitations (quartz, fluorite, barite and siderite) on the eastern border of the URG (Brockamp and Clauer, 2005).
The second group of lowest temperatures, with variable salinity, is a mix between basin brines and meteoric waters.
Fluid pathways and basin conceptual model
Different fluid signatures depending on their structural locations were identified (Fig. 12). The different steps of the paragenesis are structured as follows:
1. A first phase composed by the pre-rifting diagenetic sequence: characterised by the development of siderite (Sd1) for which signal is recorded, as matrix cement, accompanied by feldspars and quartz overgrowths. This burial sequence recorded is similar in all locations.
2. A pre-rifting deformation phase characterised by deformation bands initiation. These deformation bands are observed within the deep-seated sandstones and concentrate diagenetic apatite crystals that record diagenetic fluids’ print within deformation zones.
3. The first fracture cementations occurred, with quartz overgrowth precipitation from the formation waters enriched in Si by diagenetic illitisation.
4. Siderites (Sd2) constitute the second fracture cementation phase and are located only within deep-seated reservoirs. These siderites precipitate from composite fluids, resulting from the mix of marine brines enriched in Fe and carbonate ions, and the more alkaline formation waters.
5. The barite precipitation records the third fracture cementation event. The geochemical signature of the barite is varying depending on the location of the samples. Within the basin, fluids are colder and could be more saline than for the shoulders. Sources are meteoric and basin brines for the deep-seated mineralisations, whereas sources are mainly meteoric and geothermal for the shoulders, with a small contribution of the basin brines. The presence of hydrocarbon-loaded secondary inclusions in some barite samples outlines that these barite mineralisations pre-date the reservoir oil charge, which occurred in Miocene (Böcker et al., 2016). This fluid flow event can be linked to the previously recorded Cenozoic hydrothermal pulses recorded within the Schwarzwald and Odenwald districts (named regional phase e) (Schwinn et al., 2006; Baatartsogt et al., 2007; Pfaff et al., 2010; Walter et al., 2016, 2019).
6. Finally, the last carbonate (Sd3) mineralisation phase seals barite cemented fractures for the deep seating basin samples (regional phase e).
A U/Pb dating of these mineralisation events would help to assess more precisely the window of fracture cementation depending on the location.
We propose the following model of fluid circulations in the URG, with an evolution of the fluid pathways versus time and space (Figs. 16 and 17).
The fluid pathway system needed to account for these different fluid circulations is also composite (Fig. 16). It is formed by the fault and fracture patterns and the sedimentary layers, especially the Triassic ones.
Schematically, following the basement structuration in blocks during the rift, deep, convective cells developed in the border fault zones with probably a branch flowing deeper towards the basement, through faults and fracture corridors. In the central part of the basin, fluid circulation is restricted to the bottom of the basin.
The preferential orientation of the cemented fractures, i.e. N000-N010°E, N040°E, N140-N170°E, and of fault systems (N000°E to N030°E) is compatible with their activation or reactivation in an E-W extensional regime, mainly Oligocene in age, or/and in a left-shearing trans-tensive regime, Miocene in age (Schumacher, 2002; Lopes Cardozo and Behrmann, 2006).
The main fault damage zones and the fractures lattices are usable for meteoric water infiltration on the shoulders. Fluids are drained downward in the major border faults and toward the basin using the oblique faults network. The biggest faults, probably those with deepest roots, are usable for deep hot fluid to move upward. Secondary fault network develops horst and graben structures with several hundred meters of offset. (Boiron et al., 2011; Sanjuan et al., 2013; Dezayes and Lerouge, 2019). These faults affect the deep-seated reservoir and interact with the sedimentary aquifer. The convective cells regulating geothermal circulations is an hypothesis already proposed in previous studies (Baillieux et al., 2013; Guillou-Frottier et al., 2013; Freymark et al., 2017). According to the geometry of the structural network, the dimensions of the convective flow “loops” in the border cells, is controlled by faults segments, with a 10 km length order of magnitude. Fault zones oblique to the rift borders can also drain fluid toward the basin, which fits with flow patterns proposed by Guillou-Frottier et al. (2013), as the oblique faults and damage zones crosscut the impermeable fault cores of the border faults.
In the central part of the basin, the fluid circulation involves quartz, barite and siderite precipitation. Migration of sedimentary brine along the sedimentary cover-basement interface agrees with previous studies (Cathelineau and Boiron, 2010; Genter et al., 2010). These sedimentary brines mix with up-flowing hot low saline waters on horst structures, as in Soultz (Smith et al., 1998; Dubois et al., 2000; Guillou-Frottier et al., 2013) or Roemerberg (this study).
A hydraulic gradient must exist to allow fluid transfers. One possible way is through the hydrostatic gradient induced by the topographic effect between the shoulders and the basin (Taillefer et al., 2018; Freymark et al., 2019). Such reasoning applies to the shoulders of the Rhine graben (Fig. 17). Before the rift opening, a regional scale uplift of the Rhenish Massif started in the Late Cretaceous and Paleocene. The Rhenish Massif acted since then, as a charging area for meteoric waters. The Vosges and the Black Forest massifs were uplifted during the rift opening (Eocene-Oligocene) and reactivated during the Miocene phase. Thus, these last two hydraulic charging areas were installed during the rift opening and reinforced during the Miocene tectonic reactivation phase. The recharge of the system by these topographic highs has been advocated regionally in the Schwarzwald, in the Hunsruck-Taunus and the Vosges Massifs (Fig. 17) (Cathelineau and Boiron, 2010; Sanjuan et al., 2010; Staude et al., 2011; Loges et al., 2012; Walter et al., 2016).
In deep-seated reservoirs, fluid pressure within fault-systems would control migration by setting varying permeabilities and fracture network cementation. This effect of fluid pressure would add on the main driving force controlled by the main border fault systems parametrisation and topographic gradient. The current study reveals that an up-flow of hot geothermal fluids within the border faults led to barite mineralisations. This result emphasises the hypothesis that these border faults do not only act as a down-flow drain but acted, at least in the past, as upward fluid pathways, and is complementary to previous studies conclusions (Sanjuan et al., 2016; Walter et al., 2018, 2019; Dezayes and Lerouge, 2019; Freymark et al., 2019).
Conclusions
The petrological and geochemical analysis provides tools to decipher the sources of fluid responsible for fracture mineralisation in the deep part of the URG and along its eastern and western shoulders.
The complex thermal history and subsequent flows of fluids within the Buntsandstein Gp. sandstone reservoirs are structured in distinct phases involving several fluids from different sources that were investigated in the Schwarzwald district (Danišík et al., 2010; Pfaff et al., 2010; Staude et al., 2011, 2012; Walter et al., 2016, 2018).
Regarding the condition of barite mineralisation, a common source for sulphate is the dissolution of Triassic evaporite. The precipitation regimes are varying following previous literature data. Shoulders show the warmest and least salty fluids, while fluids with variable salinity and relatively low temperature are in the central part of the basin. These basin fluids migrate through sedimentary drains in the deep part of the basin and permeable fault zones within the basement and potentially flow up on horst structures bordering structural blocks, such as Soultz or Landau. In the central buried reservoirs, these geothermal fluids mineralising barite are preceded and followed by the fracture network’s siderite cementation. The interaction between diagenetic brines and geothermal fluids is mainly restricted to fractures oriented N000-N010°E and N150-170°E.
These fluids pathways can only be active if the blocks architecture already existed. Given the age of the observed fault systems, the cementation of fracture network could start as soon as the Paleocene, and the structures located in the central parts of the basin, then probably appeared in the Oligocene-Miocene periods. The siderite and secondary apatite are mostly linked to diagenetic mineralisation of the open fractures. However, these fluid pathways and precipitation conditions are not fully understood within the deep central part of the basin and require further investigation.
Deep infiltrations of meteoric fluids explain graben shoulder mineralisations, topographically driven, and therefore the required presence of relief in the vicinity. This meteoric fluid circulation is restricted around major faults with a high throw, as they constitute the main drains, with a rise of hot, Ba-rich fluids, which leach out a portion of the sulphate-rich formations on the buried blocks, forming the barite deposits, in fractures, at high temperature.
The influence of the basin fluids is crucial from the early diagenesis on, and continues throughout the history of the sandstones, with a matrix and fracture diagenesis whose pre-rift (Jurassic), syn-rift (Oligocene) and post-rift (Miocene) phases are marked by the role of meteoric and geothermal fluids as a vector of the material taken by dissolution in the Upper Triassic formations and mainly stored in the Triassic formations.
The up- and down-flows within border faults and the deep circulations can also still be active at present, as suggested by the presence of a positive thermal anomaly, which is related to deep geothermal fluid circulation in an area between Strasbourg in the South and Roemerberg (Speyer) in the North (Baillieux et al., 2013; Vidal and Genter, 2018). The presence of similar hot geothermal fluid pulses and structures along the northern edge of the URG (Loges et al., 2012) and the similar nature of current brines could suggest that such fluid flow pathways between the basement and sedimentary cover are still active. They should then be considered precisely in geothermal reservoir characterisation because these flow pathways are the targeted resources and their induced modifications on the fracture network’s effective permeability.
Conflicts of interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Funding statement
PhD grant from ENGIE, Neptune Energy and the Université de Lorraine supports this research and publication.
Data availability
The supplementary data file includes the geochemical and structural data used to support this study’s findings.
Acknowledgements
We want to thank the geoscience team from ENGIE Deutschland GmbH (now Neptune Energy Germany) and their partner Palatina GeoCon to supply core material, scientific support, and rich debates. All the GeoRessources laboratory, Neptune Energy and ENGIE collaborators are thanked for numerous advices, settings supplies and help with measurements, the help and experience of the collaborators from the CRPG laboratory (Nancy). We would like to thank the ChronoEnvironnement laboratory in Besançon, for the supply of cold cathode facility, and discussions on diagenesis of fractured sandstones. We are grateful to the editor and the reviewers for their time and suggestions for improving the manuscript.
Cite this article as: Bossennec C, Géraud Y, Böcker J, Klug B, Mattioni L, Bertrand L, Moretti I. 2021. Characterisation of fluid flow conditions and paths in the Buntsandstein Gp. sandstones reservoirs, Upper Rhine Graben, BSGF - Earth Sciences Bulletin 192: 35.