Parameters of the Earths orbit for the last 5 Million years in 1 kyr resolution

New values for the astronomical parameters of the Earth's orbit and rotation (eccentricity, obliquity and precession) are proposed for paleoclimatic research related to the Late Miocene, the Pliocene and the Quaternary. They have been obtained from a numerical solution of the Lagrangian system of the planetary point masses and from an analytical solution of the Poisson equations of the Earth-Moon system. The analytical expansion developed in this paper allows the direct determination of the main frequencies with their phase and amplitude. Numerical and analytical comparisons with the former astronomical solution BER78 are performed so that the accuracy and the interval of time over which the new solution is valid can be estimated. The corresponding insolation values have also been computed and compared to the former ones. This analysis leads to the conclusion that the new values are expected to be reliable over the last 5 Ma in the time domain and at least over the last 10 Ma in the frequency domain.

Stable hydrogen and carbon isotope records and concentration data for four long chain n-alkanes from core GIK16160-3 derived from the Zambezi delta in the Mozambique Channel during the last 37,000 years

Hydrogen isotope values (dD) of sedimentary terrestrial leaf wax such as n-alkanes or n-acids have been used to map and understand past changes in rainfall amount in the tropics because dD of precipitation is commonly assumed as the first order controlling factor of leaf wax dD. Plant functional types and their photosynthetic pathways can also affect leaf wax dD but these biological effects are rarely taken into account in paleo studies relying on this rainfall proxy. To investigate how biological effects may influence dD values we here present a 37,000-year old record of dD and stable carbon isotopes (d13C) measured on four n-alkanes (n-C27, n-C29, n-C31, n-C33) from a marine sediment core collected off the Zambezi River mouth. Our paleo d13C records suggest that each individual n-alkanes had different C3/C4 proportional contributions. n-C29 was mostly derived from a C3 dicots (trees, shrubs and forbs) dominant vegetation throughout the entire record. In contrast, the longer chain n-C33 and n-C31 were mostly contributed by C4 grasses during the Glacial period but shifted to a mixture of C4 grasses and C3 dicots during the Holocene. Strong correlations between dD and d13C values of n-C33 (correlation coefficient R2 = 0.75, n = 58) and n-C31 (R2 = 0.48, n = 58) suggest that their dD values were strongly influenced by changes in the relative contributions of C3/C4 plant types in contrast to n-C29 (R2 = 0.07, n = 58). Within regions with variable C3/C4 input, we conclude that dD values of n-C29 are the most reliable and unbiased indicator for past changes in rainfall, and that dD and d13C values of n-C31 and n-C33 are sensitive to C3/C4 vegetation changes. Our results demonstrate that a robust interpretation of palaeohydrological data using n-alkane dD requires additional knowledge of regional vegetation changes from which nalkanes are synthesized, and that the combination of dD and d13C values of multiple n-alkanes can help to differentiate biological effects from those related to the hydrological cycle.

The record of stable isotope composition of surface and thermocline planktic foraminifera from sediment core 64PE-174P13 for the last 460,000 years

The upper branch of the Atlantic Meridional Overturning Circulation predominantly enters the Atlantic Ocean through the southeast, where the subtropical gyre is exposed to the influence of the Agulhas leakage (AL). To understand how the transfer of Indian Ocean waters via the AL affected the upper water column of this region, we have generated new proxy records of planktic foraminifera from a core on the central Walvis Ridge, on the eastern flank of the South Atlantic Gyre (SAG). We analyzed the isotopic composition of subsurface dweller Globigerinoides ruber sensu lato, and thermocline Globorotalia truncatulinoides sinistral, spanning the last five Pleistocene glacial-interglacial (G-IG) cycles. The former displays a response to obliquity, suggesting connection with high latitude forcing, and a warming tendency during each glacial termination, in response to the interhemispheric seesaw. The d18O difference between the two species, interpreted as a proxy for upper ocean stratification, reveals a remarkably regular sawtooth pattern, bound to G-IG cyclicity. It rises from interglacials until glacial terminations, with fast subsequent decrease, appearing to promptly respond to deglacial peaks of AL. Stratification, however, bears a different structure during the last cycle, being minimal at Last Glacial Maximum, and peaking at Termination I. We suggest this to be the result of the intensified glacial wind field over the SAG and/or of the invasion of the South Atlantic thermocline by Glacial North Atlantic Intermediate Waters. The d13C time series of the two species have similar G-IG pattern, whereas their difference is higher during interglacials. We propose that this may be the result of the alternation of intermediate water masses in different circulation modes, and of a regionally more efficient biological pump at times of high pCO2.

Sea-surface reconstructions of the western Pacific Ocean during the last 5.3 million years

The late Neogene was a time of cryosphere development in the northern hemisphere. The present study was carried out to estimate the sea surface temperature (SST) change during this period based on the quantitative planktonic foraminiferal data of 8 DSDP sites in the western Pacific. Target factor analysis has been applied to the conventional transfer function approach to overcome the no-analog conditions caused by evolutionary faunal changes. By applying this technique through a combination of time-slice and time-series studies, the SST history of the last 5.3 Ma has been reconstructed for the low latitude western Pacific. Although the present data set is close to the statistical limits of factor analysis, the clear presence of sensible variations in individual SST time-series suggests the feasibility and reliability of this method in paleoceanographic studies. The estimated SST curves display the general trend of the temperature fluctuations and reveal three major cool periods in the late Neogene, i.e. the early Pliocene (4.7 3.5 Ma), the late Pliocene (3.1-2.7 Ma), and the latest Pliocene to early Pleistocene (2.2-1.0 Ma). Cool events are reflected in the increase of seasonality and meridional SST gradient in the subtropical area. The latest Pliocene to early Pleistocene cooling is most important in the late Neogene climatic evolution. It differs from the previous cool events in its irreversible, steplike change in SST, which established the glacial climate characteristic of the late Pleistocene. The winter and summer SST decreased by 3.3-5.4°C and 1.0 2.1C in the subtropics, by 0.9°C and 0.6C in the equatorial region, and showed little or no cooling in the tropics. Moreover, this cooling event occurred as a gradual SST decrease during 2.2 1.0 Ma at the warmer subtropical sites, while that at cooler subtropical site was an abrupt SST drop at 2.2 Ma. In contrast, equatorial and tropical western Pacific experienced only minor SST change in the entire late Neogene. In general, subtropics was much more sensitive to climatic forcing than tropics and the cooling events were most extensive in the cooler subtropics. The early Pliocene cool periods can be correlated to the Antarctic ice volume fluctuation, and the latest Pliocene early Pleistocene cooling reflects the climatic evolution during the cryosphere development of the northern hemisphere.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Revised age models and foraminiferal records obtained for the Last Interglacial period in six marine sediment cores

The dataset contains the revised age models and foraminiferal records obtained for the Last Interglacial period in six marine sediment cores: - the Southern Ocean core MD02-2488 (age model, sea surface temperatures, benthic d18O and d13C for the period 136-108 ka), - the North Atlantic core MD95-2042 (age model, planktic d18O, benthic d18O and d13C for the period 135-110 ka), - the North Atlantic core ODP 980 (age model, planktic d18O, sea surface temperatures, seawater d18O, benthic d18O and d13C, ice-rafted detritus for the period 135-110 ka), - the North Atlantic core CH69-K09 (age model, planktic d18O, sea surface temperatures, seawater d18O, benthic d18O and d13C, ice-rafted detritus for the period 135-110 ka), - the Norwegian Sea core MD95-2010 (age model, percentage of Neogloboquadrina pachyderma sinistral, sea surface temperatures, benthic d18O, ice-rafted detritus for the period 134-110 ka), - the Labrador Sea core EW9302-JPC2 (age model, percentage of Neogloboquadrina pachyderma sinistral, sea surface temperatures, benthic d18O for the period 134-110 ka).

Stable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from the Holocene and the Last Glacial Maximum

We measured oxygen isotopes and Mg/Ca ratios in the surface-dwelling planktonic foraminifer Globigerinoides ruber (white s.s.) and the thermocline dweller Pulleniatina obliquiloculata to investigate upper ocean spatial variability in the Indo-Pacific Warm Pool (IPWP). We focused on three critical time intervals: the Last Glacial Maximum (LGM; 18-21.5 ka), the early Holocene (8-9 ka), and the late Holocene (0-2 ka). Our records from 24 stations in the South China Sea, Timor Sea, Indonesian seas, and western Pacific indicate overall dry and cool conditions in the IPWP during the LGM with a low thermal gradient between surface and thermocline waters. During the early Holocene, sea surface temperatures increased by ~3°C over the entire region, indicating intensification of the IPWP. However, in the eastern Indian Ocean (Timor Sea), the thermocline gradually shoaled from the LGM to early Holocene, reflecting intensification of the subsurface Indonesian Throughflow (ITF). Increased surface salinity in the South China Sea during the Holocene appears related to northward displacement of the monsoonal rain belt over the Asian continent together with enhanced influx of saltier Pacific surface water through the Luzon Strait and freshwater export through the Java Sea. Opening of the freshwater portal through the Java Sea in the early Holocene led to a change in the vertical structure of the ITF from surface- to thermocline-dominated flow and to substantial freshening of Timor Sea thermocline waters.

Authigenic uranium, mass accumulation rates, benthic d13C and SST during the last interglacial period of ODP Site 177-1094

Southern Ocean sediments reveal a spike in authigenic uranium 127,000 years ago, within the last interglacial, reflecting decreased oxygenation of deep water by Antarctic Bottom Water (AABW). Unlike ice age reductions in AABW, the interglacial stagnation event appears decoupled from open ocean conditions and may have resulted from coastal freshening due to mass loss from the Antarctic ice sheet. AABW reduction coincided with increased North Atlantic Deep Water (NADW) formation, and the subsequent reinvigoration in AABW coincided with reduced NADW formation. Thus, alternation of deep water formation between the Antarctic and the North Atlantic, believed to characterize ice ages, apparently also occurs in warm climates.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) of sediment core SO188_342 from the Northern Bay of Bengal for the last 18 ka

The Indian Summer Monsoon (ISM) is a major global climatic phenomenon. Long-term precipitation proxy records of the ISM, however, are often fragmented and discontinuous, impeding an estimation of the magnitude of precipitation variability from the Last Glacial to the present. To improve our understanding of past ISM variability, we provide a continuous reconstructed record of precipitation and continental vegetation changes from the lower Ganges-Brahmaputra-Meghna catchment and the Indo-Burman ranges over the last 18,000 years (18 ka). The records derive from a marine sediment core from the northern Bay of Bengal (NBoB), and are complemented by numerical model results of spatial moisture transport and precipitation distribution over the Bengal region. The isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) are compared to results from an isotope-enabled general atmospheric circulation model (IsoCAM) for selected time slices (pre-industrial, mid-Holocene and Heinrich Stadial 1). Comparison of proxy and model results indicate that past changes in the dD of precipitation and plant waxes were mainly driven by the amount effect, and strongly influenced by ISM rainfall. Maximum precipitation is detected for the Early Holocene Climatic Optimum (EHCO; 10.5-6 ka BP), whereas minimum precipitation occurred during the Heinrich Stadial 1 (HS1; 16.9-15.4 ka BP). The IsoCAM model results support the hypothesis of a constant moisture source (i.e. the NBoB) throughout the study period. Relative to the pre-industrial period the model reconstructions show 20% more rain during the mid-Holocene (6 ka BP) and 20% less rain during the Heinrich Stadial 1 (HS1), respectively. A shift from C4-plant dominated ecosystems during the glacial to subsequent C3/C4-mixed ones during the interglacial took place. Vegetation changes were predominantly driven by precipitation variability, as evidenced by the significant correlation between the dD and d13C alkane records. When compared to other records across the ISM domain, precipitation and vegetation changes inferred from our records and the numerical model results provide evidence for a coherent regional variability of the ISM from the Last Glacial to the present.

Planktic foraminiferal Mg/Ca and oxygen isotopes, plus estimated d18Oseawater and ice volume corrected d18Oseawater for the past 35000 years of sediment core GeoB10069-3

Evidence from geologic archives suggests that there were large changes in the tropical hydrologic cycle associated with the two prominent northern hemisphere deglacial cooling events, Heinrich Stadial 1 (HS1; ~19 to 15 kyr BP; kyr BP = 1000 yr before present) and the Younger Dryas (~12.9 to 11.7 kyr BP). These hydrologic shifts have been alternatively attributed to high and low latitude origin. Here, we present a new record of hydrologic variability based on planktic foraminifera-derived d18O of seawater (d18Osw) estimates from a sediment core from the tropical Eastern Indian Ocean, and using 12 additional d18Osw records, construct a single record of the dominant mode of tropical Eastern Equatorial Pacific and Indo-Pacific Warm Pool (IPWP) hydrologic variability. We show that deglacial hydrologic shifts parallel variations in the reconstructed interhemispheric temperature gradient, suggesting a strong response to variations in the Atlantic Meridional Overturning Circulation and the attendant heat redistribution. A transient model simulation of the last deglaciation suggests that hydrologic changes, including a southward shift in the Intertropical Convergence Zone (ITCZ) which likely occurred during these northern hemisphere cold events, coupled with oceanic advection and mixing, resulted in increased salinity in the Indonesian region of the IPWP and the eastern tropical Pacific, which is recorded by the d18Osw proxy. Based on our observations and modeling results we suggest the interhemispheric temperature gradient directly controls the tropical hydrologic cycle on these time scales, which in turn mediates poleward atmospheric heat transport.

Changes in bottom water corrosiveness and carbonate ion concentrations in the sub-Antarctic Atlantic during the last glacial period

The glacial climate system transitioned rapidly between cold (stadial) and warm (interstadial) conditions in the Northern Hemisphere. This variability, referred to as Dansgaard-Oeschger variability, is widely believed to arise from perturbations of the Atlantic Meridional Overturning Circulation. Evidence for such changes during the longer Heinrich stadials has been identified, but direct evidence for overturning circulation changes during Dansgaard-Oeschger events has proven elusive. Here we reconstruct bottom water [CO3]2- variability from B/Ca ratios of benthic foraminifera and indicators of sedimentary dissolution, and use these reconstructions to infer the flow of northern-sourced deep water to the deep central sub-Antarctic Atlantic Ocean. We find that nearly every Dansgaard-Oeschger interstadial is accompanied by a rapid incursion of North Atlantic Deep Water into the deep South Atlantic. Based on these results and transient climate model simulations, we conclude that North Atlantic stadial-interstadial climate variability was associated with significant Atlantic overturning circulation changes that were rapidly transmitted across the Atlantic. However, by demonstrating the persistent role of Atlantic overturning circulation changes in past abrupt climate variability, our reconstructions of carbonate chemistry further indicate that the carbon cycle response to abrupt climate change was not a simple function of North Atlantic overturning.

Sea surface temperature estimation for the last 70 ka from IODP Site 306-U1313

During the six Heinrich Events of the last 70 ka episodic calving from the circum-Atlantic ice sheets released large numbers of icebergs into the North Atlantic. These icebergs and associated melt-water flux are hypothesized to have led to a shutdown of Atlantic Meridional Overturning Circulation (AMOC) and severe cooling in large parts of the Northern Hemisphere. However, due to the limited availability of high-resolution records the magnitude sea surface temperature (SST) changes related to the impact of Heinrich Events on the mid-latitude North Atlantic is poorly constrained. Here we present a record of UK37'-based SSTs derived from sediments of Integrated Ocean Drilling Project (IODP) Site U1313, located at the southern end of the ice-rafted debris (IRD)-belt in the mid-latitude North Atlantic (41°N). We demonstrate that all six Heinrich Events are associated with a rapid warming of surface waters by 2 to 4°C in a few thousand years. The presence of IRD leaves no doubt about the simultaneous timing and correlation between rapid surface water warming and Heinrich Events. We argue that this warming in the mid-latitude North Atlantic is related to a northward expansion of the subtropical gyre during Heinrich Events. As a wide-range of studies demonstrated that in the central IRD-belt Heinrich Events are associated with low SSTs, these results thus identify an anti-phased (seesaw) pattern in SSTs during Heinrich Events between the mid-latitude (warm) and northern North Atlantic (cold). This highlights the complex response of surface water characteristics in the North Atlantic to Heinrich Events that is poorly reproduced by fresh water hosing experiments and challenges the widely accepted view that within the IRD-belt of the North Atlantic Heinrich Events coincide with periods of low SSTs.