Alternating Southern and Northern Hemisphere climate response to astronomical forcing during the past 35 m . y .

Earth’s climate has undergone different intervals of gradual change as well as abrupt shifts between climate states. Here we aim to characterize the corresponding changes in climate response to astronomical forcing in the icehouse portion of the Cenozoic, from the latest Eocene to the present. As a tool, we use a 35-m.y.-long dObenthic record compiled from different high-resolution benthic isotope records spliced together (what we refer to as a megasplice).We analyze the climate response to astronomical forcing during four 800-k.y.long time windows. During the mid-Miocene Climatic Optimum (ca. 15.5 Ma), global climate variability was mainly dependent on Southern Hemisphere summer insolation, amplified by a dynamic Antarctic ice sheet; 2.5 m.y. later, relatively warm global climate states occurred during maxima in both Southern Hemisphere and Northern Hemisphere summer insolation. At that point, the Antarctic ice sheet grew too big to pulse on the beat of precession, and the Southern Hemisphere lost its overwhelming influence on the global climate state. Likewise, we juxtapose response regimes of the Miocene (ca. 19 Ma) and Oligocene (ca. 25.5 Ma) warming periods. Despite the similarity in dObenthic values and variability, we find different responses to precession forcing. While Miocene warmth occurs during summer insolation maxima in both hemispheres, Oligocene global warmth is consistently triggered when Earth reaches perihelion in the Northern Hemisphere summer. This pattern is in accordance with previously published paleoclimate modeling results, and suggests an amplifying role for Northern Hemisphere sea ice. INTRODUCTION Hays et al. (1976, p. 1121) conclusively identified astronomical insolation forcing as the so-called “pacemaker of the Ice Ages.” From then on, deciphering the imprint of astronomical climate forcing in geologic archives of past climates became an indispensable part of paleoclimatology and paleoceanography (Hinnov and Hilgen, 2012). Astronomical-scale climatic changes are recorded in d18O of foraminiferal calcite tests. The d18O of benthic foraminifer tests reflects a combination of temperature state and size of Earth’s cryosphere at the time of their formation. These characteristics motivated the construction of the Pliocene–Pleistocene dObenthic stack of Lisiecki and Raymo (2005; LR04, http://www.lorraine-lisiecki.com/stack. html), as well as several dObenthic compilations (Cramer et al., 2009; Friedrich et al., 2012; Zachos et al., 2001, 2008). A stack is an average of different, coeval dObenthic records, placed on a common age model, whereas a compilation is a collection of records, all in the time domain, but with their separate age models. In this study we apply a different approach to integrate paleoclimate information from globally distributed dObenthic records. We constructed a 35-m.y.-long dObenthic megasplice (Fig. 1C) of different high-resolution benthic isotope records (see the GSA Data Repository1 for details). In contrast to existing compilations and stacks that respectively mix and average time-equivalent data across sites, the megasplice comprises data from a single site at any point in time, allowing analyses that require a high-resolution record at astronomical resolution. Visualizing Different Astronomical Rhythms Using a Gaussian Process The megasplice provides a novel stratigraphic tool that is suitable as a global framework for orbital-scale correlation of paleoclimate records. Here we quantify and visualize the changing paleoclimate response to astronomical forcing 1 GSA Data Repository 2017107, detailed description of the construction of the megasplice, details and evaluation of the Gaussian process, and the R-code to generate Figures 2 and 3, is available online at http://www.geosociety.org/datarepository/2017/ or on request from editing@geosociety.org. The megasplice dataset is available on PANGAEA at doi:10.1594 /PANGAEA .869815. (i.e., different astronomical rhythms) throughout geologic time. The wavelet spectrograms in Figures 1A and 1B show a remarkable transition from a dObenthic signal dominated by eccentricity-modulated precession to an obliquity-driven signal after the mid-Miocene climatic transition (MMCT). In this paper, we focus on the dObenthic response to astronomical forcing in four specific intervals: the MMCT, the mid-Miocene climatic optimum (MMCO), and the early Miocene and late Oligocene warming intervals. We calculate dObenthic response surfaces in a three-dimensional space spanned by obliquity and eccentricity-modulated precession (Figs. 2 and 3) within 800-k.y.-long time windows (Figs. 1A, 1B). The underlying statistical model is a Gaussian process (see the Data Repository) that renders the dependence of one variable, here dObenthic, in a multidimensional space. Gaussian processes are routinely used for geospatial interpolation (i.e., kriging; Cressie, 1993). The methodology would be susceptible to circular reasoning if dObenthic records were densely tuned to obliquity and precession. To alleviate this issue, we preserved only one tie point between depth and time every ~100 k.y., so that the phase relationship between astronomical forcing and dObenthic can be evaluated independently. The four time windows discussed herein all utilize data from the equatorial Pacific (Ocean Drilling Program Sites 1137, 1138, and 1218), thus avoiding complications when comparing astronomical imprints across oceans or latitudes. The dObenthic response surfaces (Figs. 2 and 3) show heavy to light dObenthic values as a function of the astronomical configuration. Blue to red colors designate cooler to warmer global climates as well as higher to lower global ice volume. The steps between different colors are fixed to 0.02‰, so that variability can be compared among different time windows. The dObenthic gradients in Figures 2B, 2E, 3B, and 3E reveal the relative imprint of obliquity and precession. Figure 2E, for example, has a steep gradient along the y axis, indicating that the corresponding climate system was responding more strongly to eccentricity-modulated precession (e · sin ; e is eccentricity and is the longitude *E-mail: ddevleeschouwer@marum.de GEOLOGY, April 2017; v. 45; no. 4; p. 375–378 | Data Repository item 2017107 | doi:10.1130/G38663.1 | Published online 7 February 2017 © 2017 The Authors. Gold Open Access: This paper is published under the terms of the CC-BY license. 376 www.gsapubs.org | Volume 45 | Number 4 | GEOLOGY 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 Age (Ma) 5 4 3 2 1 0 B en th ic δ 18 Sit e 1 26 7 (B ell et al ., 2 01 4)


INTRODUCTION
) conclusively identified astronomical insolation forcing as the so-called "pacemaker of the Ice Ages."From then on, deciphering the imprint of astronomical climate forcing in geologic archives of past climates became an indispensable part of paleoclimatology and paleoceanography (Hinnov and Hilgen, 2012).Astronomical-scale climatic changes are recorded in d 18 O of foraminiferal calcite tests.The d 18 O of benthic foraminifer tests reflects a combination of temperature state and size of Earth's cryosphere at the time of their formation.These characteristics motivated the construction of the Pliocene-Pleistocene d 18 O benthic stack of Lisiecki and Raymo (2005; LR04, http://www.lorraine-lisiecki.com/stack.html), as well as several d 18 O benthic compilations (Cramer et al., 2009;Friedrich et al., 2012;Zachos et al., 2001Zachos et al., , 2008)).A stack is an average of different, coeval d 18 O benthic records, placed on a common age model, whereas a compilation is a collection of records, all in the time domain, but with their separate age models.In this study we apply a different approach to integrate paleoclimate information from globally distributed d 18 O benthic records.We constructed a 35-m.y.-long d 18 O benthic megasplice (Fig. 1C) of different high-resolution benthic isotope records (see the GSA Data Repository1 for details).In contrast to existing compilations and stacks that respectively mix and average time-equivalent data across sites, the megasplice comprises data from a single site at any point in time, allowing analyses that require a high-resolution record at astronomical resolution.

Visualizing Different Astronomical Rhythms Using a Gaussian Process
The megasplice provides a novel stratigraphic tool that is suitable as a global framework for orbital-scale correlation of paleoclimate records.
Here we quantify and visualize the changing paleoclimate response to astronomical forcing (i.e., different astronomical rhythms) throughout geologic time.The wavelet spectrograms in Figures 1A and 1B show a remarkable transition from a d 18 O benthic signal dominated by eccentricity-modulated precession to an obliquity-driven signal after the mid-Miocene climatic transition (MMCT).In this paper, we focus on the d 18 O benthic response to astronomical forcing in four specific intervals: the MMCT, the mid-Miocene climatic optimum (MMCO), and the early Miocene and late Oligocene warming intervals.We calculate d 18 O benthic response surfaces in a three-dimensional space spanned by obliquity and eccentricity-modulated precession (Figs. 2 and 3) within 800-k.y.-long time windows (Figs.1A, 1B).The underlying statistical model is a Gaussian process (see the Data Repository) that renders the dependence of one variable, here d 18 O benthic , in a multidimensional space.Gaussian processes are routinely used for geospatial interpolation (i.e., kriging; Cressie, 1993).
The methodology would be susceptible to circular reasoning if d 18 O benthic records were densely tuned to obliquity and precession.To alleviate this issue, we preserved only one tie point between depth and time every ~100 k.y., so that the phase relationship between astronomical forcing and d 18 O benthic can be evaluated independently.The four time windows discussed herein all utilize data from the equatorial Pacific (Ocean Drilling Program Sites 1137, 1138, and 1218), thus avoiding complications when comparing astronomical imprints across oceans or latitudes.
The d 18 O benthic response surfaces (Figs. 2 and  3) show heavy to light d 18 O benthic values as a function of the astronomical configuration.Blue to red colors designate cooler to warmer global climates as well as higher to lower global ice volume.The steps between different colors are fixed to 0.02‰, so that variability can be compared among different time windows.The d 18 O benthic gradients in Figures 2B, 2E, 3B, and 3E reveal the relative imprint of obliquity and precession.Figure 2E, for example, has a steep gradient along the y axis, indicating that the corresponding climate system was responding more strongly to eccentricity-modulated precession (e • sin ; e is eccentricity and is the longitude  of perihelion measured from the spring equinox) than to obliquity.However, detailed insight into how the climate system responded to changes in eccentricity and precession is provided in Figures 2C, 2F, 3C, and 3F.These figures should be read as polar plots, for which the azimuth represents and the radius represents e (the center represents e = 0).The month during which Earth reaches perihelion is indicated at the corresponding azimuth.Figure 2F thus represents a climate that is warmest when the Earth is in perihelion at January, under high eccentricity.We provide a well-known example in Figure DR3 in the Data Repository, in which we compare the astronomical rhythms of the 100 k.y.world (0.8-0 Ma) and the obliquity world (1.8-1.0Ma).

Contrasting Miocene Rhythms
The response of d 18 O benthic to astronomical forcing between 13.6 and 12.8 Ma is shown in Figures 2A-2C.This time window covers the first 800 k.y. after the MMCT (Fig. 1) and thus represents a climate state with an expanded East Antarctic ice sheet (Flower and Kennett, 1994;Wright et al., 1992).Figures 2D-2F show the d 18 O benthic response between 16 and 15.2 Ma, a time slice that corresponds to the culmination of the MMCO (Fig. 1).During this relatively warm period, Antarctic ice sheets were much smaller than today and presumably more dynamic (Fielding et al., 2011;Griener et al., 2015;Passchier et al., 2013;Warny et al., 2009).A large difference in East Antarctic ice volume thus exists between the two time windows depicted in Figure 2. Nevertheless, the d 18 O benthic response to obliquity is similar for both time windows, with isotopic values becoming lighter with increasing obliquity (Figs.2B and 2E).The response to climatic precession (e • sin ), however, differs markedly between the two time windows (Figs.2C, 2F).The MMCO is a period with a low global ice volume and the isotopic response to eccentricity-modulated precession is intense; the lightest d 18 O benthic values occur when eccentricity is high and Earth reaches perihelion in January (Fig. 2F).The heaviest isotope values occur at low eccentricity.In the MMCT time slice (13.6-12.8Ma), the isotopic response to precession is bimodal, with lighter d 18 O benthic values when perihelion is reached in either July or January (Fig. 2C).These significantly different responses to precession portray the effect of a highly dynamic continental ice sheet.During the MMCO, the Northern Hemisphere was free of continental ice sheets, while the Antarctic ice sheet was smaller but more dynamic than it is today.Antarctic ice cover and global climate, represented by d 18 O benthic , mainly depended on Southern Hemisphere summer insolation, with minimum Antarctic ice volume and peak warmth when Southern Hemisphere summer insolation was maximum, i.e., perihelion in January (Fig. 2F).After extensive ice growth over Antarctica ca.13.8 Ma, this ice sheet grew too big to pulse on the beat of precession, and thereby lost its overwhelming influence on the global climate state.At the same time, the influence of the Northern Hemisphere was growing.For example, through ice-albedo feedback mechanisms, which became more important in the cooling Arctic region (Darby, 2008;Krylov et al., 2008;Stein et al., 2016).The combined effect of a less dynamic Antarctic ice sheet and a dynamic Arctic cryosphere explains the bimodal precession response between 13.6 and 12.8 Ma (MMCT; Fig. 2C).

Different Warming Regimes During the Oligocene and Miocene
We compare here two intervals of long-term warming.The younger time slice (19.6-18.8Ma; Figs.3A-3C) corresponds to the warming that eventually culminated in the MMCO; the older time slice coincides with the latest Oligocene warming (Pälike et al., 2006b) (26.1-25.3 Ma;Figs. 3D-3F).Based on the absolute d 18 O benthic values and their variability, one might expect similar climate response dynamics.Broadly speaking, this is what the surface response plots show, with cooler global climate (heavy d 18 O benthic ) under low eccentricity and low obliquity (Figs.3B, 3C, 3E, 3F).Nonetheless, we observe an important difference in the response to precession between the Miocene and Oligocene warming.Whereas the Miocene time slice shows a symmetrical and bimodal d 18 O benthic response to precession (Fig. 3C), Oligocene global warmth consistently occurs when Earth reaches perihelion in May or June (Fig. 3F).The late Oligocene response is consistent with an amplification of a Northern Hemisphere summer insolation signal through a vigorous feedback mechanism.This is remarkable, considering that during the Oligocene, continental ice in the Arctic was probably limited to no more than a few isolated glaciers (Eldrett et al., 2007;Pekar et al., 2006).Arctic sea ice, however, may have played a central role in positive feedback mechanisms, through its albedo as well as through its effect on deep-water formation and ocean conveyor circulation.At the same time in the Southern Hemisphere, the Antarctic continent was still heavily glaciated.For the 26.1-25.3Ma window considered here, Pekar et al. (2006) estimated an ice volume equivalent to 75%-100% of the present-day East Antarctic ice sheet, too large to respond on a precessional time scale.Meanwhile, climate modeling results by Eldrett et al. (2009) for the Arctic show a strong decrease in Arctic winter temperatures across the Eocene-Oligocene transition, with the formation of thick largescale winter sea ice.Given this general stage of the Oligocene climate system, we suggest that the d 18 O benthic response to precession (Fig. 3F) reflects the accelerated melting of sea ice in spring and early summer, when insolation is high during those seasons.This regional response to astronomical forcing is then propagated to global climate through multiple positive feedback mechanisms (e.g., ice albedo, deep-water formation).This interpretation is reinforced by the fact that minimum d 18 O benthic occurs when the Earth reaches perihelion in May, i.e., the month during which sea-ice volume starts to decline.

CONCLUSIONS
This paper documents four different ways in which the global climate system responds to astronomical insolation forcing from the latest Eocene to the present.Shifts in key boundary conditions like continental distribution, pCO 2 , or the size and distribution of the cryosphere thus had a significant influence on climate sensitivity to astronomical forcing.The presence of a dynamic cryosphere in the Southern or Northern Hemisphere seems to exert the principal control on the response of global climate to astronomical forcing in the icehouse of the past 35 m.y.We report an alternation of the driving hemisphere from the Northern Hemisphere during the late Oligocene, to the Southern Hemisphere during the MMCO, and back to the Northern Hemisphere during the Quaternary.The intermediate time slices (early Miocene and MMCT) exhibit a bimodal response to precession and thus suggest a balanced influence of both hemispheres on global climate.

Figure 2 .
Figure 2. Benthic d 18 O response to astronomical forcing between 13.6 and 12.8 Ma (left panels: mid-Miocene climatic transition), and between 16 and 15.2 Ma (right panels: mid-Miocene Climatic Optimum).Blue to red colors designate heavy to light d 18 O benthic values, and thus indicate cooler to warmer global climates as well as higher to lower global ice volume.The e • sin plots represent a slice at e • cos while the e • sin • cos show the dependency of d 18 O benthic , averaged over the distribution of obliquity (e is eccentricity and the longitude of perihelion measured from the spring equinox; see text).NH-Northern Hemisphere; SH-Southern Hemisphere.

Figure 3 .
Figure 3. Benthic d 18 O response to astronomical forcing between 19.6 and 18.8 Ma (left panels, early Miocene warming), and between 26.1 and 25.3 Ma (right panels, late Oligocene warming).NH-Northern Hemisphere; SH-Southern Hemisphere; e is eccentricity and is the longitude of perihelion measured from the spring equinox; see text).