Next Article in Journal
Geochemistry of Carboniferous–Permian Coal from the Wujiawan Mine, Datong Coalfield, Northern China: Modes of Occurrence, Origin of Valuable Trace Elements, and Potential Industrial Utilization
Next Article in Special Issue
Editorial for Special Issue “Chemical, Mineralogical and Isotopic Studies of Diagenesis of Carbonate and Clastic Sediments”
Previous Article in Journal
Seismic Properties of a Unique Olivine-Rich Eclogite in the Western Gneiss Region, Norway
Previous Article in Special Issue
REE Characteristics of Lower Cretaceous Limestone Succession in Gümüşhane, NE Turkey: Implications for Ocean Paleoredox Conditions and Diagenetic Alteration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geochemical and Dynamic Model of Repeated Hydrothermal Injections in Two Mesozoic Successions, Provençal Domain, Maritime Alps, SE-France

1
Faculty of Science, Department of Petroleum Geosciences, Soran University, Soran, Erbil 44008, The Kurdistan Region of Iraq, Iraq
2
Scientific Research Center (SRC), Soran University, Soran, Erbil 44008, The Kurdistan Region of Iraq, Iraq
3
General Directorate of Scientific Research Center, Salahaddin University-Erbil, Erbil 44002, The Kurdistan Region of Iraq, Iraq
4
School of the Environment, University of Windsor, Windsor, ON N9B 3P4, Canada
5
Res. Grp.—Biogeochemistry & Modelling of the Earth System, Université Libre de Bruxelles, 1050 Brussels, Belgium
*
Author to whom correspondence should be addressed.
Minerals 2020, 10(9), 775; https://doi.org/10.3390/min10090775
Submission received: 9 August 2020 / Revised: 26 August 2020 / Accepted: 30 August 2020 / Published: 31 August 2020

Abstract

:
A field, petrographic and geochemical study of two Triassic–Jurassic carbonate successions from the Maritime Alps, SE France, indicates that dolomitization is related to episodic fracturing and the flow of hydrothermal fluids. The mechanism governing hydrothermal fluids has been documented with the best possible spatio-temporal resolutions specifying the migration and trapping of hydrothermal fluids as a function of depth. This is rarely reported in the literature, as it requires a very wide range of disciplines from facies analysis (petrography) to very diverse and advanced chemical methods (elemental analysis, isotope geochemistry, microthermometry). In most cases, our different recognized diagenetic phases were mechanically separated on a centimetric scale and analyzed separately. The wide range of the δ18OVPDB and 87Sr/86Sr values of diagenetic carbonates reflect three main diagenetic realms, including: (1) the formation of replacive dolomites (Type I) in the eogenetic realm, (2) formation of coarse to very coarse crystalline saddle dolomites (Types II and Type III) in the shallow to deep burial mesogenetic realm, respectively, and (3) telogenetic formation of a late calcite cement (C1) in the telogenetic realm due to the uplift incursion of meteoric waters. The Triassic dolomites show a lower 87Sr/86Sr ratio (mean = 0.709125) compared to the Jurassic dolomites (mean = 0.710065). The Jurassic calcite (C1J) shows lower Sr isotopic ratios than the Triassic C1T calcite. These are probably linked to the pulses of the seafloor’s hydrothermal activity and to an increase in the continental riverine input during Late Cretaceous and Early Cenozoic times. This study adds a new insight into the burial diagenetic conditions during multiple hydrothermal flow events.

1. Introduction

The formation of low- versus high-temperature dolomites in carbonate successions is attributed to separate diagenetic conditions and fluid origins [1,2]. Low temperature dolomitizing fluids are associated with near surface diagenetic fluids [3,4,5]. Hot dolomitizing fluids are either the result of high burial temperatures or the flux of hot hydrothermal (HT) fluids into colder, shallower buried carbonate succession environments [6,7,8,9]. The HT fluids are commonly formed by the upward flow of basinal brines [10], or deeply percolating meteoric waters [11] that interact with crystalline basement rocks [12]. In this paper, the petrographical characteristics and geochemical and isotopic compositions of dolomite and calcite in the Triassic-Jurassic successions are used to constrain the diagenetic conditions encountered within the framework of a constructed paragenetic sequence of dolomites in the Provençal Domain in SE France.

2. Regional Geology

The study area is located in the Provençal Domain close to the northeastern part of the Maritime Alps, in the southeast of France, bordered by NW Italy (Figure 1). The area, which is considered as the interior-most part of the European Proximal Margin [13], is underlain by the Argentera Massif from the southern part of the Dauphinois Domain. The Provençal Domain consists of Mesozoic rocks in shallow marine environmental settings [14] and Upper Carboniferous–Permian continental sediments overlain by Triassic rocks consisting of siliciclastics (Lower Triassic), tidal carbonates (Middle Triassic) and evaporites (Upper Triassic) [15].
The European paleo-margin experienced extensional and transtensional tectonics during the opening of the Alpine Tethys during the Jurassic and Cretaceous [16]. A subsequent regional unconformity followed the Triassic succession and encompasses the Late Triassic–Early Jurassic succession [14]. The emersion of the Provençal Domain during Early Jurassic times led to the formation of a carbonate platform followed by a continuous decrease in the sea level and the consequent establishment of a peritidal environment during the Berriasian [17,18]. The deposition of shallow-water marine carbonates continued during the Cretaceous [14]. The Mesozoic succession is separated from the Cenozoic deposits by a regional unconformity that is equivalent to the late Cretaceous-Middle Eocene interval [14]. A rapid lateral facies variation is observed in both Provençal and Dauphinois Domains and the Mesozoic succession is thinner toward the Provençal area [14].

3. Methods and Materials

Conventional and cathodoluminescence (CL) petrography was carried on 50 thin sections (out of total 98 samples) stained with Alizarin Red. The CL examinations were carried out on representative polished thin sections using CL model 8200 at operating conditions of 9–12 kV and a 300–400 μA gun current.
A total of 1.5 mg of the 43 dolomite and calcite samples were dissolved in 3 M HNO3 for the purpose of elemental analysis. The solution was diluted with 2 mL of deionized H2O (N18.2 MΩ cm−1). The element concentrations are reported in ppm. The analytical errors are given as ± %RSD. The 1σ-reproducibility for the major and trace elements of the two standard materials are ±0.18 wt. % for Ca, ±0.081 wt. % for Mg, ±22 ppm for Sr, ±17 ppm for Fe, ±1 ppm for Mn (CRM-512, n = 111), ±0.36 wt. % for Ca, ±0.002 wt. % for Mg, ±1 ppm for Sr, ±12 ppm for Fe and ±1 ppm for Mn (CRM-513, n = 111).
The oxygen and carbon stable isotopic analyses were performed by the reaction of powdered, micro-sampled calcite and dolomite with 100% phosphoric acid at 70 °C using a Gasbench II connected to a ThermoFisher Delta V Plus mass spectrometer. All values are reported in per mil relative to V-PDB. The reproducibility and accuracy were monitored by a replicate analysis of laboratory standards calibrated to international standards NBS19, NBS18 and LSVEC. Isotopic values are referred to the University of Erlangen PDB standard with reproducibilities of ±0.2‰ or better.
The 87Sr/86Sr analyses were performed at the Institute for Geology, Mineralogy and Geophysics, Ruhr-Universität Bochum (Germany). According to the ICP-OES, the Sr concentration measurements need a carbonate sample with 250 to 350 ng of Sr digested with 6 M supra-pure HCl for about 24 h at room temperature in PFA beakers and evaporated until dry. Subsequently, the dried samples were re-dissolved in 0.4 mL HNO3 3M. The Sr fraction was recovered using an ion exchange resin (Sr-resin TRISEM) with 0.05 M and 3 M HNO3 applying 2 mL of distilled water. After evaporation, the dried sample material was treated with 1 mL of a 1:1 mixture of concentrated HNO3:H2O2 removing organic remains and subsequently dried again. Later, the samples were converted to chloride-form by applying 0.4 mL 6 M HCl. The mean of 279 analyses of NIST standard NBS 987 was 0.710242 with a mean standard error of 0.000002 (±2 se) and mean standard deviation of 0.000032 (±2 sd), whereas the mean of 253 analyses of the USGS EN-1 standard was 0.709162 with a mean standard error of 0.000002 (±2 se) and mean standard deviation of 0.000026 (±2 sd).

4. Results

4.1. Field Observation

The studied carbonates display a complex fracture system with carbonate cemented breccias and zebra fabrics (Figure 2A,B and Figure 3). The brecciated fragments show angular edges cemented by coarse-crystalline dolomite. Brecciation is most extensive in the lower part of the Triassic section (Figure 2B). The development of stratabound zebra fabric in the Triassic series is associated with fracture networks with predominantly conjugated orientations extending for decimeters (Figure 2A). Dolomitization is the most common diagenetic alteration in the host carbonates in the Maritime Alps [18].

4.2. Petrography of the Triassic–Jurassic Succession

Petrographic characteristics are used to document the diagenetic fluids that influenced the Mesozoic carbonate successions in the Maritime Alps (Provençal Domain—SE France), their relative timing and burial history. The diagenesis of the carbonate is complex and includes several dolomite generations.

4.2.1. Middle Triassic Section

The micrite and replacive microdolomite (D1T) composition is more common in breccia and zebra fabrics of the Triassic than in the Jurassic carbonates. Despite strong dolomitization, the depositional features are still preserved including laminated algal/microbial mats and slightly compacted micritized ooids with darker micritic envelopes in grain-supported fabrics (Figure 4A,B). Dolomite 2 (D2T) is composed of medium to coarse grained anhedral–subhedral crystals (70–200 µm) in a nonplanar to planar-s fabric. D2T replaced the D1T dolomite and display relict ooids (Figure 4A,B).
Saddle dolomite (SD1T), which occurs as irregular compacted crystals (0.3 to 0.6 mm) with nonplanar-s fabric, replaced D2T. SD1T commonly occurs in the dolostones and, less commonly, as vug-lining cement in fractures. Coarser saddle dolomite SD2T crystals (0.5 up to 4 mm) consist of homogeneous euhedral slightly curved, sometimes planar, crystals. In contrast to SD2T, SD1T appears as a cement in irregular fenestrae and breccias and as zebra dolomite fabric. The brecciated dolostones D2T (clast sizes up to a few cm) are cemented by various dolomitic types. SD1T is distinguished from SD2T crystals by its fabric retentive, inclusion-free rims and inclusion-rich cores. SD1T and SD2T saddle dolomites are characterized by a similar growth zoning with bright red and dark red luminescence. The gray bands of D2T alternate with the white bands of the saddle dolomites (SD1, SD2), forming the zebra fabric and are thus genetically related to the fracturing of the dolostones. SD1T and SD2T are crosscut by low amplitude stylolites (Figure 3C) that were aligned parallel to fracture-filling saddle dolomite.
Light brown colored, coarse-grained, euhedral, slightly planar, Fe-rich saddle dolomite crystals (SD3T) are found in the breccias (Figure 4D–F). The breccias vary from clast-supported (jigsaw puzzle) to cement-supported textures (Figure 2B). SD3T is associated with stylolites with sets of fractures or tension gashes. The tension gashes are parallel to the columns of stylolites and terminated at the wider end of the stylolite seam. SD3T exhibited the typical saddle characteristics with a wavy extinction and curved crystal faces under a cross-polarized light, and has a uniform bright red luminescence under CL.
Euhedral, Fe-poor SD4T saddle dolomite (Figure 5A) displayed alternating inclusion-rich cores and -free rims. The margins of SD4T crystals were corroded and replaced by calcite (C1T; Figure 5B,C). Similar calcite filled the vugs and fractures and replaced D2T, SD3T and SD4T too.

4.2.2. Middle-Upper Jurassic

The peritidal Middle–Upper Jurassic carbonates in the Provençal region include ooidal packstones and grainstones, which are replaced by microcrystalline dolomite (10–300 µm) (D1J) (Figure 6A). SD1J replaced DJ (Figure 6B). A fabric destructive, second saddle dolomite generation (SD2J) with subhedral–euhedral, nonplanar crystals (>1 mm) fill fractures and vugs (Figure 6C).
Yet another phase of coarse-crystalline (several mm across), euhedral, nonplanar to planar, beige colored saddle dolomite (SD3J) (Figure 6D), which fills the fractures or replaces the dolomitic matrix, was recognized. SD3J is postdated by a dark-grey dolomite with different sizes and shapes (D4) of filling space between SD3J crystals (Figure 6D).
Another coarse-crystalline, euhedral, nonplanar, multizoned, red luminescent fracture and vug-filling saddle dolomite (SD5J) (Figure 6E,F) was observed. SD5J crystals are characterized by alternating inclusion-rich and inclusion-free zones. High-amplitude stylolites cross cut the fracture-filling saddle dolomites; no tension gashes were observed. The SD6J dolomite is distinguished from the previous saddle dolomites by its elongated crystal shapes (Figure 7A–C). This SD6J dolomite is densely packed and systematically associated with D4 (Figure 6D). CL revealed an alteration of thick red luminescence (SD5J) and thick nonluminescent SD6J (Figure 8B,C). In addition, thin bright red luminescence is associated with fractures that postdated SD6J.
A Fe-poor calcite cement (C1J) fills the central parts of the fractures and vugs, hence the latest cement (Figure 8A,B). The calcite sizes are ranged from fine to coarse crystals (>400 µm), precipitated and postdated saddle dolomite formation, and had crystals that were transparent with whitish and yellowish colors. CL revealed an alteration of thick red luminescence (SD5J) and thick nonluminescent SD6J (Figure 8B,C). In addition, thin bright red luminescence was associated with small veins and fractures that postdated SD6J.
Dolomite recrystallization (Dr) occurred after all saddle dolomitic phases, where the ghosts of the pristine oolitic facies still resisted diagenesis, and in places the coarse sized poikilotopic cement fabric (>700 μm) developed. Dr displays the same characteristics as those observed in SD6J under CL, i.e., nonluminescent thick with thin bright red luminescent bands (Figure 8E,F).
Because of the small sizes of fractures and cavities, and also the fact that sometimes two dolomite phases occur together in a fracture during mixing the carbonate phases are categorized into four groups to obtain precise data: type 1 with DJ, D1T and D2T; type II with saddle dolomites from SD1 T/J and SD2 T/J; type III with SD3 T/J, SD4T, D4J, SD5T/J, Dr, and SD6J and type IV with C1T/J. This phase succession will also make it easier to interpret the geochemical composition, especially when two carbonate phases are mixed.

4.3. Early and Late Diagenetic Stylolites

Most stylolites are used as a key of the burial depth, especially those aligned parallel to the bedding planes. The Triassic–Jurassic carbonate succession is characterized by two types of stylolites: (i) low amplitude stylolites (no more than 1 cm in length) crosscut Type I dolomites. The low amplitude stylolites (early stylolites) were the first ones produced during the Triassic–Jurassic diagenesis and are associated with a jump-up in McArthur’s graph [19]; (ii) early stylolites are followed by larger stylolites (several centimeters in length with high amplitudes) associated with a set of fractures and tension gashes. Stylolite-related fractures (tension gashes) are rare diagenetic features formed by increasing and progressive mechanical loading during burial, evidenced by the parallel orientation with the bedding plane [20,21].

4.4. Chemical Composition (Major and Trace Elements)

Nineteen samples were selected for trace element analyses from both outcrops (Table 1).
Generally, Ca (mean: 22.1 wt. %; n = 15) and Mg (mean: 12 wt. %; n = 15), Ca (mean: 39.2 wt. %; n = 2) and Mg (mean: 0.6 wt. %; n = 2) are the major elements in dolomite and calcite phases, respectively. Other elements like Sr (mean: 112 ppm; n = 15 for dolomite; mean: 330 ppm; n = 2 for calcite), Fe (mean: 580 ppm; n = 15 for dolomite; mean: 116 ppm; n = 2 for calcite), Mn (mean: 95 ppm; n = 15 for dolomite; mean: 180 ppm; n = 2 for calcite), Ba (mean: 3 ppm; n = 15 for dolomite; mean: 4 ppm; n = 2 for calcite), represent crucial “trace elements” in the carbonate diagenetic setting (Figure 9). The relative abundance of major elements and scarcity of trace elements in dolomite and calcite, are likely the recorded evidence of diagenetic processes (e.g., [22]).
The Ca/Mg ratio of the saddle dolomites (Types II and III) ranges between 1.72 and 2.14 from the Triassic samples, and 1.69 and 2.08 from the Jurassic samples. The Ca/Mg ratio of the late calcite phase in the Triassic is 68 and in the Jurassic is 85 (Figure 9). The Jurassic dolomites show a gradual increase in Fe and Mn, starting from 34 going up to 226 ppm for Mn, and 78 to 730 ppm for Fe. The Triassic dolomites exhibit high Fe and Mn values (mean: 846 and 92 ppm, respectively), with a wide range. Rarely do the Mn contents of the early dolomites in both sections overlap the Mn contents of the later dolomites. The early dolomites display the highest Fe contents. The Fe and Mn contents in the Triassic and Jurassic limestones share the same values with a few late dolomite samples. In contrast, late calcite-Type IV (C1J and C1T) contents exhibit the lowest Fe contents (C1T: Fe = 204 ppm; Mn = 146 ppm) (C1J: Fe = 27 ppm; Mn = 13 ppm). Higher Fe and Mn contents in C1T could be related to SD3T, which was overlain by the calcite cement, i.e., due to mixing of SD3T and C1T during preparation.

4.5. δ18OVPDB and δ13CVPDB Isotopes

The stable isotopic compositions of ten limestone samples (Table 2) shows that the δ18OVPDB and δ13CVPDB values vary between −0.2‰ to −4.2‰ and −0.5‰ to +3.1‰, respectively. The isotopic composition of the four dolomite types (I, II, III and IV) shows a range in δ13CVPDB values from −0.6% to +2.6‰ and δ18OVPDB values from −14.1‰ to −4.9‰. In more detail, the δ18OVPDB values of Type I (D1T/J and D2T/J) display a narrow range between −7.1‰ and −4.8% in the Jurassic and from −7.4‰ to −5.7‰ in the Triassic. The δ13CVPDB values for both sections display similar values (average +1.1‰ and +1.9‰, respectively). All δ18OVPDB values are depleted in comparison with the original isotopic compositions of the marine dolomite signature during Triassic and Jurassic times [22,23,24,25].
The δ18OVPDB values of Type II dolomite show a signature with depleted values varying between −8.3‰ and −7.9% in the Jurassic samples, and from −7.9‰ to −6.2‰ in the Triassic samples. These δ18OVPDB values of saddle dolomites decrease with the increase in crystal size. The δ13C values of Type II and Type III exhibit a narrow range in the Triassic (average +2.0‰ and +2.4‰) and in the Jurassic (average +1.7‰ and +1.9‰). Type IV late-stage calcite shows significantly depleted δ18OVPDB values in the Triassic (−14.1‰ and −12.1 ‰) and in the Jurassic (between −11.8 and −9.6‰).

4.6. 87Sr/86Sr Data

The 87Sr/86Sr data of dolomite and calcite (Table 2) in the Triassic section show a wide range (0.707735 and 0.712023), while in the Jurassic section the range is between 0.708680 and 0.71120, so the ratios are lower than those of the host limestone sample.
The two unaltered limestone samples yielded an 87Sr/86Sr ratio of 0.708766 and 0.707862 for the Triassic and Jurassic limestones, respectively. The Triassic signal displays a value different from the seawater signal when compared to the data of Korte et al. [26], which reported the 87Sr/86Sr values between 0.707562 and 0.70804. Jurassic samples have been documented by Jones et al. [27] with an Sr87/Sr86 ratio in the range of 0.707060–0.707280.

5. Interpretation and Discussion

5.1. Paragenetic Sequence Related to Dolomitizing Fluids

The petrographic and geochemical data obtained suggest that dolomitization has been accomplished by various episodes of hydrothermal flow events under various burial diagenetic conditions. The earliest dolomitization, which affected mainly the micrite mud matrix, probably occurred at near-surface conditions. This is evidenced by the anhedral, non-saddle dolomite fabric [28] and the uncompacted ooids. The nonplanar crystal shape of dolomite Type I D2T could be linked to a diagenesis at a higher temperature (e.g., ≥ 50 °C, [28]). However, these nonplanar dolomites are probably related to impurities in the growth media, which stabilize the crystal faces [29], and not to high temperature fluids. This is supported by petrographic evidence—for example, D2T crystals do not show any characteristics of hot fluid involvement “i.e., lack of curved faces and sweeping extinction”.
It is well reported in the literature that saddle dolomite forms under high temperatures typically by flow of hydrothermal fluids [30]. Type II is dissolved by low-amplitude, early stylolites, in high temperature conditions (senso, [31]) in the shallow–intermediate depth without necessarily implying deep burial conditions “i.e., eogenesis”.
Brecciation and concomitant dolomitization (Type II dolomites) suggest a flow of hydrothermal fluids along fractures (Figure 10). Brecciation was accompanied by the dissolution of the host carbonates, which resulted in the creation of new pathways for fluid of hydrothermal fluids and cementation by producing the Type III saddle dolomite. Further cementation by the latter dolomite was probably related to hydrofracturing and resetting of the zebra and breccia fabrics. Stylolite-related fractures are parallel to the amplitude of stylolite (Figure 4F), and according to Nelson [20] the tension gashes are not true tension fractures, but are extensional fractures derived from the same compressive state of stress as the fractures of stylolites [32]. Filling the tension gash fractures with SD3 confirmed that the late stylolites are coeval with Type III saddle dolomite formation (Figure 11 and Figure 12). Stylolite-related fractures (tension gashes), which were produced during deep burial mesogenesis, were cemented by Type III saddle dolomite (Figure 4F). Hence, Type III saddle dolomite formation was presumably formed by the flux of hot basinal fluids along stylolites during compressional basin tectonics [33].
Subsequent openings of the veins or fractures in Type II facies led the hot secondary brine fluid to percolate along these pathways, which resulted in the precipitation of Type III saddle dolomites. The hydrothermal activity associated with hydrofracturing affecting the studied Mesozoic series has also been reported by [34,35,36].
The homogenization temperature of saddle dolomites from the Dauphinois–Provençal Domains varied between 170 and 260 °C ([14]; Figure 10), which rules out a low temperature origin for dolomite. The final melting temperature of fluid inclusions have been reported as highly saline fluids, and the salinity was approximately equivalent to 20% to 23% CaCl2 [14], this confirms the interaction of hot circulating fluids with Argentera basement. The high temperature dolomitizing fluid is further supported by low δ18OVPDB values. The depletion of the δ18OVPDB values follows the increase in the size of saddle crystal from Type II to Type III (up to −11.6‰ in Type III dolomites). Therefore, the plausible explanation that can adequately explain the range of homogenization temperatures and considerable range of oxygen isotope values is that the precipitation of dolomite phases occurred at different diagenetic settings during which a repeated injection of hydrothermal fluid occurred (Figure 13 and Figure 14). This interpretation is further supported by the following lines of evidence (Figure 4 and Figure 6): (1) suture intergranular contacts between the ooids are engulfed by the saddle dolomite cement; (2) occurrence of sulfides and dissolution of stylolite-related fractures in type III saddle dolomites; (3) co-occurrence of Type III saddle dolomites and late stylolite-related fractures, later the cements inside the fractures were dissolved and etched; (4) coarse-crystalline poikilotopic cement fabric (>700μm) etched and corroded the Dr boundaries; (5) dissolution, compacted saddle dolomite and association with sulfides; (6) the presence of high-amplitude stylolites parallel to the overburden force.
Therefore, according to the abovementioned points, our interpretation agrees with the one of Wierzbicki et al. [37], who suggested that the partial dissolution of dolomite crystal outlines fits the deep burial mesogenesis. This late evolution could be linked to the subsidence and recrystallization processes in the Provençal Domain during Late Cretaceous–Early Eocene times, where the main extensional and transtensional tectonics affected the European paleomargin of the Alpine Tethys [38].
The calcitization of saddle dolomite and the precipitation of the fracture-filling Fe-poor calcite cement (C1T/J, Type IV; Figure 14), which cross cut the high-amplitude stylolites, are the last diagenetic event affecting the Triassic series by the incursion of meteoric waters. Furthermore, the precipitation of bladed calcite cement around saddle crystals implies that a new activity led to the uplift of the studied series and could reactivated the Alpine transcurrent shear zone during Cenozoic [39].

5.2. Geochemical Evolution of Diagenetic Fluids

A crucial question in any study of carbonate chemical composition is whether a primary marine signal is preserved or not. Carbon and oxygen isotopic compositions of the regional limestones span a considerable range with δ13C values between −0.5‰ and +3.1‰ and δ18OVPDB values between −4.2‰ and 0.2‰. Their isotopic compositions are similar to the values expected for Triassic and Jurassic marine carbonates [23,24,25,26,27].
Type I dolomite and the last event of calcite precipitation (Type IV) have considerably lighter δ18OVPDB compared to those expected for the carbonate precipitated in equilibrium with ambient seawater (Figure 13). The δ18O and δ13C values in both types of dolomite and calcite fit the inverse “J” Lohmann curve [40] pointing to a diagenetic alteration near surface conditions. The initial nucleation of dolomite in the muddy facies was therefore driven by dolomite supersaturation caused by the mixing of fresh and marine groundwater [40,41,42]. Mixed groundwater filled the pore spaces in the oolite facies and led to undersaturation of calcite, and the dissolution and precipitation of dolomite in the case of Type I dolomite. The Triassic and Jurassic rocks developed on a shelf in which the water depth was extremely shallow, probably a few meters deep, as suggested by the oolites and algae. In such a shallow depositional setting these carbonates were probably altered and dolomitized by mixed water, as no evaporitic minerals were observed in this study. In the Triassic and Jurassic sections, the δ18OVPDB shifts toward a lower value in comparison with the original depositional facies, are quite abrupt and vary around their respective average (−6.7‰ and −6.4‰) in the range of the “meteoric” [3]. The strong deviation in δ13CVPDB values from this line is related to the rock–water interaction in the Type I Triassic samples. Very weak variations or near constant δ13CVPDB values in the Type I Jurassic samples suggest decreased rock–water interactions in a distal setting from the meteoric recharge with a consequent buffering by the rocks.
Types II and III have depleted δ18O values (up to −10.4‰ for Type III in the Triassic and up to −11.6‰ for Type III in the Jurassic) with more or less constant δ13CVPDB values in the Jurassic (Types II and III) and a wider range of δ13CVPDB values in Type III in the Triassic. The δ13C values are similar to the ones of the early dolomitic facies suggesting that the saddle dolomitizing fluid was mainly buffered by interactions with early dolomite rocks (Type I). The same buffered system has been inferred from our trace element analyses and showed that positive and gradual variation in Fe and Mn contents observed between Types II and III saddle dolomites due to water–rock interactions (Figure 13).
The significant shift to lower δ18OVPDB values from Type II to Type III in both sections is probably due to the interaction of the dolomitized rocks (Type I) with hydrothermal fluids (e.g., [12,43]). Similar carbon and oxygen isotopic values (−9.0‰ to −7.3‰ and +1.8‰ to +2.4‰, respectively) were reported by Nader et al. [42] in the Jurassic rocks at Lebanon as the result of hydrothermal hot fluids under deep conditions. The tension-gashed stylolitization phase during the saddle emplacement between Type II and Type III saddle dolomites suggests that the chemical compaction or burial was coeval with the late phase of hot fluids (Type III dolomite). This phase is associated the most widespread depleted δ18OVPDB values (Figure 14). Triassic Type II facies show some similarities and overlapping δ18OVPDB values, such as those of Type I, suggesting that saddle precipitation started in the same burial conditions of the latter, but after early stylolitization. These complex scenarios of paragenetic models show that Types II and III correspond to complex polyphased dolomitizing fluids and Type III saddle dolomite was precipitated under deep burial–mesogenetic diagenesis. Their geometries (irregular, strata-discordant, zebras, even faults on the field) may suggest that the faults constituted a hydrothermal fluid conduit operating during two major steps with increasingly 18OVPDB depleted values in the saddles as evidenced by two different isotopic domains in both sections. Further evidence for deep burial setting and mesogenesis includes the dissolution of saddle dolomite along stylolites (e.g., [44]) and sutured intergranular contacts between ooids.
Type IV Fe-poor calcite records a dedolomitization and cementation by similar calcite due to the incursion of meteoric waters [40,45]. The wide variations of the δ18OVPDB and δ13CVPDB values reveals an inverse “J” Lohmann curve [40,45], at least in the Jurassic and probably in the Triassic, supporting a telogenetic meteoric origin. The δ18OVPDB and δ13CVPDB values from calcite samples overlap with the values of the HT dolomite; a similar situation has been reported from NW Spain by Shah et al. [46] and Nader et al. [42] in Lebanon on hydrothermal-hosted Jurassic rocks. These authors concluded that the wide range of depleted isotopic oxygen values with a slight shift in the isotopic carbon values of pristine facies indicate different temperatures of the dolomitizing fluids, while the extra-shifted carbon isotope values and highly depleted isotopic oxygen values in the calcite facies indicate surface-derived fluids. The abrupt Fe and Mn drops, Sr rise and weak Mn variations only reported in the Triassic samples that could be related to the SD3 overgrowth by C1T suggest a significant lack in carbonate-buffering system. In addition, the saddle dolomitic contribution to calcite oxygen contents is very low compared to the contribution of the meteoric fluids [46] due to high δ18Owater18Orock ratios. These arguments support that the calcite cement marks the end-member of the diagenetic sequence as a result of a late uplift of the carbonate platform, which initiated this late telogenetic phase (Figure 14).

5.3. Impact of Basement–Fluid Interactions, Continental Riverine and vceanic Fluxes on Strontium Isotope Signatures

Multiple parameters controlling the radiogenic strontium isotope in carbonate rocks, including: (i) the isotopic composition of seawater from which they precipitated and terrigenous (i.e., siliciclastics) sedimentary rocks; (ii) the 87Sr/86Sr ratio is mainly related to the nature of the parent rocks and their continental weathering rates (riverine) and oceanic fluxes of nonradiogenic strontium. However, diagenesis will be the main factor controlling the significant variation of the 87Sr/86Sr ratio in marine carbonates as it will be affected by hydrothermal or cold fluids that control the strontium ratio.
87Sr/86Sr in dolomites (Types I, II, III) vary from 0.707735 and 0.712023 in the Triassic and from 0.708680 up to 0.711203 in the Jurassic (Figure 15). These ratios are associated to the early and late stylolites, prior to the latest deformation during the uplift (Figure 14). The 87Sr/86Sr ratios are higher than those of ambient Triassic and Jurassic seawater [26,27]. A few of these ratios in the Triassic dolomites are lower than those of the marine carbonate obtained from host limestone samples (0.708766; Figure 15), while in the Jurassic dolomites the ratios are higher than those in the marine carbonate samples (0.707862; Figure 15).
When comparing our radiogenic values with the ones of the Triassic 87Sr/86Sr signals of the Lombardic Alps, the strontium values of the Lombardic Alps increased from 0.70805 in the Anisian to 0.708148 in the Ladinian as a result of the various rates of weathering and erosion of young volcanic rocks and old granitic rocks [47], are lower than our values in Provençal Domain (0.708766; Figure 15). Thus, it is more likely that an increase in the riverine flux was the dominant factor causing the 87Sr/86Sr seawater signature rise during the Triassic times in Provençal Domain (cf. [26]), as well as the involved Triaasic carbonates by meteoric waters.
The repeated injection of hydrothermal fluids into the Mesozoic succession in the Provençal Domain and formation of many phases of saddle dolomites might have been expelled by from overpressured lithologies and fluxed through deep-seated faults (Figure 4, Figure 6, Figure 10 and Figure 14; e.g., [30]). The high Sr isotopic ratios of saddle dolomite suggest that the hydrothermal fluids interacted with the crystalline basement and/or siliciclastic deposits [12,30].
The radiogenic enrichment of the Jurassic samples could be deviated from open ocean evaporite, however, the dolomitizing fluids are not in agreement with the involvement of the evaporites which are missing here. If a high strontium ratio represents the involvement of evaporites, then it is expected that a high strontium ratio would correlate with the depletion of oxygen values. Mostly, the dolomitizing fluids that have high strontium ratios would have been enriched with heavier oxygen isotope values, confirming the involvement of evaporites (e.g., [48]). However, this scenario is not in agreement with our Jurassic dolomite samples and the regional geology, therefore the interaction of hot dolomitizing fluid mostly interacted with the Argentera crystalline basement before flowing up to the Jurassic carbonate rocks.
In contrast to Jurassic dolomites, some Triassic dolomites have significantly lower 87Sr/86Sr values with respect to the open marine facies values. The depletion of 87Sr/86Sr ratios started from the Late Triassic up to the Turonian recording pulses of mid-oceanic hydrothermal fluxes that related to the opening and closure of the North Atlantic Ocean [27]. There were several short-term excursions of radiogenic depletion of 87Sr/86Sr ratios highlighting these pulses. The authors of [49] reported excursions of very low 87Sr/86Sr ratios during three short periods (Early Jurassic, Early Cretaceous and Late Cretaceous) that are linked to the rate of hydrothermal activity as a consequence of increased ocean-crust production at the mid-oceanic ridges. Therefore, some of our data in the Triassic could be partly influenced by the short-term excursions of strontium depletion after the deposition of Triassic sediments. In addition, following McArthur’s Graph, our data display that Type II saddle dolomites are associated with a jump-down of the 87Sr/86Sr isotope values. The repetition of the jump-down and jump-up graph of each type of dolomite follows the same trends in the Triassic and Jurassic samples (Figure 15). The graph suggests that the Triassic and Jurassic sections had the same stratigraphic position (not separated) prior to the precipitation of the calcite phase (Type IV).
The ultimate diagenetic phase (Type IV) records the highest enrichment radiogenic composition (0.712023) from the Triassic and the lowest 87Sr/86Sr ratio (0.708680) from the Jurassic samples. This phase attributed to the late stage of calcite cementation and was consistent with the very depleted δ18OVPDB values and the lowest δ13CVPDB values (cf. [50]). The meteoric water has the lowest δ 18OVPDB values and highest 7Sr/86Sr ratios (0.709630 to 0.70948, [51]). The increasing recharge from meteoric water could apply to the late stage of Fe-poor calcite precipitation and is probably related to an uplift of the buried Triassic–Jurassic rocks (Figure 11, Figure 12, Figure 13 and Figure 14). The late calcite cement C1J in the Jurassic section shifted toward decreasing Sr isotope values, and appears to be more affected by other sources besides meteoric water (in contrast to Triassic calcite value). On the basis of the tectonic graph, C1T and C1J are attributed to two different fluids, which were stratigraphically separated during the uplifting of these units (outcrops). Consequently, the lowering of Sr isotope compositions (C1J) in this study, especially in the Jurassic section, could be partly linked to the pulses of the seafloor hydrothermal activity that lowered the 87Sr/86Sr ratios. During the late calcite event, the former dolomitized rocks were strongly brecciated with an in situ accumulation of crushed fragments. This event reveals that an intense tectonic activity was ongoing during the late fluid migration [30,52,53], which could be associated to Late Cretaceous–Eocene times, where general strike-slip fault systems occurred in the European paleomargin [52,53].
The Triassic–Jurassic carbonate seawater in Provençal Domain could partly have been influenced by the input of various rates of weathering and erosion of young volcanic rocks and/or old granitic rocks. These inputs probably persisted after the deposition of Mesozoic succession in Provencal Domain. Finally, fluctuations in the rising and falling of 87Sr/86Sr seawater signatures have already been reported [26,27], therefore, the values of the 87Sr/86Sr fluctuations are the important factors and may have influence the diagenesis products beside the hydrothermal fluids and meteoric waters in Provencal Domain (Figure 15).

6. Conclusions

The main conclusions derived from this integrated field, petrographic and geochemical study of dolomitization in the Mesozoic carbonate succession, Provençal Domain, Maritime Alps, SE France show that:
  • The succession is extensively dolomitized with abundant breccia and zebra textures, and complex fracturing paths due to the influence of repeatedly injections of hydrothermal fluids.
  • Three diagenetic settings are recognized:
    (i)
    An eogenetic realm with weak diagenetic imprints. This setting is confirmed by δ18OVPDB and δ13CVPDB values that fit the inverse “J” Lohmann curve and meteoric line models.
    (ii)
    A second mesogenetic setting in two phases, the first one started with the generation of wispy stylolites (“early dolomite”), and various networks of fractures as well as zebra and breccia fabrics. The fabrics are associated with a first generation of medium- to coarse-grained dominantly euhedral saddle dolomites (Type II) formed through the migration of hot dolomitizing fluids. This occured during early Cretaceous times and was related to the extensional and transtensional tectonics of the European paleomargin of the Alpine Tethys. In the second phase (our Type III dolomite), further physical brecciation and cataclastic fractures cross cut the Type II dolomite under deeper burial conditions with late stylolite-related fractures (tension gashes). The δ18OVPDB values become more depleted with subsequent increases in the saddle crystal sizes. The δ18OVPDB saddle dolomite values (up to −8.3‰ in Type II and up to −11.6‰ in Type III) suggest that at least two pulses of hydrothermal fluids during the precipitation of saddle dolomites occurred. The chemical composition of the successive saddle dolomites was progressively modified from nonferroan to ferroan-rich in an open system and the positive co-variant trends of 87Sr/86Sr vs. δ18OVPDB strongly support an origin from from hot fluids in a deep burial-reducing mesogenetic environment.
    (iii)
    The ultimate diagenesis is related to the late phase of calcite cement precipitation, synchronous with Fe and Mn depletion and extra-negative δ18OVPDB13CVPDB values. A meteoric water under near-surface conditions during telogenesis and uplifting of the Triassic–Jurassic succession is inferred.
  • The depletion and enrichment of radiogenic strontium signals could also be linked to the pulses of mid-oceanic hydrothermal fluxes and an input of riverine fluxes during the opening and closure of the North Atlantic Ocean, beside the involvement of hydrothermal fluids and meteoric waters on Triassic–Jurassic successions.

Author Contributions

Conceptualization, N.S.; methodology, N.S.; software, N.S.; validation, N.S. and A.P.; formal analysis, N.S. and A.P.; investigation, N.S. and A.P.; resources, N.S.; data curation, N.S.; writing—original draft preparation, N.S.; writing—review and editing, H.M. and A.P.; visualization, N.S.; supervision, A.P. and H.M.; project administration, N.S., A.P. and H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors kindly thank Luca Martire, Luca Barale, and Carlo Bertok from Dipartimento di Scienze della Terra, Università di Torino, Italy for collecting samples at the field “Provençal Domain”. We acknowledge also the support of Philippe Muchez from Department of Earth and Environmental Sciences, KUL-university-Belgium for helpful revision of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Machel, H.G. Concepts and models of dolomitization: A critical reappraisal. In The Geometry and Petrogenesis of Dolomite Hydrocarbon Reservoirs; Braithwaite, C.J.R., Rizzi, G., Darke, G., Eds.; Geological Society of London: London, UK, 2004; Volume 235, pp. 7–63. [Google Scholar]
  2. Fu, Q.; Quing, H. Medium and coarsely crystalline dolomites in the Middle Devonian Ratner Formation, southern Saskatchewan, Canada: Origin and pore evolution. Carbonates Evaporites 2011, 26, 111–125. [Google Scholar] [CrossRef]
  3. Land, S. The origin of massive dolomite. J. Geol. Educ. 1985, 33, 87–110. [Google Scholar] [CrossRef]
  4. Machel, G.; Mountjoy, W. Chemistry and environment of dolomitization: A reappraisal. Earth Sci. Rev. 1986, 23, 175–222. [Google Scholar] [CrossRef]
  5. Hardie, A. Dolomitization: A critical view of some current views. J. Sediment. Petrol. 1987, 57, 166–183. [Google Scholar] [CrossRef]
  6. Qing, H.; Mounjoy, W. Formation of coarsely crystalline, hydrothermal dolomite reservoir in the Presqu’ile Barrier, western Canada Sedimentary Basin. Am. Assoc. Pet. Geol. 1994, 78, 55–77. [Google Scholar]
  7. Al-Aasm, I.S. Origin and characterization of hydrothermal dolomite in the Western Canada Sedimentary Basin. J. Geochem. Explor. 2003, 78–79, 9–15. [Google Scholar] [CrossRef]
  8. Smith, B. Origin and reservoir characteristics of Upper Ordovician Trenton-Black River hydrothermal dolomite reservoirs in New York. Bull. Am. Assoc. Pet. Geol. 2006, 90, 1691–1718. [Google Scholar] [CrossRef] [Green Version]
  9. Lavoie, D.; Chi, G. Lower Paleozoic foreland basins in eastern Canada: Tectono-thermal events recorded by faults, fluids and hydrothermal dolomites. Bull. Can. Pet. Geol. 2010, 58, 17–35. [Google Scholar] [CrossRef] [Green Version]
  10. Hardie, A. On the significance of evaporites. Annu. Rev. Earth Planet. Sci. 1991, 19, 131–168. [Google Scholar] [CrossRef]
  11. Iannace, A.; Gasparrini, M.; Gabellone, T.; Mazzoli, S. Late dolomitization in basinal limestones of the southern Apennines fold and thrust belt (Italy). Oil Gas Sci. Technol. Rev. d’IFP Energ. Nouv. 2012, 67, 59–75. [Google Scholar] [CrossRef] [Green Version]
  12. Davies, G.; Smith, L. Structurally controlled hydrothermal dolomite reservoir facies: An overview. Am. Assoc. Pet. Geol. Bull. 2006, 90, 1641–1690. [Google Scholar] [CrossRef]
  13. Mohn, G.; Manatschal, G.; Muntener, O.; Beltrando, M.; Masini, E. Unravelling the interaction between tectonic and sedimentary processes during lithospheric thinning in the Alpine Tethys margins. Int. J. Earth Sci. 2010, 99, 75–101. [Google Scholar] [CrossRef]
  14. Barale, L.; Bertok, C.; Talabani, S.N.; d’Atri, A.; Martire, L.; Piana, F.; Préat, A. Very hot, very shallow hydrothermal dolomitization: An example from the Maritime Alps (north-west Italy–south-east France). Sedimentology 2016, 63, 2037–2065. [Google Scholar] [CrossRef]
  15. Durand, M. Permian to Triassic continental successions in southern Provence (France): An overview. Geol. Soc. Italy 2008, 127, 697–716. [Google Scholar]
  16. De Graciansky, C.; Lemoine, M. Early Cretaceous extensional tectonics in the southwestern French Alps: A consequence of North-Atlantic rifting during Tethyan spreading. Bull. Société Géologique Fr. 1988, 8, 733–737. [Google Scholar] [CrossRef]
  17. Dardeau, G.; Pascal, A. La régression fin-Jurassique dans les Alpes-Maritimes: Stratigraphie, faciès, environnements sédimentaires et influence du bâti structural dans l’arc de Nice. Bull. Bur. Rech. Géologiques Minières 1982, 3, 193–204. [Google Scholar]
  18. Debrand-Passard, S.; Courbouleix, S.; Lienhart, J. Geological synthesis of south-east France. Bur. Geol. Min. Res. Orléans 1984, 125, 615. [Google Scholar]
  19. McArthur, M.; Howarth, J.; Shields, A. Strontium isotope stratigraphy. In The Geologic Time Scale 2012; Gradstein, F.M., Ogg, J.G., Schmitz, M., Ogg, G.M., Eds.; Elsevier: Boston, MA, USA, 2012; Volume 2, pp. 127–144. [Google Scholar]
  20. Nelson, A. Significance of fracture sets associated with stylolite zones. Am. Assoc. Pet. Geol. Bull. 1981, 65, 2417–2425. [Google Scholar]
  21. Eren, M. Origin of stylolite related fractures in Atoka Bank Carbonates, Eddy County, New Mexico, U.S.A. Carbonates Evaporites 2005, 20, 42–49. [Google Scholar] [CrossRef]
  22. Jones, E.; Tarney, J.; Baker, H.; Gerouki, F. Tertiary granitoids of Rhodope, Northern Greece: Magmatism related to extensional collapse of the Hellenic orogen? Tectonophysics 1992, 210, 295–314. [Google Scholar] [CrossRef]
  23. Veizer, J.; Hoefs, J. The nature of O18/O16 and C13/C12 secular trends in sedimentary carbonate rocks. Geochim. Cosmochim. Acta 1976, 40, 1387–1395. [Google Scholar] [CrossRef]
  24. Veizer, J.; Ala, D.; Azmy, K.; Bruckschen, P.; Buhl, D.; Bruhn, F.; Carden, G.A.F.; Diener, A.; Ebneth, S.; Goddéris, Y.; et al. 87Sr/86Sr, δ13C and δ18O evolution of Phanerozoic seawater. Chem. Geol. 1999, 161, 59–88. [Google Scholar] [CrossRef] [Green Version]
  25. Katz, M.E.; Wright, J.D.; Miller, K.G.; Cramer, B.S.; Fennel, K.; Falkowski, P.G. Biological overprint of the geological carbon cycle. Mar. Geol. 2005, 217, 323–338. [Google Scholar] [CrossRef]
  26. Korte, C.; Kozur, W.; Bruckschen, P.; Veizer, J. Strontium isotope evolution of Late Permian and Triassic seawater. Geochim. Cosmochim. Acta 2003, 67, 47–62. [Google Scholar] [CrossRef]
  27. Jones, E.; Jenkyns, C.; Coe, L.; Hesselbo, P. Strontium isotopic variations in Jurassic and Cretaceous seawater. Geochim. Cosmochim. Acta 1994, 58, 3061–3074. [Google Scholar] [CrossRef]
  28. Gregg, J.M.; Sibley, D.F. Epigenetic dolomitization and origin of xenotopic dolomite texture. J. Sediment. Res. 1984, 54, 908–931. [Google Scholar]
  29. Sibley, F.; Gregg, M. Classification of dolomite rock textures. J. Sediment. Res. 1987, 57, 967–975. [Google Scholar]
  30. Salih, N.; Mansurbeg, H.; Kolo, K.; Gerdes, A.; Préat, A. In situ U-Pb dating of hydrothermal diagenesis in tectonically controlled fracturing in the Upper Cretaceous Bekhme Formation, Kurdistan Region-Iraq. Int. Geol. Rev. 2019, 0020–6814, 1938–2839. [Google Scholar] [CrossRef]
  31. Kretz, R. Carousel model for the crystallization of saddle dolomite. J. Sediment. Petrol. 1992, 62, 551–574. [Google Scholar]
  32. Nelson, A. Natural fracture systems: Description and classification. Am. Assoc. Pet. Geol. Bull. 1979, 63, 2214–2232. [Google Scholar]
  33. Morad, S.; Al Suwaidi, M.; Mansurbeg, H.; Morad, D.; Ceriani, A.; Papagoni, M.; Al-Aasm, I. Diagenesis of a limestone reservoir (Lower Cretaceous), Abu Dhabi, United Arab Emirates: Comparison between the anticline crest and flanks. Sediment. Geol. 2019, 380, 127–142. [Google Scholar] [CrossRef]
  34. Carraro, F.; Dal Piaz, V.; Franceschetti, B.; Malaroda, R.; Sturani, C.; Zanella, E. Note Illustrative della Carta Geologica del Massiccio dell’Argentera alla scala 1: 50.000. Mem. Soc. Geol. Italy 1970, 9, 557–663. [Google Scholar]
  35. Malaroda, R. L’Argentera meridionale–Memoria illustrativa della “Geological Map of Southern Argentera Massif (Maritime Alps) 1:25,000”. Sci. Géologiques–Mémoires 1999, 51–52, 241–331. [Google Scholar]
  36. Bigot-Cormier, F.; Poupeau, G.; Sosson, M. Dénudations différentielles du massif cristallin externe alpin de l’Argentera (Sud-Est de la France) révélées par thermochronologie traces de fission (apatites, zircons). Comptes Rendus l’Académie Sci. Paris 2000, 330, 363–370. [Google Scholar] [CrossRef]
  37. Wierzbicki, R.; Dravis, J.; Al-Aasm, I.; Harland, N. Burial dolomitization and dissolution of Upper Jurassic Abenaki platform carbonates, Deep Panuke reservoir, Nova Scotia, Canada. Am. Assoc. Pet. Geol. 2006, 90, 1843–1861. [Google Scholar] [CrossRef]
  38. Crampton, L.; Allen, A. Recognition of forebulge unconformities associated with early-stage Foreland Basin development—Example from the North Alpine Foreland Basin. Am. Assoc. Pet. Geol. Bull. 1995, 79, 1495–1514. [Google Scholar]
  39. D’Arti, A.; Piana, F.; Barale, L.; Bertok, C.; Martire, L. Geological setting of the southern termination of Western Alps. Int. J. Earth Sci. 2016, 105, 1831–1858. [Google Scholar]
  40. Lohmann, C. Geochemical patterns of meteoric diagenetic systems and their application to studies of paleokarst. In Paleokarst; James, N.P., Choquette, P.W., Eds.; Springer: New York, NY, USA, 1988; pp. 58–80. [Google Scholar]
  41. Allan, R.; Mathhews, K. Carbon and oxygen isotopes as diagenetic and stratigraphic tools: Data from surface and subsurface Barbados, West Indies. Geology 1977, 5, 16–20. [Google Scholar] [CrossRef]
  42. Nader, H.; Swennen, R.; Ottenburgs, R. Karst-meteoric dedolomitisation in Jurassic carbonates, Lebanon. Geol. Belg. 2003, 6, 3–23. [Google Scholar]
  43. Moore, H. Carbonate Diagenesis and Porosity. Dev. Sedimentol. 1989, 46, 338. [Google Scholar]
  44. Wright, P.; Harris, P. Carbonate Dissolution and Porosity Development in the Burial (mesogenetic) Environment. In Proceedings of the Oral presentation at American Association of Petroleum Geologists, Annual Convention and Exhibition, Pittsburgh, PA, USA, 19–22 May 2013. [Google Scholar]
  45. Salih, N.; Kamal, M.; Préat, A. Classical observations and stable isotopes for identification the diagenesis of Jeribe formation at Jambour oil Fields-Kurdistan Region-Iraq. Pet. Sci. Technol. 2019, 37, 1548–1556. [Google Scholar] [CrossRef]
  46. Shah, M.; Nader, H.; Garcia, D.; Swennen, E.R. Hydrothermal dolomites in the Early Albian (Cretaceous) platform carbonates (NW Spain): Nature and origin of dolomites and dolomitizing fluids. Oil Gas Sci. Technol. 2012, 67, 97–122. [Google Scholar] [CrossRef] [Green Version]
  47. Faure, G.; Assereto, R.; Tremba, L. Strontium isotope composition of marine carbonates of Middle Triassic to Early Jur assic age, Lombardic Alps, Italy. Sedimentology 1978, 25, 523–543. [Google Scholar] [CrossRef]
  48. Li, S.; Levin, E.; Chesson, A. Continental scale variation in 17O-excess of meteoric waters in the United States. Geochim. Cosmochim. Acta 2015, 164, 110–126. [Google Scholar] [CrossRef]
  49. Jones, E.; Jenkyns, C. Seawater strontium isotopes, oceanic anoxic events, and seafloor hydrothermal activity in the Jurassic and Cretaceous. Am. J. Sci. 2001, 301, 112–149. [Google Scholar] [CrossRef] [Green Version]
  50. Salih, N.; Mansurbeg, H.; Kolo, K.; Préat, A. Hydrothermal Carbonate Mineralization, Calcretization, and Microbial Diagenesis Associated with Multiple Sedimentary Phases in the Upper Cretaceous Bekhme Formation, Kurdistan Region-Iraq. Geosciences 2019, 9, 459. [Google Scholar] [CrossRef] [Green Version]
  51. Dong, S.; Chen, D.; Qing, H.; Zhou, X.; Wang, D.; Guo, Z.; Jiang, M.; Qian, Y. Hydrothermal alteration of dolostones in the Lower Ordovician, Tarim Basin, NW China: Multiple constraints from petrology, isotope geochemistry and fluid inclusion microthermometry. Mar. Pet. Geol. 2013, 46, 270–286. [Google Scholar] [CrossRef]
  52. Sanchez, G.; Rolland, Y.; Jolivet, M.; Brichau, S.; Corsini, M.; Carter, A. Exhumation controlled by transcurrent tectonics: The Argentera-Mercantour massif (SW Alps): Exhumation controlled by transcurrent tectonics in SW Alps. Terra Nova 2011, 23, 116–126. [Google Scholar] [CrossRef]
  53. Dumont, T.; Schwartz, S.; Guillot, S.; Simon-Labric, T.; Tricart, P.; Jourdan, S. Structural and sedimentary records of the Oligocene revolution in the Western Alpine arc. J. Geodyn. 2012, 56–57, 18–38. [Google Scholar] [CrossRef] [Green Version]
Figure 1. General geological map of the south-west Alps. The red rectangle shows a part of the study area (Provençal Domain) (after [14]).
Figure 1. General geological map of the south-west Alps. The red rectangle shows a part of the study area (Provençal Domain) (after [14]).
Minerals 10 00775 g001
Figure 2. Triassic outcrops characterized by fractures arranged in horizontal and sub-vertical bands to the bedding planes, resulting in zebra-like structures (A), open fractures zones, and breccias (B).
Figure 2. Triassic outcrops characterized by fractures arranged in horizontal and sub-vertical bands to the bedding planes, resulting in zebra-like structures (A), open fractures zones, and breccias (B).
Minerals 10 00775 g002
Figure 3. Basic stratigraphic litho-log of the Triassic and Jurassic carbonate rocks in Provençal Domain. The lower most part of Triassic rock (1) is extensively brecciated and the floated fragments have been cemented by saddle dolomite (see the enlarged part of this breccia on the right part “1”). The breccia is followed by another structure, i.e., zebra dolomites.
Figure 3. Basic stratigraphic litho-log of the Triassic and Jurassic carbonate rocks in Provençal Domain. The lower most part of Triassic rock (1) is extensively brecciated and the floated fragments have been cemented by saddle dolomite (see the enlarged part of this breccia on the right part “1”). The breccia is followed by another structure, i.e., zebra dolomites.
Minerals 10 00775 g003
Figure 4. (A) The depositional texture is grain-supported with ooids (grainstones). The intergranular contacts between ooids are sutured and convexo–concave due to overburden load (arrows). (B) Transition from a mostly anhedral dolomite D2T to a euhedral saddle dolomite SD2T separated by an irregular stylolite indicating that these dolomites were formed at different time. (C) Early stylolite (arrows) cross-cuts SD1T with precipitation of sulfide minerals. (D) Light brown saddle dolomite SD3T, showing coarse crystal size and pressure dissolution. (E) Macroscopic fracture showing Fe-rich SD3T (red arrow) and Fe-poor SD2T (white arrow). A late-stage Fe-poor calcite cement appears at the top left (staining). (F) Arrows point to sets of tension gashes filled by Fe-rich SD3T fracture. The fractures and columns of stylolite (red arrows) are in the same direction and this could be related to the mechanical overburden in the deep burial diagenetic setting.
Figure 4. (A) The depositional texture is grain-supported with ooids (grainstones). The intergranular contacts between ooids are sutured and convexo–concave due to overburden load (arrows). (B) Transition from a mostly anhedral dolomite D2T to a euhedral saddle dolomite SD2T separated by an irregular stylolite indicating that these dolomites were formed at different time. (C) Early stylolite (arrows) cross-cuts SD1T with precipitation of sulfide minerals. (D) Light brown saddle dolomite SD3T, showing coarse crystal size and pressure dissolution. (E) Macroscopic fracture showing Fe-rich SD3T (red arrow) and Fe-poor SD2T (white arrow). A late-stage Fe-poor calcite cement appears at the top left (staining). (F) Arrows point to sets of tension gashes filled by Fe-rich SD3T fracture. The fractures and columns of stylolite (red arrows) are in the same direction and this could be related to the mechanical overburden in the deep burial diagenetic setting.
Minerals 10 00775 g004
Figure 5. (A) The former saddle dolomite SD3T leaving traces along the core and cortex of zoned SD4T, reveals that SD4T was formed under different diagenetic conditions. C1T filled the remaining secondary pores, and postdating the latter saddle dolomite. (B,C) Millimeter-sized fracture-filling Fe-poor, late stage calcite cement (C1T), postdating late stylolite (arrows).
Figure 5. (A) The former saddle dolomite SD3T leaving traces along the core and cortex of zoned SD4T, reveals that SD4T was formed under different diagenetic conditions. C1T filled the remaining secondary pores, and postdating the latter saddle dolomite. (B,C) Millimeter-sized fracture-filling Fe-poor, late stage calcite cement (C1T), postdating late stylolite (arrows).
Minerals 10 00775 g005
Figure 6. (A) Mostly planar-s dolomite D2J with remaining of a former replacement dolomite DJ “lithorelic” (arrows) showing dissolution and irregular compacted grain boundaries. Micrite of precursor facies on the top left corner. (B) Saddle dolomite-type SD1J matrix. The core is rich in inclusions. Thus, this dolomite facies is a key for the replacement of a former dolomite. (C) Secondary porosity filled by saddle dolomite SD2J cement. (D) Secondary porosity occluded by late dolomite D4 of different grain sizes, within SD3J. (E,F) Saddle dolomite SD5J is characterized by strongly curved crystal boundaries and a sweeping extinction under crossed nicols. D, Dolomite; SD, Saddle dolomite.
Figure 6. (A) Mostly planar-s dolomite D2J with remaining of a former replacement dolomite DJ “lithorelic” (arrows) showing dissolution and irregular compacted grain boundaries. Micrite of precursor facies on the top left corner. (B) Saddle dolomite-type SD1J matrix. The core is rich in inclusions. Thus, this dolomite facies is a key for the replacement of a former dolomite. (C) Secondary porosity filled by saddle dolomite SD2J cement. (D) Secondary porosity occluded by late dolomite D4 of different grain sizes, within SD3J. (E,F) Saddle dolomite SD5J is characterized by strongly curved crystal boundaries and a sweeping extinction under crossed nicols. D, Dolomite; SD, Saddle dolomite.
Minerals 10 00775 g006
Figure 7. (AC) More than 1mm sized-dolomite crystals SD6J with two sets of cleavages, precipitated in fractures in contact with recrystallized dolomite. SD6J shows etching and overgrowth of coarse dolomite. (D,E) Saddle dolomite rhomb with an irregular, selectively dissolved outline which predates the poikilotopic texture within the recrystallized dolomite matrix. Two sets of curvature twin laminae, and dissolution of high amplitude stylolite intersecting coarse recrystallized dolomite crystals (the upper side part of Dr and lower side part of E). (F) Bladed calcite cement (arrows) precipitated on saddle dolomite crystal: note the trace of calcite across saddle dolomite; such cement strongly refers to late meteoric alteration at near-surface condition. SD, Saddle dolomite; Dr, Dolomite recrystallization.
Figure 7. (AC) More than 1mm sized-dolomite crystals SD6J with two sets of cleavages, precipitated in fractures in contact with recrystallized dolomite. SD6J shows etching and overgrowth of coarse dolomite. (D,E) Saddle dolomite rhomb with an irregular, selectively dissolved outline which predates the poikilotopic texture within the recrystallized dolomite matrix. Two sets of curvature twin laminae, and dissolution of high amplitude stylolite intersecting coarse recrystallized dolomite crystals (the upper side part of Dr and lower side part of E). (F) Bladed calcite cement (arrows) precipitated on saddle dolomite crystal: note the trace of calcite across saddle dolomite; such cement strongly refers to late meteoric alteration at near-surface condition. SD, Saddle dolomite; Dr, Dolomite recrystallization.
Minerals 10 00775 g007
Figure 8. Photomicrograph of (A) saddle dolomite (SD2) and calcite (C1) in PPL, calcite overgrowth on the saddle dolomite is clearly highlighted under cathodoluminescence in (B). (C) Saddle dolomite (SD6) in PPL and (D) cathodoluminescence (CL) photomicrograph, illustrating the general pattern in luminescence, including thicker zones of red and nonluminescence with shades of thin bright red luminescence, reflecting formation in pore waters that chemically changed. (E) Recrystallized dolomite (Dr) in PPL, (F) the same as in (E) CL photomicrograph, the recrystallized dolomite shows thick nonluminescent zones with thin bright red luminescent bands.
Figure 8. Photomicrograph of (A) saddle dolomite (SD2) and calcite (C1) in PPL, calcite overgrowth on the saddle dolomite is clearly highlighted under cathodoluminescence in (B). (C) Saddle dolomite (SD6) in PPL and (D) cathodoluminescence (CL) photomicrograph, illustrating the general pattern in luminescence, including thicker zones of red and nonluminescence with shades of thin bright red luminescence, reflecting formation in pore waters that chemically changed. (E) Recrystallized dolomite (Dr) in PPL, (F) the same as in (E) CL photomicrograph, the recrystallized dolomite shows thick nonluminescent zones with thin bright red luminescent bands.
Minerals 10 00775 g008
Figure 9. Diagram of the main chemical elements involved in the diagenesis and of the host limestones from the study area (n = 20). For more detail about Type I, II, III, IV, see text. (a) Ca/Mg vs. Fe for Jurassic carbonate samples, (b) Ca/Mg vs. Fe for Triassic carbonate samples, (c) Mn vs. Fe for Jurassic carbonate samples, (d) Mn vs. Fe for Triassic carbonate samples, (e) Ca/Mg vs. Sr for Jurassic carbonate samples, (f) Ca/Mg vs. Sr for Triassic carbonate samples, (g) Ca/Mg vs. Mn for Jurassic carbonate samples, and (h) Ca/Mg vs. Mn for Triassic carbonate samples.
Figure 9. Diagram of the main chemical elements involved in the diagenesis and of the host limestones from the study area (n = 20). For more detail about Type I, II, III, IV, see text. (a) Ca/Mg vs. Fe for Jurassic carbonate samples, (b) Ca/Mg vs. Fe for Triassic carbonate samples, (c) Mn vs. Fe for Jurassic carbonate samples, (d) Mn vs. Fe for Triassic carbonate samples, (e) Ca/Mg vs. Sr for Jurassic carbonate samples, (f) Ca/Mg vs. Sr for Triassic carbonate samples, (g) Ca/Mg vs. Mn for Jurassic carbonate samples, and (h) Ca/Mg vs. Mn for Triassic carbonate samples.
Minerals 10 00775 g009
Figure 10. Histogram of homogenization temperatures obtained for saddle dolomite and fibrous dolomite cements [14].
Figure 10. Histogram of homogenization temperatures obtained for saddle dolomite and fibrous dolomite cements [14].
Minerals 10 00775 g010
Figure 11. Paragenetic succession illustrating the diagenetic phases of Triassic section. Provençal Domain, SE France (see text for explanation).
Figure 11. Paragenetic succession illustrating the diagenetic phases of Triassic section. Provençal Domain, SE France (see text for explanation).
Minerals 10 00775 g011
Figure 12. Paragenetic succession illustrating the diagenetic phases of Jurassic section. Provençal Domain, SE France (see text for explanation).
Figure 12. Paragenetic succession illustrating the diagenetic phases of Jurassic section. Provençal Domain, SE France (see text for explanation).
Minerals 10 00775 g012
Figure 13. Cross-plot of oxygen and carbon isotopic compositions (VPDB) of selected dolomite, calcite phases (Type I to IV: see text) and pristine limestone from Triassic–Jurassic succession (n = 80).
Figure 13. Cross-plot of oxygen and carbon isotopic compositions (VPDB) of selected dolomite, calcite phases (Type I to IV: see text) and pristine limestone from Triassic–Jurassic succession (n = 80).
Minerals 10 00775 g013
Figure 14. Petro-physical and geochemical model tracing the Triassic–Jurassic succession in Provençal Domain, SE France. The detailed microfacies in each diagenetic setting is described by digitized microphotographs under optical microscope. The initial precipitation of saddle dolomite is started in a shallow burial system (see the key evidence in digitized microphotographs). The mechanisms of hydrothermal (HT) fluid movement under deep burial conditions (mesogenesis) created another active pathway for hot fluid with different phase of saddle dolomite (see the zoned saddle dolomites, suture contact boundaries between the ooid grains, and stylolite-related fractures) from that of shallow burial conditions (see the low amplitude stylolite). The diagenesis in the Triassic–Jurassic succession ended by precipitation of Fe-poor calcite during uplifting (close up the large fracture that postdates all saddle phases on the right side of digitized microphotograph, while the one on the left side showing the digitized thin section under CL which result in a band of luminescence due to a change in fluid chemistry.
Figure 14. Petro-physical and geochemical model tracing the Triassic–Jurassic succession in Provençal Domain, SE France. The detailed microfacies in each diagenetic setting is described by digitized microphotographs under optical microscope. The initial precipitation of saddle dolomite is started in a shallow burial system (see the key evidence in digitized microphotographs). The mechanisms of hydrothermal (HT) fluid movement under deep burial conditions (mesogenesis) created another active pathway for hot fluid with different phase of saddle dolomite (see the zoned saddle dolomites, suture contact boundaries between the ooid grains, and stylolite-related fractures) from that of shallow burial conditions (see the low amplitude stylolite). The diagenesis in the Triassic–Jurassic succession ended by precipitation of Fe-poor calcite during uplifting (close up the large fracture that postdates all saddle phases on the right side of digitized microphotograph, while the one on the left side showing the digitized thin section under CL which result in a band of luminescence due to a change in fluid chemistry.
Minerals 10 00775 g014
Figure 15. The rate of change of 87Sr/86Sr according to time in order to illustrate the number of jumps between the different diagenetic conditions in both sections separately. The number of jumps (down and up) can calculate the number of fault reactivation (n = 20).
Figure 15. The rate of change of 87Sr/86Sr according to time in order to illustrate the number of jumps between the different diagenetic conditions in both sections separately. The number of jumps (down and up) can calculate the number of fault reactivation (n = 20).
Minerals 10 00775 g015
Table 1. Elemental analysis of diagenetic phases from Triassic and Jurassic samples (n = 19). Type I = D1T/J and D2T/J. Type II = SD1 T/J and SD2 T/J. Type III = SD3 T/J, SD4T, D4J, SD5T/J, DrJ, and SD6J. Type IV = CT1T/J. T = Triassic, J = Jurassic, and Lst = Host limestone. SE France sections.
Table 1. Elemental analysis of diagenetic phases from Triassic and Jurassic samples (n = 19). Type I = D1T/J and D2T/J. Type II = SD1 T/J and SD2 T/J. Type III = SD3 T/J, SD4T, D4J, SD5T/J, DrJ, and SD6J. Type IV = CT1T/J. T = Triassic, J = Jurassic, and Lst = Host limestone. SE France sections.
Sample No.TypeMn
(ppm)
Fe
(ppm)
Ba
(ppm)
Sr
(ppm)
Mg
(wt. %)
Ca
(wt. %)
Ca (wt. %)/
Mg (wt. %)
NM4Type III147433312512.721.81.7
NM8Type IV132742910.539.486.0
NM10aType II3415236012.821.91.7
NM12Type II3778619711.423.72.1
NM16Type II108317413412.122.71.9
NM23Type I110199127312.421.51.7
NM23Type III661076210611.922.71.9
NM28Type I10771848312.721.81.7
NM30Type IV14620433680.638.968.0
NM34Type III91282.328611.523.42.0
NM4Type III14044136112.220.61.7
NM21Type III226730418011.921.81.8
NM22Type III40168310312.621.81.7
NM11Type I5146328311.422.11.9
NM34Type III6129117911.421.91.9
NM32Type II6033724712.321.21.7
NM23Type III1491225825710.823.12.1
NM44Lst-Triassic113130927912.220.41.7
NM39Lst-Jurassic1015431270.337.8135.5
Table 2. δ13C (% VPDB), δ18O (% VPDB), (Sr87/Sr86) values of selected samples from Triassic–Jurassic succession (n = 100). Type I = D1T/J and D2T/J. Type II = SD1 T/J and SD2 T/J. Type III = SD3 T/J, SD4T, D4J, SD5T/J, DrJ, and SD6J. Type IV = CT1T/J. D is for dolomite, SD = saddle dolomite, C= calcite, Dr = recrystallized dolomite, T = Triassic, J = Jurassic, and Lst = Host limestone. SE France sections.
Table 2. δ13C (% VPDB), δ18O (% VPDB), (Sr87/Sr86) values of selected samples from Triassic–Jurassic succession (n = 100). Type I = D1T/J and D2T/J. Type II = SD1 T/J and SD2 T/J. Type III = SD3 T/J, SD4T, D4J, SD5T/J, DrJ, and SD6J. Type IV = CT1T/J. D is for dolomite, SD = saddle dolomite, C= calcite, Dr = recrystallized dolomite, T = Triassic, J = Jurassic, and Lst = Host limestone. SE France sections.
AgeSample No.Typeδ13Cδ18OSr87/Sr86
(‰ VPDB)(‰ VPDB)
Upper JurassicNM1Type III1.72−9.09
NM2Type II1.78−8.25
NM3Type III1.46−9.35
NM3Type III1.78−9.94
NM3Type II2.06−8.03
NM4Type I1.81−6.44
NM5Type I1.74−7.07
NM5Type II1.76−7.85
NM5Type III1.83−9.99
NM6Type I2.04−4.87
NM7Type III1.69−9.82
NM7Type III1.72−9.43
NM7Type III1.72−9.83
NM8Type III1.81−8.75
NM8Type III1.63−9.61
NM8Type IV0.38−11.830.70868
NM9Type III1.81−10.86
NM10aType II1.55−8.020.71030
NM10bType II1.44−8.180.710552
NM11Type III1.48−10.78
NM11Type III2.07−10.56
NM12Type III1.78−9.14
NM13Type III1.76−9.28
Upper JurassicNM14Type IV−0.59−11.12
NM14Type IV−0.02−9.980.710355
NM15Type III1.85−9.42
NM15Type III1.86−9.710.710154
NM16Type I1.76−6.310.708757
NM16Type III1.42−9.65
NM17Type III1.79−9.31
NM17Type III1.8−9.47
NM18Type IV0.92−10.42
NM18Type III1.59−10.26
NM18Type III1.57−9.66
NM19Type III1.57−9.720.711203
NM19Type III1.48−9.23
NM20Type IV0.98−9.62
NM20Type IV−0.15−10.52
NM20Type III1.29−9.31
NM21Type III1.61−11.630.710609
NM21Type I1.98−6.82
NM22Type III1.82−9.050.708879
NM22Type III1.84−9.45
NM22Type III1.83−9.37
TriassicNM23Type II2.3−6.2
NM23Type II2.3−6.8
NM24Type III2.3−9.98
NM24Type III1.97−9.09
NM25Type III2.56−9.86
NM27Type II2.31−7.58
NM27Type III2.21−8.280.708299
NM27Type II2.21−7.48
NM26Type I1.98−5.7
NM26Type I1.31−6.440.707735
NM28Type III2.54−8.51
NM28Type III1.46−9.35
NM29Type II2.2−6.670.708671
NM30Type IV0.65−14.140.712023
NM31Type III2.11−8.44
NM32Type III2.2−8.57
NM32Type III2.02−9.91
NM33Type III2.49−8.74
NM33Type III2.03−9.440.70942
NM34Type III1.81−9.63
NM34Type II1.97−7.93
NM34Type III2−9.160.708184
NM35Type III2.26−9.07
NM35Type IV2.09−12.08
NM36Type I0.21−7.350.708595
NM36Type I0.93−7.350.711628
NM37Type III1.16−10.350.710469
Upper JurassicNM38Limestone2.07−2.33
NM39Limestone−0.53−4.210.707862
TriassicNM40Limestone1.54−2.31
NM41Limestone1.73−2.5
NM42Limestone1.88−3.36
NM43Limestone1.98−1.67
NM44Limestone2.53−2.370.708766
NM45Limestone2.31−3.36
NM46Limestone3.05−0.23
NM47Limestone2.71−3.25

Share and Cite

MDPI and ACS Style

Salih, N.; Mansurbeg, H.; Préat, A. Geochemical and Dynamic Model of Repeated Hydrothermal Injections in Two Mesozoic Successions, Provençal Domain, Maritime Alps, SE-France. Minerals 2020, 10, 775. https://doi.org/10.3390/min10090775

AMA Style

Salih N, Mansurbeg H, Préat A. Geochemical and Dynamic Model of Repeated Hydrothermal Injections in Two Mesozoic Successions, Provençal Domain, Maritime Alps, SE-France. Minerals. 2020; 10(9):775. https://doi.org/10.3390/min10090775

Chicago/Turabian Style

Salih, Namam, Howri Mansurbeg, and Alain Préat. 2020. "Geochemical and Dynamic Model of Repeated Hydrothermal Injections in Two Mesozoic Successions, Provençal Domain, Maritime Alps, SE-France" Minerals 10, no. 9: 775. https://doi.org/10.3390/min10090775

APA Style

Salih, N., Mansurbeg, H., & Préat, A. (2020). Geochemical and Dynamic Model of Repeated Hydrothermal Injections in Two Mesozoic Successions, Provençal Domain, Maritime Alps, SE-France. Minerals, 10(9), 775. https://doi.org/10.3390/min10090775

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop