(Table S1) Age determination of sediment cores SK99_17 and MD76-131

Productivity in the Arabian Sea is one of the highest in the world. It is controlled by seasonally reversing monsoonal wind-driven upwelling of nutrient-rich deeper waters which fuel phytoplankton growth. The detailed history of upwelling-induced productivity in the eastern Arabian Sea is unknown. Here we present paleoproductivity records from a composite sediment core at the millennial scale during the last 80 kyr B.P. These records are based on relative abundance counts of planktonic foraminifera and organic carbon contents, which are shown to mainly vary in concert. The eastern Arabian Sea upwelling-induced productivity was higher in the glacial period than in the Holocene, but it fell repeatedly on millennial timescales. These productivity declines occurred during cold events in the North Atlantic region, with the most pronounced changes prevailing during the Heinrich events. Hence, seasonal monsoon winds that drive upwelling-induced productivity in the east were weak when the North Atlantic was cold. These weak winds resulted in stratification of the water column, comparable to today's Arabian Sea stratification in the intermonsoonal period. Combining the new eastern with published western Arabian Sea results shows that the entire biological factory was severely diminished during the North Atlantic Heinrich events, and the seasonal productivity change in the Arabian Sea monsoon system was reduced with year-round low productivity.

Pore water and solid-phase measurements on sediment cores from the Black Sea

The surface sediments in the Black Sea are underlain by extensive deposits of iron (Fe) oxide-rich lake sediments that were deposited prior to the inflow of marine Mediterranean Sea waters ca. 9000 years ago. The subsequent downward diffusion of marine sulfate into the methane-bearing lake sediments has led to a multitude of diagenetic reactions in the sulfate-methane transition zone (SMTZ), including anaerobic oxidation of methane (AOM) with sulfate. While the sedimentary cycles of sulfur (S), methane and Fe in the SMTZ have been extensively studied, relatively little is known about the diagenetic alterations of the sediment record occurring below the SMTZ. Here we combine detailed geochemical analyses of the sediment and pore water with multicomponent diagenetic modeling to study the diagenetic alterations below the SMTZ at two sites in the western Black Sea. We focus on the dynamics of Fe, S and phosphorus (P) and demonstrate that diagenesis has strongly overprinted the sedimentary burial records of these elements. Our results show that sulfate-mediated AOM substantially enhances the downward diffusive flux of sulfide into the deep limnic deposits. During this downward sulfidization, Fe oxides, Fe carbonates and Fe phosphates (e.g. vivianite) are converted to sulfide phases, leading to an enrichment in solid phase S and the release of phosphate to the pore water. Below the sulfidization front, high concentrations of dissolved ferrous Fe (Fe2+) lead to sequestration of downward diffusing phosphate as authigenic vivianite, resulting in a transient accumulation of total P directly below the sulfidization front. Our model results further demonstrate that downward migrating sulfide becomes partly re-oxidized to sulfate due to reactions with oxidized Fe minerals, fueling a cryptic S cycle and thus stimulating slow rates of sulfate-driven AOM (~ 1-100 pmol/cm**3/d) in the sulfate-depleted limnic deposits. However, this process is unlikely to explain the observed release of dissolved Fe2+ below the SMTZ. Instead, we suggest that besides organoclastic Fe oxide reduction, AOM coupled to the reduction of Fe oxides may also provide a possible mechanism for the high concentrations of Fe2+ in the pore water at depth. Our results reveal that methane plays a key role in the diagenetic alterations of Fe, S and P records in Black Sea sediments. The downward sulfidization into the limnic deposits is enhanced through sulfate-driven AOM with sulfate and AOM with Fe oxides may provide a deep source of dissolved Fe2+ that drives the sequestration of P in vivianite below the sulfidization front.

Age, 232Th, 230Th, 238U, CaCO3, organic carbon, lithogenic fraction and opal concentrations of 4 sediment cores from the southwest Pacific, with calculated 230Th-normalized mass flux

No clear scenario has yet been able to explain the full carbon drawdown that occurred during the Last Glacial Maximum (LGM); however, increased export production (EP) in the Subantarctic Zone (SAZ) of the Southern Ocean due to iron (Fe) fertilisation has been proposed to have provided a key mechanism affecting the air-sea partitioning of carbon. We chronicle changes in marine EP based on four sediment cores in Subtropical Waters (STW) and SAZ around New Zealand since the LGM. For the first time in this region, we present 230-Thorium normalised fluxes of biogenic opal, carbonate (CaCO3), excess Barium (xsBa), and organic Carbon (Corg). In STW and SAZ, these flux variations show that EP did not change markedly since the LGM. The only exception was a site in the SAZ close to the STF, where we suggest the STF shifted over the core site, driving increased EP. To understand why EP was mostly low and constant we investigated dust deposition changes by measuring lithogenic fluxes at the four sites. These data are coherent with an increased dust deposition in the southwest Pacific during the LGM. Additionally, we infer an increased lithogenic material discharge from erosion and glacier melts during the deglaciation, limited to the Campbell Plateau. Therefore, we propose that even though increased glacial dust deposition may have relieved Fe limitation within the SAZ, the availability of silicic acid (Si(OH)4) limited any resultant increase in carbon export during the LGM. Consequently, we infer low Si(OH)4 concentrations in the SAZ that have not significantly changed since the LGM. This result suggests that both Si(OH)4 and Fe co-limit EP in the SAZ around New Zealand, which would be consistent with modern process studies.

Bulk sediment parameter of three sediment cores from permafrost lowlands, Alaska Arctic Coastal Plain

Arctic lowland landscapes have been modified by thermokarst lake processes throughout the Holocene. Thermokarst lakes form as a result of ice-rich permafrost degradation and they may expand over time through thermal and mechanical shoreline erosion. We studied proximal and distal sedimentary records from a thermokarst lake located on the Arctic Coastal Plain of northern Alaska to reconstruct the impact of catchment dynamics and morphology on the lacustrine depositional environment and to quantify carbon accumulation in thermokarst lake sediments. Short cores were collected for analysis of pollen, sedimentological and geochemical proxies. Radiocarbon and Pb/Cs dating, as well as extrapolation of measured historic lake expansion rates, were applied to estimate a minimum lake age of ~ 1,400 calendar years BP. The pollen record is in agreement with the young lake age as it does not include evidence of the "alder high" that occurred in the region ~ 4.0 cal ka BP. The lake most likely initiated from a remnant pond in a drained thermokarst lake basin (DTLB) and deepened rapidly as evidenced by accumulation of laminated sediments. Increasing oxygenation of the water column as shown by higher Fe/Ti and Fe/S ratios in the sediment indicate shifts in ice regime with increasing water depth. More recently, the sediment source changed as the thermokarst lake expanded through lateral permafrost degradation, alternating from redeposited DTLB sediments, to increased amounts of sediment from eroding, older upland deposits, followed by a more balanced combination of both DTLB and upland sources. The characterizing shifts in sediment sources and depositional regimes in expanding thermokarst lakes were therefore archived in the thermokarst lake sedimentary record. This study also highlights the potential for Arctic lakes to recycle old carbon from thawing permafrost and thermokarst processes.

The permafrost carbon inventory on the Tibetan Plateau : a new evaluation using deep sediment cores

Acknowledgements We are grateful for Dr. Jens Strauss and the other two anonymous reviewers for their insightful comments on an earlier version of this MS, and appreciate members of the IBCAS Sampling Campaign Teams for their assistance in field investigation. This work was supported by the National Basic Research Program of China on Global Change (2014CB954001 and 2015CB954201), National Natural Science Foundation of China (31322011 and 41371213), and the Thousand Young Talents Program.

Orbitally tuned time scale and astronomical forcing in the middle Eocene to early Oligocene determined on ODP and IODP sediment cores

Deciphering the driving mechanisms of Earth system processes, including the climate dynamics expressed as paleoceanographic events, requires a complete, continuous, and high-resolution stratigraphy that is very accurately dated. In this study, we construct a robust astronomically calibrated age model for the middle Eocene to early Oligocene interval (31-43 Ma) in order to permit more detailed study of the exceptional climatic events that occurred during this time, including the Middle Eocene Climate Optimum and the Eocene/Oligocene transition. A goal of this effort is to accurately date the middle Eocene to early Oligocene composite section cored during the Pacific Equatorial Age Transect (PEAT, IODP Exp. 320/321). The stratigraphic framework for the new time scale is based on the identification of the stable long eccentricity cycle in published and new high-resolution records encompassing bulk and benthic stable isotope, calibrated XRF core scanning, and magnetostratigraphic data from ODP Sites 171B-1052, 189-1172, 199-1218, and 207-1260 as well as IODP Sites 320-U1333, and -U1334 spanning magnetic polarity Chrons C12n to C20n. Subsequently we applied orbital tuning of the records to the La2011 orbital solution. The resulting new time scale revises and refines the existing orbitally tuned age model and the Geomagnetic Polarity Time Scale from 31 to 43 Ma. Our newly defined absolute age for the Eocene/Oligocene boundary validates the astronomical tuned age of 33.89 Ma identified at the Massignano (Italy) global stratotype section and point. Our compilation of geochemical records of climate-controlled variability in sedimentation through the middle-to-late Eocene and early Oligocene demonstrates strong power in the eccentricity band that is readily tuned to the latest astronomical solution. Obliquity driven cyclicity is only apparent during very long eccentricity cycle minima around 35.5 Ma, 38.3 Ma and 40.1 Ma.

Clay mineralogy of sediment cores from the South-East Atlantic off Namibia

Independent proxies were assessed in two Late Quaternary sediment cores from the eastern South Atlantic to compare deep-water changes during the last 400 kyr. Two cores were recovered from beneath North Atlantic Deep Water (NADW) at approximately 3 000 m depth. Late Quaternary presence of NADW is indicated by the Cibicidoides wuellerstorfi assemblage on the Walvis Ridge (Core GeoB 1214) and the Bulimina alazanensis assemblage on the Namibian continental slope (Core GeoB 1710). The propagation of NADW is exclusively observed during interglacials, with maximum factor loadings in Stages 1, 5, 7, 9 and 11. These maxima are consistent with peaks in kaolinite/chlorite ratios and maxima of poorly crystalline smectite in the clay-mineral record. Kaolinite and poorly crystalline smectite are products of intense chemical weathering. They are injected into the NADW at low latitudes, north of the study area, and advected south. Chlorite, which is stable under cold weathering regimes, is a characteristic mineral of water masses of southern origin. During glacial stages, it is advected north with Southern Component Water (SCW). Above the NADW/SCW depths, kaolinite/chlorite ratios vary only slightly without a significant glacial-interglacial pattern, as measured in a core (GeoB 1712) from 1 000 m deep on the same profile of the Namibian continental slope off Walvis Bay.

(Table 2) Age determination of Irminger Basin sediment cores

High-resolution records of coarse lithic content and oxygen isotope have been obtained in a piston core from the Irminger Basin. The last glacial period is characterized by numerous periods of increased iceberg discharges originating partly from Iceland and corresponding to millennial-scale instabilities of the coastal ice sheets and ice shelves in the Nordic area. A comparison with midlatitude sediment cores shows that ice-rafted material corresponding to the Heinrich events was deposited synchronously from 40° to 60°N. There are thus two oscillating systems: every 5-10 kyr massive iceberg armadas are released from large continental ice caps, whereas more frequent instabilities of the coastal ice sheets in the high latitude regions occur every 1.2-3.8 kyr. At the time of the Heinrich events the synchroneity of the response from all the northern hemisphere ice sheets attests the existence of strong interactions between the two systems.

Palaeoenvironment records of six sediment cores collected in the Gulf of Carpentaria

The Gulf of Carpentaria is an epicontinental sea (maximum depth 70 m) between Australia and New Guinea, bordered to the east by Torres Strait (currently 12 m deep) and to the west by the Arafura Sill (53 m below present sea level). Throughout the Quaternary, during times of low sea-level, the Gulf was separated from the open waters of the Indian and Pacific Oceans, forming Lake Carpentaria, an isolation basin, perched above contemporaneous sea-level with outlet channels to the Arafura Sea. A preliminary interpretation is presented of the palaeoenvironments recorded in six sediment cores collected by the IMAGES program in the Gulf of Carpentaria. The longest core (approx. 15 m) spans the past 130 ka and includes a record of sea-level/lake-level changes, with particular complexity between 80 and 40 ka when sea-level repeatedly breached and withdrew from Gulf/Lake Carpentaria. Evidence from biotic remains (foraminifers, ostracods, pollen), sedimentology and geochemistry clearly identifies a final marine transgression at about 9.7 ka (radiocarbon years). Before this transgression, Lake Carpentaria was surrounded by grassland, was near full, and may have had a surface area approaching 600 km-300 km and a depth of about 15 m. The earlier rise in sea-level which accompanied the Marine Isotopic Stage 6/5 transgression at about 130 ka is constrained by sedimentological and biotic evidence and dated by optical- and thermoluminescence and amino acid racemisation methods.

Laminated sediment cores in the Bering Sea

During the last glacial termination, the upper North Pacific Ocean underwent dramatic and rapid changes in oxygenation that lead to the transient intensification of oxygen minimum zones (OMZs), recorded by the widespread occurrence of laminated sediments on circum-Pacific continental margins. We present a new laminated sediment record from the mid-depth (1100 m) northern Bering Sea margin that provides insight into these deglacial OMZ maxima with exceptional, decadal-scale detail. Combined ultrahigh-resolution micro-X-ray-fluorescence (micro-XRF) data and sediment facies analysis of laminae reveal an alternation between predominantly terrigenous and diatom-dominated opal sedimentation. The diatomaceous laminae are interpreted to represent spring/summer productivity events related to the retreating sea ice margin.We identified five laminated sections in the deglacial part of our site. Lamina counts were carried out on these sections and correlated with the Bølling-Allerød and Preboreal phases in the North Greenland Ice Core (NGRIP) oxygen isotope record, indicating an annual deposition of individual lamina couplets (varves). The observed rapid decadal intensifications of anoxia, in particular within the Bølling-Allerød, are tightly coupled to short-term warm events through increases in regional export production. This dependence of laminae formation on warmer temperatures is underlined by a correlation with published Bering Sea sea surface temperature records and d18O data of planktic foraminifera from the Gulf of Alaska. The rapidity of the observed changes strongly implies a close atmospheric teleconnection between North Pacific and North Atlantic regions.We suggest that concomitant increases in export production and subsequent remineralization of organic matter in the Bering Sea, in combination with oxygen-poor waters entering the Being Sea, drove down oxygen concentrations to values below 0.1ml/l and caused laminae preservation. Calculated benthic-planktic ventilation ages show no significant variations throughout the last deglaciation, indicating that changes in formation rates or differing sources of North Pacific mid-depth waters are not prime candidates for strengthening the OMZ at our site. The age models established by our correlation procedure allow for the determination of calendar age control points for the Bølling-Allerød and the Preboreal that are independent of the initial radiocarbon-based chronology. Resulting surface reservoir ages range within 730-990 yr during the Bølling-Allerød, 800-1100 yr in the Younger Dryas, and 765-775 yr for the Preboreal.

(Table 1) Age determination of sediment cores EW9302-JPC1 and EW9302-JPC2

Results from two deep sea cores from northeast of Newfoundland at 1251 and 2527 m water depth, respectively, indicate that during the time period from 160,000 to 10,000 years BP, ice rafting events in the Labrador Sea were accompanied by rapid variations in deep and surface water circulation. Twelve ice-rafting events occurred, each coinciding with high concentrations of detrital carbonate and oxygen isotopic depletion of both surface and bottom waters. Eleven of these can be correlated with the North Atlantic Heinrich events H1-H11. The remaining very conspicuous ice-rafting event took place early in MIS substage 5e, at a time when the planktic faunal assemblage suggests marked warming of the sea surface. In the shallower core, benthic d13C values rise from a minimum during the deglaciation to peak substage 5e values following the last ice-rafting event, indicating that the ventilation of intermediate depths was renewed after the deglaciation was complete and continued throughout substage 5e. The benthic foraminifera suggest that this well-ventilated water mass was comparable to the modern Labrador Sea Water (LSW). The benthic faunas suggest that a relatively warm intermediate water mass entered the SE Labrador Sea during Heinrich events. Generally low benthic d13C values indicate that this water mass was poorly ventilated and rich in inorganic nutrients. Isotope data and benthic faunal distributions indicate that North Atlantic Deep Water (NADW) formed in the Norwegian-Greenland Sea reached the SE Labrador Sea between the Heinrich events.

(Table 1) AMS 14C dates from sediment cores RC12-10 and GYRE97-6PV-20

Proxy records from two piston cores in the Gulf of Mexico (GOM) provide a detailed (50-100 year resolution) record of climate variability over the last 14,000 years. Long-term (millennial-scale) trends and changes are related to the transition from glacial to interglacial conditions and movement of the average position of the Intertropical Convergence Zone (ITCZ) related to orbital forcing. The d18O of the surface-dwelling planktic foraminifer Globigerinoides ruber show negative excursions between 14 and 10.2 ka (radiocarbon years) that reflect influx of meltwater into the western GOM during melting of the Laurentide Ice Sheet. The relative abundance of the planktic foraminifer Globigerinoides sacculifer is related to transport of Caribbean water into the GOM. Maximum transport of Caribbean surface waters and moisture into the GOM associated with a northward migration of the average position of the ITCZ occurs between about 6.5 and 4.5 ka. In addition, abundance variations of G. sacculifer show century-scale variability throughout most of the Holocene. The GOM record is consistent with records from other areas, suggesting that century-scale variability is a pervasive feature of Holocene climate. The frequency of several cycles in the climate records is similar to cycles identified in proxy records of solar variability, indicating that at least some of the century-scale climate variability during the Holocene is due to external (solar) forcing.

(Table 2) Radiocarbon content and 14C AMS ages of planktonic foraminifera from sediment cores Louis_1900 and Louis_2023

North American freshwater runoff records have been used to support the case that climate flickers were caused by shutdowns of the ocean thermohaline circulation (THC) resulting from reversals of meltwater discharges. Inconsistencies in the documentation of these meltwater switches, however, continue to fuel the debate on the cause/s of the oscillatory nature of the deglacial climate. New oxygen and carbon isotope records from the northern Gulf of Mexico depict in exceptional detail the succession of meltwater floods and pauses through the southern routing during the interval 16 to 8.9 ka (14C years BP; ka, kiloannum). The records underscore the bimodal role played by the Gulf of Mexico as a destination of meltwater discharges from the receding Laurentide Ice Sheet. The evidence indicates that the Gulf of Mexico acted as the principal source of superfloods at 13.4, 12.6, and 11.9 ka that reached the North Atlantic and contributed significantly to density stratification, disruption of ocean ventilation, and cold reversals. Gulf of Mexico lapsed into a "relief valve" position in post-Younger Dryas time, when meltwater discharges were rerouted south at 9.9, 9.7, 9.4, and 9.1 ka, thus temporarily interrupting North Atlantic-bound freshwater discharges from Lake Agassiz. The history of meltwater events in the Gulf of Mexico contradicts the model that meltwater flow via the eastern outlets into the North Atlantic disrupted the ocean THC, causing cooling, while diversions to the Gulf of Mexico via the Mississippi River enhanced THC and warming.