# Climate-controlled organic matter accumulation as recorded in the Upper Jurassic Argiles de Châtillon Formation, a shallow-marine counterpart of the Kimmeridge Clay Formation

I. van der Hoeven, R. Verreussel, A. Riboulleau, Nicolas Tribovillard, B. van de Schootbrugge

► To cite this version:

I. van der Hoeven, R. Verreussel, A. Riboulleau, Nicolas Tribovillard, B. van de Schootbrugge. Climate-controlled organic matter accumulation as recorded in the Upper Jurassic Argiles de Châtillon Formation, a shallow-marine counterpart of the Kimmeridge Clay Formation. Geological Magazine, 2023, 160 (3), pp.579-600. 10.1017/s0016756822001121. hal-04219714

HAL Id: hal-04219714

https://hal.science/hal-04219714v1

Submitted on 9 Feb 2024

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
## Original Article

**Cite this article:** van der Hoeven IC, Verreussel RMCH, Riboulleau A, Tribovillard N, and van de Schootbrugge B (2023) Climate-controlled organic matter accumulation as recorded in the Upper Jurassic Argiles de Châtillon Formation, a shallow-marine counterpart of the Kimmeridge Clay Formation. *Geological Magazine* **160**: 579–600. https://doi.org/10.1017/S0016756822001121

Received: 24 May 2022
Revised: 29 September 2022
Accepted: 3 October 2022
First published online: 2 December 2022

**Keywords:** carbon isotopes; geochemistry; Jurassic; Milankovitch cyclicity; palynology; sulfurization

**Author for correspondence:** B. van de Schootbrugge, Email: B.vanderSchootbrugge@uu.nl

*Current address: Department of Ocean Systems, NIOZ Royal Netherlands Institute for Sea Research, Landsdiep 4, 1797 SZ 't Horntje, The Netherlands

# Climate-controlled organic matter accumulation as recorded in the Upper Jurassic Argiles de Châtillon Formation, a shallow-marine counterpart of the Kimmeridge Clay Formation

I. C. van der Hoeven¹,² Ⓓ, R. M. C. H. Verreussel¹, A. Riboulleau³, N. Tribovillard³ and B. van de Schootbrugge²,* Ⓓ

¹TNO Geological Survey of the Netherlands, Princetonlaan 6, 3584 CB, Utrecht, The Netherlands; ²Department of Earth Sciences, Utrecht University, Princetonlaan 8A, 3584CB, Utrecht, The Netherlands and ³Laboratory of Oceanology & Geosciences, Université de Lille, UMR LOG 8187 Univ Lille-CNRS-ULCO-IRD, 59000, Lille, France

## Abstract

Mudstones from the Argiles de Châtillon Formation exposed in the Boulonnais region of Northern France represent a proximal lateral equivalent of the organic-rich Kimmeridge Clay Formation. The Argiles de Châtillon Formation is composed of two subunits that straddle the Kimmeridgian–Tithonian boundary. Each subunit contains an organic-rich interval. The two conspicuous organic-rich intervals have been linked to either periods of high sea level or greenhouse warming. Here, we use palynology to further understand climate and environmental mechanisms that drove organic matter enrichment. We use bulk organic carbon isotope records (δ¹³Corg) to correlate the Boulonnais sections with those of the Kimmeridge Clay Formation. The palynological results suggest that the stratigraphically lower organic-rich interval (Kimmeridgian) was deposited under suboxic to anoxic stratified conditions. A large-scale climate shift from cooler/humid to warmer/arid conditions marked the Kimmeridgian–Tithonian boundary, influencing organic matter enrichment in the stratigraphically higher organic-rich interval (Tithonian). In contrast with the lower organic-rich interval, there are no indications of stratified conditions for the higher organic-rich interval. Within this thicker organic-rich interval, cyclic variations in amorphous organic matter distribution, total organic carbon and δ¹³Corg trends on a 2 m scale are observed. They co-occur with fluctuations of the palynological assemblages, indicative of more humid versus arid climate conditions, likely alternating on a ~100 kyr eccentricity timescale. Our results show that under the most humid phases of these overall arid climate conditions, sulfurization of carbohydrates was the dominant control on organic matter preservation. This climate-controlled process that drives organic matter enrichment in the Tithonian can be recognized on a basin-wide scale.

## 1. Introduction

The coast near Boulogne-sur-Mer in NW France offers a well-exposed succession of Upper Jurassic sandstone units alternating with organic-rich mudstones. This study focuses on one of these organic-rich mudstone units, the Argiles de Châtillon Formation (ACF). The ACF is a proximal analogue for the Kimmeridge Clay Formation (KCF) that crops out in southern England and is an important source rock in the North Sea Basin (Herbin et al. 1995; Proust et al. 1995; Tribovillard et al. 2005; Cornford, 2018). The KCF consists predominantly of laminated mudstones, intercalated with limestones and dolostones, indicating a relatively distal marine depositional environment. The less complete and more proximal ACF is characterized by more bioturbated mudstones, thin sandstones and storm-induced bioclastic limestones. Based on ammonite biostratigraphy, the age of the ACF is latest Kimmeridgian to earliest Tithonian (Geyssant et al. 1993; Herbin et al. 1995), which corresponds to the middle part of the KCF.

The ACF is stratigraphically sandwiched between two prominent sandstone units, the Grès de Châtillon Formation below and the Grès de la Crèche Formation above (Fig. 1). Halfway through the 20 m thick succession of the ACF, an erosive unconformity, interpreted as a sequence boundary (Wignall & Newton, 2001), subdivides the formation into a lower and an upper member (Fig. 1). Two distinct organic-rich intervals (ORIs) occur in the ACF, one near the base of the lower member and one in the lower part of the upper member. These display total organic carbon (TOC) contents of a maximum of 8 %, but TOC contents are laterally variable (Herbin et al. 1995; Tribovillard et al. 2001). Sequence stratigraphic studies allow the correlation of the two ORIs to periods of rising and high sea level (Proust et al. 2001; Wignall &

© The Author(s), 2022. Published by Cambridge University Press. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.

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

CAMBRIDGE
UNIVERSITY PRESS
<figure><img src="image_3.png" /></figure>

(1) Proust et al. (1995); (2) Proust et al. (2001); (3) Herbin et al. (1995);

(4) Braaksma et al. (2006); (5) Tribovillard et al. (2005)

<caption>**Fig. 1.** Stratigraphic framework of the Upper Jurassic succession from the Boulonnais area, displaying the lithostratigraphy, chronostratigraphy sequence stratigraphy, and organic matter rich intervals.</caption>

Newton, 2001). However, according to others, sea level is not the sole factor controlling organic matter (OM) enrichment. Changes in global climate, driving changes in primary productivity, may have played an equally important role (Herbin et al. 1995). Moreover, it was recently proposed that the two ORIs in the ACF are the result of different driving forces, with the lower interval deposited under calm and restricted conditions, and the upper interval deposited under oxic and more energetic conditions (Tribovillard et al. 2019). The apparent contradiction of OM enrichment under oxic conditions is explained by assuming a system of storm-induced OM concentration. According to these authors, sulfurization occurred in more proximal environments and storms remobilized the OM to the oxygenated, shallow-marine setting. The chemically recalcitrant structure of the sulfurized OM made it resistant against oxidation after deposition.

The ACF spans the transition from the Kimmeridgian to the Tithonian (Gallois et al. 2019), a period that records a climate change in NW Europe, evolving from warm and humid conditions in late Kimmeridgian time (Eudoxus and Autissiodorensis ammonite zones) to more arid conditions in early Tithonian time (Elegans ammonite Zone; Abbink et al. 2001; Deconinck & Baudin, 2008; Hesselbo et al. 2009). Smaller-scale Milankovitch cyclicity-controlled climate variations have been recorded in the ACF as well as in the KCF, with alternations on a decimetre-metre scale (e.g. Waterhouse, 1995, 1999; Weedon et al. 1999; Huang et al. 2010).

The present study aims to understand the relationship between variations in Late Jurassic climate, regional carbon isotope trends and OM distribution of the Upper Jurassic ACF. Because of its relatively proximal marine setting, the ACF is adequately suited for palynological analysis, as it provides a balanced mix between land-derived and marine palynomorphs. The land-derived palynomorphs are expected to reflect trends in vegetation and, by inference, climate. We aim to investigate if and how climate change controlled the OM distribution in the ACF. In addition, the ACF and the base of the overlying Grès de la Crèche Fm have been sampled in high resolution for stable isotope analysis (δ¹³Corg) in order to correlate the Boulonnais outcrops to the regional record from the KCF (Morgans-Bell et al. 2001). An accurate correlation between the ACF and KCF will allow comparison between OM accumulation processes in the proximal and distal domains and establish an OM accumulation model for the latest Kimmeridgian – earliest Tithonian period that is applicable on a basin-wide scale.

## 2. Geological setting

During Kimmeridgian and Tithonian times, Western Europe was composed of an archipelago of low-lying massifs set in a broad and relatively shallow epicontinental sea, the Laurasian Seaway (Ziegler, 1990; Fig. 2). The Laurasian Seaway connected the Panboreal Superrealm in the north to the Tethys Ocean to the south. This area comprises a series of interconnected basins, such as the Paris, Wealden, Cleveland and Wessex basins, that evolved as a result of the break-up of Pangaea that started in the Triassic Period (Underhill & Stoneley, 1998; Mansy et al. 2003; Zanella & Coward, 2003). Estimations of the water depth in these basins varied between several tens of metres (e.g. Hallam, 1975) and hundreds of metres (e.g. Herbin et al. 1991).

The sediments of the ACF were deposited at the easternmost edge of the Wessex Basin, bounded in the north by the London–Brabant Massif (Fig. 2). The Wessex Basin was periodically anoxic, alternating with shorter or longer oxygenation events (Oschmann, 1988; van Kaam-Peters et al. 1997; Waterhouse, 1999). As part of a clastic-dominated ramp environment, the Boulonnais area accumulated sand-dominated deposits alternating with silt and clay-dominated deposits. Sedimentation was controlled by relative sea level (Proust et al. 1995; Braaksma et al. 2006) and affected by synsedimentary faulting (Hatem et al. 2014). Multiple studies of the Boulonnais sections have provided a depositional and sequence stratigraphical framework that is useful to correlate and compare the different outcrops and to explain variations in OM content (Herbin et al. 1995; Proust et al. 1995; Deconinck et al. 1996; Wignall et al. 1996; Tribovillard et al. 2001; Wignall & Newton, 2001; Braaksma et al. 2006; Deconinck & Baudin, 2008; Hatem et al. 2017).
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** (Colour online) Palaeogeographic map of NW Europe for the Late Jurassic period showing the locations of the Boulonnais area in relation to other basins and emergent land masses. B – Boulonnais; WB – Wessex Basin; PB – Paris Basin; CB – Cleveland Basin; BFB – Broad Fourteens Basin; CG – Central Graben; VG – Viking Graben; LBM – London–Brabant Massif; AM – Armorican Massif. Modified after Coward *et al.* (2003).</figcaption></figure>

The most complete and best accessible outcrops of the ACF and the underlying Grès de Châtillon Fm occur between the village of Audresselles and Cap Gris Nez (Fig. 3). The 5 m thick Grès de Châtillon Fm consists of thin carbonate cemented sandstones with clay partings, common cross-bedding, wave ripples and trace fosils such as *Rhizocorallium*, *Ophiomorpha* and *Thalassinoides* (Wignall *et al.* 1996; Schlirf, 2003). The sandstones were deposited in a tidally influenced, shallow-marine setting and are shallowing-upwards (Proust *et al.* 1995; Braaksma *et al.* 2006; Angus *et al.* 2020). A conspicuous 1 m thick uncemented and heterogeneous sandstone with bright orange colours marks the top of the formation and is laterally continuous over several kilometres (Fig. 5b). Crocodile teeth (pers. obs.), charcoal fragments, bioclastic gravel strings (Angus *et al.* 2020) and features indicating synsedimentary faulting (Hatem *et al.* 2014, 2016) occur in this layer.

The 18–25 m thick ACF consists of bioturbated mudstones and siltstones, limestones and laminated paper shale (Tribovillard *et al.* 2001; Wignall & Newton, 2001; Hatem *et al.* 2016). The base of the unit consists of 35–50 cm thick bioturbated sandy limestone, locally associated with oyster patch reefs. The oyster patch reefs and cone-in-cone structures, occurring on concretions in the ACF, indicate fluid expulsion at the time of deposition, related to synsedimentary faulting (Hatem *et al.* 2014, 2016; Tribovillard *et al.* 2018). The shales of the ACF generally represent a low-energy shelf facies, deposited below wave base. The lithological succession of the ACF differs slightly from outcrop to outcrop, with the main differences being variations in overall thickness and a laterally variable amount of storm beds (Wignall & Newton, 2001). Storm beds occur as thin, hummocky cross-stratified

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** (Colour online) Overview of the Boulonnais area and location of the three studied sections, the Cran aux Œufs, Cran du Noirda and Cap de la Crèche sections. For the palaeogeographic position of the Boulonnais in NW Europe/Wessex Basin, see Figure 2.</figcaption></figure>

(HCS) sandstones (Fig. 5d), or as shell-bearing limestones (Fig. 5e), typically consisting almost exclusively of the small oyster species *Nanogyra nana* (Tribovillard *et al.* 2019). These shell beds are referred to as coquina beds, and this oyster type mainly occurs in brackish bays and lagoons (Fürsich, 1981; Fürsich & Oschmann, 1986). Particularly around the Kimmeridgian–Tithonian boundary, storm beds occur commonly (Fig. 4). Apart from the storm beds, the overall hydrodynamic regime was quite low, based on the abundant occurrence of the shallow-water foraminifer genus *Ammobaculites* (Colpaert *et al.* 2021). The calcareous foraminifer
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** (Colour online) Lithological log of the Cran aux Œufs section. The subsections bearing sample codes OEB and OEF are correlated via a marker bed on top of a small unexposed interval, marked with a cross. The position of the stage, ammonite zone and member boundaries are based on the δ¹³Corg correlations. Outcrop pictures are displayed in Figure 5.</figcaption></figure>

genus *Lenticulina* dominates the foraminifera assemblages in the ORIs of the ACF, pointing to temporal hypoxia (Colpaert *et al.* 2021). The ORIs of the ACF have been compared with the more distal time-equivalent mudstones of the KCF (Herbin *et al.* 1991, 1995; Proust *et al.* 1993; Deconinck *et al.* 1996), but the TOC values, ranging from <1 to 8%, are much lower in comparison. The sulfur contents, on the other hand, are up to three times higher in the ACF than the sulfur contents at equivalent TOC

values for the KCF in Dorset (Tribovillard *et al.* 2001). The hydrogen index (HI) of the ORIs from the upper ACF varies between 50 and 550 mg HC/g TOC, with an extrapolated maximum of 600 mg HC/g TOC and exhibits a strong positive correlation (r² = 0.97) with the TOC (Tribovillard *et al.* 2001). The OM of the laminated intervals can be classified as oil-prone, type II-S kerogen, and with a relatively low average Tmax value of 423 °C, it is considered immature with respect to hydrocarbon generation.

At Cap Gris Nez, the ACF is succeeded by sandstones of the Grès de la Crèche Fm (Figs 1, 3). The Grès de la Crèche Fm displays yellow and brownish colours that contrast strongly with the dark grey shales of the ACF (Fig. 5a). The Grès de la Crèche Fm is 10 to 15 m thick and consists of a lower and an upper sandstone unit, separated by siltstones and mudstones of the Marnes Intercalaires (Proust *et al.* 1995, 2001). Cross-bedding occurs throughout the Grès de la Crèche Fm, as well as abundant bivalves and trace fossils (Proust *et al.* 1995; Schlirf, 2003). The most conspicuous feature of the Grès de la Crèche Fm is the occurrence of large spheroidal concretions, locally referred to as 'boules', often occurring in specific stratigraphic levels and in some places amalgamated into cemented layers (Fig. 5a). These concretions are the result of early diagenesis (Hatem *et al.* 2016; Tribovillard *et al.* 2018).

## 2.a. Late Jurassic climate in NW Europe

The Jurassic Period is generally considered to be a greenhouse world without any major glaciations (Donnadieu *et al.* 2011). Hothouse conditions in Early Jurassic time were induced by high atmospheric CO₂ levels as a result of volcanic outgassing related to the Central Atlantic Magmatic Province (Korte & Hesselbo, 2011). The greenhouse climate and high sea level favoured the development of tropical cyclones up to 40°N (Agustsdottir *et al.* 1999; Colombié *et al.* 2018). In Late Jurassic to Early Cretaceous times, the Boulonnais region occupied a palaeolatitude between 30 and 35°N under a warm, humid and possibly monsoonal climate (Hallam, 1984; Sellwood & Valdes, 2008; Hesselbo *et al.* 2009; Armstrong *et al.* 2016; Angus *et al.* 2020). The predominantly humid climate was interrupted by a stepwise aridification in Late Jurassic time, based on palynological (Abbink *et al.* 2001; Schneider *et al.* 2018) and clay mineral records (Hallam, 1993; Schnyder *et al.* 2006). Climate conditions in the latest Kimmeridgian *Eudoxus-Autissiodorensis* ammonite zones, on the other hand, are thought to have been relatively humid based on a comparison of clay mineral data (Hesselbo *et al.* 2009). The most humid conditions prevailed during the *Autissiodorensis* Zone, followed by a shift towards more arid conditions in the succeeding ammonite zones.

The change towards warmer and more arid conditions has been recorded in the palynological record from the UK, Denmark and the Netherlands and is referred to as the Scitulus climate shift (Abbink, 1998; Abbink *et al.* 2006; Verreussel *et al.* 2018). The climate shift towards drier conditions is reflected by the sudden increase in *Classopollis* at the base of the *Elegans* ammonite Zone (Abbink *et al.* 2001). The maximum aridity in NW France in the studied interval is reached in the *Scitulus* Zone (Deconinck *et al.* 1983). The warm and arid phase continued across the Jurassic-Cretaceous boundary and gradually changed into a more humid climate in Berriasian time (Abbink *et al.* 2001; Hesselbo *et al.* 2009; Schneider *et al.* 2018).

An increasing abundance of storm deposits around the Kimmeridgian stage in Western Europe has been described as a
<figure><img src="image_7.png" /><figcaption>**Fig. 5.** (Colour online) Outcrop pictures from the Cran aux Œufs section. (a) Subsection OEF, representing the upper part of the Cran aux Œufs section that covers the upper ACF. The base of subsection OEF is a distinct marker horizon consisting of a 75 cm thick shell bed. Two correlatable storm beds are indicated with arrows, the lower one being displayed in close-up in (e). The star marks the base of the Grès de la Crèche Fm. Note persons in circle for scale. (b) Subsection OEB, representing the lower part of the Cran aux Œufs section that covers the Grès de Châtillon Fm and the lower ACF. In the background and masked with a transparent overlay, the upper ACF is visible. Correlatable storm beds are indicated by arrows. The white star in the far background on the left marks the base of the Grès de Châtillon Fm, as portrayed in (a). (c) Up to 10 cm thick yellow-orange weathering level at the base of the lower ORI near the base of the ACF. (d) Thin hummocky cross-stratified (HCS) sandstone (upper ACF) near the base of the upper ORI. (e) Thick storm bed with coquina shells. Note the cross-stratified set at the top. (f) Laminated organic-rich mudstones, upper ACF, upper ORI. Pen for scale is 14 cm long.</figcaption></figure>
consequence of a change from a more humid to drier climate, coeval with deep seawater warming, that resulted in the expansion of the tropical zone (e.g. Hallam, 1985; Valdes & Sellwood, 1992; Baniak et al. 2014; Colombié et al. 2018). Recently, it was proposed that an expanded Hadley Cell, with an intensified but alternating hydrological cycle and associated increased storminess, heavily influenced sedimentation and TOC enrichment by promoting primary productivity and OM burial in the Boreal and Laurasian seaway (Armstrong et al. 2016; Colombié et al. 2018; Atar et al. 2019). Armstrong et al. (2016) argued that OM-lean sediments are deposited under relatively dry conditions while OM-rich sediments are deposited under monsoonal conditions, with increased freshwater and nutrient input, all taking place on a 100 kyr eccentricity timescale (Huang et al. 2010).

# 3. Materials and methods

## 3.a. Sections and samples

All sample material is derived from cliff outcrops of Upper Jurassic deposits located along the 25 km long coastline of the Boulonnais region (Fig. 3). Three different sets of samples are used in this study. The main sample set is derived from a field campaign in 2020 and includes 62 samples from the ACF (Fig. 4). These samples have been analysed for palnology and stable isotopes (Table 1; online Supplementary Material). The ACF was logged and sampled in the Cran aux Œufs section, situated in the north of the Boulonnais area (Fig. 3), GPS coordinates 50° 51′ 1.23" N and 1° 34′ 52.37" E. This location was favoured over other outcrops of the ACF because of its easy access and the lack of large apparent hiatuses (Wignall & Newton, 2001). Other outcrop localities of the ACF, such as Cran du Noirda or Cap de la Crèche, are frequently showcased in publications, but the ACF is described as being less complete (Wignall & Newton, 2001).

The Cran aux Œufs section covers the entire interval from the Grès de Châtillon Fm to the Grès de la Crèche Fm and includes the complete succession of the ACF. The section was sampled in two nearby locations, situated ~50 m apart (Fig. 5). The sub-section from the southernmost outcrop locality is coded OEB and covers 10 ~ m of the basal part of the ACF (Fig. 5b), including the boundary with the underlying Grès de Châtillon Fm (Figs 4, 5b). The other outcrop locality is coded OEF and covers the upper 15 m of the ACF, including the boundary with the overlying Grès de la Crèche Fm (Fig. 5a). The two subsections are easily correlated via a laterally traceable marker horizon, a 75 cm thick shell bed. Between the top of the OEB subsection and this shell bed, an interval of ~2 m is poorly exposed and therefore not sampled.

Near the base of the ACF in subsection OEB, an organic-rich interval of c. 1 m thick occurs containing thin ‘rusty’ layers (10 cm) with yellow and orange weathering colours (Fig. 5c). This interval is succeeded by predominantly silty shale, including two 20 cm thick sandstone layers. Shell fragments gradually decrease in abundance towards the top of the OEB subsection. The base of the OEF subsection is characterized by a 75 cm thick coquina bed, followed by silty shale with some shell fragments. The middle part, corresponding to the second organic-rich interval, is dominated by storm deposits, either consisting of centimetre- to decimetre-thick coquina beds, in some cases with HCS (Fig. 5e), or consisting of thin sandstone layers with HCS (Fig. 5d). In this interval yellow and orange (rusty) weathering levels occur (Fig. 5c), as well as laminated mudstones (Fig. 5f). Towards the top, the lithology becomes increasingly more silty and sandy. Storm beds are

<table><caption>**Table 1.** Overview of sections, samples and applied methods used in this study</caption><thead><tr><th rowspan="2">Sample set</th><th rowspan="2">Outcrop locality</th><th rowspan="2">Number of samples</th><th colspan="2">Analyses</th></tr><tr><th>δ<sup>13</sup>C<sub>org</sub></th><th>Palnology</th></tr></thead><tbody><tr><td>2020 field campaign</td><td>Cran aux Œufs</td><td>62</td><td>x</td><td>x</td></tr><tr><td rowspan="2">2019 duplicates Tribovillard</td><td>Cap de la Crèche</td><td>27</td><td>x</td><td></td></tr><tr><td>Cran du Noirda</td><td>28</td><td>x</td><td></td></tr><tr><td>2016 field campaign</td><td>Cap de la Crèche</td><td>13</td><td>x</td><td></td></tr></tbody></table>

absent in this part. The transition from the ACF to the Grès de la Crèche Fm appears to be very sharp from a distance; the colour changes abruptly from dark grey to yellow-brown (Fig. 5a). In terms of lithology, the transition is more subtle; sandy siltstones change into pure sandstones.

A second set of samples that was used in this study is a duplicate set derived from a geochemical study on the ACF that was published by Tribovillard et al. (2019). These samples are constrained to relatively small ORIs from two outcrop localities, Cap de la Crèche and Cran du Noirda (Fig. 3; Table 1). These samples were analysed for stable isotopes only (δ¹³Corg). Note, however, that supplementary clay mineral and geochemical data are available for these samples (Tribovillard et al. 2019) that allow comparison with the results of the present study.

A third sample set, collected in a fieldwork campaign by TNO in 2016, covers the transition from the ACF to the Grès de la Crèche Fm at outcrop locality Cap de la Crèche (Fig. 3; Table 1). This set was analysed for stable isotopes with the sole purpose of complementing the stable isotope curve from the ACF into the base of the Grès de la Crèche Fm (Fig. 6).

## 3.b. TOC and stable isotope analysis

In total, 130 samples were processed for stable isotope and TOC content analysis. Samples were powdered and subsequently dispatched to Iso-Analytical Ltd, Cheshire, UK, for TOC and stable isotope analysis. Carbonate was removed with hydrochloric acid, after which the samples were washed using distilled water. Subsequently, stable isotopes and TOC contents were measured using an elemental analyser isotope ratio mass spectrometer (EA-IRMS). Wheat flour (IA-R001, δ¹³C<sub>V-PDB</sub> = of -26.43 ‰) was used as a reference material. For quality control purposes, check samples of beet sugar (IA-R005, δ¹³C<sub>V-PDB</sub> = -26.03 ‰) and cane sugar (IA-R006, δ¹³C<sub>V-PDB</sub> = -11.64 ‰) were analysed during batch analysis of the samples.

## 3.c. Palynological processing and analysis

The palynological sample processing was carried out by Palynological Laboratory Services Ltd (PLS), in Anglesey, Wales, UK, according to standard palynological processing procedures. The standard processing routine includes treatment with hydrochloric acid to digest the carbonate and hydrofluoric acid in order to destroy the silicate minerals and release the acid-resistant OM. The organic residues were sieved using a 15 µm sieve. The organic residues were not treated with oxidizing agents. Norland cement was used as a permanent mounting medium for the sieved residue.
<figure><img src="image_8.png" /><figcaption>**Fig. 6.** (Colour online) Proposed correlation of the stable isotope (δ¹³Corg) and TOC trends from the Upper Jurassic of the Boulonnais area and the KCF from Kimmeridge Bay, southern England. The KCF record is modified after Morgans-Bell *et al.* (2001); it is a composite δ¹³Corg record through the KCF of the Swanworth Quarry 1 and Metherhills bore-holes. The Kimmeridge Bay record represents the more distal depositional setting. The stable isotope-based correlation allows for improved age constraints on the Upper Jurassic of the Boulonnais area. The numbers in the lithological column of the KCF refer to bed group numbers. The duration of the ammonite zones is derived from Hesselbo *et al.* (2020). Absolute ages must be made with reservation; the duration of the ammonite zones varies over the different editions of the geological timescale.</figcaption></figure>

For the palynological analysis, three categories were distinguished and counted separately: (1) palynofacies, (2) miscellaneous palynomorphs and (3) sporomorphs. The palynofacies and miscellaneous palynomorph analyses were carried out on all 62 available samples. A total of 43 samples was selected for the sporomorph analysis (see online Supplementary Material). Dinoflagellate cysts have not been counted on a species level on account of the rather poor state of preservation.

For the palynofacies analysis, between 200 and 350 OM particles were counted per sample. Five particle types were distinguished: (1) large particles of amorphous organic matter (AOM)

(AOM heavy), (2) small particles of AOM (AOM light), (3) palynomorphs, (4) wood, and (5) wood elongate. AOM heavy differs from AOM light in being larger (>100 µm) and darker. Microphotographs of particle types are displayed in Figure 13. The results are displayed in closed sum and frequency diagrams.

For the miscellaneous palynomorph analysis, between 100 and 150 palynomorph specimens were counted. Five different palynomorph types were distinguished: (1) acritarchs, (2) dinoflagellate cysts, (3) foraminiferal organic linings, (4) prasinophyte algae and (5) sporomorphs (Fig. 13). The results are displayed in closed sum and frequency diagrams.
Boulonnais

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

Kimmeridge Clay Formation

<figure><img src="image_10.png" /><figcaption>**Fig. 7.** δ¹³Corg against TOC plots for the ACF and KCF. (a) Lower and (b) upper ACF in the Cran aux Œufs section, and (c) the upper ACF in the Cap de la Crèche section. There is a clear positive correlation between δ¹³Corg and TOC in the upper ACF (b, c), which is absent in the lower ACF. In the KCF, δ¹³Corg against TOC plots for the (d) *Eudoxus*, (e) *Autissiodorensis*, (f) *Elegans* and (g) *Scitulus* ammonite zones display an increasing positive correlation with time. Plots for the age equivalent of the (h) lower and (i) upper ACF in the KCF based on δ¹³Corg correlations (Fig. 8) show, respectively, no correlation and a positive correlation.</figcaption></figure>

For the sporomorph analysis, between 100 and 200 pollen and spore specimens were counted. If possible, the specimens were taxonomically determined to species level, but afterwards all sporomorphs were combined into nine groups (Fig. 12): (1) *Callialasporites*, containing the species *Araucariacites australis*, *Callialasporites turbatus*, *C. dampieri* and *C. trilobatus*, (2)

*Cerebropollenites*, containing only *Cerebropollenites mesozoicus*, (3) *Classopollis*, containing only *Classopollis torosus*, (4) *Exesipollenites*, containing only *Exesipollenites tumulus*, (5) *Perinopollenites*, containing only *Perinopollenites elatoides*, (6) psilatrilete spores, containing only psilatrilete (smooth-walled) spores, (7) spores indifferent, containing *Cicatricosiporites* spp.,
<figure><img src="image_11.png" /><figcaption>**Fig. 8.** (Colour online) Local stable isotope correlation of outcrops from the Boulonnais region.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Fig. 9.** (Colour online) Compilation of the stable isotope and palynological results from the Cran aux Œufs section. Three closed sum diagrams display the palynological results: the miscellaneous palynomorph counts, the palynofacies counts and the sporomorph counts (excluding the bisaccate pollen and round indet. spore groups). The results of the stable isotope analyses are plotted in green (TOC) and red (δ¹³Corg) as an overlay on the palynofacies panel.</figcaption></figure>

Ischyosporites spp., Retitriletes spp., trilete granulate spores, trilete verrucate spores, Kraeuselisporites spp., Trachysporites spp., Taurocusporites spp., Gleicheniidites spp. and Contignisporites spp., (8) bisaccate pollen and (9) round indeterminable sporomorphs: these are likely either poorly preserved *Classopollis*, *Exesipollenites* or *Perinopollenites* specimens. The results of the sporomorph analysis are displayed in closed sum and frequency diagrams. The interpretations in terms of palaeoclimatic trends are established based on botanical affinities and ecological preferences of sporomorphs (Abbink, 1998 and references therein; Abbink et al. 2001, 2004; Schneider et al. 2018). *Classopollis*, for example, is derived from the extinct conifer family Cheirolepidiaceae that thrived under warm, seasonal arid conditions (Francis, 1984; Abbink et al. 2001). This pollen type is considered to be indicative of arid conditions, generally associated with coastal environments (Francis, 1984; Abbink et al. 2001). *Callialasporites* is associated with araucarian trees, which grew in forests not far from the sea and are considered indicative of warm climates and coastal habitats. Therefore, an assemblage with

dominant *Classopollis*, *Callialasporites* and *Araucariacites* is an indication of relatively warm and arid environments. *Perinopollenites* is considered an extinct representative of the Taxodiaceae and is associated with wet lowlands, probably indicating more humid environments (Abbink, 1998; Abbink et al. 2001). Spores are mainly produced in low-lying swampy deltaic-alluvial plain areas. Smooth trilete spores are from ferns representing a host of different parent plant affinities, but, as a rule of thumb, are associated with rivers and as such indicative of wet, lowland environments (Traverse, 1988). Bisaccate pollen are gymnosperm pollen from hinterland vegetation and are mostly transported into the marine realm by air.

Bisaccate pollen are excluded from the closed sum; this pollen type is generally overrepresented in marine environments due to their efficient wind and floating dispersal properties across long distances (Traverse, 1988) and are not very suitable to use as a proxy for palaeoclimate (Abbink et al. 2004). Round indeterminable sporomorphs are not included in the closed sum; the botаниcal affinity of this group is unclear.
<figure><img src="image_13.png" /><figcaption>**Fig. 10.** (Colour online) Sawtooth percentage diagrams of the palynofacies counts and the miscellaneous palynomorph counts. On the left, the lithological column of the Cran aux Œufs section is displayed together with the individual samples and with the units that are used for the description of the results. AOM – amorphous organic matter.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Fig. 11.** (Colour online) Closed sum and sawtooth diagrams showing the results of the sporomorph counts. The panel on the left is a closed sum diagram of seven sporomorph groups. The closed sum diagram does not include the bisaccate and round indet. spores groups; these are figured separately on the right as percentage of total sporomorphs. On the left-hand side the lithological column of the Cran aux Œufs section is displayed, together with the individual samples and with the units that are used for the description of the results.</figcaption></figure>

# 4. Results

## 4.a. TOC and stable isotope trends

The TOC contents in the Cran aux Œufs section range between 0.3 and 7% (Fig. 6; online Supplementary Material). The Cran aux Œufs section is subdivided into seven units for the description of the results. For completeness, a composite record for the ACF was constructed. It includes the entire ACF in the Cran aux Œufs section (solid line) and the transition from the upper ACF into the Grès de la Crèche Fm, from the Cap de la Crèche section (dashed line, Figs 6, 8). Two intervals with elevated TOC values occur, one in Unit II at the base of the lower ACF and one in Unit V in the upper ACF. The high TOC interval from Unit II displays a sharp increase in TOC, from 0.5 to 7%, followed by a more gradual decrease that continues into Unit III and stabilizes at a value of 1% in Unit IV. The high TOC interval of Unit V displays TOC values between 3 and 4% that remain stable throughout the 4 m thick Unit V. In the succeeding Unit VI, the TOC values remain elevated with an average of 2%. Towards the top of the ACF, in Unit VII, the TOC values drop to less than 1% and remain relatively stable across the boundary into the Grès de la Crèche Fm (Fig. 6, dashed line).

The measured δ¹³Corg values of the samples from the Cran aux Œufs section range between -26.25 ‰ and -23.94 ‰, with an average of -25.02 ‰ and a standard deviation of 0.46 ‰ (Fig. 6; online Supplementary Material). The degree of scatter is low, indicating that the results are reliable. In Unit II at the base of the ACF, the stable isotope trend starts off with a shift towards more positive values, followed by a sharp decline to the most negative values of the entire ACF record. After that, the stable isotope trend gradually increases and stabilizes to -25 ‰ in Unit III. In units V and VI of the upper ACF, more variation in the isotope signature is observed. Most conspicuous is the occurrence of several smaller-scale,

apparently rhythmic, trends. The stable isotope trend displays a pattern of shifts towards more positive values, followed by returns to more negative values. These trends appear to be rhythmic, with a vertical extent of ~2 m. In the uppermost part of the ACF the stable isotope trend stabilizes at a value of 24.7 ‰.

In the interval that covers units II and III, the high TOC interval in Unit II displays a positive correlation between TOC and δ¹³Corg (Fig. 6), but the remainder of the interval (Unit III) shows no correlation (R² of 0.0052; Fig. 7a). Conversely, in the interval that covers units IV–VII, a good fit exists between the TOC and the δ¹³Corg isotope trend: high TOC content correlates with positive shifts in δ¹³Corg values, with a higher R² of 0.35 and a slope of 0.25 (Fig. 7b). The positive correlation between TOC and the δ¹³Corg values is most evident in the 2019 sample set from the Cap de la Crèche section, with an R² of 0.92 and a comparable slope of 0.37 (Figs 7, 8).

## 4.b. Palynological results

The palynological assemblages from the Cran aux Œufs sections vary considerably in terms of richness and preservation. In the ORIs, the preservation of the palynomorphs is distinctly poorer than in the other intervals. For the richness, the same applies: the ORIs are less rich in palynomorphs compared to the low TOC intervals. The palynological results reveal both large- and small-scale trends. The results are discussed per unit in order to compare the palynology-based trends with the variations observed in lithology and in TOC (Fig. 9).

### 4.b.1. Palynofacies trends

The results of the palynofacies analyses are displayed in a closed sum diagram (Fig. 9) and in a sawtooth diagram (Fig. 10). The dominant groups are AOM and palynomorphs. In Unit II, 75 to 90% AOM is recorded. The very high AOM values drop to about
<figure><img src="image_15.png" /><figcaption>**Fig. 12.** (Colour online) Microphotographs of palynomorphs from the Argiles de Châtillon Formation. (a) *Classopollis* OEF-19. (b) *Classopollis* sample OEB-17. (c) *Perinopollenites* sample OEB-02. (d) *Perinopollenites* sample OEF-19. (e) *Exesipollenites* sample OEF-19. (f) *Exesipollenites* sample OEB-21. (g) *Callialasporites dampieri* sample OEB-22. (h) *Callialasporites dampieri* sample OEB-02. (i) *Cerebropollenites* sample OEB-02. (j) Trilete spore sample OEB-02. (k) *Gleicheniidites* sample OEF-19. (l) *Ischyosporites* sample OEF-19. (m) *Osmundacidites* sample OEB-21. (n) *Osmundacidites* sample OEB-21. (o) Dinocyst sample OEF-10. (p) Foraminifera remains sample OEB-02. Scale bars = 50 µm.</figcaption></figure>

20 % in the lower part of Unit III, after which the values increase again to 40 to 60 % in the middle part of Unit III. The entire upper ACF displays high abundances of AOM, but the highest levels are reached in Unit V, with an average value of 80 %. Lower percentages of AOM are observed in Unit VI, with an average of 45 %. Not surprisingly, the intervals rich in AOM correspond to the ORIs. In units V and VI of the upper ACF, the AOM distribution displays a rhythmic trend. These rhythms, or cyclic patterns, are typically 2 m long and occur superimposed on an overall decreasing trend towards the top of the section. Note that heavy AOM co-varies with AOM, but in much lower frequencies, and is most prominently represented in Unit V.

### 4.b.2. Miscellaneous palynomorph trends

The closed sum diagram of the miscellaneous palynomorphs (Fig. 9) shows an overall dominance of dinoflagellate cysts throughout the entire Cran aux Œufs section. Sporomorphs are the second most abundant palynomorph type. Foraminiferal organic linings are quite rare except for units II and V where they are common. The foraminiferal organic linings are mainly of the coiled *Globothalamea* type (Tyszka et al. 2021). Acritarchs are rare

and prasinophyte algae are in general very rare. Note however that prasinophyte algae are consistently present in Unit V and in the lower part of Unit VI. The dinoflagellate cysts increase from c. 40 % at the base of the ACF to c. 90 % at the top. Superimposed on the overall increasing trend, a higher order frequency pattern is displayed in which sporomorphs and dinoflagellate cysts alternate in abundance on a 50–70 cm scale. This feature is best expressed in the more densely sampled units II and III of the lower ACF.

### 4.b.3. Sporomorph trends

The results of the sporomorph analysis are displayed in closed sum and in sawblade diagrams (Figs 9, 11). Relevant sporomorph genera are displayed in Figure 12. The three most abundant groups are *Classopollis*, *Perinopollenites* and psilatrilete spores. The other groups, *Callialasporites*, *Cerebropollenites*, *Exesipollenites* and ‘spores indifferent’ are in general quite common. *Classopollis* and *Callialasporites* show a marked increase from the base to the top of the section. A distinct step towards higher percentages of *Classopollis* is recorded at the base of Unit IV. *Perinopollenites* and psilatrilete spores, on the other hand, show a concomitant
<figure><img src="image_16.png" /><figcaption>**Fig. 13.** (Colour online) Microphotographs of palynomorphs from the Argiles de Châtillon Formation. (a) Palynofacies from the lower OM-rich interval (sample OEB-03) showing dominant AOM; (b) same but with higher magnification. (c) Palynofacies from the upper OM-rich interval (sample OEF-10) showing dominant AOM. (d) Palynofacies from sample OEB-12, showing abundant wood particles. (e) Palynofacies from sample OEB-19, showing abundant sporomorphs. (f) Palynofacies from sample OEF-13 with a foraminiferal organic lining. (g) Palynofacies from sample OEF-03, showing heavy AOM. (h) Heavy AOM particle from sample OEF-03. (i) AOM particle from sample OEB-03; (j) the same AOM particle with incident UV light, showing an acritarch inside the particle. (k) Palynofacies from the lower OM-rich interval (sample OEB-03); (l) the same photograph but with incident UV light, showing a previously 'hidden' dinoflagellate cyst. (m) AOM particle from sample OEF-03; (n) the same particle with incident UV light. (o) An unidentifiable dinoflagellate cyst from sample OEF-13; (p) the same specimen with incident UV light. Scale bars = 50 µm.</figcaption></figure>

decrease from base to top. The lower ACF displays elevated abundances of spores compared to the upper ACF. The other groups, *Cerebropollenites*, *Exesipollenites* and indifferent spores, show variation but no distinct trends. The bisaccate pollen show an overall, gradual increasing trend in the lower ACF, and an increasing trend in the upper ACF that displays smaller-scale variation. The round indeterminable sporomorphs are most abundant in Unit II and in Unit V (Fig. 11).

### 4.b.4. *Trends and observations from the combined palynological categories*

Unit II is dominated by AOM, including the presence of heavy AOM. The upper and lower boundaries are sharp: the rise and the decline of AOM occurs almost instantaneously. Dinoflagellate cysts display the lowest recorded abundances in this interval. Sporomorphs, on the other hand, are abundant, and a peak in foraminiferal organic linings is present. Within the sporomorphs, a rapid increase in psilatrilete spores is recorded at the

base. The psilatrilete spores retain high percentages throughout the interval. *Classopollis* and *Callialasporites* display low relative abundances.

In Unit III, a gradual increase and a subsequent sharp decrease in palynomorphs occurs. In the middle of Unit III, the highest values of palynomorphs (75 %) are reached. AOM displays the opposite trend: a gradual decrease followed by a rapid increase. Low values of AOM (<15 %) are recorded in this interval. Within the palynomorphs, an overall increase in dinoflagellate cysts, concomitant with an overall decrease in sporomorphs, is observed. Within the sporomorphs, *Perinopollenites* and psilatrilete spores are equally abundant and together dominate the assemblages.

The base of Unit IV is characterized by a peak abundance of wood and palynomorphs, which in turn is immediately succeeded by a peak in AOM, including heavy AOM. Within the sporomorphs, a shift towards higher relative abundances of *Classopollis* is recorded at the base of this interval. *Callialasporites* and *Cerebropollenites* show a similar trend.
Indifferent spores reach their highest abundance at the base of Unit IV.

Unit V is characterized by the absolute dominance of AOM (up to 95 %). The highest relative abundances of heavy AOM and foraminiferal organic linings are recorded in this interval. Prasinophyte algae, which are mostly absent from the ACF record, are consistently present in this interval. Within the sporomorphs, *Classopollis*, *Cerebropollenites* and *Callialasporites* display parallel trends. Combined, these groups make up 25 to 50 % of the assemblages. The groups alternate in dominance with *Perinopollenites*, which accounts for 20 to 60 % of the assemblages.

Unit VI is characterized by a repeating pattern of an increase in AOM up to 70 %, followed by a decrease in AOM to 30 %. Wood remains relatively low. A general trend towards more palynomorphs is visible. Within the palynomorphs, the relative amounts of foraminiferal organic linings decrease, whereas the dinoflagellates increase. Within the sporomorphs, *Classopollis* displays an increase and becomes the most abundant group in this interval. *Perinopollenites* is in general quite low in abundance but shows an increase in the AOM-rich samples.

Unit VII is characterized by low abundances of AOM, values as low as 10 % being observed towards the top. Wood, on the contrary, shows a very rapid increase towards the top of the section, while palynomorphs also increase. Within the palynomorphs, dinoflagellate cysts are most abundant, at the cost of the sporomorphs. Foraminiferal organic linings are low in abundance. The sporomorph assemblage is relatively similar to that of the underlying interval, but an increase in *Callialasporites* and a decrease in *Classopollis* is noted towards the top.

# 5. Discussion

## 5.a. Isotope-based correlations to the KCF and chronostratigraphy

For the first time, a comprehensive stable isotope record for the Kimmeridgian–Tithonian transition of the Boulonnais area in northern France is established (Fig. 6). This stable isotope record can be correlated with the high-resolution stable isotope curve from the more distally marine setting of the KCF from southern England that was deposited in the same basin (Morgans-Bell et al. 2001; Jenkyns et al. 2021; Fig. 6). A correlation between the composite records of the Boulonnais and the KCF will allow for improved chronostratigraphic constraints in the ACF (Fig. 6), which is essential for the matching of climate and OM preservation trends through the Wessex Basin.

The negative trend towards the most negative peak in the upper part of Unit II is used as an anchor point for the lower part of the section and is correlated to the negative trend and peak between 357 and 354 m borehole depth in the KCF. The trend to more positive and stable values of around -25 ‰ in Unit III is roughly correlated to the interval between 355 and 340 m depth in the KCF. An important anchor point for the correlation of the ACF to the KCF is the shift towards more positive values in the ACF (Unit V). This shift is correlated with a similar shift in the uppermost part of the *Autissiodorensis* ammonite Zone of the KCF. In the units V and VI, the cyclic pattern of the stable isotope curve mimics the cyclic stable isotope pattern of the KCF curve. In the Boulonnais area, this pattern continues throughout the upper part of the ACF and into the lowermost part of the succeeding Grès de la Crèche Fm. In the KCF, the pattern continues until level 275 m in the *Scitulus* ammonite Zone.

Conversely to the KCF, the age constraints on the formations deposited in the Boulonnais area are relatively poor. This problem is due to the widespread presence of erosional sequence boundaries and the lateral variability of the succession (Wignall & Newton, 2001). The biostratigraphic work on the ACF from Geyssant et al. (1993) also indicates the scarcity of well-preserved ammonite remains in these formations.

The isotope correlation to the KCF suggests that units IV and V correlate to the upper part of the *Autissiodorensis* ammonite Zone (Irius Subzone), in the upper part of bed group 35 in the record of Morgans-Bell et al. (2001). The implied chronostratigraphy contradicts the available biostratigraphic data on the ACF (Geyssant et al. 1993; Herbin et al. 1995; Proust et al. 1995), which all indicate that the *Irius* Subzone is absent in the Boulonnais region. *Aulacostephanus volgensis* (Vischn.) (Upper *Autissiodorensis* ammonite Zone; Morgans-Bell et al. 2001) and *Pectinatites* (Arkellites) bleicheri (Loriol) (Elegans ammonite Zone; Morgans-Bell et al. 2001) were found, respectively, at the base (unit 15) and at the top (unit 18) of a 3 m thick interval in the Cran Barbier section (Geyssant et al. 1993). This latter section is located near the Cran aux Œufs section, but separated by several synsedimentary faults. As the stable isotope samples were taken at a relatively higher resolution, and in the absence of biostratigraphic constraints on the Cran aux Œufs section, it is not excluded that the equivalent of the lower part of bed group 35 from the KCF, the *Irius* Subzone, is present in the interval between 11 and 14 m in the Cran aux Œufs section. The presence of the *Irius* Subzone is likely limited to the region around the Cran aux Œufs section. As indicated by Wignall & Newton (2001), the most complete record of the ACF in terms of preserved parasequences lies near the Cran aux Œufs section. The stable isotope correlation panel between Cran aux Œufs, Cran du Noirda and Cap de la Crèche shows a similar result, with the thickest lower ACF sequence preserved at Cran aux Œufs. Apparently the depocentre shifted to this area during the deposition of the upper part of the lower ACF (Fig. 8). Furthermore, from the correlation to the KCF and the Geological Time Scale (Hesselbo et al. 2020), it is concluded that the Kimmeridgian–Tithonian boundary, which lies within the *Autissiodorensis* ammonite Zone, is located approximately at the base of the OEF subsection, on top of the three coquina beds.

## 5.b. Palynological response to variations in depositional environment

In shallow-marine, clastic-dominated Jurassic deposits, such as the ACF, palynology is generally considered to be an excellent tool for relative sea level interpretation. In particular, the marine–terrestrial ratio is often used as a proxy for relative sea level change (Partington et al. 1993; Duxbury et al. 1999; Fraser et al. 2003). The marine–terrestrial ratio is the ratio between marine-derived palynomorphs such as acritarchs, dinoflagellate cysts and foraminiferal organic linings, and land-derived palynomorphs, such as sporomorphs. The closed sum diagram (Fig. 9) shows a weak upward increase in dinoflagellate cysts and a concomitant weak upward decrease in sporomorphs (Figs 9, 10). Normally, this would be interpreted as a gradual deepening trend from a shallow-marine shelf setting to a more distal marine setting. However, this interpretation must be rejected because compelling evidence from other proxies points to a different sea level scenario. Based on lithology and OM content (Proust et al. 1993, 1995, 2001; Herbin et al. 1995; Tribovillard et al. 2001, 2019; Braaksma et al. 2006; Hatem et al. 2017), grain size distribution (Williams et al.
2001; Braaksma *et al.* 2006; Hatem *et al.* 2017), sedimentary struc-
tures such as HCS (Wignall & Newton, 2001; Hatem *et al.* 2017;
Tribovillard *et al.* 2019) and on the presence of an erosional uncon-
formity halfway up the ACF (Wignall & Newton, 2001), it is evi-
dent that the relative sea level scenario involves two transgressive-
regressive cycles, including two maximum flooding surfaces (MFS)
and one sequence boundary. The two MFS more or less correspond
to the organic-rich levels of units II and V (Fig. 6). Surprisingly,
these intervals display low abundances of dinoflagellate cysts.
Apparently, another mechanism than relative sea level was in play
here that determined the relative abundance of dinoflagellate cysts.

The absence of dinoflagellate cysts and the occurrence of abundant prasinophyte algae are indicative of stratification of the water column and bottom-water anoxia during oceanic anoxic events (Palliani et al. 2002; van de Schootbrugge et al. 2013; Houben et al. 2021). The decrease in dinoflagellate cysts and calcareous nannofossils (microplankton) and increase in green algae in the photic zone is interpreted as a reaction to palaeoenvironmental stress. In the Cran aux Œufs section, dinoflagellate cysts are always present and display a steady increase towards the top, whereas prasinophyte algae never reach substantial percentages. This distribution pattern is probably caused by a change from a stratified dysoxic or anoxic setting to a more oxygenated setting, which is in line with data from Tribovillard et al. (2019). Both sources of evidence indicate that fully stratified water-column conditions combined with bottom-water anoxia did not exist. Despite this, conditions were likely not permanently fully oxic either, as indicated by a general lack of bioturbation in the organic-rich laminated shales (Fig. 5f) and dysaerobic foraminifera assemblages, dominated by *Lenticulina* and *Ammobaculites* (Colpaert et al. 2021).

**5.c. Palaeoclimate reconstructions**

5.c.1. *Sporomorph-based climate trends*

Information about Late Jurassic climate in the ACF can be
obtained from the seven most-important sporomorph groups
and their associated climate indications, as explained in the meth-
ods section. There are trends visible on a larger (section) and
smaller (metre) scale. An increase in *Classopollis*, *Callialasporites*
and *Cerebropollenites* and a decrease in *Perinopollenites* is recorded
at the transition from the OEB (Kimmeridgian) to OEF
(Tithonian) subsection (Unit IV). This stratigraphical distribution
shows that an important change towards warmer and more arid
climate conditions occurred during the deposition of the ACF,
in agreement with the ‘Scitulus climate shift’ (Abbink, 1998)
and the clay mineral results of Hesselbo *et al.* (2009).

Smaller-scale climate trends are also recorded in the sporomorph dataset of the upper ACF. We see alternations of higher and lower abundances of the warmer and more arid against the cooler and more humid sporomorph groups, reflected in the sporomorph closed sum (red/left versus green/right; Fig. 9). The warm and arid sporomorphs fluctuate between ~50 % and ~30 %, but remain higher than in the OEB subsection, where they do not exceed ~20 %. In addition, the levels with dominance of the more humid sporomorphs correspond to the higher AOM and TOC contents in the storm bed dominated intervals, whereas the AOM-lean levels show higher proportions of arid sporomorphs.

5.c.2. Clay mineral data linked to climate variations

By integrating the clay mineral data from the Cap de la Crèche and
Cran du Noirda sections (Tribovillard *et al.* 2019; Fig. 8) with our

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

<figure><img src="image_18.png" /><figcaption>**Fig. 14.** (Colour online) Upper: Correlation of TOC and climate trends in clay mineral assemblages (kaolinite/illite and kaolinite/illite-smectite) of Tribovillard *et al.* (2019) for the Cap de la Crèche section. There is a clear large-scale positive correlation visible between TOC and humidity. Lower: Large-scale shift from humid climate below the sequence boundary to arid climate above the sequence boundary in Cran du Noirda. Note the higher values for the clay mineral ratios compared to those in Cap de la Crèche. The Scitulus climate shift is located within the sequence boundary, based on observations in the Cran aux Œufs record.</figcaption></figure>

stable isotope, TOC and palynological data for the same sample set,
we can place them in a context of climate and OM distribution.
Kaolinite commonly forms in soils under a tropical humid climate.
Therefore, the kaolinite/illite and kaolinite/illite-smectite mixed-
layer ratios are calculated (Fig. 14), where increased ratios are
indicative of more humid conditions on the continent (Thiry,
2000). In addition to the previously recorded large-scale variations
in the clay mineral assemblages (Deconinck et al. 1996; Hatem
et al. 2017; Tribovillard et al. 2019), we observe 2 m scale variations
in the distribution of the clay minerals as well.

In the Cran du Noirda record, there is a clear shift visible
towards lower ratios around the base of the upper ACF. The large
shift to lower values is indicative of a rapid shift towards a more
arid climate, which we correlate to the early Tithonian climate shift
that was recorded in the sporomorph dataset. For the upper ORI at
Cap de la Crèche, the data show a trend towards more humid con-
ditions with increasing TOC values. This pattern is very similar to
the increase in humid sporomorph assemblages with increasing
TOC in the equivalent interval in the Cran aux Œufs section.

The clay mineral record from Hatem *et al.* (2017) shows alter-
nations in smectite abundance on a similar 2 m scale in the upper
ACF, most evident in the most proximal section of Cap Gris Nez.
Our results shows that 2 m scale variability in the sporomorph dis-
tribution reflects humid–arid alternations. Climate fluctuations
might have similarly influenced the clay mineral signal via climate-
controlled runoff variations, with enhanced smectite in more
proximal sites, closer to the source area (Hatem *et al.* 2017).
The supply of clay minerals from the continent suggests runoff
as a primary mechanism controlling nutrient supply. In other
work, nutrient supply was also suggested as the link between clay
mineral variations and matching variations in OM content on an
orbital eccentricity timescale (Armstrong *et al.* 2016).

5.c.3. Milankovitch cycles in the ACF and KCF

In the Tithonian interval of the ACF (units V and VI), several
asymmetrical, yet repetitive, trends towards increasing values of
δ¹³C<sub>org</sub>, TOC and AOM are observed on a 2 m scale (Fig. 9). In
the lower-resolution sporomorph record, repetitive trends occur
that seem to line up with the palynofacies and stable isotope trends.
On the other hand, these alternations are clearly absent in units II
and III (Kimmeridgian). The strong correlation between the TOC,
palynofacies and sporomorph variability suggests a link between
OM accumulation and climate cyclicity for the upper ACF interval.
It is not uncommon to observe Milankovitch cyclicity in palyno-
logical data (e.g. Waterhouse, 1999; Bonis *et al.* 2010). Besides, sim-
ilar-scale cyclic variations have been described in the literature for
the ACF and KCF, although with varying interpretations on the
duration and underlying mechanism (Herbin *et al.* 1995;
Waterhouse, 1995, 1999; Weedon *et al.* 1999, 2004; Armstrong
*et al.* 2016). For example, Waterhouse (1999) observed palynofa-
cies cycles of 60–70 cm and 30–40 cm in the ACF. The ratio
between these cycles suggest obliquity and precession control. It
supports a roughly 100 kyr eccentricity control for the 2 m scale
alternations in the ACF.

The absolute age control on the duration of the *Elegans* Zone is
rather weak, but is in the order of 0.2–0.5 Myr (Weedon *et al.* 2004;
Ogg, 2005; Huang *et al.* 2010; Hesselbo *et al.* 2020). The presence of
three 2 m scale δ¹³Corg shifts and TOC fluctuations suggests a
roughly 100 kyr eccentricity control, in agreement with similar
palynofacies cycles (e.g. Waterhouse, 1995). Our results thus link
the variations in OM content on a 100 kyr scale to terrestrial veg-
etation changes that are characteristic of a changing Late Jurassic
climate.

5.d. Organic matter deposition

5.d.1. Spatiotemporal variations in δ¹³Corg and TOC content

The stable isotope values from the Cran aux Œufs section range
between -26.25 ‰ and -23.94 ‰, whereas the values for the same
time interval of the KCF (Morgans-Bell et al. 2001) reach more
positive values, ranging between -26.34 ‰ and -22.13 ‰. For
the TOC content, the same applies; the overall absolute values
for the KCF are considerably higher than the values of the ACF
from the Cran aux Œufs section. In addition, the TOC values of
the KCF display a larger variation. Interestingly, the difference
in absolute values is also observed within the smaller scale of
the Boulonnais area. The absolute TOC and δ¹³Corg values for
the Cap de la Crèche section are higher than the values of the
Cran aux Œufs section (Figs 7, 8). For the TOC values, this was
already noted by Herbin et al. (1995) and Tribovillard et al.
(2001). The Cap de la Crèche section is generally considered to
represent a slightly more distally marine depositional setting than
the Cran aux Œufs section, and the depositional setting of the
Kimmeridge Clay section is the most distally marine.

When looking at the variation of the TOC and δ¹³Corg values
within the ACF, a conspicuous difference is noted between the

Tithonian interval of the OEF subsection and the Kimmeridgian interval of the OEB subsection. A strong positive correlation between δ¹³Corg and TOC is recorded in units IV–VI of Tithonian age, but not in units II and III of Kimmeridgian age. A similar conclusion can be drawn from the record of the KCF, as indicated in Figure 7. Cross-plots between δ¹³Corg and TOC for the different ammonite stages ranging from *Eudoxus* to *Scitulus* show that the correlation between TOC and δ¹³Corg becomes more evident after the Kimmeridgian–Tithonian boundary. In fact, for the *Eudoxus* and *Autissiodorensis* ammonite zones the correlation is weak, whereas the R² for the *Elegans* and *Scitulus* zones are higher in both the Boulonnais sections and the KCF (Fig. 7d–g). A closer analysis of the δ¹³Corg and TOC records of the KCF (Fig. 6) shows that the positive correlation between TOC and δ¹³Corg starts at the base of bed group 35 (Upper *Autissiodorensis*, *Irius*), which correlates to the base of Unit IV and the base of the Tithonian. The base of the Tithonian thus marks an important shift in the geochemical signature recorded in the Wessex Basin (Fig. 7h, i).

5.d.2. Basin-wide controls on OM accumulation

Based on the similar pattern in the δ¹³C<sub>org</sub> records between the ACF and KCF, it is evident that the δ¹³C<sub>org</sub> record can be used to make high-resolution correlations on a regional scale, despite the variations in absolute values (Figs 6, 8). The stable isotope record from the Boulonnais area shows great similarities to that of the KCF, and it is therefore considered that both records were controlled by the same (basin-scale) mechanism. In fact, δ¹³C<sub>org</sub> correlations by Turner *et al.* (2019) showed that a water mass exchange with the North Sea Basin existed, and this correlation can potentially even be extended north of the Greenland–Norwegian seaway. The detailed correlation between the ACF and KCF (Fig. 6) shows that in the latest *Autissiodorensis* and *Elegans* ammonite zones (Tithonian), a mechanism was in play that led to a positive correlation between TOC and δ¹³C<sub>org</sub> and a cyclic pattern of the TOC curve, which was absent during the Kimmeridgian. Apparently, local factors that could have the influenced the δ¹³C<sub>org</sub> signal, such a sediment supply from the hinterland, climate and vegetation, were quite comparable through the basin, as well as on a larger regional scale. However, the differences between these formations in terms of lithology, depositional setting and absolute TOC are considerable (Proust *et al.* 1993; Herbin *et al.* 1995; Morgans-Bell *et al.* 2001; Wignall & Newton, 2001; Braaksma *et al.* 2006; Gallois *et al.* 2019). Therefore, a basin-wide controlling mechanism is required to explain the enhanced accumulation of OM in the basin.

Common boundary conditions for the formation of organic-
rich sediments are enhanced productivity of OM or enhanced pres-
ervation of OM. Most productivity estimates for the KCF and ACF
indicate that average productivity was in both cases not increased
(Tribovillard et al. 2001, 2005, 2019; Weedon et al. 2004; van
Dongen et al. 2006). One of the most intriguing differences is that
indications for photic zone anoxia and euxinic bottom-water con-
ditions are plentiful for the KCF (van Kaam-Peters et al. 1997,
1998; Sælen et al. 2000; van Dongen et al. 2006), but these indica-
tions are virtually absent for the ACF (Tribovillard et al. 2005,
2019; Hatem et al. 2017). Despite these differences, sedimentolog-
ical, palynological and palaeontological evidence shows that con-
ditions were likely not permanently fully oxic either. So, the
question arises how the enhanced OM accumulation could have
worked on a basin-wide scale, by assuming a mechanism that
can be applied also under predominant oxic conditions that are
generally not prone to OM enrichment.

5.d.3. Preservation of OM via sulfurization of carbohydrates

The positive correlation between TOC and δ¹³Corg has been
observed previously in the lower Tithonian of the KCF from the
Wessex Basin (Huc et al. 1992; Sinninghe Damsté et al. 1998;
van Kaam-Peters et al. 1998; Sælen et al. 2000; Morgans-Bell
et al. 2001; van Dongen et al. 2006; Fig. 7) and the Yorkshire basin
(Atar et al. 2019). Such correlation has been ascribed to the early
sulfurization of carbohydrates (Sinninghe Damsté et al. 1998; van
Kaam-Peters et al. 1998; van Dongen et al. 2006). In most environ-
ments, carbohydrates are rapidly remineralized and therefore con-
tribute to a very small proportion of the OM preserved in
sediments. However, carbohydrates are highly reactive towards
sulfides (Kok et al. 2000). Free sulfides are mostly formed in the
sediment by organoclastic sulfate-reduction and sulfate-driven
anaerobic oxidation of methane. Sulfate-reduction reactions can
take place both in the sediment and in the water column (Raven
et al. 2018). In the presence of free sulfides, carbohydrates
cross-link via the formation of sulfide bonds, or incorporate sulfur,
and thus become more resistant to remineralization. This abiotic
reaction of OM with sulfides, which also affects lipids, is called sul-
furization (Sinninghe Damsté & De Leeuw, 1990). The δ¹³Corg sig-
nature of carbohydrates is estimated to be 4–10 ‰ heavier than
that of lipids (Tyson, 1995; Sinninghe Damsté et al. 1998), which
explains why varying amounts of preserved carbohydrates may
result in δ¹³Corg shifts of up to 6 ‰. Since carbohydrates are also
rich in hydrogen, early sulfurization of carbohydrates contributes
to the preservation of high amounts of hydrogen-rich OM. This
leads to the general correlation of TOC, HI and δ¹³Corg values
as observed in the KCF.

Based on the positive correlation between TOC, organic sulfur and HI, Tribovillard *et al.* (2001, 2019) demonstrated that the OM of the ORIs from the upper ACF was sulfurized. The observed correlation between TOC and sulfurized OM is similar to what was described for the KCF (Huc *et al.* 1992; van Kaam-Peters *et al.* 1997, 1998). By analogy with the KCF, we propose that the correlation between AOM, TOC, HI and δ¹³Corg values observed in the Tithonian (units IV–VII) of the ACF results more specifically from the early sulfurization of carbohydrates.

It is important to realize that such large shifts in δ¹³Corg, up to 6 ‰, cannot be achieved by a change in the ratio between terrestrial and marine-derived kerogen alone (see discussion in Atar et al. 2019). In the Jurassic, marine type II OM had an average δ¹³Corg value of -28 ‰ (Tyson, 1995; van Dongen et al. 2006). Type III terrestrial OM had a δ¹³Corg signature ranging between -22 ‰ and -25 ‰ (Tyson, 1995; Pearce et al. 2005; van Dongen et al. 2006). The difference in δ¹³Corg of these kerogens implies that varying amounts can account for 3 ‰ at most (Tyson, 1995; van Dongen et al. 2006). The Jurassic type III OM had a more positive δ¹³Corg signature than that of type II OM. This is in contrast with its present-day signature, which is more negative than type II OM. Quaternary δ¹³Corg values for type III OM (Leng & Lewis, 2017) cannot be safely extrapolated to be used as a marine-terrestrial proxy for the Mesozoic. Furthermore, OM with a low HI and high oxygen index, plotting as type III OM and generally corresponding to terrestrial OM, can in fact alternatively correspond to altered type II marine OM (Hare et al. 2014). In addition, the terrestrial contribution to the total OM is weak in both the KCF (Tyson, 1989; Tribovillard et al. 1994; Van Kaam-Peters et al. 1998; Turner et al. 2019) and the ACF, where it is

mostly dominated by marine AOM and the contribution of terres-
trial OM to the palynofacies and palynomorphs rarely exceeds
35 % (Fig. 8). This implies that the δ¹³Corg signature is unlikely
to be influenced significantly by changes in the source of the
OM alone.

A remarkable outcome of this study is thus that during
Tithonian time, the apparently eccentricity-driven process of selec-
tive preservation of OM via the sulfurization of carbohydrates acts
simultaneously in the basin centre (southern England) as well as in
the more marginal areas of the basin (Boulonnais), as indicated by
the correlation based on stable isotopes (Fig. 6).

5.e. Depositional model: climate controls on OM deposition in the ACF

The present study shows that there are two ORIs in the ACF that
formed in contrasting settings, which agrees with the previous
work of Tribovillard et al. (2019). The differences between the
two ORIs are reflected in the palynological results and δ¹³Corg-
to-TOC correlations. The lower ORI is thinner, has sharply defined
boundaries, shows less variation in AOM content and reaches
higher TOC contents. The lowest abundance of dinoflagellate cysts
in the interval with highest TOC contents points towards moderate
stratification, suggesting that preservation played an important
role in the development of the lower ORI in the ACF, in agreement
with redox-sensitive trace elements of Tribovillard et al. (2005,
2019). The deposition occurred under stable, cooler and more
humid climate, during rising and high relative sea level.

The upper black shales have formed under contrasting condi-
tions. They were deposited after a climate shift to warmer and more
arid conditions in early Tithonian time, which is also documented
in other parts of NW Europe (Abbink et al. 1998; Hesselbo et al.
2009). There is a clear positive correlation between TOC and
δ¹³Corg. Here, there is no strong palynological evidence for strati-
fication of the water column, in agreement with the trace-element
analysis of Tribovillard et al. (2005, 2019). Despite this, the absence
of bioturbation in the paper shale intervals and the presence of
solely hypoxic foraminifera assemblages (Colpaert et al. 2021)
points in the direction of scarce benthic life and the (temporary)
absence of fully oxygenized conditions. Under these circumstan-
ces, sulfurization of OM, in particular carbohydrates, took place,
as evident from the positive correlation between TOC and
δ¹³Corg. This ORI is thicker but shows subordinate variations on
a 2 m scale. AOM and TOC content are likely controlled by the
100 kyr eccentricity cycle that controls alternations between more
humid and more arid phases, reflected in the sporomorph assem-
blage and the clay mineral ratios. An eccentricity-modulated shift
and expansion of the subtropical Hadley Cell was recently hypoth-
esized to have affected the Late Jurassic climate, resulting in alter-
nating humid and arid climate conditions (Armstrong et al. 2016;
Colombié et al. 2018; Atar et al. 2019). According to these model
simulations, monsoonal conditions, an intensified hydrological
cycle with more storminess, more continental weathering and
increased runoff occurred during the humid climate conditions.

Our results show that under the most humid phases of these
overall arid climate conditions, sulfurization of carbohydrates
was the dominant control on OM preservation. There are several
studies where climate has been described to affect the sulfurization
and preservation of OM (van Kaam-Peters et al. 1998; Riboulleau
et al. 2003; Atar et al. 2019). Iron supply is a key factor in sulfu-
rization (Hartgers et al. 1997; Tribovillard et al. 2015; Shawar
et al. 2018; Raven et al. 2021). As proposed for the KCF, the
<figure><img src="image_19.png" /><figcaption>**Fig. 15.** Depositional model for the ACF and KCF, explaining the (climate-controlled) mechanisms behind organic matter (OM) enrichment in the Kimmeridgian and Tithonian, and the process leading to OM sulfurization under a humid climate. The humid-arid alternations affected OM distribution on a basin-wide scale, whereas the ventilation caused by storminess resulted in oxygenated conditions in the proximal domain. PP – primary productivity; TOC – total organic carbon.</figcaption></figure>

availability of reactive iron likely depends on changing climate, probably driven by Milankovitch cycles, which can have a large impact on the amount of sulfurization via the supply of reactive iron from the continent (aeolian dust, runoff) (van Kaam-Peters et al. 1998). A similar link between increasing aridity, humid-arid alternations and OM preservation during Late Jurassic time was recorded by Riboulleau et al. (2003). Productivity was increased as a consequence of both aeolian and fluvial nutrient supply, combined with enhanced mixing during storms. The high primary productivity increased the oxygen demand, leading to anoxic conditions in the sediment, and less efficient reoxygenation during fair-weather periods.

Based on the similarity with these examples, we propose a possible model whereby the large-scale shift to a more arid and dry climate in Tithonian time triggered eccentricity-controlled humid-arid alternations, and the onset of sulfurization as the dominant control on OM preservation in these humid periods (Fig. 15). During the most humid phases, induced by the eccentricity-modulated shift and expansion of the subtropical Hadley Cell, monsoonal conditions and an intensified hydrological cycle increased continental weathering. The terrigenous material (smectite, Fe, nutrients) was efficiently transferred to the sea because of increased runoff (Armstrong et al. 2016; Colombié et al. 2018; Atar et al. 2019). This input of nutrients promoted primary productivity in the basin. As sulfurized OM was abundant in the humid intervals in the Tithonian, there should have been a (relative) reactive iron shortage (Tribovillard et al. 2002). The balance between aeolian and fluvial iron supply and consumption through iron sulfide precipitation can be subtle; small changes in primary productivity can drive large changes in OM preservation in environments near the redox threshold (Raven et al. 2018). Iron limitation can occur owing to high production of reduced sulfur in the sediment, induced by high inputs of OM and high sulfate-reduction rates (Tribovillard et al. 1994; Riboulleau et al. 2003). If OM supply increases and fosters sulfate-reduction reactions, the created sulfide can react with iron, while an excess of free sulfides can react with OM. Thus, a higher flux of OM leads to oxygen consumption in the sediment, favouring sulfurization of carbohydrates. Nevertheless, the increased storminess leads to the recurrent mixing of the water column, so that the upper sediment is never fully anoxic, at least in the marginal part of the basin. In the more distal

setting, the basin is deeper, so that even if storms are present, the water column might not be entirely mixed. Anoxia of the sediment and possibly water column can develop better, leading to even better preservation of OM.

Independent of how exactly climate drives the sulfurization of OM, different mechanisms characterize OM accumulation in the Kimmeridgian versus the Tithonian. The eccentricity-controlled humid-arid alternations were initiated during Tithonian time, simultaneously with and potentially linked to the aridification of the NW European continent. This regional, climate-driven mechanism is recorded in both the distal facies of the KCF as well as the proximal facies of the ACF of Tithonian age.

## 6. Conclusions

Palynological, TOC and $\delta^{13}C_{org}$ analyses were combined to study the OM distribution in relation to palaeoclimate for the Kimmeridgian–Tithonian interval of the ACF from the Boulonnais area. The ACF contains two intervals with elevated TOC, one in the uppermost Kimmeridgian, corresponding to the lower part of the *Autissiodorensis* ammonite Zone, and one in the lowermost Tithonian, corresponding to the upper *Autissiodorensis* (Irius) and *Elegans* ammonite zones. Because the ACF was deposited in a shallow-marine setting, the acquired quantitative pollen and spore record is considered to reflect actual changes in vegetation patterns that can be used as a proxy for palaeoclimate. The climate-driven changes in terrestrial vegetation are identified in the shallow-marine domain and are correlated to a more distal marine setting. The pollen and spore record of the ACF reveals a change towards warmer and more arid climate conditions around the Kimmeridgian–Tithonian boundary. Superposed on the large-scale climate trend, small-scale cycles observed in the palynological record of the upper ACF were interpreted to reflect eccentricity-driven increases in humidity and match alternations in OM content, TOC and $\delta^{13}C_{org}$. Parallel to the palynological record, a high-resolution $\delta^{13}C_{org}$ record was established that allows direct comparison and correlation of the proximal mudstones from the ACF to the distal mudstones of the KCF. The detailed correlation leads to the conclusion that the underlying mechanism that dictated the trends in TOC and $\delta^{13}C_{org}$ of the KCF is also applicable to the shallower marine ACF record, albeit in a more moderate fashion. The cyclic increases in TOC and $\delta^{13}C_{org}$ reflect sulfurization of carbohydrates, which we related to the eccentricity-driven humid-arid alternations. The ORI from the uppermost Kimmeridgian *Autissiodorensis* Zone does not exhibit this relationship: the TOC values and the $\delta^{13}C_{org}$ do not show a correlation and small-scale climate trends are absent in the palynological record. Here, OM accumulation is thus solely related to stratification of the water column. In conclusion, an important change occurred in the NW European climate across the Kimmeridgian–Tithonian boundary, towards warmer and more arid conditions with superimposed eccentricity-controlled changes in humidity. The latter appears to be a crucial driver for sulfurization of OM, which controlled OM accumulation in the Tithonian period on a basin-wide scale.

**Supplementary material.** To view supplementary material for this article, please visit https://doi.org/10.1017/S0016756822001121

**Acknowledgements.** This work was financially supported by TNO Geological Survey of the Netherlands. We thank Nico Janssen, Sander Houben and Bart van der Kwaak for scientific discussions and/or assistance during field work.
We thank Andrew Aplin and one anonymous reviewer for their comments that improved the manuscript.

**Conflict of interest.** None.

## References

Abbink O (1998) Palynological investigations in the Jurassic of the North Sea region. *LPP Contributions Series* 8, 192 pp.

Abbink OA, Mijnlieff HF, Munsterman DK and Verreussel RMCH (2006) New stratigraphic insights in the Late Jurassic of the Southern Central North Sea Graben and Terschelling Basin (Dutch Offshore) and related exploration potential. *Netherlands Journal of Geosciences* 85, 221–38.

Abbink O, Targarona J, Brinkhuis H and Visscher H (2001) Late Jurassic to earliest Cretaceous palaeoclimatic evolution of the southern North Sea. *Global and Planetary Change* 30, 231–56.

Abbink O, Van Konijnenburg-Van Cittert JHA and Visscher H (2004) A sporomorph ecogroup model for the Northwest European Jurassic-Lower Cretaceous: concepts and framework. *Netherlands Journal of Geosciences* 83, 17–31.

Agustsdottir AM, Barron EJ, Bice KL, Colarusso LA, Cookman JL, Cosgrove BA, De Lurio JL, Dutton JF, Frakes BJ, Frakes LA, Moy CJ, Olszewski TD, Pancost RD, Poulsen CJ, Ruffner CM, Sheldon DG and White TS (1999) Storm activity in ancient climates 2. An analysis using climate simulations and sedimentary structures. *Journal of Geophysical Research: Atmospheres* 104, 27295–320.

Angus L, Hampson GJ, Palci F and Fraser AJ (2020) Characteristics and context of high-energy, tidally modulated, barred shoreface deposits: Kimmeridgian–Tithonian sandstones, Weald Basin, southern UK and northern France. *Journal of Sedimentary Research* 90, 313–35.

Armstrong HA, Wagner T, Herringshaw LG, Farnsworth AJ, Lunt DJ, Harland M, Imber J, Loptson C and Atar EF (2016) Hadley circulation and precipitation changes controlling black shale deposition in the Late Jurassic Boreal Seaway. *Paleoceanography* 31, 1041–53.

Atar E, März C, Aplin AC, Dellwig O, Herringshaw LG, Lamoureux-Var V, Leng MJ, Schnetger B and Wagner T (2019) Dynamic climate-driven controls on the deposition of the Kimmeridge Clay Formation in the Cleveland Basin, Yorkshire, UK. *Climate of the Past* 15, 1581–601.

Baniak GM, Gingras MK, Burns BA and George Pemberton S (2014) An example of a highly bioturbated, storm-influenced shoreface deposit: Upper Jurassic Ula Formation, Norwegian North Sea. *Sedimentology* 61, 1261–85.

Bonis NR, Ruhl M and Kürschner WM (2010) Milankovitch-scale palynological turnover across the Triassic–Jurassic transition at St. Audrie's Bay, SW UK. *Journal of the Geological Society, London* 167, 877–88.

Braaksma H, Drijkoningen GG, Filippidou N, Kenter JAM and Proust JN (2006) The origin and nature of seismic reflections of sharp-based shoreface deposits (upper Jurassic Siliciclastics, northern France). *Geophysical Prospecting* 54, 211–36.

Colombié C, Carcel D, Lécuyer C, Ruffel A and Schnyder J (2018) Temperature and cyclone frequency in Kimmeridge Greenhouse period (Late Jurassic). *Global and Planetary Change* 170, 126–45.

Colpaert CM, Nikitenko B, Danelian T, Elena F and Igor K (2021) Upper Kimmeridgian to Lower Tithonian foraminifers and ostracods of the Boulonnais region, Northern France. Stratigraphic and palaeoenvironmental implications and inter-regional correlations with the Siberian arctic. *Marine Micropaleontology* 164, 101976. doi: 10.1016/j.marmicro.2021.101976.

Cornford C (2018) Petroleum systems of the South Viking Graben. In *Rift-Related Coarse-Grained Submarine Fan Reservoirs; the Brae Play, South Viking Graben, North Sea* (eds CC Turner and BT Cronin). American Association of Petroleum Geologists Memoir vol. 115.

Coward MP, Dewey JF, Hempton M and Holroyd J (2003) Tectonic evolution. In *The Millennium Atlas: Petroleum Geology of the Central and Northern North Sea* (eds D Evans, C Graham, A Armour and P Bathurst), pp. 17–33. London: Geological Society of London.

Deconinck J-F and Baudin F (2008) Les dépôts du Kimmeridgien et du Tithonien du nord-ouest du bassin de Paris (Haute-Normandie, Boulonnais). *Annales de la Société Géologique du Nord* 15, 77–90.

Deconinck J-F, Chamley H, Debrabant P and Colebeaux J-P (1983) Le Boulonnais au Jurassique supérieur: donnees de la mineralogie des argiles et de la geochemie. *Annales de la Societe Geologie Nord* 102, 145–52.

Deconinck J-F, Geyssant JR, Proust J-N and Vidier JP (1996) Sédimentologie et biostratigraphie des dépôts Kimméridgiens et Tithoniens du Boulonnais. *Annales de la Société Géologique du Nord* 4, 157–70.

Donnadieu Y, Dromart G, Goddéris Y, Pucéat E, Brigaud B, Dera G, Dumas C and Olivier N (2011) A mechanism for brief glacial episodes in the Mesozoic greenhouse. *Paleoceanography* 26, PA3212. doi: 10.1029/2010PA002100

Duxbury SDSR, Kadolsky D and Johansen S (1999) Sequence stratigraphic subdivision of the Humber Group in the Outer Moray Firth area (UKCS, North Sea). In *Biostratigraphy in Production and Development Geology* (eds RW Jones and MD Simmons), pp. 23–54. Geological Society of London, Special Publication no. 152.

Francis JE (1984) The seasonal environment of the Purbeck (Upper Jurassic) fossil forests. *Palaeogeography, Palaeoclimatology, Palaeoecology* 48, 285–307.

Fraser SI, Robinson AM, Johnson HD, Underhill JR, Kadolsky DGA, Connell R, Johannesson P and Ravnás R (2003) Upper Jurassic. In *The Millennium Atlas: Petroleum Geology of the Central and Northern North Sea* (eds D Evans, C Graham, A Armour and P Bathurst), pp. 157–89. London: Geological Society of London.

Fürsich FT (1981) Salinity-controlled benthic associations from the Upper Jurassic of Portugal. *Lethaia* 14, 203–23.

Fürsich FT and Oschmann W (1986) Storm shell beds of *Nanogyra virgula* in the Upper Jurassic of France. *Neues Jahrbuch Geologische Paläontologische Abhandlungen* 172, 141–61.

Gallois R, Vadet A and Etches S (2019) Correlation of the Kimmeridgian–Tithonian (Jurassic) boundary beds exposed in the Boulonnais, France with those at Kimmeridge, Dorset, UK. *Proceedings of the Geologists’ Association* 130, 187–95.

Geyssant JR, Vidier JP, Herbin JP, Proust JN and Deconinck JF (1993) Biostratigraphie et paléoenvironnement des couches de passage Kimméridgien/Tithonien du Boulonnais (Pas-de-Calais): nouvelles données paléontologiques (ammonites), organisation séquentielle et contenu en matière organique. *Géologie de la France* 4, 11–24.

Hallam A (1975) *Jurassic Environments*. Cambridge: Cambridge University Press, 269 pp.

Hallam A (1984) Continental humid and arid zones during the Jurassic and Cretaceous. *Palaeogeography, Palaeoclimatology, Palaeoecology* 47, 195–223.

Hallam A (1985) A review of Mesozoic climates. *Journal of the Geological Society, London* 142, 433–45.

Hallam A (1993) Jurassic climates as inferred from the sedimentary and fossil record. *Philosophical Transactions of the Royal Society of London, Series B: Biological Sciences* 341, 287–96.

Hare AA, Kuzyk ZZA, Macdonald RW, Sanei H, Barber D, Stern GA and Wang F (2014) Characterization of sedimentary organic matter in recent marine sediments from Hudson Bay, Canada, by Rock-Eval pyrolysis. *Organic Geochemistry* 68, 52–60.

Hartgers WA, Lopez JF, Sinninghe Damsté JS, Reiss C, Maxwell J and Grimalt JO (1997) Sulfur-binding in recent environments: II. Speciation of sulfur and iron and implications for the occurrence of organo-sulfur compounds. *Geochimica et Cosmochimica Acta* 61, 4769–88.

Hatem E, Tribovillard N, Averbuch O, Bout-Rouazeilles V, Trentesaux A, Deconinck JF, Baudin F and Adatte T (2017) Small-scaled lateral variations of an organic-rich formation in a ramp-type depositional environment (the Late Jurassic of the Boulonnais, France): impact of the clastic supply. *Bulletin de la Société géologique de France* 188, 31–47.

Hatem E, Tribovillard N, Averbuch O, Sansjofre P, Adatte T and Guillot F (2016) Early diagenetic formation of carbonates in a clastic dominated ramp environment impacted by synsedimentary faulting-induced fluid seepage – evidence from the Late Jurassic Boulonnais Basin (N France). *Marine and Petroleum Geology* 72C, 12–29.

Hatem E, Tribovillard N, Averbuch O, Vidier D, Sansjofre P, Birgel D and Guillot F (2014) Oyster patch reefs as indicators of fossil hydrocarbon seeps
induced by synsedimentary faults. *Marine and Petroleum Geology* 55, 176-85.

Herbin JP, Fernandez-Martinez JL, Geyssant JR, Albani AE, Deconinck JF, Proust JN, Colbeaux JP and Vidier JP (1995) Sequence stratigraphy of source rocks applied to the study of the Kimmeridgian/Tithonian in the north-west European shelf (Dorset/UK, Yorkshire/UK and Boulonnais/ France). *Marine and Petroleum Geology* 12, 177-94.

Herbin JP, Muller C, Geyssant J, Melieres F and Penn I (1991) Heterogeneity of organic matter distribution in relation to a transgressive systems tract: Kimmeridge Clay (Jurassic), England, *American Association of Petroleum Geologists Bulletin* 75, 593-4.

Hesselbo SP, Deconinck JF, Huggett JM and Morgans-Bell HS (2009) Late Jurassic palaeoclimatic change from clay mineralogy and gamma-ray spectrometry of the Kimmeridge Clay, Dorset, UK. *Journal of the Geological Society, London* 166, 1123-33.

Hesselbo SP, Ogg JG, Ruhl M, Hinnov LA and Huang CJ (2020) The Jurassic Period. In *Geologic Time Scale 2020* (eds FM Gradstein, JG Ogg, MD Schmitz and GM Ogg), pp. 955–1021. Amsterdam: Elsevier.

Houben AJ, Goldberg T and Slomp CP (2021) Biogeochemical evolution and organic carbon deposition on the Northwestern European Shelf during the Toarcian Ocean Anoxic Event. *Palaeogeography, Palaeoclimatology, Palaeoecology* 565, 110191. doi: 10.1016/j.palaeo.2020.110191.

Huang C, Hesselbo SP and Hinnov L (2010) Astrochronology of the late Jurassic Kimmeridge Clay (Dorset, England) and implications for Earth system processes. *Earth and Planetary Science Letters* 289, 242–55.

Huc AY, Lallier-Vergès E, Bertrand P, Carpentier B and Hollander DJ (1992) Organic matter response to change of depositional environment in Kimmeridgian Shales, Dorset, U.K. In *Organic Matter: Productivity, Accumulation, and Preservation in Recent and Ancient Sediments* (eds JK Whelan and JW Farrington), pp. 469–86. New York: Columbia University Press.

Jenkyns HC, Coe AL, Gallois RW, Hesselbo SP, Marshall JEA, Morgans-Bell HS, Sellwood BW, Tyson RV, Weedon GP and Williams CJ (2021) Wireline Logging and Core Analysis Datasets Collected from Swanworth Quarry and Metherhills Boreholes during the NERC Funded Rapid Global Geological Events (RGGE) Project. NERC EDS National Geoscience Data Centre. (Dataset). doi: 10.5285/a8e86b5e-f341-419e-be98-0463efe3201f.

Kok MD, Schouten S and Damsté JSS (2000) Formation of insoluble, nonhydrolyzable, sulfur-rich macromolecules via incorporation of inorganic sulfur species into algal carbohydrates. *Geochimica et Cosmochimica Acta* 64, 2689-99.

Korte C and Hesselbo SP (2011) Shallow marine carbon and oxygen isotope and elemental records indicate icehouse-greenhouse cycles during the Early Jurassic. *Paleoceanography* 26, PA4219. doi: 10.1029/2011PA002160.

Leng MJ and Lewis JP (2017) Bulk C/N ratios and carbon isotopes in estuarine environments. In *Applications of Paleoenvironmental Techniques in Estuarine Studies* (eds K Weckström, K Saunders, P Gell and G Skilbeck), pp. 213–37. Developments in Palaeoenvironmental Research 20. Dordrecht, Netherlands: Springer.

Mansy JL, Manby GM, Averbuch O, Everaerts M, Bergerat F, van Vliet-Lanoë B, Lamarche J and Vandycke S (2003) Dynamics and inversion of the Mesozoic Basin of the Weald-Boulonnais area: role of basement reaction. *Tectonophysics* 373, 161–79.

Morgans-Bell HS, Coe AL, Hesselbo SP, Jenkyns HC, Weedon GP, Marshall JE, Tyson RV and Williams CJ (2001) Integrated stratigraphy of the Kimmeridge Clay Formation (Upper Jurassic) based on exposures and boreholes in south Dorset, UK. *Geological Magazine* 138, 511–39.

Ogg J (2005) The Jurassic Period. In *A Geologic Time Scale 2004* (eds F Gradstein, J Ogg and A Smith), pp. 307–43. Cambridge: Cambridge University Press.

Oschmann W (1988) Kimmeridge Clay sedimentation — a new cyclic model. *Palaeogeography, Palaeoclimatology, Palaeoecology* 65, 217–51.

Palliani RB, Mattioli E and Riding JB (2002) The response of marine phytoplankton and sedimentary organic matter to the early Toarcian (Lower Jurassic) Oceanic Anoxic Event in northern England. *Marine Micropaleontology* 46, 223–45.

Partington MA, Copestake P, Mitchener BA and Underhill JR (1993) Biostratigraphic calibration of genetic stratigraphic sequences in the

Jurassic-lowermost Cretaceous (Hettangian to Ryazanian) of the North Sea and adjacent areas. In *Petroleum Geology of Northwest Europe: Proceedings of the 4th Conference* (ed. JR Parker), pp. 371–86. Geological Society of London, Petroleum Geology Conference Series vol. 4.

Pearce CR, Hesselbo SP and Coe AL (2005) The mid-Oxfordian (Late Jurassic) positive carbon-isotope excursion recognised from fossil wood in the British Isles. *Palaeogeography, Palaeoclimatology, Palaeoecology* 221, 343–57.

Proust JN, Deconinck JF, Geyssant JR, Herbin JP and Vidier JP (1993) Nouvelles données sédimentologiques dans le Kimméridgien et le Tithonien du Boulonnais (France). *Comptes Rendus de l'Académie des Sciences Paris* 316, 363-9.

Proust JN, Deconinck JF, Geyssant JR, Herbin JP and Vidier JP (1995) Sequence analytical approach to the Upper Kimmeridgian/Lower Tithonian storm-dominated ramp deposits of the Boulonnais. A landward time equivalent to offshore marine source rocks (Northern France). *Geologische Rundschau* 84, 225-71.

Proust J, Mahieux G and Tessier B (2001) Field and seismic images of sharp-based shoreface deposits: implications for sequence stratigraphic analysis. *Journal of Sedimentary Research* 71, 944-57.

Raven MR, Fike DA, Gomes ML, Webb SM, Bradley AS and McClelland HLO (2018) Organic carbon burial during OAE2 driven by changes in the locus of organic matter sulfurization. *Nature Communications* 9, 1-9.

Raven MR, Keil RG and Webb SM (2021) Rapid, concurrent formation of organic sulfur and iron sulfides during experimental sulfurization of sinking marine particles. *Global Biogeochemical Cycles* 35, e2021GB007062. doi: 10.1029/2021GB007062.

Riboulleau A, Baudin F, Deconinck JF, Derenne S, Largeau C and Tribovillard N (2003) Depositional conditions and organic matter preservation pathways in an epicontinental environment: the Upper Jurassic Kashpir Oil Shales (Volga Basin, Russia). *Palaeogeography, Palaeoclimatology, Palaeoecology* 197, 171-97.

Sælen G, Tyson RV, Telnæs N and Talbot MR (2000) Contrasting watermass conditions during deposition of the Whitby Mudstone (Lower Jurassic) and Kimmeridge Clay (Upper Jurassic) formations, UK. *Palaeogeography, Palaeoclimatology, Palaeoecology* 163, 163-96.

Schlurf M (2003) Palaeoecologic significance of Late Jurassic trace fossils from the Boulonnais, N France. *Acta Geologica Polonica* 53, 123-42.

Schneider AC, Heimhofer U, Heunisch C and Mutterlose J (2018) From arid to humid—the Jurassic–Cretaceous boundary interval in northern Germany. *Review of Palaeobotany and Palynology* 55, 57–69.

Schnyder J, Ruffell A, Deconinck JF and Baudin F (2006) Conjunctive use of spectral gamma-ray logs and clay mineralogy in defining late Jurassic–early Cretaceous palaeoclimate change (Dorset, UK). *Palaeogeography, Palaeoclimatology, Palaeoecology* 229, 303–20.

Sellwood BW and Valdes PJ (2008) Jurassic climates. *Proceedings of the Geologists' Association* 119, 5-17.

Shawar L, Haley I, Said-Ahmad W, Feinstein S, Boyko V, Kamyshny A and Amrani A (2018) Dynamics of pyrite formation and organic matter sulfurization in organic-rich carbonate sediments. *Geochimica et Cosmochimica Acta* 241, 219-39.

Sinninghe Damsté JS and De Leeuw JW (1990) Analysis, structure and geochemical significance of organically-bound sulphur in the geosphere: state of the art and future research. *Organic Geochemistry* 16, 1077-101.

Sinninghe Damsté JS, Kok MD, Köster J and Schouten S (1998) Sulfurized carbohydrates: an important sedimentary sink for organic carbon? *Earth and Planetary Science Letters* 164, 7-13.

Thiry M (2000) Palaeoclimatic interpretation of clay minerals in marine deposits: an outlook from the continental origin. *Earth-Science Reviews* 49, 201-21.

Traverse A (1988) *Paleopalynology*. Boston: Unwin Hyman, 600 pp.

Tribovillard N, Averbuch O, Bialkowski A and Deconinck JF (2002) Early diagenesis of marine organic-matter and magnetic properties of sedimentary rocks: the role of iron limitation and organic-matter source organisms. *Bulletin de la Société géologique de France* 173, 295-306.

Tribovillard N, Bialkowski A, Tyson RV, Lallier-Verges E and Deconinck JF (2001) Organic facies variation in the Late Kimmeridgian of the Boulonnais area (northernmost France). *Marine and Petroleum Geology* 18, 371-89.

Tribovillard NP, Desprairies A, Lallier-Verges E, Bertrand P, Moureau N, Ramdani A and Ramanampisoa L (1994) Geochemical study of organic-
rich cycles from the Kimmeridge Clay Formation of Yorkshire (UK): produc-
tivity versus anoxia. *Palaeogeography, Palaeoclimatology, Palaeoecology* **108**,
165–81.

Tribovillard N, Hatem E, Averbuch O, Barbecot F, Bout-Rouazeilles V and Trentesaux A (2015) Iron availability as a dominant control on the primary composition and diagenetic overprint of organic-matter-rich rocks. *Chemical Geology* **401**, 67–82.

Tribovillard N, Koched H, Baudin F, Adatte T, Delattre M, Abraham R and Ferry JN (2019) Storm-induced concentration of sulfurized, marine-origin, organic matter as a possible mechanism in the formation of petroleum source-rock. *Marine and Petroleum Geology* **109**, 808–18.

Tribovillard N, Petit A, Quijada M, Riboulleau A, Sansjofre P, Thomazo C, Huguet A, Birgel D and Averbuch O (2018) A genetic link between synsedimentary tectonics-expelled fluids, microbial sulfate reduction and cone-in-cone structures. *Marine and Petroleum Geology*, **93**, 437–50.

Tribovillard N, Ramdani A and Trentesaux A (2005) Controls on organic accumulation in Upper Jurassic shales of northwestern Europe as inferred from trace-metal geochemistry. In *The Deposition of Organic-Carbon-Rich Sediments: Models, Mechanisms, and Consequences* (ed. N Harris), pp. 145–64. SEPM Special Publication vol. 82.

Turner HE, Batenburg SJ, Gradstein FM and Gale AS (2019) The Kimmeridge Clay Formation (Upper Jurassic–Lower Cretaceous) of the Norwegian Continental shelf and Dorset, UK: a chemostratigraphic correlation. *Newsletters on Stratigraphy* **52**, 1–32.

Tyson RV (1989) Late Jurassic palynofacies trends, Piper and Kimmeridge Clay formations, UK onshore and northern North Sea. In *Northwest European Micropalaeontology and Palynology* (eds DJ Batten and MC Keen), pp. 135–72. Chichester: Ellis Horwood.

Tyson RV (1995) *Sedimentary Organic Matter. Organic Facies and Palynofacies*. London: Chapman Hall.

Tyszka J, Godos K, Goleń J and Radmacher W (2021) Foraminiferal organic linings: functional and phylogenetic challenges. *Earth-Science Reviews* **220**, 103726. doi: 10.1016/j.earscirev.2021.103726.

Underhill JR and Stoneley R (1998) Introduction to the development, evolution and petroleum geology of the Wessex Basin. In *Development, Evolution and Petroleum Geology of the Wessex Basin* (ed. JR Underhill), pp. 1–18. Geological Society of London, Special Publication no. 133.

Valdes PJ and Sellwood BW (1992) A palaeoclimate model for the Kimmeridgian. *Palaeogeography, Palaeoclimatology, Palaeoecology*, **95**, 47–72.

van de Schootbrugge B, Bachan A, Suan G, Richoz S and Payne JL (2013) Microbes, mud and methane: cause and consequence of recurrent Early Jurassic anoxia following the end-Triassic mass extinction. *Palaeontology* **56**, 685–709.

van Dongen BE, Schouten S and Sinninghe Damsté JS (2006) Preservation of carbohydrates through sulfurization in a Jurassic euxinic shelf sea:

examination of the Blackstone Band TOC cycle in the Kimmeridge Clay
Formation, UK. *Organic Geochemistry* **37**, 1052–73.

van Kaam-Peters HM, Schouten S, de Leeuw JW and Sinninghe Damsté JS (1997) A molecular and carbon isotope biogeochemical study of biomarkers and kerogen pyrolysates of the Kimmeridge Clay Facies: palaeoenvironmental implications. *Organic Geochemistry* **27**, 399–422.

van Kaam-Peters HM, Schouten S, Köster J and Sinninghe Damsté JS (1998) Controls on the molecular and carbon isotopic composition of organic matter deposited in a Kimmeridgian euxinic shelf sea: evidence for preservation of carbohydrates through sulfurisation. *Geochimica et Cosmochimica Acta* **62**, 3259–83.

Verreussel RMCH, Bouroullec R, Munsterman DK, Dybkjær K, Geel CR, Houben AJP, Johannessen PN and Kerstholt-Boegehold SJ (2018) Stepwise basin evolution of the Middle Jurassic-Early Cretaceous rift phase in the Central Graben area of Denmark, Germany and The Netherlands. In *Mesozoic Resource Potential in the Southern Permian Basin* (eds B Kilhams, PA Kukla, S Mazur, T McKie, HF Mijnlieff and K van Ojik), pp. 305-40. Geological Society of London, Special Publication no. 469.

Waterhouse HK (1995) High-resolution palynofacies investigation of Kimmeridgian sedimentary cycles. In *Orbital Forcing Timescales and Cyclostratigraphy* (eds MR House and AS Gale), pp. 75-114. Geological Society of London, Special Publication no. 85.

Waterhouse HK (1999) Orbital forcing of palynofacies in the Jurassic of France and the United Kingdom. *Geology* **27**, 511–4.

Weedon GP, Coe AL and Gallois RW (2004) Cyclostratigraphy, orbital tuning and inferred productivity for the type Kimmeridge Clay (Late Jurassic), Southern England. *Journal of the Geological Society, London* **161**, 655–66.

Weedon GP, Jenkyns HC, Coe AL and Hesselbo SP (1999) Astronomical calibration of the Jurassic time-scale from cyclostratigraphy in British mudrock formations. *Philosophical Transactions of the Royal Society of London, Series A: Mathematical, Physical and Engineering Sciences* **357**, 1787–813.

Wignall PB and Newton R (2001) Black shales on the basin margin: a model based on examples from the Upper Jurassic of the Boulonnais, northern France. *Sedimentary Geology* **144**, 335–56.

Wignall PB, Sutcliffe OE, Clemson J and Young E (1996) Unusual shoreface sedimentology in the Upper Jurassic of the Boulonnais, northern France. *Journal of Sedimentary Research* **66**, 577–86.

Williams CJ, Hesselbo SP, Jenkyns HC and Morgans-Bell HS (2001) Quartz silt in mudrocks as a key to sequence stratigraphy (Kimmeridge Clay Formation, Late Jurassic, Wessex Basin, UK). *Terra Nova* **13**, 449–55.

Zanella E and Coward MP (2003) Structural framework. In *The Millennium Atlas: Petroleum Geology of the Central and Northern North Sea* (eds D Evans, C Graham, A Armour and P Bathurst), pp. 45–59. London: Geological Society of London.

Ziegler PA (1990) *Geological Atlas of Western and Central Europe*. The Hague: Shell Internationale Petroleum Maatschappij B.V., 239 pp.