High-resolution planktic stable isotope and trace element record of ODP Hole 161-976

Multi-decadal to centennial planktic d18O and Mg/Ca records were generated at ODP976 in the Alboran Sea. The site is in the flow path of Atlantic inflow waters entering the Mediterranean and captured North Atlantic signals through the surface inflow and the atmosphere. The records reveal similar climatic oscillations during the last two glacial-to-interglacial transitions, albeit with a different temporal pacing. Glacial termination 1 (T1) was marked by Heinrich event 1 (H1), post-H1 Bolling/Allerod (B/A) warming and Younger Dryas (YD) cooling. During T2 the H11 d18O anomaly was twice as high and lasted 30% longer than during H1. The post-H11 warming marked the start of MIS5e while the subsequent YD-style cooling occurred during early MIS5e. The post-H11 temperature increase at ODP976 matched the sudden Asian Monsoon Termination II at 129 ka BP. Extending the 230Th-dated speleothem timescale to ODP976 suggests glacial conditions in the Northeast Atlantic region were terminated abruptly and interglacial warmth was reached in less than a millennium. The early-MIS5e cooling and freshening at ODP976 coincided with similar changes at North Atlantic sites suggesting this was a basin-wide event. By analogy with T1 we argue that this was a YD-type event that was shifted into the early stages of the last interglacial period. This scenario is consistent with evidence from northern North Atlantic and Nordic Sea sites that the continuing disintegration of the large Saalian Stage (MIS6) ice sheet in Eurasia delayed the advection of warm North Atlantic waters and full-strength convective overturn until later stages of MIS5e.

Age model, stable isotopes, major element proportions and endmember unmixing results from four western tropical Atlantic cores

We investigate changes in the delivery and oceanic transport of Amazon sediments related to terrestrial climate variations over the last 250 ka. We present high-resolution geochemical records from four marine sediment cores located between 5 and 12° N along the northern South American margin. The Amazon River is the sole source of terrigenous material for sites at 5 and 9° N, while the core at 12° N receives a mixture of Amazon and Orinoco detrital particles. Using an endmember unmixing model, we estimated the relative proportions of Amazon Andean material ("%-Andes", at 5 and 9° N) and of Amazon material ("%-Amazon", at 12° N) within the terrigenous fraction. The %-Andes and %-Amazon records exhibit significant precessional variations over the last 250 ka that are more pronounced during interglacials in comparison to glacial periods. High %-Andes values observed during periods of high austral summer insolation reflect the increased delivery of suspended sediments by Andean tributaries and enhanced Amazonian precipitation, in agreement with western Amazonian speleothem records. Increased Amazonian rainfall reflects the intensification of the South American monsoon in response to enhanced land-ocean thermal gradient and moisture convergence. However, low %-Amazon values obtained at 12° N during the same periods seem to contradict the increased delivery of Amazon sediments. We propose that reorganizations in surface ocean currents modulate the northwestward transport of Amazon material. In agreement with published records, the seasonal North Brazil Current retroflection is intensified (or prolonged in duration) during cold substages of the last 250 ka (which correspond to intervals of high DJF or low JJA insolation) and deflects eastward the Amazon sediment and freshwater plume.

The Antarctic ice core chronology (AICC2012)

Abstract of Bazin et al. (2013): An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120-800 ka. In this framework, new measurements of d18Oatm over Marine Isotope Stage (MIS) 11-12 on EDC and a complete d18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from d18Oatm, dO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one. Abstract of Veres et al. (2013): The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120-800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on d15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores.

(Table 1) Age determination of stalagmites of the Marcello Arévalo Cave in Chile

Stalagmites are important palaeo-climatic archives since their chemical and isotopic signatures have the potential to record high-resolution changes in temperature and precipitation over thousands of years. We present three U/Th-dated records of stalagmites (MA1-MA3) in the superhumid southern Andes, Chile (53°S). They grew simultaneously during the last five thousand years (ka BP) in a cave that developed in schist and granodiorite. Major and trace elements as well as the C and O isotope compositions of the stalagmites were analysed at high spatial and temporal resolution as proxies for palaeo-temperature and palaeo-precipitation. Calibrations are based on data from five years of monitoring the climate and hydrology inside and outside the cave and on data from 100 years of regional weather station records. Water-insoluble elements such as Y and HREE in the stalagmites indicate the amount of incorporated siliciclastic detritus. Monitoring shows that the quantity of detritus is controlled by the drip water rate once a threshold level has been exceeded. In general, drip rate variations of the stalagmites depend on the amount of rainfall. However, different drip-water pathways above each drip location gave rise to individual drip rate levels. Only one of the three stalagmites (MA1) had sufficiently high drip rates to record detrital proxies over its complete length. Carbonate-compatible element contents (e.g. U, Sr, Mg), which were measured up to sub-annual resolution, document changes in meteoric precipitation and related drip-water dilution. In addition, these soluble elements are controlled by leaching during weathering of the host rock and soils depending on the pH of acidic pore waters in the peaty soils of the cave's catchment area. In general, higher rainfall resulted in a lower concentration of these elements and vice versa. The Mg/Ca record of stalagmite MA1 was calibrated against meteoric precipitation records for the last 100 years from two regional weather stations. Carbonate-compatible soluble elements show similar patterns in the three stalagmites with generally high values when drip rates and detrital tracers were low and vice versa. d13C and d18O values are highly correlated in each stalagmite suggesting a predominantly drip rate dependent kinetic control by evaporation and/or outgassing. Only C and O isotopes from stalagmite MA1 that received the highest drip rates show a good correlation between detrital proxy elements and carbonate-compatible elements. A temperature-related change in rainwater isotope values modified the MA1 record during the Little Ice Age (~0.7-0.1 ka BP) that was ~1.5 °C colder than today. The isotopic composition of the stalagmites MA2 and MA3 that formed at lower drip rates shows a poor correlation with stalagmite MA1 and all other chemical proxies of MA1. 'Hendy tests' indicate that the degassing-controlled isotope fractionation of MA2 and MA3 had already started at the cave roof, especially when drip rates were low. Changing pathways and residence times of the seepage water caused a non-climatically controlled isotope fractionation, which may be generally important in ventilated caves during phases of low drip rates. Our proxies indicate that the Neoglacial cold phases from ~3.5 to 2.5 and from ~0.7 to 0.1 ka BP were characterised by 30% lower precipitation compared with the Medieval Warm Period from 1.2 to 0.8 ka BP, which was extremely humid in this region.

Iberian Margin 420 kyr geochemical and color variations record

Here we report 420 kyr long records of sediment geochemical and color variations from the southwestern Iberian Margin. We synchronized the Iberian Margin sediment record to Antarctic ice cores and speleothem records on millennial time scales and investigated the phase responses relative to orbital forcing of multiple proxy records available from these cores. Iberian Margin sediments contain strong precession power. Sediment "redness" (a* and 570-560 nm) and the ratio of long-chain alcohols to n-alkanes (C26OH/(C26OH + C29)) are highly coherent and in-phase with precession. Redder layers and more oxidizing conditions (low alcohol ratio) occur near precession minima (summer insolation maxima). We suggest these proxies respond rapidly to low-latitude insolation forcing by wind-driven processes (e.g., dust transport, upwelling, precipitation). Most Iberian Margin sediment parameters lag obliquity maxima by 7-8 ka, indicating a consistent linear response to insolation forcing at obliquity frequencies driven mainly by high-latitude processes. Although the lengths of the time series are short (420 ka) for detecting 100 kyr eccentricity cycles, the phase relationships support those obtained by Shackleton []. Antarctic temperature and the Iberian Margin alcohol ratios (C26OH/(C26OH + C29)) lead eccentricity maxima by 6 kyr, with lower ratios (increased oxygenation) occurring at eccentricity maxima. CO2, CH4, and Iberian SST are nearly in phase with eccentricity, and minimum ice volume (as inferred from Pacific d18Oseawater) lags eccentricity maxima by 10 kyr. The phase relationships derived in this study continue to support a potential role of the Earth's carbon cycle in contributing to the 100 kyr cycle.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Paleoclimate record combinatio of speleothems and of sediment core GeoB3910-2

A substantial strengthening of the South American monsoon system (SAMS) during Heinrich Stadials (HS) points toward decreased cross-equatorial heat transport as the main driver of monsoonal hydroclimate variability at millennial time-scales. In order to better constrain the exact timing and internal structure of HS1 over tropical South America we assessed two precisely dated speleothem records from central-eastern and northeastern Brazil in combination with two marine records of terrestrial organic and inorganic matter input into the western equatorial Atlantic. During HS1 we recognize at least two events of widespread intensification of the SAMS across the entire region influenced by the South Atlantic Convergence Zone (SACZ) at 16.11-14.69 kyr BP and 18.1-16.66 kyr BP (labeled as HS1a and HS1c, respectively), separated by a dry excursion from 16.66-16.11 kyr BP (HS1b). In view of the spatial structure of precipitation anomalies, the widespread increase of monsoon precipitation over the SACZ domain was termed 'Mega-SACZ'.

Sedimentology and geochemistry of ODP Site 184-1144

We analyzed sediment from Ocean Drilling Program (ODP) Site 1144 in the northern South China Sea to examine the weathering response of SE Asia to the strengthening of the East Asian Monsoon (EAM) since 14 ka. Our high-resolution record highlights the decoupling between continental chemical weathering, physical erosion and summer monsoon intensity. Mass accumulation rates, Ti/Ca, K/Rb, hematite/goethite and 87Sr/86Sr show sharp excursions from 11 to 8 ka, peaking at 10 ka. Clay minerals show a shorter-lived response with a higher kaolinite/(illite + chlorite) ratio at 10.7-9.5 ka. However, not all proxies show a clear response to environmental changes. Magnetic susceptibility rises sharply between 12 and 11 ka. Grain-size becomes finer from 14 to 10 ka and then coarsens until ~7 ka, but is probably controlled by bottom current flow and sealevel. Sr and Nd isotopes show that material is dominantly eroded from Taiwan with a lesser flux from Luzon, while clay mineralogy suggests that the primary sources during the Early Holocene were reworked via the shelf in the Taiwan Strait, rather than directly from Taiwan. Erosion was enhanced during monsoon strengthening and caused reworking of chemically weathered Pleistocene sediment largely from the now flooded Taiwan Strait, which was transgressed by ~8 ka, cutting off supply to the deep-water slope. None of the proxies shows an erosional response lasting until ~6 ka, when speleothem oxygen isotope records indicate the start of monsoon weakening. Although more weathered sediments were deposited from 11 to 8 ka when the monsoon was strong these are reworked and represent more weathering during the last glacial maximum (LGM) when the summer monsoon was weaker but the shelves were exposed.

(Table 1) Age control points of ODP Site 160-968

The astronomical timescale of the Eastern Mediterranean Plio-Pleistocene builds on tuning of sapropel layers to Northern Hemisphere summer insolation maxima. A 3000-year precession lag has become instrumental in the tuning procedure as radiocarbon dating revealed that the midpoint of the youngest sapropel, S1, in the early Holocene occurred approximately 3000 years after the insolation maximum. The origin of the time lag remains elusive, however, because sapropels are generally linked to maximum African monsoon intensities and transient climate modeling results indicate an in-phase behavior of the African monsoon relative to precession forcing. Here we present new high-resolution records of bulk sediment geochemistry and benthic foraminiferal oxygen isotopes from ODP Site 968 in the Eastern Mediterranean. We show that the 3000-year precession time lag of the sapropel midpoints is consistent with (1) the global marine isotope chronology, (2) maximum (monsoonal) precipitation conditions in the Mediterranean region and China derived from radiometrically dated speleothem records, and (3) maximum atmospheric methane concentrations in Antarctica ice cores. We show that the time lag relates to the occurrence of precession-paced North Atlantic cold events, which systematically delayed the onset of strong boreal summer monsoon intensity. Our findings may also explain a non-stationary behavior of the African monsoon over the past 3 million years due to more frequent and intensive cold events in the Late Pleistocene.

Point clouds of measurements in the Dechen Cave near Iserlohn, Germany

Modeling natural phenomena from 3D information enhances our understanding of the environment. Dense 3D point clouds are increasingly used as highly detailed input datasets. In addition to the capturing techniques of point clouds with LiDAR, low-cost sensors have been released in the last few years providing access to new research fields and facilitating 3D data acquisition for a broader range of applications. This letter presents an analysis of different speleothem features using 3D point clouds acquired with the gaming device Microsoft® Kinect. We compare the Kinect sensor with terrestrial LiDAR reference measurements using the KinFu pipeline for capturing complete 3D objects (< 4m**3). The results demonstrate the suitability of the Kinect to capture flowstone walls and to derive morphometric parameters of cave features. Although the chosen capturing strategy (KinFu) reveals a high correlation (R2=0.92) of stalagmite morphometry along the vertical object axis, a systematic overestimation (22% for radii and 44% for volume) is found. The comparison of flowstone wall datasets predominantly shows low differences (mean of 1 mm with 7 mm standard deviation) of the order of the Kinect depth precision. For both objects the major differences occur at strongly varying and curved surface structures (e.g. with fine concave parts).

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