Abridged English version
1 Introduction
During the last two decades, geochemistry has turned towards chemostratigraphic studies [4,8,19]. By facing geochemical data of trace elements and stable isotopes (δ18O, δ13C) with sedimentological data, it is possible to propose sequence stratigraphic interpretations. These chemostratigraphic studies are well developed in basin settings where diagenesis is slightly influent, because of a homogeneous stable LMC (Low Magnesian Calcite) mineralogy. Such chemostratigraphic approaches in platform settings are recent [12,16,21,22,24,25], but more complex and discussed because of: (1) an intermittent and diversified sedimentary record, (2) numerous initial mineralogy (aragonite, LMC, HMC or High Magnesian Calcite), and (3) a highly diversified diagenesis. The aim of this work is to understand the meaning of the isotopic signals of bulk neritic carbonates. Through a confrontation with different data sets (sedimentological data, microdrilling isotopic data, petrophysical data), the main scope is to distinguish diagenetic and original parts within the bulk signal.
2 Geological setting
This study is performed on the HTM102 borehole, which crosses the Upper-Jurassic carbonates of Lorraine (eastern edge of the Paris Basin, Fig. 1) from Middle Oxfordian to Lower Kimmeridgian (Fig. 2; [9]).
2.1 Depositional environments
Facies heterogeneity along the core (Fig. 2; legend of Fig. 3) is interpreted as the vertical and lateral evolution of five major depositional environments from subtidal to supratidal [23]. The environmental variations illustrate the stages of evolution of the Lorraine platform that begin with reefal installation above the Callovo-Oxfordian clays in Early Middle Oxfordian, and end at Lower/Upper Kimmeridgian with the return of a terrigeneous sedimentation [6,7,10,11].
2.2 Sequence stratigraphy
Six middle-frequency stratigraphic cycles (Lo1 to Lo6, each of about 1 Ma, 3rd order; Fig. 2) are arranged into a lower frequency cycle (2nd order; Fig. 2) bounded by two maximum of deepening located in the underlying Callovo-Oxfordian clays [3,18] and in the overlying Exogyra Nana Marls [13].
3 Isotopic signal: acquisition and results
3.1 Bulk rocks
A small amount of each sample powder (∼0.5 g) was treated by phosphoric acid attack to extract the CO2. Measurements were performed with a mass spectrometer within the ‘Laboratoire de géologie des bassins sédimentaires’ of the Paris-6 University. All data are reported to the PDB, and the standard deviations proposed for δ18O and δ13C are respectively 0.08‰ and 0.04‰.
The curves are constructed by analysing 100 samples along the core. All the breaks in the curves correspond to sedimentological boundaries (Fig. 2). The less negative δ18O are associated with micritic levels (mudstone to packstone), whereas the more negative values correspond to sparitic levels (grainstone). This is confirmed by a simple bivariate scatter of the isotopic data according to their exact texture (Fig. 4). Fig. 4, however, shows that there is no particular arrangement with the δ13C.
3.2 Separated products
Sampling is performed with a microdrill and concerns: (1) initial LMC of shells insensitive to diagenesis (brachiopods, oysters) that reflect the initial marine isotopic composition, and (2) blocky LMC cements that reflect the intraformational fluid compositions during burial. Measurement procedure and characteristics are identical to those used for bulk rocks.
The mean δ13C of the two sets of products are almost equal, whereas the mean δ18O are different with blocky calcite cement values more negative than initial LMC values (Fig. 5). The negative bulk δ18O of granular carbonates are due to the very negative δ18O of the cements (as already shown by [15]), which here point to parent fluids with mixed marine/meteoric compositions [23], whereas the δ18O of the muddy carbonates come close to the initial marine values but remain more negative (contrary to [20]).
4 Discussion
4.1 Relationship between bulk isotopic signal and texture
The statistical range of δ18O and δ13C for each texture is proposed in Figs. 6 and 7.
An ANOVA confirms that there is a highly significant effect of the texture of the carbonates on their δ18O, and a second ANOVA shows that this is due to the presence/absence of blocky cements (Fig. 6). On the other hand, another set of two ANOVA shows that the significant effect of the texture of carbonates on their δ13C is not due to the presence/absence of blocky cement (Fig. 7). This effect of texture on the δ13C is only ascribed to the differentiation of the wackestone values (Fig. 7), which all correspond to either (1) exposure levels, or (2) encrusted bioconstructions (Fig. 2). In the first case, the occurrence of pedogenetic features can explain the decrease of the δ13C [1]. In the second case, the weak δ13C values may illustrate the influence of a CO2 enriched in 12C generated by organic activity.
4.2 Influence of diagenesis: alteration of micrites
The bulk δ18O variations are not only governed by the textures of the carbonates, as shown by the similarity between the δ18O curve and the evolution of porosity along the borehole (Fig. 2). This similarity can be explained by the study of a specific interval (220–285 m; Fig. 2), which displays an alternation of porous micrite-dominant levels and less porous sparite-dominant levels. According to the relationship between δ18O and texture, the δ18O curve should be irregular along this interval. However, this is not the case and the curve is smoothed.
This is due to the diagenetic history of the porous Oxfordian micrites that have undergone telogenetic dissolutions and reprecipitations as shown by particular nanostructural indications (relic grains, micritic bridges; Fig. 8). The δ18O data obtained after microsampling of these micrites clearly point to a diagenetic alteration, whereas the dense micrites data already reflect marine compositions (Fig. 9). It is possible that the products of reprecipitation between grains, observed within porous micrites, can keep a negative isotopic record of the meteoric fluids that have flowed through the carbonate formations of Lorraine since the Tertiary. Some of these fluids lead to large dissolutions [23]. All these observations go against the widespread idea that considers porous micrites as initial, and dense micrites as the most modified [14,17].
Unlike the δ18O, the bulk δ13C signal is not so modified by these diagenetic alterations as its curve really differs from the porosity evolution (Fig. 2). This is expected since the amount of new C coming with the diagenetic fluids is low compared to the amount of C in the sedimentary rocks. Thereby, the diagenetic fluids need less rock–water interactions to be in equilibrium with surrounding carbonates according to the δ13C rather than according to the δ18O [2].
5 Conclusions
Even if the breaks in the δ18O curve correspond to sequential boundaries (maximum of regression and deepening; Fig. 2), this signal in carbonate platform settings is a single lithologic counterpart only reflecting facies variations, and is moreover modified by diagenetic alteration such as telogenetic dissolutions–reprecipitations. The bulk δ18O does not retain any initial component.
On the other hand, the correspondence between the breaks of the δ13C curve and the sequential boundaries is neither the result of a lithological effect nor the result of diagenesis. This signal could be considered as a promising tool in proximal realm, since it may reflect pedogenesis and original organic activity.
1 Introduction
Au cours des 20 dernières années, la géochimie sédimentaire s'est focalisée sur le développement d'études chimiostratigraphiques à partir de l'analyse de la phase carbonatée des roches totales (bulk analysis ; [4,8,19]). Des comparaisons entre les données géochimiques élémentaires et isotopiques (δ18O, δ13C) et les données sédimentologiques permettent d'effectuer des interprétations en termes de stratigraphie séquentielle.
Ces études chimiostratigraphiques concernent essentiellement les environnements de bassins, où les phénomènes diagénétiques sont peu influents du fait d'une minéralogie homogène en LMC (Low Magnesian Calcite) stable. Les alternances typiques marnes–calcaires reflètent, de plus, une très faible variabilité des milieux de dépôt et limitent les circulations de fluides durant l'enfouissement.
La transposition de cette approche chimiostratigraphique aux plates-formes carbonatées est plus récente [12,16,21,22,24,25], mais l'interprétation des signaux géochimiques est plus complexe et débattue, car : (1) l'enregistrement sédimentaire est plus discontinu et plus diversifié que dans les bassins (hiatus sédimentaires) ; (2) les minéralogies initiales sont multiples (aragonite, LMC, HMC ou High Magnesian Calcite) ; (3) la diagenèse est bien plus développée et diversifiée que dans les bassins. Les données géochimiques obtenues pour les plates-formes sont donc potentiellement le reflet de signaux multiples.
Le but de cette étude est de comprendre la signification des signaux isotopiques de l'oxygène et du carbone des carbonates néritiques en roche totale. Il s'agit d'acquérir ces signaux le long d'une coupe de référence et de les confronter avec : (1) les données lithologiques, (2) des données isotopiques obtenues sur des fractions séparées caractéristiques (calcites initiales, ciments de blocage) ; (3) des données pétrophysiques. L'objectif final est de tenter de déterminer la nature et l'influence des diverses composantes qui contribuent au signal isotopique en roche totale des carbonates néritiques.
2 Cadre de l'étude
L'objet principal de cette étude est le forage HTM 102, réalisé par l'Andra (Agence nationale pour la gestion des déchets radioactifs) en Lorraine (Fig. 1), et qui traverse notamment les formations carbonatées du Jurassique supérieur. L'extension stratigraphique de l'intervalle d'étude débute dans l'Oxfordien moyen, zone à Plicatilis, et se termine dans le Kimméridgien inférieur, zone à Baylei (Fig. 2 ; [9]).
2.1 Environnements de dépôt
L'hétérogénéité des faciès observés (Fig. 2 ; légende Fig. 3) peut être traduite en termes d'évolution verticale et latérale de cinq environnements de dépôt principaux s'étalant du subtidal au supratidal [23]. Les variations environnementales illustrent les grands stades d'évolution de la plate-forme de Lorraine qui repose au nord sur le massif de Londres–Brabant (massif des Ardennes) et est limitée au sud par la faille de Vittel [6,7,10,11]. Cette plate-forme se développe au début de l'Oxfordien moyen avec l'installation de récifs au-dessus des argilites callovo-oxfordiennes ( « argiles de la Woëvre »), et disparaı̂t au Kimméridgien inférieur/supérieur, avec le retour d'une sédimentation terrigène ( « Marnes à Exogyra nana »).
2.2 Interprétation séquentielle
Le motif d'empilement des environnements de dépôt est interprété en terme de cycles stratigraphiques de régression/transgression. Six cycles à moyenne fréquence sont identifiés (Lo1 à Lo6, chacun d'environ 1 Ma, i.e. 3e ordre ; Fig. 2) et sont agencés dans un cycle de plus basse fréquence (2e ordre ; Fig. 2), délimité par deux maxima d'approfondissement situés dans les « argiles de la Woëvre » sous-jacentes [3,18] et dans les « marnes à Exogyra nana » sus-jacentes [13].
3 Signal isotopique (18O, 13C) : acquisition et résultats
3.1 Analyses sur roche totale
Les échantillons bruts sont sciés, lavés, séchés et broyés finement. Une quantité de poudre de 50 mg subit une attaque à l'acide orthophosphorique, durant 6 h ou 15 min, suivant que la réaction se déroule à 25 °C ou 50 °C, respectivement. Le CO2 produit est séparé des autres gaz grâce à un jeu d'extraction impliquant de l'azote liquide (−270 °C), puis de la neige carbonique (−50 °C). Les rapports isotopiques du CO2 ainsi purifié sont ensuite mesurés à l'aide d'un spectromètre de masse. Les analyses ont été effectuées au laboratoire de géologie des bassins sédimentaires (FRE 2400) de l'université Paris-6. Les données sont toutes exprimées par rapport au PDB et les écarts types pour les mesures du δ18O et du δ13C sont respectivement de 0,08 et 0,04‰.
Les courbes d'évolution du δ18O et du δ13C ont été établies le long du forage par l'analyse de 100 échantillons prélevés en fonction des données sédimentologiques. Les ruptures observées sur les deux courbes isotopiques coı̈ncident avec les limites sédimentologiques (Fig. 2). Les intervalles à δ18O les moins négatifs (−2 à −4‰) sont associés aux niveaux lithologiques à matrice micritique (mudstone à packstone), alors que les intervalles à δ18O les plus négatifs (⩽−5‰) correspondent aux niveaux lithologiques à ciments sparitiques (grainstone). La courbe d'évolution du δ13C ne permet pas de faire une analyse du même type. L'utilisation d'un diagramme bivarié des données isotopiques habillées en fonction de leur texture (Fig. 4) montre qu'il n'y a pas d'agencement particulier en fonction des variations du δ13C, contrairement à celles du δ18O (Fig. 4).
3.2 Analyses sur fractions séparées
Les fractions séparées sont prélevées mécaniquement à l'aide d'une fraise de dentiste sous une loupe binoculaire. Il s'agit de LMC originelles de tests d'organismes (brachiopodes, huı̂tres), donc non sensibles à la diagenèse, et de ciments de blocage également en LMC (pores, fractures). Les premières reflètent les compositions isotopiques marines initiales et les secondes reflètent les compositions des fluides baignant les carbonates durant l'enfouissement.
Les moyennes des δ13C des deux ensembles de produits sont proches (Fig. 5). En revanche, leurs δ18O diffèrent, avec des ciments de blocage nettement plus négatifs que les LMC originelles (Fig. 5). Les δ18O très négatifs des ciments de blocage traduisent des solutions parentes de compositions mixtes marine/météorique [23], et contribuent donc à alléger les δ18O en roche totale des carbonates à textures granulaires [15]. Les δ18O des carbonates à textures micritiques sont plus proches des valeurs initiales, mais restent plus négatifs (contrairement aux données de [20], notamment).
4 Discussion
4.1 Relation signal isotopique roche totale–faciès
Les Figs. 6 et 7 illustrent les étendues statistiques (moyenne ± écart type) des mesures de δ18O et δ13C effectuées pour chaque type de texture.
- • La relation entre le faciès et les δ18O apparaı̂t plus clairement et se traduit par l'agencement des textures, depuis les plus fines jusqu'aux plus grossières, en fonction de la décroissance des valeurs de δ18O (Fig. 6). Une ANOVA (ANalysis Of VAriance) réalisée sur ces valeurs montre qu'il existe un effet hautement significatif de la texture des carbonates sur leur valeur en δ18O (Fig. 6). Mais cette relation est associée à la présence/absence de ciment de blocage dans les échantillons, comme le montre une seconde ANOVA, qui atteste un effet hautement significatif de ce paramètre sur les δ18O (Fig. 6).
- • Aucune relation n'est visuellement attestée entre faciès et δ13C. Même s'il existe un effet significatif de la texture sur les δ13C (Fig. 7), il n'est pas lié à la présence de ciment (Fig. 7), mais seulement à la différenciation des valeurs des wackestones (Fig. 7). Il apparaı̂t d'ailleurs que tous les échantillons répertoriés comme wackestone correspondent, (1) soit à des niveaux émersifs (cotes 215 et 119 m, Fig. 2), (2) soit à des niveaux bioconstruits encroûtés (microbialites). Dans le premier cas, les diminutions du δ13C sont dues à une faible pédogenèse associée à ces niveaux émersifs (Fig. 2 ; [1]). Dans le second cas, les faibles δ13C mesurés reflètent probablement l'influence d'un CO2 enrichi en 12C généré par l'activité organique.
4.2 Influence de la diagenèse : modifications des micrites
Une comparaison de la courbe d'évolution des δ18O avec les variations de porosité le long du forage HTM102 (Fig. 2) montre une très grande similarité entre les deux signaux, ce qui nous indique que les variations du δ18O en roche totale ne sont pas uniquement contrôlées par les variations texturales. Cette similitude peut être expliquée par l'étude détaillée d'un intervalle situé entre les cotes 220 et 285 m (Fig. 2), au sein duquel alternent des niveaux poreux à textures micritiques dominantes et des niveaux moins poreux à textures sparitiques dominantes. Si l'on se réfère à la relation δ18O/texture mise en évidence précédemment, la courbe d'évolution des δ18O le long de cet intervalle devrait être très fluctuante, avec des niveaux micritiques nettement moins négatifs que les niveaux sparitiques. Or, la courbe semble lissée et les niveaux micritiques montrent des δ18O aussi négatifs que les niveaux sparitiques (Fig. 2).
Ceci est dû à une particularité des micrites poreuses de l'Oxfordien de Lorraine. Au MEB (microscope électronique à balayage), ces micrites montrent des grains reliques à aspect fondu, qui caractérisent des dissolutions de calcite (Fig. 8), et des ponts micritiques entre les grains qui correspondent à des produits issus de phases de reprécipitation (Fig. 8). Des mesures de 18O ont été réalisées après microprélèvements dans des micrites poreuses (7 échantillons) et compactes (9 échantillons). Les valeurs isotopiques des micrites poreuses sont proches de celles observées pour les ciments de blocage (autour de −7‰) et caractérisent une altération diagénétique, tandis que les δ18O des micrites compactes se rapprochent des valeurs obtenues sur les calcites marines (Fig. 9). Il est probable que les produits issus de reprécipitations aux contacts entre grains, observés dans les micrites poreuses, gardent la trace des fluides météoriques, à rapport isotopique très négatif, qui circulent dans les formations carbonatées de Lorraine dès le Tertiaire, et dont les plus agressifs entraı̂nent des dissolutions importantes [23]. Malgré le fait qu'ils soient volumétriquement moins abondants que les grains reliques, ils paraissent néanmoins susceptibles d'entraı̂ner le δ18O des micrites poreuses vers un pôle plus négatif. Ces données vont à l'encontre de l'hypothèse, classiquement admise, qui considère les nanostructures poreuses lâches comme reflétant la structure quasi initiale des micrites, et les nanostructures plus denses (Fig. 8) comme résultant de modifications diagénétiques plus intenses au cours de l'enfouissement [14,17].
À l'inverse du δ18O, le signal du δ13C en roche totale ne montre aucune similarité avec l'évolution des porosités (Fig. 2), et ne semble donc pas perturbé par ces altérations diagénétiques. Ceci n'est pas étonnant, dans la mesure où la quantité de carbone apportée par les fluides diagénétiques est infime par rapport au réservoir quasi infini que représentent les formations traversées. De ce fait, le rééquilibrage de ces fluides par rapport aux carbonates encaissants est nettement plus rapide pour le δ13C que pour le δ18O [2].
5 Conclusions
Même s'il apparaı̂t que les ruptures du signal isotopique de l'oxygène coı̈ncident avec des limites séquentielles (maxima de régression et d'approfondissement ; Fig. 2), ceci est logique, puisque toutes ces limites correspondent à des variations lithologiques (c'est-à-dire texturales) majeures qui, de plus, délimitent des horizons diagénétiques distincts (Fig. 2). Le δ18O dans les carbonates néritiques n'est qu'une illustration des variations de faciès qui, de surcroı̂t, peut être modifiée par les altérations diagénétiques les plus récentes (dissolutions–reprécipitations télogénétiques). Il ne possède donc aucun souvenir d'une composante originelle.
Pour le δ13C, la correspondance entre limites séquentielles et ruptures du signal (Fig. 2) n'est ni le reflet d'un effet de la lithologie, ni celui d'un effet de la diagenèse. Au contraire, il semble que certaines variations en δ13C puissent refléter l'influence de phénomènes tels que la pédogenèse et l'activité organique. Le possible reflet d'une composante originelle dans ce signal est prometteur, mais d'autres travaux s'avèrent dès à présent nécessaires pour mieux comprendre l'ensemble des processus qui engendrent ses variations.
Remerciements
Nous tenons à remercier les rapporteurs, dont les précieuses remarques ont permis d'améliorer cette note. Ce travail a été soutenu par une bourse de thèse Andra.