# A Magnetic Geothermometer in Moderately Buried Shales

Charles Aubourg, Myriam Kars, Jean-Pierre Pozzi, Martin Mazurek, O. Grauby

► To cite this version:

Charles Aubourg, Myriam Kars, Jean-Pierre Pozzi, Martin Mazurek, O. Grauby. A Magnetic Geothermometer in Moderately Buried Shales. Minerals, 2021, 11, 10.3390/min11090957. [https://doi.org/10.3390/min11090957](https://doi.org/10.3390/min11090957)

HAL Id: hal-03629197

https://amu.hal.science/hal-03629197

Submitted on 4 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# A Magnetic Geothermometer in Moderately Buried Shales

Charles Aubourg <fnref n="1" />*, Myriam Kars <fnref n="2" />, Jean-Pierre Pozzi <fnref n="3" />, Martin Mazurek <fnref n="4" /> and Olivier Grauby <fnref n="5" />

**Abstract:** Shales contain magnetic minerals generally at very low concentrations. In the early stages of diagenesis, the inherited magnetic minerals are altered, while magnetic nanominerals are formed. In this study, we proposed a study of shales over a stratigraphic thickness of 1.3 km from a borehole in the Paris basin (Borehole EST 433, France), and shales from the same formation (Opalinus Clay) collected in seven boreholes in the Jura molasse basin (Swiss). Magnetic measurements at experimental temperatures <30 K allowed the formation of a proxy of magnetite nanograins named PM. We showed that some of these nanograins formed around the pyrite grains, probably under the action of temperature and organic matter. PM was then compared to the maturity values of the organic matter. We found a correlation between PM and the percentage of reflectance of vitrinite. The shales from both Paris and molassic Swiss basins showed very comparable magnetic characteristics for a given maturity level. The magnetic study therefore provided constraints on the maturity level of the shales in the oil window area. Our study showed that PM can be used as a geothermometer in shales in which CaCO₃ is lower than 60%.

**Keywords:** shales; magnetite; goethite; oil window; geothermometer

**Citation:** Aubourg, C.; Kars, M.; Pozzi, J.; Mazurek, M.; Grauby, O. A Magnetic Geothermometer in Moderately Buried Shales. *Minerals* **2021**, *11*, 957. https://doi.org/10.3390/min11090957

Academic Editor: Alvar Soesoo

Received: 13 August 2021

Accepted: 30 August 2021

Published: 31 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Shales, which account for two-thirds of the volume of sedimentary rocks [1], are of significant economic interest for geological storage, unconventional resources, and fluid transfers. The diagenesis of shales during burial modifies their petrophysical properties [2], and it is therefore relevant to evaluate the maximum burial conditions in shales. Many geothermometers are used to determine the peak temperature of burials. For burial temperatures below 150 °C, common geothermometers in shales include RockEval pyrolysis [3], vitrinite reflectance [4], fluid inclusions [5], illite crystallinity [6], and apatite fission tracks [7]. In the Paris basin (France), Blaise et al [8] proposed a comparative study of several techniques, including an original study based on the magnetic properties of shales, initially proposed by Aubourg and Pozzi [9]. In this study, we propose to compare the use of the geothermometer according to two approaches; one consisting of describing, in a more complete way than in the study of Blaise et al. [8], the evolution of the magnetic parameters through different sedimentary formations (Drilling EST433); and the other consisting of studying the same formation (Opalinus Clay) at different burial depths in the Swiss Molasse Basin.

<footnote n="1">¹ Laboratoire des Fluides Complexes et leurs Réservoirs, UMR 5150 CNRS TOTAL, Université de Pau et des Pays de l'Adour, Avenue de l'Université, CEDEX, 64013 Pau, France</footnote>

<footnote n="2">² Center for Advanced Marine Core Research, Kochi University, B200 Monobe, Nankoku 783-8502, Japan; mkars@kochi-u.ac.jp</footnote>

<footnote n="3">³ Laboratoire de Géologie, CNRS-Ecole Normale Supérieure PSL University, 75013 Paris, France; pozzi@biotite.ens.fr</footnote>

<footnote n="4">⁴ Rock-Water Interaction, Institute of Geological Sciences, University of Bern, Baltzerstrasse 3, 3012 Bern, Switzerland; martin.mazurek@geo.unibe.ch</footnote>

<footnote n="5">5 Cinam, UMR 7325 CNRS, Campus de Luminy, Aix-Marseille Université, CEDEX 9, 13288 Marseille, France; olivier.grauby@univ-amu.fr</footnote>

<footnote>* Correspondence: charles.aubourg@univ-pau.fr</footnote>
## 2. Rock Magnetism of Shales

Aubourg and Pozzi [9] initially suggested that the magnetic minerals of argillaceous rocks may be representative of modest burial; i.e., less than 10 km (assuming a thermal gradient of 25 °C/km). By analyzing the magnetic properties of shales at low analytical temperature, in the range of 10 K to 300 K (~263 °C to 27 °C), they reported a diagnostic signature, from which a PM parameter could be extracted (Figure 1A). When a strong magnetic field (typically on the order of 2.5 Tesla) is applied to a rock, an isothermal remanent magnetization will be recorded in the ferromagnetic minerals. The remanent magnetization is named the saturated isothermal remanent magnetization (SIRM). If this magnetic field is applied at room temperature (RT = 300 K), a certain fraction of magnetic minerals will record this magnetic field. This magnetization is called an RT-SIRM. This fraction includes single-domain (SD) and multidomain (MD) magnetic grains. For example, for magnetite, these are grains with a size greater than 20 nm [10]. If the magnetic field is applied to the rock at a very low temperature, typically LT = 10 K, another category of grains will also be able to record this magnetic field. In the case of magnetite, grains <20 nm, called superparamagnetic grains (SP), will be able to record a remanent magnetization called LT-SIRM. In summary, if we assume that the ferromagnetic mineral is magnetite, the RT-SIRM is a proxy for SD and MD, whereas the LT-SIRM integrates a larger fraction with SP, SD, and MD. If a rock fragment is heated from 10 K very quickly, the LT-SIRM will lose its intensity if the fraction of SP is important. This decrease is therefore a proxy for the concentration of SP in the rock powder. This is precisely what is measured with the PM parameter: the decrease of LT-SIRM from 10 K to 35 K (Figure 1A). Here, we used the definition of Aubourg and Pozzi [9], where PM = (LT-SIRM<sub>10K</sub> - LT-SIRM<sub>35K</sub>)/LT-SIRM<sub>10K</sub>

Aubourg and Pozzi [9] studied two categories of shales, Opalinus Clay (Dogger) from the Jura mountains and Bure Clay (Callovo-Oxfordian) from the Paris Basin. The peak burial temperature of Opalinus and Bure shales has been estimated at ~85 °C and ~40 °C, respectively. They showed a PM close to 0.97 in Opalinus Clay and 0.64 in Bure Clay (Figure 1A). This indicated that: (1) the rock could contain a very large fraction of nanometric magnetic grains, and (2) this parameter could vary in significant proportions—almost 97% of LT-SIRM was reduced from 10 K to 35 K in the case of the Opalinus Clay. Aubourg and Pozzi [9] performed burial-like heating experiments using a gold capsule in an autoclave at 150 °C (53 days), 200 °C (31 days) and 250 °C (27 days). Hydrostatic pressure was maintained at 100 MPa. Based on rock magnetic analysis of natural claystones and heated claystones, Aubourg and Pozzi [9] suggested that PM displayed an evolution consistent with the peak burial temperature (T<sub>B</sub>) in the range of 50 °C to 250 °C (Figure 1B). The PM evolution followed two branches: a PM<sub>up</sub> branch (50–90 °C) and a PM<sub>down</sub> branch (90–250 °C). The maximum value of PM (PM<sub>max</sub> = 0.97) was observed for a burial temperature near ~90 °C.

This high proportion of nanograins in shales has been observed in different basins. Kars et al. [11–13] observed the diagnostic signal in samples (cores and cuttings) taken from boreholes. Abdelmalak et al. [14] reported the diagnostic signal in argillaceous rocks from volcanic margin basins (East Greenland). Chou et al. [15] reported the diagnostic signal in siltstones from a fold-and-thrust belt in Taiwan. One of the reasons for this ubiquity is that inherited magnetic iron oxides are altered during the early stage of burial in anoxic conditions [16]. This means that the essential, if not all, ferromagnetic minerals are neoformed in argillaceous rocks. A decisive experiment showed that the association of pyrite, organic matter, and a heating at 90 °C allowed the production of magnetite [17]. These components are classically present in shales.

More generally, in the context of shale burial, Aubourg et al. [18,19] proposed the existence of a window in which magnetite would be the main ferromagnetic mineral (Figure 2). The magnetite window covers a burial temperature range from ~60 °C to 330 °C [19]. At ~60 °C, magnetite can coexist with greigite. From ~340 °C, magnetite, then at a very low concentration (a few ppmv), coexists with monoclinic pyrrhotite SD and MD.
Aubourg et al. [19] showed that the concentration of magnetite decreases by one order of magnitude from the oil window entrance (T~B~60 °C) to the metamorphic domain (T~B~330 °C). Above this temperature of 330 °C, pyrrhotite becomes the predominant mineral in the shales (Figure 2). Pyrrhotite is formed from magnetite and pyrite [20].

<figure><img src="image_4.png" /></figure>

**Figure 1.** (A) Representative warming monitoring of artificial isothermal remanence imparted at 10 K at 2.5 Tesla (LT-SIRM) of Opalinus Clay from Mont Terri (Swiss). Note the log scale and the huge dropdown of LT-SIRM, quantified by the PM parameter. This indicated a large predominance of nanoparticles in the superparamagnetic state. (B) Published PM parameter versus peak burial temperature of natural shales (Opalinus and Bure) and laboratory-heated shales. The oil window is indicated. The maximum value of PM is localized in the oil window. Two branches, PM-up and PM-down, were drawn. This indicated that a PM value corresponded to two possible burial temperatures (redrawn from Aubourg and Pozzi [9]).

<figure><img src="image_5.png" /></figure>

**Figure 2.** Artificial remanence (RT-SIRM) of shales from the Pyrenees and Taiwan mountain belt. Data from Aubourg et al. [19]. The black circles correspond to shales in which magnetite was the main ferromagnetic mineral. For the white circles, pyrrhotite was the main mineral. The peak burial temperature (T~B~) of shales was derived from vitrinite analysis and from Raman spectroscopy of carboniferous material at 50 °C to ~330 °C. Pyrrhotite >1 µm was formed at the expense of magnetite and pyrite. Redrawn from Aubourg et al. [19].
### 3. Geological Setting and Sampling

In this study, we combined analysis of Borehole EST433 from the Basin of Paris, [8,21] and Opalinus Clay samples from the Swiss Molasse Basin [22].

The thermal subsidence of the sedimentary basin of Paris began during the late Triassic until the early Cretaceous time [23] (Figure 3). During the Jurassic, shales and carbonates were deposited in a successive sequence of transgression and regression. Borehole EST433 was drilled in 2008, and ~1.2 km of cores of Jurassic rocks have been recovered [24]. Blaise et al. [8] examined the evolution of vitrinite reflectance within rocks from Borehole EST433. The vitrinite reflectance was evolving from 0.35% (~500 m) to 0.75% (~1900 m) (Figure 4). This indicated that the organic matter was immature to early mature in the petroleum sense. This was confirmed by a low Tmax (413–436 °C) in the RockEval analysis. The total organic matter content was generally <1%, except in the Bathonian and lower Toarcian, where content could exceed 10%. The Toarcian claystones had marine organic matter (type II kerogen) and were the only potential oil-prone source rock of the section. Otherwise, the organic matter was mainly terrestrial (type III kerogen) [21]. Peak burial temperatures derived from a set of different techniques ranged between ~40 and ~80 °C from the lower Cretaceous to the Triassic [8]. We collected 33 core samples (Figure 4, Table 1) located essentially in the Callovo-Oxfordian shales and the Lower Jurassic shales. A few months after the sampling party, we collected fragments from the cores. The rock fragments were reduced manually into powder (~500 mg) and sealed in a gel cap. Kars [25] showed that aging rock powder for few months may drastically reduce the PM parameter. For this reason, we measured the different rock magnetic parameters within a few hours after crushing.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** The Paris Basin, showing the locations of Boreholes EST433 and EST211. The white rectangle frames the location of boreholes where Opalinus Clay was sampled.</figcaption></figure>

Opalinus Clay has long been the subject of intense investigation in the Mont Terri project, where the feasibility of geological disposal of nuclear waste is being studied [26]. Further to the east, this formation is also present at various depths in the Swiss Molasse Basin (Figures 3 and 5). Opalinus Clay is a laterally continuous, lithologically homogeneous formation with a regionally comparable diagenetic evolution (Figure 5A) [22,27]. Mazurek et al. [22] used data from apatite fission track analysis, vitrinite reflectance, and biomarker isomerization to constrain peak burial temperatures from boreholes. Peak burial
temperatures were estimated from modeling and available vitrinite data (RO%). When RO% was provided, peak burial temperature was estimated from Easy%R0 thermal modeling by Mazurek et al. [22]. In addition, we proposed the empirical calibration of vitrinite data provided by Barker and Goldstein [28] (Table 2). The highest burial temperature (~115 °C) was reached eastward, for the Herdern, Berlingen, and Kreuzlingen boreholes. Lower peak burial temperature (~85 °C) was estimated for the Opalinus Clay in the Weiach and Benken boreholes. We collected fresh fragments from seven boreholes (Table 2). Rock magnetism was performed right after crushing the fresh fragments to avoid an aging effect.

<table><caption>**Table 1.** Data from Borehole EST433, including depth location of the samples in the borehole. RT-SIRM: room-temperature saturated isothermal remanent magnetization; LT-SIRM: low-temperature saturated remanent magnetization (10 K). The PM value was inferred from the LT-SIRM drop between 10 K and 35 K. PM = (LT-SIRM<sub>10K</sub> - LT-SIRM<sub>35K</sub>)/LT-SIRM<sub>10K</sub>. K: low-field magnetic susceptibility; CaCO<sub>3</sub>: carbonate content.</caption><thead><tr><th></th><th>Depth</th><th>RT-SIRM</th><th>LT-SIRM</th><th>PM</th><th>K</th><th>CaCO3</th></tr><tr><th></th><th>m</th><th>μAm2/kg</th><th>μAm2/kg</th><th></th><th>μSI</th><th>%</th></tr></thead><tbody><tr><td rowspan="15">Callovo-Oxfordian</td><td>531.95</td><td>223</td><td>530</td><td>0.66</td><td>24</td><td>83</td></tr><tr><td>531.95</td><td>26</td><td>141</td><td>0.66</td><td>70</td><td>51</td></tr><tr><td>542.10</td><td>n.m.</td><td>477</td><td>0.51</td><td>8</td><td>35</td></tr><tr><td>547.95</td><td>127</td><td>733</td><td>0.43</td><td>80</td><td>53</td></tr><tr><td>556.10</td><td>14</td><td>126</td><td>0.71</td><td>31</td><td>19</td></tr><tr><td>566.00</td><td>73</td><td>1040</td><td>0.54</td><td>31</td><td>30</td></tr><tr><td>573.11</td><td>131</td><td>893</td><td>0.47</td><td>84</td><td>25</td></tr><tr><td>587.81</td><td>64</td><td>1050</td><td>0.59</td><td>115</td><td>26</td></tr><tr><td>598.00</td><td>84</td><td>1130</td><td>0.50</td><td>138</td><td>30</td></tr><tr><td>609.25</td><td>78</td><td>1130</td><td>0.50</td><td>198</td><td>30</td></tr><tr><td>625.04</td><td>82</td><td>1660</td><td>0.56</td><td>196</td><td>46</td></tr><tr><td>636.85</td><td>48</td><td>1100</td><td>0.55</td><td>132</td><td>24</td></tr><tr><td>654.29</td><td>133</td><td>945</td><td>0.55</td><td>93</td><td>43</td></tr><tr><td>676.29</td><td>50</td><td>706</td><td>0.64</td><td>93</td><td>96</td></tr><tr><td>725.06</td><td>77</td><td>158</td><td>0.27</td><td>12</td><td>100</td></tr><tr><td>756.10</td><td>29</td><td>78</td><td>0.38</td><td>43</td><td>100</td></tr><tr><td>796.50</td><td>7</td><td>n.m</td><td>0.71</td><td>78</td><td>49</td></tr><tr><td rowspan="2">Dogger</td><td>860.14</td><td>25</td><td>570</td><td>0.81</td><td>63</td><td>94</td></tr><tr><td>891.97</td><td>15</td><td>71</td><td>0.51</td><td>19</td><td>5</td></tr><tr><td rowspan="5">Toarcian</td><td>966.38</td><td>58</td><td>2100</td><td>0.86</td><td>364</td><td>5</td></tr><tr><td>968.98</td><td>57</td><td>2680</td><td>0.92</td><td>331</td><td>43</td></tr><tr><td>1032.44</td><td>32</td><td>627</td><td>0.89</td><td>30</td><td>16</td></tr><tr><td>1035.55</td><td>26</td><td>538</td><td>0.89</td><td>57</td><td>49</td></tr><tr><td>1034.45</td><td>15</td><td>58</td><td>0.94</td><td>451</td><td>3</td></tr><tr><td rowspan="4">Domerian Carixian</td><td>1128.47</td><td>92</td><td>5800</td><td>0.93</td><td>257</td><td>3</td></tr><tr><td>1129.63</td><td>n.m.</td><td>7750</td><td>0.94</td><td>266</td><td>13</td></tr><tr><td>1186.05</td><td>57</td><td>3260</td><td>0.97</td><td>296</td><td>46</td></tr><tr><td>1188.40</td><td>57</td><td>3620</td><td>0.97</td><td>225</td><td>10</td></tr><tr><td rowspan="4">Rhetian</td><td>1247.96</td><td>39</td><td>250</td><td>0.61</td><td>75</td><td>11</td></tr><tr><td>1261.14</td><td>n.m.</td><td>68</td><td>0.59</td><td>189</td><td>1</td></tr><tr><td>1263.02</td><td>n.m.</td><td>5920</td><td>0.94</td><td>477</td><td>1</td></tr><tr><td>1270.95</td><td>322</td><td>314</td><td>0.87</td><td>226</td><td>32</td></tr><tr><td>Trias</td><td>1889.89</td><td>34</td><td>288</td><td>0.60</td><td>247</td><td>3</td></tr></tbody></table>

n.m = not measured.
<table><caption>**Table 2.** Opalinus Clay data, including depth location of the samples in the borehole. RT-SIRM: room-temperature saturated isothermal remanent magnetization; LT-SIRM: low-temperature saturated remanent magnetization (10 K); K: low-field magnetic susceptibility; G%: goethite proxy derived by the slope of RT-SIRM during cooling. The PM value was inferred from the LT-SIRM drop between 10 K and 35 K. PM = (LT-SIRM<sub>10K</sub> - LT-SIRM<sub>35K</sub>)/LT-SIRM<sub>10K</sub>. RO%: percentage of vitrinite reflectance. Peak burial temperature derived from an equation by Barker and Goldstein [28]. EasyRO% and modeling are data from Mazurek et al. [22].</caption><thead><tr><th rowspan="2">Borehole</th><th>Depth</th><th>RT-SIRM</th><th>LT-SIRM</th><th>K</th><th>G%</th><th>PM</th><th>RO%</th><th colspan="3">Peak Burial Temperature</th></tr><tr><th>m</th><th>μAm2/kg</th><th>μAm2/kg</th><th>μSI</th><th></th><th></th><th>%</th><th>B&G 89</th><th>EasyRO%</th><th>Modeling</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>°C</th><th>°C</th><th>°C</th></tr></thead><tbody><tr><td>Riniken</td><td>385.16</td><td>52</td><td>1250</td><td>267</td><td>23</td><td>0.83</td><td>0.40-0.44</td><td>42-54</td><td>60-71</td><td></td></tr><tr><td>Schafisheim</td><td>1015.67</td><td>61</td><td>2010</td><td>234</td><td>25</td><td>0.88</td><td>~0.48</td><td>~65</td><td></td><td></td></tr><tr><td>Weiach</td><td>581.86</td><td>72</td><td>1660</td><td>168</td><td>29</td><td>0.78</td><td>0.53-0.57</td><td>77-86</td><td>87-92</td><td>84</td></tr><tr><td>Benken</td><td>635.67</td><td>35</td><td>1330</td><td>299</td><td>10</td><td>0.92</td><td>0.52-0.58</td><td>75-88</td><td>85-93</td><td>85</td></tr><tr><td>Herdern</td><td>1732.38</td><td>80</td><td>2300</td><td>227</td><td>28</td><td>0.86</td><td>n.d.</td><td></td><td></td><td>115</td></tr><tr><td>Berlingen</td><td>1912.78</td><td>61</td><td>2610</td><td>256</td><td>16</td><td>0.91</td><td>n.d.</td><td></td><td></td><td>110</td></tr><tr><td>Kreuzlingen</td><td>2175.58</td><td>59</td><td>1280</td><td>299</td><td>16</td><td>0.85</td><td>n.d.</td><td></td><td></td><td>115</td></tr></tbody></table>

n.d. = not determined.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Peak burial temperature proxies from Borehole EST433. The reflectance of vitrinite (RO%) from ~0.3 to ~0.8 indicated that organic matter maturity was within the oil window. The numbers on the left refer to sampling for magnetic analysis. Data from Blaise et al. [8].</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Location of boreholes from which Opalinus Clay was sampled. (A) Map view; coordinates are on the Swiss km grid. (B) Cross-section; the Herdern borehole location is indicated. Map and geological cross-section adapted from [27].</figcaption></figure>

**4. Methods**

We measured the magnetic susceptibility (MFK1 Kappabridge, Agico, Inc. Brno, Czech Republic) and monitored the remanence at low temperature (<300 K) using a Magnetic Properties Squid Magnetometer (MPMS-XL, Quantum Design, California, CA, USA). Low-temperature experiments allowed us to limit considerably the magnetic transformations, with a detection threshold on the order of one part per billion [29]. This technique is particularly suitable for detecting diagnostic magnetic phase transitions in magnetite [30], pyrrhotite [31,32], and hematite [33]. Referring to the low-temperature measurements, Kars [34] demonstrated that the internal reproducibility (same powder measured several times), as well as external reproducibility (different powders from the same rock fragment), were better than 1%. In this study, some stratigraphic levels were measured several times (different fragments taken at about few centimeters from the core) to check the consistency of the results.

For the low-temperature measurements, we monitored a remanence on cooling or warming obtained either at room temperature (300 K, RT) or at low temperature (10 K, LT). This remanence was acquired with the application of a 2.5 T magnetic field. We expected the saturation achieved and the resulting remanence was then a saturated isothermal remanent magnetization (SIRM). The LT-SIRM was monitored through the warming (10 K to 300K) in a null magnetic field (trapped field < 100 nT). The RT-SIRM was monitored on cooling under a 50 µT magnetic field. This procedure, initially proposed by Aubourg and Pozzi [9], aimed to exacerbate some magnetic transitions of antiferromagnetic minerals [35], and to illuminate an input of induced magnetization that is believed to be carried by nanosized pyrrhotite and magnetite. This is the P-behavior as defined by Kars et al. [35]. For a few specimens, we also performed the monitoring of RT-SIRM on cooling and warming. It should be noted that a few samples displayed a noisy RT-SIRM monitoring curve on cooling due to low magnetization.

We carried out scanning electronic microscopy (SEM) on selected specimens to obtain information on iron sulfides. SEM observation was performed on representative samples from Borehole EST433 (ANDRA, Bure, France). To obtain information on neoformed magnetic minerals, we performed laboratory heating at 95 °C for tens of days according
to the protocol proposed by Aubourg et al. [36]. We then performed a focused ion beam
(with gallium) section of a pyrite framboid. This 0.2 × 5 µm section was analyzed with a
MET JEOL 2000 FX 200 Kv. Chemical analysis of sulfide and oxygen was performed using
a mono-crystal analyzer (Si/Li EDX Bruker).

**5. Results**

5.1. *Bure Samples*

5.1.1. Rock Magnetism

The calcite content of the samples was obtained from ANDRA (Agence nationale pour la gestion des déchets radioactifs). It varies in large proportion, from 100% to near 1%, and can reach almost 100% in marls. By contrast, in shales, the calcite content can be as low as 1%. The low-field magnetic susceptibility (K) of our samples ranged from 8 to 477 µSI. Hysteresis loops (not shown) displayed straight lines, which indicated a dominant paramagnetic behavior. Both low magnetic susceptibility and relatively strong paramagnetic susceptibility suggested that a minute amount of ferromagnetic minerals was present in our samples. Assuming that magnetite was the main magnetic mineral, with the knowledge that the saturation remanent magnetization of magnetite is SIRM < 20 Am²/kg [37], a value for RT-SIRM on the order of 10⁻⁴ Am²/kg (Table 1) gives a magnetite volumetric concentration of >2 ppmv [19]. We checked for possible correlations between the calcite content, the TOC content, and the magnetic parameters. We did not observe significant correlations. However, the magnetic susceptibility tended to be lower when calcite content was elevated. The LT-SIRM was on the order of 10⁻³ Am²/kg; i.e. one order of magnitude more than for the RT-SIRM. The ratio of LT-SIRM/RT-SIRM evolved from 1 (Rhetian sample) to 64 (Carixian sample).

Representative remanence curves of two end members are displayed in Figure 6. The depth between the two chosen samples was ~600 m. In both samples, the calcite content was lower than 30%, and the TOC was about the same (~0.6%). When cycling the RT-SIRM of the Early Oxfordian sample upon cooling and warming (Figure 6A), we noted two main behaviors: (1) near 120 K, a nonreversible Verwey transition pointed to the existence of stoichiometric magnetite (Figure 6A); and (2) upon cooling, the RT-SIRM increased regularly. When larger than ~10%, this increase can be attributed to the goethite contribution [9]. We derived from this behavior a slope that we named G%: G% = 100 × (RT-SIRM<sub>1-50K</sub>RT - SIRM<sub>300K</sub>)/RT-SIRM<sub>150K</sub>. For the Early Oxfordian sample, G% was 11. When warming the LT-SIRM (Figure 6B), there was a steady decrease of remanence from 10 K to 300 K. At 300 K, the residual LT-SIRM was similar to the RT-SIRM. For the Early Oxfordian sample, 40% of the LT-SIRM was lost at 35 K (PM = 0.4). This example was similar to the Callovo-Oxfordian claystones presented by Aubourg and Pozzi [9] from a different bore-hole (EST211) (Figure 1B).

The Carixian sample behaved differently. When cycling the RT-SIRM, the Verwey transition was much more pronounced (Figure 6C), and the parameter G% was less than 10. However, the most remarkable difference was the development of P-behavior. This P-behavior resulted from a superimposition of an induced magnetization, which paralleled the applied field of 5 µT, and the RT-SIRM [35]. Note that the Early Oxfordian sample displayed embryonic P-behavior (Figure 6A). For Carixian claystones, the LT-SIRM on warming displayed a diagnostic feature, with a PM = 0.9 and a marked Verwey transition (Figure 6D). This Carixian sample's rock magnetic signature can be compared to the Early Dogger Opalinus claystones from Aubourg and Pozzi [9]. For selected samples, the cooling of RT-SIRM and the warming of LT-SIRM are shown together (Figure 6E,F). The RT-SIRM pattern displayed P-behavior and a different goethite contribution (Figure 6E). There was a trend toward high goethite contribution and low P-behavior at low depths, and very low goethite contribution and high P-behavior at higher depths. The LT-SIRM pattern displayed a large amplitude of loss between 10 and 35 K (Figure 6F). The shallowest drop corresponded to the shallowest depth.
The evolution of G% and PM are displayed in Figure 7. From 531 m to 637 m depth, G% increased from 0 to a maximum value of 60 (Figure 7A). At larger depths, G% was lower than a few %. The value of G% was not correlated to the content of calcite. The parameter PM displayed a range of 0.27 to 0.97 (Figure 7B). We distinguished three ranges of calcite content: 0–30%, 30–60%, and >60%. When samples with low calcite content (<30%) were taken, a general pattern was observed, with a quite regular increase from 0.4 at lower depths to a maximum value (0.97) at 1186 m (Figure 7B). After the maximum value of PM at ~1030 m, we noticed a sharp decrease in PM, from ~0.97 to ~0.59 at about 60 m depth. Samples fragments from the same horizon provided quite the same PM. For example, at depths of 1035.55 and 1034.45 m (lower Toarcian), five measurements were taken, and the value of PM was remarkably defined (PM = 0.91 ± 0.02). The samples with intermediate content of calcite (30–60%) displayed about the same pattern; however, calcite-rich (>60%) samples were off-trend.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Cooling (A,C,E) and warming (B,D,F) monitoring of artificial remanence acquired at room temperature (RT-SIRM) and 10 K (LT-SIRM), respectively. All data are from Borehole EST433. Some data were published in Blaise et al. [19]. Early Oxfordian shales (A,B) displayed characteristics of immature rocks in the sense of organic matter maturity. The Carixian sample (C,D) displayed characteristics of early mature rocks. Selected samples taken at different depths (E,F); note the gradual increase of PM with depth (F). M is the Verwey transition at ~120 K, and was indicative of stoichiometric magnetite. The P-behavior and PM were indicative of fine magnetic particles in the superparamagnetic state.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** (A) Proxy of goethite (G%) versus depth for Borehole EST433. For depths >700 m, G% was negligible. (B) PM parameter versus depth; calcite content is indicated.</figcaption></figure>

5.1.2. Microscopic Observation

We inspected fragments from Borehole EST433 using SEM, and focused our analysis on pyrite grains. Indeed, Brothers et al. [17] showed that magnetite formed at the expense of pyrite at 90 °C. It is then likely to track magnetite near pyrite framboids. The most common observation all along the borehole were framboid aggregates of euhedral pyrites of different sizes (Figure 8A). However, in Toarcian claystones (968.38 m depth), we observed some rounded grains within framboid aggregates (Figure 8B). SEM inspection of a rounded aggregate of framboids of Oxfordian claystones revealed the presence of oxidized framboids (Figure 8C), compatible with magnetite.

To obtain more information on oxidized pyrite, we ran laboratory experiments on Callovo-Oxfordian claystones (Borehole EST211, depth 675 m). This procedure was designed to increase our chances to collect material for TEM observation. We performed laboratory heating for 67 days following the experimental protocol of Aubourg et al. [36] (atmospheric pressure, half-confined environment).

From the heated sample, we obtained a 5 × 0.2 µm FIB section across a framboid (Figure 9). The framboid slice showed the pyrites contoured by a gel shape and pores. To better understand the nature of this gel, we conducted a transmission electron microscopy (TEM) study coupled with EDX analysis (Figure 10). A section across two pyrite grains displayed S/Fe and O/Fe ratios (Figure 10A). The gel had an O/Fe ratio close to ~1, consistent with magnetite (Fe₂O₃). The pyrites had a S/Fe ratio close to 2, consistent with the formulation FeS₂. There was no evidence of oxidation in pyrite. When mapping sulfide (S) and oxygen (O) (Figure 10B), it was apparent that unaltered pyrite was framed by iron oxides, likely magnetite. A close inspection of the boundary of pyrite showed euhedral iron oxides of less than 20 nm in size embedded in a gel-like material (Figure 11). The nature of the gel itself remains enigmatic, and requires further investigation. Assuming that these grains encapsulated in the gel were magnetite, this size implied that they were in a superparamagnetic state (<20 nm). It is therefore plausible that these grains were responsible for the characteristic magnetic behavior (parameter PM and P-behavior).
<figure><img src="image_11.png" /><figcaption>B) SEM Toarcian claystones (968.38 m)</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>C) SEM of Oxfordian claystones (676.29 m)</figcaption></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

**Figure 8.** Scanning electron microscopy (SEM) observation of shales from Borehole EST433. Framboids of pyrite (A) and partially oxidized pyrite (B) observed in Toarcian claystones. (C) Spherule of former pyrite framboids oxidized observed in Oxfordian claystones and the EDS analysis.

<figure><img src="image_15.png" /></figure>

**Figure 9.** Focused ion beam (FIB) section of a pyrite framboid in Callovo-Oxfordian claystone from Borehole EST211. This sample was heated for 67 days at 95 °C. Lateral view (right) shows that pyrite (dark grey) is framed by a gel (light grey). Porosity (dark) is present within the framboids.
<figure><img src="image_16.png" /><figcaption>**Figure 10.** Transmitted electronic microscopy (TEM) microphotograph of pyrite from the FIB section (see Figure 9). (A) EDX profile between two pyrite grains. In between, iron oxide was detected. This was presumably magnetite. (B) EDX map of oxygen and sulfide. The pyrite grain was not oxidized, and iron oxide—presumably magnetite—framed the pyrite.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 11.** Transmitted electronic microscopy (TEM) microphotograph of pyrite framed by a gel in which nanometric euhedral crystals can be seen. These minerals were probably superparamagnetic (SP) magnetite.</figcaption></figure>

## 5.2. *Samples from the Opalinus Clay*

Seven samples of Opalinus Clay were analyzed magnetically (Table 2). The magnetic susceptibility ranged from 168 to 299 µSI. Both RT-SIRM and LT-SIRM were homogenous in magnitude (average $6 \times 10^{-5}$ Am²/kg and $1.7 \times 10^{-3}$ Am²/kg, respectively). The ratio of LT-SIRM/RT-SIRM ranged between 22 and 43. These values were comparable to those obtained for the Opalinus Clay from Mont Terri studied by Aubourg and Pozzi [9]. G% ranged from 10 to 29, and PM ranged from 0.78 to 0.92. The RT-SIRM and LT-SIRM monitoring curves are shown in Figure 12. RT-SIRM cooling curves showed a Verwey transition at ~120 K, a potential goethite contribution, and P-behavior (Figure 12A). LT-SIRM warming curves displayed a large drop of intensity in the range of 10–35 K (Figure 12B). Those curves were similar to those identified at Mont Terri [9], and also comparable to Bure EST433 shales (Figure 6F). The PM value contrasted with PM values from Borehole EST433, which evolved from 0.40 to 0.94. We distinguished boreholes to the West (Benken, Weiach, Riniken, Schafisheim) from those to the East (Berlingen, Herdern, Kreuzlingen), which experienced a higher burial temperature [22]. Samples from Weiach and Riniken had the same magnetic signatures, with PM values of ~0.8; they also were located along an E–W trend. Further south of those two boreholes, the sample from Schafisheim developed a large P-behavior and a PM of ~0.9. Further east, the sample from Berlingen also displayed a large PM of ~0.9. South of these two boreholes, samples from the Herdern and Kreuzlingen boreholes had a PM of 0.85.
<figure><img src="image_18.png" /><figcaption>**Figure 12.** Magnetic results for Opalinus Clay. Cooling (A) and warming (B) of artificial remanence imparted at room temperature (RT-SIRM) and 10 K (LT-SIRM). The values of RT-SIRM and LT-SIRM were normalized. See Table 2 for the magnitude of the remanence.</figcaption></figure>

# 6. Discussion

## 6.1. *Rock Magnetism*

Our rock magnetic investigation of shales from the Basin of Paris (Borehole EST433) and the Swiss Molassic Basin (Opalinus Clay) showed the presence of a magnetic assemblage constituted essentially of stoichiometric magnetite and potential goethite. Magnetite is commonly observed in moderately buried shales [18,19]. By contrast, goethite is rarely reported in shales because standard magnetic measurements at room temperature or at high temperature (>300 K) are not appropriate [25]. It should be noted that we did not detect hematite or pyrrhotite, characterized by the Morin (260 K) transition and the Besnus (34 K) transition, respectively. However, very small pyrrhotite in the superparamagnetic state was suggested by the P-behavior transition [9].

Observations of samples from the Callovo-Oxfordian at the micron- and nanoscales demonstrated that some magnetite formed at the expense of pyrite, leaving the core of the pyrite preserved from oxidation (Figures 9 and 10). The observation of pyrite grains rimmed with magnetite was in line with the laboratory experiments of Brothers et al. [17], who showed that a moderate temperature elevation of 90 °C of aqueous solutions of ferric-ligand complexes and pyrite resulted in the formation of magnetite. The replacement of pyrite with magnetite is believed to act during diagenesis. Thanks to the FIB section of a framboid previously heated at 95 °C for 67 days (Figure 10), we were able to specify that the magnetite taken in a gel form could have a size lower than 20 nm (Figure 11). We therefore proposed that in the burial processes, at least part of the magnetite in a superparamagnetic state was formed in contact with pyrites.

Aubourg and Pozzi [9] were the first to identify goethite within the Callovo-Oxfordian shales from the Bure borehole. This finding was confirmed in immature shales from Baffin Bay [14] and in shales from the Grès d'Annot formation (France) [11]. Kars et al.
[25] suggested that goethite in Callovo-Oxfordian shales from Bure is not likely an alteration by-product. Goethite occurs as nanoparticles dispersed in the clayey matrix, and not enclosed in other minerals or in organic matter. In Borehole EST433, we found goethite using the G% parameter (Figure 7A). Assuming G% is a proxy of goethite concentration, the trend of G% suggested an increase in concentration from 532 m to 637 m. This goethite window took place in the less-buried shales of the EST233 borehole. Whether this increase was due to burial processes or was related to early diagenesis remains to be demonstrated. The goethite seems to be less present for depths >650 m. Goethite occurrence also was detected in Opalinus Clay from the seven boreholes of the Swiss Molassic Basin. The G% value ranged from 10 (barely detectable) to 29 (easily detectable) (Table 2). This observation suggested that goethite occurrence was more related to rock composition rather than the burial temperature.

## 6.2. Rock Magnetism and Burial

Aubourg and Pozzi [9] and Blaise et al. [8] suggested that the PM parameter derived from low-temperature warming monitoring of artificial remanence is sensitive to burial temperature. Data from Borehole EST433 indicated that the PM was potentially erroneous when a calcite fraction greater than 60% was present (Figure 6B). Moreau and Ader [38] observed in the Paris Basin that calcite-rich samples ($CaCO_3 > 20\%$) preserved inherited magnetic minerals better than clay-rich facies. Sulfate reduction was more efficient in clayey facies, leading to greater dissolution of inherited magnetic minerals. When removing these samples with $CaCO_3 > 60\%$, the PM values displayed up and down paths, with a maximum at a depth of 1186 m. This agreed with the PM-up and PM-down branches initially proposed by Aubourg and Pozzi [9] (Figure 1B). The PM-up branch implied that magnetite formation was continuous during the diagenesis process in the studied temperature range below ~100 °C. This was suggested experimentally by Kars et al. [13] on shales heated from 50 to 130 °C. The data from Borehole EST433 thus supported the experimental work. It was the production of magnetite, probably in framboidal pyrites, that governed the PM value, as the vast majority of these magnetites were below the blocking volume (<20 nm) (Figure 11).

Comparing the EST433 data (Figure 7B) with the initial model of Aubourg and Pozzi [9] (Figure 1B), we noticed two differences: (1) the PM-up branch was not linear in shape, and had a distinguishable convex shape; and (2) the PM-down branch was much more abrupt. The convex shape of the PM-up evolution deserves some comments. It was interesting to compare this evolution to the evolution of the reflectance of the vitrinite (RO%) according to the Easy%RO model [22,39]. Assuming a constant sedimentation rate, a constant heat flow, and a similar conductivity of rocks, the evolution of RO% should draw a convex shape. The reason for this is that the reflectance of vitrinite is a function of both temperature and time to the exposition of this temperature. That the PM-up branch displayed a convex shape has a possible profound implication regarding the mechanism of magnetic mineral formation. Similar to vitrinite, the production of magnetic nanograins is controlled by both temperature and the time of exposure to temperature. At imparted experimental temperature, Aubourg et al. [36] showed that magnetite production continued to increase over time as long as water was available in the system. These conditions must be met during burial. Kars et al. [13] complemented this observation by performing successive steps of increasing temperature. This experiment, even if very simplified compared to the natural burial mechanism, seemed to indicate that at a given temperature, a production of nanomagnetite took place, stabilized, then reactivated when the temperature increased.

Selected PM with calcite content <60% were plotted with reflectance of vitrinite data provided by Blaise et al. [8] in semilog scale (Figure 13A). Both PM and RO% varied linearly with the log of depth. PM ranged from 0.43 to 0.97, while RO% ranged from 0.35 to 0.62. PM therefore covered a wider range of values than RO% for a comparable burial.
From the linear fitting curves, we could deduce an equivalent value of the reflectance of vitrinite (RO%c), following the equation:

$$R0\%_c = (PM + 3.6681) \times 0.4444 - 1.5081 \quad (1)$$

Then, PM was plotted against burial temperature estimated for Borehole EST433 (Figure 13B) [8]. The following semiempirical equation could then be established:

$$T(^{\circ}\text{C}) = 81 \times PM^2 - 74.7PM + 66 \quad (2)$$

where T is the peak burial temperature in the range of 45–70 °C. This equation works for the entry of the so-called “oil window: of shales in the petroleum sense; i.e., the diagenesis of organic matter. We do not claim that this equation is universal. This empirical relationship was the result of a localized burial temperature model in EST433, based on the joint use of multiproxies.

<figure><img src="image_19.png" /><figcaption>**Figure 13.** (A) Vitrinite reflectance data RO% (with standard deviation at one sigma) and PM magnetic parameter as a function of depth (m). Note the log scale of depth. (B) PM magnetic parameter for EST433 samples with CaCO₃ < 60% as a function of burial temperature, as proposed by Blaise et al. [8]. The PM instrumental uncertainty was below 1%. Best fittings were calculated for curves.</figcaption></figure>

It then became interesting to compare the Bure data to those from the Opalinus Clay. For this, we proposed to compare available RO% and PM values. The Opalinus Clay's PM data ranged from 0.78 to 0.92, giving vitrinite equivalents RO%c between 0.47 and 0.53 from Equation (1). From Table 2; the range of published RO% of Opalinus Clay was 0.40 to 0.58, including Mont Terri. Within the uncertainty range of vitrinite reflectance data (generally larger than 10%) and PM values (within 1% from instrumental error), we noticed that there was a good agreement between the range of RO%c deduced from PM and the measured RO%. This first-order observation showed a good agreement between the magnetic data and the vitrinite data from the Callovo-Oxfordian shales and the Opalinus Clay. The three boreholes for which the Opalinus Clay Formation reached the highest burial temperature, close to ~115°C (Herdern, Berlingen, and Kreuzlingen), had PM between 0.85 and 0.91 (Table 2). This temperature did not fit into the calibration range established with Equation (2), and these points should rather be located on the PM-down branch. However, data from EST433 suggested a rapid drop in PM after reaching a maximum (Figure 7B). Additional data are needed to better constrain this trend, as data from
the Paris Basin appeared to show a much faster decrease in PM than data from boreholes
in the Swiss Molasse Basin.

However, this downward trajectory implied that the magnetic nanograin population
reduced beyond a burial temperature between 70 and 90 °C. This may have been either a
dissolution or an increase in size above 20 nm. Kars et al. [13] showed experimentally that
continuous magnetite production predominated up to an experimental temperature of
~130 °C. It should be noted that these experiments were done without pressure, and that
the temperature was probably overestimated compared to natural conditions. In any case,
after the magnetite production phase, most likely associated with the oxidation of pyrites,
it would seem that this process stopped or even reversed. Indeed, Aubourg et al [19] ob-
served a general decrease in the concentration of magnetite down to the metamorphic
domain (Tburial > 330 °C) (Figure 2). A possible mechanism for the gradual disappearance
of magnetites is their progressive replacement by pyrrhotites in the superparamagnetic
state, as initially suggested by Rochette [20] and discussed by Aubourg et al. [19].

This comparison between the Opalinus Clay and the Callovo-Oxfordian from Bore-
hole EST433 showed that it is more necessary than ever to have several markers to cor-
rectly model the burial temperature. For a burial temperature range below ~100 °C, the
magnetic geothermometer has advantages over vitrinite reflectance and the RockEval
technique. Beyond the nature of the organic matter, which can significantly influence the
value of the proxies, it is undoubtedly the possible inheritance of detrital organic matter
that can pose the most problems [40]. The organic matter of a rock in an anoxic environ-
ment can be inherited from eroded rock that has undergone higher burial temperatures
than autochthonous organic matter. In the case of shales, in an anoxic environment, it was
demonstrated that a large part of the detrital magnetic minerals was dissolved in the first
few meters of the sediment, mainly by sulfate-reduction mechanisms [16] This dissolution
of magnetic minerals during the early stages of diagenesis is of paramount importance for
the use of the magnetic geothermometer. Moreau and Ader [38] observed in the Paris Ba-
sin that the dissolution of inherited magnetic minerals in carbonate facies (CaCo₃ > 20%)
was much lower than in carbonate-poor facies, which are generally richer in organic mat-
ter. Our study confirmed this observation, since the magnetic parameters deviated from
the geothermometer trends for carbonate values higher than 60%. However, it should be
noted that satisfactory parameters were obtained for carbonate levels between 30 and 60%
(Figure 7). One of the main flaws of the magnetic geothermometer is the existence of two
solutions for the burial temperatures for an observed PM value (PM-up and PM-down
curves). This is why we recommend its use for burial temperatures lower than ~100°C.
The PM values can be advantageously compared to the values of RO% and Tmax, the prox-
ies of the vitrinite reflectance and RockEval.

7. Conclusions

This study concerns two series of shales; the shales from Borehole EST433 in the Paris
Basin taken from different formations to trace the evolution of magnetic properties with
depth over 1.3 km, and the Opalinus Clay taken from seven boreholes in the Swiss Mo-
lasse Basin, brought to variable burial temperatures (~85 °C to ~115 °C).

These shales showed very comparable magnetic mineral assemblages, with stoichio-
metric nanometric magnetites, and goethite restricted in some formations. Micrometric
and nanometric observations identified nanometric magnetites around the periphery of
pyrite grains in framboid aggregates. It was the concentration of magnetites smaller than
20 nm (superparamagnetic state) that dictated the evolution of magnetic properties at low
temperatures (10–300 K). The PM parameter was used as a proxy to plot this evolution.

In Borehole EST433, the PM values of the shales, in which the calcite concentration
was less than 60%, showed a very clear evolution with depth, with values of 0.4 to 0.9.
Vitrinite values showed, in the same depth interval, values from 0.3 to 0.6%. The convex
shape of PM as a function of depth reflected a temperature and time coupling effect, in a
manner comparable to that of vitrinite. A calibration curve could be obtained, the first in
this range of burial temperatures (~40–70 °C).

Opalinus Clay showed PMs that were consistent with published RO% values. In the
warmest part of the Swiss Molasse Basin, it was likely that decreasing PM values were
consistent with an evolution suggested by laboratory experiments in an anterior study.

Overall, it was noticeable that shales from two different basins, sharing similar rock
magnetic signatures with PM ~0.9, and with the same vitrinite reflectance RO% in the
range of 0.5–0.6, have been the matter of different evaluations of peak burial temperatures
(~70 °C in the Paris Basin, ~90 °C in the Swiss Molasse Basin).

**Author Contributions:** Conceptualization and methodology, C.A., M.K., and J.P. Microscopy: O.G., Opalinus Samples: M.M., C.A.; writing—review and editing, All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partly funded by ANDRA (France) and grants from the Institute of Rock Magnetism (Minneapolis, MN, USA).

**Acknowledgments:** The authors wish to thank all the people who, by their advice and remarks, made it possible to identify this magnetic property of shales. The stays at the Institute of Rock Magnetism (Minneapolis, MN, USA), with the precious advice of M. Jackson, as well as the experimental work of D. Janots (ENS, Paris, France), have also favored this work. The authors would like to thank the editor O. Zhu and the two anonymous reviewers for improving this article.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**References**

1. Shales and Mudstones, Vols. I and II; Schieber, J., Zimmerle, W., Sethi, P.S., Eds.; Schweizerbart Science Publishers: Stuttgart, Germany, 1998; ISBN 9783510651832.

2. Ma, L.; Slater, T.; Dowey, P.J.; Yue, S.; Rutter, E.H.; Taylor, K.G.; Lee, P.D. Hierarchical integration of porosity in shales. *Sci. Rep.* **2018**, *8*, 1–14, doi:10.1038/s41598-018-30153-x.

3. Hart, B.S.; Steen, A.S. Programmed pyrolysis (Rock-Eval) data and shale paleoenvironmental analyses: A review. *Interpretation* **2015**, *3*, SH41–SH58, doi:10.1190/INT-2014-0168.1.

4. Mukhopadhyay, P.K. Vitrinite Reflectance as Maturity Parameter. *ACS 1994*, *1*, 24, doi:10.1021/bk-1994-0570.ch001.

5. Chi, G.; Diamond, L.W.; Lu, H.; Lai, J.; Chu, H. Common problems and pitfalls in fluid inclusion study: A review and discussion. *Minerals* **2021**, *11*, 1–23, doi:10.3390/min11010007.

6. Kübler, B.; Jaboyedoff, M. Illite crystallinity. *Comptes Rendus l'Académie Des Sci.-Ser. IIA-Earth Planet. Sci.* **2000**, *331*, 75–89, doi:10.1016/S1251-8050(00)01395-1.

7. Donelick, R.A.; O'Sullivan, P.B.; Ketcham, R.A. Apatite fission-track analysis. *Rev. Mineral. Geochemistry* **2005**, *58*, 49–94, doi:10.2138/rmg.2005.58.3.

8. Blaise, T.; Barbarand, J.; Kars, M.; Ploquin, F.; Aubourg, C.; Brigaud, B.; Cathelineau, M.; El Albani, A.; Gautheron, C.; Izart, A.; et al. Reconstruction of low temperature (<100 °C) burial in sedimentary basins: A comparison of geothermometer in the intra-continental Paris Basin. *Mar. Pet. Geol.* **2014**, *53*, 71–87, doi:10.1016/j.marpetgeo.2013.08.019.

9. Aubourg, C.; Pozzi, J.-P. Toward a new <250 °C pyrrhotite-magnetite geothermometer for claystones. *Earth Planet. Sci. Lett.* **2010**, *294*, 47–57, doi:10.1016/j.epsl.2010.02.045.

10. Özdemir, Ö.; Dunlop, D.J. Effect of crystal defects and internal stress on the domain structure and magnetic properties of magnetite. *J. Geophys. Res. Solid Earth* **1997**, *102*, 20211–20224, doi:10.1029/97jb01779.

11. Kars, M.; Aubourg, C.; Labaume, P.; Berquó, T.S.; Cavailhes, T. Burial diagenesis of magnetic minerals: New insights from the Grès d'Annot transect (SE France). *Minerals* **2014**, *4*, 667–689, doi:10.3390/min4030667.

12. Kars, M.; Aubourg, C.; Suárez-Ruiz, I. Neoformed magnetic minerals as an indicator of moderate burial: The key example of middle Paleozoic sedimentary rocks, West Virginia. *Am. Assoc. Pet. Geol. Bull.* **2015**, *99*, 389–401, doi:10.1306/06301413006.

13. Kars, M.; Aubourg, C.; Pozzi, J.P.; Janots, D. Continuous production of nanosized magnetite through low grade burial. *Geochem. Geophys. Geosyst.* **2012**, *13*, 1–12, doi:10.1029/2012GC004104.

14. Abdelmalak, M.M.; Aubourg, C.; Geoffroy, L.; Laggoun-Defarge, F. A new oil-window indicator? The magnetic assemblage of clay stones from the Baffin Bay volcanic margin (Greenland). *Am. Assoc. Pet. Geol. Bull.* **2012**, *96*, 205–215, doi:10.1306/07121111008.

15. Chou, Y.M.; Song, S.R.; Aubourg, C.; Lee, T.Q.; Boullier, A.M.; Song, Y.F.; Yeh, E.C.; Kuo, L.W.; Wang, C.Y. An earthquake slip zone is a magnetic recorder. *Geology* **2012**, *40*, 551–554, doi:10.1130/G32864.1.

16. Roberts, A.P. Magnetic mineral diagenesis. *Earth-Sci. Rev.* **2015**, *151*, 1–47, doi:10.1016/j.earscirev.2015.09.010.
17. Brothers, L.A.; Engel, M.H.; Elmore, R.D. The late diagenetic conversion of pyrite to magnetite by organically complexed ferric iron. *Chem. Geol.* **1996**, *130*, 1–14, doi:10.1016/0009-2541(96)00007-1.

18. Aubourg, C.; Pozzi, J.P.; Kars, M. Burial, claystones remagnetization and some consequences for magnetostratigraphy. *Geol. Soc. Spec. Publ.* **2012**, *371*, 181–188, doi:10.1144/SP371.4.

19. Aubourg, C.; Jackson, M.; Ducoux, M.; Mansour, M. Magnetite-out and pyrrhotite-in temperatures in shales and slates. *Terra Nov.* **2019**, *31*, 534–539, doi:10.1111/ter.12424.

20. Rochette, P. Metamorphic control of the magnetic mineralogy of black shales in the Swiss Alps: Toward the use of "magnetic isogrades." *Earth Planet. Sci. Lett.* **1987**, *84*, 446–456, doi:10.1016/0012-821X(87)90009-4.

21. Blaise, T.; Izart, A.; Michels, R.; Suarez-Ruiz, I.; Cathelineau, M.; Landrein, P. Vertical and lateral changes in organic matter from the Mesozoic, eastern Paris Basin (France): Variability of sources and burial history. *Int. J. Coal Geol.* **2011**, *88*, 163–178, doi:10.1016/j.coal.2011.09.001.

22. Mazurek, M.; Hurford, A.J.; Leu, W. Unravelling the multi-stage burial history of the Swiss Molasse Basin: Integration of apatite fission track, vitrinite reflectance and biomarker isomerisation analysis. *Basin Res.* **2006**, *18*, 27–50, doi:10.1111/j.1365-2117.2006.00286.x.

23. Carpentier, C.; Brigaud, B.; Blaise, T.; Vincent, B.; Durlet, C.; Boulvais, P.; Pagel, M.; Hibsch, C.; Yven, B.; Lach, P.; et al. Impact of basin burial and exhumation on Jurassic carbonates diagenesis on both sides of a thick clay barrier (Paris Basin, NE France). *Mar. Pet. Geol.* **2014**, *53*, 44–70, doi:10.1016/j.marpetgeo.2014.01.011.

24. Landrein, P.; Vigneron, G.; Delay, J.; Lebon, P.; Pagel, M. Lithologie, hydrodynamisme et thermicité dans le système sédimentaire multicouche recoupé par les forages Andra de Montiers-sur-Saulx (Meuse). *Bull. La Soc. Geol. Fr.* **2013**, *184*, 519–543, doi:10.2113/gssgbull.184.6.519.

25. Kars, M.; Lerouge, C.; Grangeon, S.; Aubourg, C.; Tournassat, C.; Madé, B.; Claret, F. Identification of nanocrystalline goethite in reduced clay formations: Application to the Callovian-Oxfordian formation of Bure (France). *Am. Mineral.* **2015**, *100*, 1544–1553, doi:10.2138/am-2015-5096.

26. Bossart, P.; Thury, M. Research in the Mont Terri Rock laboratory: Quo vadis? *Phys. Chem. Earth Parts A B C* **2007**, *32*, 19–31, doi:10.1016/J.PCE.2006.04.031.

27. Aschwanden, L.; Diamond, L.W.; Adams, A. Effects of progressive burial on matrix porosity and permeability of dolostones in the foreland basin of the Alpine Orogen, Switzerland. *Mar. Pet. Geol.* **2019**, *100*, 148–164, doi:10.1016/j.marpetgeo.2018.10.055.

28. Barker, C.E.; Goldstein, R.H. Fluid-inclusion technique for determining maximum temperature in calcite and its comparison to the vitrinite reflectance geothermometer. *Geology* **1990**, *18*, 1003–1006, doi:10.1130/0091-7613(1990)018<1003:FITFDM>2.3.CO;2.

29. Dunlop, D.J.; Üzdemir, Ü. *Rock Magnetism*; AA(University of Toronto), AB(University of Toronto): Toronto, ON, Canada, 2001.

30. Jackson, M.; Moskowitz, B.; Bowles, J. The magnetite Verwey transition. *IRM Q.* **2011**, *20*, 1–11.

31. Rochette, P.; Fillion, G.; Mattéi, J.L.; Dekkers, M.J. Magnetic transition at 30-34 Kelvin in pyrrhotite: Insight into a widespread occurrence of this mineral in rocks. *Earth Planet. Sci. Lett.* **1990**, *98*, 319–328, doi:10.1016/0012-821X(90)90034-U.

32. Wolfers, P.; Fillion, G.; Ouladdiaf, B.; Ballou, R.; Rochette, P. The Pyrrhotite 32 K Magnetic Transition. *Solid State Phenom.* **2011**, *170*, 174–179, doi:10.4028/www.scientific.net/SSP.170.174.

33. Özdemir, Ö.; Dunlop, D.J.; Berquó, T.S. Morin transition in hematite: Size dependence and thermal hysteresis. *Geochem. Geophys. Geosyst.* **2008**, *9*, doi:10.1029/2008GC002110.

34. Kars, M. *Calibration and Application of the MagEval Geothermometer in Sedimentary Rocks*; Université of Pau and Pays de l'Adour: Pau, France, 2012.

35. Kars, M.; Aubourg, C.; Pozzi, J.P. Low temperature magnetic behaviour near 35 K in unmetamorphosed claystones. *Geophys. J. Int.* **2011**, *186*, 1029–1035, doi:10.1111/j.1365-246X.2011.05121.x.

36. Aubourg, C.; Pozzi, J.P.; Janots, D.; Sahraoui, L. Imprinting chemical remanent magnetization in claystones at 95 °C. *Earth Planet. Sci. Lett.* **2008**, *272*, 172–180, doi:10.1016/j.epsl.2008.04.038.

37. Maher, B.A.; Karloukovski, V.V.; Mutch, T.J. High-field remanence properties of synthetic and natural submicrometre haematites and goethites: Significance for environmental contexts. *Earth Planet. Sci. Lett.* **2004**, *226*, 491–505, doi:10.1016/j.epsl.2004.05.042.

38. Moreau, M.G.; Ader, M. Effects of diagenesis on magnetic mineralogy in a Jurassic claystone-limestone succession from the Paris Basin. *J. Geophys. Res. Solid Earth* **2000**, *105*, 2797–2804, doi:10.1029/1999jb900319.

39. Sweeney, J.J.; Burnham, A.K. Evaluation of a Simple Model of Vitrinite Reflectance Based on Chemical Kinetics.1. *Am. Assoc. Pet. Geol. Bull.* **1990**, *74*, 1559–1570, doi:10.1306/0C9B251F-1710-11D7-8645000102C1865D.

40. Nzoussi-Mbassani, P.; Copard, Y.; Disnar, J.R. Vitrinite recycling: Diagnostic criteria and reflectance changes during weathering and reburial. *Int. J. Coal Geol.* **2005**, *61*, 223–239, doi:10.1016/j.coal.2004.08.002.