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.

Sea surface temperature reconstruction based on alkenones of sediment core M35003-4

Evidence for abrupt climate changes on millennial and shorter timescales is widespread in marine and terrestrial climate records (Dansgard et al., 1993, doi:10.1038/364218a0; Bond et al., 1993, doi:10.1038/365143a0; Charles et al., 1996, doi:10.1016/0012-821X(96)00083-0, Bard et al., 1997, doi:10.1038/385707a0). Rapid reorganization of ocean circulation is considered to exert some control over these changes (Broecker et al., 1985, doi:10.1038/315021a0), as are shifts in the concentrations of atmospheric greenhouse gases (Broecker, 1994, doi:10.1038/372421a0). The response of the climate system to these two influences is fundamentally different: slowing of thermohaline overturn in the North Atlantic Ocean is expected to decrease northward heat transport by the ocean and to induce warming of the tropical Atlantic (Crowley, 1992, doi:10.1029/92PA01058; Manabe and Stouffer, 1997, doi:10.1029/96PA03932), whereas atmospheric greenhouse forcing should cause roughly synchronous global temperature changes (Manabe et al., 1991, doi:10.1175/1520-0442(1991)004<0785:TROACO>2.0.CO;2). So these two mechanisms of climate change should be distinguishable by the timing of surface-water temperature variations relative to changes in deep-water circulation. Here we present a high-temporal-resolution record of sea surface temperatures from the western tropical North Atlantic Ocean which spans the past 29,000 years, derived from measurements of temperature-sensitive alkenone unsaturation in sedimentary organic matter. We find significant warming is documented for Heinrich event H1 (16,900-15,400 calendar years bp) and the Younger Dryas event (12,900-11,600 cal. yr bp), which were periods of intense cooling in the northern North Atlantic. Temperature changes in the tropical and high-latitude North Atlantic are out of phase, suggesting that the thermohaline circulation was the important trigger for these rapid climate changes.

(Table 1) Age determination of western North Atlantic sediment cores

Foraminiferal abundance, 14C ventilation ages, and stable isotope ratios in cores from high deposition rate locations in the western subtropical North Atlantic are used to infer changes in ocean and climate during the Younger Dryas (YD) and Last Glacial Maximum (LGM). The d18O of the surface dwelling planktonic foram Globigerinoides ruber records the present-day decrease in surface temperature (SST) of ~4°C from Gulf Stream waters to the northeastern Bermuda Rise. If during the LGM the modern d18O/salinity relationship was maintained, this SST contrast was reduced to 2°C. With LGM to interglacial d18O changes of at least 2.2 per mil, SSTs in the western subtropical gyre may have been as much as 5°C colder. Above ~2.3 km, glacial d13C was higher than today, consistent with nutrient-depleted (younger) bottom waters, as identified previously. Below that, d13C decreased continually to -0.5 per mil, about equal to the lowest LGM d13C in the North Pacific Ocean. Seven pairs of benthic and planktonic foraminiferal 14C dates from cores >2.5 km deep differ by 1100 ± 340 years, with a maximum apparent ventilation age of ~1500 years at 4250 m and at ~4700 m. Apparent ventilation ages are presently unavailable for the LGM < 2.5 km because of problems with reworking on the continental slope when sea level was low. Because LGM d13C is about the same in the deep North Atlantic and the deep North Pacific, and because the oldest apparent ventilation ages in the LGM North Atlantic are the same as the North Pacific today, it is possible that the same water mass, probably of southern origin, flowed deep within each basin during the LGM. Very early in the YD, dated here at 11.25 ± 0.25 (n = 10) conventional 14C kyr BP (equal to 12.9 calendar kyr BP), apparent ventilation ages <2.3 km water depth were about the same as North Atlantic Deep Water today. Below ~2.3 km, four YD pairs average 1030 ± 400 years. The oldest apparent ventilation age for the YD is 1600 years at 4250 m. This strong contrast in ventilation, which indicates a front between water masses of very different origin, is similar to glacial profiles of nutrient-like proxies. This suggests that the LGM and YD modes of ocean circulation were the same.

(Table S1) Age determination of sediment core RAPiD-21-03K

A decadal resolution time series of sea surface temperature (SST) spanning the last two millennia is reconstructed by combining a proxy record from a new sediment sequence with previously published data from core MD99-2275, north of Iceland. The alkenone based SST reconstruction is validated with historic observational data and compared to a new similar temporal resolution reconstruction obtained from sediment core RAPiD21-3K, in the subpolar North Atlantic. The two SST paleorecords show consistent multidecadal scale coolings throughout the interval and similar expressions during the contrasted climatic periods 'colloquially known' as the Medieval Climatic Anomaly (MCA) and Little Ice Age (LIA). In order to further understand the temporal and spatial SST variations and investigate the influence of natural forcings on the observed SST changes during the last millennium, we compare our time series to simulations using the Institut Pierre-Simon Laplace IPSLCM4-v2 climate model. This comparison highlights the potential importance of volcanism as a natural forcing driving coherent abrupt cooling events captured in the subpolar North Atlantic records.

(Table S1) Age determination of ODP Site 165-1002

A series of 14C measurements in Ocean Drilling Program cores from the tropical Cariaco Basin, which have been correlated to the annual-layer counted chronology for the Greenland Ice Sheet Project 2 (GISP2) ice core, provides a high-resolution calibration of the radiocarbon time scale back to 50,000 years before the present. Independent radiometric dating of events correlated to GISP2 suggests that the calibration is accurate. Reconstructed 14C activities varied substantially during the last glacial period, including sharp peaks synchronous with the Laschamp and Mono Lake geomagnetic field intensity minimal and cosmogenic nuclide peaks in ice cores and marine sediments. Simulations with a geochemical box model suggest that much of the variability can be explained by geomagnetically modulated changes in 14C production rate together with plausible changes in deep-ocean ventilation and the global carbon cycle during glaciation.

Sea-surface temperature reconstruction for sediment core GIK18287-3 in the tropical South China Sea

The timing and magnitude of sea-surface temperature (SST) changes in the tropical southern South China Sea (SCS) during the last 16,500 years have been reconstructed on a high-resolution, 14C-dated sediment core using three different foraminiferal transfer functions (SIMMAX28, RAM, FP-12E) and geochemical (Uk'37) SST estimates. In agreement with CLIMAP reconstructions, both the FP-12E and the Uk'37 SST estimates show an average late glacial-interglacial SST difference of 2.0°C, whereas the RAM and SIMMAX28 foraminiferal transfer functions show only a minor (0.6°C) or no consistent late glacial-interglacial SST change, respectively. Both the Uk'37 and the FP-12E SST estimates, as well as the planktonic foraminiferal delta18O values, indicate an abrupt warming (ca. 1°C in <200 yr) at the end of the last glaciation, synchronous (within dating uncertainties) with the Bølling transition as recorded in the Greenland Ice Sheet Project 2 (GISP2) ice core, whereas the RAM-derived deglacial SST increase appears to lag during this event by ca. 500 yr. The similarity in abruptness and timing of the warming associated with the Bølling transition in Greenland and the southern SCS suggest a true synchrony of the Northern Hemisphere warming at the end of the last glaciation. In contrast to the foraminiferal transfer function estimates that do not indicate any consistent cooling associated with the Younger Dryas (YD) climate event in the tropical SCS, the Uk'37 SST estimates show a cooling of ca. 0.2-0.6°C compared to the Bølling-Allerød period. These Uk'37 SST estimates from the southern SCS argue in favor of a Northern Hemisphere-wide, synchronous cooling during the YD period.

(Table 2) Radiocarbon age determination from the polar planktonic foraminifera Neogloboquadrina pachyderma s for sediment core HU90-013-029

An additional Heinrich ice-rafting event is identified between Heinrich events 5 and 6 in eight cores from the Labrador Sea and the northwest Atlantic Ocean. It is characterized by sediment rich in detrital carbonate (40% CaCO3) with high concentration of floating dropstones, high coarse-fraction (% > 150 µm) content, and has a sharp contact with the underlying but grades into the overlying hemipelagic sediment. It also shows lighter d18ONpl values, indicating freshening due to iceberg rafting and/or meltwater discharge. This event is correlated with Dansgaard-Oeschger event 14 and interpreted as an additional Heinrich event, H5a. The thickness of H5a in the Labrador Sea reaches up to 220 cm. This additional Heinrich event has also been reported in cores PS2644 and SO82-5 from the northern North Atlantic. With the recognition of H5a the temporal spacing between Heinrich events 1 to 6 becomes more uniform (~7 ka).