Alexandre Braconnier

Stany Gallier

Christian Chauveau

Fabien Halter

Combustion of a single aluminum droplet burning in O2/CO2/N2 mixtures

Background

Metal combustion has an interesting energetic potential and an elevated temperature flame. According to its high energetic density, aluminum micrometric particles are currently used in industrial solid propellants, especially to increase thrust of solid rocket motors (SRM). However, aluminum particles combustion is a complex process sensitive to many parameters. Indeed it is already proved that combustion time depends on particle diameter. Thus, the initial particle size impacts combustion dynamic and larger particles burn as vapor flame in a diffusion mode whereas smaller particles combustion occurs in a kinetic-limited regime. Relying on experiments, different other factors like pressure or ambient composition change the aluminum burning time and determination of dependence with influencing parameters is an ongoing effort to elaborate accurate correlations, especially in less documented atmospheres such as CO, CO 2 or H 2 O which are oxidizers released by solid propellant combustion. Combustion of a single aluminum droplet burning in O2/CO2/N2 mixtures A. BRACONNIER (1,2)* , S. GALLIER [START_REF] Dreizin | Experimental study of stages in aluminium particle combustion in air[END_REF] , C. CHAUVEAU [START_REF] Legrand | Ignition and combustion of levitated magnesium and aluminum particles in carbon dioxide[END_REF] , F. HALTER [START_REF] Legrand | Ignition and combustion of levitated magnesium and aluminum particles in carbon dioxide[END_REF] • Further investigations to understand and quantify the occurrence of asymmetric phase and its influence on burning parameters varying surrounding composition and pressure. • Autonomous combustion is initiated by using a separated CO 2 laser beam.

Explain and quantify influence of ambient SRM conditions on combustion process analyzing single burning aluminum particle

• Optical signature of burning aluminum particle is recorded using filtered photomultipliers and a high-speed camera.

Experimental set-up

CO 2 laser
Photomultipliers (540 nm / 488 nm)

High speed camera (25000 fps)

HP chamber + levitator

HP chamber

Optical access

Electrodynamic levitator

Ambient Gases :

 O 2  N 2  Ar  CO 2  H 2 O  CO
 Varying pressure in a low range seems to modify related parameters of asymmetric vapor phase. Influence of higher pressure is then less considerable.

 Existence of a symmetric and asymmetric vapor phase (avp) regimes  Introduction of related parameters : t avp and D avp , respectively time and droplet diameter before asymmetric phase occurrence.

:•

  Development of an experimental technique with high-speed image processing • Determination of combustion time dependence to improve numerical simulations • Improvement of phenomenological understanding Objectives • Recording of optical emissions and laser signal during combustion sequences Results (1) ICARE-CNRS, UPR3021, Université d'Orléans, 1C avenue de la Recherche Scientifique , 45071 Orléans, France (2) ArianeGroup, Centre de Recherches du Bouchet, 9 Rue Lavoisier, 91710 Vert le Petit, France *E-mail: alexandre.braconnier@cnrs-orleans.fr

Future works 

 works Combustion time  Temporal identification of singular emissions  Confirmation of laser extinction • Combustion in various mixtures at atmospheric pressure : O 2 -N 2 -Ar  Varying diluent environment, both in composition and proportion, modify combustion times and related parameters of asymmetric vapor phase.  N 2 and Ar species have similar influences on combustion dynamic but N 2 has a more considerable effect, favoring the apparition of asymmetric vapor phase in comparison with Ar. • Highlighting the existence of different stages [2] of combustion • Combustion in air at various pressure : 1 to 31 bars  Combustion time t c  Initial droplet diameter D 0  Evolution of droplet diameter D g  Evolution of flame diameter D f  Identification of phenomena • Recording of dynamic combustion sequences • Single aluminum particle (30-150µm) is isolated in an electrodynamic levitator confined in a high pressure chamber supplied by various gases [1].

  1 bars • D 0 = 90 µm


Xavier Mangenot
email: xavier-mangenot@club.fr

Marta Gasparrini

Virgile Rouchon

Magali Bonifacie

Basin scale thermal and fluid-flow histories revealed by carbonate clumped isotopes ( 47 ) -Middle Jurassic carbonates of the Paris Basin depocentre

The reconstruction of past diagenetic conditions in sedimentary basins is often under-constrained. This results from both the analytical challenge of performing the required analyses on the minute sample amounts available from diagenetic mineral phases and the lack of tracers for some of the This article is protected by

INTRODUCTION

Much of the current knowledge on the thermal, fluid-flow and diagenetic histories of sedimentary basins is based on data, acquired on rocks and minerals, by means of a variety of tools [START_REF] Tucker | Carbonate Sedimentology[END_REF]Moore, 2001). In basins including carbonate successions, the main approach to investigate the thermal and chemical conditions of ancient diagenetic fluids is to examine their fossil archives, by means of the chemical and isotopic compositions of diagenetic minerals and their fluid inclusions [START_REF] Morse | Geochemistry of Sedimentary Carbonates[END_REF][START_REF] Goldstein | Systematics of Fluid Inclusions[END_REF]. In particular, the fluid inclusion microthermometry (leading to both temperature and salinity estimates) and the oxygen isotope compositions of carbonates ( 18 O carbinforming on either temperature or the oxygen isotope composition of the parent fluid δ 18 O water ) are widely used for carbonate successions. These methods are well-established and have been largely applied to many carbonate diagenesis case studies, in different geological contexts [START_REF] Tucker | Carbonate Sedimentology[END_REF]Moore, 2001). Nevertheless, some drawbacks may limit their successful application on carbonate minerals. Indeed, fluid inclusions in carbonate phases are not always present or may not be large enough to be reliably analyzed.

In addition, post-entrapment processes such as stretching, leakage and refill can also affect the reliability of fluid inclusion microthermometry. On the other hand, the application of the δ 18 O carb thermometer requires that the oxygen isotope composition of the parent fluid (δ 18 O water ) from which the mineral formed is known (McCrea, 1950;[START_REF] Epstein | Carbonate-water isotopic temperature[END_REF]. In most cases, δ 18 O water values are under-constrained and strong assumptions are needed because the direct measurements of  18 O water of palaeo-fluids are difficult to determine on a systematic basis (for example, via isotopic composition of O and/or H in crushed fluid inclusions; [START_REF] Affolter | New online method for water isotope analysis of speleothem fluid inclusions using laser absorption spectroscopy (WS-CRDS)[END_REF].

This article is protected by copyright. All rights reserved.

The precise and independent estimations of both the absolute precipitation temperature and the δ 18 O water composition of the carbonate parent fluid can now be acquired from a single method, the carbonate clumped isotopes (denoted as  47 ) thermometry [START_REF] Eiler | Clumped-isotope' geochemistry-The study of naturally-occurring, multiplysubstituted isotopologues[END_REF], thus overcoming the main limitations of the more conventional techniques. The  47 thermometry is based on the thermally-dependant preference of the heavy isotopes of carbon and oxygen 13 C and 18 O to bond with each other in the carbonate lattice at low temperatures. It can be applied to various carbonate minerals (for example, calcite and dolomite) and in a wide range of growth temperatures (0 to 300°C; [START_REF] Kluge | Laboratory calibration of the calcium carbonate clumped isotope thermometer in the 25-250°C temperature range[END_REF][START_REF] Bonifacie | Correction for the 17 O interference in δ( 13 C) measurements when analyzing CO 2 with stable isotope mass spectrometry[END_REF]. The  47 thermometry has been applied to a broad array of studies, mainly focused on the reconstruction of Earth surface palaeotemperatures [START_REF] Passey | High temperature environments of human evolution in East Africa based on bond reordering in paleosol carbonates[END_REF][START_REF] Snell | Hot summers in the Bighorn Basin during the early Paleogene[END_REF] with only few applications to subsurface problems, such as burial and structural diagenesis, water and hydrocarbon fluid-flows, and basin thermal evolution [START_REF] Bristow | A hydrothermal origin for isotopically anomalous cap dolostone cements from south China[END_REF][START_REF] Huntington | Use of clumped-isotope thermometry to constrain the crystallization temperature of diagenetic calcite[END_REF]Loyd et al., 2012;[START_REF] Budd | Syndepositional deformation features in high-relief carbonate platforms: long-lived conduits for diagenetic fluids[END_REF][START_REF] Dale | Time-capsule concretions: unlocking burial diagenetic processes in the Mancos Shale using carbonate clumped isotopes[END_REF]Loyd et al., 2015;[START_REF] Huntington | Carbonate clumped isotope thermometry in continental tectonics[END_REF][START_REF] Swart | Evidence for high temperature and 18 O-enriched fluids in the Arab-D of the Ghawar Field, Saudi Arabia[END_REF].

The aim of the present study is to evaluate the contribution of the Δ 47 tool to the reconstruction of a complete paragenetic sequence, including multiple cement phases with different mineralogy, precipitation temperatures and parent fluid composition, and for which available thermal and fluidflow histories are relatively well-established. For this purpose, the Middle Jurassic subsurface carbonate reservoirs of the central part of the Paris Basin are an ideal geological target, mainly because:

1 They host multiple petrographically well differentiated carbonate phases (i.e. sedimentary biogenic and abiogenic, different burial diagenetic cements) [START_REF] Worden | Cross formational flow in the Paris Basin[END_REF][START_REF] Gaumet | Fondements g ologi ues pour la mod lisa on stra graphi ue des systèmes carbonat s. e Jurassi ue moyen de l'Angleterre la M diterran e[END_REF] This article is protected by copyright. All rights reserved.

2 They belong to a stratigraphic succession for which the geodynamic evolution [START_REF] Guillocheau | Meso-Cenozoic geodynamic evolution of the Paris Basin: 3D stratigraphic constraints[END_REF], formation water geochemistry [START_REF] Worden | Cross formational flow in the Paris Basin[END_REF] and burial histories [START_REF] Uriarte | Maturité thermique des sédiments de la bordure sud-est du Bassinde Paris[END_REF] have been extensively studied.

3 Δ 47 thermometry and fluid inclusion microthermometry acquired on diagenetic carbonates from the same rock unit (five calcitic or dolomitic cement specimens) show very consistent results (Mangenot et al., in press), allowing the Δ 47 data to be confidently used as a solid palaeotemperature (and palaeoδ 18 O water ) proxy for this study.

This study includes detailed petrographic observations (including UV light and cathodoluminescence microscopy) and isotopic analyses ( 18 O,  13 C and  47 ) of several samples of bioclasts, micritic matrix and four different calcitic and dolomitic burial cements. All samples were collected from four different Middle Jurassic cores located tens of kilometres apart in the central part of the Paris Basin.

After examining the consistency of the petrographic and isotopic dataset, which revealed the thermal and palaeofluids composition of each carbonate phase with unprecedented precision, this study discusses the significance of  47 for guiding geological interpretations in the light of the extended knowledge pre-existing for the Paris Basin.

GEOLOGICAL AND HYDROLOGICAL CONTEXT

The Paris Basin is one of the most documented sedimentary basins in the world. It has been studied since the 18 th Century [START_REF] Guettard | Mémoire et carte minéralogique sur la nature et la situation des terrains qui traversent la France et l'Angleterre[END_REF] and is considered a typical example of intracratonic basin [START_REF] Brunet | Subsidence of the Paris Basin[END_REF][START_REF] Pomerol | Naissance et développement de la Géologie du Bassin Parisien au cours des XVIIe et XIXe siècles[END_REF][START_REF] Perrodon | Paris Basin. In: IntraCratonic Basins[END_REF]. It is mainly composed of Mesozoic sediments lying uncomfortably on a Palaeozoic basement [START_REF] Perrodon | Paris Basin. In: IntraCratonic Basins[END_REF] This article is protected by copyright. All rights reserved. [START_REF] Guillocheau | Meso-Cenozoic geodynamic evolution of the Paris Basin: 3D stratigraphic constraints[END_REF]. It has been widely studied as exploration target for oil and gas, hydrological resource, nuclear waste disposal and potential gas storage (natural gas and CO 2 ). The sedimentary column is ca 3000 m thick in the central part of the basin and includes various lithologies spanning from Upper Triassic to Tertiary (Fig. 1A andB). Besides minor uplifts expressed by corresponding unconformities, the Paris Basin experienced a simple burial history punctuated by periods of rapid subsidence in the Jurassic and Late Cretaceous, followed by a major uplift of variable extent during Tertiary tectonism. The four major unconformities have ages corresponding to: the Jurassic-Cretaceous transition (late Cimmerian phase); the Late Aptian (Austrian phase); and the Late Cretaceous to Early Tertiary [START_REF] Guillocheau | Meso-Cenozoic geodynamic evolution of the Paris Basin: 3D stratigraphic constraints[END_REF]. These unconformities have been associated with the early stage of the rise of the thermal dome in the central North Sea, a further stage of Atlantic Ocean opening during the Early Cretaceous and Aptian, and Africa-Europe convergence during the Late Cretaceous and Early Palaeogene [START_REF] Ziegler | Geological Atlas of Western and central Europe[END_REF]. Within the central part of the Paris Basin, these tectonic events led to major compressive flexure and fault reactivation (Fig. 1C). Recent apatite fission track analyses in the southern margin of the Paris Basin revealed that the main tectonic inversion occurred at the Cretaceous-Palaeogene boundary and led to the main uplift phase causing the erosion of 200 to 700 m of Late Cretaceous chalk [START_REF] Barbarand | Lower Paleogene denudation of Upper Cretaceous cover of the Morvan Massif and southeastern Paris Basin (France) revealed by AFT thermochronology and constrained by stratigraphy and paleosurfaces[END_REF].

The last main tectonic event is related to the Late Eocene to Early Oligocene east-west extension, occurring in the eastern part of the basin (Rhine Graben opening), whereas a compressive regime was dominant in the western part [START_REF] Guillocheau | Meso-Cenozoic geodynamic evolution of the Paris Basin: 3D stratigraphic constraints[END_REF]Andre et al., 2010).

The Middle Jurassic carbonates extend across the entire Paris Basin, from exposed outcrops in the eastern and southern margins up to and beyond the English Channel [START_REF] Purser | Plateformes carbonatées: exemple du Jurassique Moyen du Bassin de Paris. Dynamique et méthodes d'étude des bassins sédimentaires[END_REF][START_REF] Javaux | La plateforme parisienne et bourguignonne au Bathonien terminal et au Callovien[END_REF][START_REF] Gaumet | Contrôle stratigraphique des faciès, géométries et profils de dépôt de la plate-forme carbonatée bourguignonne au Bathonien-Callovien[END_REF][START_REF] Gaumet | Fondements g ologi ues pour la mod lisa on stra graphi ue des systèmes carbonat s. e Jurassi ue moyen de l'Angleterre la M diterran e[END_REF][START_REF] Delmas | Paris Basin. Rapport r gional d' valuation p trolière[END_REF]. In the depocentre area, the Middle Jurassic carbonates are documented oil and geothermal water reservoirs, which have experienced a complex and still debated hydrogeological history (Fontes & This article is protected by copyright. All rights reserved. [START_REF] Fontes | Geochemistry and origin of formation brines from the Paris Basin, France[END_REF]Matray et al., 1994;[START_REF] Worden | Cross formational flow in the Paris Basin[END_REF]. The present-day temperature of the Middle Jurassic strata in the central part of the basin is around 60 to 80°C (Ménétrier et al. 2005;[START_REF] Garibaldi | Détermination des températures profondes du bassin du Sud-Est de la France et relations entre anomalies thermiques, géologie et circulations hydrothermales par modelisation 3[END_REF]. However, organic maturity data from Lower Jurassic (mainly Hettangian and Toarcian) shales (Fig. 1F), as well as fluid inclusions and clay minerals thermometric data from Middle Jurassic rocks of the depocentre consistently indicate that they have experienced higher burial temperatures, up to 80 to 100°C [START_REF] Espitalie | Aspects of hydrocarbon migration in the Mesozoic in the Paris Basin as deduced from an organic geochemical survey[END_REF][START_REF] Guilhaumou | Détermination de paléotempératures dans les roches-mères du bassin de Paris: etude d'inclusions fluides et implications pour l'histoire thermique du bassin[END_REF][START_REF] Gaulier | Modeling present and past thermal regimes in the Paris Basinpetroleum implications[END_REF][START_REF] Clauer | Clay minerals as records of temperature conditions and duration of thermal anomalies in the Paris Basin, France[END_REF][START_REF] Worden | Cross formational flow in the Paris Basin[END_REF][START_REF] Spötl | Clay minerals as records of temperature conditions and duration of thermal anomalies in the Paris Basin, France[END_REF][START_REF] Gonçalvès | Fluid inclusions as constraints in a three-dimensional hydro-thermo-mechanical model of the Paris basin, France[END_REF].

The other main reservoirs in the central part of the Paris Basin are the Carnian Chaunoy and the Rhaetien fluvial sandstone units, here grouped together and referred to as the Upper Triassic reservoirs. A Triassic (Keuper) halite body sits to the east and is separated from the Upper Triassic reservoirs by interbedded sands and anhydrite-rich mudstones [START_REF] Worden | Basin scale evolution of formation waters: a diagenetic and formation water study of the Triassic Chaunoy Formation, Paris Basin[END_REF]. Sandwiched between the Upper Triassic sandstone reservoirs and the Middle Jurassic carbonate reservoirs, are Lower Jurassic Hettangian and Toarcian shales, which represent the main source rocks within the basin [START_REF] Poulet | Hydrocarbon migration in the Paris Basin. In: Migration of Hydrocarbons in Sedimentary Basins[END_REF][START_REF] Espitalie | Aspects of hydrocarbon migration in the Mesozoic in the Paris Basin as deduced from an organic geochemical survey[END_REF]. The Lower Jurassic source rocks of the Paris Basin only experienced oil window temperatures (90 to 115°C) for the Toarcian shales (Ménétrier et al., 2005). Hydrocarbon maturation started during the Early Cretaceous and reached the peak generation at the beginning of the Late Cretaceous (90 to 70 Ma; Monticone et al., 2011). Oil migration possibly started just after, though the main migration seems to have occurred in the Oligocene during an east-west extensional phase [START_REF] Tremolières | Mécanismes de la deformation en zones de plate-frome. Méthode et application au Bassin de Paris[END_REF].

This article is protected by copyright. All rights reserved.

SAMPLES AND METHODS

Sample location and sampling procedure

Four sub-surface rock cores were selected from the depocentre of the Paris Basin, corresponding to the following exploration and production wells (present depth of 1500 to 1900 m): Baulne en Brie, Rigny la Nonneuse, Villeperdue and Fossoy. The studied core intervals belong to the same Middle Jurassic lithostratigraphic unit (carbonates of the Comblanchian Formation, Upper Bathonian), and consists of similar sedimentary facies (see details in the Petrographic analysis section below).

Location of the cores, present-day burial depth of the Middle Jurassic sampled intervals and T max values (from Rock-Eval pyrolysis) of the underlying Lower Jurassic shales are reported in Fig. 1. The studied area covers more than 5000 km 2 . A total of 70 rock samples were collected from the four cores to create thin sections (30 m thick) for petrographic analysis. Geochemical analyses were accomplished on different carbonate phases, manually sampled with a dental drill from the highly polished rock slabs used for thin section preparation. A total of 46 samples from both sedimentary and diagenetic carbonates (39 cements of calcite or dolomite, three bioclasts and four micrites) were selected and drilled for isotopic analyses.

Petrographic analysis

All of the prepared thin sections were carefully analyzed using conventional transmitted light and cathodoluminescence microscopy (CL). This was essential to select homogeneous specimens, and avoid phase mixing during manual sampling, prior to geochemical analysis. CL apparatus used in the Δ 47 ,  18 O and  13 C measurements and data processing Measurements of Δ 47 ,  18 O and  13 C compositions of carbonate samples were performed at Institut de Physique du Globe de Paris (IPGP, Stable Isotope Team) with a Thermo Scientific MAT253 gassource mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) after digestion of carbonate powder in 104% phosphoric acid. About 5 mg of carbonate were reacted at 90°C in a common acid bath for 20 minutes for calcites and 1 hour for dolomites. The methods used for carbonate digestion, CO 2 purification and isotope measurements follow the procedure detailed in [START_REF] Bonifacie | Correction for the 17 O interference in δ( 13 C) measurements when analyzing CO 2 with stable isotope mass spectrometry[END_REF]. Each measurement consisted of 70 cycles of comparison between the CO 2 extracted from the sample against a working internal reference CO 2 gas (Oztech; with  13 C = -3.72‰ and  18 O = -6.06‰, versus VPDB, verified with the international carbonate reference material NBS19). The signal integration time was 26 seconds (i.e. total integration time of 1820 seconds for each CO2 sample) for a signal of 12V on m/z = 44. Each carbonate sample was analyzed two or three times to determine its  47 ,  18 O and  13 C compositions in eight discrete analytical sessionsTraditional δ 18 O and δ 13 C carbonate compositions were acquired as part of each Δ 47 analysis and 17 O corrections were made using the parameters from [START_REF] Santrock | Isotopic analyses based on the mass spectrum of carbon dioxide[END_REF]. In order to account This article is protected by copyright. All rights reserved.

for the temperature dependence of oxygen isotope fractionation between CO 2 gas and carbonate, resulting from the reaction with phosphoric acid at 90°C, a fractionation factor of 1,00811 for calcite and 1,0093 for dolomite was used (following Rosenbaum & Sheppard, 1989;[START_REF] Katz | Laboratory-grown coccoliths exhibit no vital effect in clumped isotope (∆ 47 ) composition on a range of geologically relevant temperatures[END_REF]. The  18 O and  13 C of the carbonate samples are expressed in ‰ relative to the VPDB standard.

For constructing correction frames for processing raw ∆ 47 values, CO 2 gases driven to isotopologue equilibrium both at 1000°C and 25°C were used. The equilibrated CO 2 gas standards in this study have bulk isotopic compositions spanning the entire range of measured cement samples, were purified and analyzed in the same way as carbonate samples or carbonate standards, and were typically run every four to five analyses (Table S1). The raw ∆ 47 data were first corrected for linearity effects using a fixed common Equilibrated Gas Line slope fitted to the equilibrated CO2 gases at both 1000°C and 25°C. Subsequently, as recommended by [START_REF] Dennis | Defining an absolute reference frame for 'clumped' isotope studies of CO2[END_REF], the raw ∆ 47 values (expressed relative to the working gas) were transferred into the Carbon Dioxide Equilibrated Scale (CDES) using the CO2 gases driven to isotopologue equilibrium both at 1000°C and 25°C with theoretically predicted  47 values of 0.0266‰ and 0.9252‰, respectively (after Wang et al., 2004).

Finally, ∆ 47 data were projected into the 25°C acid digestion reference frame for easier comparison with previously published calibration data. For this, an acid fractionation value of +0.092‰, determined by [START_REF] Henkes | Carbonate clumped isotope compositions of modern marine mollusk and brachiopod shells[END_REF] was added.

To ensure accuracy of the entire data reduction process, and ultimately the accuracy of the ∆ 47 data presented here, two carbonate reference materials, also reported by [START_REF] Dennis | Defining an absolute reference frame for 'clumped' isotope studies of CO2[END_REF] and many other studies, were routinely analyzed (IPGP-Carrara marble and 102-GCAZ01b). One of these two carbonate standards was analyzed typically every five analyses and distributed along the diagenetic cement samples in all runs in order to check for analytical stability/accuracy of the whole procedure This article is protected by copyright. All rights reserved.

(including carbonate digestion, CO 2 purification, stability of the conditions for analyses of CO 2 inside the mass spectrometer and/or accuracy of the correction frames constructed with standards of equilibrated CO 2 gas -namely the accuracy of the equilibrated gas lines and empirical transfer function lines), as well as long-term reproducibility of the ∆ 47 measurements. The Δ 47 values obtained for these carbonate reference materials over the period of this study (July 2014to December 2015) are: Δ 47 = 0.413 ± 0.016‰ (1SD, n =34) for IPGP-Carrara and Δ 47 = 0.719 ± 0.019‰ (1SD, n=27) for 102-GC-AZ01b. These ∆ 47 values are indistinguishable from those obtained in previous studies (e.g. [START_REF] Dennis | Defining an absolute reference frame for 'clumped' isotope studies of CO2[END_REF][START_REF] Henkes | Carbonate clumped isotope compositions of modern marine mollusk and brachiopod shells[END_REF] and on a longer timescale at IPGP [START_REF] Bonifacie | Correction for the 17 O interference in δ( 13 C) measurements when analyzing CO 2 with stable isotope mass spectrometry[END_REF].

Noticeably, the recent work of [START_REF] Daëron | Absolute isotopic abundance ratios and the accuracy of  47 measurements[END_REF] and [START_REF] Schauer | Choice of 17 O correction affects clumped isotope (Δ 47 ) values of CO 2 measured with mass spectrometry[END_REF] underline that when the 17 O parameters from Brand et al. (2010) are used to process isotopic data (instead of the 17 O parameters from [START_REF] Santrock | Isotopic analyses based on the mass spectrum of carbon dioxide[END_REF] commonly used in the clumped isotope community to date) small differences in  47 values may arise, with direction and magnitude depending on the bulk isotopic composition of the samples and of the standards of equilibrated CO 2 gas used to transfer the data into the CDES reference frame. However, it should be noted that such reprocessing only led to small differences in  47 values of the standards and samples (i.e. less than ±0.010‰), within the analytical reproducibility on  47 measurements, and thus does not significantly change the reported temperatures (and conclusions) of this study.

Finally, the corrected  47 values were then converted into temperatures (T 47 ) using the interlaboratory composite  47 -T calibration recently determined for all carbonate minerals for the 0 to 300°C temperature range (i.e. equation 3 from [START_REF] Bonifacie | Correction for the 17 O interference in δ( 13 C) measurements when analyzing CO 2 with stable isotope mass spectrometry[END_REF] that compiles 103 mean  47 data from seven different laboratories and with proper error propagation -that is  47 = 0.0422*10 6 /T + 0.2182 in the 25°C acid digestion reference frame). The oxygen isotopic compositions of the water [ 18 O water, expressed relative to VSMOW(Vienna Standard Mean Ocean Water)] from which the carbonates precipitated were calculated for each estimated T 47 using the This article is protected by copyright. All rights reserved.

 18 O carb values measured for the carbonate, as well as the oxygen isotopes fractionation between the carbonate and water from O' Neil et al. (1969) for calcites and Horita (2014) for dolomites.

RESULTS

Petrography and paragenesis

The four studied cored intervals correspond to the Upper Bathonian (Middle Jurassic) Comblanchian Formation [START_REF] Purser | Plateformes carbonatées: exemple du Jurassique Moyen du Bassin de Paris. Dynamique et méthodes d'étude des bassins sédimentaires[END_REF][START_REF] Gaumet | Contrôle stratigraphique des faciès, géométries et profils de dépôt de la plate-forme carbonatée bourguignonne au Bathonien-Callovien[END_REF][START_REF] Gaumet | Fondements g ologi ues pour la mod lisa on stra graphi ue des systèmes carbonat s. e Jurassi ue moyen de l'Angleterre la M diterran e[END_REF]. The reconstructed depositional profile includes, from proximal to distal: (i) a grain-dominated lithoclast rich intertidal deposit with bird's-eye and fenestral porosities; (ii) a mud-dominated succession deposited within a protected lagoon rich in green algae (Cayeuxias and Dasycladaceas) and benthic foraminifera (Meyendorfinas); and (iii) an open platform deposit rich in crinoids, isolated corals and shell fragments from brachiopods and bivalves. The intertidal to lagoonal facies locally display bird's-eye structures with geopetal infill by internal micritic sediment and calcite spars. The porosity evaluated from petrographic observations is low (<5%) and was modified through time by a broad spectrum of diagenetic processes, including multiple cementation episodes, dissolution, neomorphism, mechanical and chemical compaction and micro-fracturing. In this study, only two brachiopods and one bivalve, with no hints of diagenetic alteration (i.e. showing non-luminescent pristine fibres; see Fig. 2A andB), were used for the geochemical characterization of syn-sedimentary carbonates.

Presently, most of the primary and secondary pores are filled by calcite and dolomite cements (Figs 2 to 4). Petrographic relationships between carbonate cements and other diagenetic products allowed the establishment of a complete paragenetic sequence, composed of 14 stages (Fig. 5), ranging from syn-sedimentary to late burial processes in the current study.

This article is protected by copyright. All rights reserved.

A more detailed petrographic study was focused on the five main generations of calcite and dolomite cements, named Cal0, Cal1, Dol1, Cal2 and Dol2. The petrographic features of these cements are illustrated in Figs 2 to 4. Although these cements do not always occur in each of the four investigated cores, a generalized mineral paragenesis, based on CL cement stratigraphy, was reconstructed (Fig. 5).

The Cal0 cementation event includes two petrographically distinct non-ferroan calcites, lining abiotic grains and fenestral or mouldic pores (Fig. 2C). The former is defined by a rim of non-luminescent bladed calcite crystals (10 to 40 m) with zoned bright yellow CL. The latter is characterized by a rim of equant granular crystals with a uniform dull brown to red CL. Due to their habit and CL response these calcite cements possibly precipitated during eogenesis, respectively from meteoric and marine fluids (as also reported by [START_REF] Tucker | Carbonate Sedimentology[END_REF][START_REF] Tucker | Carbonate Sedimentology[END_REF]Flügel, 2004). They represent a volumetrically minor cementation event and could not be sampled for geochemical analyses. However, they possibly played a role in promoting early induration of the muddy lagoonal sediments, that reduced the impact of successive mechanical compaction during burial.

Cal1 is composed of blocky crystals of non-ferroan calcite, that filled a variety of pore types: interparticle, fenestral, mouldic and vuggy. This cement induced an important reduction of the remaining pore space and seems to be the main cause for the present-day relatively low porosity of the studied carbonate reservoirs. It can be subdivided into two petrographic types, named here respectively Cal1(a) and Cal1(b). Type Cal1(a) is rarely observed and consists of crystals (10 to 100 m) displaying concentric zonation under CL, with bright orange and dull orange zones. In contrast, Cal1(b) is more common and consists of crystals (100 m to 2 mm) which display a homogeneous bright orange CL (Fig. 2C to E). Type Cal1(b) is referred simply as Cal1 in the following text, figures and tables, becuaes it is the only petrographic type of Cal1 that could be further investigated with geochemical tools.

This article is protected by copyright. All rights reserved.

Type Dol1 is composed of dolomite crystals (200 m to 2 mm) which display curved crystal faces (Fig. 4B andC) and an undulose extinction, typical of saddle dolomites [START_REF] Searl | Saddle dolomite : a new view of its nature and origin[END_REF]. The crystals exhibit a concentric zoned CL, dominantly non-luminescent with some thin red zones (Fig. 4B). This dolomite cement is non-ferroan and mostly occurs within fenestral pores, commonly showing geopetal features. Primary hydrocarbon bearing fluid inclusions were observed in Dol1 crystals (Fig. 4C), in the Villeperdue and Fossoy cores. These inclusions display shapes (rectangular to negative crystal) controlled by crystallographic planes and are assembled along crystal growth zones.

Type Cal2 is composed of non-ferroan blocky calcite crystals (100 m to 1 mm), mostly precipitated in pore spaces created by diagenetic modifications occurring after Cal1 precipitation (for example, stylolite planes opening, micro-fracturing and partial dissolution of Cal1). Cal2 cement displays a uniform dull red to brown CL (Fig. 3) with only minor sectorial zoning. Type Dol2 is composed of saddle dolomite crystals (100 to 500 m) precipitated as overgrowths on Dol1 with a uniform bright red CL (Fig. 4A). It represents the last carbonate cement of the reconstructed mineral paragenesis. Dol2 cement is rarely observed and never voluminous. Although it was petrographically identified in the four studied cores, it could be sampled for geochemical purposes only from the Rigny la Nonneuse core. Although most of the above described cements have been petrographically recognized in the four studied cores, they were not always voluminous enough to be sampled for geochemical analyses.

This article is protected by copyright. All rights reserved. S1). Remarkably, each of the six investigated carbonate phases shows a discrete combination of their  18 O carb ,  13 C carb , T 47 and  18 O water values with only limited variations among samples of each phase (although the investigated samples are from four different core sections, tens of kilometres apart). Table 2 reports, for each of the six carbonate phases, their average  18 O carb ,  13 C carb ,  47 , T 47 and  18 O water values (obtained from two to 21 samples, and that represent a total of three to 51 isotopic measurements for a given carbonate phase) together with their respective one standard deviation of the mean.

The carbon isotope compositions of all the different calcite cements (average  13 C carb = 1.81 ± 0.34‰, n = 32 samples for Cal1 and Cal 2; Table 2) fall in the same range as the micrites (average  13 C carb = 1.75 ± 0.27‰, n = 4 samples) but are significantly lower than for the dolomite cements (average  13 C carb = 3.2 ± 0.4‰, n = 6 samples for Dol1 and Dol2) (Fig. 6). In contrast, the six discrete carbonate phases investigated show a significant scatter in their oxygen isotope compositions the three analyzed bioclasts fall within the range of expected values for biogenic calcites, precipitated in equilibrium with Middle Jurassic seawater [START_REF] Veizer | 87 Sr/ 86 Sr, δ 13 C and δ 18 O evolution of Phanerozoic seawater[END_REF].

Similarly to the  18 O carb values, the different carbonate phases investigated show a significant scatter in their  47 compositions (i.e.  47 varying from 0.524 to 0.678‰ in all the investigated samples), with however very clustered  47 values when only one carbonate phase is considered (with one standard deviation of the mean averaging at ± 0.009‰). Remarkably, each of the six carbonate phases shows a discrete combination of tightly clustered  18 O carb and  47 values (Fig. 7). As such, marine bioclasts yield average  47 values of 0.663 ± 0.013‰ (1 standard deviation; Table 2), corresponding to average temperature of 29 ± 4°C. Samples of the surrounding micrite display strictly lower  47 values (mean  47 = 0.636±0.007‰, 1 standard deviation), corresponding to significantly higher temperatures of 45 ± 3°C. The other groups correspond to different phases of diagenetic cements with average  47 values of 0.587 ± 0.009‰ for Cal1, 0.537 ± 0.009‰ for Dol1, 0.557 ± 0.008‰ for Cal2, and 0.582 ± 0.007‰ for Dol2 (with all uncertainties reported as 1 standard deviation with n = 21, 4, 9 and 2, respectively). Those very homogeneous  47 values inside each carbonate phase can be converted into precise temperature estimates, with T 47 of 64 ± 4°C for Cal1, 90 ± 6°C for Dol1, 78 ± 4°C for Cal2, and 67 ± 3°C for Dol2 (with all uncertainties reported as 1 standard deviation of the mean). For each investigated carbonate phase, the reported standard deviation of the mean for  47 data (between 0.007‰ and 0.009‰, typically resulting from a total number of three to 51 analyses depending on the generation studied) are lower than those obtained on homogeneous carbonate standards (with typical 1 standard deviation of ±0.017‰ on n >50 measurements in the present study). This further suggests the high homogeneity of the carbonates investigated here at both the sub-sample scale and between the different co-genetic cement phases collected in different core sections located tens of kilometres apart.

This article is protected by copyright. All rights reserved.

The oxygen isotopic compositions of the parent fluids for the six carbonate phases investigated exhibit variable values, with respective average δ 18 O water and their associated one standard deviation of the mean of: +1.8 ± 1.4‰ for the bioclasts, +2.3 ± 0.7‰ for the micrites, +2.5 ± 1‰ for Cal1, 0.5 ± 0.4‰ for Dol1, -4.1 ± 0.5‰ for Cal2, and -7.0 ± 0.8‰ for Dol2. As illustrated in Figs 7B and8, the carbonate parent fluids show slightly positive O isotopic compositions for the syn-sedimentary carbonates (bioclasts and micrites) and Cal1 blocky calcites. The parent fluids become then slightly depleted in 18 O going from Cal1, to Dol1 and Cal2, to reach the most negative  18 O water values for Dol2, which is the last phase of the mineral paragenesis. Because each of the six carbonate phases investigated here are highly homogeneous in their  18 O carb ,  13 C carb and  47 data, all the T 47 ,  18 O water as well as  18 O carb and  13 C carb values reported and discussed hereafter are values averaged by phase.

PREVIOUSLY PUBLISHED MICROTHERMOMETRIC DATA

Microthermometry analyses have been accomplished by previous authors (Matray et al. 1989;[START_REF] Demars | volu on diag n ue, pal o uides et pal othermicit dans les r servoirs du Keuper et du Dogger du bassin de Paris[END_REF][START_REF] Worden | Cross formational flow in the Paris Basin[END_REF]Mangenot et al., in press) for two-phase aqueous fluid inclusions in carbonate cements from the Middle Jurassic carbonates of the Paris Basin depocentre.

The results, including homogenization temperatures (Th) and total salinities, are compiled in Fig. 9. These former microthermometry data (from Matray et al. 1989;[START_REF] Demars | volu on diag n ue, pal o uides et pal othermicit dans les r servoirs du Keuper et du Dogger du bassin de Paris[END_REF][START_REF] Worden | Cross formational flow in the Paris Basin[END_REF] This article is protected by copyright. All rights reserved.

Fluid inclusions from undifferentiated cements measured by Matray et al. (1989), [START_REF] Demars | volu on diag n ue, pal o uides et pal othermicit dans les r servoirs du Keuper et du Dogger du bassin de Paris[END_REF] and [START_REF] Worden | Cross formational flow in the Paris Basin[END_REF] display a large range of homogenization temperatures (55 to 120°C), with a broad mode and salinities (0 to 25 wt% NaCl eq), with a multi-modal distribution (Fig. 9A.1 and 9B.1). Although the microthermometry data from Mangenot et al. (in press) globally overlap those reported from these previous authors, more precision is gained since the three carbonate cements are considered separately (Fig. 9A.2 and 9B.2). Indeed, both homogenization temperatures and salinities, acquired with such petrographic control, are more consistent (Fig. 9A.2 and 9B.2).

Mangenot et al. (in press

) also evaluated the crystallization temperatures for the same five carbonate samples via clumped isotope (Δ 47 ) analysis. This study reveals a very good consistency between the thermal information acquired from these two independent palaeothermometers. The usefulness of these previous studies for the present survey is two-fold:

1

DISCUSSION

Carbonates formed in surface and sub-surface settings record the temperature and isotope composition of the fluids from which they precipitated (i.e. the parent fluids). The samples investigated here recorded exceptional petrographic homogeneity and discrete average Δ 47 values for the six carbonate phases identified. The temperature homogeneity suggests similar burial temperatures characterizing the host-rocks at the time of cementation, whereas the compositional homogeneity possibly reflects rapid fluid circulation episodes accompanied by little interaction with the host-rocks. This allowed reconstruction of a detailed thermal and compositional evolution of the fluids circulating in the Paris Basin through time, from depositional to current burial settings. The following sections interpret the findings of this survey in the framework of the known geological and hydrological history of the Paris Basin. This provided constraints on the timing of carbonate cementations, the parent fluid sources and the evolution of the water/rock systems during the main subsidence and uplift stages of the basin history.

Thermal conditions of carbonate crystallization

Based on CL microscopy, consistent petrographic features and a relative paragenetic sequence were deduced at basin scale for the different cement phases (i.e. Cal1, Dol1, Cal2 and Dol2), and a comprehensive cement stratigraphy was proposed (Fig. 5). The paragenesis was then updated considering the  47 temperatures measured on some samples of these cements and is shown in Fig. 10. Overall, the data suggest that the investigated Middle Jurassic carbonates (from four stratigraphically equivalent core intervals) recorded in a similar manner several fluid-flow episodes, leading to corresponding cementation, and spanning from Jurassic to present time. However, in This article is protected by copyright. All rights reserved. some examples the CL pattern of carbonate cements has revealed not to be fully reliable to establish cement stratigraphy at basin scale (e.g. [START_REF] Frank | Chronostratigraphic significance of cathodoluminescence zoning in syntaxial cement: Mississippian Lake Valley Formation, New Mexico[END_REF]. Therefore, it should be stressed that such basin scale cement stratigraphy correlations should preferentially rely also on geochemical data ( 18 O carb ,  13 C carb and  47 ), when possible.

Preservation of the initial T 47 and  18 O water values

As rocks experience elevated temperatures during deep burial and subsequent cooling during uplift, their  47 composition may be altered by solid-state reordering of C-O bonds within the mineral lattice [START_REF] Passey | Carbonate clumped isotope bond reordering and geospeedometry[END_REF][START_REF] Henkes | Temperature limits for preservation of primary calcite clumped isotope paleotemperatures[END_REF][START_REF] Stolper | The kinetics of solid-state isotope-exchange reactions for clumped isotopes in carbonates: a study of inorganic calcites and apatites from natural and experimental samples[END_REF] and thus compromise the determination of the δ 18 O water of the parent fluid. For the carbonates studied here, the possibility of solid-state diffusion alteration can be ruled out given the three following arguments. First, current knowledge on solid-state diffusion predicts that carbonate could start to lose its original  47 composition only if at ca 80 to 100°C for more than 100 Ma [START_REF] Passey | Carbonate clumped isotope bond reordering and geospeedometry[END_REF][START_REF] Stolper | The kinetics of solid-state isotope-exchange reactions for clumped isotopes in carbonates: a study of inorganic calcites and apatites from natural and experimental samples[END_REF], conditions that are well above those experienced by the studied carbonate unit (as independently determined by Ménétrier et al., 2005;[START_REF] Garibaldi | Détermination des températures profondes du bassin du Sud-Est de la France et relations entre anomalies thermiques, géologie et circulations hydrothermales par modelisation 3[END_REF]Goncalvès et al., 2010).

Second, the excellent consistency between data from fluid inclusion microthermometry and  47 thermometry (Mangenot et al., in press), two independent methods based on different principles, strongly suggests that the measured  47 values of the carbonates were acquired at the time of crystallization. Third, the bioclasts investigated here show T 47 of ca 29°C, that agrees with growth conditions in a lagoonal setting. Consequently, the δ 18 O water of the carbonate parent fluids can be meaningfully estimated.

This article is protected by copyright. All rights reserved.

Syn-depositional and early diagenetic carbonates

As suggested by the preservation of primary prismatic structures of the bivalve shell and the lack of luminescence of the brachiopod shells (pristine fibres; Fig. 2A andB), the three marine bioclasts analyzed do not display petrographic evidences of post-depositional modifications. However, the possibility that a minor amount of syn-depositional marine cement occurs in the punctae of the brachiopod shell is not excluded. The growth temperatures of 29 ± 4°C calculated from  47 also agree with a primary marine origin, possibly reflecting Middle Jurassic seawater temperatures in a restricted lagoon environment. Seawater temperatures of 15 to 27°C were independently deduced for the Tethys ocean during Bathonian time from carbonates located at latitudes similar to those of the investigated samples [START_REF] Picard | δ 18 O values of coexisting brachiopods and fish: temperature differences and estimates of paleo-water depths[END_REF][START_REF] Lécuyer | Carbonate Reservoirs: Porosity, Evolution and Diagenesis in a Sequence Stratigraphic Framework[END_REF].

The micrite matrix surrounding the bioclasts record T 47 of 45 ± 3°C, significantly above the 15 to 27°C surface seawater temperature expected for the Middle Jurassic time [START_REF] Picard | δ 18 O values of coexisting brachiopods and fish: temperature differences and estimates of paleo-water depths[END_REF][START_REF] Lécuyer | Carbonate Reservoirs: Porosity, Evolution and Diagenesis in a Sequence Stratigraphic Framework[END_REF]. These high temperatures and the micrite negative  18 O carb values (-4.09 ± 0.61‰that are too depleted in 18 O compared to Middle Jurassic marine carbonates of [START_REF] Veizer | 87 Sr/ 86 Sr, δ 13 C and δ 18 O evolution of Phanerozoic seawater[END_REF] could be inherited from micrites carbonate mud recrystallization (for example, mineralogical stabilization) occurring under shallow burial conditions (as also suggested by [START_REF] Vincent | Geodynamic control on carbonate diagenesis: petrographic and isotopic investigation of the Upper Jurassic formations of the Paris Basin (France)[END_REF][START_REF] Carpentier | Origin of micropores in Late Jurassic (Oxfordian) micrites of the Eastern Paris Basin, France[END_REF] for Upper Jurassic micrites in the eastern Paris Basin). An alternative hypothesis lies on the non-straightforward interpretation of  47 values from bulk samples of very fine-grained sediments, because they could reflect the natural mixing of micrite with differently sourced material, such as detrital carbonates, and early to late diagenetic cements (refer to Defliese & Lohman, 2015, for mixing effects). In spite of these uncertainties, the recorded T 47 of 45°C dictate that at least some components of the micrite formed at temperatures higher than depositional conditions. This article is protected by copyright. All rights reserved.

Burial diagenetic carbonates

Clumped isotope temperature data on 21 Cal1 cement specimens, collected from three different cores and in a wide variety of pore types and depositional environments, indicate a very homogeneous precipitation event, occurring at the basin scale under burial temperatures of 64 ± 4°C (Fig. 7). The homogeneity of  13 C carb and  18 O carb values, observed in Cal1 sub-samples over 2 cm across a representative Cal1 cement specimen (cf. Mangenot et al., in press), suggests uniform conditions for precipitation. Furthermore, the large volume of Cal1 cements precipitated in the different cores indicates a massive and pervasive fluid-flow through the porous lithologies of the studied reservoir unit, mainly represented by the grain-supported facies of the Comblanchian Formation. This flow and consequent Cal1 cementation induced a significant porosity loss. This points to a good regional connectivity and an interconnected pore network, possibly controlled by the more porous intertidal facies, between the different wells, and thus throughout the central part of the Paris Basin. Compared to Cal1 cements, the four Dol1 cement samples record drastic changes in the CL pattern (Fig. 4B), more negative  18 O carb compositions, and the highest T 47 of the whole paragenetic sequence (i.e. 90 ± 6°C). Such results conclusively indicate that this saddle dolomite cement was precipitated during further burial and by hotter fluids. The occurrence of possibly primary oil inclusions in some Dol1 samples (Fig. 3C) leads to the conclusion that hydrocarbons migrated synchronously with the dolomite parent fluid, at least at the Villeperdue and Fossoy core localities. The high T 47 and the oil inclusions lead to the conclusion that precipitation of Dol1 occurred close to maximum burial conditions (as deduced from [START_REF] Espitalie | Aspects of hydrocarbon migration in the Mesozoic in the Paris Basin as deduced from an organic geochemical survey[END_REF][START_REF] Guilhaumou | Détermination de paléotempératures dans les roches-mères du bassin de Paris: etude d'inclusions fluides et implications pour l'histoire thermique du bassin[END_REF][START_REF] Gaulier | Modeling present and past thermal regimes in the Paris Basinpetroleum implications[END_REF][START_REF] Clauer | Clay minerals as records of temperature conditions and duration of thermal anomalies in the Paris Basin, France[END_REF][START_REF] Worden | Cross formational flow in the Paris Basin[END_REF][START_REF] Spötl | Clay minerals as records of temperature conditions and duration of thermal anomalies in the Paris Basin, France[END_REF][START_REF] Gonçalvès | Fluid inclusions as constraints in a three-dimensional hydro-thermo-mechanical model of the Paris basin, France[END_REF], also broadly corresponding to the main hydrocarbon migration phase from the source rocks (as deduced from [START_REF] Tremolières | Mécanismes de la deformation en zones de plate-frome. Méthode et application au Bassin de Paris[END_REF][START_REF] Poulet | Hydrocarbon migration in the Paris Basin. In: Migration of Hydrocarbons in Sedimentary Basins[END_REF][START_REF] Poulet | Hydrocarbon migration in the Paris Basin. In: Migration of Hydrocarbons in Sedimentary Basins[END_REF].

This article is protected by copyright. All rights reserved.

The two last cements of the paragenetic sequence, Cal2 and Dol2, exhibit both negative  18 O carb values and cooler T 47 of 78 ± 4°C and 67 ± 3°C, respectively. This cooling trend is interpreted to be related to the Tertiary uplift, which is extensively documented in the Paris Basin [START_REF] Guillocheau | Meso-Cenozoic geodynamic evolution of the Paris Basin: 3D stratigraphic constraints[END_REF][START_REF] André | Chronology of fracture sealing under a meteoric fluid environment: microtectonic and isotopic evidence of major Ceinozoic events in the eastern Paris Basin (France)[END_REF][START_REF] Barbarand | Lower Paleogene denudation of Upper Cretaceous cover of the Morvan Massif and southeastern Paris Basin (France) revealed by AFT thermochronology and constrained by stratigraphy and paleosurfaces[END_REF]. Because these two carbonate phases are volumetrically minor, the corresponding cementation events did not play an important role in the porosity reduction of the studied Middle Jurassic reservoir units.

Timing and duration of carbonate cementation

In the studied area of the Paris Basin depocentre, the thermal history of the Middle Jurassic reservoirs was already modelled (Uriarte, 1997) using BasinMod ® 1D software and a palaeo-heat flow of 62MW•m 2 at the Ambreville well locality (Figs 1E,1F and 11). This model was constructed to force the best compromise between stratigraphic data and palaeotemperatures issued from different palaeothermometers, such as vitrinite reflectance, fluid inclusions microthermometry, Rock-Eval pyrolysis and clay mineralogy. The key features of the modelled thermal history are the increase in heating from deposition until the beginning of the Tertiary, during the subsidencedominated stage of the basin, followed by a continuous cooling during the main uplift stage and consequent erosion [START_REF] Uriarte | Maturité thermique des sédiments de la bordure sud-est du Bassinde Paris[END_REF]. Because the Ambreville well is geographically close to the cores studied here (ca 50 to 80 km; Fig. 1E andF), located at the present-day burial depth of 1500 to 1900 m and falling within the 435 to 445°C T max maturity zone (Fig. 1E andF; [START_REF] Delmas | Paris Basin. Rapport r gional d' valuation p trolière[END_REF], it may be considered that this well experienced burial and thermal histories similar to the four cores investigated in this study.

This article is protected by copyright. All rights reserved.

The thermal history reconstructed by [START_REF] Uriarte | Maturité thermique des sédiments de la bordure sud-est du Bassinde Paris[END_REF] for the basin depocentral area can thus be used to convert the calculated T 47 into relative dates for the different carbonate cementation events, assuming a thermal equilibrium between the cement parent fluid and the host rock. The average TΔ 47 defined for each cement phase yields uncertainties for the average ages of about ±10 Ma.

Whilst the exact cementation ages may be biased, they have been derived using the same methodology and are therefore comparable. Cementation ages of 135 to 110 Ma (Lower Cretaceous) for Cal1, 75 to 65 Ma (Upper Cretaceous) for Dol1, 60 to 55 Ma (Palaeocene-Eocene) for Cal2 and 45 to 30 Ma (Eocene-Oligocene) for Dol2 (Fig. 11) were thus assumed. The ages obtained by this way confirm the paragenetic succession proposed for the different cementation events, as determined by petrographic analysis only (Figs 5 and 10). The time-temperature information derived for each carbonate cement (Fig. 11) also confirms that two of these cements (Cal1 and Dol1) precipitated during the main Cretaceous subsidence stage, whereas the two others (Cal2 and Dol2) precipitated during the main Tertiary uplift stage.

Carbonate cementation has been interpreted both as a prolonged process, that occurs slowly over a long span of time and as a rapid process, in which short bursts of activity are followed by long periods of quiescence [START_REF] Tucker | Carbonate Sedimentology[END_REF]Moore, 2001). Here, the inferred cementation ages (Fig. 11) indicate that cementation occurred as a rapid (<10 Ma) and episodic process, in response of punctual bursts of fluid-flow activity. With the present knowledge, it is however impossible to characterize precise cementation rates.

Parent fluid geochemistry and evolution

Coupling the independent temperature estimates out of  47 data and the carbonate oxygen isotope values ( 18 O carb ) allows a direct calculation of the oxygen isotope composition of the fluid ( 18 O water ) from which the carbonate formed. This provides straightforward constraints on water-rock This article is protected by copyright. All rights reserved.

interactions that are not achievable using  18 O carb values only. Furthermore, when TΔ 47 are interpreted together with fluid inclusions salinity, they can provide even better constraints on the isotope and compositional evolution of subsurface water-rock systems [see Mangenot et al. (in press) for a practical description of the approach], and therefore yield insights into regional groundwater fluid-flows. The fluid inclusions dataset produced from five samples of Cal1, Dol1 and Cal2 (Mangenot et al., in press; see synthesis in Fig. 12) are also integrated and discussed here.

Marine and early diagenetic pore-water geochemistry

In the studied Middle Jurassic carbonates, the earliest phases investigated are syn-depositional marine in origin; these are biogenic carbonates (from brachiopod and bivalve shells) and the surrounding micrites. Both carbonate groups display similar  8 O water values both averaging 1.8‰; Table 1 and Fig. 8), that are in agreement with the expected isotopic composition of Jurassic seawater in restricted lagoon settings [START_REF] Picard | δ 18 O values of coexisting brachiopods and fish: temperature differences and estimates of paleo-water depths[END_REF][START_REF] Lécuyer | Carbonate Reservoirs: Porosity, Evolution and Diagenesis in a Sequence Stratigraphic Framework[END_REF]. When interpreted together with the TΔ 47 data (Fig. 7), the  8 O water values indicate that the marine bioclast carbonates possibly derived from Middle Jurassic lagoon seawater at 29 ± 4°C, whereas the micrites seem to have been recrystallized at about 45 ± 3°C during burial (Fig. 11) from primary marine fluids trapped within the initial micrite pores (for example, connate waters).

Basinal brines associated with mesogenetic cements

The Cal1 mineralizing fluids are enriched in both 18 O and salt (with average δ 18 O water of 2.5‰ and salinity mode of 14 wt% NaCl eq; Figs 8 and12) compared to modern seawater (δ 18 O water between -1‰ and +1‰ and salinities of 0.5 to 5 wt% NaCl eq; [START_REF] Rohling | Oxygen isotope composition of seawater[END_REF]. Such features conjointly argue for a basinal brine origin. The Dol1 cement recorded δ 18 O water (Fig. 8) close to marine values (δ 18 O water This article is protected by copyright. All rights reserved. of 0.5‰), while fluid inclusion salinities are about three times higher than seawater (from 10.1 to 11.5 wt% NaCl eq; Fig. 12). Such parent fluid compositions could be explained by three different scenarios: (i) a primary brine (evaporated seawater); (ii) a normal seawater which got enriched in salt and 18 O by water-evaporite interaction during burial; and (iii) emplacement of an externally sourced brine circulating through the sedimentary succession during burial. The first two hypotheses can be rejected based on the accepted sedimentological models for the studied Middle Jurassic reservoirs, which support the lack of both evaporitic seawater conditions and massive evaporite bodies [START_REF] Purser | Plateformes carbonatées: exemple du Jurassique Moyen du Bassin de Paris. Dynamique et méthodes d'étude des bassins sédimentaires[END_REF][START_REF] Javaux | La plateforme parisienne et bourguignonne au Bathonien terminal et au Callovien[END_REF][START_REF] Gaumet | Contrôle stratigraphique des faciès, géométries et profils de dépôt de la plate-forme carbonatée bourguignonne au Bathonien-Callovien[END_REF]. The data thus support the third scenario, implying the flow of externally sourced fluids through the Middle Jurassic reservoirs.

External sources of salty water in the central part of the Paris Basin could be the Upper Triassic aquifers (i.e. Chaunoy and Rhetien sandstone reservoirs) that display present-day salinity up to 20 wt% NaCl eq and  18 O water up to 1.1‰ [START_REF] Worden | Cross formational flow in the Paris Basin[END_REF][START_REF] Worden | Basin scale evolution of formation waters: a diagenetic and formation water study of the Triassic Chaunoy Formation, Paris Basin[END_REF]Millot et al., 2011). Furthermore, the occurrence of a cross-formational flow between Triassic and Middle Jurassic aquifers has already been proposed in diagenetic models from previous authors (Matray et al., 1994;[START_REF] Worden | Cross formational flow in the Paris Basin[END_REF][START_REF] Bril | Effects of the 'Pays De Bray' fault on fluid paleocirculations in the Paris Basin Dogger Reservoir[END_REF]. Most importantly, this hypothesis is also supported by the geochemistry dataset for Cal1 and Dol1 revealing parent fluid salinity and δ 18 O water compositions similar to the present-day groundwater of the Upper Triassic aquifers (Fig. 13; Worden & Maray, 1995;[START_REF] Worden | Basin scale evolution of formation waters: a diagenetic and formation water study of the Triassic Chaunoy Formation, Paris Basin[END_REF]. Further geochemical investigations ( 3 He, 34 S, 14 C and trace elements geochemistry) conclusively showed the occurrence of a Triassic saline component within the present-day Middle Jurassic aquifers (Marty et al., 1988;[START_REF] Fouillac | Sulphur and oxygen isotopes of dissolved sulphur species in formation waters from the Dogger geothermal aquifer, Paris Basin, France[END_REF][START_REF] Wei | Regional modelling of groundwater flow and salt and environmental tracer transport in deep aquifers in the Paris Basin[END_REF]Marty et al., 1993;[START_REF] Pinti | Noble gases in crude oils from the Paris Basin, France: implications for the origin of fluids and constraints on oil-water-gas interactions[END_REF][START_REF] Pinti | Atmosphere-derived noble gas evidence for the preservation of ancient waters in sedimentary basins[END_REF]. All of the studied subsurface cores exhibit voluminous Cal11 and Dol1 cements, with similar petrographic and geochemical features. This implies that this upward fluid migration, leading to the invasion of the Middle Jurassic reservoirs by hypersaline brines, was a widespread process that occurred at basin scale. Fluids with similar or This article is protected by copyright. All rights reserved. even higher salinities than those recorded in Cal1 and Dol1 (up to 25 wt% NaCl eq; Fig. 9B.1) were indeed recorded from other cores of the basin depocentre (Matray et al., 1989;[START_REF] Demars | volu on diag n ue, pal o uides et pal othermicit dans les r servoirs du Keuper et du Dogger du bassin de Paris[END_REF][START_REF] Worden | Cross formational flow in the Paris Basin[END_REF], allowing the extrapolation of the dataset here presented at larger scale. Some Dol1 crystals contain hydrocarbon bearing fluid inclusions (Fig. 4C), suggesting that oil was trapped at the time of Dol1 precipitation. Such hydrocarbon bearing fluid inclusions were not reported in previous diagenetic studies of the Paris Basin and attest that Dol1 parent fluids and oil have migrated together through the Middle Jurassic reservoirs. Because the most probable source rocks for these hydrocarbons are the underlying Lower Jurassic shales [START_REF] Poulet | Hydrocarbon migration in the Paris Basin. In: Migration of Hydrocarbons in Sedimentary Basins[END_REF][START_REF] Espitalie | Aspects of hydrocarbon migration in the Mesozoic in the Paris Basin as deduced from an organic geochemical survey[END_REF], the dolomitizing brine probably passed through these source rocks during the upward cross-formational flow described above. This hypothesis also satisfies the supposed mineralization age of Dol1 at 75 to 65 Ma (Fig. 11), which corresponds to the timing of peak-oil generation for these source rocks (90 to 70 Ma; Monticone et al., 2011).

Meteoric freshening associated with telogenetic cements

The negative δ 18 O water values determined for Cal2 and Dol2 parent fluids (averaging at -4.1‰ for Cal2 and at -7‰ for Dol2; Fig. 8), the salinity dilution pattern observed in Cal2 fluid inclusions (Fig. 12), as well as the different estimated ages (Fig. 11), suggest that two distinct pulses of diluted meteoric waters occurred through the Middle Jurassic reservoirs during the Tertiary uplift of the basin. These diluted fluid pulses led to the progressive freshening of the pre-existing brines described above (refer to Basinal brines associated with mesogenetic cements section above). Fluids with similar low salinity (Fig. 9B.1) were indeed recorded in other cores of the basin depocentre [START_REF] Demars | volu on diag n ue, pal o uides et pal othermicit dans les r servoirs du Keuper et du Dogger du bassin de Paris[END_REF][START_REF] Worden | Cross formational flow in the Paris Basin[END_REF] indicating that the influx of meteoric recharge fluids was not a local process.

This article is protected by copyright. All rights reserved.

These meteoric freshening pulses recorded by the last two cements of the paragenesis (i.e. the telogenetic Cal2 and Dol1) would explain why the present-day formation water in the Middle Jurassic reservoirs displays a very low saline load (around 3% TDS) and a negative δ 18 O water value (Fig. 13). Lateral recharges of meteoric waters were previously documented from the study of diagenetic cements in regions closer to the eastern margin of the basin [START_REF] Buschaert | Widespread cementation induced by inflow of continental water in the eastern part of the Paris Basin: O and C isotopic study of carbonate cements[END_REF][START_REF] Vincent | Geodynamic control on carbonate diagenesis: petrographic and isotopic investigation of the Upper Jurassic formations of the Paris Basin (France)[END_REF][START_REF] Brigaud | The origin and timing of multiphase cementation in carbonates: impact of regional scale geodynamic events on the Middle Jurassic Limestones diagenesis (Paris Basin, France)[END_REF][START_REF] André | Chronology of fracture sealing under a meteoric fluid environment: microtectonic and isotopic evidence of major Ceinozoic events in the eastern Paris Basin (France)[END_REF][START_REF] Carpentier | Impact of basin burial and exhumation on Jurassic carbonates diagenesis on both sides of a thick clay This article is protected by copyright. All rights reserved. barrier[END_REF]). In the central part of the basin, meteoric water infiltrations were already suggested by hydrogeological modelling studies [START_REF] Wei | Regional modelling of groundwater flow and salt and environmental tracer transport in deep aquifers in the Paris Basin[END_REF][START_REF] Gonçalvès | 3-D modelling of salt and heat transport during the 248 m.y. evolution of the Paris Basin; diagenetic implications[END_REF][START_REF] Violette | Reactive transport modelling of carbonate cementation in a deep saline aquifer, the Middle Jurassic isotopologues of molecular gases[END_REF], chlorine concentrations, hydrogen and oxygen isotope compositions (Matray et al., 1994) and noble gas geochemistry [START_REF] Castro | Noble gases as natural tracers of water circulation in the Paris Basin: 2. Calibration of a groundwater flow model using noble gas isotope data[END_REF]. However, such multiple meteoric water recharges were so far not being assessed from carbonate diagenetic cements of the central part of the basin. This meteoric water, initially cold at the recharge area (i.e. surface temperature), must have re-equilibrated with the ambient rock temperatures during circulation at depth to explain the palaeotemperatures of 78 ± 4°C and 67 ± 3°C recorded in Cal2 and Dol2 cements, respectively. In such burial conditions, the preservation of This article is protected by copyright. All rights reserved.

observations suggest that periods of active tectonics in the basin favoured fluid-flows, leading to the precipitation of carbonate cements in the remaining pore space of the Middle Jurassic reservoirs.

This geodynamic control on cementation was previously demonstrated in the Upper Jurassic carbonates of the eastern Paris Basin, located at ca 500 km from the studied area, which were mainly affected by meteoric cementation [START_REF] Buschaert | Widespread cementation induced by inflow of continental water in the eastern part of the Paris Basin: O and C isotopic study of carbonate cements[END_REF][START_REF] Gonçalvès | Contribution of a three-dimensional regional scale basin model to the study of the past fluid flow evolution and the present hydrology of the Paris basin, France[END_REF][START_REF] Vincent | Geodynamic control on carbonate diagenesis: petrographic and isotopic investigation of the Upper Jurassic formations of the Paris Basin (France)[END_REF][START_REF] Brigaud | The origin and timing of multiphase cementation in carbonates: impact of regional scale geodynamic events on the Middle Jurassic Limestones diagenesis (Paris Basin, France)[END_REF][START_REF] André | Chronology of fracture sealing under a meteoric fluid environment: microtectonic and isotopic evidence of major Ceinozoic events in the eastern Paris Basin (France)[END_REF][START_REF] Carpentier | Impact of basin burial and exhumation on Jurassic carbonates diagenesis on both sides of a thick clay This article is protected by copyright. All rights reserved. barrier[END_REF][START_REF] Blaise | Evolution of porewater composition through time in limestone aquifers: Salinity and D/H of fluid inclusion water in authigenic minerals (Jurassic of the eastern Paris Basin, France)[END_REF].

Previous hydrogeological models agreed with the occurrence of a brine migration from the Upper Triassic to the Middle Jurassic aquifers of the studied depocentre area, though the relative timing of these upward fluid-flows is still debated. Early studies (Matray et al., 1989) suggested an Early to Late Cretaceous timing for this cross-formational fluid-flow, whereas more recent ones [START_REF] Worden | Cross formational flow in the Paris Basin[END_REF][START_REF] Bril | Effects of the 'Pays De Bray' fault on fluid paleocirculations in the Paris Basin Dogger Reservoir[END_REF][START_REF] Gonçalvès | 3-D modelling of salt and heat transport during the 248 m.y. evolution of the Paris Basin; diagenetic implications[END_REF][START_REF] Violette | Reactive transport modelling of carbonate cementation in a deep saline aquifer, the Middle Jurassic isotopologues of molecular gases[END_REF] suggested that it was a Cenozoic phenomenon, occurring during the compressional tectonic phase of the Alpine Orogeny. This article is protected by copyright. All rights reserved.

of the proposed fluid-flow scenario. The Cretaceous timing of brine migration, inferred for Cal1 and Dol1 parent fluids in the central part of the basin, can also be temporally related to carbonate cementation events previously identified in the eastern and southern margins of the Paris Basin [START_REF] Brigaud | The origin and timing of multiphase cementation in carbonates: impact of regional scale geodynamic events on the Middle Jurassic Limestones diagenesis (Paris Basin, France)[END_REF][START_REF] Gigoux | Evidence for an Early Cretaceous mineralizing event above the basement/sediment unconformity in the intracratonic Paris Basin: paragenetic sequence and Sm-Nd dating of the world-class Pierre-Perthuis stratabound fluorite deposit[END_REF], as well as in the Wessex Basin of southern England [START_REF] Hendry | Geological control on regional subsurface carbonate cementation: an isotopic paleohydrologic investigation of middle Jurassic limestones in Central England. In: Diagenesis and Basin Development[END_REF]Worden et al., 2015), the latter being the extension of the Paris Basin to the North.

Indeed, [START_REF] Gigoux | Evidence for an Early Cretaceous mineralizing event above the basement/sediment unconformity in the intracratonic Paris Basin: paragenetic sequence and Sm-Nd dating of the world-class Pierre-Perthuis stratabound fluorite deposit[END_REF] reported a major fluid-flow leading to fluorine mineralization in the southern Paris Basin at 130 ± 15Ma (from Sm-Nd dating), whereas [START_REF] Brigaud | The origin and timing of multiphase cementation in carbonates: impact of regional scale geodynamic events on the Middle Jurassic Limestones diagenesis (Paris Basin, France)[END_REF] reported a contemporaneous meteoric water input, leading to carbonate cements, at the Aptian-Albian transition in Upper Jurassic strata of the eastern margin of the basin. Cretaceous fluid-flows with meteoric affinity have also been documented in the Middle Jurassic strata of the Wessex Basin in southern England [START_REF] Hendry | Geological control on regional subsurface carbonate cementation: an isotopic paleohydrologic investigation of middle Jurassic limestones in Central England. In: Diagenesis and Basin Development[END_REF]Worden et al., 2015). When all these studies are integrated together, they provide strong evidence of multiple and tectonically controlled fluid-flows during Cretaceous time, that led to massive movements of differently sourced waters through the different carbonate units of the Paris Basin.

Successively, tectonic inversion occurred in the Paris Basin, mainly linked to the beginning of the Pyrenean compression at the Late Cretaceous to Paleogene boundary [START_REF] Barbarand | Lower Paleogene denudation of Upper Cretaceous cover of the Morvan Massif and southeastern Paris Basin (France) revealed by AFT thermochronology and constrained by stratigraphy and paleosurfaces[END_REF].

These geodynamic changes induced the exposure of the basin margins and a drastic modification in the fluid-flow mechanisms that would have favoured lateral infiltrations of surface derived waters at depth (i.e. meteoric recharge; see for instance [START_REF] Garven | Theoretical analysis of the role of groundwater flow in the genesis of strata bound ore deposits[END_REF]. The findings here agree with this hypothesis and demonstrate that Cal2 and Dol2 (the last two cement phases of the paragenesis)

have been precipitated by meteoric-derived fluids during Palaeocene-Eocene (60 to 55 Ma) and

Eocene-Oligocene (45 to 30 Ma) times, respectively. The mixing of these meteoric waters with the pre-existing hypersaline groundwater implies a topography-driven fluid-flow, with a sufficiently high hydraulic gradient, and a well-interconnected pore network to channel the fluids from surface to depth (more than 1500 m) along several hundreds of kilometres. [START_REF] Dorobek | Petrography, geochemistry, and origin of burial diagenetic facies, Siluro-Devonian Helderberg Group (carbonate rocks), Central Appalachians[END_REF] This article is protected by copyright. All rights reserved.

away from the inferred recharge area, in tectonically uplifted basin margins, buried beneath some 250 m of low permeable strata. The broadly similar scale of this process suggests that in the Paris Basin, Tertiary telogenetic meteoric waters would have been capable of rapidly extending through the Middle Jurassic strata at basin scale, to reach the deepest parts of the basin center. However, the results of isotope analysis on present-day formation waters (Fig. 13) show that the Middle Jurassic aquifers were not completely flushed away by the infiltrated Tertiary meteoric waters (Matray et al., 1994). This suggests that the hydraulic head, triggering the meteoric recharge fluidflow, was not sufficiently high to cause the complete flushing of the pre-existing brines, possibly due to the lack of a high-relief highland.

To conclude, when the whole Paris Basin hydrogeological scenario is considered, it shares numerous similarities with the fluid-flow histories recorded in other intracratonic basins with comparable architecture, such as the Illinois Basin [START_REF] Stueber | Glacial recharge and paleohydrologic flow systems in the Illinois Basin: evidence from chemistry of Ordovician carbonate (Galena) formation waters[END_REF][START_REF] Labotka | Isotopic and geochemical characterization of fossil brines of the Cambrian Mt. Simon Sandstone and Ironton-Galesville Formation from the Illinois Basin, USA[END_REF] or the Michigan Basin [START_REF] Wilson | Geochemistry and isotope chemistry of CaNaCl brines in Silurian strata, Michigan Basin[END_REF]. Indeed, these basins also experienced late stage infiltration of meteoric fluids into more ancient saline evaporative waters. Therefore, the combination of clumped isotopes thermometry with conventional petrography and fluid inclusion studies on carbonate burial cements could serve as a powerful investigation approach to achieve a more detailed reconstruction of the hydrogeological histories and water-rock interactions in other sedimentary basins.

CONCLUSIONS

Middle Jurassic carbonates from the central part of the Paris Basin were investigated based on a sample collection from four sub-surface cores (present depth of 1500 to 1900 m). The thermal and diagenetic evolution of this stratigraphic interval was precisely reconstructed using petrographic observations combined with extensive  47 analyses for each identified carbonate phase.

This article is protected by copyright. All rights reserved.

The implications of the collected data can be summarized as follows:

 Sedimentation occurred from slightly evaporated seawater in a lagoonal setting ( 18 O water of 1.8 ± 1.4‰) at temperatures of 29 ± 4°C.

 Syn-sedimentary micrites record recrystallization during burial at higher temperatures (close to 45°C) from connate waters also enriched in 18 O ( 18 O water of 2.3 ± 0.7‰).

 A first major pore-filling calcite cementation phase (Cal1) precipitated at 64 ± 4°C, from high salinity fluids (14 wt.% NaCl eq) with a  18 O water close to 2.5‰, pointing to a primary brine.

 A second cementation phase induced dolomite (Dol1) precipitation at 90 ± 6°C, from similar fluids with moderate salinity (11 wt.% NaCl eq) and a  18 O water value of 0.5 ± 0.4‰, coeval with hydrocarbon reservoir charging.

 A third cementation phase resulted in calcite precipitation (Cal2) at 78 ± 4°C associated with the circulation of meteoric derived fluids, with low salinity (from 12 to 0 wt.% NaCl eq) and a

 18 O water value of -4.1 ± 0.5‰. This article is protected by copyright. All rights reserved.

Overall, the results demonstrate how a careful petrographic characterization is an important prerequisite for identifying individual carbonate cementation phases at basin scale, and to apply the carbonate clumped isotope (Δ 47 ) thermometry at its highest level of performance (in terms of both accuracy and precision on temperature and  18 O water of the carbonate parent fluids). Such study permits to refine constraints on the fluid-flow and water-rock system evolution at the basin scale.

More broadly, this study also illustrates how Δ 47 thermometry provides information of primary importance to constrain carbonate diagenesis that largely extends what is currently achievable with conventional isotopic and petrographic data alone. The use of Δ 47 thermometry on diagenetic cements is a powerful tool to investigate fluid-flows and palaeo-thermicity, thereby addressing a longstanding conundrum in basin research. Jurassic carbonates (modified after [START_REF] Delmas | Paris Basin. Rapport r gional d' valuation p trolière[END_REF]. (F) Map of Rock-Eval T max from pyrolysis of organic matter (OM) of the Lower Toarcian source rocks (modified after [START_REF] Delmas | Paris Basin. Rapport r gional d' valuation p trolière[END_REF]. In (E) and (F), the location of the studied reservoirs is shown with grey dots whereas the white square represents a well used for thermal history modelling (Uriatre, 1997). This article is protected by copyright. All rights reserved. et al., 1989;[START_REF] Demars | volu on diag n ue, pal o uides et pal othermicit dans les r servoirs du Keuper et du Dogger du bassin de Paris[END_REF][START_REF] Worden | Cross formational flow in the Paris Basin[END_REF]. (A.2 and B.2)

FIGURE CAPTIONS

Histograms of homogenization temperatures and salinities acquired for the investigated cores after careful petrography and phase differentiation prior to analysis (Mangenot et al., in press). Analyzed

This article is protected by copyright. All rights reserved.

phases are from: BEBJ12 and BEBJ8 samples (Cal1, Baulne en Brie), FOS1610 sample (Dol1, Fossoy), and BEBJ2 and BEBJ samples (Cal2, Baulne en Brie).

  (average  18 O carb values varying by more than 13‰) although with homogeneous  18 O carb values when only one carbonate phase is considered (with one standard deviation values averaging at ±0.5‰). As such, the average  18 O carb values are clustered at  18 O carb = -1.59 ± 0.97‰ (n = 3 samples) for the bioclasts;  18 O carb = -4.09 ± 0.61‰ (n = 4 samples) for the micrites;  18 O carb = -6.97 ± 0.72‰ (n = 21 samples) for Cal1;  18 O carb = -9.23 ± 0.39‰ (n = 4 samples) for Dol1;  18 O carb = -15.27 ± 0.49‰ (n = 11 samples) for Cal2; and  18 O carb = 13.50 ± 0.37‰ (n = 2 samples) for Dol2. The  18 O carb values of This article is protected by copyright. All rights reserved.

  are distinguished from the most recent data (from Mangenot et al., in press), acquired on five samples from the Baulne en Brie and Fossoy cores. The former authors did not differentiate the various carbonate phases bearing the measured fluid inclusions, thus their microthermometric data probably record the mixing of different generations of fluid inclusions (Fig. 9A.1 and 9B.1). Conversely, prior to microthermometry analysis, Mangenot et al. (in press) differentiated three geochemically and petrographically different carbonate phases (transmitted light, UV light and CL), which correspond to Cal1 (two samples) Dol1 (one sample) and Cal2 (two samples) of the present study (Fig. 9A.2 and 9B.2).

2

  The high resolution microthermometry data fromMangenot et al. (in press), which globally overlap the data of previous authors on petrographically undifferentiated carbonate cements from several sub-surface cores, allow discussion of the diagenetic evolution of the Middle Jurassic reservoirs at a larger scale. The good consistency between the thermal information acquired by fluid inclusion and clumped isotope thermometry on the same samples of Cal1, Dol1 and Cal2 (five samples in total) in the 60 to 100°C temperature range demonstrates the accuracy of the TΔ 47 estimates, which were achieved in this study on additional 35 samples collected from the same carbonate rock unit.

  and -7.0‰ (typical of meteoric fluids in the Paris Basin;Matray et al. 1994), could reflect high water-rock ratios, pointing to an open system behaviour (as also suggested by[START_REF] Lavastre | Water circulation control on carbonate-δ 18 O records in a low permeability clay formation and surrounding limestones: the Upper Dogger-Oxfordian sequence from the eastern Paris Basin, France[END_REF] in the eastern part of the Paris Basin). In such context, the factors governing the switch from calcite to dolomite in the cementation sequence (Cal1-Dol1-Cal2-Dol2) cannot be easily explained.Geodynamic control on fluid-flowsIn the central part of the Paris Basin, the four cementation dates estimated above (135 to 110 Ma for Cal1, 75 to 65 Ma for Dol1, 60 to 55 Ma for Cal2 and 45 to 30 Ma for Dol2; Fig.11) correspond to well-documented periods of active tectonic stresses (i.e. Late Aptian unconformity, Late Cretaceous unconformity, Palaeocene north-south compression and Oligocene extension; Fig.11). Such

  The novel contribution of this work, by giving access to the temperature,  18 O water , salinity and relative timing of Cal1 and Dol1 brine migrations(Figs 7, 8 and 11 to 13), largely supports the first hypothesis(Matray et al., 1989). Indeed, the findings of this study reveal that the saline load of the Middle Jurassic aquifers was due to the externally sourced parent fluids of Cal1 and Dol1, which precipitated respectively during Early Cretaceous (135 to 110 Ma) and Late Cretaceous (75 to 65 Ma) time, in a dominantly subsiding context. During subsidence stages of basin evolution, pore fluid pressure tends to increase and may favour upward fluid-flows. Permeable fault zones may then act as drains for fluids to move. In this respect, Cal1 and Dol1 parent fluids could have been channeled from the Upper Triassic donor aquifers towards the Middle Jurassic receiver aquifers by such upward flow along faults (Fig.13B). Interestingly,Guillaucheau et al. (2000) have demonstrated that the Variscan Bray fault system (Fig.1C), cross-cutting the basin depocentre, experienced different episodes of reactivation during Lower and Upper Cretaceous times, largely supporting the reliability



  A fourth and last cementation phase led to the formation of dolomite (Dol2) at temperatures of 67 ± 3°C by meteoric derived fluids with a negative  18 O water value of -7.0 ± 0.8‰. The thermal information and fluid properties obtained by Δ 47 thermometry is consistent with the burial, hydrogeological and tectonic evolution known for the Paris Basin and allows refinement of previous diagenetic models. Findings in this study imply two major brine circulation episodes in the Early (Cal1) and Late (Dol1) Cretaceous times, during the main subsidence phase, followed by two diluted water influxes (Cal2 and Dol2), during the Tertiary basin inversion and consequent uplift. The  18 O water and salinity for Cal1 and Dol1 parent fluids indicate a likely contribution of Triassic reservoir fluids that reached the Middle Jurassic lithologies through a cross-formational flow along major faults, such as the Bray fault system. Conversely,  18 O water and salinity for Cal2 and Dol2 parent fluids indicate a freshening of the Middle Jurassic reservoirs by a topography-driven flow during the Tertiary, from an unknown meteoric recharge area.

Fig. 1 .

 1 Fig. 1. (A) Synthetic stratigraphic column of the Paris basin with location of the studied stratigraphic

Fig. 2 .

 2 Fig. 2. Macrophotography, thin-section scans and microphotography (transmitted light and CL) of

Fig. 3 .

 3 Fig. 3. Macrophotography, thin-section scans and microphotography (transmitted light and CL)

Fig. 5 .

 5 Fig. 5. Paragenetic sequence, composed of 14 diagenetic events, reconstructed for the Middle

Fig. 6 .

 6 Fig. 6. Stable oxygen and carbon isotope composition of the sedimentary and diagenetic carbonates

Fig. 9 .

 9 Fig. 9. Compilation of fluid inclusions microthermometry results for carbonate cements from Middle

Fig. 10 .

 10 Fig.10. (A) Schematic drawing illustrating the diagenetic evolution of the studied Middle Jurassic

Fig. 13 .

 13 Fig. 13. Conceptual diagram showing the relative fluid-flows chronology, deduced from carbonate cements geochemistry. On the left is illustrated the salinity versus  18 O water cross-plot, summarizing

Stable isotopes data ( 18 O carb ,  13 C carb ,  47 ,  18 O water )

  Table 1 reports the average  18 O carb ,  13 C carb ,  47 , T 47 and  18 O water values determined from two to three replicate isotopic measurements of the same powder for each investigated carbonate sample, belonging to six genetically different carbonate phases. Detailed information on each single measurement, correction and standardization can be found in the supplementary material of Appendix S1 (Table

  described Silurian meteoric cements in basinal limestones of the central Appalachians precipitated at least 150 km

  symbol size. The dashed box represents  18 O and  13 C values of biogenic calcites precipitated in equilibrium with Middle Jurassic seawater (from[START_REF] Veizer | 87 Sr/ 86 Sr, δ 13 C and δ 18 O evolution of Phanerozoic seawater[END_REF]. Fig. 7. (A) Cross-plot of the carbonate precipitation temperatures versus conventional  18 O carb estimated from  47 thermometry. Reported uncertainties on each data point are one standard deviation of the mean for replicate measurements of the same carbonate powder for T 47 and are included in the symbol size for  18 O carb values (Table1). The long-term external  47 reproducibility measured on more than 50 carbonate standards over the course of this study (i.e. ±0.017‰) is also reported on the left bottom side for comparison. The six investigated carbonate phases are represented with different colours. Average temperatures calculated from all of the samples of each carbonate phase are also reported in bold with their respective one standard deviation of the mean (Table2). (B) Cross-plot of the carbonate precipitation temperatures estimated from  47 thermometry versus  18 O water values back-calculated from  47 and  18 O carb compositions. Fig. 8. Average oxygen isotope composition of the carbonate parent fluids ( 18 O water ) for each investigated carbonate phase (Table2). The length of the coloured rectangles represents one standard deviation of the mean  18 O water value calculated for each carbonate phase (reported as vertical black bars and on the left of each box, together with the number of samples considered).

Table 1 .

 1 Stable isotope ( 47 ,  18 O and  13 C) results from individual diagenetic cements.

Table 2 .

 2 Average composition for each cement population.

ACKNOWLEDGMENTS

We thank Carine Chaduteau and Amandine Katz for help with  47 measurements at IPGP. We are also grateful to Magali Ader, Pascal Houel, Geneviève Bessereau, Benjamin Brigaud, Jean-Pierre Girard and Christophe Durlet for the constructive discussions about the thermicity, the diagenesis and fluid-flow histories of the Paris Basin. The subsurface cores of Baulne en Brie, Rigny la Nonneuse and Fossoy were available from the IFPEN storage collection of the BEPH (Bureau Exploration-Production d'Hydrocarbures). The Villeperdue subsurface core was available from the University of Dijon. We are also grateful to Alain Buisson and Pierre André Duboin from Lundin Petroleum for sharing their knowledge on the petroleum reservoirs of the Paris Basin.

evidence from δ 18 O values of marine fish teeth. Paleoceanography, 18, 1076. This article is protected by copyright. All rights reserved. Loyd, S.J., Corsetti, F.A., Eiler, J.M. and Tripati, A.K. (2012) Determining the diagenetic conditions of concretion formation: assessing temperatures and pore waters using clumped isotopes. J. Sediment. Res, 82, 1006-1016. Loyd S. J., Corsetti F. A., Eagle R. A., Hagadorn J. W., Shen Y., Zhang X., Bonifacie M. and Tripati A. K. (2015) Evolution of Neoproterozoic Wonoka-Shuram anomaly-aged carbonates: evidence from clumped isotope paleothermometry. Precambrian Res, 264, 179-191. Mangenot, X., Bonifacie, M. Gasparrini, M., Goetz, A., Chaduteau, C., Ader, M., Rouchon, V. (in press) Coupling  47 and fluid inclusion thermometries on carbonate cements to precisely reconstruct the temperature, salinity and  18 O of circulating paleowater in sedimentary basins. Chemical Geology. Marty, B., Criaud, A. and Fouillac, C. (1988) Low enthalpy geothermal fluids from the Paris sedimentary basin-1. Characteristics and origin of gases. Geothermics, 17, 619-633. M T g s n T M ni V O'Ni ns R K and De Marsily, G. (1993) Helium isotope fluxes and groundwater ages in the Dogger Aquifer, Paris Basin. Water Resour. Res, 29, 1025-1035. Matray J.M., Meunier, A., Thomas, M. and Fontes, J.C. (1989) The Triassic and Dogger formation waters from the Paris Basin: history and diagenetic effects upon the reservoirs. Soc. Nat. Elf-Aquitaine -Internal Report, 13, 484-504. Matray, J.M., Lambert, M. and Fontes, J.C. (1994) Stable isotope conservation and origin of saline waters from the Middle Jurassic aquifer of the Paris Basin, France. Appl. Geochemistry, 9, 297-309. McCrea, J.M. (1950) On the isotopic chemistry of carbonates and a paleotemperature. The Journal of Chemical Physics, 18, 849 Ménétrier, C., Élie, M., Martinez, L. and Le, A. (2005) Estimation of the maximum burial palaeotemperature for Toarcian and Callovo-Oxfordian samples in the central part of the Paris Basin using organic markers. Comptes Rendus Geosci, 337, 1323-1330.

This article is protected by copyright. All rights reserved. Millot, R., Guerrot, C., Innocent, C., Philippe, N. and Sanjuan, B. (2011) Chemical , multi-isotopic (Li-B-Sr-U-O) and thermal characterization of Triassic formation waters from the Paris Basin. Chemical Geology, 283, 226-241. Monticone, B., Duval, M., Knispel, R., Wojciak, P. and Dubille, M. H. (2011) Shale oil potential of the Paris Basin, France. AAPG Int. Conf. Exhib, 10384.