# Rapid light carbon releases and increased aridity linked to Karoo-Ferrar magmatism during the early Toarcian oceanic anoxic event

Eric Font, Luís Vítor Duarte, Mark J. Dekkers, Celine Remazeilles, Ramon Egli, Jorge E. Spangenberg, Alicia Fantasia, Joana Ribeiro, Elsa Gomes, José Mirão, et al.

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

## To cite this version:

Eric Font, Luís Vítor Duarte, Mark J. Dekkers, Celine Remazeilles, Ramon Egli, et al. Rapid light carbon releases and increased aridity linked to Karoo-Ferrar magmatism during the early Toarcian oceanic anoxic event. Scientific Reports, 2022, 12, 10.1038/s41598-022-08269-y. insu-03691309

HAL Id: insu-03691309

https://insu.hal.science/insu-03691309

Submitted on 9 Jun 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_2.png" /></figure>

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

OPEN

# Rapid light carbon releases and increased aridity linked to Karoo–Ferrar magmatism during the early Toarcian oceanic anoxic event

Eric Font<fnref n="1" />,²✉, Luís Vítor Duarte<fnref n="1" />,³, Mark J. Dekkers⁴, Celine Remazeilles⁵, Ramon Eglí⁶, Jorge E. Spangenberg⁷, Alicia Fantasia⁸, Joana Ribeiro<fnref n="1" />,⁹, Elsa Gomes<fnref n="1" />,¹⁰, José Mirão¹¹ & Thierry Adatte¹²

Large-scale release of isotopically light carbon is responsible for the carbon isotope excursion (CIE) of the Toarcian Oceanic Anoxic Event during the Lower Jurassic. Proposed sources include methane hydrate dissociation, volcanogenic outgassing of carbon dioxide and/or thermogenic methane release from the Karoo-Ferrar magmatic province (southern Africa). Distinct small-scale shifts superimposed on the long-term CIE have been interpreted as rapid methane pulses linked to astronomically forced climate changes. In the Peniche reference section (Portugal), these small-scale shifts correspond to distinct brownish marly layers featuring markedly high mercury (Hg) and magnetic mineral concentration. Total organic carbon and Hg increase are uncorrelated, which suggests input of Hg into the atmosphere, possibly released after the intrusion of the Karoo-Ferrar sills into organic-rich sediments. Enhanced magnetic properties are associated with the presence of martite, washed-in oxidized magnetite, inferred to be due to increased aridity on the continental hinterland. This study provides strong evidence for a direct link between the Karoo-Ferrar magmatism, the carbon-isotope shifts and the resulting environmental changes.

The early Toarcian (early Jurassic) was marked by an episode of global warming (seawater temperature rise of 4–7 °C<fnref n="1" />,²), increase in organic carbon burial, ocean anoxia and mass extinction known as the Toarcian Oceanic Anoxic Event (T-OAE, ~183 Ma)³, also referred to as the Jenkyns Event e.g.⁴,⁵. The T-OAE is characterized by one of the largest negative carbon-isotope excursions (CIE; ~–2 to –8‰) of the last 250 Myr⁶ and represents a global perturbation of the carbon cycle e.g.<fnref n="1" />,⁷⁻¹¹. This negative CIE has been ascribed to the release of thousands of gigatons of isotopically light carbon into the ocean–atmosphere system⁷,⁸,¹²⁻¹⁵. Proposed sources of isotopically light carbon include Karoo-Ferrar magmatism (southern Africa; Fig. 1) and associated release of mantle-derived carbon dioxide¹⁶⁻²¹, thermogenic methane emissions from the intrusion of Karoo-Ferrar sills into Gondwanan coal deposits¹³,²², as well as methane release from marine clathrates and terrestrial environments⁷,⁸,²³,²⁴. The T-OAE carbon-cycle perturbation corresponds to a period of low carbonate productivity and is thought be associated with enhanced hydrological cycling and runoff, resulting in an increased delivery of nutrients and

<footnote n="1">¹Department of Earth Sciences, University of Coimbra, 3030 790 Coimbra, Portugal. ²Instituto Dom Luís (IDL), Faculdade de Ciências, Universidade de Lisboa, 1749-026 Lisboa, Lisboa, Portugal. ³MARE and Department of Earth Sciences, University of Coimbra, 3030 790 Coimbra, Portugal. ⁴Department of Earth Sciences, Paleomagnetic Laboratory ‘Fort Hoofddijk’, Utrecht University, Princetonlaan 8a, 3584 CB Utrecht, The Netherlands. ⁵Laboratoire Des Sciences de L’Ingénieur Pour L’Environnement, UMR CNRS 7356 La Rochelle Université, Pôle Sciences et Technologie, Avenue Michel Crépeau, 17042 La Rochelle Cedex 1, France. ⁶Central Institute for Meteorology and Geodynamics (ZAMG), 1190 Wien, Austria. ⁷Institute of Earth Surface Dynamics, University of Lausanne, 1015 Lausanne, Switzerland. ⁸Univ Lyon, UCBL, ENSL, UJM, CNRS, LGL-TPE, 69622 Villeurbanne, France. ⁹Institute of Earth Science - Porto Pole, Rua do campo alegre 687, 4169-007 Porto, Portugal. ¹⁰CITEUC—Centre for Earth and Space Research of the University of Coimbra, 3040-004 Coimbra, Portugal. ¹¹HERCULES Laboratory, School of Sciences and Technology-Geosciences Department, University of Évora, 7000-671 Évora, Portugal. ¹²ISTE, Institut des Sciences de la Terre, Université de Lausanne, GEOPOLIS, 1015 Lausanne, Switzerland. ¹³email: efont@uc.pt</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (A) Global Toarcian palaeogeographic reconstruction⁹⁰ with the location of the Karoo-Ferrar large igneous province in southern Gondwana, the position of palaeoclimatic belts after¹⁰², and the location of the Peniche section (yellow star). (B) Simplified geological map of the Lusitanian Basin, Portugal. (C) Field photographs of the Peniche section. The interval of the Carbon Isotope Excursion (CIE) of the Toarcian Ocean Anoxic Event (T-OAE) is indicated by the white arrow. White dashed lines indicate the position of the three brownish layers. (D,E) Field photographs of the second (labelled MC2 samples) (D) and third (labelled MC3 samples) brownish layers (E).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Peniche section: Mass specific magnetic susceptibility (χ) and Isothermal Remanent Magnetization (IRM) parameters after unmixing using the Kruiver et al.<sup>96</sup> software. SIRM is the saturation of isothermal remanent magnetization. Isotope composition of carbonate (δ<sup>13</sup>C<sub>carb</sub>) and organic carbon (δ<sup>13</sup>C<sub>org</sub>) are from Hesselbo et al.<sup>103</sup> and Fantasia et al.<sup>31</sup>, respectively. δ<sup>13</sup>C<sub>org</sub> shifts identified at Yorkshire (UK, Ref.<sup>5</sup>) are shown for comparison. MC1, MC2, and MC3 correspond to the three brownish layers studied here (see also Fig. 3). D1 and D3 correspond to the discontinuities proposed by Pittet et al.<sup>44</sup>.</figcaption></figure>

terrestrial organic matter to the ocean<sup>15,25-31</sup>. These processes would have caused an increase in the primary productivity and the development of oxygen-deficient conditions in the oceans leading to an increased burial of organic matter<sup>9,30,32,33</sup>.

The long-term CIE (300–500 kyr, Ref.<sup>34–36</sup>) is punctuated by at least three distinct stepwise abrupt δ<sup>13</sup>C negative shifts, with an amplitude of 2–3‰ (vs. VPDB), which are observed in many T-OAE sections, including the Yorkshire section in UK<sup>8</sup> (Fig. 2), the Sancerre borehole in the Paris Basin<sup>14,24</sup>, France, the Belluno Trough section in Italy<sup>37</sup>, the Mochras Farm section in Wales<sup>37,38</sup> and the Sakuraguchi-dani section in Japan<sup>29,39</sup>. These distinct shifts have been interpreted as being due to astronomically controlled millennial-scale carbon injections from methane hydrate or terrestrial permafrost sources<sup>8,12,23,29,40,41</sup>. However, the link with magmatic activity from the Karoo-Ferrar province, as well as the climate response to these rapid carbon injections, are unclear. For instance, some authors argue that the Karoo-Ferrar magmatism would have only contributed to the initial decrease of the CIE, whereas the posterior stepwise negative δ<sup>13</sup>C shifts would have resulted from astronomically forced climate change<sup>8,12,40</sup>. To unravel the nature and origin of these short-scale stepwise negative δ<sup>13</sup>C shifts, we provide new high-resolution (cm-scale) records of the organic carbon-isotope composition (δ<sup>13</sup>C<sub>org</sub>), mercury (Hg) content, whole-rock mineralogy, organic matter geochemistry, magnetic properties, and micro-Raman spectra analysis of three stratigraphic levels associated to the distinct δ<sup>13</sup>C shifts of the lowermost portion of the CIE at Peniche (Lusitanian Basin, Portugal), that includes the Toarcian Global boundary Stratotype Section and Point (GSSP)<sup>42</sup> (Fig. 1).

## Results

Low-field mass specific magnetic susceptibility (χ) of samples from the entire Peniche section (labelled PT samples) varies from 1.14 × 10<sup>-8</sup> to 1.90 × 10<sup>-7</sup> m<sup>3</sup>/kg, typical of marine sedimentary rocks (Fig. 2; Table S1 in the supplementary data). Lowest values are observed at the base of the section, in the limestones and marlstones of the topmost Pliensbachian (uppermost part of the Lemede Formation<sup>43</sup>). The Pliensbachian-Toarcian (Pl/To) boundary interval features a rapid shift in χ, which possibly reflects the presence of a condensed interval (the so-called D1 discontinuity described by Pittet et al.<sup>44</sup>). Field observations indicate that the lithology abruptly changes from limestone to marl-dominated sediments. Moving upsection, χ increases gradually along the *poly-morphum* Zone and reaches maximum values in the lowermost CIE interval of the *levisoni* Zone (Fig. 2). Further
upsection, $\chi$ gradually decreases. Particularly noteworthy in the Peniche section are three atypical prominent $\chi$ peaks that occur at the same level as rapid shifts in $\delta^{13}C_{org}$ (Fig. 2) previously determined by Fantasia et al.³¹. These peaks correspond to discrete ~5–10 cm-thick distinct brownish marly layers (Fig. 1D,E).

Two principal ferromagnetic phases are identified in the early Toarcian of the Peniche sediments based on the unmixing of Isothermal Remanent Magnetization (IRM) acquisition curves (Fig. S1 in the supplementary data). The low-coercivity phase is characterized by a skewed coercivity distribution modelled by two Gaussian functions (components 1 and 2) with median acquisition fields (noted as $B_{1/2}$) of ~15–25 mT and ~45–55 mT and dispersion parameter (DP) of 0.2–0.4 log mT and 0.27–0.34 log mT, respectively. First-order reversal curves (FORC) (Fig. S2) show that this low-coercivity phase is composed by ~97% of >0.1 µm detrital (titano-)magnetite with the typical signature found in sedimentary rocks, and for the remaining ~3% by <0.1 µm magnetite, possibly of secondary origin⁴⁵⁻⁴⁸. The lack of biogenic (secondary) magnetite, typically identified by DP < 0.2 and a narrow central ridge in FORC diagrams, and which is generally abundant in pelagic sediments⁴⁹, is common in hemipelagic environments with high sedimentation rate and low TOC (Peniche: ~1 m/kyr during the early Toarcian cf.⁵⁰). The ratio of components 1 and 2 is relatively constant at ~1:19 throughout the section (Fig. S3), confirming their common origin from a single mineral phase with fixed properties. The high-coercivity phase is represented by a single Gaussian function (component 3) with $B_{1/2}$ of ~200–400 mT and DP of 0.30–0.50 (log mT), typical of hematite. The relative content in magnetite (components 1 and 2) and hematite (component 3) is expressed by the values of the IRM at saturation (SIRM) of the corresponding components in Figs. 2 and 3. SIRM of magnetite and hematite is low in most samples but typically one order of magnitude higher in the brownish levels. In these levels hematite's contribution to the SIRM increases from a base-line value of ~15 to ~37%.

To uncover the environmental significance of these magnetically enhanced brownish levels associated with the rapid $\delta^{13}C$ shifts, we measured their $\chi$, IRM, $\delta^{13}C_{org}$, TOC and bulk Hg content at a cm-scale resolution (MC transects; Fig. 3). Both the base and top of the brownish layers feature low $\chi$ and SIRM values with a gradual increase towards maximum values more to their centers (Fig. 3). Brownish layers show $B_{1/2}$ values of 350–400 mT for the high-coercivity component 3 while surrounding sediments have lower values of ~200–280 mT (Fig. S3 and Table S3 in the supplementary data), suggesting that the former contain a distinct population of hematite. $\delta^{13}C_{org}$ shifts of 2–3‰ previously obtained by³¹ in Peniche and by⁸ in Yorkshire are confirmed here by our new measurements illustrated in Fig. 2. These shifts, named A, B and C in Yorkshire⁸ (Fig. 2), are shown to be abrupt in the case of the first (MC1) and third (MC3) brownish layer at Peniche, occurring within a 5–7 cm-thick interval. In these levels concentration-dependent magnetic properties ($\chi$ and SIRM) change abruptly as well (Fig. 3; Table S4 in the supplementary data). The second brown layer (MC2) shows more gradual shifts in measured parameters within a 14 cm-thick interval. The three brownish layers also feature lower TOC contents of ~0.1–0.2% compared to the over- and underlying sediments which have twice these values (~0.2–0.6%, Table S4). Notably, Hg varies from 20–30 ppb in the over- and underlying sediments up to 50–85 ppb in the brownish layers. This Hg enrichment is coeval with the increase in $\chi$ and SIRM and the decrease in $\delta^{13}C_{org}$ observed in these brown beds (Fig. 3, Table S4).

Optical microscopy coupled to micro-Raman analyses of six MC samples revealed the presence of non-magnetic Ti-oxides (mainly anatase and rutile) in all samples (Fig. S5 in the supplementary data), possibly of a detrital origin or as the product of reductive dissolution of primary iron oxides in a sulfidic anoxic environment. All samples contain abundant framboidal pyrite, mainly associated to organic matter (Fig. S6). Pyrite framboids of heterogeneous grainsize associated with wood fragments suggest that they did not nucleate in the water column but formed after sedimentation, by alteration of organic matter. Well-preserved pyrite and the overall absence of alteration to hematite or goethite in the Raman spectra indicate that the samples did not suffer late diagenetic alteration or oxidation (Fig. S5-S6). The brownish layers contain abundant Ti-bearing iron oxides (Fig. S7), which are exceedingly rare or completely absent in other samples. Crystals exhibit a cubic to octahedral habit typical of magnetite (Fig. S7.A-E) with different degree of alteration expressed by dissolution and exsolution features (Fig. S7). Semi-hexagonal habits typical of hematite or ilmenite are also observed (Fig. S7.G-H). Micro-Raman analyses reveal that some crystals are made of a mixed composition of magnetite and hematite, whereas others are entirely hematite, suggesting that the magnetite host underwent significant martitization (Fig. 4). Relics of magnetite in hematite are rare but were identified under reflected light microscopy (plane and crossed polarized light) in sample MC3-33 (Fig. S7.I). Some hematite crystals have subhedral forms (Fig. S7.J), typical of magnetite, but, occasionally, tabular hematite crystals are observed (Fig. S7.K).

## Discussion

A significant magmatic emplacement event has been linked to major changes in the T-OAE paleoenvironment. Specifically, the timing of the emplacement of the Karoo-Ferrar large igneous province (LIP) appears to coincide with the Pl/To mass extinction, recorded in marine and terrestrial sedimentary settings¹⁶,¹⁷,⁵¹. Paleomagnetism and geochronology (U–Pb and Ar-dating) of intrusive and extrusive magmatic rocks from the Karoo and Ferrar LIPs bracket the main emplacement event to 182–183 Ma, with a duration of 300–500 kyr¹⁹,²⁰,⁵²,⁵³. The age of the Pl/To boundary is placed at 183.6 + 1.7/– 1.1 Ma based on U–Pb geochronology and biostratigraphy¹⁷. Based on the correlation of U–Pb dates of zircon from ash beds deposited during the early Toarcian in Palquilla (southern Peru) with previously published ages of Karoo intrusives, Sell et al.⁵⁴ place the Pl/To boundary at no younger than 183.5 Ma. Starting from an age of 183.5 Ma for the Pl/To boundary, Sell et al.⁵⁴ and Burgess et al.⁵⁵ argue that Ferrar magmatism begins within uncertainty at the onset of the T-OAE and the CIE, and continues through this interval. Based on high-precision U–Pb and ⁴⁰Ar/³⁹Ar geochronology, Greber et al.⁵⁶ suggest that Karoo-LIP magmatism occurred after the Pl/To boundary and therefore postdated the extinction, but support the onset of Karoo-LIP activity and the T-OAE and the early Toarcian warming being synchronous. However, the resolution of radiometric methods does not allow for a more precise correlation between the discrete stepwise
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Three brownish layers (MC1, MC2, MC3). From left to right: Magnetic susceptibility (χ), saturation isothermal remanent magnetization (SIRM) of magnetite and hematite after unmixing using the Maxunmix software⁹⁷, organic carbon isotope composition (δ¹³Corg), total organic carbon (TOC), mercury content and Hg/TOC ratio. These distinct layers are coeval to the three rapid shifts in carbon isotope composition (shifts A–C) illustrated in Fig. 2.</figcaption></figure>

negative δ¹³C shifts within the CIE and their possible linkage to the emplacement and eruptive history of the Karoo-Ferrar LIP. Also, the age estimate of the Pl/To boundary is not well-constrained enough to definitively link it to the Karoo-Ferrar LIP.

The stepwise pattern of the negative CIE and presence of distinct δ¹³C shifts have been tentatively linked to Karoo-Ferrar LIP magmatism following two model scenarios. The first scenario invokes the liberation of huge amounts of isotopically light carbon to the ocean-atmosphere system through methane hydrate dissociation controlled by astronomically forced climate variations⁸,¹²,⁴⁰. In this model scenario, volcanogenic CO₂ from the Karoo-Ferrar magmatic events would not have effused fast enough to cause the δ¹³C shifts. Instead, magmatism would have only contributed to the slower initial ~2‰ decrease in δ¹³C that began before the first methane hydrate dissociation pulse inferred from the Yorkshire record⁸ (Fig. 2). Global warming from volcanic outgassing
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Typical reflected light microscopic photographs of martite and corresponding micro-Raman spectra observed in the brownish layers. Particles show the typical habit of magnetite with a mixed composition of magnetite (M) and hematite (H). I, II and III in the top right panel belong to spots I to III.</figcaption></figure>

might then have acted in concert with astronomical climate forcing to drive episodic dissociation of methane hydrate, resulting in distinct stepped negative δ¹³C shifts. This hypothesis is supported by spectral analysis of carbon-isotope and carbonate abundance that revealed a regular cyclicity, compatible with precession⁸,⁴⁰ or obliquity¹². Assuming that the astronomical forcing is controlled by orbital precession cycles (~21 kyr in the early Jurassic), Kemp et al.⁸ estimated the duration of the three pulses of methane hydrate dissociation to <2 kyr each.

The second scenario implies a direct link between the Karoo-Ferrar eruptions and the stepped δ¹³C shifts, either by repeated pulses of massive degassing of volcanogenic CO₂¹,¹⁵ or by the release of thermogenic CO₂ and CH₄ from organic-rich sediments in southern Gondwana induced by the intrusion of magma¹³,²². These magmatic processes are not predicted to exhibit astronomical periodicity that would explain the cyclic nature of
the stepwise δ¹³C shifts. Both hypotheses can be tested with high-resolution Hg chemostratigraphy. Hg has been used to infer the causal relationship between LIP activity and extinction events recorded in marine and terrestrial sedimentary rocks e.g.,⁵⁷⁻⁶². Hg enrichments recorded during the T-OAE interval at several sites from Europe and the southern hemisphere have been interpreted as evidence for a causal link between Karoo-Ferrar volcanism and coeval environmental changes³¹,⁶². In the Peniche section abrupt and short-lived positive excursions in Hg concentrations and Hg/TOC ratios coincide with the PI/To boundary and with the onset of the negative CIE during the T-OAE⁶². However, a precise correlation of Hg content and the distinct δ¹³C shifts is beyond the sample resolution of those studies (10–50 cm in⁶²; 20 cm in³¹). The resolution in our data set (2–5 cm sample spacing) allows for these two scenarios to be investigated. Findings suggest a direct link between the distinct stepwise pattern that is superposed on the negative CIE and magmatic pulses associated with Karoo-Ferrar activity. We show that the distinct discrete brownish marly layers correspond to the stratigraphic intervals where rapid shifts in δ¹³C org (called Shift A, B and C by analogy with the Yorkshire section, Fig. 2) were documented previously³¹. An absence of significant variations in oxygen (OI) and hydrogen (HI) values (Table S4) suggests that the decrease in δ¹³C org values cannot be explained by a change in the nature and composition of the organic matter (see also⁶³). Nonetheless, each shift is associated with increased Hg content (Fig. 3). For example, shift B shows a rapid and striking increase up to 85 ppb Hg, coeval with the onset of the δ¹³C org decrease. Hg contents then gradually decrease concomitantly with the decrease in δ¹³C org values (Fig. 3) suggesting that these events are linked. Due to the fact that Hg has a strong affinity for organic matter, correlated Hg and TOC values in modern and ancient sediments are generally taken to reflect the input of Hg bound to organic matter or detrital particles (clays) from terrestrial run-off⁵⁷,⁶⁰,⁶⁴. Alternatively, some authors⁶⁴ infer that Hg enrichment within the T-OAE interval may be caused by local redox conditions and/or remobilization of terrestrial Hg. At Peniche, the general increase in Hg background values suggest a link between the T-OAE and the onset of Karoo-Ferrar volcanic activity³¹,⁶². We noted that localized factors, such as the presence of pyritized wood fragments could also have affected the Hg record³¹. However, in the case of the brown layers, Hg does not correlate with TOC (Fig. 3) or clay contents (Fig. S7 in the supplementary data). In the Hg enriched brownish layers, TOC are at least two times lower than that in the over- and underlying sediments (Fig. 3). Because pyrite is mainly associated with organic matter (Fig. S4 in the supplementary data), lower TOC would also point to lower pyrite contents, suggesting that post-depositional Hg remobilization and fixation by pyrite is not the cause of the Hg enrichment. After normalizing for varying TOC levels, Hg/TOC peaks remain a prominent feature within the brownish marls, with values five times higher than background levels (Fig. 3). Note that Grasby et al.⁵⁷ suggest that very low TOC contents close to the limit of detection may lead to spurious results and misinterpretation and proposed that Hg/TOC values should not be calculated for <0.2wt% TOC. For this reason, we only consider the bulk, non-normalized, values of Hg. With or without normalization to TOC, these data suggest that a rapid Hg increase is coeval with a rapid δ¹³C shift, which we infer to be from direct atmospheric deposition of Hg sourced from Karoo-Ferrar magmatic events. Whether this Hg increase results directly from volcanic outgassing or corresponds to liberation of thermogenic Hg after the intrusion of Karoo-Ferrar magma into organic-rich sediments remains to be elucidated.

As volcanogenic carbon has δ¹³C < -6‰, previous workers suggested that the Karoo-Ferrar volcanic eruptions would have not effused fast enough to directly cause the largest T-OAE excursion recorded worldwide⁷,⁸,⁶⁵. Volcanic degassing is a potentially large but poorly constrained source of mercury⁶⁶. Basaltic rocks with minimal crustal contamination have generally low Hg concentrations ranging from 0.2 to 7 ppb⁶⁷⁻⁶⁹. For example, the Karoo dolerite sill (chilled margins) have average Hg concentrations of 2.5 ppb⁶⁷. The low concentration of Hg in mantle-derived magmatic rocks is due to the fact that Hg is a chalcophile element and behaves as an incompatible element with silicate mineral assemblages. The amount of gaseous elemental mercury released during the eruptions of continental flood basalts will thus mainly depend on the Hg contamination from the crust. Most of the Karoo-Ferrar LIP is represented by sills and dykes intruded into Carboniferous to early Jurassic organic-rich sediments²². Based on mass balance calculations constrained by contact aureole around the sills, Svensen et al.²² suggested that 27,400 Gt of CO₂ may have been released from the entire Karoo basin during the intrusive event. The released gas would have a δ¹³C isotopic composition close to that of bulk oil and organic matter (~30 to ~25‰)⁷⁰; at the same time significant amounts of Hg are released from the organic-rich sediments. Decomposition of gas hydrate deposits associated to hydrothermal activity, such as in the eastern shelf of Sakhalin (Russia), may release significant amounts of mercury into shelf waters⁷¹. In this case, Hg is supplied by hydrothermal fluids from deeper fault zones. However, no evidence of hydrothermalism is observed at Peniche. In addition, the main argument against the hypothesis that methane hydrate destabilization was the driver for the three discrete δ¹³C shifts discussed here, relies on the time needed to recharge the methane hydrate reservoir. The duration of the entire Toarcian CIE has been estimated to be comprised between 300 and 900 kyr³⁶,⁷². A more conservative estimate of 300–500 kyr has been calculated for Peniche, Sancerre (Paris), Yorkshire (UK) and the Talghemt section (High Atlas in Morocco)³⁴,³⁵. With a thickness of ~10 m for the negative CIE at Peniche (Fig. 2), a sedimentation rate of ~2 to ~3.3 cm/kyr would result. Considering that each brownish layer is spaced by 1–3 m, the time span between each δ¹³C shift is ~50 to ~150 kyr. These durations are way too short for the methane hydrate reservoir to be recharged as explained next. Based on coupled climate model simulations of orbitally induced changes in ocean circulation and intermediate water temperature during the Paleocene-Eocene Thermal Maximum, Lunt et al.⁷³ calculated that the submarine gas hydrate stores recharge on a timescale of 500 kyr, i.e. 3 to 10 times longer than the interval at Peniche. Using the GEOCLIM carbon cycle model, Heimdal et al.⁷⁴ show that a total of 20,500 Gt C replicates the Toarcian pCO₂ and δ¹³C proxy data, and that thermogenic carbon represents a plausible source for the observed negative CIEs. In this scenario, the extremely isotopically depleted carbon signature of methane clathrates is not required to replicate the negative CIEs. Therefore, we suggest that the Hg enrichment associated with the distinctively stepwise abrupt δ¹³C negative shifts observed at Peniche is better explained by the rapid, sporadic and massive release of thermogenic carbon and Hg from intrusive events of the Karoo-Ferrar LIP.
Interestingly, the three distinctive brownish beds straddling the δ¹³C org shifts are characterized by higher Hg levels and show a strong increase in magnetite and hematite concentration (Figs. 2 and 3). Micro-Raman analyses and optical microscopy reveal that magnetite and hematite correspond to martites (Fig. 4), hematite pseudomorphic after magnetite that typically forms under oxidizing conditions (e.g.,⁷⁵). Provenance changes (composition of the sediment source, aeolian contribution, etc.) and sea-level changes may result in fluctuations of the magnetic signal in marine sediments. This pertains to the relative contribution of diamagnetic calcite from carbonate productivity and dia-/para-/ferromagnetic materials washed-in from terrestrial detrital sources. Also, diagenesis and precipitation of new magnetic minerals at the redox front of the sediment/water interface may yield magnetic anomalies⁷⁶. This goes as well for aeolian input of dusts rich in magnetic particles or climate induced production of magnetic soils, both of which result in a magnetically enhanced detrital input. In the Peniche section, phyllosilicate and calcite are uncorrelated to χ in the brown bed which indicates that the enhanced magnetization is not primarily controlled by lithology in these layers, or by sea-level fluctuations (Fig. S8). In addition, the stratigraphic positions of the brownish layers do not correspond to stratigraphic discontinuities, like those described by Pittet et al.⁴⁴ at the PI/To boundary (D1) and at the *polymorphum-levisoni* zonal transition (D3), that likely reflect transgression-related sediment starvation (Fig. 2). The close association of the magnetically enhanced levels with the distinct δ¹³C shifts while the remainder of the section features background magnetic levels, together with the presence of well-preserved pyrite also argue against a late diagenetic imprint due to burial, fluid circulation or oxidation. Rather, the close relationship between the sediment color, increased Hg contents and δ¹³C org shifts suggests that magnetic enhancement and the presence of martite are the consequence of climate and environmental changes triggered by magmatic activity of the Karoo-Ferrat LIP. Martitization is a widespread phenomenon in various geological settings, however the processes involved in its formation are not fully understood (e.g.,⁷⁵,⁷⁷). Martite occurs in banded iron formations, ore deposits and continental red beds but has been rarely described in marine hemipelagic sediments such as those studied here. The transformation of magnetite into hematite has been considered a redox reaction linked to increased oxygen fugacity. However, the oxygen fugacity in many geologic environments is deemed too low to account for hematite formation via classic redox reactions, leading some authors to consider non-redox formation reactions involving the leaching of Fe²⁺ ions by acidic solutions⁷⁵,⁷⁷,⁷⁸. Although oxygen fugacity is relatively low in hemipelagic and pelagic settings, redox transformation of magnetic minerals may occur in some organic-rich marine sediments. It has, for example, been documented in so called "burnt down" fronts of Mediterranean sapropels⁷⁶,⁷⁹⁻⁸³. Sapropels are deposited under anoxic conditions, and they can be partially or completely oxidized upon re-establishment of oxygenated bottom water ventilation⁸⁴,⁸⁵. During this oxidation, authigenic iron oxides are precipitated, which leads to enhanced magnetizations in the oxidation front. Upon re-establishing the oxygen diffusion into the sediment, it reacts with Fe²⁺ liberated from pyrite and also with upward-diffusing Fe²⁺ from the organic rich sapropel sediments below the oxidation front. However, in contrast to the present study, magnetite is the main magnetic mineral formed in the sapropel oxidation front: hematite is barely or not at all formed⁷⁹⁻⁸². When present, the minute hematite contribution is thought to reflect variations in the input of aeolian (Saharan) dust⁸⁶. The entire Peniche section is characterized by low concentrations of redox sensitive elements (V and Mo) indicating that the depositional conditions were not favorable to the development of long-term anoxic conditions, but was rather dysoxic³¹. Pyrite framboids are formed in situ and precipitated in close association with organic matter (Fig. S6). They do not show any signs of oxidation or dissolution under scanning electron and optical microscopes, or compositional changes in the Raman spectra (Fig. S6). Periods of intermittent bottom-water oxygenation can be created by turbidites, which are rather frequent in the Peniche section during the *levisoni Zone*⁸⁷,⁸⁸. However, no turbidites are associated with the brownish layers; they are documented at least a few decimeters away from them (Fig. 1D,E). In-situ martitization by non-redox reactions involving acid solutions is also unlikely given the buffering capacity of the prevalent carbonates in the sediment. Therefore, we argue that martitization took place in the hinterland, followed by martite deposition in the sedimentary basin. For instance, martite has been found in the C-horizon of loessic soils, where it forms as weathering product of parent minerals⁸⁹. The increased relative concentration of the hematite magnetic component in the brown levels can thus be ascribed to episodes of increased soil erosion. We note here that Peniche is located at the boundary between warm temperate and arid climate during the Toarcian⁹⁰ (Fig. 1A). Oxic and arid conditions on the continent during the *levisoni Zone* were also suggested by Rodrigues et al.⁹¹ based on the characterization of the terrestrial organic matter. The region is thus sensitive to record rapid and brief shifts in paleoclimatic belts, which provides an explanation for the distinctive magnetic signature of the brownish layers. Increased aridity would have resulted in enhanced martitization, as well as in a higher level of mineralization of the organic matter and/or lower productivity, explaining the lower TOC values observed in these samples. The brownish layers also appear to contain the lowest amounts of calcite (Fig. S8), which indicates reduced carbonate production. Lower carbonate production during the T-OAE is hypothesized to have been at least partly triggered by ocean acidification linked to increased atmospheric pCO₂⁹²,⁹³. The boron isotopic composition (δ¹¹B) of brachiopod shells from Peniche shows that the seawater pH decreased ~ 0.6 units just prior to and within the T-OAE, in conjunction with a lower CaCO₃ content in the sediments, concurring with changes in atmospheric pCO₂ derived from stomatal indices⁹⁴,⁹⁵. In this context, martitization on the continent could have been possibly mediated by more acid conditions due to a higher pCO₂ and SO₂ release triggered by Karoo-Ferrat magmatic activity.

## Conclusions

We show that stepwise abrupt δ¹³C negative shifts observed in many T-OAE section worldwide and featuring the onset of the long-term early Toarcian carbon isotope excursion, are represented at Peniche by discrete (~5–10 cm-thick) layers with a marked brownish colour, and distinct geochemical and mineralogical properties. These brownish layers exhibit increased Hg contents and Hg/TOC ratios concomitant with the rapid decrease of
$\delta^{13}$C of ~2–3‰, suggesting a direct link between the Karoo-Ferrar magmatism and the rapid injection of light carbon into the atmosphere. The time span between the deposition of each brownish layer is too short for the discrete shift to be due to destabilization of methane clathrates but can, instead, be explained by sporadic and massive release of thermogenic carbon dioxide and methane, together with Hg, from the intrusion of Karoo-Ferrar sills into organic-rich sedimentary host rocks. In addition, we show that the increased magnetic properties of these brownish layers result from the presence of martite, i.e., hematite pseudomorphs after magnetite that typically form under conditions of increased oxygen fugacity. The presence of these martites is an indicator of palaeoenvironmental and palaeoclimate change, specifically increased aridity, caused by the early Toarcian global warming event.

## Methods

**Study locations.** The Pliensbachian and Toarcian of the Peniche section crops out on the western coast of Portugal and is generally composed of an alternation of hemipelagic marl and limestone deposited on a homoclinal ramp system in the distal part of the Lusitanian Basin, e.g.,⁴³ (Fig. 1). During the early Toarcian, the Peniche area received siliciclastic sediments from the emerged Berlenga-Farilhões Horst located to the west⁸⁷,⁸⁸. The studied interval is about 35 m thick and includes the *emaciatum* ammonite-Zone of the upper Pliensbachian, and the *polymorphum* and *levisoni* Zones of the lower Toarcian⁴³ (Fig. 2). The upper ammonite Zone belongs to the Lemede Formation (Pliensbachian) and consists of decimetric-thick alternations of bioturbated limestone and marlstone. The beds often show shell accumulations and are capped by a condensed interval just above the PL/To boundary (topmost portion of the Lemede Formation). The lower Toarcian (Cabo Carvoeiro Formation) is locally characterized by bioturbated marl-dominated deposits. The lower-middle part of the *levisoni* Zone contains mixed-carbonate siliciclastic turbidites and abundant wood fragments. Detailed mineralogical and geochemical data of this interval of the Peniche section, including clay mineralogy, organic matter geochemistry, major and trace elements, and carbon isotope stratigraphy are published in Fantasia et al.³¹.

**Sample collection and processing.** Magnetic susceptibility ($\chi$) and Isothermal Remanent Magnetization (IRM) data shown in Fig. 2 were measured on 210 rock fragments (labelled PT) collected with a sampling spacing of ~20 cm along the entire ~35 m-thick of the Peniche section. These samples correspond to the same samples previously studied by Fantasia et al.³¹. Samples from the three levels of high magnetic susceptibility (labelled MC1 with 20 samples, MC2 with 33 samples and MC3 with 28 samples) corresponding to the brownish layers associated with the rapid $\delta^{13}$C$_{org}$ shifts A, B and C (Fig. 3) were collected at a later stage with a higher, cm-scale, sampling resolution. The MC samples were crushed in the laboratory using an agate mortar.

Magnetic measurements were performed at the Paleomagnetism Laboratory of the Instituto Dom Luís (IDL) of the University of Lisbon and the Department of Earth Sciences of the University of Coimbra, both in Portugal. Magnetic properties of the MC samples include the measurement of mass specific magnetic susceptibility ($\chi$, in m³/kg) and determination of IRM acquisition curves. Stepwise IRM (typically 30 steps) was acquired with an impulse magnetizer (model IM-30; ASC Scientific) and measured with a JR6 magnetometer (AGICO; measurement range up to 12,500 A/m; noise level and sensitivity of 2.4 × 10⁻⁶ A/m for a standard-sized (10 cc) paleomagnetic sample). The maximum applied field was 1.2 T. IRM acquisition curves were unmixed into several components based on cumulative log-Gaussian (CLG) functions with the Kruiver et al.⁹⁶ software or on skewed generalized (log) Gaussian functions with the Max UnMix program (Maxbauer et al.⁹⁷) to isolate the contributions of magnetite and hematite (denoted here as SIRM of component 1–2 and component 3, respectively). One sample from MC3 was further characterized with high-resolution FORC measurements performed with a Lake Shore 8600 Vibrating Sample Magnetometer at the material magnetism laboratory of the Central Institute of Meteorology and Geodynamics in Austria (13 stacks of 851 curves each, acquired in regular field steps of 0.35 mT) and processed with VARIFORC⁹⁸,⁹⁹.

Carbon isotope compositions of bulk organic matter ($\delta^{13}$C$_{org}$, ‰ VPDB) were determined on decarbonated (10% HCl treatment) powdered whole-rock samples by flash combustion, using an elemental analyser Carlo Erba 1108 elemental analyser connected to a Thermo Fisher Scientific Delta V Plus isotope ratio mass spectrometer at the Institute of Earth Surface Dynamics of the University of Lausanne (IDYST–UNIL), Switzerland. The repeatability and intermediate precision of the $\delta^{13}$C$_{org}$ values was better than 0.05‰.

The C$_{org}$ concentrations were obtained by Rock–Eval pyrolysis following the procedure described by Behar et al.¹⁰, at the Institute of Earth Science of the University of Lausanne (ISTE-UNIL). The determined parameters are total organic carbon (TOC), the Hydrogen Index (HI as mg HC/g TOC) and the Oxygen Index (OI as mg CO$_2$/g TOC), which permit an overall characterization of the sedimentary organic matter.

Mercury (Hg) content of the MC samples was measured with a Zeeman R-915F (Lumex, St-Petersburg, Russia) high-frequency atomic absorption spectrometer set at Mode 1 (700 °C) at the ISTE-UNIL. Analyses were performed on whole-rock powdered samples. The measurements were done in duplicates, and a certified external standard (GSD-11, Chinese alluvium, Hg: 72 ± 6 ppb) was used for calibration purposes ($r = 0.99$, for measured vs certified values) and to guarantee the analytical quality.

Micro-Raman analyses of iron oxides were conducted on six samples (MC1-33, MC2-31, MC3-10, MC3-33, MC3-37, MC3-60) at the LaSIE laboratory of La Rochelle University (France) using a Horiba High Resolution Raman spectrometer (LabRAM HR Evolution) equipped with a confocal microscope and a Peltier-based cooled charge coupled device (CCD) detector. Excitation was provided by a solid-state diode pumped green laser at 532 nm. The spot under the × 100 lens has a diameter of ~3 µm. Spectra were recorded with the acquisition LabSpec6 software at room temperature with a resolution of approximately 2 cm⁻¹ and corresponded to an accumulation of 2 × 60 s acquisitions. The laser power was reduced to 0.03 mW in order to prevent the transformation of Fe-compounds into hematite (α-Fe$_2$O$_3$) due to an excessive heating.
Whole-rock mineralogy of MC samples was analyzed at the ISTE-UNIL by X-Ray Diffraction (XRD), using a
Thermo Scientific ARL X-TRA diffractometer. The whole-rock mineralogy was determined by a semiquantitative
method, using XRD peak intensities of the main minerals compared to external standards¹⁰¹.

For optical microscopy observations whole rock polished pellets of selected samples were prepared in accord-
ance with standard procedures (ISO 7404-2, 2009). Samples were embedded in a mixture of resin and hardener
and polished after drying using a Buehler grinder-polisher. Petrographic observations illustrated in Figure S6-D
were carried out at Institute of Earth Sciences—Porto (Portugal) in a Leica 4000-M microscope equipped with
a Discus-Fossil system under standard conditions, and a 50 × oil immersion objective.

Received: 4 October 2021; Accepted: 3 March 2022

Published online: 14 March 2022

References

1. Suan, G., Mattioli, E., Pittet, B., Mailliot, S. & Lecuyer, C. Evidence for major environmental perturbation prior to and during the Toarcian (Early Jurassic) oceanic anoxic event from the Lusitanian Basin, Portugal. *Paleoceanography* https://doi.org/10.1029/2007pa001459 (2008).

2. Bailey, T. R., Rosenthal, Y., McArthur, J. M., van de Schootbrugge, B. & Thirlwall, M. F. Paleooceanographic changes of the Late Pliensbachian-Early Toarcian interval: A possible link to the genesis of an Oceanic Anoxic Event. *Earth Planet Sci. Lett.* **212**, 307–320 (2003).

3. Jenkyns, H. C. The early Toarcian (Jurassic) anoxic event—stratigraphic, sedimentary, and geochemical evidence. *Am. J. Sci.* **288**, 101–151 (1988).

4. Müller, T. et al. New multiproxy record of the Jenkyns Event (also known as the Toarcian Oceanic Anoxic Event) from the Mecsek Mountains (Hungary): Differences, duration and drivers. *Sedimentology* **64**, 66–86 (2017).

5. Reolid, M., Mattioli, E., Duarte, L. V. & Marok, A. The Toarcian oceanic anoxic event and the Jenkyns event (IGCP-655 final report). *Episodes* **43**, 833–844 (2020).

6. Jenkyns, H. C. Geochemistry of oceanic anoxic events. *Geochem. Geophys. Geosyst.* **11**, 2 (2010).

7. Hesselbo, S. P. et al. Massive dissociation of gas hydrate during a Jurassic oceanic anoxic event. *Nature* **406**, 392–395. https://doi.org/10.1038/35019044 (2000).

8. Kemp, D. B., Coe, A. L., Cohen, A. S. & Schwark, L. Astronomical pacing of methane release in the Early Jurassic period. *Nature* **437**, 396–399 (2005).

9. Fantasia, A. et al. The Toarcian Oceanic Anoxic Event in southwestern Gondwana: An example from the Andean Basin, northern Chile. *J. Geol. Soc. Lond.* **175**, 883–902 (2018).

10. French, K. L., Sepulveda, J., Trabucho-Alexandre, J., Grocke, D. R. & Summons, R. E. Organic geochemistry of the early Toarcian oceanic anoxic event in Hawsker Bottoms, Yorkshire, England. *Earth Planet Sci. Lett.* **390**, 116–127. https://doi.org/10.1016/j.epsl.2013.12.033 (2014).

11. van Breugel, Y., Baas, M., Schouten, S., Mattioli, E. & Damste, J. S. S. Isorenieratane record in black shales from the Paris Basin, France: Constraints on recycling of respired CO2 as a mechanism for negative carbon isotope shifts during the Toarcian oceanic anoxic event. *Paleoceanography* **21**, 2. https://doi.org/10.1029/2006pa001305 (2006).

12. Ruebsam, W., Mayer, B. & Schwark, L. Cryosphere carbon dynamics control early Toarcian global warming and sea level evolution. *Global Planet Change* **172**, 440–453. https://doi.org/10.1016/j.gloplacha.2018.11.003 (2019).

13. McElwain, J. C., Wade-Murphy, J. & Hesselbo, S. P. Changes in carbon dioxide during an oceanic anoxic event linked to intrusion into Gondwana coals. *Nature* **435**, 479–482 (2005).

14. Hermoso, M. et al. Dynamics of a stepped carbon-isotope excursion: Ultra high-resolution study of Early Toarcian environmental change. *Earth Planet Sci. Lett.* **319**, 45–54. https://doi.org/10.1016/j.epsl.2011.12.021 (2012).

15. Brazier, J. M. et al. Calcium isotope evidence for dramatic increase of continental weathering during the Toarcian oceanic anoxic event (Early Jurassic). *Earth Planet Sci. Lett.* **411**, 164–176 (2015).

16. Duncan, R. A., Hooper, P. R., Rehacek, J., Marsh, J. S. & Duncan, A. R. The timing and duration of the Karoo igneous event, southern Gondwana. *J. Geophys. Res.-Solid Earth* **102**, 18127–18138 (1997).

17. Palfy, J. & Smith, P. L. Synchrony between Early Jurassic extinction, oceanic anoxic event, and the Karoo-Ferrat flood basalt volcanism. *Geology* **28**, 747–750 (2000).

18. Jourdan, F., Feraud, G., Bertrand, H., Watkeys, M. K. & Renne, P. R. The Ar-40/Ar-39 ages of the sill complex of the Karoo large igneous province: Implications for the Pliensbachian-Toarcian climate change. *Geochem. Geophys. Geosyst.* **9**, 2. https://doi.org/10.1029/2008gc001994 (2008).

19. Svensen, H., Corfu, F., Polteau, S., Hammer, O. & Planke, S. Rapid magma emplacement in the Karoo Large Igneous Province. *Earth Planet Sci. Lett.* **325**, 1–9 (2012).

20. Moulin, M. et al. Eruptive history of the Karoo lava flows and their impact on early Jurassic environmental change. *J. Geophys. Res.-Solid Earth* **122**, 738–772. https://doi.org/10.1002/2016JB013354 (2017).

21. Suan, G. et al. Polar record of Early Jurassic massive carbon injection. *Earth Planet Sci. Lett.* **312**, 102–113 (2011).

22. Svensen, H. et al. Hydrothermal venting of greenhouse gases triggering Early Jurassic global warming. *Earth Planet Sci. Lett.* **256**, 554–566. https://doi.org/10.1016/j.epsl.2007.02.013 (2007).

23. Them, T. R. et al. High-resolution carbon isotope records of the Toarcian Oceanic Anoxic Event (Early Jurassic) from North America and implications for the global drivers of the Toarcian carbon cycle. *Earth Planet Sci. Lett.* **459**, 118–126 (2017).

24. Hermoso, M., Le Callonnec, L., Minoletti, F., Renard, M. & Hesselbo, S. P. Expression of the Early Toarcian negative carbon-isotope excursion in separated carbonate microfractions (Jurassic, Paris Basin). *Earth Planet Sci. Lett.* **277**, 194–203 (2009).

25. Them, T. R. et al. Evidence for rapid weathering response to climatic warming during the Toarcian Oceanic Anoxic Event. *Sci. Rep.-Uk* **7**, 2 (2017).

26. Cohen, A. S., Coe, A. L., Harding, S. M. & Schwark, L. Osmium isotope evidence for the regulation of atmospheric CO2 by continental weathering. *Geology* **32**, 157–160 (2004).

27. Percival, L. M. E. et al. Osmium isotope evidence for two pulses of increased continental weathering linked to Early Jurassic volcanism and climate change. *Geology* **44**, 759–762. https://doi.org/10.1130/G37997.1 (2016).

28. van Acken, D., Tutken, T., Daly, J. S., Schmid-Rohl, A. & Orr, P. J. Rhenium-osmium geochronology of the Toarcian Posidonia Shale SW Germany. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **534**, 2. https://doi.org/10.1016/j.palaeo.2019.109294 (2019).

29. Izumi, K., Kemp, D. B., Itamiya, S. & Inui, M. Sedimentary evidence for enhanced hydrological cycling in response to rapid carbon release during the early Toarcian oceanic anoxic event. *Earth Planet Sci. Lett.* **481**, 162–170 (2018).

30. Montero-Serrano, J. C. et al. Continental weathering and redox conditions during the early Toarcian Oceanic Anoxic Event in the northwestern Tethys: Insight from the Posidonia Shale section in the Swiss Jura Mountains. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **429**, 83–99 (2015).
31. Fantasia, A. et al. Global versus local processes during the Pliensbachian-Toarcian transition at the Peniche GSSP, Portugal: A multi-proxy record. *Earth-Sci. Rev.* **198**, 2 (2019).

32. Wignall, P. B., Newton, R. J. & Little, C. T. S. The timing of paleoenvironmental change and cause-and-effect relationships during the early Jurassic mass extinction in Europe. *Am. J. Sci.* **305**, 1014–1032. https://doi.org/10.2475/ajs.305.10.1014 (2005).

33. Bodin, S. et al. Toarcian carbon isotope shifts and nutrient changes from the Northern margin of Gondwana (High Atlas, Morocco, Jurassic): Palaeoenvironmental implications. *Palaeogeogr., Palaeoclimatol. Palaeoecol.* **297**, 377–390 (2010).

34. Boulila, S. & Hinnov, L. A. A review of tempo and scale of the early Jurassic Toarcian OAE: Implications for carbon cycle and sea level variations. *Newsl. Stratigr.* **50**, 363–389 (2017).

35. Boulila, S., Galbrun, B., Sadki, D., Gardin, S. & Bartolini, A. Constraints on the duration of the early Toarcian T-OAE and evidence for carbon-reservoir change from the High Atlas (Morocco). *Global Planet Change* **175**, 113–128 (2019).

36. Suan, G. et al. Duration of the Early Toarcian carbon isotope excursion deduced from spectral analysis: Consequence for its possible causes. *Earth Planet Sci. Lett.* **267**, 666–679. https://doi.org/10.1016/j.epsl.2007.12.017 (2008).

37. Jenkyns, H. C., Grocke, D. R. & Hesselbo, S. P. Nitrogen isotope evidence for water mass denitrification during the early Toarcian (Jurassic) oceanic anoxic event. *Paleoceanography* **16**, 593–603. https://doi.org/10.1029/2000pa000558 (2001).

38. Xu, W. M. et al. Evolution of the Toarcian (Early Jurassic) carbon-cycle and global climatic controls on local sedimentary processes (Cardigan Bay Basin, UK). *Earth Planet Sci. Lett.* **484**, 396–411 (2018).

39. Kemp, D. B. & Izumi, K. Multiproxy geochemical analysis of a Panthalassic margin record of the early Toarcian oceanic anoxic event (Toyora area, Japan). *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **414**, 332–341 (2014).

40. Kemp, D. B., Coe, A. L., Cohen, A. S. & Weedon, G. P. Astronomical forcing and chronology of the early Toarcian (Early Jurassic) oceanic anoxic event in Yorkshire, UK. *Paleoceanography* **26**, 2. https://doi.org/10.1029/2011pa002122 (2011).

41. Hesselbo, S. P. & Pienkowski, G. Stepwise atmospheric carbon-isotope excursion during the Toarcian Oceanic Anoxic Event (Early Jurassic, Polish Basin). *Earth Planet Sci. Lett.* **301**, 365–372. https://doi.org/10.1016/j.epsl.2010.11.021 (2011).

42. da Rocha, R. B. et al. Base of the Toarcian Stage of the Lower Jurassic defined by the Global Boundary Stratotype Section and Point (GSSP) at the Peniche section (Portugal). *Episodes* **39**, 460–481 (2016).

43. Duarte, L. V. Lithostratigraphy, sequence stratigraphy and depositional setting of the Pliensbachian and Toarcian series in the Lusitanian Basin (Portugal). *In R. B. Rocha (Ed.) – The Peniche section (Portugal). Contributions to the definition of the Toarcian GSSP. International Subcommission on Jurassic Stratigraphy* ISBN **978-972-8893-14-9**, 17–23 (2007).

44. Pittet, B., Suan, G., Lenoir, F., Duarte, L. V. & Mattioli, E. Carbon isotope evidence for sedimentary discontinuities in the lower Toarcian of the Lusitanian Basin (Portugal): Sea level change at the onset of the Oceanic Anoxic Event. *Sediment. Geol.* **303**, 1–14. https://doi.org/10.1016/j.sedgeo.2014.01.001 (2014).

45. Egli, R. Characterization of individual rock magnetic components by analysis of remanence curves, 1 Unmixing natural sediments. *Stud. Geophys. Geod.* **48**, 391–446. https://doi.org/10.1023/B:Sgeg.0000020839.45304.6d (2004).

46. Egli, R., Chen, A. P., Winklhofer, M., Kodama, K. P. & Horng, C. S. Detection of noninteracting single domain particles using first-order reversal curve diagrams. *Geochem. Geophys. Geosyst.* https://doi.org/10.1029/2009gc002916 (2010).

47. Ludwig, P. et al. Characterization of primary and secondary magnetite in marine sediment by combining chemical and magnetic unmixing techniques. *Global Planet Change* **110**, 321–339. https://doi.org/10.1016/j.gloplacha.2013.08.018 (2013).

48. Lascu, I. et al. Magnetic unmixing of first-order reversal curve diagrams using principal component analysis. *Geochem. Geophys. Geosyst.* **16**, 2900–2915 (2015).

49. Roberts, A. P., Chang, L., Heslop, D., Florindo, F. & Larrasoana, J. C. Searching for single domain magnetite in the "pseudo-single-domain" sedimentary haystack: Implications of biogenic magnetite preservation for sediment magnetism and relative paleointensity determinations. *J. Geophys. Res. Solid Earth* **117**, 2 (2012).

50. Ruhl, M. et al. Astronomical constraints on the duration of the early Jurassic Hettangian stage and recovery rates following the end-Triassic mass extinction (St Audrie's Bay/East Quantoxhead, UK). *Earth Planet Sci. Lett.* **295**, 262–276. https://doi.org/10.1016/j.epsl.2010.04.008 (2010).

51. Courtillot, V. E. & Renne, P. R. On the ages of flood basalt events. *Cr Geosci.* **335**, 113–140. https://doi.org/10.1016/S1631-0713(03)00006-3 (2003).

52. Burgess, S. D., Bowring, S. A., Fleming, T. H. & Elliot, D. H. High-precision geochronology links the Ferrar large igneous province with early-Jurassic ocean anoxia and biotic crisis. *Earth Planet Sci. Lett.* **415**, 90–99 (2015).

53. Ivanov, A. V. et al. Timing and genesis of the Karoo-Ferrar large igneous province: New high precision U-Pb data for Tasmania confirm short duration of the major magmatic pulse. *Chem. Geol.* **455**, 32–43 (2017).

54. Sell, B. et al. Evaluating the temporal link between the Karoo LIP and climatic–biologic events of the Toarcian stage with high-precision U-Pb geochronology. *Earth Planet Sci. Lett.* **408**, 48–56 (2014).

55. Burgess, S. D. & Bowring, S. A. High-precision geochronology confirms voluminous magmatism before, during, and after Earth's most severe extinction. *Sci. Adv.* **1**, 2 (2015).

56. Greber, N. D. et al. New high precision U-Pb ages and Hf isotope data from the Karoo large igneous province; implications for pulsed magmatism and early Toarcian environmental perturbations. *Res. Geochem. 1*, 100005. https://doi.org/10.1016/j.ringeo.2020.100005 (2020).

57. Grasby, S. E., Them, T. R., Chen, Z. H., Yin, R. S. & Ardakani, O. H. Mercury as a proxy for volcanic emissions in the geologic record. *Earth-Sci. Rev.* **196**, 2 (2019).

58. Font, E. & Bond, D. Volcanism and Mass Extinction. In *Encyclopedia of Geology* 2nd edn (eds Alderton, D. & Elias, S. A.) 596–606 (Academic Press, 2021).

59. Font, E. et al. Mercury anomaly, Deccan volcanism, and the end-Cretaceous mass extinction. *Geology* **44**, 171–174. https://doi.org/10.1130/g37451.1 (2016).

60. Sanei, H., Grasby, S. E. & Beauchamp, B. Latest Permian mercury anomalies. *Geology* **40**, 63–66. https://doi.org/10.1130/G32596.1 (2012).

61. Sial, A. N. et al. Mercury as a proxy for volcanic activity during extreme environmental turnover: The Cretaceous-Paleogene transition. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **387**, 153–164. https://doi.org/10.1016/j.palaeo.2013.07.019 (2013).

62. Percival, L. M. E. et al. Globally enhanced mercury deposition during the end-Pliensbachian extinction and Toarcian OAE: A link to the Karoo-Ferrar Large Igneous Province. *Earth Planet Sci. Lett.* **428**, 267–280 (2015).

63. Rodrigues, B., Duarte, L. V., Silva, R. L. & Mendonca, J. G. Sedimentary organic matter and early Toarcian environmental changes in the Lusitanian Basin (Portugal). *Palaeogeogr., Palaeoclimatol. Palaeoecol.* **554**, 2 (2020).

64. Them, T. R. et al. Terrestrial sources as the primary delivery mechanism of mercury to the oceans across the Toarcian Oceanic Anoxic Event (Early Jurassic). *Earth Planet Sci. Lett.* **507**, 62–72. https://doi.org/10.1016/j.epsl.2018.11.029 (2019).

65. Dickens, G. R., Oneil, J. R., Rea, D. K. & Owen, R. M. Dissociation of oceanic methane hydrate as a cause of the carbon-isotope excursion at the end of the paleocene. *Paleoceanography* **10**, 965–971. https://doi.org/10.1029/95pa02087 (1995).

66. Edwards, B. A., Kushner, D. S., Outridge, P. M. & Wang, F. Y. Fifty years of volcanic mercury emission research: Knowledge gaps and future directions. *Sci. Total Environ.* **757**, 2 (2021).

67. Zintwana, M. P., Cawthorn, R. G., Ashwal, L. D., Roelofse, F. & Cronwright, H. Mercury in the Bushveld Complex, South Africa, and the Skaergaard Intrusion Greenland. *Chem. Geol.* **320**, 147–155 (2012).

68. Canil, D., Crockford, P. W., Rossin, R. & Telmer, K. Mercury in some arc crustal rocks and mantle peridotites and relevance to the moderately volatile element budget of the Earth. *Chem. Geol.* **396**, 134–142 (2015).
69. Coufalik, P. et al. Model of mercury flux associated with volcanic activity (vol 101, pg 549, 2018). *B Environ. Contam. Tox.* **101**, 554–555 (2018).

70. Hunt, J. M. Petroleum Geochemistry and Geology, W. H. Freeman and Company, San Francisco, California. 617 (1979).

71. Luchsheva, L., Konovalov, Y. & Obzhirov, A. Gas hydrate deposits are a real source of mercury in marine ecosystems. *Norwegian J. Dev. Int. Sci.* **6**, 10–14 (2017).

72. Boulila, S. et al. Astronomical calibration of the Toarcian Stage: Implications for sequence stratigraphy and duration of the early Toarcian OAE. *Earth Planet Sci. Lett.* **386**, 98–111 (2014).

73. Lunt, D. J. et al. A model for orbital pacing of methane hydrate destabilization during the Palaeogene. *Nat. Geosci.* **4**, 775–778 (2011).

74. Heimdal, T. H., Godderis, Y., Jones, M. T. & Svensen, H. H. Assessing the importance of thermogenic degassing from the Karoo Large Igneous Province (LIP) in driving Toarcian carbon cycle perturbations. *Nat. Commun.* **12**, 2 (2021).

75. Otake, T., Wesolowski, D. J., Anovitz, L. M., Allard, L. F. & Ohmoto, H. Experimental evidence for non-redox transformations between magnetite and hematite under H-2-rich hydrothermal conditions. *Earth Planet Sci. Lett.* **257**, 60–70 (2007).

76. Roberts, A. P. Magnetic mineral diagenesis. *Earth Sci. Rev.* **151**, 1–47 (2015).

77. Orberger, B. et al. Micro- to nano-scale characterization of martite from a banded iron formation in India and a lateritic soil in Brazil. *Phys. Chem. Miner.* **41**, 651–667 (2014).

78. Ohmoto, H. Nonredox transformations of magnetite-hematite in hydrothermal systems. *Econ. Geol. Bull. Soc.* **98**, 157–161 (2003).

79. Passier, H. F., de Lange, G. J. & Dekkers, M. J. Magnetic properties and geochemistry of the active oxidation front and the youngest sapropel in the eastern Mediterranean Sea. *Geophys. J. Int.* **145**, 604–614 (2001).

80. Larrasoana, J. C., Roberts, A. P., Stoner, J. S., Richter, C. & Wehausen, R. A new proxy for bottom-water ventilation in the eastern Mediterranean based on diagenetically controlled magnetic properties of sapropel-bearing sediments. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **190**, 221–242. https://doi.org/10.1016/S0031-0182(02)00607-7 (2003).

81. Passier, H. F. & Dekkers, M. J. Iron oxide formation in the active oxidation front above sapropel S1 in the eastern Mediterranean Sea as derived from low-temperature magnetism. *Geophys. J. Int.* **150**, 230–240 (2002).

82. Kruiver, P. P. & Passier, H. F. Coercivity analysis of magnetic phases in sapropel S1 related to variations in redox conditions, including an investigation of the S ratio. *Geochem. Geophys. Geosyst.* **2**, 2 (2001).

83. Larrasoana, J. C., Roberts, A. P., Hayes, A., Wehausen, R. & Rohling, E. J. Detecting missing beats in the Mediterranean climate rhythm from magnetic identification of oxidized sapropels (Ocean Drilling Program Leg 160). *Phys. Earth Planet. Inter.* **156**, 283–293 (2006).

84. De Lange, G. J. et al. Synchronous basin-wide formation and redox-controlled preservation of a Mediterranean sapropel. *Nat. Geosci.* **1**, 606–610 (2008).

85. van Santvoort, P. J. M. et al. Active post-depositional oxidation of the most recent sapropel (S1) in sediments of the eastern Mediterranean Sea. *Geochim. Cosmochim.* **Ac** **60**, 4007–4024 (1996).

86. Larrasoana, J. C., Roberts, A. P., Rohling, E. J., Winklhofer, M. & Wehausen, R. Three million years of monsoon variability over the northern Sahara. *Clim. Dynam.* **21**, 689–698 (2003).

87. Duarte, L. V. Facies analysis and sequential evolution of the Toarcian-lower Aalenian series in the Lusitanian Basin (Portugal). *Comunicações do Instituto Geológico e Mineiro* **83**, 65–94 (1997).

88. Wright, V. P. & Wilson, R. C. L. A carbonate submarine-fan sequence from the Jurassic of Portugal. *J. Sediment. Petrol.* **54**, 394–412 (1984).

89. Maher, B. A., Alekseev, A. & Alekseeva, T. Magnetic mineralogy of soils across the Russian Steppe: climatic dependence of pedogenic magnetite formation. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **210**, 321–341 (2003).

90. Dera, G. et al. Distribution of clay minerals in Early Jurassic Peritethyan seas: Palaeoclimatic significance inferred from multi-proxy comparisons. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **271**, 39–51 (2009).

91. Rodrigues, B., Duarte, L. V., Mendonca, J. G., Santos, L. G. & de Oliveira, A. D. Evidence of terrestrial organic matter deposition across the early Toarcian recorded in the northern Lusitanian Basin Portugal. *Int. J. Coal. Geol.* **168**, 35–45 (2016).

92. Trecalli, A., Spangenberg, J., Adatte, T., Follmi, K. B. & Parente, M. Carbonate platform evidence of ocean acidification at the onset of the early Toarcian oceanic anoxic event. *Earth Planet Sci. Lett.* **357**, 214–225. https://doi.org/10.1016/j.epsl.2012.09.043 (2012).

93. Ettinger, N. P. et al. Ocean acidification and photic-zone anoxia at the Toarcian oceanic anoxic event: Insights from the Adriatic carbonate platform. *Sedimentology* **68**, 63–107 (2021).

94. Müller, T. et al. Ocean acidification during the early Toarcian extinction event: Evidence from boron isotopes in brachiopods. *Geology* **48**, 1184–1188. https://doi.org/10.1130/G47781.1 (2020).

95. Ruebsam, W., Reolid, M. & Schwark, L. delta C-13 of terrestrial vegetation records Toarcian CO2 and climate gradients. *Sci. Rep.-Uk* **10**, 2 (2020).

96. Kruiver, P. P., Dekkers, M. J. & Heslop, D. Quantification of magnetic coercivity components by the analysis of acquisition curves of isothermal remanent magnetisation. *Earth Planet Sci. Lett.* **189**, 269–276. https://doi.org/10.1016/s0012-821x(01)00367-3 (2001).

97. Maxbauer, D. P., Feinberg, J. M. & Fox, D. L. MAX UnMix: A web application for unmixing magnetic coercivity distributions. *Comput. Geosci.-Uk* **95**, 140–145. https://doi.org/10.1016/j.cageo.2016.07.009 (2016).

98. Egli, R. VARIFORC: An optimized protocol for calculating non-regular first-order reversal curve (FORC) diagrams. *Global Planet Change* **110**, 302–320. https://doi.org/10.1016/j.gloplacha.2013.08.003 (2013).

99. Egli, R. Magnetic characterization of geologic materials with first-order reversal-curves. In *Magnetic Measurement Techniques for Materials Characterization* (eds Franco, V. & Dodrill, B.) (Springer Nature Publishing Group, 2021).

100. Behar, F., Beaumont, V. & Penteado, H. L. D. Rock-Eval 6 technology: Performances and developments. *Oil Gas Sci. Technol.* **56**, 111–134 (2001).

101. Adatte, T., Stinnesbeck, W. & Keller, G. Lithologic and mineralogic correlations of near K/T boundary clastic sediments in NE Mexico: Implication for origin and nature of deposition. *Spec. Publ. Geol. Soc. Am.* **307**, 211–226 (1996).

102. Scotese, C.R. Paleomap Project, www.scotese.com. (2013).

103. Hesselbo, S. P., Jenkyns, H. C., Duarte, L. V. & Oliveira, L. C. V. Carbon-isotope record of the Early Jurassic (Toarcian) oceanic anoxic event from fossil wood and marine carbonate (Lusitanian Basin, Portugal). *Earth Planet Sci. Lett.* **253**, 455–470. https://doi.org/10.1016/j.epsl.2006.11.009 (2007).

## Acknowledgements

This work was financially supported by the project FCT/UIDB/50019/2020 - IDL funded by FCT. Magnetic analyses was supported by EPOS-TNA-MSL 2018-001. We thank Glen Snyder, Matias Reolid and Peter Barry for their fruitful reviews.
**Author contributions**

E.F. conducted sampling and rock magnetic analysis, analysed the results and wrote the manuscript. L.V.D. conducted sampling and field work. M.J.D. conducted rock magnetic analysis. C-R. conducted micro-Raman analysis. R.E. conducted FORC diagram. J.E.S conducted carbon isotope analysis. A.F. contributed to sampling and field work. J.R. conducted petrographic analysis of organic matter. E.G. conducted petrographic observations. J.M. conducted SEM observations. T.A. provided mineralogical and geochemical data. All authors reviewed and contributed to the manuscript.

**Funding**

Fundação Ciencia e Tecnologia (FCT/UIDB/50019/2020—IDL).

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-08269-y.

**Correspondence** and requests for materials should be addressed to E.F.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022