# High-resolution record reveals climate-driven environmental and sedimentary changes in an active rift

Lisa Mcneill, Donna Shillington, Gareth Carter, Jeremy Everest, Robert Gawthorpe, Clint Miller, Marcie Phillips, Richard Collier, Aleksandra Cvetkoska, Gino de Gelder, et al.

► To cite this version:

Lisa Mcneill, Donna Shillington, Gareth Carter, Jeremy Everest, Robert Gawthorpe, et al. High-resolution record reveals climate-driven environmental and sedimentary changes in an active rift. Scientific Reports, 2019, 9 (1), pp.3116. 10.1038/s41598-019-40022-w. hal-02133929

HAL Id: hal-02133929

https://hal.umontpellier.fr/hal-02133929

Submitted on 18 May 2020

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
OPEN

# High-resolution record reveals climate-driven environmental and sedimentary changes in an active rift

Received: 28 June 2018

Accepted: 28 January 2019

Published online: 28 February 2019

Lisa C. McNeill ID <fnref n="1" />, Donna J. Shillington ID ², Gareth D. O. Carter ³, Jeremy D. Everest ³, Robert L. Gawthorpe ID ⁴, Clint Miller ⁵, Marcie P. Phillips ⁶, Richard E. LI. Collier ID ⁷, Aleksandra Cvetkoska ⁸, Gino De Gelder ⁹, Paula Diz ¹⁰, Mai-Linh Doan ID ¹¹, Mary Ford ¹², Maria Geraga ¹³, Jack Gillespie ID ¹⁴, Romain Hemelsdaël ¹⁵, Emilio Herrero-Bervera ¹⁶, Mohammad Ismaiel ¹⁷, Liliane Janikian ¹⁸, Katerina Kouli ID ¹⁹, Erwan Le Ber ²⁰, Shunli Li ²¹, Marco Maffione ²², Carol Mahoney ⁷, Malka L. Machlus ID ²,²³, Georgios Michas ²⁴, Casey W. Nixon ⁴, Sabire Asli Oflaz ID ²⁵, Abah P. Omale ²⁶, Kostas Panagiotopoulos ID ²⁷, Sofia Pechlivanidou ID ⁴, Simone Sauer ID ²⁸, Joana Seguin ²⁹, Spyros Sergiou ¹³, Natalia V. Zakharova ID ³⁰ & Sophie Green ³

Young rifts are shaped by combined tectonic and surface processes and climate, yet few records exist to evaluate the interplay of these processes over an extended period of early rift-basin development. Here, we present the longest and highest resolution record of sediment flux and paleoenvironmental changes when a young rift connects to the global oceans. New results from International Ocean Discovery Program (IODP) Expedition 381 in the Corinth Rift show 10s–100s of kyr cyclic variations in basin paleoenvironment as eustatic sea level fluctuated with respect to sills bounding this semi-isolated basin, and reveal substantial corresponding changes in the volume and character of sediment delivered into the rift. During interglacials, when the basin was marine, sedimentation rates were lower

<footnote n="1">¹School of Ocean and Earth Science, University of Southampton, Southampton, S014 3ZH, United Kingdom.
²Lamont-Doherty Earth Observatory of Columbia University, 61 Route 9W, Palisades, NY, 10964, USA. ³British Geological Survey, The Lyell Centre, Edinburgh, United Kingdom. ⁴Department of Earth Science, University of Bergen, Bergen, Norway. ⁵Department of Earth, Environmental and Planetary Sciences, Rice University, Houston, USA. ⁶Institute for Geophysics, University of Texas at Austin, Austin, USA. ⁷School of Earth and Environment, The University of Leeds, Leeds, United Kingdom. ⁸Department of Animal Ecology and Systematics, Justus Liebig University, Giessen, Germany. ⁹Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, Paris, France. ¹⁰Departamento Geociencias Marinas v Ordenación del Territorio, Facultad de Ciencias del Mar, Universidad de Vigo, Vigo, Spain. ¹¹Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, IRD, IFSTTAR, and ISTERre, Le Bourget-du-Lac, France. ¹²CRPG, UMR 7358, France. Also at: Université de Lorraine, ENSG, INP, Nancy, France. ¹³Department of Geology, University of Patras, Patras, Greece. ¹⁴Center for Tectonics, Resources, and Exploration (TRaX), Department of Earth Sciences, School of Physical Sciences, University of Adelaide, Adelaide, Australia. ¹⁵Géosciences Montpellier, Université de Montpellier, Montpellier, France. ¹⁶University of Hawai'i at Manoa, Hawai'i Institute of Geophysics and Planetology, Honolulu, USA. ¹⁷University Centre for Earth and Space Sciences, University of Hyderabad, Hyderabad, India. ¹⁸Universidade Federal de São Paulo, Departamento de Ciências do Mar, São Paulo, Brazil. ¹⁹Department of Geology and Geoenvironment, National and Kapodistrian University of Athens, Athens, Greece. ²⁰School of Geography, Geology and the Environment, University of Leicester, Leicester, United Kingdom. ²¹School of Energy Resources, China University of Geosciences, Beijing, China. ²²School of Geography Earth and Environmental Sciences, University of Birmingham, Birmingham, United Kingdom. ²³Department of Physical Sciences, Kingsborough Community College, City University of New York, New York, USA. ²⁴Laboratory of Geophysics and Seismology, Technological Educational Institute of Crete, Irakleio, Greece. ²⁵Graduate School "Human Development in Landscapes", Christian-Albrechts-Universität zu Kiel, Kiel, Germany. ²⁶Department of Geology and Geophysics, Louisiana State University, Baton Rouge, USA. ²⁷Institute of Geology and Minerology, University of Cologne, Cologne, Germany. ²⁸Ifremer, Department of Marine Geosciences, Centre Bretagne, Plouzané, France. ²⁹Institute for Ecosystem Research, Christian-Albrechts-Universität zu Kiel, Kiel, Germany. ³⁰Department of Earth and Atmospheric Sciences, Central Michigan University, Mount Pleasant, USA. Correspondence and requests for materials should be addressed to L.C.M. (email: lcmn@noc.soton.ac.uk)</footnote>
(excepting the Holocene), and bioturbation and organic carbon concentration higher. During glacials, the basin was isolated from the ocean, and sedimentation rates were higher (~2–7 times those in interglacials). We infer that reduced vegetation cover during glacials drove higher sediment flux from the rift flanks. These orbital-timescale changes in rate and type of basin infill will likely influence early rift sedimentary and faulting processes, potentially including syn-rift stratigraphy, sediment burial rates, and organic carbon flux and preservation on deep continental margins worldwide.

Active continental rift zones generate rapidly subsiding basins with significant accumulations of sediments. These settings are thought to be highly sensitive to the interplay of extensional tectonics, sedimentary processes, climate and sea level change¹,². Rift basins and other basins in active environments close to sea level are known to experience dramatic changes in environment due to glacio-eustatic fluctuations³,⁴. Both pronounced environmental change and active faulting are thought to control sediment delivery and accumulation in active basins¹⁻³. This includes changes in temperature, precipitation and amount and type of vegetation cover affecting erosion and sediment flux, and tectonics driving changes in catchment relief and area, erosion of different lithologies and shifting loci of sedimentation within a basin. However, the paucity of high-resolution constraints makes it difficult to isolate each of these contributions. Erosion and deposition of sediments also appear to have profound effects on rift localisation, fault longevity, crustal creation, and thermal structure⁵⁻⁹. However, the record of these processes is deeply buried and inaccessible in most ancient and mature rifts; high resolution data on age and rates of process are not available for deep and old sediments, the spatial resolution of seismic data is low, and sections are deeply buried and often overprinted by later deformational phases so do not preserve information on early fault history and its link to sedimentary processes. Basins in active environments tend to only have very short (mostly <15–25 ka) records from piston cores. Rare longer temporal records exist (Black Sea, Lake Malawi²,¹⁰) but in environments of lower tectonic activity or without ocean connection. Therefore, our ability to constrain in sufficient resolution temporal and spatial variations in basin paleoenvironment, sedimentation patterns or rift development, or the role of climate and sea level in moderating these factors is severely restricted. As a result, many hypotheses tend to be model derived and remain untested.

The Corinth Rift (Fig. 1) is a region of rapid, localised extension and high seismic activity. Current extension rates reach 10–15 mm/yr¹¹⁻¹³, some of the highest in the world. Corinth's high rates of tectonic activity, high sediment fluxes, closed drainage system and preservation of the syn-rift record make it a unique laboratory for the study of extension, sedimentation and paleoenvironment in a young rift.

Rifting began ~5 Ma with 3 main phases identified by integration of onshore depositional records and offshore seismic stratigraphy¹,¹⁴⁻²¹. Following an early continental and lacustrine phase (phase 1), depocentre deepening occurred ~2 Ma (phase 2) resulting in increased subsidence and sediment supply at the location of the modern Gulf of Corinth. At ~0.5–0.8 Ma, fault activity generally stepped northward, establishing the dominant N-dipping rift fault system along the southern boundary (Fig. 1) and development of fan deltas along the western part of the southern Gulf margin (phase 3). The transition from Phase 2 to 3 may mark the onset of repeated connection to the open ocean, although earlier marine incursions are recognised onshore²¹.

From onshore deposits and seismic imaging of offshore stratigraphy, the rift basin environment during phase 3 has been interpreted to alternate between marine conditions during interglacial/highstand periods and isolated from the open ocean during glacial/lowstand periods, as eustatic sea level fluctuated relative to the boundaries of the basin in the Late Quaternary¹⁷⁻¹⁹,²¹⁻²⁵. Sills at the western mouth of the Corinth Basin (Rion sill: ~60 m; Acheloos-Cape Pappas sill: ~50–55 m) currently control the connection with the Ionian Sea (Fig. 1) and are interpreted to control connection for at least the last 200 kyr²²,²⁶,²⁷. Prior to this, the cyclical connection to the ocean was either controlled by these sills in the west or by one or more sills in the east, primarily the Corinth Isthmus (Fig. 1), currently above sea level but submerged in the past²⁸. Prior to IODP Expedition 381, the only data available to constrain paleoenvironment evolution, and thus to examine alternating conditions, were shallow piston cores recovering sediments no deeper than 30 metres below seafloor (mbsf) and no older than 50 ka²⁹,³⁰. These data showed a change in sediment facies and micropaleontological assemblages at ~12 ka, marking the timing of transgression above the basin sill following the last glacial maximum. To date, these short records and similar records from other active basins do not provide a sufficiently long-duration record to evaluate sea-level driven changes in environment and their implications for sedimentation and rift processes.

IODP Expedition 381 drilled and cored the most recent ~1–2 Myr of syn-rift sediments to a depth of 705 mbsf, in October-December, 2017 (Fig. 1)³¹. This is the longest and highest resolution record of its kind in a young extensional basin at the point of connection to the global oceans, and it provides the first constraints on the age of the full rift sequence, syn-rift stratigraphy, rates and timings of rift tectonic processes, sediment fluxes and basin environmental conditions. Here we use this new record to test the hypothesis that glacial-interglacial timescale (10's-100 kyr) climatic and environmental change strongly influences the nature and volume of accumulating sediment within the Corinth basin. If such change occurs, we explore how sediment delivery and accumulation are affected, including effects on sediment fluxes, grain size and lithostratigraphy. Alternatively, the basin is insensitive to regional and local short-term environmental and climatic change being, instead, dominated by longer term tectonic-driven changes. Finally, we hypothesise that regional climatic change at 10–100 kyr timescales results in major basin environmental change with a corresponding change in aquatic biota.
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Map of the Gulf of Corinth and Corinth Rift system. Includes primary active rift faults (after ref. ¹⁷), positions of bounding sills to the Gulf, IODP drill sites, and primary catchment areas and rivers feeding the offshore basin (from ref. ⁴²). Inset shows regional setting.</figcaption></figure>

## Results

Three sites were drilled during Expedition 381 (Fig. 1). Site M0079, the focus of observations in this paper, is in the central Corinth basin, the primary depocentre of the Gulf of Corinth (Fig. S4). The nature of the site, within a primary locus of sedimentation, is based on analysis of syn-rift sediment thicknesses from integration of seismic profiles around the rift¹⁷. This site samples an expanded section of the most recent rift phase (phase 3), and thus provides a high-resolution record of extension, sedimentation and sea level change over the last ~750 ka. The site contains a thick, continuous succession of fine-grained distal facies and no faults. Thus, the depositional history including sediment accumulation rates (both absolute and variations through time) from this borehole are representative of the Corinth primary depocentre. The upper lithostratigraphic Unit 1, the focus of this paper, is subdivided into 16 subunits and has its base at 677 mbsf (Fig. 2), with a thin Unit 2 section beneath.

**Evidence for cyclical rift basin paleoenvironment.** Microfossil assemblages in Unit 1 alternate primarily between “marine” assemblages composed exclusively of marine microfossils, and assemblages with significant amounts of both non-marine and marine microfossils that reflect a range of complex conditions that are neither fully marine nor freshwater, with different degrees of marine influence (hereafter, “isolated”). During marine intervals, both marine and terrestrial microfossils are abundant, including calcareous nannofossils, marine diatoms, planktic and benthic foraminifera, dinoflagellate cysts, foraminifera test linings, pollen and spores (Figs 2 and 3, Table S1). During isolated intervals, non-marine diatoms are typically present in moderate to high abundances (Figs 2 and 3, Table S1), and observed in combination with moderate to high abundances of green algae coenobia and spores, marine or brackish dinoflagellate cysts, pollen and spores and low abundances of marine microfossils. The variability of microfossil and palynomorph assemblages within and between marine and isolated intervals directly indicates changes in aquatic basin paleoenvironment (e.g., salinity) and in terrestrial environment around the rift basin, paleoclimate and depositional setting. Based on these assemblages, the 16 subunits of Unit 1 were identified as marine or isolated intervals as a function of connection to or isolation from the open ocean. A combination of lithology, microfossil assemblages and core physical properties data provided the basis for pinpointing the subunit boundaries³¹.

**Mid-Late Quaternary Chronostratigraphy.** Biostratigraphy and magnetostratigraphy provide age control that allows us to link changing paleoenvironment to the eustatic sea level curve (Fig. 2). Calcareous nannofossils in the marine intervals confirm that Unit 1 is Middle Pleistocene through Recent based on (1) the Last Occurrence (LO) of *Pseudoemiliania lacunosa* (0.43 Ma) at 496 mbsf; (2) the First Occurrence (FO) of *Emiliania huxleyi* (0.29 Ma) at 405 mbsf; and (3) the crossover in dominance between *E. huxleyi* and *Gephyrocapsa* “small” (<3 µm) within subunit 1–2 (Fig. 2), documented at ~70 ka in the Mediterranean³². Due to the fluctuating environment within the basin, the syn-rift stratigraphy does not contain a continuous marine section, thus the LO of *P. lacunosa* and FO of *E. huxleyi* may not mark the true respective LO and FO (see Methods). However, the observed specimens are well preserved and moderately to highly abundant, and support other age markers. Magnetostratigraphy indicates the Brunhes-Matuyama chron boundary (0.773 Ma) occurs at ~665 mbsf at the base of Unit 1 (Fig. 2).

The overall pattern of downhole cyclicity of the Unit 1 subunits (marine vs isolated) strongly supports the palaeoenvironment being predominantly controlled by eustatic sea level fluctuation; tectonic control of this cyclical
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Results from Site M0079. (A) Age constraints from calcareous nannofossils and magnetostratigraphy (black text) and interpreted ages from correlation of basin environment with eustatic sea level curve (grey dotted lines and text, see Methods). (B–M) Brunhes-Matuyama chron boundary; E. Hux.: first occurrence (FO) of *E. huxleyi*; P.Lac: last occurrence (LO) of *P. lacunosa*; E.Hux/Gp. crossover: crossover in dominance between *E. huxleyi* and *Gephyrocapsa*. Arrows indicate potential uncertainty of true depth position of age marker (see text and Methods), and age is given in parentheses; (B) Oxygen isotope stages (OIS); (C) Lithostratigraphic units and subunits, coloured by basin environment interpretation (blue: marine; green: isolated/semi-isolated; grey: undetermined), subunit 1–11 is slumped interval (orange), yellow denotes Unit 2; (D) Abundance counts of calcareous nannofossils (light blue), benthic foraminifer (dark blue, dotted) and non-marine diatoms (green, dashed); (E) Estimated sand percent; (F) Bioturbation intensity; (G) Decompressed sedimentation rate (see Methods); and H) Sediment lithology coloured according to facies association, longer bars representing coarser intervals (Table S3; ref. ³¹).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Core images. (a) Core section examples (from Site M0079) of typical sedimentary facies in marine (left; Facies Association 1, homogeneous mud) and isolated (right; Facies Association 4, laminated greenish grey to grey mud with mud beds) intervals; (b) Typical microfossil assemblages from marine (above) and isolated (below) intervals (from Site M0079). Marine assemblage image shows planktic forams from >125 µm fraction; isolated assemblage (at 1000x magnification) shows non-marine planktonic diatom taxon *Pantocsekiella ocellata*.</figcaption></figure>

nature can be ruled out. In addition, the marine phases include short-lived isolated intervals, and vice versa, as expected from details of the eustatic fluctuations (e.g., within subunit 1–2, Fig. 2). We correlated the transitions between marine and isolated intervals to the eustatic sea level curve of Spratt and Lisecki³³ assuming a bounding sill depth of –60 m (present day Rion sill depth²⁷) and that flooding of this sill marks the transition between isolated and marine conditions (Fig. 2). This correlation and a sill depth of –60 m agree with age constraints from bio- and magnetostratigraphy and work exceptionally well to ~545 mbsf and the base of subunit 1–13, which is correlated with oxygen isotope stage (OIS) 13 at ~535 ka (Fig. 2). Below this depth, marine fauna occur in thin stratigraphic intervals within subunits 1–15 and 1–16 (Fig. 2). This may be explained by the basin-controlling sill being at a shallower level at this earlier stage, resulting in a limited record of earlier highstands in the basin.

Currently the primary age uncertainty derives from the presumption of a constant sill depth. As above, the nature of the marine record in the deeper part of Unit 1 suggests the sill was potentially shallower at this time. However, the persistence of the marine intervals and their relatively consistent thickness until we reach the deeper part of Unit 1 gives strong support to the controlling sill remaining close to this depth, supported by previous studies²⁷. A significant shift of the sill to shallower depths (e.g., 0–20 m) would severely restrict marine incursions, and this is not observed (unless we invoke very significant differences in sedimentation rate between marine intervals with time). Shallowing the sill level by up to 30 m does not significantly change the inferred age model, and the relative differences in sedimentation rates between isolated and marine intervals persist (see Methods, Fig. S3). In the absence of additional constraints, a constant sill depth is the most conservative assumption. As further research is undertaken on these drill cores, we will be able to add further absolute age constraints, including from tephra that are starting to be identified. This will reduce uncertainties of the age model and in turn provide constraints on sill height.

**Syn-rift sedimentary stratigraphy and processes.** Fine-grained (mud dominated), carbonate-rich turbiditic and hemipelagic sediments dominate the syn-rift succession of the Corinth Rift basin. However, we observe significant differences in lithology and depositional processes between the marine and isolated intervals of lithostratigraphic Unit 1. Marine subunits of Site M0079 are moderately to highly bioturbated and dominated by homogeneous to poorly-bedded greenish grey muds (Figs 2 and 3a). Biogenic material is common and
comprises calcareous nannofossils, foraminifera and marine diatoms, and fragments of gastropods and bivalves. Marine subunits are also characterised by increased total organic carbon (TOC) concentration with respect to the isolated subunits³⁴. Isolated subunits are characterised by laminated to thinly bedded and homogeneous grey and greenish grey muds (Fig. 3a), some with black, organic-rich laminations and beds, and are generally lacking or having sparse bioturbation (Fig. 2). Overall, isolated intervals contain a higher proportion of relatively coarser grained lithologies (silts and sands). This is particularly well expressed in the deeper subunits of Unit 1, for example compare sand content of marine subunits 1–5, 1–7 and 1–9 with isolated subunits 1–6 and 1–8 (Fig. 2).

We calculated sedimentation rates for each marine and isolated subunit using decompacted thicknesses of total sediment (hemipelagic plus gravity flow deposit) derived using borehole porosity data (see Methods)³¹. The primary uncertainties in these calculations include the sill depth, which controls the precise estimated timing of the transition between marine and isolated intervals and their duration, potential differences in seafloor porosity in the past, and the applicability of the decompaction curve used, particularly for the marine versus isolated intervals, which may compact differently (see Methods). Unit 1 has an average sedimentation rate of 1.1 mm/yr, but a striking result is the clear difference in sedimentation rate between the marine and isolated intervals/subunits within Unit 1 (Fig. 2). Sedimentation rates in isolated intervals/subunits (range of averages per interval: 1.7–3.3 mm/yr) are generally higher than those of the marine intervals/subunits (range of averages per interval: 0.3–0.7 mm/yr). For each pair of successive marine-isolated intervals (e.g., between subunits 1–3 and 1–2), excluding the Holocene subunit 1–1, the rates are 1.9–6.7 times greater in the isolated than the marine intervals. This result is supported by indicators³⁵ of enhanced fluvial input into the basin during glacial periods. Sedimentation rates for the 0–12 ka (Holocene) marine interval (subunit 1-1) are unusually high (average 2.9 mm/yr) compared to older marine intervals and are more similar to those of the isolated intervals. This suggests that the Holocene section recovered by piston cores in the Gulf of Corinth cannot be used to generalise conditions farther back in time. Overall there is also a slight pattern of increased rates from the past to present for both isolated and marine intervals. Deviating from this pattern, the sedimentation rate is particularly high in subunit 1–6 (MIS 8) at 3.3 mm/yr, and this is considerably higher than the prior and subsequent isolated subunits (Fig. 2).

Translating results from Site M0079 to the thickest Unit 1 section of the whole rift basin (1.5 km versus 0.68 km at Site M0079) using seismic stratigraphic unit thickness¹⁷ generates a maximum average sedimentation rate for Unit 1 in the basin of ~2 mm/yr and rates up to 3–4 mm/yr over individual ~100 kyr time periods. Overall these results indicate high sediment accumulation rates overall and large differences in sediment accumulation between marine/interglacial and isolated/glacial intervals.

## Discussion

**Changing environment in the Late Quaternary rift basin.** The changing environmental conditions reflected by the microfossil assemblages are interpreted to arise from fluctuating eustatic sea level with respect to the bounding sills of the Gulf of Corinth. Many previous studies have hypothesised that alternating marine and isolated conditions in the Gulf are driven by sea level change¹⁵,¹⁸,²⁵, but previous observations from this and other semi-isolated basins (e.g., Sea of Marmara, Black Sea) primarily came from piston cores and thus only sampled part of the most recent glacial cycle and the Holocene³,⁴,²⁹,³⁰. The cores recovered during IODP Expedition 381 provide the first direct evidence of sea-level driven changes in paleoenvironment in an active rift basin over hundreds of thousands of years (0–750 ka).

Observations from Site M0079 demonstrate profound climate and sea-level driven changes in the basin paleoenvironment, sediment accumulation and geochemistry³¹,³⁴ of this young rift over the last ~750 kyr. The microfossil assemblages are highly complex and variable – this is particularly true within the isolated intervals and within the marine-isolated transitions. We also note variations between the successive marine and isolated intervals. This variability suggests a much wider range of basin environments than simple marine versus freshwater end members for the marine and isolated intervals, respectively. The primary controls on the distribution of microfossil assemblages present are salinity, nutrient availability, pH, light, turbidity, temperature, and transport. Salinity, apparently the primary driver, is likely to be controlled by the combined effects of connection to the open ocean, freshwater influx/dilution, precipitation, evaporation, sill overspill from the Mediterranean Sea, and water body stratification and overturning that together result in a variety of conditions as eustatic sea level and climate fluctuated. Further research will investigate the specific environmental conditions indicated by the full suite of microfossils and how these vary through time.

Distinct differences in sedimentary processes and sedimentation rate are also observed between the marine and isolated intervals. The marine intervals are characterised by more homogeneous mud sequences, reduced coarse-grained sediment, and increased TOC concentration³⁴. The changes in accumulation of organic carbon between marine and isolated intervals result from the interplay between basin productivity, terrestrial carbon input, preservation and sediment flux, all driven by cycling climate and basin environment³⁴. Finally, the correspondence of increased bioturbation and homogeneity of sediments in the marine intervals indicates increased benthic organism activity. In summary, the clear increase in sedimentation rate during glacial periods when sea level is low and the basin is isolated (Fig. 2) implies that sediment fluxes (rates, volumes and grain size) into the basin from the subaerial rift flanks are increased during glacial times, probably leading to reduced benthic faunal activity within the basin.

**Influences on sediment flux into the basin.** Compared to other basins at similar evolutionary stages and over comparable time periods, sedimentation rates in Corinth (up to ~3–4 mm/yr) are high (e.g., Sea of Marmara, Lake Baikal, Malawi, Northern Gulf of California²,³,³⁶,³⁷). The Northern Gulf of California and Sea of Marmara have comparable sedimentation rates (up to 2–4 mm/yr³,³⁷). In Corinth, sediment flux is primarily controlled by combinations of high uplift rates of the rift flanks (Late Quaternary rates of 1–5 mm/yr¹⁵,³⁸–⁴¹) and erosion of weakly consolidated/lithified materials into accommodation space created through rapid extension
and subsidence. The ability of the fluvial systems to transport sediment is also an important control on delivery of sediments to the basin. However, in the Corinth system, fluvial systems do not appear to be limited by their ability to transport sediment³⁸, and thus the increased sediment fluxes during glacial periods reflect increased sediment production and supply.

Many of the most significant catchments are located on the southern rift margin, i.e., the N. Peloponnese (with the exception of the Mornos on the NW rift margin), coincident with highest uplift rates in the footwall of the southern rift boundary N-dipping fault system¹ (Fig. 1). In spite of the exceptionally high sedimentation rates shown here, the Corinth basin is presently underfilled, and it is likely that it has been underfilled for many 100's of kyr and potentially at least 1.5 Myr based on the height of delta foresets on the Gulf margins and those exposed onshore²¹. This highlights that rift basin subsidence rates exceed sediment accumulation rates, even during glacial/lowstand periods when sediment flux into the basin is enhanced by a factor of 2–7. This may be amplified by catchment averaged erosion rates being insufficient to keep pace with uplifting topography.

Although tectonic subsidence and uplift clearly control the larger-scale development of the basin, they would not be expected to fluctuate cyclically on timescales of 10's kyr, as sedimentation rates fluctuate here. Therefore, we propose that basin sediment accumulation, including variations between glacial and interglacial periods, is a function of one or more of the following: (a) climate (temperature, precipitation) controlling erosion rates and sediment flux²⁹,⁴²; (b) vegetation type and cover changing runoff and sediment retention/erosion rates; c) changing basin salinity driving more efficient hyperpycnal flows during isolated/glacial intervals¹. Because the shelf is narrow in much of the Gulf of Corinth, enhanced slope failure due to direct sediment supply to the shelf edge is unlikely to contribute and catchments are unlikely to significantly increase in area, during sea level lowstands. However, the fluvial and deltaic material exposed on the shelf during lowstands may be more easily erodable, and the partially exposed, relatively unstable margin slopes may be more susceptible to failure²⁷. A previous study proposed that a combination of cool, wet winters in the last glacial and reduction in tree cover led to increased sediment availability and flux during glacials estimated over the last ~100 kyr in the eastern Corinth rift²⁹. However, pollen-based precipitation reconstructions from the region⁴³⁻⁴⁵ indicate a decrease in precipitation during glacials, not an increase. Long pollen records across the Balkan Peninsula suggest interglacial forested landscapes alternating with open vegetation during glacials⁴⁴⁻⁴⁶, and our new drill cores in the Corinth basin show the same results. Our new cores also have lower pollen and higher corroded pollen grain concentrations in isolated intervals, the latter indicating increased reworking attributed to increased soil erosion, often linked to open vegetation⁴⁷. These pollen and sedimentation rate results from the new IODP borehole record and previous studies lead us to conclude that both reduction and change in type of vegetation cover resulted in increased soil erosion and higher sediment flux into the basin during the isolated/glacial intervals. Analysis of bedrock fault scarp weathering shows evidence for enhanced physical weathering during glacial periods in the Mediterranean region⁴⁸, which may add to the observed Corinth sediment flux increase. Contrasts in seasonality and storminess between glacials and interglacials may have also contributed to erosion rates, but this parameter has not yet been resolved from our new sedimentary record, or previous studies in the region. Sedimentation rates from piston cores and drilling have also revealed variations in sedimentation rates over glacial/interglacial cycles in the Sea of Marmara and Black Sea³,¹⁰.

The Holocene (last marine/interglacial interval) appears to have significantly higher sedimentation rates than earlier marine/interglacial intervals and is comparable to isolated/glacial intervals (Fig. 2). Forest clearance and agriculture are visible in pollen records from Greece since the mid-Holocene (~6 ka)⁴⁹, therefore it is likely that human activity in the study region is responsible for the high sediment fluxes observed in Holocene subunit 1–1⁵⁰.

**Implications for the early history of rifts.** The Corinth Rift basin provides a unique window into an active rift basin environment at the point in time where the rift basin connects to the global ocean system due to rifting and resulting subsidence. It therefore is intermediate between terrestrial rift basins that have yet to achieve a connection to the open oceans (e.g., Lake Baikal, the majority of the East African Rift system) and rifts that have fully connected and no longer oscillate between open and isolated basin environment (e.g., Gulf of California, passive rifted margins). The new Corinth boreholes provide an extended record with unprecedented resolution over multiple glacio-eustatic cycles and will add insight to how the environment of other semi-isolated basins will have developed over time, such as the Sea of Marmara and Black Sea. Our new results illustrate the significant impact of regional climate at orbital timescales (10 s–100 s kyr) on sediment accumulation, and how this in turn affects patterns of organic matter accumulation, preservation and burial history (see also ref.³⁴). We would expect to see similar cyclical patterns of organic carbon to those observed here in other rift basins regularly connecting to the oceans. Observed sedimentation patterns, controlled by sea level and climate, could have implications for the distribution of potential source and reservoir rocks within the early phase of sediment build-up on rifted margins, in particular where icehouse conditions prevailed. This is because sediment build up impacts the burial history, preservation and thermal evolution of deep rifted margin sediments and organic matter. Cyclical fluctuations in sediment accumulation combined with feedbacks on fault activity will have important implications for how tectonics and climate interact to control the stratigraphy of mature rifts and rifted margins worldwide⁵¹.

The observation of high average rates of sedimentation that vary substantially over ~10–100 kyr timescales could have important consequences for tectonic processes in Corinth and early rifting in general. Surface processes can impact lithospheric extensional processes, in particular, the redistribution of mass from the uplifted footwalls into the basin promotes greater extension on rift faults before they are abandoned in favour of new faults⁹. Thicker sediments could also reduce differences in buoyancy forces resulting from thinning, and thus promote narrow rifting⁵. Many of these inferences are based on modelling or larger scale structure but with limited temporal constraints. For the relatively high slip rates on faults bounding the Corinth Rift (2–10 mm/yr)¹⁵, the temporal changes in erosion onshore and deposition offshore observed here could modulate fault evolution.
Further application of the new drilling results in the Corinth rift to analysis of the fault activity will allow us to test whether sedimentation changes, even on relatively short timescales, can impact rift faulting.

## Conclusions

A new scientific drill core record from the syrift sequence of the active Corinth rift provides the first direct observations that basin environment fluctuated between marine conditions during eustatic highstands and isolated conditions during eustatic lowstands when the basin is cut off by bounding sills over the last ~700 ky. Sedimentation rates in the basin show significant variations on 10's-100 kyr timescales and are markedly increased during glacial/isolated periods. In contrast, bioturbation and organic carbon concentrations are reduced during glacial/isolated periods. The new borehole data, supported by other studies, suggest sedimentation rate changes are a function of a decrease and change in type of vegetation cover in glacial periods, resulting in increased erosion and basin sediment flux. The aquatic basin environment is clearly influenced by 10-100 kyr sea level fluctuation and climate, with the microfossil assemblages sampled indicating much greater complexity than simple alternations between fully marine to freshwater conditions. These results thus reveal the dominant role of climate and sealevel change in generating ~100-ky variations in sedimentation rate and basin environment in this active rift basin.

## Methods

**Micropaleontology and Basin Environment.** Microfossils (calcareous nannofossils, marine and non-marine diatoms, planktic and benthic foraminifer, dinoflagellate cysts, foraminifer test linings, freshwater algae coenobia and spores, and aquatic pollen and spores) were used to distinguish basin paleoenvironment, principally differentiating “marine” or “isolated”. Supplementary Table S1 gives abundances of key microfossil groups: calcareous nannofossils and benthic foraminifer (marine indicators) and non-marine diatoms (isolated indicators) with depth in Hole M0079A. Qualitative counts of calcareous nannofossils and non-marine diatoms are based on the Cascading Count Method⁵². The numerical approximations associated with abundance in Fig. 2 (abundance: Barren (B), Very Rare (VR), Rare (R), Few (F), Common (C) and Abundant (A)) for calcareous nannofossils and non-marine diatoms are outlined below:

Calcareous Nannofossils: $B = 0$; $VR = 1-5$; $R = 6-100$; $F = 101-1500$; $C = 1501-5000$; $A = 5001-10,000$+

Non-marine Diatoms: $B = 0$; $VR = 1$; $R = 2$; $F = 3-10$; $C = 11-20$; $A = 21-50$+

The abundance of benthic foraminifer is represented by the number of individuals found in ~10 cc of wet sediment. Benthic foraminifer were counted in the >125 µm fraction. In the isolated intervals, where abundance in the >125 µm is generally lower than 10 individuals, the 63-125 µm fraction was also screened.

**Age Model.** Age constraints were from shipboard biostratigraphic (calcareous nannofossils) and magnetostratigraphic analyses. Biozonation of calcareous nannofossils was applied using existing studies³,⁵,⁴. The website www.mikrotax.org aided identification of calcareous nannofossils. Calcareous nannofossils provide three age markers (shown on Fig. 2). Due to the fluctuating environment within the basin, the syn-rift stratigraphy does not contain a continuous marine section, thus the LO of *P. lacunosa* and FO of *E. huxleyi* may not mark the true respective LO and FO. We note that *P. lacunosa* was identified in a coherent interval within a large-scale slump that defines subunit 1–11 (Fig. 2)³¹ surrounded by intervals interpreted as isolated. The stratigraphic interval containing this marine calcareous nannofossil assemblage is intact with no evidence for reworking. Therefore, this coherent interval of slumped sediments containing *P. lacunosa* is interpreted to represent a part of an older marine interval, most likely the time equivalent of the underlying marine subunit 1–13. Preliminary magnetostratigraphy analysis provides one age marker in Hole M0079A, the Brunhes-Matuyama chron boundary at 0.773 Ma at a depth of 665 mbsf (Fig. 2). A total of 532 discrete sediment cubic samples were collected from working halves at intervals of ~1.5 m throughout the borehole. All samples were demagnetised using alternating field (AF) treatment in 14 progressive field steps from 5 to 40 mT (with 5-mT increments) and from 40 to 100 mT (with 10-mT increments). Remanent magnetisation direction and intensity were measured before and after each demagnetization step using the horizontal pass-through super-conducting cryogenic rock magnetometer (SRM 755–4000, 2G Enterprises) at the University of Bremen. The inclination of the remanence after AF demagnetisation at 40 mT was used to determine the polarity of each sample (i.e., normal or reversed) and build a magnetostratigraphy downhole. Magnetozones identified from the data were correlated to the Geomagnetic Polarity Time Scale – GPTS⁵⁴,⁵⁵.

In addition to the above age markers, the age model was developed by tying the Unit 1 subunit boundaries between the marine and isolated intervals to eustatic sea level (Supplementary Fig. S1 and Table S2³³). A sill depth of ~60 m (the current depth of the Rion sill at the western end of the Gulf of Corinth²⁷) was used to determine the transition timing between marine and isolated. See Supplementary Information for further discussion of sill depth and impact of its depth change on sedimentation rates. All ages between these transitions were extrapolated linearly. Below a depth of ~545 mbsf, marine intervals were thin or absent (Fig. 2). Therefore, the age model developed from microfossil-based sea level correlation could not be applied from here to the base of the hole; however, the Brunhes-Matuyama chron boundary at 665 mbsf provides absolute age constraint near the base of the hole. Ages and depths in Hole M0079A are shown in Supplementary Table S2.

**Facies Associations and Bioturbation.** The lithostratigraphy of the syn-rift sediments drilled during Expedition 381 were categorised into facies associations (FA³¹) defined by physical and biogenic features of the sediment, including bedding and lamination style, grain size, colour, body and trace fossils. The facies associations used in this paper are defined in Supplementary Table S3 and a simplified version is depicted in Fig. 2. The degree of bioturbation applied to the cores is a semi-quantitative assessment ranging from 0 (no bioturbation) to 6 (completely bioturbated)⁵⁶.
**Sedimentation Rate.** Sedimentation rates were calculated using ages described above and decompacted sediment thicknesses. Decompression was based on porosities measured using the "moisture and density" technique on discrete 6 cm³ samples spaced at ~1.5 m on cores from Hole M0079A. The wet and dry masses of these samples were measured before and after being dried in a convection oven at 60° ± 5°C for 24 hrs, respectively. The volume of dried sample was measured with a helium-displacement pycnometer. These measurements were then used to calculate the volume and mass of water originally in the samples, and the porosity of the samples (Supplementary Table S4³¹).

Next, we determined a smooth porosity function to use for decompaction. We first removed outliers by fitting a 2nd order polynomial to the measured porosities and discarding values with residuals greater than 1.5, and then fit a 35th order polynomial to the remaining points. A high-order polynomial was required to capture the observed variations in porosity between marine and isolated subunits (Supplementary Fig. S2 and Table S5).

The decompacted thickness of sediments from a given depth interval, $T_i^*$, was determined (Supplementary Table S5) assuming that there is not alteration of the grains:

$$T_i^* = \frac{T_i(1 - \phi_i)}{(1 - \phi_i^*)}$$

where $T_i$ and $\phi_i$ are the compacted thickness and porosity of a given interval, respectively, where the porosity is taken from the smoothed function. $\phi_i^*$ is the initial porosity, and is assumed to be 56%, the porosity at the modern seafloor in Hole M0079A (Supplementary Table S4³¹).

## Data Availability

All data and material pertinent to this study is contained within the manuscript and Supplementary Information, and/or the IODP Expedition 381 Preliminary Report (Shillington et al., 2019; ref.³¹). The full Expedition Report and dataset from IODP Expedition 381 will become openly available on March 1, 2019 (McNeill et al., 2019; ref.³⁴) via the IODP website (https://www.iodp.org/resources/access-data-and-samples), including access to core materials and logging data.

## References

1. Gawthorpe, R. L. et al. Tectono-sedimentary evolution of the Plio-Pleistocene Corinth rift, Greece. *Basin Research*, https://doi.org/10.1111/bre.12260 (2017).

2. Lyons, R. et al. Continuous 1.3-million-year record of East African hydroclimate, and implications for patterns of evolution and biodiversity. *PNAS* 51, 15568–15573 (2015).

3. McHugh, C. M.G. et al. The last reconnection of the Marmara Sea (Turkey) to the World Ocean: A paleoceanographic and paleoclimatic perspective. *Marine Geol.*, https://doi.org/10.1016/j.margeo.2008.1007.1005 (2008).

4. Mudie, P., Rochon, A., Aksu, A. & Gillespie, H. Dinoflagellate cysts, freshwater algae and fungal spores as salinity indicators in Late Quaternary cores from Marmara and Black seas. *Marine Geol.* 190, 203–231 (2002).

5. Bialas, R. W. & Buck, W. R. How sediment promotes narrow rifting: Application to the Gulf of California. *Tectonics* 28, https://doi.org/10.1029/2008TC002394 (2009).

6. Dorsey, R. & Umhoefer, P. J. In *Tectonics of Sedimentary Basins: Recent Advances*. (eds Busby, C. & Azor, A.) 209–225 (2012).

7. Han, L. et al. Seismic imaging of the metamorphism of young sediment into new crystalline crust in the actively rifting Imperial Valley, California. *Geochem. Geophys. Geosys.* 17, 4566–4584 (2016).

8. Lizarralde, D. et al. Variation in styles of rifting in the Gulf of California. *Nature* 448, 466–469 (2007).

9. Olive, J.-A., Behn, M. D. & Malatesta, L. C. Modes of extensional faulting controlled by surface processes. *Geophys. Res. Lett.* 41, 6725–6733 (2014).

10. Hsü, K. J. et al. In *Deep Sea Drilling Project Initial Reports*, 42, 509–524 (1978).

11. Bernard, P. et al. Seismicity, deformation and seismic hazard in the western rift of Corinth: New insights from the Corinth Rift Laboratory (CRL). *Tectonophysics* 426, 7–30 (2006).

12. Briole, P. et al. Active deformation of the Corinth rift, Greece: Results from repeated Global Positioning System surveys between 1990 and 1995. *J. Geophys. Res.* 105, 25605–25625 (2000).

13. Clarke, P. et al. Crustal strain in central Greece from repeated GPS measurements in the interval 1989–1997. *J. Geophys. Res.* 135, 195–214 (1998).

14. Armijo, R., Meyer, B., King, G., Rigo, A. & Papanastassiou, D. Quaternary evolution of the Corinth Rift and its implications for the Late Cenozoic evolution of the Aegean. *Geophys. J. Intl.* 126, 11–53 (1996).

15. Bell, R. E. et al. Fault architecture, basin structure and evolution of the Gulf of Corinth, central Greece. *Basin Research*, https://doi.org/10.1111/j.1365-2117.2009.00401.x (2009).

16. Higgs, B. Syn-sedimentary structural controls on basin deformation in the Gulf of Corinth, Greece. *Basin Research* 1, 155–165 (1988).

17. Nixon, C. et al. Rapid spatiotemporal variations in rift structure during development of the Corinth Rift, central Greece. *Tectonics* 35, 1225–1248, https://doi.org/10.1002/2015TC004026 (2016).

18. Sachpazi, M., Clement, C., Laigle, M., Hirn, A. & Roussos, N. Rift structure, evolution, and earthquakes in the Gulf of Corinth, from reflection seismic images. *Earth Planet Sci. Lett.* 2003, 243–257 (2003).

19. Taylor, B. et al. The structures, stratigraphy and evolution of the Gulf of Corinth rift, Greece. *Geophys. J. Intl.* 185, 1189–1219 (2011).

20. Ford, M., Hemelsdael, R., Mancini, M. & Palyvos, N. Rift migration and lateral propagation: evolution of normal faults and sediment-routing systems of the western Corinth rift (Greece). In: *The Geometry of Normal Faults* (Ed. by Childs, C., Holdsworth, R. E., Jackson, C.A.-L., Manzocchi, T., Walsh, J. J., Yielding, G.) *Geol. Soc. London, Spec. Publ.*, (439) London (2016).

21. Ford, M. et al. Tectono-sedimentary evolution of the western Corinth Rift (Central Greece). *Basin Res.* 25, 3–25 (2013).

22. Bell, R. E., McNeill, L. C., Bull, J. M. & Henstock, T. J. Evolution of the offshore western Gulf of Corinth. *GSA Bulletin* 120, 156–178 (2008).

23. Leeder, M. et al. Tectono-sedimentary processes along an active marine/lacustrine half-graben margin: Alkyonides Gulf, E. Gulf of Corinth, Greece. *Basin Research* 14, 25–41 (2002).

24. Lykousis, V., Sakellariou, D., Moretti, I. & Kaberi, H. Late Quaternary basin evolution of the Gulf of Corinth: Sequence stratigraphy, sedimentation, fault-slip and subsidence rates. *Tectonophysics* 440, 29–51 (2007).

25. Heezen, B. C., Ewing, M. & Johnson, G. L. The Gulf of Corinth floor. *Deep Sea Research* 13, 381–411 (1966).

26. Beckers, A. et al. Influence of bottom currents on the sedimentary processes at the western tip of the Gulf of Corinth, Greece. *Marine Geol.* 378, 312–332 (2016).
27. Perissoratis, C., Piper, D. J. W. & Lykousis, V. Alternating marine and lacustrine sedimentation during late Quaternary in the Gulf of Corinth rift basin, central Greece. *Marine Geol.* **167**, 391–411 (2000).

28. Collier, R. E. L. & Thompson, J. Transverse and linear dunes in an upper pleistocene marine sequence, Corinth basin, Greece. *Sedimentology* **38**, 1021–1040 (1991).

29. Collier, R. et al. High sediment yields and cool, wet winters: Test of last glacial paleoclimates in the northern Mediterranean. *Geology* **28**, 999–1002 (2000).

30. Moretti, I. et al. Sedimentation and subsidence rate in the Gulf of Corinth: what we learn from the Marion Dufresne's long-piston coring. *Tectonics* **336**, 291–299 (2004).

31. Shillington, D. J., McNeill, L. C. & Carter, G. D. O. & the Expedition 381 Participants. Expedition 381 Preliminary Report: Corinth Active Rift Development. *College Station, TX, International Ocean Discovery Program.*, https://doi.org/10.14379/iodp.pr.381.2018 (2019).

32. Anthonissen, D. E. & Ogg, J. *The Geologic TimeScale 2012*, 1083–1127 (2012).

33. Spratt, R. M. & Lisecki, L. E. A Late Pleistocene sea level stack. *Climate of the Past* **12**, 1079–1092 (2016).

34. McNeill, L. C., Shillington, D. J., Carter, G. D. O. & Expedition 381 Participants. Corinth Active Rift Development. *Proceedings of the International Ocean Discovery Program*, **381**, College Station, TX, International Ocean Discovery Program., https://doi.org/10.14379/iodp.proc.381.2019 (2019).

35. Toucanne, S. et al. Tracking rainfall in the northern Mediterranean borderlands during sapropel deposition. *Quat. Sci. Rev.* **129**, 178–195 (2015).

36. Edginton, D. N. et al. Sedimentation rates, residence times and radionuclide inventories in Lake Baikal from 137Cs and 210Pb in sediment cores. *Nature* **350**, 601–604 (1991).

37. Dorsey, R. Sedimentation and crustal recycling along an active oblique-rift margin: Salton Trough and northern Gulf of California. *Geology* **38**, 443–446 (2010).

38. Demoulin, A., Beckers, A. & Hubert-Ferrari, A. Patterns of Quaternary uplift of the Corinth rift southern border (N Peloponnese, Greece) revealed by fluvial landscape morphometry. *Geomorph.* **246**, 188–204 (2015).

39. Michas, G., Vallianatos, F. & Sammonds, P. Statistical mechanics and scaling of fault populations with increasing strain in the Corinth Rift. *Earth Planet Sci. Lett.* **431**, 150–163 (2015).

40. Palyvos, N., Pantosti, D. & Zabci, C. Paleoseismological evidence of recent earthquakes on the 1967 Mudurnu Valley earthquake segment of the North Anatolian fault. *Bull. Seis. Soc. Am.* **97**, 1646–1661 (2007).

41. Pirazzoli, P. A., Stiros, S. C., Fontugne, M. & Arnold, M. Holocene and Quaternary uplift in the central part of the southern coast of the Corinth Gulf (Greece). *Marine Geol.* **212**, 35–44 (2004).

42. Watkins, S. E. et al. Are landscapes buffered in high frequency climate change? A comparison of sediment fluxes and depositional volumes in the Corinth Rift, central Greece, over the past 130 kyrs. *Geol. Soc. Am. Bull.*, https://doi.org/10.1130/B31953.1 (2018).

43. Dormoy, I. et al. Terrestrial climate variability and seasonality changes in the Mediterranean region between 15 000 and 4000 years BP deduced from marine pollen records. *Climate of the Past* **5**, 615–632 (2009).

44. Sadori, L. et al. Climate, environment and society in southern Italy during the last 2000 years. A review of the environmental, historical and archaeological evidence. *Quat. Sci. Rev.* **136**, 173–188 (2016).

45. Tzedakis, P., Hooghiemstra, H. & Palike, H. The last 1.35 million years at Tenaghi Philippon: revised chronostratigraphy and long-term vegetation trends. *Quat. Sci. Rev.* **25**, 3416–3430 (2006).

46. Tzedakis, P. et al. Ecological thresholds and patterns of millennial-scale climate variability: The response of vegetation in Greece during the last glacial period. *geology* **32**, 109–112 (2004).

47. Wilmshurst, J. M. & McGlone, M. S. Corroded pollen and spores as indicators of changing lake sediment sources and catchment disturbance. *J. Paleolimnol.* **34**, 503–517 (2005).

48. Tucker, G. E. et al. 2011, Geomorphic significance of postglacial bedrock scarps on normal-fault footwalls. *J. Geophys. Res.* **116**, F01022 (2011).

49. Weiberg, E. et al. The socio-environmental history of the Peloponnese during the Holocene: Towards an integrated understanding of the past. *Quat. Sci. Rev.* **136**, 40–65 (2016).

50. Butzer, K. Environmental history in the Mediterranean world: cross-disciplinary investigation of cause-and-effect for degradation and soil erosion. *J. Archaeological Sci.* **32**, 1773–1800 (2005).

51. Balázs, A., Granjeon, D., Matenco, L., Sztanó, O. & Cloetingh, S. Tectonic and Climatic Controls on Asymmetric Half-Graben Sedimentation: Inferences From 3-D Numerical Modeling. *Tectonics* **36**, 2123–3141, https://doi.org/10.1002/2017TC004647 (2017).

52. Stzyzen, M. J. Cascading Counts of Nannofossil Abundance. *J. Nannoplankton Res.*, **19** (1) (1997).

53. Backman, J., Raffi, I., Rio, D., Fornaciari, E. & Pálike, H. Biozonation and biochronology of Miocene through Pleistocene calcareous nannofossils from low and middle latitudes. *Newsletters on Stratigraphy* **45**(3), 221–244 (2012).

54. Ogg, J., Ogg, G. & Gradstein, F. *A Concise Geologic Time Scale 2016*. Elsevier B.V., ISBN 978-0-444-63771-0 (2016).

55. Cande, S. C. & Kent, D. V. Revised calibration of the geomagnetic polarity timescale for the Late Cretaceous and Cenozoic. *J. Geophys. Res.* **100**(B4), 6093–6095, https://doi.org/10.1029/94JB03098 (1995).

56. Taylor, A. M. & Goldring, R. Description and analysis of bioturbation and ichnofabric. *J. Geol. Soc. London* **150**, 141–148, https://doi.org/10.1144/gsjgs.150.1.0141 (1993).

## Acknowledgements

We sincerely thank all involved with the successful completion of IODP Expedition 381, including ECORD Science Operator staff, ship and drilling crew of the D/V *Fugro Synergy*, and staff at MARUM, University of Bremen. We thank Bill Ryan for useful comments on an earlier version of the manuscript. We thank two anonymous reviewers for their very helpful reviews and comments that improved the paper.

## Author Contributions

McNeill, Shillington, Gawthorpe, Miller and Phillips synthesised results and prepared the publication. Phillips, Cvetkoska, Diz Ferreriro, Geraga, Kouli, Oflaz, Panagiotopoulos characterised paleoenvironment based on microfossil assemblages. Gawthorpe, Collier, De Gelder, Ford, Gillespie, Hemelsdaël, Janikian, Li, Nixon, Pechlivanidou and Sergiou characterised sediment lithology and assigned facies associations. Phillips undertook biostratigraphy with calcareous nannofossils, and Herrero-Bervera and Maffione undertook magnetostratigraphy. Le Ber, Doan, Ismaiel, Machlus, Michas, Omale and Zakharova acquired and analysed physical properties data, including density and porosity. Miller, Mahoney, Sauer and Seguin acquired and analysed geochemical data, including TOC. Miller and McNeill established ties to the sea level curve, and Miller and Shillington calculated decompaction and sediment rates. McNeill, Shillington, Carter, Everest and Green organised the expedition that produced these results. All authors generated and analysed core and log data during Expedition 381 and contributed to the preparation of this paper.
**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-40022-w.

**Competing Interests:** The authors declare no competing interests.

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

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

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

© The Author(s) 2019