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.

Eocene sedimentary calcium carbonate contents and stable isotope composition of benthic foraminifera

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

Coral oxygen isotope and Sr/Ca data from the Northern Gulf of Aqaba (Red Sea)

The last interglacial period (about 125,000 years ago) is thought to have been at least as warm as the present climate (Kukla et al., 2002, doi:10.1006/qres.2001.2316). Owing to changes in the Earth's orbit around the Sun, it is thought that insolation in the Northern Hemisphere varied more strongly than today on seasonal timescales (Berger, 1987, doi:10.1175/1520-0469(1978)035<2362:LTVODI>2.0.CO;2), which would have led to corresponding changes in the seasonal temperature cycle (Montoya et al., 2000, doi:10.1175/1520-0442(2000)013<1057:CSFKBW>2.0.CO;2). Here we present seasonally resolved proxy records using corals from the northernmost Red Sea, which record climate during the last interglacial period, the late Holocene epoch and the present. We find an increased seasonality in the temperature recorded in the last interglacial coral. Today, climate in the northern Red Sea is sensitive to the North Atlantic Oscillation (Felis et al., 2000 doi:10.1029/1999PA000477; Rimbu et al., 2001, doi:10.1029/2001GL013083), a climate oscillation that strongly influences winter temperatures and precipitation in the North Atlantic region. From our coral records and simulations with a coupled atmosphere-ocean circulation model, we conclude that a tendency towards the high-index state of the North Atlantic Oscillation during the last interglacial period, which is consistent with European proxy records (Zagwijn, 1996, doi:10.1016/0277-3791(96)00011-X; Aalbersberg and Litt, 1998, doi:10.1002/(SICI)1099-1417(1998090)13:5<367::AID-JQS400>3.0.CO;2-I; Klotz et al., 2003, doi:10.1016/S0921-8181(02)00222-9), contributed to the larger amplitude of the seasonal cycle in the Middle East.

Occurrence and stratigraphy of Coccolithus pelagicus in ODP Hole 120-747A

A record based on counts of the relative abundance of the dominant calcareous nannofossil taxa Coccolithus pelagicus and Reticulofenestra spp. in sediments recovered from Ocean Drilling Program Hole 747A (Kerguelen Plateau, Southern Indian Ocean) is established in this paper. This record (17 m.y. long) virtually spans the entire Miocene. Broad, steplike variations in the abundance of C. pelagicus range between 0% and 96%. Based on these variations, five stratigraphic units characterized by high abundance in C. pelagicus are delineated. We suggest that these variations are caused by water-mass movements (such as the north/south shifting of a front). This pronounced signal is compared with paleoceanographic events revealed by isotopic (d18O and d13C) studies. The five defined units are tentatively correlated to well-known global isotopic events. In particular, Units A and D correlate respectively with the Oligocene/Miocene boundary glaciation and the middle Miocene cooling event. Time-series analysis indicates the presence of the three main periodic components of the eccentricity of the Earth's orbit. A 200-k.y. cycle is also present. The stratigraphic and paleoceanographic significance of this record is discussed.

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

Radioisotopes and He accumulation in deep sea sediments of the Ontong Java Plateau

The accumulation of extraterrestrial 3He, a tracer for interplanetary dust particles (IDPs), in sediments from the Ontong Java Plateau (OJP; western equatorial Pacific Ocean) has been shown previously to exhibit a regular cyclicity during the late Pleistocene, with a period of ~100 ka. Those results have been interpreted to reflect periodic variability in the global accretion of IDPs that, in turn, has been linked to changes in the inclination of Earth's orbit with respect to the invariable plane of the solar system. Here we show that the accumulation in OJP sediments of authigenic 230Th, produced by radioactive decay of 234U in seawater, exhibits a 100-ka cyclicity similar in phase and amplitude to that evident in the 3He record. We interpret the similar patterns of 230Th and 3He accumulation to reflect a common origin within the ocean-climate system. Comparing spatial and temporal patterns of sediment accumulation against regional patterns of biological productivity and against the well-established pattern of CaCO3 dissolution in the deep Pacific Ocean leads to the further conclusion that a common 100-ka cycle in accumulation of biogenic, authigenic and extraterrestrial constituents in OJP sediments reflects the influence of climate-related changes in sediment focusing, rather than changes in the rate of production or supply of sedimentary constituents.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Sedimentary records of ODP Site 198-1209 on Shatsky Rise

The accumulation of wind blown (eolian) dust in deep-sea sediments reflects the aridity/humidity conditions of the continental region supplying the dust, as well as the "gustiness" of the climate system. Detailed studies of Pleistocene glacial-interglacial dust fluxes suggest changes in accumulation rates corresponding to orbital variations in solar insolation (Milankovitch cycles). While the orbital cycles found in sedimentary archives of the Pleistocene are intricately related to glacial growth and decay, similar global orbital signals recognized in deep-sea sediments of early Paleogene age, the last major greenhouse interval ~65-45 million years ago, could not have been linked to the waxing and waning of large ice sheets. Thus orbital signals recorded in early Paleogene sediments must reflect some other climate response to changes in solar insolation. To explore the potential connection between orbital forcing and the climate processes that control dust accumulation, we generated a high-resolution dust record for ~58 Myr old sediments from Shatsky Rise (ODP Site 1209, paleolatitude ~15°N-20°N). The dust accumulation data provide the first evidence of a correlation between dust flux to the deep sea and orbital cyclicity during the early Paleogene, indicating dust supply responded to insolation forcing during the last major interval of greenhouse climate. Furthermore, the relative amplitude of the dust flux response during the early Paleogene greenhouse was comparable to that during icehouse climates. Thus, subtle variations in solar insolation driven by changes in Earth's orbit about the Sun may have had a similar impact on climate during intervals of overall warmth as they did during glacial-interglacial states.

Dating analysis from DSDP Hole 90-594

Date-32 is a fast and easily used computer program developed to date Quaternary deep-sea cores by associating variations in the earth's orbit with recurring oscillations in core properties, such as carbonate content or isotope composition. Starting with known top and bottom dates, distortions in the periodicities of the core properties due to varying sedimentation rates are realigned by fast Fourier analysis so as to maximise the spectral energy density at the orbital frequencies. This allows age interpolation to all parts of the core to an accuracy of 10 kyrs, or about 1.5% of the record duration for a typical Brunhes sequence. The influence of astronomical forcing is examined and the method is applied to provide preliminary dates in a high-resolution Brunhes record from DSDP Site 594 off southeastern New Zealand.

Stable carbon and oxygen isotope compositions of Eocene benthic foraminifera from ODP Site 207-1258

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

(Supplementary Table 1) Iceberg-rafted debris (IBRD) mass accumulation rates (MAR) of IODP site 318-U1361

The Pliocene and Early Pleistocene, between 5.3 and 0.8 million years ago, span a transition from a global climate state that was 2-3 °C warmer than present with limited ice sheets in the Northern Hemisphere to one that was characterized by continental-scale glaciations at both poles. Growth and decay of these ice sheets was paced by variations in the Earth's orbit around the Sun. However, the nature of the influence of orbital forcing on the ice sheets is unclear, particularly in light of the absence of a strong 20,000-year precession signal in geologic records of global ice volume and sea level. Here we present a record of the rate of accumulation of iceberg-rafted debris offshore from the East Antarctic ice sheet, adjacent to the Wilkes Subglacial Basin, between 4.3 and 2.2 million years ago. We infer that maximum iceberg debris accumulation is associated with the enhanced calving of icebergs during ice-sheet margin retreat. In the warmer part of the record, between 4.3 and 3.5 million years ago, spectral analyses show a dominant periodicity of about 40,000 years. Subsequently, the powers of the 100,000-year and 20,000-year signals strengthen. We suggest that, as the Southern Ocean cooled between 3.5 and 2.5 million years ago, the development of a perennial sea-ice field limited the oceanic forcing of the ice sheet. After this threshold was crossed, substantial retreat of the East Antarctic ice sheet occurred only during austral summer insolation maxima, as controlled by the precession cycle.

Pollen abundances, and temperature and precipitation reconstruction for the Eocene to Oligocene transition in northern high latitudes

A profound global climate shift took place at the Eocene-Oligocene transition (~33.5 million years ago) when Cretaceous/early Palaeogene greenhouse conditions gave way to icehouse conditions (Zachos et al., 2001, doi:10.1126/science.1059412; Coxall et al., 2005, doi:10.1038/nature03135; Lear et al., 2008, doi:10.1130/G24584A.1). During this interval, changes in the Earth's orbit and a long-term drop in atmospheric carbon dioxide concentrations (Pagani et al., 2005, doi:10.1126/science.1110063; Pearson and Palmer, 2000, doi:10.1038/35021000; DeConto and Pollard, 2003, doi:10.1038/nature01290) resulted in both the growth of Antarctic ice sheets to approximately their modern size (Coxall et al., 2005, doi:10.1038/nature03135; Lear et al., 2008, doi:10.1130/G24584A.1) and the appearance of Northern Hemisphere glacial ice (Eldrett et al., 2007, doi:10.1038/nature05591; Moran et al., 2006, doi:10.1038/nature04800). However, palaeoclimatic studies of this interval are contradictory: although some analyses indicate no major climatic changes (Kohn et al., 2004, doi:10.1130/G20442.1; Grimes et al., 2005, doi:10.1130/G21019.1), others imply cooler temperatures (Zanazzi et al., 2007, doi:10.1038/nature05551), increased seasonality (Ivany et al., 2000, doi:10.1038/35038044; Terry, 2001, doi:10.1016/S0031-0182(00)00248-0) and/or aridity (Ivany et al., 2000, doi:10.1038/35038044; Terry, 2001, doi:10.1016/S0031-0182(00)00248-0; Sheldon et al., 2002, doi:10.1086/342865; Dupont-Nivet et al., 2007, doi:10.1038/nature05516). Climatic conditions in high northern latitudes over this interval are particularly poorly known. Here we present northern high-latitude terrestrial climate estimates for the Eocene to Oligocene interval, based on bioclimatic analysis of terrestrially derived spore and pollen assemblages preserved in marine sediments from the Norwegian-Greenland Sea. Our data indicate a cooling of ~5 °C in cold-month (winter) mean temperatures to 0-2 °C, and a concomitant increased seasonality before the Oi-1 glaciation event. These data indicate that a cooling component is indeed incorporated in the d18O isotope shift across the Eocene-Oligocene transition. However, the relatively warm summer temperatures at that time mean that continental ice on East Greenland was probably restricted to alpine outlet glaciers.

Age model of sediment core V20-120 (Table 2)

Time series analyses of atmospheric and oceanic variables in a late Pleistocene record from the northwest Pacific show the complex relationship of the response of various segments of the climate system to changes in the earth's orbit. Most variance spectra of time series from this subarctic record contain frequency peaks with periods corresponding to at least one of the major orbital components of eccentricity, obliquity, or precession. Although the radiolarian faunal (water mass) assemblages have prominent spectral peaks with 41,000-year periods which are coherent with obliquity at this frequency, only the Transitional faunal assemblage contains variance focused at a frequency corresponding to the 100,000-year period of eccentricity. Three of these faunal time series also show variance concentrated at a frequency with a 20,000-year period. These three time series are not coherent at a 20,000-year frequency with either of the dominant spectral peaks of precession. They are coherent, however, with variations in the second harmonic of the obliquity cycle. Changes in obliquity apparently affect siliceous faunal abundances in the northwest Pacific region of this high-latitude site more than variations in eccentricity or precession. Maxima in the time series of quartz abundance occur during low values of eccentricity and high glacial ice volume. Because atmospheric winds serve as the major source of supply of quartz to the sediments at this site, these high quartz values reflect increased aridity at the source region. Except for short periods during interglacials, the characteristics of the surface waters in this region of the subarctic Pacific during much of the last 460,000 years were similar to those which exist today in the Sea of Okhotsk. The spectrum of winter sea surface temperature estimates, derived from siliceous microfaunal abundances, contains dominant frequency peaks at periods of 100,000, 41,000, and 23,000 years which are coherent with eccentricity, obliquity, and precession, respectively. Based upon the relationship of the Subarctic Front with the dominance of specific faunal asemblages, the front was positioned south of its present-day location throughout much of the late Pleistocene.