Middle to Late Quaternary sedimentation and rock-magnetic of ODP Site 105-646

We examine rock-magnetic, carbonate, and planktonic foraminiferal fluxes to identify climatically controlled changes of terrigenous and pelagic sedimentation at Ocean Drilling Program (ODP) Site 646 (the Labrador Sea). Terrigenous sediments are brought to the site principally by bottom currents. We use a rock-magnetic parameter sensitive to changes in magnetic mineral grain size, the ratio of anhysteretic susceptibility to low-field magnetic susceptibility (XARM/X), to monitor changes in bottom-current intensity over time, with large values of XARM/X (finer-grained magnetic minerals) indicating weaker bottom currents. A second rock-magnetic parameter, magnetic mineral accumulation rate (KaT) was used to indicate variations in terrigenous flux. Planktonic foraminiferal and carbonate accumulation rates (Pfar and CaC03ar) are used as indicators of pelagic flux. Absolute age assignments are based on correlation between the planktonic foraminiferal oxygen-isotope variations for Site 646 and the SPECMAP master oxygen-isotope curve. Cross-correlation analyses of the parameters that we studied with respect to the SPECMAP curve suggest that from oxygen-isotope stages 21 to 11, sedimentation rate, KaT, X, CaCO3ar, and Pfar were at their maximums, whereas XARM/X was at its minimum during peak interglacials (i.e., 0 k.y. lag time with respect to minimum ice volume). However, all parameters we examined lag behind minimum ice volume from stages 11 to 1, indicating a change in timing of both pelagic and terrigenous fluxes at approximately 400 k.y. BP. The negative correlation coefficient between XARM/X and the SPECMAP curve further suggest that finer-grained magnetic minerals are deposited during glacial periods, which probably reflects weaker bottom currents. The shift observed in the lag times of parameters examined with respect to the SPECMAP record is attributed to a change in significance of orbital parameters. Spectral results exhibit strong power in eccentricity (about 100 k.y.) throughout the record. Kap X, CaCO3flr, and Pfar show significant power in obliquity (about 41 k.y.), whereas XARM/X shows significant power at 73 k.y. from stages 21 to 11. The 73-k.y. period in XARM/X is near the difference tone of obliquity and eccentricity: 1/43-1/102 = 1/69. Kar and XARM/X show power only in eccentricity from stages 11 to 1. X and Pfar show significant power in precession (about 18 and 22 k.y.) whereas CaC03ar has power at 34 k.y, which could be a combination of precession and obliquity. The shift in power of orbital parameters may by attributed to the effect of the about 413-k.y. signal of eccentricity.

(Appendix 1) Coral analysis and isostatic rebound effects from different Holes of IODP Expedition 310

Past sea-level records provide invaluable information about the response of ice sheets to climate forcing. Some such records suggest that the last deglaciation was punctuated by a dramatic period of sea-level rise, of about 20 metres, in less than 500 years. Controversy about the amplitude and timing of this meltwater pulse (MWP-1A) has, however, led to uncertainty about the source of the melt water and its temporal and causal relationships with the abrupt climate changes of the deglaciation. Here we show that MWP-1A started no earlier than 14,650 years ago and ended before 14,310 years ago, making it coeval with the Bølling warming. Our results, based on corals drilled offshore from Tahiti during Integrated Ocean Drilling Project Expedition 310, reveal that the increase in sea level at Tahiti was between 12 and 22 metres, with a most probable value between 14 and 18 metres, establishing a significant meltwater contribution from the Southern Hemisphere. This implies that the rate of eustatic sea-level rise exceeded 40 millimetres per year during MWP-1A.

Sedimentation rates calculated on surface sediment samples from different site of the Atlantic and Pacific Oceans (Table 1)

Radiocarbon ages on CaCO3 from deep-sea cores offer constraints on the nature of the CaCO3 dissolution process. The idea is that the toll taken by dissolution on grains within the core top bioturbation zone should be in proportion to their time of residence in this zone. If so, dissolution would shift the mass distribution in favor of younger grains, thereby reducing the mean radiocarbon age for the grain ensemble. We have searched in vain for evidence supporting the existence of such an age reduction. Instead, we find that for water depths of more than 4 km in the tropical Pacific the radiocarbon age increases with the extent of dissolution. We can find no satisfactory steady state explanation and are forced to conclude that this increase must be the result of chemical erosion. The idea is that during the Holocene the rate of dissolution of CaCO3 has exceeded the rain rate of CaCO3. In this circumstance, bioturbation exhumes CaCO3 from the underlying glacial sediment and mixes it with CaCO3 raining from the sea surface.

Radiolarian-based bransfer function for the estimation of sea surface temperatures in the Southern Ocean (Atlantic Sector)

A new radiolarian-based transfer function for sea surface temperature (SST) estimations has been developed from 23 taxa and taxa groups in 53 surface sediment samples recovered between 35° and 72°S in the Atlantic sector of the Southern Ocean. For the selection of taxa and taxa groups ecological information from water column studies was considered. The transfer function allows the estimation of austral summer SST (December-March) ranging between -1 and 18°C with a standard error of estimate of 1.2°C. SST estimates from selected late Pleistocene squences were sucessfully compared with independend paleotemperature estimates derived from a diatom transfer function. This shows that radiolarians provide an excellent tool for paleotemperature reconstructions in Pleistocene sediments of the Southern Ocean.

(Table 1) Sea surface temperature reconstruction for eastern equatorial Pacific surface sediment samples

Significant uncertainties persist in the reconstruction of past sea surface temperatures in the eastern equatorial Pacific, especially regarding the amplitude of the glacial cooling and the details of the post-glacial warming. Here we present the first regional calibration of alkenone unsaturation in surface sediments versus mean annual sea surface temperatures (maSST). Based on 81 new and 48 previously published data points, it is shown that open ocean samples conform to established global regressions of Uk'37 versus maSST and that there is no systematic bias from seasonality in the production or export of alkenones, or from surface ocean nutrient concentrations or salinity. The flattening of the regression at the highest maSSTs is found to be statistically insignificant. For the near-coastal Peru upwelling zone between 11-15°S and 76-79°W, however, we corroborate earlier observations that Uk'37 SST estimates significantly over-estimate maSSTs at many sites. We posit that this is caused either by uncertainties in the determination of maSSTs in this highly dynamic environment, or by biasing of the alkenone paleothermometer toward El Niño events as postulated by Rein et al. (2005).

Last Interglacial synthesis of high-latitude temperature: temperature anomalies and associated errors for 4 time slices

The Last Interglacial (LIG, 129-116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain. Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG. We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

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

Database of the GLAMAP project

A uniform chronology for foraminifera-based sea surface temperature records has been established in more than 120 sediment cores obtained from the equatorial and eastern Atlantic up to the Arctic Ocean. The chronostratigraphy of the last 30,000 years is mainly based on published d18O records and 14C ages from accelerator mass spectrometry, converted into calendar-year ages. The high-precision age control provides the database necessary for the uniform reconstruction of the climate interval of the Last Glacial Maximum within the GLAMAP-2000 project.

Modes of sapropel formation in the eastern Mediterranean: some constraints based on pyrite properties

Pyrite formation within and directly below sapropels in the eastern Mediterranean was governed by the relative rates of sulphide production and Fe liberation and supply to the organic-rich layers. At times of relatively high [SO4]2- reduction, sulphide could diffuse downward from the sapropel and formed pyrite in underlying sediments. The sources of Fe for pyrite formation comprised detrital Fe and diagenetically liberated Fe(II) from sapropel-underlying sediments. In organic-rich sapropels, input of Fe from the water column via Fe sulphide formation in the water may have been important as well. Rapid pyrite formation at high saturation levels resulted in the formation of framboidal pyrite within the sapropels, whereas below the sapropels slow euhedral pyrite formation at low saturation levels occurred. d34S values of pyrite are -33 per mil to -50 per mil. Below the sapropels d34S is lower than within the sapropels, as a result of increased sulphide re-oxidation at times of relatively high sulphide production and concentration when sulphide could escape from the sediment. The percentage of initially formed sulphide that was re-oxidized was estimated from organic carbon fluxes and burial efficiencies in the sediment. It ranges from 34% to 80%, varying significantly between sapropels. Increased palaeoproductivity as well as enhanced preservation contributed to magnified accumulation of organic matter in sapropels.

(Table 1) Ratio of hydroxy-GDGT vs. the total core GDGT was calculated based on the detection of hydroxy-GDGT

Hydroxylated glycerol dialkyl glycerol tetraethers (hydroxy-GDGTs) were detected in marine sediments of diverse depositional regimes and ages. Mass spectrometric evidence, complemented by information gleaned from two-dimensional (2D) 1H-13C nuclear magnetic resonance (NMR) spectroscopy on minute quantities of target analyte isolated from marine sediment, allowed us to identify one major compound as a monohydroxy-GDGT with acyclic biphytanyl moieties (OH-GDGT-0). NMR spectroscopic and mass spectrometric data indicate the presence of a tertiary hydroxyl group suggesting the compounds are the tetraether analogues of the widespread hydroxylated archaeol derivatives that have received great attention in geochemical studies of the last two decades. Three other related compounds were assigned as acyclic dihydroxy-GDGT (2OH-GDGT-0) and monohydroxy-GDGT with one (OH-GDGT-1) and two cyclopentane rings (OH-GDGT-2). Based on the identification of hydroxy-GDGT core lipids, a group of previously reported unknown intact polar lipids (IPLs), including the ubiquitously distributed H341-GDGT (Lipp J. S. and Hinrichs K. -U. (2009) Structural diversity and fate of intact polar lipids in marine sediments. Geochim. Cosmochim. Acta 73, 6816-6833), and its analogues were tentatively identified as glycosidic hydroxy-GDGTs. In addition to marine sediments, we also detected hydroxy-GDGTs in a culture of Methanothermococcus thermolithotrophicus. Given the previous finding of the putative polar precursor H341-GDGT in the planktonic marine crenarchaeon Nitrosopumilus maritimus, these compounds are synthesized by representatives of both cren- and euryarchaeota. The ubiquitous distribution and apparent substantial abundance of hydroxy-GDGT core lipids in marine sediments (up to 8% of total isoprenoid core GDGTs) point to their potential as proxies.

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.

GDGT-based proxies of sediment core YD0903

We reconstruct the environmental evolution of the East China Sea in the past 14 kyr based on glycerol dialkyl glycerol tetraethers (GDGTs) in a sediment core from the subaqueous Yangtze River Delta. Two primary phases are recognized. Phase I (13.8-8 cal kyr BP) reflects a predominantly continental influence, showing distinctly higher concentrations of branched GDGTs (averaged 143 ng/g dry sediment weight, dsw) than isoprenoid GDGTs (averaged 36 ng/g dsw), high BIT index (branched vs. isoprenoid tetraethers) values (>0.78) and a fluctuating GDGT-0/crenarchaeol ratio (R0/5, varied from 0.52 to 3.81). Within this interval, temporal increases of terrestrial and marine influence are attributed to Younger Dryas (YD) (ca. 12.9-12.2 cal kyr BP) cold event and melt-water pulse (MWP) -1B (11.5-11.1 cal kyr BP), respectively. The prominent transition from 8 to 7.9 cal kyr BP shows a sharp decrease in BIT index value (<0.4) and increase in crenarchaeol, which marks the beginning of phase II. Afterwards, the proxies remain relatively constant, which indicates that phase II (7.9 cal kyr BP-present) is a shelf sedimentary environment with high stand of sea level. Overall, the BIT index in our record serves as a good marker for terrestrial influence at the site, and likely reflects the flooding history of the region. The TEX86 (TetraEther Index of tetraethers consisting of 86 carbons) proxy is not applicable in phase I because of an excess terrestrial influence; but it seems to be valid for revealing the annual SST in phase II (21.6±0.9°C, n=49). In contrast, the MBT'/CBT (Methylation of Branched Tetraethers and Cyclization of Branched Tetraethers) proxy appears to faithfully record the annual mean air temperature (MAT) (14.3±0.63°C, n=68) and presents an integrated signal over the middle and lower Yangtze River drainage basin.