The giant Jurassic-aged pockmark field of Beauvoisin developed in a 800 m wide depression for over 3.4 Ma during the Oxfordian; it formed below about 600 m water depth. It is composed of sub-sites organized in clusters and forming vertically stacked carbonate lenses encased in marls . This fine-scale study is focused on a detailed analysis of petrographical organization and geochemical signatures of crystals that grew up in early to late fractures of carbonate lenses, surrounding nodules, and tubes that fed them. The isotopic signature (C, O and Sr) shows that at least three different episodes of fluid migration participated to the mineralization processes. Most of the carbonates precipitated when biogenic seepage was active in the shallow subsurface during the Oxfordian. The second phase occurred relatively soon after burial during early Cretaceous and thermogenic fluids came probably from underlying Pliensbachian, Late Toarcian or Bajocian levels. The third phase is a bitumen-rich fluid probably related to these levels reaching the oil window during Mio-Pliocene. The fluids migrated through faults induced by the emplacement of Triassic-salt diapir of Propiac during the Late Jurassic and that remained polyphased drain structures over time.

Le pockmark géant de Beauvoisin s’est développé dans une dépression de 800 m de diamètre pendant plus de 3,4 Ma à l’Oxfordien. Il est composé de sous-sites formant des lentilles carbonatées empilées verticalement et inter-stratifiées avec des marnes. Cette étude est basée sur une analyse détaillée de l’organisation pétrographique et des signatures géochimiques des cristaux qui se sont développés dans les fractures précoces à tardives des lentilles de carbonate, des nodules environnants et des tubes qui les alimentaient. La signature isotopique (C, O et Sr) montre qu’au moins trois épisodes différents de migration de fluides ont participé aux processus de minéralisation. La plupart des carbonates ont précipité grâce aux apports de fluides biogéniques provenant des marnes sous-jacentes pendant l’Oxfordien. La deuxième phase s’est produite relativement peu de temps après l’enfouissement au début du Crétacé et ce sont des fluides thermogéniques provenant des niveaux plus profonds du Pliensbachien, du Toarcien supérieur ou du Bajocien qui l’ont traversé, sans savoir si ceux-ci ont atteint le paléo-fond de mer. La troisième phase est un fluide riche en bitume probablement lié à ces mêmes niveaux atteignant la fenêtre à huile pendant le Mio-Pliocène. Les fluides ont migré par des failles induites par la mise en place du diapir de sel du Trias de Propiac, mais ces conduits initiaux restent ouverts et sont donc réutilisés par des fluides plus tardifs liés à l’enfouissement de cette structure au court du temps.

At seep sites, sulfide and methane release results in authigenic carbonate precipitation commonly forming lenses encased in surrounding limestones or marlstones (e.g., Campbell et al., 2002). Seeping fluids sustain chemosynthetic bacteria and organisms at the seafloor, which have been recognized both in modern and ancient seep sites (Clari et al., 1994; Campbell and Bottjer, 1995b; Campbell, 2006; Kiel and Little, 2006; Gay et al., 2006; Amano et al., 2010; Kiel, 2010; Teichert and van de Schootbrugge, 2013; Vrijenhoeck, 2013). The fluid seeps are all related to an underlying focused fluid flow system connecting a reservoir at depth through a network of preferential pathways (Bohrmann et al., 1998; Aloisi et al., 2000; Aloisi et al., 2002; Agirrezabala et al., 2013). Depending on the type of reservoir and/or trap (i.e., turbiditic channels, salt diapirs, anticline structures, glaciomarine tills), the nature (marine, lacustrine) and depth of the source rock, and the residence time in the sedimentary column, the fluids can have contrasting signatures once they reach the seafloor (Gay et al., 2007). It also means that a carbonate body, such as a carbonate lens or a fluid pipe, can record several generations of fluids during burial. For example, this is the case in the Mariana forearc where calcite and aragonite precipitated at different stages during the formation of a chimney (Tran et al., 2014).

The giant fossil pockmark of Beauvoisin (Drôme department, Haute-Provence area, France) is one of the few examples worldwide being equivalent to Cenozoic chemosymbiotic-related fossil sites (Campbell and Bottjer, 1995a). It was first interpreted as a hydrothermal vent (Rolin, 1987; Gaillard et al., 1996) and then reinterpreted as a cold seep site (Peckmann et al., 1999; Kiel, 2013). It is now considered as a long-lasting (3.44 Ma) pockmark in the vicinity of Triassic salt diapirs, active during the Oxfordian (Tribovillard et al., 2013; Gay et al., 2019). The carbonate constructions have been also termed “pseudo-bioherms” by the pioneering studies of authors from the University of Lyon (France; e.g., Artru, 1972; Bourseau, 1977; Gaillard et al., 1985).

The isotopic signature of all minerals found in carbonate lenses and masses, nodules and tubes will provide information on the conditions of formation and the origin of fluids. In particular, the Strontium Isotope Stratigraphy (SIS) methodology (McArthur et al., 2012) was recently applied to fossil seep-deposits providing reliable formation ages (Ge and Jiang, 2013; Kiel et al., 2014; Kiel and Hansen, 2015). In this study, the mineralogical assemblages were investigated in details to date the different phases of fluid flow in the giant fossil pockmark of Beauvoisin. The aim of this paper is to provide the fluid flow history during burial of a seep site, based on petrographical organization and geochemical signatures in carbonate lenses, marls, nodules and tubes.

The carbonate, lucinid and gastropod-bearing, constructions of Beauvoisin developed under synsedimentary conditions, during the Jurassic in the South-Eastern Basin of France (Artru, 1972; Bourseau, 1977; Rolin, 1987; Gaillard et al., 1996; Peckmann et al., 1999; Kiel, 2013; Tribovillard et al., 2013; Gay et al., 2019). Prior to the development of the pseudo-bioherms, from the Lias to the Dogger the shallow-water carbonate platform underwent frequent subaerial exposure or erosion indicating a stable margin with a slow rate of subsidence (Lemoine, 1985). During the Bathonian, the moderate subsidence of the southern platform was compensated by higher rates of sedimentation, but the environments remained shallow. The Bathonian sequence is characterized by marls alternating with thin mudstone layers. However, the northern platform was already deepening as indicated by cherts and organic-rich marls facies (Dardeau, 1988). A general subsidence affected the entire platform from the Callovian and during the Oxfordian leading to deep marl facies deposition alternating with pluri-cm to m-thick dolomudstones layers in the northern part, locally remobilized by submarine slides. The giant fossil pockmark of Beauvoisin developed during the Oxfordian on the slope of the Liguro-Tethyan Ocean where the subsidence was at a maximum, leading to the deposition of up to 2000–2500 m of dark-grey organic-rich marls called the “Terres Noires” Formation (Gaillard et al., 1985; Gaillard and Rolin, 1988; Caillaud, 2018). The Late Oxfordian is characterized by the so-called Argovian sequence composed of thick amber-colored dolomudstone layers alternating with light brown marls in a general thickening upward sequence. Within the giant pockmark of Beauvoisin the thickness of the Argovian sequence varies by a factor of 3 to 4 which could be due to syn-to-post depositional erosion or slump processes on slopes of the depression (Gay et al., 2019). Subsidence was mostly due to basement faults with kinematics controlled by salt withdrawal in the extensional domain of the margin (Mascle, 1988). The giant pockmark of Beauvoisin is bounded to the west by major salt-rooted faults, which facilitated growth of the salt diapir of Propiac composed of Triassic evaporites (Fig. 1). The faults and the correlative salt diapir movements were active during the deposition of the Terres Noires in the area of Beauvoisin and some other places in the South-eastern Basin (Artru, 1972; Debrand-Passard, 1984; Tribovillard et al., 1987).

The site of Beauvoisin is characterized by 19 fossil-rich carbonate lenses embedded in fine-grained marls (sub-sites A to T, see Fig. 2 for log stratigraphy). Since its first recognition in the 1970s (Artru, 1972; Bourseau, 1977; Gaillard et al., 1985; Rolin, 1987, Gaillard and Rolin, 1988; Rolin et al., 1990), erosion has exposed new sub-sites that have been recently reported in a precise stratigraphical log (Gay et al., 2019) following the initial nomenclature established 20 years ago by Rolin (1987). Some key layers, such as the few cm-thick bentonite layers documented by Pellenard et al. (2003) or the R1 to R20 layers of the Argovian sequence (Fig. 2), are used as reference layers and correlated to the Tethyan sequences (Gradstein, 2012). This study is based on the analysis of the macroscopic and microscopic fabric of the carbonate lenses encased in the surrounding marls of the “Terres Noires” Formation. Twenty-four uncovered vertical 30 μm thin-sections regarding the carbonate lenses were prepared to perform the petrographic study, prior to geochemical analyses.

Analyses of C and O isotopes were conducted at the University of Barcelona on calcites extracted from uncovered thin-sections using a micro-drill to determine the carbon and oxygen stable isotope ratio of the different cements using the standard technique of Craig and Gordon (1965) and Claypool et al. (1980). The CO2 was extracted from 60 ± 10 μg of powdered carbonate samples, which were reacted with 103% phosphoric acid for 2 min at 70° C for calcites. The CO2 was analyzed using an automated Kiel Carbonate Device attached to a Thermal Ionization Mass Spectrometer Thermo Electron (Finnigan) MAT-252. The results are precise to ±0.02‰ for δ13C and ±0.04‰ for δ18O.

Strontium isotopic analyses were performed at University of Montpellier on calcites extracted from uncovered thin-sections using a micro-drill. Samples were cleaned during three cycles, in alternating rinsing of alcohol and MilliQ high purity de-ionized water for half an hour in an ultrasonic bath. Once cleaned and dried in a laminar flow hood, the samples were weighed in pre-cleaned Teflon beakers. The samples were then dissolved in closed Teflon beakers with a mixture of concentrated HF48% and HCl12N (0.5:2) and put on a hot plate at 110°C during 24 hours. After evaporation to dryness, HNO3 was added to the residue and kept at 110 °C for six hours before another evaporation step to dryness. Strontium isotopic compositions of representative samples were analyzed on unspiked samples. Chemical separation of Sr was done following two main steps. The first chemistry was done using the modified procedure of Richard et al. (1976) using a cationic exchange resin (Biorad AG50WX8). The objective was to concentrate only the alkaline earth elements and to remove the other elements from the matrix. The second chemistry was done using the Eichrom Sr-Spec ion-exchange resin in order to isolate Sr following a procedure modified from Pin et al. (1994). Total chemistry blanks were less than 20 pg and thus negligible for this study. Strontium isotopes were measured by thermal ionization mass spectrometry using a Triton Finnigan Mat spectrometer at Labogis (University of Nîmes). The samples were alternatively run with international NBS 987 standards using measurement procedure wherein standards were run every six unknowns. The 88Sr beam intensities for all standards and samples ranged from 10 V to 20 V. 87Sr/86Sr isotopic ratios were internally corrected from the instrumental bias using a value of 0.1194 for the 86Sr/88Sr ratio. Then the corrected ratios were normalized to the NBS 987 standards that gave a mean value of 0.710249 with a reproducibility of ±0.000006 (2 s, n = 5) during the course of this study. All ages are provided following the SIS (Strontium Isotope Stratigraphy) curve (McArthur et al., 2012) and the revised Jurassic strontium curve (Wierzbowski et al., 2017) as shown in the Figure 3.

The thin-sections were studied under optical microscopy, cathodoluminescence microscopy (CL), scanning electron microscopy (SEM) and Raman spectroscopy:

We used an Olympus BX 60 microscope, which allows polarized transmitted light and reflected light microscopic observation on the thin-sections. Conventional optical microscopy was supplemented by cathodoluminescence imaging, using a CITL MK4 optical Cold Cathode Luminescence stage fitted on an Olympus BX41 microscope. The accelerating voltage was 8–12 kV and the current of the electron gun was 200–300 μA. The images were acquired with a Spot Flex color camera of Diagnostic Instrument using a Kodak scientific grade CCD sensor of 4 megapixel. This cooled-high resolution-low light camera allows to image carbonate cements using a rather short exposure time.

SEM observations and analyses were conducted using a FEI Quanta 200 Environmental Scanning Electron Microscope (ESEM) fitted with an X-ray Energy Dispersive System (EDS). The tungsten electron source was at 20 kV. The spatial resolution of the EDS probe, a Roentec Single Drift Detector, is about 1 μm3. Pictures were taken using the back-scattered electrons (BSE) mode. The ESEM was used under high vacuum conditions with carbon-coated thin-sections to obtain better accuracy on the chemical analyses of minerals using the EDS system.

A Raman spectrometer (JY Horiba, Labram HR800UV) has been used to identify the mineral phases on some polished thin-sections. This apparatus is equipped with an electronically cooled charge-coupled device (CCD) detector, using 532 nm laser excitation. Spectra were obtained using 100× magnification objective. The acquisition time ranges between 20 and 60 s and the number of iteration ranges between 1 and 3 for each analysis.

The giant pockmark of Beauvoisin is characterized by 19 sub-sites (A to T) forming local unconformities in the Terres Noires Formation (see Gay et al., 2019 and references therein). Each sub-site is composed of sub-vertically stacked 0.5 to 2 m thick carbonate lenses surrounded by cm to dm nodules in a 100–120 m wide smooth depression 4 to 6 m deeper than the seabed (Gay et al., 2019). The nodules are encased in the carbonate lenses in the center of the depression whereas their size varies depending on the distance from the carbonate lenses. The sub-sites are organized in clusters grouped together in the same stratigraphic interval and the same geographic zone within the pockmark. The coalescence of all clusters (and their associated depressions) leads to the formation of a 800 m wide general depression interpreted as a giant fossil pockmark (Gay et al., 2019).

In this context, the sub-site E1 is very well exposed allowing a pseudo-3D view of the internal organization of an individual carbonate lens (see Fig. 2a for location). It developed during the Middle Oxfordian and it belongs to the Plicatilis zone and Antecedens sub-zone (Fig. 2b), as previously defined by Gaillard et al. (1985) and Rolin (1987), referring to the Tethyan sequences and age model from Gradstein (2012). The sub-site E1 forms a 4 m high edifice composed of sub-vertically stacked, 2–10 m wide and 0.5 to 1.5 m thick, massive carbonate lenses that are in contact (i.e., the top of an underlying lens is in sharp contact with the bottom of the overlying lens) or interbedded with nodule-rich marls (Fig. 4a). The carbonate lenses are composed of a micritic matrix rich in Mg-calcite, aragonite and dolomite, which are the common authigenic minerals of cold-seep sites (Roberts et al., 1993). Within an individual carbonate lens the nodule size is generally 1 to 2 cm at the base, decreasing to mm at the top whereas they can reach 15 cm in diameter in surrounding marls. Bitumen is frequent within the nodules encased in marls and not in nodules encased in carbonate lenses. As for sub-site F (see Gay et al., 2019 for details), the carbonate lenses are affected by breccia composed of mm to dm sub-angular clasts. The top of the upper lens of many sub-sites is affected by pervasive corrosion forming an irregular pluri-cm thick Mn-Fe-rich crust with pyrite inclusions. Remnant pyrite-encrusted micrite nodules can be found in overlying marls for a few cm to a few dm above the oxidized crust.

The carbonate lenses are dominated by endofauna, such as lucinid bivalves Beauvoisina carinata (Gaillard et al., 1985; Rolin et al., 1990; Gaillard et al., 1992), where they can form dense groups with some specimens reaching 18 cm in diameter (Rolin, 1987). Some other species have been identified in various abundance, such as gastropods, ammonites, crustacean exoskeletons, fragments and coprolites, fecal pellets, fish teeth, sponge spicules, echinoids, holothuroids, crinoids, benthic foraminifers, planktonic foraminifers, ostracods, radiolarians and dinoflagellates (Gaillard et al., 1985; Gaillard and Rolin, 1988; Rolin et al., 1990; Gaillard et al., 1992; Senowbari-Daryan et al., 2007; Kiel, 2010; Gaillard et al., 2011; Kiel, 2013). This biota is common in the Jurassic seafloor in the area suggesting a bathyal environment estimated around 600 m water depth (Tribovillard et al., 2013).

Most of the carbonate lenses are rich in macro-fauna, mainly lucinids, with a few exceptions. This is particularly the case towards the top of the sub-site E1. The carbonate lenses are lenticular in shape, meaning that they are laterally passing to the marls. On a vertical section (left view or front view on Fig. 4b) the thickest part of an overlying carbonate lens is systematically located where the underlying carbonate lens is thinning or absent (Fig. 4b). This suggests that the sub-site E1 is formed by the progressive lateral shift of the carbonate lenses in any direction through time with respect to the general stratigraphy (S0).

The sub-site R (see Fig. 2a for location) developed during the Late Oxfordian and belongs to the Bifurcatis zone and Stenocycloides sub-zone (Fig. 2b) referring to the Tethysian sequences and age model from Gradstein (2012). As for sub-site E1, the sub-site R is composed of vertically stacked carbonate masses that are lenticular in shape (Fig. 5). The front view (Fig. 5a1–a2) shows that the thickness of the carbonate lenses is laterally compensated, meaning that the sub-site R is formed by the successive lateral shift of carbonate lenses through time. On the left view the irregular bases and tops of some carbonate lenses laterally pass to more massive carbonate beds (Fig. 5c1–c2) concordant to the general stratigraphy (S0). These lenses are generally thicker than the other carbonates lenses, suggesting more favorable environmental conditions linked to better conditions of carbonate precipitation and/or preservation.

The top of the uppermost lens is affected by pervasive corrosion creating vugs and irregular surfaces. The corrosion surface is associated to a thick and massive Mn-Fe-rich carbonate crust (Fig. 6a and see Fig. 5c1 for location). This oxidized crust is in close association with numerous pyrite-rich nodules and vertical tubes, a few mm to 1 cm in diameter or height encased in the overlying marls (Fig. 6b and see Fig. 5c1 for location). The latter were never found elsewhere (i.e., laterally at the same stratigraphic level), suggesting a genetic relationship with the sub-site R.

Lucinid bivalves are the most prominent biota associated with the carbonate lenses of the sub-site R where they can locally form dense aggregates (Fig. 7a and see Fig. 4a1 for location). The shells can reach a diameter of 5 to 8 cm, suggesting that they fed on relative abundant food and nutrients (Rolin, 1987; Gay, 2002; Gay et al., 2019), but probably less abundant than on sub-site E1 where shells are larger. A thin-section realized in a lucinid-rich carbonate lens displays pluri-cm clasts made of yellowish or light-brown, microbial, micritic matrix, separated by veins (Fig. 7b, c). The micritic matrix is composed of detrital grains (mainly quartz) and pyrite framboids. it also contains some well-preserved shells and numerous mm to pluri-mm dark nodules coating body fossils, such as ammonites, bivalves, spicules and any biodetritus or burrows (Gaillard et al., 1985; Rolin, 1987). Preserved shells and more densely packed nodules are located in the vicinity of conduits and veins (Fig. 7b, c). The δ13C values can reach a minimum of −25‰ V-PDB on average in encased nodules (Peckmann et al., 1999), similar to the micritic matrix (Gay et al., 2019), whereas the δ13C values in shells are +5‰ V-PDB on average (Peckmann et al., 1999).

The basal lenses are systematically brecciated, forming various-sized, cemented, sub-angular clasts in a jigsaw-like structure, as for most of the other sub-sites in the area. However, they lack from macro-fauna, such as lucinids, or the shells are smaller than in overlying lenses. The brecciated facies is a matrix-supported cemented macro-breccia made up of cm to dm-size sub-angular micritic to microsparitic clasts within the carbonate matrix associated with abundant pyrite framboids (Fig. 8a, see Fig. 5a1 for sampling location). Thin-section observations in the basal lens shows that the brecciated facies is in close association with very sinuous mineralized veins filled by several generations of carbonate cements (Fig. 8b). In plane polarized light, a first generation of iron-rich micro-sparitic calcite (C1) is generally found in shell fragments or thin veins (Fig. 8c). This first generation of C1 is then followed by at least three generations of high Mg-calcite or aragonite (C2–C4) filling the veins (Fig. 8c).

The sample R-a was taken in a shell fragment containing the first generation of micro-sparitic calcite C1 (Fig. 8c). The isotopic analyses display a δ13C value of −14.82‰ V-PDB and a δ18O value of −4.69‰ V-PDB (Tab. 1). The measured 87Sr/86Sr ratio of 0.706855 ± 10.10−6 is consistent with an Early Oxfordian or Late Oxfordian age (Fig. 3).

Two samples (R-b and R-c) were taken at two different locations, in the middle of a thick vein hosting the following generation of a brown drusy calcite C2 (Fig. 8c). The calcite C2 is in contact with the micritic matrix and is usually coating clasts. The isotopic analyses (δ13C, δ18O and 87Sr/86Sr) are significantly different from sample R-a. Samples R-b and R-c display very close isotopic values with, respectively, δ13C values of +4.35‰ and +4.03‰ V-PDB and δ18O values of, respectively, +2.44‰ and +2.04‰ V-PDB (Tab. 1). The measured 87Sr/86Sr ratios for the two samples are homogeneous at 0,707069 ± 14.10−6 and are consistent with an Early Bathonian or Early Tithonian age (Fig. 3).

The sample R-d was taken in the middle of a thick vein hosting equant or blocky calcite filling all remaining space (Fig. 8c). Two successive generations of calcite were identified in plane polarized light and they were attributed to C3 and C4. Unfortunately the size of the micro-drill did not allow recovering samples at two different locations and the Rd sample is probably a mix between C3 and C4. The isotopic analyses display a δ13C value of +8.46‰ V-PDB and a δ18O value of +1.8‰ V-PDB (Tab. 1). The measured 87Sr/86Sr ratio of 0,707281 ± 14.10−6 is crosscutting the Sr curve several times (Pliensbachian, Late Toarcian, Early Bajocian, Late Berriasian, Late Aptian and even Late Turonian) as shown on Figure 3.

All δ13C and δ18O isotopic values obtained from samples R-a, R-b, R-c and R-d are in the range of values measured in previous studies (Peckmann et al., 1999; Tribovillard et al., 2013). However, no bitumen was detected in any samples of the sub-site R.

In the Beauvoisin area the sub-site A (See Fig. 2a for location) was first described as a columnar structure (Gaillard et al., 1985; Rolin, 1987). However, a recent study has shown that the sub-site A is composed of lenticular carbonate lenses with a general organization very similar to the other sub-sites B to T, and particularly to the sub-sites E1 and R previously described in this study (Fig. 9) (Gay et al., 2019). Three units A1, A2 and A3, respectively from base to top, were individualized as main sub-sites. They are located in the Early Oxfordian interval and thus correspond to the first bio-constructions occurring in the area (See Fig. 2b for log dating). They all belong to the Cordatum zone and the Cardioceras sub-zone after Gradstein (2012).

Three samples were collected in carbonate masses at various vertical positions related to the sub-site A1 (Fig. 8).

Basal carbonate lens: sample A1

The sample A1 was taken in the basal carbonate lens (see Fig. 9a and 9b for location), which is dominated by breccia and mineralized veins, at the same position as for sub-site R. In thin-section, the sample A1 displays a micritic to microsparitic matrix, rich in detrital grains (quartz, feldspar), shell fragments (mainly lucinids and gastropods), and pyrites (about 95% framboidal and 5% pseudo euhedral), crossed by mm- to pluri-mm wide mineralized veins (Fig. 10a). Depending on the concentration of bio-detritus, the matrix can locally be considered as biomicritic to biomicro-sparitic. The samples A1-a and A1-b were taken in the calcite C1 replacing shell fragments (Fig. 10a). In plane-polarized light, the calcite C1 is always crosscut by all other generations of calcite (Fig. 10b) supporting that the calcite C1 is the first generation of calcite occurring in the basal carbonate lens. The samples A1-a and A1-b have δ13C values of −9.82‰ and −12.01‰ V-PDB respectively and δ18O values of −4.16‰ and −6.08‰ V-PDB respectively (Tab. 1). The measured 87Sr/86Sr ratios of 0.706896 ± 10.10−6 and 0.706802 ± 14.10−6 for both samples are consistent with an Early Oxfordian or Middle Oxfordian age. The samples A1-c to A1-g were taken in the light brown drusy calcite C2 which is in contact with the micritic matrix in most places along the major vein or locally coating clasts. In details, the calcite C2 in contact with the matrix is microsparitic and becomes more drusy towards the center of the vein as shown in the plane-polarized section (Fig. 10b). Cathodoluminescence microscopic image does not exhibit any clear evidences for the limit between the matrix and the calcite C2 (Fig. 10b). The average value of δ13C is +16.54 V-PDB and the average value of δ18O is −0.51‰ V-PDB. The measured 87Sr/86Sr ratio of 0.707045 ± 10.10−6 for A1-c gives an age of Early Bathonian or Early Tithonian. The samples A1-h and A1-i were taken in the equant calcite C3 in the major vein (Fig. 10a). The calcite C3 is always in contact with the calcite C2 but never in contact with the micritic matrix as shown in the plane polarized section (Fig. 10b). Contrary to the natural light showing a difference in color between the light brown calcite C2 and the light yellowish calcite C3, the cathodoluminescence shows that the limit is not as clear (Fig. 10b). The drusy calcite C2 is locally passing to a more equant calcite C3. The average value of δ13C is +9.43 V-PDB and the average value of δ18O is 1.62‰ V-PDB. The average measured 87Sr/86Sr ratio of 0.707195 ± 14.10−6 crosscut the Sr curve several times (Pliensbachian, Middle Toarcian, Middle Bajocian, Early Berriasian and Late Aptian) as shown on Figure 3. The samples A1-j to A1-l were taken in the white calcite C4 in the middle of the main vein (Fig. 10a). Where the calcite C4 is in contact with the calcite C3 it appears fibrous (Fig. 10b). As shown in the cathodoluminescence, the calcite C4 may locally be the result of the recrystallization of C3 in thin veins where it becomes more microsparitic (Fig. 10b). Towards the center of the major vein, the calcite C4 becomes more equant or blocky, filling all remaining space. The average value of δ13C is +9.03‰ V-PDB and the average value of δ18O is 1.73‰ V-PDB. Only the sample A1-l was analyzed to get the 87Sr/86Sr ratio yielding a value of 0.707112 ± 14.10−6. This value crosscuts the Sr curve several times (Pliensbachian, Early Toarcian, Late Bajocian and Late Tithonian). In the sample A1, a last generation of a white microsparitic calcite (C5) was identified in thin veins crosscutting all previous generations of calcite (C1–C4) and the micritic matrix.

Carbonate nodule: sample N

The sample N was taken in a 15 cm wide carbonate mass at about 5 m beneath the basal lens of the sub-site A1 (see Fig. 9a and 9c for location). Its shape is very irregular and not well-rounded. It is part of an aggregate of several masses as shown in Gay et al. (2019). In cross section the carbonate mass displays a homogeneous micrite encasing some 1 to 3 cm size well-preserved shells and dispersed pluri-mm size nodules (Fig. 11a). Its composition is very similar to that of sample A1 but it appears richer in pyrite based on surface distribution. Very close to the edge of the mass the micritic matrix is brecciated by a 2 cm wide mineralized vein forming a jigsaw-like structure. Some clasts were incorporated into the mineralized vein as shown in thin-section (Fig. 11b). Even nodules in the close vicinity of the vein are brecciated. The samples N-a to N-d were analyzed in respect to the calcites C1 to C5 (Fig. 11b). For the sample N-a (C1) the value of δ13C is −9,68‰ V-PDB and the value of δ18O is −5,36‰ V-PDB. The measured 87Sr/86Sr ratio of 0.706830 ± 10.10−6 gives an estimated age of Early Oxfordian to Middle Oxfordian. The sample N-b (C2) gives a δ13C value of −4,12‰ V-PDB and a δ18O value of −0,54‰ V-PDB. The measured 87Sr/86Sr ratio of 0.707040 ± 10.10−6 gives an estimated age of Middle Bathonian or Early Tithonian. The sample N-c (C3) gives a δ13C value of −5,22‰ V-PDB and a δ18O value of −1,52‰ V-PDB. The measured 87Sr/86Sr ratio of 0.707264 ± 10.10−6 crosscuts the Sr curve several times (Pliensbachian, Late Toarcian, Early Bajocian, Late Berriasian and Late Aptian) as shown on Figure 3. Unfortunately the sample N-d, taken in the calcite C4 did not provide any data due to technical issues. As for the sample A1 in the overlying basal lens, the white microsparitic calcite C5 is crosscutting previous phases of calcite precipitation (Fig. 11b). However, it also appears as blocky calcite with 10 to 15 mm wide crystals filling all remaining space in the central conduit. Hydrocarbon has been locally found in the connected porosity between calcite C5 crystals. The sample N-e (C5) gives a δ13C value of −22,05‰ V-PDB and a δ18O value of 0,65‰ V-PDB. The measured 87Sr/86Sr ratio of 0.707629 ± 10.10−6 gives an estimated age of Trias, Early Sinemurian or Late Campanian (Fig. 3).

Carbonate tube: sample T

The sample T corresponds to a tubular carbonate mass, 30 cm high and 20 cm in diameter, found about 10 m beneath the basal lens of the sub-site A1 and about 5 m beneath the sample N (see Fig. 9a and 9d for location). It contains numerous nodules encased in a micritic matrix (Fig. 12a). Rare lucinid shells of 1 to 3 cm maximum size were also identified. It is characterized by a large number of mineralized veins and a 1 to 2 cm wide central vein in its center. A void is clearly visible both in thin-section and in cross-polarized illumination (Fig. 12b, c). This sample is affected by intense fracturing and defining the phases remained unclear in some places as the clasts and veins are crosscut several times by other veins. However, the surrounding marls are not affected by the fractures. In this context, the sample T-a was taken in a calcite C1 surrounded by the calcite C2. This sample T-a (C1) gives a δ13C value of −8,59‰ V-PDB and a δ18O value of −1,8‰ V-PDB. Locally, strontianite has been identified between the calcites C1 and C2. The sample T-b was taken in the drusy calcite C2, which is in direct contact with the micritic matrix. It gives a δ13C value of +4,55‰ V-PDB and a δ18O value of +0,47‰ V-PDB. The measured 87Sr/86Sr ratio of 0.707038 ± 7.10−6 gives an estimated age of Middle Bathonian or Late Kimmeridgian. In the tube, it is very difficult to distinguish the calcite C3 from the calcite C4 as they are not very well developed. However, both phases can be found in the middle part of veins in contact with C2. The sample T-c (C3–C4) has a δ13C value of +12,5‰ V-PDB and a δ18O value of +0,47‰ V-PDB. The measured 87Sr/86Sr ratio of 0.707250 ± 7.10−6 could be related to several ages (Pliensbachian, Late Toarcian, Early Bajocian, Late Berriasian, Late Aptian) as shown on Figure 3. As for the nodule N, a calcite (C5) has been identified in the middle of the conduit and in some few mm wide thin veins. In fact, this phase displays large crystals of 10 to 15 mm wide equant calcite growing towards the center where the remaining void is present (Fig. 12c). The sample T-d (C5) has a δ13C value of −19,63‰ V-PDB and a δ18O value of +0,42‰ V-PDB. The measured 87Sr/86Sr ratio of 0.707602 ± 7.10−6 gives an estimated age of Trias, Early Sinemurian or Late Campanian. The calcite C5 can locally contain traces of hydrocarbons or bitumen identified both on Raman and SEM, as previously shown by Peckmann et al. (1999) in this area. In two Early and Middle Oxfordian samples, the gas chromatograms (C13+) yielded n-alkanes with chain lengths from 12 to 33 carbons, indicating that the fluids originated from the thermal maturation of organic matter. They also show isotopic signatures typical of petroleum hydrocarbons derived from planktonic or detrital sources with δ13C values of −67,5‰ and −75,7‰ V-PDB (Peckmann et al., 1999).

Sub-site organization

A recent study on the same area has shown that a single sub-site (A to T) is generally organized as a set of vertically stacked carbonate lenses interbedded with marls or claystones (Gay et al., 2019). Each sub-site is located in the center of a concentric area of nodules forming a 100–120 m wide depression, 4 to 6 m in depth compared to the surrounding palaeo-seafloor. Although the carbonate lenses contain reef-building organisms it is now considered that they did not form any significant relief on the seafloor (Gaillard et al., 1985; Rolin, 1987; Peckmann et al., 1999; Gay, 2002) even if some examples on modern sites show that they can form a positive relief (up to 10 m) above the seabed (Teichert et al., 2005; Haas et al., 2010; Himmler et al., 2015). It has been shown that carbonate pavements formed in situ into the sediments (Bayon et al., 2009).

The sub-sites correspond to seven clusters of pseudo-bioherms and only one cluster was active at a given time. The seepage activity was spread over time, cluster after cluster. The spatio-temporal 3D reconstruction of the position of these clusters shows that the coalescence of the clusters led to the formation of a 800 m wide giant fossil pockmark in about 3.4 My. However, periods of fluid seeping alternated with periods of apparent quiescence, about 200 ky each (Gay et al., 2019). During the period of apparent quiescence, the reduced fluid flux in the shallow sub-surface could not sustain massive carbonate precipitation, but it does not mean that the fluid flux was also reduced at depth, with possible lateral migrations perturbating the zone beneath the giant pockmark, which gave birth to a new cluster of sub-sites at a different location in the general depression of the pockmark (Gay et al., 2019). To date all sub-sites were considered as individual pockmarks (Gaillard et al., 1985; Rolin, 1987) and carbonate lenses alternating with marls were considered to be the result of varying fluid fluxes as observed in modern sites (Feng et al., 2010). However, the well exposed sub-site E1 clearly shows that the carbonate lenses are lenticular in shape, laterally passing to the marls where nodules are abundant close to the carbonate lens. The thickest part of an overlying carbonate lens is systematically located where the underlying carbonate lens is thinning or absent, suggesting that the fluids allowing the massive carbonate precipitation were, at some point, deviated to be expelled laterally in any possible direction beneath the sub-site acting as a seal and leading to a new carbonate lens at seafloor. The basal lenses are composed of a breccia made of mineralized veins, clasts-rich conduits crosscutting limestones and embedded nodules whereas the overlying lenses are richer in macrofauna, such as lucinids, but are less affected by brecciation. This facies is very similar to the “brecciated limestone” from Italy (Ladanza et al., 2013). It means that the continuous fluid supply from below the sub-site could not shift laterally anymore. The overpressured fluids were trapped and could get through the basal lenses after having brecciated them.

Paragenesis

The detailed study of thin-sections evidenced similar paragenesis for sub-sites A1 and R, nodules N and tubes T (Fig. 13). They are all characterized by a dark brown micritic matrix containing well-preserved shells of various sizes, depending on their position in the carbonate lens, shell fragments made of calcite C1 (or aragonite) and disseminated pyrite (or marcasite) identified on Raman and on thin-sections. The isotopic analyses conducted on both sub-sites have shown that the isotopic signature of the calcite C1 in shell fragments are very similar with depletion in δ13C (−10,92 and −14,82‰ respectively) and in δ18O (−5,12 and −4,69‰ respectively). These values are consistent with recent geochemical analyses conducted on several Beauvoisin sub-sites that have led to the conclusion that they were dependent on chemosynthesis through AOM (Peckmann et al., 1999; Tribovillard et al., 2013). This conclusion is also supported by specific enrichments in Mo, As and Sb in carbonate lenses evidenced in the sub-site F (Tribovillard et al., 2013), or such as in the modern mud-volcano of the Malta Plateau (Cangemi et al., 2010).

Filamentous bacterial mats are of primary importance in cold seep sites where they form laminations, mineralized biofilms, crusts, tufs, peloidal fabrics, pyrite framboids, aragonite botryoids and crystal aggregates (Cavagna et al., 1999; Aharon, 2000; Clari and Martire, 2000; Peckmann et al., 2001; Barbieri and Cavalazzi, 2005). At the Beauvoisin sub-sites the activity of sulfate-reducing bacteria is also indicated by the dispersively-distributed framboidal pyrite identified both on thin-sections and in SEM (Scanning Electron Microscopy). Their concentration is particularly high within nodules N (Gaillard et al., 1992; Gay, 2002) and they are considered to be of diagenetic origin (Bréheret, 1994). These observations have led to the conclusion that pyrite and calcite C1 are contemporaneous and formed soon after the deposition of the micritic matrix.

The pyrite-C1 phase is then followed by a period of brecciation as shown by broken shells included in clasts. The shape of clasts is concordant with the massive matrix of the carbonate lenses suggesting a jigsaw-like structure (Gay et al., 2019). A dark brown cement-supported micro-breccia is more concentrated towards the veins and the clasts in the brecciated basal lenses are systematically clotted with a light brown drusy calcite C2. The calcite C2 displays a very different isotopic signature with positive values of δ13C (+ 17,04‰ at maximum for A1 and +4,55‰ at maximum for R) and less depleted values or slightly positive values of δ18O (−0,51‰ on average for A1 and +2,47‰ at maximum for R). These values are much closer to the isotopic signature of marls in the area both for site A1 and T respectively (Peckmann et al., 1999; Tribovillard et al., 2013). This not necessarily means that the fluid that initiated the first carbonate lenses and precipitated in shells to form the calcite C1 was different from the fluid that led to the brecciation and then the precipitation of the drusy calcite C2 around clasts, but this could mean that the second phase came from long-lasting fluids trapped beneath the sub-site and that became more in equilibrium with surrounding marls before they were expelled due to overpressure causing the brecciation. Locally strontianite was identified in replacement of calcite C2 in the nodules N and in the tube T. This suggests that an additional Sr-rich fluid may have caused the brecciation beneath the sub-site A1. However, this is not the case for the sub-site R where strontianite was not evidenced.

The subsequent yellowish equant calcite C3 and microsparitic to equant calcite C4 (locally dolomite or aragonite) grew up in contact with C2 towards the middle of the veins and they are never in direct contact with the micritic matrix or the calcite C1. The size and shape of crystals suggest slow processes of precipitation, probably due 1) to slow motion of fluids (compared to C2) within the veins in sub-sites A1 and R, the nodules in carbonate masses N and conduits in the tube T and 2) to the continuous supply of fluids of the same nature.

The last phase corresponds to a white microsparitic to blocky calcite C5. The calcite C5 was not identified in the sub-site R. In the carbonate lenses of the sub-site A1 the calcite C5 is only microsparitic in small straight veins crosscutting all previous phases C1 to C4. It is more blocky in the carbonate masses forming the nodules where it always filled all remaining space in veins. It means that the calcite C5 is always surrounded by the calcite C4. In some places, dolomite and aragonite have been identified through SEM observations, suggesting that the fluid was enriched in Mg at the end of the migration process. Barite and celestite were not identified in the nodules beneath the sub-site A1. However, they are present in nodules closer to sub-sites B to H corresponding to the clusters 3 and 4 in the new nomenclature (Gay et al., 2019). They often are in close association with bi-pyramidal quartz in nodules the type of which is called septarium (Guilhaumou et al., 1996) suggesting that the temperature at which they were formed was higher than 90 °C. These quartz often contain hydrocarbon inclusions (Touray and Barlier, 1975) indicating that hydrocarbon migration and quartz precipitation have occurred contemporaneously at burial depths greater than 3000 m given a standard heatflow of 30 °C/km in passive margins (Ranalli and Rybach, 2005). However, the heatflow may have been higher due to the vicinity of salt diapirs (Mello et al., 1995), in particular the diapir of Propiac only 2 km west of the studied area (See Fig. 1 for location). In the tubes T, few mm to cm wide voids are commonly interconnected, defining a very high porosity which is channelized within the conduit. In that case, the calcite C5 is always identified on the wall of the voids/conduits where it appears more equant.

Hydrocarbons, analyzed as bitumen (Peckmann et al., 1999), are trapped in the porosity between C5 crystals in the tube T, in the nodule N and in the carbonate lenses of sub sites A1. They are also trapped as inclusions in the calcite C5 and in quartz within nodules (Guilhaumou et al., 1996; Gay, 2002). However, no hydrocarbons were detected above the sub-site E (see Fig. 2 for location). Regionally, the P-T conditions have been estimated from clay mineral assemblages, vitrinite reflectance determinations (Barlier, 1974) and from aqueous and hydrocarbon-bearing fluid inclusions hosted by quartz crystals in the nodules and septaria (Touray and Barlier, 1975). There is a progressive transition from diagenetic evolution, with temperatures around 150 °C and pressures lower than 0.2 kbars, to low grade metamorphism with temperatures up to 250 °C and pressures near 0.64 kbars (Gaidon, 1988). The growth of single quartz crystals in the septaria was always synchronous with the thermal culmination. This probably relates to Oligocene sedimentary burial in the west of the basin, and to tectonic overburial caused by Alpine thrust-emplacement in the east (i.e. this study area) which resulted in the generation and the migration of hydrocarbons (Roure et al., 1994).

Generations of fluids

All samples (T-N-A1-R) of calcite C1 display isotopic δ13C and δ18O values between −16‰ and −7‰ V-PDB and between −1‰ and −8‰ V-PDB respectively (Fig. 14a), very similar to the isotopic values found in botryoidal aragonite from previous studies (Peckmann et al., 1999; Tribovillard et al., 2013). The C1 phase, including yellowish calcite C1 and botryoidal aragonite can be interpreted as bacterial in origin (Roberts et al., 1993) and negative isotopic values are a marker of the microbial mediation in the building of carbonate lenses (Peckmann et al., 1999; Louis-Schmid et al., 2007). The isotopic values of C2 to C4 phases are more dispersed with positive δ13C and δ18O values except for samples N-b and N-c that are closer to the C1 signature (Fig. 14a). The calcite C5 in T-d and N-e samples displays very negative values of δ13C between −24‰ and −19‰ V-PDB (Fig. 14a). These values are concordant with isotopic values measured in nodules from previous studies (Peckmann et al., 1999; Louis-Schmid et al., 2007), indicating that fluids that migrated through the tube and nodule had the same signature, probably related to hydrocarbon-bearing phase. However, the isotopic composition taken alone does sign neither the origin of fluids nor their relative occurrence in the history of the palaeo-pockmark.

Seafloor seep carbonates usually form in equilibrium with the Sr isotope composition of surrounding seawater and pore fluids (Naehr et al., 2000). However, methane seeps pose the potential risk that the seeping fluid carries Sr with a 87Sr/86Sr ratio quite different from that of ambient seawater (Paull et al., 1992), which could corrupt the utility of the Sr isotope signature of the seep carbonate for stratigraphic purposes (Kiel et al., 2014). Deviation from the expected seawater 87Sr/86Sr value possibly indicates mixing with deep fluids carrying a different isotopic signature, which can be modified during the upward migration by fluid-rock interactions (Teichert et al., 2005; Joseph et al., 2012). Recent studies conducted on Phanerozoic samples have shown that diagenetic alteration, not initial contamination, is the main issue in this approach (Kiel et al., 2014) and new examples of Caribbean fossil seep communities were successfully dated using strontium isotope stratigraphy (Kiel and Hansen, 2015).

In the Beauvoisin area, at least three groups can be identified using their signature on plots of 87Sr/86Sr versus δ18O and versus δ13C (Fig. 14b and 14c). Using the revised version of the Middle-Upper Jurassic strontium isotopic curve (Wierzbowski et al., 2017) as well as LOWESS 5 abacuses from McArthur et al. (2012), several ages were provided for each sample. The fluids expelled at the giant palaeo-pockmark of Beauvoisin hypothetically had a biogenic origin, namely the decomposition of the organic matter contained in the Terres Noires Formation (Gaillard et al., 1996). It is now confirmed by the Sr isotopic signature of the group C1. It means that the age of their formation corresponds to the biostratigraphic age of their host formation: Early to Middle Oxfordian for the nodule N and carbonate lenses in the sub-site A1 and Late Oxfordian for the sub-site R. It also means that they all grew up rapidly into sediments, avoiding any further interaction (and mixing) with later seawater even during shallow burial. This could be due to clay-rich sediments deposited in the Terres Noires Formation avoiding or limiting fluid exchanges. Contrary to hydrothermal sites where the source of fluids can be deep (Tran et al., 2014), the methane-rich fluids were originate locally from the host formation, interacting with sulfate coming from seawater through the anaerobic oxidation of methane (AOM). This first generation of fluids can be therefore considered to be biogenic.

The groups C2 and C3–C4 have very close Sr isotopic signatures (Fig. 14b and 14c). The group C2 is dated of the Early/Middle Bathonian or Early Tithonian whereas the group C3–C4 shows very similar ages of the Pliensbachian, Middle/Late Toarcian, Early/Late Bajocian, Late Tithonian, Early/Late Berriasian, Late Aptian or even Late Turonian. Given an average rate of sedimentation of 10 cm per 1000 years during the Late Oxfordian in the Beauvoisin area (Gay et al., 2019) and a 30% of compaction for shallow burial, the sub-site A1, nodule N and tube T were buried at about 1000 m below seafloor and the sub-site R was buried at about 750 m below seafloor when the Cretaceous started. This means that the Cretaceous seawaters could not interact with the buried features. However, deeper fluids may have provided with a Pliensbachian, Late Toarcian or Bajocian signatures. As shown by the paragenesis and the mineralogical assemblage, the groups C2 and C3–C4 were created shortly after the brecciation while the carbonate lenses and surrounding marls were still unconsolidated. This means that this episode of fluid migration probably occurred at a relative shallow burial (<1000 m) with no interactions with the overlying seawaters. Due to their isotopic signatures and the mineralogical continuum from light brown drusy calcite to yellowish equant calcite, the group C2 has evolved to group C3–C4 during the same event of fluid venting.

The group C5 has a Sr isotopic signature of Triassic, Sinemurian or Campanian (Fig. 14b and 14c). For similar reasons as for groups C2 and C3–C4, these samples were deeply buried during the Sinemurian or Campanian, meaning that the fluids may have a signature from Triassic or a mixed signature depending on the levels the fluids had to get through. Paleohydrothermal fluid circulation were initiated during the Lias and potentially occurred until the Bathonian along the Ardèche margin due to salt tectonics of the sedimentary cover over the Triassic salt (Lemoine et al., 1982; Mascle, 1988; Bonijoly et al., 1996; Aquilina et al., 1997; Aquilina and De Dreuzy, 2011). In the Diois and Baronnies geographic areas (including the Beauvoisin area) the main salt diapiric rising is dated of the Mio-Pliocene by means of K/Ar techniques (Bellon and Perthuisot, 1980). In the Propiac diapir (see Fig. 1 for location), the CO2/N2 ratio are the same in the primary inclusions measured in neoformed quartz from the diapiric breccia (Perthuisot and Guilhaumou, 1983) as in the inclusions measured in the quartz redeposited from brownish cap-rock, raft units, or even Triassic material which occur near the top of the diapir (Desmaison et al., 1988). These fluid inclusions were probably formed during a single fluid injection, not necessarily contemporaneous, however, with the thermal culmination of 300 °C (Edon, 1983; Edon et al., 1991). Fluorescence conducted on fluid inclusions coupled to stable carbon isotope analyses have shown the presence of oil containing n-alkanes generated from thermal maturation of organic matter from the sedimentary pile (Peckmann et al., 1999). The latter is consistent with the total thickness of the Terres Noires Formation which clearly reached a burial depth through the oil window in the lower part and with the Toarcian or Pliensbachian signature of fluids for C3 and C4.

Evolution of a fluid seep structure

In marine sediments, methane generated at greater depth and migrating upward interacts with sulfate at the base of the sulfate reduction zone and the top of the methanogenic zone (Martens and Berner, 1974). This sulfate–methane transition zone (SMTZ) is generally found at a depth of about 10 m below seafloor (t1 in Fig. 15), where the concentrations of sulfate and methane reach their minimum values (Martens and Berner, 1977; Iversen and Jørgensen, 1985; Alperin et al., 1988). Due to an increase in the methane supply from deeper levels during Early Oxfordian, the sulfate is depleted due to anaerobic oxidation of methane (AOM) (Borowski et al., 1999) and the SMI is shifted upward near the center of a pockmark (Gay et al., 2006). At this stage, the byproducts are sulfides in the form of H2S migrating above the SMI and shifting upward the oxic zone (t2 in Fig. 15). The center of the pockmark is characterized by reduced sediments and H2S combined with Fe ions to precipitate framboidal pyrite through bacterial mediation. Pyrite is disseminated in the micritic sediments coming from the hemipelagic sedimentation. Another byproduct of AOM is HCO3 ions that combined with Ca of seawater to form aragonite or Mg-rich calcite around nucleus such as sponge spicules or lithic fragments in growing nodules (Gaillard et al., 1992; Gay, 2002; Gay et al., 2019).

As the flux of methane is increasing, the SMTZ is shifted upward until the fluid can be expelled directly at the seabed (t3 in Fig. 15). In that case, methane seepage results in carbonate precipitation forming lenses or pavements encased in surrounding limestones or marlstones as shown in many fossil (Campbell and Bottjer, 1995b; Campbell et al., 2002; Kiel and Little, 2006; Kiel, 2010; Vrijenhoeck, 2013) and modern (Campbell, 2006; Amano et al., 2010; Teichert and van de Schootbrugge, 2013) basins worldwide. Sedimentary facies at the bottom are concentrically arranged, with methane-dependent features in the center of the pockmark and sulfide-dependent features at the periphery of the depression as shown in the modern Regab giant pockmark in the Lower Congo Basin (Gay et al., 2006).

It was shown that carbonate pavements formed in situ and grew up downward into the sediments (Bayon et al., 2009). It means that all porosity will be clogged by a microsparitic calcitic cement C1 forming the carbonate lens. This is also the case beneath the sub-site where the feeding conduits are progressively clogged forming tubes. The fluids are then trapped beneath the carbonate lens playing the role of a permeability barrier, as shown in modern seep sites (Bayon et al., 2009), and the fluids migrated laterally in any direction forming a new carbonate lens at the seafloor. As shown for the sub-site E1 (Fig. 4b), a sub-site is formed by compensation of the carbonate lenses due to the lateral shift of the fluid seep point (t4 and t5 in Fig. 15). All stacked carbonate lenses and their associated nodules formed in a general 120 m wide depression of about 4 to 6 m compared to the surrounding seafloor (Gay et al., 2019), and probably 5 to 8 m given a compaction rate of 30%.

Once a significant number of carbonate lenses had been formed, the fluids accumulated beneath the sub-site increasing the pore fluid pressure and the unconsolidated basal lenses were brecciated by hydraulic fracturing, allowing the fluids to be expelled at the seabed. The clasts in the breccia often host strontianite, which is always crystalized in contact with the micritic matrix of the clasts and the calcite C1 in shell fragments (t6 in Fig. 15).

At sites where AOM is active, methane concentration varies strongly depending on spatial and temporal availability of methane supply and microbial turnover rates (Knittel and Boetius, 2009; Deusner et al., 2014). Sulfate depletion is driven by methane oxidation in sediment pore-waters through bacterially-mediated reactions, inducing a change in redox conditions (Feng and Roberts, 2011) and favoring carbonate precipitation (Bayon et al., 2007; Ge et al., 2010; Vanneste et al., 2012, and references therein). However, once the focused biogenic fluid flow was reduced at depth, or because of the lateral shift of the fluid seeps during a period of relative quiescence (Gay et al., 2019), methane supply did not sustain the AOM at the seabed and both the SMTZ and the oxic zone were shifted downward, leading to the formation of an oxidized Mn-Fe-rich crust on top of the last (i.e., youngest and not buried yet) carbonate lens of a sub-site (t7 in Fig. 15), as previously shown for sub-sites F and G in the area (Gay et al., 2019). Finally the giant pockmark of Beauvoisin is the result of repeated cycles of fluid expulsion from t1 to t7 leading to the 19 stacked sub-sites A to T until Late Oxfordian (t8 in Fig. 15).

Based on the C and O isotopic signatures it was shown that a second phase of fluid flow came from long-lasting fluids that became more in equilibrium with surrounding marls. It is evidenced by continuous growing of crystals from light-brown drusy calcite C2 to yellowish microsparitic to equant calcites C3 and C4. The Sr ratios show that the fluid that induced the precipitation of the calcite C2 has a Sr isotope composition identical to Pliensbachian or Late Toarcian seawater. The fluids then migrated through pre-existing pipes, conduits and fractures (t9 in Fig. 15). However, no other fossil seep sites were identified in the surrounding Cretaceous intervals in the Beauvoisin area. So, the fluids may have been trapped leading to slow processes of crystallization for C2 and C3–C4.

A third phase of bitumen-rich fluids was identified in the Beauvoisin palaeo-pockmark (t10 in Fig. 15). It is probably related to deep fluids derived from the same Pliensbachian and/or Toarcian source rocks that have reached the oil window and the produced fluids have migrated through fault systems developed around the Triassic salt diapir of Propiac, which is located only 2 km westward. This has implications for fluid remobilization during burial as shown in the Adriatic foreland basin where salt withdrawal may be the main cause for oil migration (Ladanza et al., 2015).

Previous studies conducted in the Beauvoisin area showed that all sub-sites A to R developed in a 800 m wide depression for over 3.4 Ma, meaning that it can now be considered to be a giant fossil pockmark (Gay et al., 2019 and references therein). However, the origin of fluids remained unclear since the first identification of the Beauvoisin seep site in the 1980s (Gaillard et al., 1985). The detailed petrographic analyses of all mineralogical phases identified in a tube, in nodules and in the sub-sites A1, E1 and R have led to the conclusion that the evolution of a seep site is entirely recorded from its birth to its death and long after its burial:

  • the increase of methane-rich fluids flow to the seafloor leads to an upward shift of diagenetic horizons, such as the sulfate-methane interface and the oxic zone, allowing carbonate precipitation at the seafloor through AOM;

  • a carbonate lens probably represents a permeability barrier and fluids have to laterally shift to find a new seep pathway, forming a new carbonate lens by compensation. A sub-site is composed of all compensated carbonate lenses. If the fluids are trapped beneath a sub-site, the pore-fluid pressure increases and leads to hydraulic fractures in the basal lenses;

  • once the focused fluid flow is reduced, the diagenetic horizons are shifted downward to their common position below the seafloor and the upper carbonate lens is oxidized, marking the end of a sub-site before its burial;

  • even if a seep site is buried, the conduits, pipe and fractures remain open in unconsolidated and shallow buried sediments, allowing a new phase of deeper fluids to be precipitated. In that case, the mineralogy shows that it is a long-lasting process and there is no evidence that they reached the seafloor. A last phase of bitumen-rich fluids can be recorded in the conduits, meaning that a source rock has reached the oil window at depth and followed the same pathway at the expulsion stage.

The study of the giant fossil pockmark of Beauvoisin shows that several episodes of fluid migration have led to its formation. However, most of the precipitated mass of carbonates was created early, when it was active and cropping out at the seafloor. The other phases remained relict and one (or more) fluid phase may have been missed (not recorded). Some results suggest that the salt tectonic activity may have caused the late two fluid migrations in the area, one coming from the Late Toarcian and/or Bajocian levels during early Cretaceous, one coming from the Triassic evaporites, such as the Propiac diapir, during Mio-Pliocene. It shows that fluid migration in fossil cold seeps can be used as a marker of the regional geodynamic activity of the SE basin of France.

The authors would like to thank all Masters students involved in this project by order: Elie Boidin, Yolaine Rubert, Morgane Bizeray, Alexandra Gueguen, Amandine Castillo, Charly Poitevin, Stanislas Delivet, Maëva Evesque, Laure Matiakh, Aloïse Chabbert-Gondart and Audrey Laplanche. The authors acknowledge all colleagues who participated in discussion in the field and greatly improved our understanding of the area. They also thank Action Marges (INSU) that supported the project for many years since 2009.

Cite this article as: Gay A, Favier A, Potdevin J-L, Lopez M, Bosch D, Tribovillard N, Ventalon S, Cavailhes T, Neumaier M, Revillon S, Travé A, Bruguier O, Delmas D, Nevado C. 2020. Poly-phased fluid flow in the giant fossil pockmark of Beauvoisin, SE basin of France, BSGF - Earth Sciences Bulletin 191: 35.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.