A Deep Alteration and Oxidation Profile in a Shallow Clay Aquitard: Example of the Tégulines Clay, East Paris Basin, France

Theoxidation profileof a surficial clay aquitard was studied on a 35-meter borecore fromthe Albian T´egulines Clay near Brienne-le-Chˆateau (Paris Basin, France). Mineralogical, geochemical, and petrophysical data showed evidences of gradual oxidation taking place down to a depth of 20 m. Below 20m, the clay material was nonplastic and nonfractured, and it inherited reduced redox conditions from bacterial sulfate reduction that occurred after sediment deposition. Above 20 m, the clay material was plastic. Up to a depth of 10-11m, only rare yellowish aggregates of glauconite attested to limited oxidation, and pore water chemistry was unmodified. The 5–11m depth interval was characterized by intensive pyrite oxidation, calcite dissolution, and formation of sulfate and iron hydroxide minerals. The upper 2-3 m was ochrous and entirely oxidized. These mineralogical changes were mirrored with pore water chemistry modifications such as an increase of alkalinity and sulfate concentration in the upper part of the profile. The presence of siderite at ∼ 11m evinced the reactivity of Fe(II) in the structure of clay minerals with dioxygen from meteoric waters that infiltrated into the T´egulines Clay through vertical fractures.


Introduction
Clay-rich geological formations act as natural low permeability protective barriers preventing contamination of groundwater resources due to extremely long time scales for groundwater flow and solute transport and the capacity of clay minerals to fix cationic contaminants [1][2][3][4][5][6]. Determining transport mechanisms (advection/diffusion) is first important for evaluating and modelling how contaminants migrate through a geological formation. One way of characterizing the transport mechanisms is to study the vertical profiles of natural conservative tracers in pore waters, such as oxygen and hydrogen isotopes, and dissolved halogens such as chlorine [4,5]. When tracer concentrations do not vary vertically, the dominant transport mechanism may be advection [3], whereas curvilinear concentration profiles suggest that diffusion dominates [1][2][3][4]. In addition to the relevant transport regime, how speciation occurs during transport

Geological Setting
The studied area of Brienne-le-Château (Aube department) is located in the eastern part of the Paris Basin. The Gault clay formation in the studied area consists of siliciclastic shales deposited in an open marine environment at the maximum flooding surface [25] from Middle to Upper Albian (Lower Cretaceous) on the Greensands formation [26] (Figure 1). Temperatures of deep marine waters during Albian were estimated at around 13-18 ∘ C [27]. The limit between Greensands and Gault clay is a major truncative unconformity. An important hiatus between the Aptian and the Albian formations exists in the Aube stratotype due to a general low sea level of eustatic origin [28]. The area that we have studied is the most subsident one in the Paris Basin. The stratotype of the Gault clay defined in the Aube department consists of Argiles Tégulines de Courcelles (82 m) (in the following referred to as Tégulines Clay) overlaid by the Marnes de Brienne (43 m).
The Tégulines Clay consists of three successive facies which are from bottom to top: (1) phosphatic nodule beds, (2) silty clays, and (3) clays and limestone beds. The Marnes de Brienne consists of (1) clay marls, (2) gaize, and (3) clay marls. Phosphatic nodule beds [31], hardgrounds, and glauconitic and sandy beds and layers rich in fossils [32] are currently described in the Tégulines Clay [25]. Several phosphatic nodule beds are continuous over more than 130 km, to Aube in the southwest and Perthois in the northeast. They are interpreted as Flooding Surfaces of 3rd-order cycles [25].
Thick successions of chalk deposited over the Gault clay during Cenomanian-Turonian and Campanian-Maastrichtian transgressions of Late Cretaceous times [33]. Then the basin emerged at the Cretaceous-Tertiary limit, allowing the progressive exhumation of the underlying Gault clay sediments. The retreat of the Upper Cretaceous Sea defined the primitive river network to the east of Paris. The rivers developed on the Cretaceous marine regression surface retreating to the north. Therefore, the rivers appeared sooner in the south and later in the north [34]. At −23 Ma (Early Miocene), the Paris Basin was eroded and similar to the present day. Nowadays the Gault clay formation outcrops as an 8-10 km wide and 80-km long band of terranes oriented NE-SW through the Aube department in the eastern part of the Paris Basin (Figure 1(a)). The thickness increases from east to west between 60 and 120 m, indicating that the total thickness of Upper Cretaceous and Tertiary sediments overlying the Gault clay, and a part of the Gault clay itself, was eroded.
The maximum burial in the Paris Basin was reached during the Upper Cretaceous [35], but the total thickness of the sediments that overlaid the Gault clay in the studied area remains approximate (Figure 1(b)). About 700 m of Late Cretaceous chalk and Tertiary sediments were crosscut by drilling at Poigny, southeast of Paris [36], and 700 m is also the thickness evaluated in the Brie Plaine [33] or in Champagne [35], 600 m in the eastern Paris Basin [37], and 400 ± 300 m in south-eastern Paris Basin [38]. Thermal models indicate thickness of 300-350 m of Late Cretaceous chalk near Troyes based on fluid inclusion microthermometry [39] and 400-600 m near Saint Dizier based on organic matter [40]. Assuming a thickness of about 300-600 m of sediments overlying the Gault clay and a surface temperature of ∼ 23 ± 2 ∘ C [37], the maximum burial temperature did not exceed 30-36 ∘ C. During the Late Tertiary and Quaternary most of the continent emerged, and the disappearance of the sea influence induced general climatic cooling. During the Quaternary the glacial and interglacial periods alternated and influenced erosion and weathering processes.

Drilling Campaigns, Lithology of Tégulines Clay, and
Sampling. Albian Gault clay formation in the area at east of Brienne-le-Château was studied for its mineralogy, its chemistry, and its petrophysical properties. More than 70 samples were collected from six different boreholes (referenced AUB111, AUB121, AUB131, AUB141, AUB151, and AUB161) from the 2013 drilling campaign and of one borehole (referenced AUB1010) from the 2015 drilling campaign (Figure 1(c)). All boreholes were drilled through surficial formations and the Gault clay down to Greensands. Surficial formations were between 0.7 and 6 m thick and consisted of sands, silty clay, and reworked calcareous materials. The Gault clay was only represented by the Tégulines Clay; the Brienne marls were absent. The thickness of the Tégulines Clay crosscut by the seven boreholes ranged between 0 and 80 m, which was influenced by the low dipping of 1 ∘ from southeast to northwest. The thickness was highest in the AUB111 and AUB121 boreholes (65-75 m) located in the northwest of the studied area near Juzanvigny, of ∼35-40 m in AUB141 and AUB1010 boreholes, ∼25-30 m in AUB131 and AUB151 boreholes, and minimal in the AUB161 borehole (<10 m) in the southeast. Aptian Greensands were uncompacted glauconitic quartz sands; however at their upper levels in contact with Tégulines Clay, 0.35 to 1 meters of sand was consolidated. It had quite even thickness of ∼6-10 m in the seven boreholes.
The lithology of the Tégulines Clay was established on the basis of the observations of the AUB111 and AUB121 boreholes for which the formation thickness was the highest (Andra 2015) (Figure 1(d)). The bottom of the clay formation (defined as the UAQ unit) consisted of a 4-5 m thick green silty layer rich in coarse quartz grains and in glauconite at the bottom, followed by 5-7 m of successive small sedimentary cycles of black silty claystones and glauconite-rich silty claystone limited by a bioturbated surface at the top. This was overlain by ∼35-40 m thick dark claystones rich in bioturbations (defined as the UA unit) and containing some layers rich in phosphate nodules. The top of the formation (defined as the UAC unit) only present in boreholes AUB111 and AUB121 consisted of a 20-25 m thick dark grey claystone containing decimeter-thick carbonate layers near the transition with the Brienne marls. Numerous mm to cm sized phosphate nodules occurred at different levels of the clay formation crosscut in boreholes ( Figure 1(d)), which is consistent with literature [31]. Pyrite was not macroscopically identified in the core samples because of its small size, but rare centimeterthick nodules and bioturbations were observed in borehole AUB111, at depth 22.0 m and in borehole AUB1010 at depths 6.63 m and 8.25 m (Figure 1(d)).
Weathering of the Tégulines Clay was visible macroscopically on the core samples, marked by the presence of iron hydroxides and gypsum, and a yellowish rock color. These changes affected the surficial part of the clay formation at variable depths in boreholes: ∼3-4 m for boreholes AUB111, AUB121, AUB131, and AUB141, ∼10 m for boreholes AUB151 and AUB1010, and ∼22 m for borehole AUB161. In this last borehole, the Tégulines Clay was thin, and the weathering also affected the underlying Greensands. In the AUB1010 borehole at 9.25 m, iron hydroxides were evidenced as filling both a subvertical mm thick fracture perpendicular to the bedding, as well as coating planes in the clay-rock bedding (Figure 1(e)).

Sample Conditioning and Analytical
Techniques. The cores were conditioned immediately after being removed from the core sampler to minimize atmospheric contact. The list of analyses and the type of conditioning adapted to the analyses are given in Table 1 and the analytical methods are detailed in Supplementary Appendix 1. Samples conditioned in a Dewar flask containing liquid N 2 (77 K) on the field were directly lyophilized at laboratory. An intact piece of rock of each sample was embedded in a resin for thin section preparation and the remaining parts were crushed into fine-grained powder (<80 m) and stored in a glove box under N 2 atmosphere. Samples conditioned in double-layered heat-sealed plasticized aluminum bags under vacuum. These bags were stored in glove boxes under N 2 atmosphere in the laboratory facilities until analyses.

Rock Characterization
All the major mineralogical, chemical, and isotopic data of clay-rocks are given in Supplementary Appendix 2 (xls file).

Lithology and Mineralogy.
Main mineralogy based on XRD analysis and bulk major element chemistry was obtained for the 77 Tégulines Clay samples. The clay formation consisted dominantly of clay fraction (35-66 wt.% with a mean value of 51 wt.%) associated with quartz-feldspar silty fraction (18-58 wt.% with a mean value of 36 wt.%) and a calcite-dominant carbonate fraction (0-30 wt.% with a mean value of 9 wt.%). Accessory minerals representing less than 5 wt.% of the rock were dolomite, ankerite, siderite, pyrite, and gypsum. Siderite was only a major mineral in sample AUB1010-11.32 m. Total organic carbon was very low and almost homogeneous (0.5 ± 0.1 wt.%). Bulk iron content in the Tégulines Clay was homogeneous at the scale of the formation (∼5 wt.% equivalent Fe 2 O 3 ). The Fe 2+ /(Fe 2+ + Fe 3+ ) ratios of bulk clay-rocks from the AUB1010 borehole are ∼0.8 below 11.32 m and ∼0.2 above.
Mineral variability in the Tégulines Clay was mainly vertical without major lateral variations. We compared the overall vertical variability on the basis of the quantified mineralogy of samples from the AUB111 and AUB121 boreholes for which the Tégulines Clay was thickest, and the three main mineral units are represented (Figure 2).
At the bottom of the clay formation, the 9-12 m thick clay-quartz rich unit (UAQ) was characterized by the highest quartz-feldspar silty content (32-62% with a mean value of 46%) and the lowest carbonate content (0-12 wt.% with a mean value of 3.5 wt.%), a value of almost 0 at its bottom. At the top of the clay formation, the 20-25 m thick carbonateclay-rich unit (UAC) was characterized by the highest carbonate content (9-25 wt.% with a mean value of 17 wt.%) and the lowest quartz-feldspar silty content

Clay Minerals.
Toward the surface, in the first meters of the Tégulines Clay, the ochre color replaced the greenish color, because of the oxidation of iron in clay minerals. The clay fraction of the Tégulines Clay was composed of a quite homogeneous fine-grained matrix of illite-micas (12 ± 5 wt.%), glauconite (8 ± 6 wt.%), illite-smectite mixed layers minerals (16 ± 4 wt.%), and kaolinite (16 ± 5 wt.%), with minor chlorite (<2 wt.%) (Figure 3(a)). The few samples poor in kaolinite and one sample poor in illite/micas/glauconite were from the bottom part of the clay formation (UAQ clayquartz unit). The glauconite content was high and varied a lot at the extremities of the Tégulines Clay (Figure 3 Chemical compositions of the <2 m clay matrix reported in the 4Si-M + -R 2+ diagram modified from Meunier and Velde [19] corresponded dominantly to the field of illite and illite-smectite mixed layers; some others were aligned between the montmorillonite endmember and the chlorite endmember suggesting smectite-chlorite mixing ( Figure 3(e)). Some chemical compositions of the <2 m clay matrix in the AUB131-5 m sample near the surface evolve to the beidellite endmember [19]. This chemical evolution corresponds to a depletion of divalent cations in the clays, probably due to the iron oxidation. Flakes of white mica were on the line between muscovite and phengite. Chlorite flakes generally contained traces of titanium and potassium, strongly suggesting that chlorite derived from biotite alteration. Chemical compositions of glauconite through the Tégulines Clay from the AUB1010 borehole were reported in the M + ⋅(Si/4) −1 versus VI Fe(II + III)/Σ VI cations diagram (Figure 3(g)). Rounded aggregates of wellpreserved green glauconite had homogeneous compositions belonging to the field of highly evolved Fe-poor glauconite according to [4,20,21]. Rounded aggregates of glauconite with yellowish changing color had more heterogeneous chemical compositions giving evidence of a loss of iron and potassium and of an increasing smectite component.

The Quartz-Feldspar Silty
Fraction. This fraction consisted of dominant quartz (34±9 wt.%) with minor K-feldspar (3 ± 0.6 wt.%) generally disseminated in the clay matrix; plagioclase was absent. Quartz was essentially present as detrital angular <100 m sized grains. Dissolution features affecting detrital grains of K-feldspar were observed in samples toward the surface.

Carbonates.
The carbonate fraction consisted of bioclasts and of diagenetic carbonates, represented by major calcite with minor dolomite and siderite. Calcitic bioclasts were locally abundant. They were foraminifer tests of few microns and clasts of shells and tests of variable sizes ranging from few microns to centimeters (Figure 4(a)).
The largest carbonate shells observed by SEM had almost preserved structures. Diagenetic calcite occurred as micrite disseminated in the clay matrix, often associated with pyrite, and as microsparite that filled residual porosity in bioclasts ( Figure 4 . Some nodules contained pyrite-or calcite-infilled septarian veining very similar to previous observations [31]. From XRD of a phosphate nodule we identified dominant fluorapatite with minor quartz, pyrite, an amorphous phase, and traces of kaolinite and muscovite (Supplementary Appendix 3). The diffractogram was slightly different from those of the fluorapatite reference. The lattice parameters = were estimated at 9.326Å (instead of 9.3684Å) and at 6.898Å (instead of 6.881Å). These variations were consistent with the partial substitution of the (PO 4 ) 3− groups by (CO 3 ) 2− groups involving the reduction of the tetrahedral coordination in the a-b plane and low variation along the axis. According to [41], the phosphate is so-called francolite and the (CO 3 ) 2− content was estimated at 9.2 wt.%.

Sulfides and Sulfates.
Pyrite was homogeneously present through the Tégulines Clay, except in samples from near the surface. Its content was estimated to 0.8-1.1 wt.% in the different boreholes, with an average value of 0.9±0.3 wt.%. It occurred as isolated framboids, as elongated clusters of framboids disseminated in the clay matrix (Figure 4(e)), as replacement/infillings of bioclasts (Figure 4(b)), and as rare pyrite nodules. Celestite in association with pyrite was identified by SEM in only few samples ( Figure 5(h)). Pyrite associated with micrite and dolomite is common early diagenetic mineral formed in the bacterial sulfate reduction (BSR) zones of marine clay formations that evolved in reducing conditions [7,8]. Rare celestite associated with framboidal pyrite was probably formed from residual sulfates at the end of the bacterial activity.
Toward the surface, in the first meters of the Tégulines Clay, pyrite was partially to entirely destabilized and replaced by iron hydroxides and gypsum ( Figures 5(a), 5(b), 5(d), and 5(e)). The depth limit of visible pyrite oxidation in the AUB1010 borehole corresponded to the change of bulk Fe 2+ /(Fe 2+ + Fe 3+ ) ratio, which was ∼0.8 at 11.32 m and deeper, and decreased down to 0.2 when shallower. Gypsum occurred as coronas surrounding pyrite (AUB1010-8.25 m; Figure 5(h)) and euhedral grains with sizes ranging from 500 m to 1 mm ( Figure 5(g)). A second generation of celestite was observed as filling secondary porosity due to pyrite dissolution ( Figure 5(f)). The presence of iron sulfates (probable rhomboclase) was suspected in a cm thick oxidized pyrite analyzed by EDS in SEM and Raman spectrometry. Further investigations are needed to confirm it.  [42] allowed us to determine goethite content in the samples and confirmed the goethite content increase between ∼0.8 and 2 wt.% (Table 2).

Stable Isotopes of Diagenetic Minerals.
We analyzed carbon, oxygen, and strontium isotopes of diagenetic minerals including carbonates and phosphates to determine the origin of the diagenetic fluids from which they precipitated.
The 87 Sr/ 86 Sr ratio of a calcite bivalve (0.70736) used to characterize pristine marine fluids is consistent with Cretaceous seawater [43]. The 87 Sr/ 86 Sr ratio of a phosphate nodule (0.70737) is very similar to the signature of initial Albian seawater, indicating the early deposition of phosphates from Albian seawater, in agreement with previous work [31,44].
Carbon and oxygen isotopic analyses of carbonates were carried out on bulk rock due to the impossibility of separating minerals in clay-rock and consequently represented mixing between different generations of diagenetic carbonates.
The calcite 13 C ranged between 0.4 and 2.9‰ PDB consistent with the field of marine sedimentary carbonates. The 13 C of all the dolomites (−3.5 to 0.6‰ PDB ) confirmed a marine origin, except the dolomite in sample AUB1010-11.32 m which had a lower 13 C of −5.6‰ PDB . The siderite 13 C ranged between −11.3 and −7.6‰ PDB so it was systematically lower than marine values. The 13 C lower than marine      [22], dolomite-water [23], and siderite-water [24].
values measured in the dolomite in sample AUB1010-11.32 m and in all the siderite could be rather due to contribution of an organic component.  (Figure 6(a)). Dolomite occurred as early diagenetic mineral present as small amounts and had almost homogeneous 18 O and 13 C through the clay formation, although dolomite/ankerite exhibited dissolution/oxidation in the first ten meters near surface. In unweathered Tégulines Clay, dolomite 18 O of +29/+30‰ SMOW indicated precipitation from marinederived waters at maximum burial of 35 ∘ C (Figure 6(b)). Siderite occurred as flattened crystals disseminated in unweathered samples of the AUB1010 borehole at 13 and 16 m and had 18 O of +29/ + 30‰ SMOW . This 18 O range was consistent with a precipitation from marine-derived waters at maximum burial of 35 ∘ C (Figure 6(c)). The sample AUB1010-11.32 m was peculiar because it contained a network of massive euhedral siderite (16 wt%) and because calcite, dolomite, and siderite exhibited the lowest 18

Physical Properties.
Porosity values deduced from the gravimetric water content (porosity 1) and from density measurements (porosity 2) are consistent at the scale of the studied area. They mainly ranged between 25 and 40% with a trend increasing toward the surface (Figure 7). No correlation was evidenced between mineralogy and porosity. In boreholes AUB111 and AUB121 crosscutting the highest thickness of Tégulines Clay, the clay-rock porosity was quite homogenous below 35-40 m, with values of ∼25-30%. Above 35-40 m and in all other boreholes for which the thickness of Gault clay was lower than 35-40 m, the clay-rock porosity regularly increased from ∼25% up to 35% at 20 m and up to 40% near surface.
The pore size distribution based on mercury intrusion (results not shown) was characterized by dominant pore diameters between 30 and 80 nm in agreement with the clayrich nature of the samples. The highest values of dominant pore sizes were identified in the samples located near the surface (<20 m) and corresponded also to the highest total porosity. Mercury intrusion porosity was about 20% lower than the density and water content based porosity values. This discrepancy could be explained by the presence of pores inaccessible to mercury (e.g., pores with diameter below 3 nm), mainly interlayer pores related to smectite clay phases [45].  alkalinity values and the lowest ionic strengths and may be classified as Ca-Mg-sulfate-Cl type waters. Pore waters in weathered samples of the upper part of the clay formation were characterized by increased alkalinity, increased ionic strength values, and evolved to Ca-Mg-sulfate type waters toward the surface. S 2 O 3 2− species was only measured in pore waters in samples deeper than 20 m, attesting of reduced conditions. Iron and Al were systematically below the detection limit (0.01 mmol/L).

5.2.
Leaching. The C leached,Cl and C leached,SO4 were measured in leaching of sixteen samples with MilliQ water at three solid/water ratios (10, 20, and 50 g/L). Among them, five samples from the first drilling campaign and eleven samples from the AUB1010 borehole were preserved in liquid nitrogen on the field. C leached,Cl and C leached,SO4 were correlated with solid/liquid ratios, suggesting they were not controlled by mineral dissolution/precipitation and adsorption mechanisms.
The C poro,Cl and C poro,SO4 deduced from leaching data from the samples from the AUB1010 borehole and the C bulk,Cl and C bulk,SO4 obtained by squeezing on few samples were in the same range of values, showing that the anion accessible porosity was almost similar to the total porosity ( Figure 8). Assuming C poro,i deduced from leaching data as representative of the pore water chemistry (C bulk,i ), we carried out additional leaching with solid/liquid ratio of 100 g/L on more than thirty samples preserved in aluminum bags under vacuum to complete the anion profiles as a function of depth ( Figure 8).
In the lower part of the Tégulines Clay, the C poro,Cl deduced from the leaching data from all the samples were low (<10 mmol/L PW) and almost consistent, whereas in the upper part of the formation toward the surface, they were more heterogeneous and higher (up to 15 mmol/L PW).
The C poro,SO4 systematically increased toward the surface. Two samples AUB1010-6.63 m and AUB1010-8.25 m exhibit very high C poro,SO4 values (∼520 mmol/L PW), out of range of the other samples. After verification, these two values were due to the presence of gypsum that dissolved during sample leaching; consequently the C leached,SO4 measured in these two samples cannot be considered as representative of the C poro,SO4 , although C leached,SO4 were correlated with solid/liquid ratios. In the lower part of the Tégulines Clay, the C poro,SO4 deduced from leaching of samples conditioned in liquid nitrogen was lower than 10 mmol/L PW, whereas C poro,SO4 deduced from leaching of samples conditioned in aluminum bags were systematically higher (9-30 mmol/L PW), confirming the importance of the sample conditioning [46,47]. These values were not further discussed.

Cation Exchange Capacity and Distribution of Cations
on the Clay Exchanger. The CEC measured by colorimetry (CEC col ) ranged between 10 and 24.9 meq/100 g dry rock with a mean value of 15.0 ± 3 meq/100 g dry rock . The rare highest values were measured either in weathered samples near surface or near the interface with Greensands or in boreholes in which thickness of the Tégulines Clay was reduced due to erosion (AUB161 borehole). This CEC increase could be due (1) to a slight increase of the smectite component in the clay fraction and/or (2) to a slight contribution of iron hydroxides. Smectite component and iron hydroxides increased in weathered samples near surface, whereas smectite component only increased in samples near the interface with Greensands.
The CEC col was almost similar to the sum of exchangeable cations (CEC Σ = Ca + Mg + K + Na) in weathered samples, indicating that CEC measurements were not perturbed by mineral dissolution, except the two samples containing gypsum (AUB1010-6.63 m and AUB1010-8.25 m). The CEC col were systematically higher than the sum of exchangeable cations (CEC Σ = Ca + Mg + K + Na) in reduced samples. The NH 4 + contents measured in reduced samples from the AUB1010 borehole below 10 m (data not shown) were surprisingly high; this excess of NH 4 + corresponded to cobalt hexammine trichloride releasing NH 4 + in response to reduction of its Co(III) and could attest to the good preservation of the initial reduced state of the clay samples [48]. Depth (m) Exchangeable cation concentrations were quite homogeneous through the different boreholes crosscutting the Tégulines Clay (Supplementary Appendix 2). Ca was the dominant adsorbed cation, following by Mg, K, and Na. Toward the surface (<10 m), the CEC and the distribution of exchangeable cations were more heterogeneous. Some samples had high CEC col values associated with higher cation concentrations. High anomalous Ca values were measured in the two samples containing gypsum (AUB1010-6.63 m and AUB1010-8.25 m).

Exchangeable Strontium: Concentrations and 87 Sr/ 86 Sr
Ratios. When coupled with Sr concentrations, Sr-isotope systematics can be used to investigate mixing of different groundwaters [49]. The 87 Sr/ 86 Sr ratios of exchangeable strontium, considered as representative of 87 Sr/ 86 Sr of present-day pore waters [50,51], were measured in Tégulines Clay of the AUB1010 borehole. The exchangeable strontium contents decreased from the bottom of the Tégulines Clay toward the surface; this decrease could be due to dilution of marine-derived pore waters by surficial waters, but also to slight precipitation of celestite. The 87 Sr/ 86 Sr ratios of samples below 15 m ranged between 0.707416 and 0.707450, that is, slightly higher values than the signature of Cretaceous seawater [43]. When shallower (0-15 m), 87 Sr/ 86 Sr ratios slightly increased toward the surface up to 0.708004. The slight enrichment in 87 Sr could originate from the clay formation itself due to weathering processes of clay minerals but could be also due to support of external source such as meteoric waters [52], atmospheric deposition [53], anthropogenic sources [54], and/or mineral dissolution in surficial formations [49].  [29], and Greensand aquifer (Bougligny) [30] are reported.

Oxygen and Hydrogen
heterogeneous and evolve to local present-day meteoric waters. The pore waters 18 O and D of sample at the bottom of the Tégulines Clay at interface with Greensands (35.7 m) are also close to present-day pore waters, suggesting local support of fresh meteoric waters in the Greensands aquifer. The values systematically lower than the meteoric water line and the few values clearly different from the others by their significant 18 O and D enrichment (for example, sample 11.32 m) might be interpreted in terms of climate conditions (rather humid conditions) and evaporation processes. These data need to be validated.

Oxidation Profile.
Redox fronts based on mineralogical (redox-sensitive minerals) and geochemical (redox-sensitive elements such as Fe, Ce, and U) investigations in fractures were highly studied in granitic environments [12,14,55] and in some types of sedimentary rocks [12,56] but less in clay-rocks [13]. For the Tégulines Clay, we favored general macroscopic and mineralogical observations of the clay rather than a targeted study of fractures, as fractures were rare and difficult to identify due to the plasticity of the rock. We were able to collect much "macroscopic" evidence of weathering processes, including oxidation. First, the color changed from ochrous to greenish grey below 5 m. Second, gypsum and iron hydroxides were macroscopically observed down to 9-10 m as in situ replacement of pyrite and iron hydroxides were seen to fill a vertical fracture and planes subparallel to the bedding. Third, we saw a change of consistency from a plastic to a nonplastic clay below ∼20 m.
Combined macroscopic and microscopic observations and chemical, isotopic, and petrophysical data rather indicate that oxidation of the clay formation is characterized more by progressive reactions than by a sharp redox reaction front. In Tégulines Clay, extension of oxidation reactions can be described according to various related effects on the rock properties as observed on the well-preserved reduced clay at ∼21 m to the oxidized clay at surface (Figure 10).
Deeper than 21 m, Tégulines Clay samples were dark green, hard, and reduced. The general texture of the Tégulines Clay and the low burial depth (about 300-600 m) indicated low diagenetic process intensity without extensive precipitation of minerals (such as micritic carbonates). However observations of septarian veining in phosphate nodules attested to compaction processes. The low diagenesis process intensity resulted in larger pore diameters and porosity (28 ± 3%) in comparison to clay-rock strongly influenced by chemical compaction, such as the Callovo-Oxfordian or Opalinus clay formations for instance [57][58][59]. Unweathered Tégulines Clay consisted of more than 75% of detrital minerals including quartz, K-feldspar, kaolinite, chlorite, illitesmectite and chlorite-smectite mixed layers, detrital calcite and phosphate bioclasts, and infratraces of goethite and magnetite. Detrital minerals did not provide evidence of dissolution features. Phosphate bioclasts were often enclosed in nodules of diagenetic phosphate (francolite), indicating they were not stable, whereas calcite bioclasts (shells, tests) were relatively well preserved. Diagenetic minerals consisted of glauconite (8 ± 6%) with variable amounts of calcite and minor dolomite, siderite, pyrite, and phosphates. All these diagenetic minerals are commonly formed at early stage of diagenesis in the bacterial sulfate reduction (BSR) zones of marine clay formations that evolved in reducing conditions [7,8]. The microsparite and celestite, contemporaneous of the compaction, filled the residual porosity. The bulk Fe 2+ /(Fe 2+ + Fe 3+ ) ratio of ∼0.8 characterized the reduced state of the clay-rocks. The dark green color and the chemical composition of glauconite aggregates are consistent with highly evolved Fe-poor glauconite formed in a context of high sedimentation rate according to [20].
From ∼21 m up to ∼10-11 m, Tégulines Clay samples remained dark green and reduced and became plastic. The bulk Fe 2+ /(Fe 2+ + Fe 3+ ) ratio remained close to ∼0.7-0.8 attesting of the stability of the reduced state of the clayrocks. The oxidation degree was low and mainly detectable by the change of consistency due to the progressive increase of water content and porosity (30 ± 3%) toward shallower depths, by slight increase of goethite content, disappearance of magnetite, and observations under the microscope of some yellowish-greenish rounded aggregates of glauconite (transmitted light).
From ∼10-11 m up to 5 m, Tégulines Clay samples were green and plastic. Plasticity is always correlated with the increase of water content and porosity (33 ± 4%) toward shallower depths. The degree of oxidation was high, characterized by observations of some yellowish-brownish rounded aggregates of glauconite, substantial oxidation of pyrite into gypsum and iron hydroxides (major goethite), total oxidation of ankerite into iron hydroxides and partial dissolution of K-feldspar and calcite. The high oxidation is marked by the significant change of the Fe 2+ /(Fe 2+ + Fe 3+ ) ratio from 0.8 to ∼0.2. From 5 m deep and up to the surface, the Tégulines Clay color changed from green to ochrous and was highly plastic. The porosity attained value of 35 ± 4%. Pyrite was entirely oxidized into gypsum and iron hydroxides (goethite ∼2%). The Fe 2+ /(Fe 2+ + Fe 3+ ) ratio was closed to 0-0.1.
Above 9-11 m, the plasticity and porosity changes were the most developed and could be associated with several processes: increase of water interacting with clay particles, nonswelling clay minerals transforming into smectite, minerals dissolving, formation of expansive minerals, and rock physically decompacting due to progressive decrease in the lithostatic pressure. Among them, the pyrite oxidation is the major process involved in the deterioration of clay-rock, due to the formation of expansive phases such as gypsum and goethite [60]. No extensive clay minerals to smectite transformations were evidenced by X-ray diffraction, CEC measurements, and strontium isotopes of exchangeable strontium. The depth of ∼9-11 m is the transition between the highly reactive clay and the slightly oxidized clay. This transition could be also interpreted as the maximum depth of the structural modifications (vertical fracture) of the Tégulines Clay due to weathering processes (e.g., climate changes, saturation/desaturation cycles). Structural modifications might have favored infiltration of oxidizing meteoric fluids at depth and the network of siderite in the AUB1010-11.32 m sample might have been an exceptional example of the reactivity of reduced Tégulines Clay with surficial fluids enriched in iron and in dissolved bicarbonates diffusing through these vertical fractures. The process of siderite precipitation and extension of the vertical fractures need to be validated by further investigations. A distribution of iron (Fe 2+ /(Fe 2+ + Fe 3+ )) among the mineral phases of the Tégulines Clay was approximatively estimated for the unweathered (reduced) Tégulines Clay and weathered clay using a mass balance approach based on bulk iron contents, Fe 2+ /(Fe 2+ + Fe 3+ ) ratios, mineral contents, and mineral chemistry (Table 3). Due to its low content compared to other Fe 2+ -bearing phases, ankerite was not considered in the calculations. Siderite oxidation was assumed to be negligible compared to pyrite oxidation [61]. Glauconite was assumed to be constant through the Tégulines Clay, although it slightly decreased. The oxidation of other Febearing clay minerals was assumed in a first approximation to be negligible compared to glauconite oxidation, although the presence of iron hydroxides was detected in rare large flakes of chlorite, for example. Oxidation essentially affected pyrite (∼90%), glauconite (not determined), and magnetite (100%) and mobilized about ∼7 mmoles of Fe per 100 g of rock, that is, ∼9-10% of the total iron. This liberated iron essentially precipitated essentially in situ in the form of goethite, according to magnetic remanence measurements. Oxidation of rare pyrite-rich layers can mobilize more iron and sulfur.

Consequence of Weathering on Pore Water Chemistry.
Present-day pore waters have strongly evolved from the Albian seawater during early diagenesis, compaction, burial, and uplift due to erosion up to outcrop. Their global history is complex and must consider the past influences of external waters coming from surrounding rocks (Overlying Cretaceous chalk and underlying Greensands), the percolation of meteoric waters, and the water/rock interactions. The previous rock characterization of this study highlighted several zones of mineral reaction due to weathering in the upper part of the Tégulines Clay that must be taken into account to better understand the present-day chemical variations in the pore waters.
The 18 O and D of present-day Tégulines pore waters show that surficial waters have diffused through the clay formation ( Figure 11). In the lower part of the Tégulines Clay (below 20 m in the AUB1010 borehole but not close to the Greensands), early reduced diagenetic conditions are preserved and all the early mineralogy remain stable. Consequently pore waters of this zone can be considered as the most representative of diagenetic pore waters that have been little modified by weathering ( Figure 11). They are characterized by low C poro,Cl (<10 mmol/L) and C poro,SO4 (<10 mmol/L), low ionic strengths, and a cation exchange population dominated by Ca and Mg due to exchange preference for divalent cations over monovalent cations at low ionic strength [62]. 87 Sr/ 86 Sr ratios close to Cretaceous seawater combined with exchangeable strontium contents strongly confirm dilution by meteoric fluids but low water/clay interactions concerning strontium. By comparison, the pore waters in the upper part of the layer (above 20 m in the AUB1010 borehole) are characterized by heterogeneous and higher C poro,Cl and C poro,SO4 , higher alkalinity, higher ionic strengths (up to 175 mmol/L), and a modified cation exchange population. The chemical variations in pore waters toward the surface essentially result from mineral changes due to oxidation. Observations combined with high 87 Sr/ 86 Sr ratios suggest instead that 87 Sr support in pore waters is mainly due to external support but also to dissolution/alteration of Srbearing minerals such as glauconite and K-feldspar. The change of sulfur (S(0) and S(−2) to S(+6)) and iron (Fe(2) to Fe(3)) speciation toward the surface is characterized by the partial to total dissolution of Fe 2+ bearing minerals (pyrite, ankerite, siderite, and magnetite). The consequences of pyrite oxidation following the reaction equation FeS 2 + 7/2 O 2 + H 2 O → Fe 2+ + 2SO 4 2− + 2 H + are the most important [60]. The reaction liberates sulfates, iron ions, and protons. According to mineral characterization, pyrite oxidation in 100 g of Tégulines Clay mobilizes about 7 mmoles of Fe, 14 mmoles of sulfates, and 14 mmoles of H + .
Sulfates partially precipitate as gypsum and minor celestite. Iron precipitates as iron hydroxides, dominantly goethite identified by magnetic remanence. Liberation of protons decreases the pH and by reactivity dissolves the carbonates to produce CO 2 , with the liberation of calcium and possibly magnesium.
The highest variations of pore water chemistry were observed in the 5-10 m zone which is the most reactive ( Figure 11). Above 5 m, pore waters are diluted by infiltrating meteoric waters and controlled by the mineral assemblage formed in oxidized conditions. Below 10-11 m, pore waters result from the slow diffusion of pore waters descending from the 5-10 m reactive zone. At this stage of the work, mineral changes may explain high sulfate concentrations in pore waters, the increasing alkalinity, and the heterogeneities of Ca and Mg of the clay exchanger due to partial dissolution of carbonates (e.g., calcite, dolomite, and ankerite). The varying Cl concentrations in pore waters remain unexplained. They could be due to external support and/or the dissolution of Clbearing phases that have not been identified.

Reactivity of Tégulines Clay to Weathering: Key Solute
Transport Parameters. Understanding the chemical composition of pore waters and its regulation is important for predicting contaminant mobility in rock formations. The speciation, solubility, and retention properties of contaminants are strongly influenced by the pore water chemistry. Spatial evolution in pore water chemistry in Tégulines Clay depends on the present-day hydrogeological context, topography, erosion, and weathering processes. Even interactions between rock and diffusive surficial oxidizing waters remain limited, so the sorption capacity of the clay minerals is preserved. The effect of oxidation on the pore water chemistry can impact the speciation of redox-sensitive contaminants (uranium for instance), so further investigations are required in this field.
Chloride concentrations acquired both by aqueous leaching and by squeezing reveal the low anion exclusion, phenomenon classically described in other clay-rock formations (Callovo-Oxfordian mudstones, Opalinus clay) [63]. The low anion exclusion in the Tégulines Clay can be explained by the geometry of the pore space (and limited diagenetic processes) characterized by pore diameters mainly higher than 30-50 nm. The effect of electrostatic interactions on solute transport appears to be relatively restricted in the Tégulines Clay, potentially easing the undersetting of contaminant transport.
Beyond the uncertainties generated by oxidation on contaminant mobility, weathering-related effects can also help us to better understand and quantify the transport of solute at the scale of the rock formation. Natural tracers ( 18 O and D of pore waters, 87 Sr/ 86 Sr ratios of exchangeable strontium) highlighted evidence of diffusion of meteoric waters through Tégulines Clay, involving the dilution of pristine marine-derived waters. Coupled reactive transport modelling is now required to properly interpret the chemistry of the pore waters and extract quantitative parameters on solute transport, if possible.

Conclusions
Combined macroscopic and microscopic observations and geochemical and petrophysical data allow us to propose a first sketch model of the redox changes in the Tégulines Clay and their consequences on the clay pore water chemistry ( Figure 12). This first sketch model is based on the major ironbearing redox-sensitive minerals such as pyrite, glauconite, siderite, and iron hydroxides, on the valence of iron and on the petrophysical properties. Other redox proxies such as Mn, Ce, and U [12,13], the partial pressure of CO 2 which is a classic parameter used to control pH or alkalinity in clayrocks [64][65][66], trace element behavior, and isotopes sensitive to clay weathering such as Li and Sr are currently investigated. This first sketch model shows that oxidation of Tégulines Clay is characterized more by progressive reactions than by a sharp redox reaction front: (1) The extension of oxidation reactions can be described as observed on the well-preserved reduced clay at ∼20 m to the oxidized clay at surface. This zone corresponds to the transition between a nonplastic to a plastic rock and to glauconite oxidation.
(2) The thickness of the entirely oxidized clay is limited to the first 2-3 meters.
(3) The following 7-8 m corresponds to a highly reactive clay influenced essentially by pyrite oxidation and affected by fracturation.
(4) Siderite network might represent the reaction front between oxidizing meteoric fluids and the reduced clay between a highly reactive zone and reduced clay.
The present-day Tégulines pore water chemistry is heterogeneous and results from the dilution of modified marine waters by oxidizing meteoric fluids diffusing through the clay formation and by water-rock interactions and more specifically by the pyrite oxidation in the highly reactive transition zone.
All these data constitute a set of parameters that should allow us to constrain vertical solute reactive transport modelling through the Tégulines Clay firstly at the scale of a borehole and secondly at the scale of the studied area taking into account erosion and the hydrological network that has been in place for 23 Ma.

Conflicts of Interest
The authors declare that they have no conflicts of interest.