Geochemistry of sediments from the upwelling off Morocco

The high-productive upwelling area off Morocco is part of one of the four major trade-wind driven continental margin upwelling zones in the world oceans. While coastal upwelling occurs mostly on the shelf, biogenic particles derived from upwelling are deposited mostly at the upper continental slope. Nutrient-rich coastal water is transported within the Cape Ghir filament region at 30°N up to several hundreds of kilometers offshore. Both upwelling intensity and filament activity are dependent on the strength of the summer Trades. This study is aimed to reconstruct changes in trade wind intensity over the last 250,000 years by the analysis of the productivity signal contained in the sedimentary biogenic particles of the continental slope and beneath the Cape Ghir filament. Detailed geochemical and geophysical analyses (TOC, carbonate, C/N, delta13Corg, delta15N, delta13C of benthic foraminifera, delta18O of benthic and planktic foraminifera, magnetic susceptibility) have been carried out at two sites on the upper continental slope and one site located further offshore influenced by the Cape Ghir filament. A second offshore site south of the filament was analyzed (TOC, magnetic susceptibility) to distinguish the productivity signal related to the filament signal from the general offshore variability. Higher productivity during glacial times was observed at all four sites. However, the variability of productivity during glacial times was remarkably different at the filament-influenced site compared to the upwelling-influenced continental slope sites. In addition to climate-related changes in upwelling intensity, zonal shifts of the upwelling area due to sea-level changes have impacted the sedimentary productivity record, especially at the continental slope sites. By comparison with other proxies related to the strength and direction of the prevailing winds (Si/Al ratio as grain-size indicator, pollen) the productivity record at the filament-influenced site reflects mainly changes in trade-wind intensity. Our reconstruction reveals that especially during glacial times trade-wind intensity was increased and showed a strong variability with frequencies related to precession.

Susceptibility measurements at Gubbio sections

The complex interplay between extraterrestrial events and earth-bound processes that triggered one of the greatest biological crises of the Phanerozoic requires a high resolution timescale. Detailed magnetic susceptibility measurements at the Contessa Highway and Bottaccione sections (Italy) span the Cretaceous-Paleogene boundary and reveal clear orbital signatures in the sedimentary record. Identification of precession and 405 kyr eccentricity cycles allows an estimate of 324+/-40 kyr for the duration of the Maastrichtian part of Chron C29r. We present in the same high resolution time frame sites in Spain and the North and South Atlantic and bio-horizons, biotic changes, stable isotopic excursions and the decrease in Osmium isotopes recorded in these sections. The onset of 187 Os/ 188 Os decrease coincides with the d13 C negative excursion K-PgE1, thus suggesting a first pulse in Deccan volcanism at 66.64 Ma. The K-PgE3 d13 C negative excursion is possibly the expression of a second pulse at 66.26 Ma. Late Maastrichtian d13 C negative excursions are of low intensity and span durations of one to two eccentricity cycles, whereas early Danian excursions are brief (about 30 kyr) and acute. In Biotic response to late Maastrichtian perturbations occurred with a delay of ca. 200 kyr after the beginning of K-PgE1 shortly before K-PgE3. The biotic perturbation could be thus either a delayed response to K-PgE1, or a direct response to K-PgE3, and possibly, a threshold response to the stepwise buildup of CO2 atmospheric injections. No delay is evident in response to early Danian hyperthermal events. These differences suggest that short-lived, volcanically-derived environmental perturbations were buffered within the stable late Maastrichtian oceanic realm whereas they were amplified by the more sensitive and highly disturbed early Danian oceanic ecosystem.

(Table 2) Age model of ODP Site 184-1146

Orbital tuning of benthic d18O is a common approach for assigning ages to ocean sediment records. Similar environmental forcing of the northern South China Sea and the southeast Asian cave regions allows for transfer of the speleothem d18O radiometric chronology to the planktonic and benthic d18O records from Ocean Drilling Program Site 1146, yielding a new chronology with 41 radiometrically calibrated datums, spanning the past 350 kyr. This approach also provides for an independent assessment of the accuracy of the orbitally tuned benthic d18O chronology for the last 350 kyr. The largest differences relative to the latest chronology occur in marine isotope stages (MIS) 5.4, 5.5, 6, 7, and 9.3. Prominent suborbital-scale structure believed to be global in nature is identified within MIS 5.4 and MIS 7.2. On the basis of the radiometrically calibrated chronology, the time constant of the ice sheet is found to be 5.4 kyr at the precession band (light d18O lags precession minima by -55.4°) and 10.4 kyr at the obliquity band (light d18O lags obliquity maxima by 57.4°). These values are significantly shorter than the single 17 kyr time constant originally estimated by Imbrie et al. (1984), based primarily on the timing of terminations I and II and the 15 kyr time constant used by Lisiecki and Raymo (2005, doi:10.1029/2004PA001071).

(Table 1) Age determination of sediment core SO139-74KL

A multiproxy record has been acquired from a piston core (SO139-74KL) taken offshore southern Sumatra, an area which is situated in the southwestern sector of the tropical Indo-Pacific Warm Pool. The high-resolution data sets (X-ray fluorescence, total organic carbon, and C37 alkenones) were used to track changes in paleoproductivity, freshwater budget, and sea surface temperature (SST) of the tropical climate system at orbital time scales over the past 300 ka. Our paleoclimatic data show that enhanced marine paleoproductivity was directly related to strengthening of coastal upwelling during periods of increased boreal summer insolation and associated SE monsoon strength with a precessional cyclicity. Changes in freshwater supply were primarily forced by precession-controlled changes in boreal NW winter monsoon rainfall enclosing an additional sea level component. SST variations of 2°-5°C occurred at eccentricity and precessional cyclicity. We suggest that the sea surface temperature variability off southern Sumatra is predominantly related to three major causes: (1) variations in upwelling intensity; (2) an elevated freshwater input into the southern Makassar Strait leading to reduced supply of warmer surface waters from the western Pacific and increased subsurface water transport via the Indonesian Throughflow into the Indian Ocean; and (3) long-term changes in the intensity or frequency of low-latitude climate phenomena, such as El Niño-Southern Oscillation.

Total organic carbon (TOC) on ground bulk and carbonate-free sediment samples of ODP Hole 151-911A

A multiproxy analysis of Hole 911A (Ocean Drilling Program (ODP) Leg 151) drilled on the Yermak Plateau (eastern Arctic Ocean) is used to investigate the behaviour of the Svalbard/Barents Sea ice sheet (SBIS) during late Pliocene and early Pleistocene (~3.0-1.7 Ma) climate changes. Contemporary with the 'Mid-Pliocene (~3 Ma) global warmth' (MPGW), a warmer period lasting ~300 kyr with seasonally ice-free conditions in the marginal eastern Arctic Ocean is assumed to be an important regional moisture source, and possibly one decisive trigger for intensification of the Northern Hemisphere glaciation in the Svalbard/Barents Sea area at ~2.7 Ma. An abrupt pulse of ice-rafted debris (IRD) to the Yermak Plateau at ~2.7 Ma reflects distinct melting of sediment-laden icebergs derived from the SBIS and may indicate the protruding advance of the ice sheet onto the outer shelf. Spectral analysis of the total organic carbon (TOC) record being predominantly of terrigenous/fossil-reworked origin indicates SBIS and possibly Scandinavian Ice Sheet response to incoming solar radiation at obliquity and precession periodicities. The strong variance in frequencies near the 41 kyr obliquity cycle between 2.7 and 1.7 Ma indicates, for the first time in the Arctic Ocean, a close relationship of SBIS growth and decay patterns to the Earth's orbital obliquity amplitudes, which dominated global ice volume variations during late Pliocene/early Pleistocene climate changes.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

(Table 1) Elemental mean values of the upper Paleocene and lower Eocene in IODP Hole 302-M0004A

We reconstruct the latest Paleocene and early Eocene (~57-50 Ma) environmental trends in the Arctic Ocean and focus on the Paleocene-Eocene thermal maximum (PETM) (~55 Ma), using strata recovered from the Lomonosov Ridge by the Integrated Ocean Drilling Program Expedition 302. The Lomonosov Ridge was still partially subaerial during the latest Paleocene and earliest Eocene and gradually subsided during the early Eocene. Organic dinoflagellate cyst (dinocyst) assemblages point to brackish and productive surface waters throughout the latest Paleocene and early Eocene. Dinocyst assemblages are cosmopolitan during this time interval, suggesting warm conditions, which is corroborated by TEX86'-reconstructed temperatures of 15°-18°C. Inorganic geochemistry generally reflects reducing conditions within the sediment and euxinic conditions during the upper lower Eocene. Spectral analysis reveals that the cyclicity, recorded in X-ray fluorescence scanning Fe data from close to Eocene thermal maximum 2 (~53 Ma, presence confirmed by dinocyst stratigraphy), is related to precession. Within the lower part of the PETM, proxy records indicate enhanced weathering, runoff, anoxia, and productivity along with sea level rise. On the basis of total organic carbon content and variations in sediment accumulation rates, excess organic carbon burial in the Arctic Ocean appears to have contributed significantly to the sequestration of injected carbon during the PETM.

(Table 1) Age determination of sediment core MD01-2404

We analyzed the high-resolution foraminifer isotope records, total organic carbon (TOC), and opal content from an Okinawa Trough core MD012404 in order to estimate the monsoon hydrography and productivity changes in the East China Sea (ECS) of the tropical western Pacific over the past 100,000 years. The variability shown in the records on orbital time scales indicates that high TOC intervals coincide with the increases of boreal May-September insolation driven by precession cycles (~21 ka), implying a strong connection to the variations in monsoons. We also observed possibly nearly synchronous, millennial-scale changes of the ECS surface hydrography (mainly driven by salinity changes but also by temperature effects) and productivity coincident with monsoon events in the Hulu/Dongge stalagmite isotope records. We found that increased freshening and high productivity correlate with high monsoon intensity in interstadials. This study suggests that the millennial-scale changes in monsoon hydrography and productivity in the ECS are remarkable and persistent features over the past 100,000 years.

Cyclostratigraphy and eccentricity tuning of the early Oligocene through early Miocene record from Walvis Ridge

Few astronomically calibrated high-resolution (<=5 kyr) climate records exist that span the Oligocene-Miocene time interval. Notably, available proxy records show responses varying in amplitude at frequencies related to astronomical forcing, and the main pacemakers of global change on astronomical time-scales remain debated. Here we present newly generated X-ray fluorescence core scanning and benthic foraminiferal stable oxygen and carbon isotope records from Ocean Drilling Program Site 1264 (Walvis Ridge, southeastern Atlantic Ocean). Complemented by data from nearby Site 1265, the Site 1264 benthic stable isotope records span a continuous ~13-Myr interval of the Oligo-Miocene (30.1-17.1 Ma) at high resolution (~3.0 kyr). Spectral analyses in the stratigraphic depth domain indicate that the largest amplitude variability of all proxy records is associated with periods of ~3.4 m and ~0.9 m, which correspond to 405- and ~110-kyr eccentricity, using a magnetobiostratigraphic age model. Maxima in CaCO3 content, d18O and d13C are interpreted to coincide with ~110 kyr eccentricity minima. The strong expression of these cycles in combination with the weakness of the precession- and obliquity-related signals allow construction of an astronomical age model that is solely based on tuning the CaCO3 content to the nominal (La2011_ecc3L) eccentricity solution. Very long-period eccentricity maxima (~2.4-Myr) are marked by recurrent episodes of high-amplitude ~110-kyr d18O cycles at Walvis Ridge, indicating greater sensitivity of the climate/cryosphere system to short eccentricity modulation of climatic precession. In contrast, the responses of the global (high-latitude) climate system, cryosphere, and carbon cycle to the 405-kyr cycle, as expressed in benthic d18O and especially d13C signals, are more pronounced during ~2.4-Myr minima. The relationship between the recurrent episodes of high-amplitude ~110-kyr d18O cycles and the ~1.2-Myr amplitude modulation of obliquity is not consistent through the Oligo-Miocene. Identification of these recurrent episodes at Walvis Ridge, and their pacing by the ~2.4-Myr eccentricity cycle, revises the current understanding of the main climate events of the Oligo-Miocene.

Magnetic susceptibility dataset from the Danian Contessa Highway section, Gubbio, Italy

We studied a high-resolution multiproxy data set, including magnetic susceptibility (MS), CaCO3 content, and stable isotopes (d18O and d13C), from the stratigraphic interval covering the uppermost Maastrichtian and the lower Danian, represented by the pelagic limestones of the Scaglia Rossa Formation continuously exposed in the classic sections of the Bottaccione Gorge and the Contessa Highway near Gubbio, Italy. Variations in all the proxy series are periodic and reflect astronomically forced climate changes (i.e., Milankovitch cycles). In particular, the MS proxy reflects variations in the terrigenous dust input in this pelagic, deep-marine environment. We speculate that the dust is mainly eolian in origin and that the availability and transport of dust are influenced by variations in the vegetation cover on the Maastrichtian-Paleocene African or Asian zone, which were respectively located at tropical to subtropical latitudes to the south or far to the east of the western Tethyan Umbria-Marche Basin, and were characterized by monsoonal circulation. The dynamics of monsoonal circulation are known to be strongly dependent on precession-driven and obliquity-driven changes in insolation. We propose that a threshold mechanism in the vegetation coverage may explain eccentricity-related periodicities in the terrigenous eolian dust input. Other mechanisms, both oceanic and terrestrial, that depend on the precession amplitude modulated by eccentricity, can be evoked together with the variation of dust influx in the western Tethys to explain the detected eccentricity periodicity in the d13C record. Our interpretations of the d18O and MS records suggest a warming event ~400 k.y. prior to the Cretaceous-Paleogene (K-Pg) boundary, and a period of climatic and environmental instability in the earliest Danian. Based on these multiproxy phase relationships, we propose an astronomical tuning for these sections; this leads us to an estimate of the timing and duration of several late Maastrichtian and Danian biostratigraphic and magnetostratigraphic events.

Age determination of North Atlantic sediment cores

Available overwash records from coastal barrier systems document significant variability in North Atlantic hurricane activity during the late Holocene. The same climate forcings that may have controlled cyclone activity over this interval (e.g., the West African Monsoon, El Niño-Southern Oscillation (ENSO)) show abrupt changes around 6000 yrs B.P., but most coastal sedimentary records do not span this time period. Establishing longer records is essential for understanding mid-Holocene patterns of storminess and their climatic drivers, which will lead to better forecasting of how climate change over the next century may affect tropical cyclone frequency and intensity. Storms are thought to be an important mechanism for transporting coarse sediment from shallow carbonate platforms to the deep-sea, and bank-edge sediments may offer an unexplored archive of long-term hurricane activity. Here, we develop this new approach, reconstructing more than 7000 years of North Atlantic hurricane variability using coarse-grained deposits in sediment cores from the leeward margin of the Great Bahama Bank. High energy event layers within the resulting archive are (1) broadly correlated throughout an offbank transect of multi-cores, (2) closely matched with historic hurricane events, and (3) synchronous with previous intervals of heightened North Atlantic hurricane activity in overwash reconstructions from Puerto Rico and elsewhere in the Bahamas. Lower storm frequency prior to 4400 yrs B.P. in our records suggests that precession and increased NH summer insolation may have greatly limited hurricane potential intensity, outweighing weakened ENSO and a stronger West African Monsoon-factors thought to be favorable for hurricane development.

Magnetic suceptibility, carbon and oxygen stable isotope ratios and reflectance from the Bottaccione Gorge and Furlo section, Italy

The oceans at the time of the Cenomanian-Turonian transition were abruptly perturbed by a period of bottom-water anoxia. This led to the brief but widespread deposition of black organic-rich shales, such as the Livello Bonarelli in the Umbria-Marche Basin (Italy). Despite intensive studies, the origin and exact timing of this event are still debated. In this study, we assess leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world, by providing a 6-Myr-long astronomically-tuned timescale across the Cenomanian-Turonian boundary. We procure insights in the relationship between orbital forcing and the Late Cretaceous carbon cycle by deciphering the imprint of astronomical cycles on lithologic, geophysical, and stable isotope records, obtained from the Bottaccione, Contessa and Furlo sections in the Umbria-Marche Basin. The deposition of black shales and cherts, as well as the onset of oceanic anoxia, is related to maxima in the 405-kyr cycle of eccentricity-modulated precession. Correlation to radioisotopic ages from the Western Interior (USA) provides unprecedented age control for the studied Italian successions. The most likely tuned age for the Livello Bonarelli base is 94.17 ± 0.15 Ma (tuning #1); however, a 405-kyr older age cannot be excluded (tuning #2) due to uncertainties in stratigraphic correlation, radioisotopic dating, and orbital configuration. Our cyclostratigraphic framework suggests that the exact timing of major carbon cycle perturbations during the Cretaceous may be linked to increased variability in seasonality (i.e. a 405-kyr eccentricity maximum) after the prolonged avoidance of seasonal extremes (i.e. a 2.4-Myr eccentricity minimum). Volcanism is probably the ultimate driver of oceanic anoxia, but orbital periodicities determine the exact timing of carbon cycle perturbations in the Late Cretaceous. This unites two leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world.

Age model and sedimentation rate of core M174/KL11 (Table 1)

Stable isotope and faunal records from the central Red Sea show high-amplitude oscillations for the past 380,000 years. Positive delta18O anomalies indicate periods of significant salt buildup during periods of lowered sea level when water mass exchange with the Arabian Sea was reduced due to a reduced geometry of the Bab el Mandeb Strait. Salinities as high as 53 per mil and 55 per mil are inferred from pteropod and benthic foraminifera delta18O, respectively, for the last glacial maximum. During this period all planktonic foraminifera vanished from this part of the Red Sea. Environmental conditions improved rapidly after 13 ka as salinities decreased due to rising sea level. The foraminiferal fauna started to reappear and was fully reestablished between 9 ka and 8 ka. Spectral analysis of the planktonic delta18O record documents highest variance in the orbital eccentricity, obliquity, and precession bands, indicating a dominant influence of climatically - driven sea level change on environmental conditions in the Red Sea. Variance in the precession band is enhanced compared to the global mean marine climate record (SPECMAP), suggesting an additional influence of the Indian monsoon system on Red Sea climates.

Alkenone accumulation rates, stable carbon isotope record, and sea surface temperature reconstruction for IODP Site 306-U1313

Here we present orbitally-resolved records of terrestrial higher plant leaf wax input to the North Atlantic over the last 3.5 Ma, based on the accumulation of long-chain n-alkanes and n-alkanl-1-ols at IODP Site U1313. These lipids are a major component of dust, even in remote ocean areas, and have a predominantly aeolian origin in distal marine sediments. Our results demonstrate that around 2.7 million years ago (Ma), coinciding with the intensification of the Northern Hemisphere glaciation (NHG), the aeolian input of terrestrial material to the North Atlantic increased drastically. Since then, during every glacial the aeolian input of higher plant material was up to 30 times higher than during interglacials. The close correspondence between aeolian input to the North Atlantic and other dust records indicates a globally uniform response of dust sources to Quaternary climate variability, although the amplitude of variation differs among areas. We argue that the increased aeolian input at Site U1313 during glacials is predominantly related to the episodic appearance of continental ice sheets in North America and the associated strengthening of glaciogenic dust sources. Evolutional spectral analyses of the n-alkane records were therefore used to determine the dominant astronomical forcing in North American ice sheet advances. These results demonstrate that during the early Pleistocene North American ice sheet dynamics responded predominantly to variations in obliquity (41 ka), which argues against previous suggestions of precession-related variations in Northern Hemisphere ice sheets during the early Pleistocene.

Geochemical and stable isotope composition of Middle Campanian marl-limestone rhythmites of the Lehrte West Syncline near Hannover (Lower Saxony, northern Germany)

A cyclic marl-limestone succession of Middle-Late Campanian age has been investigated with respect to a Milankovitch-controlled origin of geochemical data. In general, the major element geochemistry of the marl-limestone rhythmites can be explained by a simple two-component mixing model with the end-members calcium carbonate and 'average shale'-like material. Carbonate content varies from 55 to 90%. Non-carbonate components are clay minerals (illite, smectite) and biogenic silica from sponge spicules, as well as authigenically formed zeolites (strontian heulandite) and quartz. The redox potential suggests oxidizing conditions throughout the section. Trace element and stable isotopic data as well as SEM investigations show that the carbonate mud is mostly composed of low-magnesium calcitic tests of planktic coccolithophorids and calcareous dinoflagellate cysts (calcispheres). Diagenetic overprint results in a decrease of 2% d18O and an increase in Mn of up to 250 ppm. However, the sediment seems to preserve most of its high Sr content compared to the primary low-magnesium calcite of co-occurring belemnite rostra. The periodicity of geochemical cycles is dominated by 413 ka and weak signals between 51 and 22.5 ka, attributable to orbital forcing. Accumulation rates within these cycles vary between 40 and 50 m/Ma. The resulting cyclic sedimentary sequence is the product of (a) changes in primary production of low-magnesium calcitic biogenic material in surface waters within the long eccentricity and the precession, demonstrated by the CaCO3 content and the Mg/Al, Mn/Al and Sr/Al ratios, and (b) fluctuations in climate and continental weathering, which changed the quality of supplied clay minerals (the illite/smectite ratio), demonstrated by the K/Al ratio. High carbonate productivity correlates with smectite-favouring weathering (semi-arid conditions, conspicuously dry and moist seasonal changes in warmer climates). Ti as the proxy indicator for the detrital terrigenous influx, as well as Rb, Si, Zr and Na, shows only low frequency signals, indicating nearly constant rates of supply throughout the more or less pure pelagic carbonate deposition of the long-lasting third-order Middle-Upper Campanian sedimentary cycle.