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

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 (δ13Corg) 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 δ13Corg 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.


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 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 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 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 & 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 . 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, 1995Waterhouse, , 1999Weedon 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 (δ 13 C org ) in order to correlate the Boulonnais outcrops to the regional record from the KCF ). 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 Kimmeridgianearliest Tithonian period that is applicable on a basin-wide scale.

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).

580
IC van der Hoeven et al.
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 fossils 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 shallowingupwards 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(Hatem et al. , 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(Hatem et al. , 2016Tribovillard 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 (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   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. 1991Proust 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 2 = 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 T max 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. , 2001. Cross-bedding occurs throughout the Grès de la Crèche Fm, as well as abundant bivalves and trace fossils 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 2 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    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.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 palynology and stable isotopes (Table 1; online Supplementary Material). The ACF was logged and sampled in the Cran aux OEufs section, situated in the north of the Boulonnais area (Fig. 3), GPS coordinates 50°51 0 1.23″ N and 1°34 0 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 OEufs 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 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 (δ 13 C org ). 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, δ 13 C V-PDB = of −26.43 ‰) was used as a reference material. For quality control purposes, check samples of beet sugar (IA-R005, δ 13 C V-PDB = −26.03 ‰) and cane sugar (IA-R006, δ 13 C V-PDB = −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. 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.   2001); it is a composite δ 13 C org record through the KCF of the Swanworth Quarry 1 and Metherhills boreholes. 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.
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 botanical affinity of this group is unclear.

4.a. TOC and stable isotope trends
The TOC contents in the Cran aux OEufs section range between 0.3 and 7 % ( Fig. 6; (Fig. 6, dashed line). The measured δ 13 C org values of the samples from the Cran aux OEufs 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 δ 13 C org (Fig. 6), but the remainder of the interval (Unit III) shows no correlation (R 2 of 0.0052; Fig. 7a). Conversely, in the interval that covers units IV-VII, a good fit exists between the TOC and the δ 13 C org isotope trend: high TOC content correlates with positive shifts in δ 13 C org values, with a higher R 2 of 0.35 and a slope of 0.25 (Fig. 7b). The positive correlation between TOC and the δ 13 C org values is most evident in the 2019 sample set from the Cap de la Crèche section, with an R 2 of 0.92 and a comparable slope of 0.37 (Figs 7, 8).

4.b. Palynological results
The palynological assemblages from the Cran aux OEufs 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)

590
IC van der Hoeven et al.
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 OEufs 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 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. (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.

592
IC van der Hoeven et al.
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.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 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 Herbin et al. 1995;Proust et al. 1995), which all indicate that the Irius Subzone is absent in the Boulonnais region.  . This latter section is located near the Cran aux OEufs 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 OEufs 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 OEufs section. The presence of the Irius Subzone is likely limited to the region around the Cran aux OEufs section. As indicated by Wignall & Newton (2001), the most complete record of the ACF in terms of preserved parasequences lies near the Cran aux OEufs section. The stable isotope correlation panel between Cran aux OEufs, Cran du Noirda and Cap de la Crèche shows a similar result, with the thickest lower ACF sequence preserved at Cran aux OEufs. 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. , 2001Herbin et al. 1995;Tribovillard et al. 2001Tribovillard et al. , 2019Braaksma et al. 2006;Hatem et al. 2017), grain size distribution (Williams et al. Climate-controlled organic matter accumulation in the Argiles de Châtillon Formation 593 2001; Braaksma et al. 2006;Hatem et al. 2017), sedimentary structures such as HCS (Wignall & Newton, 2001;Hatem et al. 2017;Tribovillard et al. 2019) and on the presence of an erosional unconformity halfway up the ACF (Wignall & Newton, 2001), it is evident that the relative sea level scenario involves two transgressiveregressive 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 OEufs 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.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 methods 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 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 mixedlayer 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 conditions 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 OEufs section. The

594
IC van der Hoeven et al.
might have similarly influenced the clay mineral signal via climatecontrolled 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 δ 13 C org , 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 palynological data (e.g. Waterhouse, 1999;Bonis et al. 2010). Besides, similar-scale cyclic variations have been described in the literature for the ACF and KCF, although with varying interpretations on the duration and underlying mechanism Waterhouse, 1995Waterhouse, , 1999Weedon et al. 1999Weedon et al. , 2004Armstrong et al. 2016). For example, Waterhouse (1999) observed palynofacies 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 δ 13 C org 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 vegetation changes that are characteristic of a changing Late Jurassic climate. 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 OEufs 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 δ 13 C org values for the Cap de la Crèche section are higher than the values of the Cran aux OEufs 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 OEufs section, and the depositional setting of the Kimmeridge Clay section is the most distally marine.
When looking at the variation of the TOC and δ 13 C org 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 δ 13 C org 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 δ 13 C org and TOC for the different ammonite stages ranging from Eudoxus to Scitulus show that the correlation between TOC and δ 13 C org becomes more evident after the Kimmeridgian-Tithonian boundary. In fact, for the Eudoxus and Autissiodorensis ammonite zones the correlation is weak, whereas the R 2 for the Elegans and Scitulus zones are higher in both the Boulonnais sections and the KCF (Fig. 7d-g). A closer analysis of the δ 13 C org and TOC records of the KCF (Fig. 6) shows that the positive correlation between TOC and δ 13 C org 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 δ 13 C org records between the ACF and KCF, it is evident that the δ 13 C org 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, δ 13 C org 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 δ 13 C org and a cyclic pattern of the TOC curve, which was absent during the Kimmeridgian. Apparently, local factors that could have the influenced the δ 13 C org 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 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 organicrich sediments are enhanced productivity of OM or enhanced preservation of OM. Most productivity estimates for the KCF and ACF indicate that average productivity was in both cases not increased (Tribovillard et al. 2001(Tribovillard et al. , 2005(Tribovillard et al. , 2019Weedon 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 conditions are plentiful for the KCF (van Kaam-Peters et al. 1997, 1998Saelen et al. 2000;van Dongen et al. 2006), but these indications are virtually absent for the ACF (Tribovillard et al. 2005(Tribovillard et al. , 2019Hatem et al. 2017). Despite these differences, sedimentological, palynological and palaeontological evidence shows that conditions 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 δ 13 C org 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;Saelen 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 van Kaam-Peters et al. 1998;van Dongen et al. 2006). In most environments, carbohydrates are rapidly remineralized and therefore contribute 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 sulfurization (Sinninghe Damsté & De Leeuw, 1990). The δ 13 C org signature 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 δ 13 C org 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 δ 13 C org values as observed in the KCF.
Based on the positive correlation between TOC, organic sulfur and HI, Tribovillard et al. (2001Tribovillard et al. ( , 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 δ 13 C org 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 δ 13 C org , 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 δ 13 C org value of −28 ‰ (Tyson, 1995;van Dongen et al. 2006). Type III terrestrial OM had a δ 13 C org signature ranging between −22 ‰ and −25 ‰ (Tyson, 1995;Pearce et al. 2005;van Dongen et al. 2006). The difference in δ 13 C org 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 δ 13 C org 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 δ 13 C org 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 terrestrial OM to the palynofacies and palynomorphs rarely exceeds 35 % (Fig. 8). This implies that the δ 13 C org 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 selective 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 δ 13 C orgto-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. (2005Tribovillard et al. ( , 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 conditions. 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 δ 13 C org . Here, there is no strong palynological evidence for stratification of the water column, in agreement with the trace-element analysis of Tribovillard et al. (2005Tribovillard et al. ( , 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 circumstances, sulfurization of OM, in particular carbohydrates, took place, as evident from the positive correlation between TOC and δ 13 C org . 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 assemblage and the clay mineral ratios. An eccentricity-modulated shift and expansion of the subtropical Hadley Cell was recently hypothesized to have affected the Late Jurassic climate, resulting in alternating 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 sulfurization (Hartgers et al. 1997;Tribovillard et al. 2015;Shawar et al. 2018;Raven et al. 2021). As proposed for the KCF, the 596 IC van der Hoeven et al.
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.

Conclusions
Palynological, TOC and δ 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 δ 13 C org . Parallel to the palynological record, a high-resolution δ 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 δ 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 δ 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 δ 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 basinwide 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.  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. PPprimary productivity; TOCtotal organic carbon.