CPClimate of the PastCPClim. Past1814-9332Copernicus PublicationsGöttingen, Germany10.5194/cp-14-2071-2018Climate evolution across the Mid-Brunhes TransitionClimate evolution across the Mid-Brunhes TransitionBarthAaron M.abarth@ehc.eduhttps://orcid.org/0000-0003-2042-6095ClarkPeter U.BillNicholas S.HeFenghttps://orcid.org/0000-0002-3355-6406PisiasNicklas G.Department of Geography and Earth and Environmental Sciences, Emory and Henry College, Emory, VA 24327, USACollege of Earth, Ocean, and Atmospheric Sciences, Oregon State University, Corvallis, OR 97331, USADepartment of Geoscience, University of Wisconsin – Madison, Madison, WI 53706, USACenter for Climatic Research, Nelson Institute for Environmental Studies, University of Wisconsin – Madison, Madison, WI 53706, USAAaron M. Barth (abarth@ehc.edu)21December201814122071208711March201827March20185November201826November2018This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/This article is available from https://cp.copernicus.org/articles/14/2071/2018/cp-14-2071-2018.htmlThe full text article is available as a PDF file from https://cp.copernicus.org/articles/14/2071/2018/cp-14-2071-2018.pdf
The Mid-Brunhes Transition (MBT) began ∼ 430 ka with an increase in
the amplitude of the 100 kyr climate cycles of the past 800 000 years. The
MBT has been identified in ice-core records, which indicate interglaciations
became warmer with higher atmospheric CO2 levels after the MBT, and
benthic oxygen isotope (δ18O) records, which suggest that
post-MBT interglaciations had higher sea levels and warmer temperatures than
pre-MBT interglaciations. It remains unclear, however, whether the MBT was a
globally synchronous phenomenon that included other components of the climate
system. Here, we further characterize changes in the climate system across
the MBT through statistical analyses of ice-core and δ18O
records as well as sea-surface temperature, benthic carbon isotope, and dust
accumulation records. Our results demonstrate that the MBT was a global event
with a significant increase in climate variance in most components of the
climate system assessed here. However, our results indicate that the onset of
high-amplitude variability in temperature, atmospheric CO2, and sea
level at ∼430 ka was preceded by changes in the carbon cycle, ice
sheets, and monsoon strength during Marine Isotope Stage (MIS) 14 and MIS 13.
Introduction
The last 800 kyr of the Pleistocene epoch is characterized by the emergence
of dominant ∼100 kyr glacial–interglacial climate cycles (Pisias and
Moore, 1981; Imbrie et al., 1993; Raymo et al., 1997; Clark et al., 2006).
These climate cycles typically have long glacial periods punctuated by short
interglaciations. Since ∼430 ka (i.e., starting with Marine Isotope
Stage (MIS) 11), interglaciations have experienced warmer temperatures
(Jouzel et al., 2007) and higher concentrations of atmospheric CO2
(Lüthi et al., 2008) relative to earlier interglaciations of the last
800 kyr (Fig. 1). The transition to higher-amplitude interglaciations has
also been recognized in deep-sea records of δ18O measured in
benthic foraminifera (Lisiecki and Raymo, 2005) that identify lesser ice
volume and/or warmer deep-ocean temperatures (Fig. 1).
Jansen et al. (1986) originally described this change in amplitude of
interglaciations as a singular Mid-Brunhes Event, but Yin (2013) argued that
it is more appropriately considered as a transition between two distinct
climate states, thus referring to it as the Mid-Brunhes Transition (MBT).
The change from low-amplitude to high-amplitude 100 kyr variability at
∼430 ka occurs during an interval of reduced eccentricity and
corresponding precession (Fig. 1), but similar orbital forcing occurred at
times before and after the onset of the MBT with no comparable response,
suggesting that the MBT was an unforced change internal to the climate
system. Mechanisms proposed for the MBT include a latitudinal shift in the
position of the Southern Hemisphere westerlies that increased upwelling of
respired carbon in the post-MBT Southern Ocean (Kemp et
al., 2010) and a change in Antarctic Bottom Water (AABW) formation through
insolation-induced feedbacks on sea ice and surface water density (Yin,
2013). However, several questions remain. (1) How and when was the MBT
expressed in other components of the climate system? (2) Was the MBT a
global or regional transition? (3) Did components expressing a transition
change synchronously? Here, we address these questions by providing a
statistical characterization of changes occurring over the last 800 kyr as
recorded by a variety of paleoclimatic proxies with broad spatial coverage.
Orbital forcing and climate records for the last 800 kyr.
(a, b, c) Precession, obliquity, and eccentricity (Laskar et al.,
2004). (d) Deuterium-derived temperature from the European Project
for Ice Coring in Antarctica (EPICA) Dome C ice core in Antarctica (Jouzel et
al., 2007). (e) Atmospheric CO2 from EPICA Dome C
(EPICA-community-members, 2004; Lüthi et al., 2008).
(f) Atmospheric CH4 from EPICA Dome C
(EPICA-community-members, 2004). (g) Global benthic oxygen isotope
stack (Lisiecki and Raymo, 2005).
Site locations. Map indicating the locations of the cores used in
this research and modern sea-surface temperature values. Each symbol
represents a different proxy record. Diamonds indicate sea-surface temperatures.
Circles indicate benthic ∂13C. Squares indicate dust.
MethodsData collection
We compiled all available published records of sea-surface temperature
(SST), benthic marine carbon isotopes ratios (δ13C), and dust
accumulation (dust) that met our selection criteria and closely represented
a global distribution as attainable (Fig. 2). Each data set has an average
temporal resolution of <5 kyr, does not include any large age gaps,
and spans much or all of the entire time period of consideration to limit
biasing of the younger parts of the record. Lisiecki (2014) placed all of
the δ13C records on the LR04 age model. Published SST records
that were not on the LR04 age model were placed on it in one of two ways. If
the original data had depth and benthic δ18O data, the SST
record was placed on LR04 using the ager script in MATLAB as part of the
ARAND software package (Howell et al., 2006). When only benthic
δ18O records were available, the SST records were placed on
LR04 by selecting corresponding tie points in the δ18O data
series using the AnalySeries version 2.0 software (Paillard et al.,
1996). Because some dust records could not be placed on the LR04 age model,
certain statistical analyses of them (e.g., phase–lag relationships) are
likely not robust, but the overall variance in them is preserved. Each
record was then interpolated to a time step (Δt) of 2 kyr. With each
record having an average resolution <5 kyr, this Δt allows
for the preservation of higher-frequency variability while limiting the
number of interpolated data points.
We used empirical orthogonal function (EOF) analysis to characterize the
dominant modes of variability and robustly demonstrate global and regional
signals of the SST, δ13C, and dust records. We then used
spectral analyses of each resulting principal component (PC) to characterize
their periodicity, phase, and amplitude.
Sea-surface temperatures
We used 11 SST records that span the entire 800 kyr time period and four
additional records that span 8–758 ka. Inclusion of these four shorter
records does not change our conclusions. The SST records cover the Pacific
(n=9), Atlantic (n=5), and Indian (n=1) oceans (Fig. 2, Table 1).
We note that Shakun et al. (2015) reconstructed a global SST stack for the
last 800 kyr using 49 records, but only seven of these spanned the entire
800 kyr. Comparison of our SST PC1 based on 15 records to the Shakun SST
stack shows excellent agreement (Fig. S1 in the Supplement).
Carbon isotopes (δ13C)
We analyzed the global set of δ13C records compiled by Lisiecki (2014) (n=26; Fig. 2), and separately analyzed the records in the
Atlantic (n=14) and the Pacific (n=4) basins, thus distinguishing
between the dominant water masses within each basin and removing the muting
effect of the more negative Pacific values on the more positive Atlantic.
Similar to SSTs, Lisiecki (2014) reconstructed a global δ13C
stack for the last 3 Myr using 46 records, but only 18 of these spanned the
last 800 kyr. Comparison of our δ13C PC1 to the Lisiecki
δ13C stack shows excellent agreement (Fig. S2).
We then looked at regional and depth stacks of the δ13C records
in the Atlantic basin to characterize changes in the dominant water masses
on orbital timescales. Regional stacks were broken into North Atlantic
(>20∘ N; n=4), equatorial Atlantic
(20∘ S to 20∘ N; n=14), and South
Atlantic (>20∘ S; n=8). We also created
stacks for the deep North Atlantic (depth >2000 m; n=4) and
intermediate North Atlantic (depth <2000 m; n=3). All included
records were averaged to create the stack and each stacked record was
interpolated to a 2 kyr time step. Stacking improves the signal-to-noise
ratio of the δ13C records, making regional stacks useful in
identifying circulation changes and comparing circulation responses with
other climate records (Lisiecki, 2014).
Data compilation. All data sets are used in these analyses with
associated locations, proxy type, references, and digital object identifier
when available.
CoreLocationLatitudeLongitudeΔt (ka)ProxyDOI/URL (last access: 15 September 2018)ReferenceDSDP 607Northeastern Atlantic41.0012-32.95732.8SST – Transfer function10.1594/PANGAEA.701229Ruddiman et al. (1989)ODP 846Eastern equatorial Pacific-3.0949-90.8182.3SST – Uk3710.1126/science.1185435Herbert et al. (2010)ODP 982North Atlantic57.5165-15.86674.7SST – Uk37https://www.ncdc.noaa.gov/paleo/study/8624Lawrence et al. (2009)ODP 1143Western equatorial Pacific9.3619113.28511.9SST – Uk3710.1594/PANGAEA.786444Li et al. (2011)ODP 1082South Atlantic-21.094111.82054.5SST – Uk3710.1594/PANGAEA.786701Etourneau et al. (2009)ODP 1313North Atlantic41.0011-32.95731.4SST – Uk3710.1594/PANGAEA.744483Naafs et al. (2012)ODP 722Arabian Sea16.621859.79532.0SST – Uk3710.1126/science.1185435Herbert et al. (2010)ODP 1146South China Sea19.4567116.27271.6SST – Uk3710.1126/science.1185435Herbert et al. (2010)ODP 846Eastern equatorial Pacific-3.0949-90.8181.3SST – Uk3710.1038/nature02338Liu and Herbert (2004)MD97-2140Western Pacific warm pool2.1792141.45813.9SST – Mg / Cahttps://www.ncdc.noaa.gov/paleo/study/6266de Garidel-Thoron et al. (2005)MD06-3018Tropical western Pacific-22.5977166.86285.2SST – Mg / Cahttps://www.ncdc.noaa.gov/paleo/study/11188Russon et al. (2010)ODP 806Western equatorial Pacific0.319159.3612.4SST – Mg / Ca10.1594/PANGAEA.772015Medina-Elizalde and Lea (2005)DSDP 594Southwest Pacific-45.5235174.9482.6SST – Modern analog10.1594/PANGAEA.691478Schaefer et al. (2005)V22-174South Atlantic-10.0667-12.81674.1SST – Transfer function10.1594/PANGAEA.52228Mix et al. (1999)RC13-110Eastern equatorial Pacific0-964.9SST – Transfer functionPisias and Mix (1997)ODP 659Eastern equatorial Atlantic18.0772-21.02623.6Dust flux10.1594/PANGAEA.696121Tiedemann et al. (1994)ODP 1090Subantarctic Atlantic-42.91378.89970.3Dust MAR10.1594/PANGAEA.767460Martinez-Garcia et al. (2011)ODP 1090Subantarctic Atlantic-42.91378.89970.3Fe MAR10.1594/PANGAEA.767460Martinez-Garcia et al. (2011)CLPChinese Loess Plateau1.0Grain size10.1029/2006GC001287Sun and An (2005)Lake BaikalSouthern Russia0.5Silica %https://www.ncdc.noaa.gov/paleo/study/6068Prokopenko et al. (2006)PS75-074Pacific Southern Ocean-56.4696-142.99540.4Fe counts10.1594/PANGAEA.826600Lamy et al. (2014)PS75-076Pacific Southern Ocean-56.4696-142.99540.4Fe wt %10.1594/PANGAEA.826600Lamy et al. (2014)ODP 663Eastern equatorial Atlantic-1.1978-11.87852.7Terrestrial %10.1594/PANGAEA.208129deMenocal et al. (1993)EPICA Dome CAntarctica-75.06123.210.2Dust10.1594/PANGAEA.695995Lambert et al. (2008)ODP 982North Atlantic57.5-15.92.5∂13C10.1594/PANGAEA.700897Venz et al. (1999)ODP 983North Atlantic60.4-23.61.0∂13Chttps://www.ncdc.noaa.gov/paleo/study/2543McIntyre et al. (1999)ODP 984North Atlantic61-243.5∂13Chttps://www.ncdc.noaa.gov/paleo/study/5897Raymo et al. (2004)DSDP 607North Atlantic41-334.1∂13C10.1594/PANGAEA.52379Ruddiman et al. (1989)ODP 658North Atlantic20.8-18.71.6∂13C10.1594/PANGAEA.68570Tiedmann et al. (1994)U 1308North Atlantic49.9-24.20.3∂13Chttps://www.ncdc.noaa.gov/paleo/study/10250Hodell et al. (2008)ODP 980/981North Atlantic55.5-14.71.6∂13C10.1594/PANGAEA.698998Oppo et al. (1998)DSDP 502Equatorial Atlantic11.5-79.4∂13C10.1594/PANGAEA.701470deMenocal et al. (1992)ODP 664Equatorial Atlantic0.1-23.23.2∂13Chttps://www.ncdc.noaa.gov/paleo/study/2529Raymo et al. (1997)ODP 925Equatorial Atlantic4.2-43.54.3∂13CBickert et al. (1997)ODP 926Equatorial Atlantic3.7-42.92.7∂13CLisiecki et al. (2008)ODP 927Equatorial Atlantic5.5-44.5∂13CBickert et al. (1997)ODP 929Equatorial Atlantic5.5-44.54.9∂13CBickert et al. (1997)ODP 928Equatorial Atlantic5.5-44.82.5∂13CLisiecki et al. (2008)ODP 1090South Atlantic-42.98.92.8∂13C10.1594/PANGAEA.696106Venz and Hodell (2002)GeoB 1032South Atlantic-22.963.7∂13C10.1594/PANGAEA.54655Wefer et al. (1996)GeoB 1035South Atlantic-21.653.9∂13C10.1594/PANGAEA.58766Bickert and Wefer (1996)ODP 1089South Atlantic-47.99.90.4∂13C10.1594/PANGAEA.701432Hodell et al. (2003)GeoB 1211South Atlantic-24.57.54.9∂13C10.1594/PANGAEA.103634Bickert and Wefer (1996)GeoB 1214South Atlantic-24.77.24.5∂13C10.1594/PANGAEA.103635Bickert and Wefer (1996)
Continued.
CoreLocationLatitudeLongitudeΔt (ka)ProxyDOI/URL (last access: 15 September 2018)ReferenceRC13-229South Atlantic-25.511.33.8∂13C10.1594/PANGAEA.701361Oppo et al. (1990)TN 576South Atlantic-42.98.91.5∂13Chttps://www.ncdc.noaa.gov/paleo/study/2576Hodell et al. (2000)ODP 1143Pacific9.4-246.73.8∂13C10.1594/PANGAEA.784150Cheng et al. (2004)ODP 677Pacific4.2-83.72.8∂13C10.1594/PANGAEA.701316Shackleton et al. (1990)ODP 846Pacific-3.1-90.82.2∂13C10.1594/PANGAEA.808207Mix et al. (1995a)ODP 849Pacific0.2-110.53.5∂13C10.1594/PANGAEA.701400Mix et al. (1995b)ODP 1123Southwest Pacific-41.7862-171.4990.8Mg / Ca10.1594/PANGAEA.786205Elderfield et al. (2012)EPICA Dome CAntarctica-75.06123.210.4CH4https://www.ncdc.noaa.gov/paleo/study/6093Loulergue et al. (2008)EPICA Dome CAntarctica-75.06123.210.4CO2https://www.ncdc.noaa.gov/paleo/study/6091Lüthi et al. (2008)EPICA Dome CAntarctica-75.06123.213.0Deuteriumhttps://www.ncdc.noaa.gov/paleo/study/6080EPICA-community-members (2004)Dust
We analyzed eight proxy records of dust that span the entire 800 kyr time
period, and then separated them by hemisphere (Northern Hemisphere had three;
Southern Hemisphere had five) to characterize hemispheric differences (Fig. 2). The various proxies
for dust include Fe mass accumulation rates, weight percent of terrigenous
material and Fe, flux of lithogenic grains, and grain size analysis. We
standardized each record before analysis to account for these various proxy
types and their differing range in values, thus allowing for comparison of
their relative amplitudes of variation.
Empirical orthogonal function (EOF) analysis
We used EOF analysis to objectively characterize the climate variability
recorded by the proxies across the MBT. Analyses of covariance between the
data were conducted using the EOF script as part of the ARAND software
package (Howell et al., 2006). The results provide both the dominant
variability as a time series (principal component) and a spatial
distribution of variance contribution (factor loadings). The records for SST
and δ13C were kept in their original values of degrees and per
mil, respectively, to preserve the original variance. Dust records were
standardized to a mean value of zero and unit variance so that each record
provided equal weight to the EOF. Statistical significance of all EOFs was
determined through segmented linear regression analysis. All resulting break
points occur on or after the second EOF and are thus considered significant.
Spectral analysis
We used the Blackman–Tukey technique in the ARAND software package for
spectral analysis of each PC (Howell et al., 2006). Analyses were conducted
using all data points within the time interval of interest, boxcar windowing
of the input data, and hamming spectral filter. Multiple tests were conducted
for the 8–800, 450–800, and 8–350 ka time slices. These intervals
characterize the dominant frequency of variability over the entire 800 kyr
record, and for the pre- and post-MBT intervals, respectively. The removal of
the 350–450 ka interval limited the influence of MIS 11, MIS 12, and
Termination V (T5), as these were shown to potentially bias the spectral
power. Furthermore, these selected intervals result in time series of equal
length to limit biasing of longer records. Additional tests were conducted
using wavelet analyses that characterize the change in spectral power as a
time series. Complementary spectral analyses were conducted on CO2
and CH4 records from the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core
(EPICA-community-members, 2004; Jouzel et al., 2007), and benthic
δ18O using the LR04 stack (Lisiecki and Raymo, 2005).
Cross-spectral analyses were conducted for the PCs against mean insolation
values to determine phase and coherency of each PC. Mean insolation values were
calculated for each of the dominant periodicities (eccentricity, obliquity,
and precession) with the data derived from AnalySeries (Laskar et al., 2004;
Paillard et al., 1996).
Variance tests
We used f tests to test for variance changes across the MBT for each
principal component from the EOF analysis as well as for CO2, CH4,
and the LR04 δ18O records. Analyses were conducted in MATLAB
using the vartest2 script. This approach assumes the null hypothesis that the pre-
and post-MBT distributions of the time series of each climate component have
the same normally distributed variance. If the resulting variance values
reject this hypothesis of no statistical difference, then the pre- and
post-MBT time series are determined to have undergone a significant change
in variance across the MBT. We interpret the change in variance to reflect a
change in the amplitude of each climate signal.
ResultsCO2, CH4, and benthic δ18O
Time series of the greenhouse gases CO2 and CH4 and of the LR04
stack of benthic δ18O suggest an increase in their interglacial
values across the MBT (Fig. 1). Spectral analyses of the LR04 stack and
atmospheric CO2 indicate a small post-MBT increase in the 100 kyr band,
whereas results for CH4 indicate a decrease (Fig. S3). All three
records show an increase in the precessional band (19–23 kyr). Variance
tests suggest that δ18O and CO2 have a statistically
significant increase in variance across the MBT, while CH4 variance
decreases (Table S1).
Sea-surface temperatures
EOF analysis of global SSTs over the last 758 kyr identifies two
statistically significant principal components (Fig. 3a). The first and
second principal components (PC1 and PC2, respectively) account for 69 %
of the total variance, with PC1 explaining 49 % alone. While some degree of
regional variability in each record exists, factor loadings indicate that
each record positively contributed to PC1 with a larger contribution coming
from high-latitude records. Thus, PC1 is representative of a global SST
signal. SST PC1 demonstrates a stepwise increase in variance starting at 436 ka, with an increase of interglacial temperatures, while showing no
significant change in the lower limit glacial values, which is one of the
defining characteristics of the MBT. The highest spectral density is in the
100 kyr frequency band throughout the entire time period (Fig. S3d).
Wavelet analysis (Fig. 4a) shows a significant increase in the
100 kyr frequency band at 580 ka that reaches its maximum spectral power during
MIS 11 and persists throughout most of the remaining interval, albeit with
decreasing intensity after ∼250 ka. Variance f tests reveal a
significant increase in amplitude from the pre- to post-MBT SSTs (Table S1).
These results thus confirm that there was a stepwise global transition of
SSTs from lower- to higher-amplitude interglaciations as previously inferred
from individual records.
Principal components. Plots of the first (PC1; blue) and second
(PC2; red) principal components from our EOF analysis of each climate
variable. Percent variance is explained by each PC represented by the numbers
with the corresponding color. (a) Sea-surface temperatures.
(b) Dust records. (c) Global ∂13C.
(d)∂13C of the Atlantic.
(e)∂13C of the Pacific.
Wavelet analysis. Wavelets of four of the first principal
components. (a) Sea-surface temperatures. (b) Dust records.
(c) Global ∂13C.
(d)∂13C of the Atlantic. Red colors represent
higher spectral power. Blue colors represent lower spectral power.
Statistical significance is highlighted by the thin black line. Milankovitch
periods are highlighted by the dashed horizontal lines.
Variance calculations on proxies of bottom water temperature (Elderfield et al., 2012) and on the Antarctic EPICA ice-core
deuterium record (EPICA-community-members, 2004), a measure of Antarctic
atmospheric temperature, also indicate statistically significant increases
in variance across the MBT (Table S1). In both proxies, the time series
indicate an increase of interglacial temperature values while showing no
significant change to the lower limit glacial values, similar to PC1 of SSTs
(Fig. 5).
Dust
The EOF analysis of the global dust records identifies two statistically
significant principal components, with PC1 representing 56 % of the total
variance and PC2 15 % (Fig. 3b). All records but the one from the
Chinese Loess Plateau (CLP) reflect increased dust accumulation due to
increased aridity and/or wind strength during glaciations, whereas higher
dust accumulation in the CLP record reflects increased summer Asian monsoon
strength, which is an interglacial signal (Sun and An, 2005).
Accordingly, factor loadings for the dust records are all positive for PC1
except for the CLP.
Temperature records. (a) Deuterium-based temperature record
from EPICA Dome C in Antarctica (light yellow; Jouzel et al., 2007). The
darker yellow line is a 15-point moving average. (b) The first
principal component of our sea-surface temperature analysis (red).
(c) Bottom water temperature derived from Mg / Ca measurements
at ODP 1123 (light blue; Elderfield et al., 2012). The dark blue line is a
15-point moving average.
In contrast to the change in variance seen in temperature, CO2, and
CH4 during MIS 11, variance tests of the dust PC1 suggest a stepwise
increase in variance during MIS 12, with subsequent glaciations having higher
amplitudes (Table S1). Separating the records by hemisphere shows that the
increase in glacial amplitude starting at MIS 12 occurs in the southern PC1
but not in the northern PC1 (Fig. 6). Similarly, the signal during MIS 14
present in the global PC1 is absent in the northern PC1, suggesting that the
northern control on dust accumulation was skipped during that glacial.
Spectral analysis of the global PC1 indicates dominant power in the
100 kyr frequency band that increases in spectral power across the MBT
(Fig. S3b). Furthermore, wavelet analysis of PC1 demonstrates an increase
in the spectral power of the 100 kyr band at ∼600 ka with its
highest power during MIS 11 (Fig. 4b), similar to the SST PC1. The 100 kyr
frequency remains statistically significant throughout the 100–600 ka interval.
δ13C
The first principal component of the global δ13C (δ13CG; PC1) explains 58 % of the total variance (Fig. 3c). EOF
analysis of δ13C records from the Atlantic basin (δ13CATL) yields two statistically significant PCs, with PC1 and PC2
explaining 58 % and 13 % of the total variance, respectively (Fig. 3d). EOF analysis of δ13C records from the Pacific (δ13CPAC) yields one statistically significant principal component
(PC1 is 81 % total variance) (Fig. 4e).
Dust principal components. The first principal components of our
dust analysis for the global (yellow), north (red), and south (blue) records.
Vertical gray boxes highlight specific glacial (dark gray) and interglacial
(light gray) periods. The numbers indicate the associated MIS of each box.
Both the global and Atlantic PC1 exhibit a strong 100 kyr frequency that is
persistent from 680 to 180 ka (Fig. 4c, d). Unlike SST and dust, however,
δ13CG and δ13CATL
demonstrate a stronger 100 kyr power prior to MIS 11 with its highest power
throughout MIS 13 and 12 (510–460 ka). Spectral analysis shows a decrease
in power of the 100 kyr frequency band from pre- to post-MBT (Fig. S3f, g).
Variance tests show that the pre- and post-MBT intervals for
δ13CG and δ13CATL are
statistically different with higher variance during the pre-MBT (Table S1).
Spectral analyses and variance tests of δ13CPAC
PC1 are similar to δ13CG and
δ13CATL PC1s. The only difference between the
three PC1s is that there is less variance recorded in
δ13CPAC (Fig. 3e). We interpret this muted signal
to be a result of three factors: the large size of the Pacific relative to
the Atlantic, less mixing between water mass end members such as the positive
NADW and more negative AABW, and ocean circulation aging the carbon isotopes
over time leading to more homogenized water masses in the Pacific.
Factor loadings for δ13CATL PC1 are all positive,
suggesting that the time series is representative of the entire Atlantic
basin. In contrast, δ13CATL PC2 yields negative values for
all but the intermediate North Atlantic records and does not show strong
100 kyr spectral power. As such, these results suggest that PC2 exhibits the
dominant mode of variability recorded in the benthic δ13C of
North Atlantic waters shallower than 2000 m depth. Curry and Oppo (2005)
show that NADW formation to below ∼2000 m is reduced in the
North Atlantic during glacial times. The sites with positive factor loadings
in PC2 are located at depths <2000 m, and therefore each site
should remain consistently bathed in NADW through glacial–interglacial
cycles. We thus interpret PC2 as a record of changes in the isotopic values
of the North Atlantic carbon reservoir rather than circulation changes.
During MIS 13, all three δ13C PC1s (global, Atlantic, and
Pacific) demonstrate high positive values. This excursion, first recognized
in individual records by Raymo et al. (1997), clearly stands out
relative to other δ13C interglacial values recorded throughout
the last 800 kyr. The MIS 13 excursion is even more apparent when compared
against other proxy records such as atmospheric CO2, SST, and CH4
(Fig. 7). This high-amplitude change in δ13C values is
similar to the changes recorded in other proxies during MIS 11, yet precedes
the MBT by one glacial cycle. Removal of the MIS 13 interval from variance
tests results in no statistical difference in variance before and after the
MBT, suggesting a large effect of the carbon isotope excursion on these
calculations.
Global ∂13C proxy comparison. Comparison of the
global ∂13C first principal component (PC1; black) to
(a) EPICA Dome C CO2 (yellow; EPICA-community-members,
2004; Lüthi et al., 2008), (b) sea-surface temperature PC1 from
this research (blue), and (c) EPICA Dome C CH4 (red;
EPICA-community-members, 2004).
δ13C gradients
Figure 8 shows regional stacks of δ13C from the deep (>2000 m) and intermediate (<2000 m) North Atlantic and the deep South
Atlantic. As discussed, the intermediate North Atlantic (INA) signal is
predominantly controlled by changes in the carbon reservoir over orbital
timescales. In contrast, the deep North Atlantic (DNA) is controlled by
changes in the relative influence of isotopically more positive NADW and
isotopically more negative AABW, as well as any δ13C
changes to the reservoir that feeds the deep basin from shallower and
surficial waters (i.e., INA). Subtracting the INA from the DNA record (i.e.,
depth gradient) removes the influence of reservoir changes, with the residual
time series reflecting only the relative influences of AABW and NADW on the
isotopic values of carbon in the deep North Atlantic. This is supported by
comparing the North Atlantic depth gradient time series against the South
Atlantic stack (Fig. S4). Both time series demonstrate good correlation for
the entire time interval (r2=0.58), but even more striking is the
similarity in δ13C values, with both time series showing
similar variability and range in δ13C space. The isotopic
similarity between the two records suggests adequate removal of reservoir
influences with the North Atlantic depth gradient, thus reflecting changes in
dominant water mass influence (i.e., circulation). We also note that the
correlation between the two records increases starting at MIS 15 (∼530 ka).
Regional
∂13C stacks. Stacked records of benthic
∂13C separated into three regions: intermediate North
Atlantic (orange), deep North Atlantic (red), and deep South Atlantic (blue).
All plots are shown in ∂13C space to highlight different
isotopic values.
MIS 13 and 5e contour plots of ∂13C. Contour
plots of the ∂13C values in the North Atlantic basin for
the interglacials MIS 13 and MIS 5e. Red colors represent more positive,
enriched values. Blue colors represent lower, depleted values. The plot was created
using Ocean Data Viewer.
Latitudinal ∂13C gradient. (a) North
Atlantic regional ∂13C stack plotted in
∂13C space (red) authigenic εNd (yellow; Howe
et al., 2017). (b) Latitudinal gradient of Atlantic
∂13C regional stacks (North Atlantic minus South Atlantic;
blue). Lower values demonstrate increased similarity between the records.
(c) South Atlantic regional ∂13C stack plotted
in ∂13C space (black). Vertical gray bars indicate glacial
periods. Numbers represent Marine Isotope Stage numbers for interglacials.
The depth gradient does not show the prominent MIS 13 excursion that was
present in the original DNA stack (Fig. 8), suggesting that the excursion
is likely due to a change in the carbon reservoir (represented by the INA)
and not related to ocean circulation. Figure 9 shows contour δ13C plots of the Atlantic basin for MIS 13 and MIS 5e. Although there
is some uncertainty in the these plots due to limited spatial coverage, they
show a clear enrichment of the entire basin during MIS 13 relative to
average post-MBT interglacial conditions, as represented here by MIS 5e. The
global and Pacific δ13C PC1s also show the MIS 13 δ13C excursion, suggesting a change in the global carbon reservoir.
Next, we evaluate the latitudinal gradient between the South Atlantic signal
and the DNA signal in order to further assess the relative influence of the
more negative AABW δ13C values on North Atlantic δ13C values (Fig. 10). Lisiecki (2014) interpreted weaker gradients
during glaciations to reflect shoaling of NADW and greater penetration of
AABW, which could result from reduced NADW formation or stronger AABW
formation. Figure 10b shows a stepwise drop in mean values beginning in MIS 12 (∼436 ka), suggesting a weakening of the gradient due to
greater similarity between North Atlantic and South Atlantic glacial and
interglacial δ13C values.
Discussion
Our new analyses demonstrate that there was a statistically significant
increase in variance in atmospheric CO2, Antarctic temperature, global
SSTs, and bottom water temperature at 436 ka. These changes are consistent
with a transition between two distinct climate states associated with
higher-amplitude interglaciations starting with MIS 11, supporting the notion of a
MBT as defined by Yin (2013). The same climate variables mentioned above
also show an increase in spectral power in the 100 kyr frequency band after
the MBT. On the other hand, the dust analyses suggest that the transition to
greater variability was experienced in the Southern Hemisphere in the
glacial periods starting with MIS 12.
MIS 13 carbon isotope excursion
The PC1 of δ13CG shows a strong correlation with the
CO2 record for most of the last 800 kyr (Fig. 7a). The exception is
during MIS 13, when CO2 levels were still at pre-MBT levels, while
δ13CG shows an anomalously high enrichment relative to
other interglacial values. This is further illustrated by δ13C
contour plots showing that the Atlantic basin was enriched in δ13C during MIS 13 relative to the MIS 5e (Fig. 9).
We evaluated records of biologic activity in various locations of the
Atlantic and Pacific oceans to assess potential sources and sinks in the
carbon system during MIS 13. Ba / Fe from the Antarctic zone (AZ) records the
sedimentary concentration of biogenic Ba and is thus a proxy of organic
matter flux to the deep ocean south of the Polar Front (Jaccard et al., 2013), whereas alkenone
concentrations from the Subantarctic zone (SAZ) indicate export productivity
to the deep ocean in the region north of the Polar Front
(Martínez-Garcia et al., 2009). Based on these proxies, Jaccard
et al. (2013) argued that there were two modes of export productivity in the
Southern Ocean (SO), where high/low export occurs in the AZ during
interglaciations/glaciations, and low/high export occurs in the SAZ during
interglaciations/glaciations. They attributed the increase in SAZ export
productivity to iron fertilization from increased dust accumulation in the
SAZ associated with intensified SO westerlies during glacial periods. Our
Southern Hemisphere dust PC1 record supports this hypothesis in showing that
high values of dust accumulation correlate with increased values of SAZ
export productivity over the last 800 kyr (Fig. S5). We note, however,
that the increase in dust starting at MIS 12 does not have an associated
decrease in glacial CO2 values, suggesting that if iron fertilization
contributed to lower CO2 levels, it had an upper limit beyond which
additional dust fluxes had little effect.
The antiphase relationship between export productivity between the SAZ and
AZ requires a mechanism to increase organic matter productivity in the AZ
during interglaciations as suggested by the Ba / Fe signal (Fig. S5c). In
the modern SO, vertical mixing and upwelling drive the delivery of
nutrient-rich waters necessary for biologic activity to the surface ocean.
Wind-driven upwelling is associated with SO westerlies which shift poleward
during interglaciations (Toggweiler et al., 2006). Thus, any
reduction of upwelling would result from a more northerly position or
decrease in strength of the westerlies; a further decrease in nutrient-rich
surface waters in the AZ during glaciations likely resulted from increased
SO stratification (Sigman et al., 2010;
Jaccard et al., 2013). We note, however, that Jaccard et al. (2013) find no
AZ export productivity during MIS 13, whereas all other interglaciations over
the last 800 kyr show some evidence for it (Fig. S5c). This skipped
interglaciation in export productivity suggests some combination of a change
in the position/strength of the SO westerlies or stratification of the AZ
that limited the delivery of nutrient-rich deep waters to the surface as
compared to other interglaciations of the last 800 kyr.
Marine Isotope Stages 15 to 13 and the carbon isotope excursion.
(a) Summer insolation at 65∘ N (red; Laskar et al., 2004).
(b) First principal component of Atlantic ∂13C
(black). (c) EPICA Dome C CO2 (yellow;
EPICA-community-members, 2004; Lüthi et al., 2008).
(d) Detrended sea-level equivalent from Shakun et al. (2015) (blue).
Derived from ∂18Osw calculations. Negative
numbers indicate lower sea level and increased ice volume.
(e) Chinese Loess Plateau grain size indicating relative Asian
summer monsoon strength (brown; Sun and An, 2005).
(f) Quartz/calcite ratios from site U1313 in the North Atlantic as a
measure of ice-rafted debris (light blue; Naafs et al., 2012). Dark gray bars
highlight the interglacials (MIS 15 and MIS 13) between ∼630 and ∼470 ka. The light gray bar highlights MIS 14.
The PC1s of δ13C (global, Atlantic, and Pacific) demonstrate
that the global ocean was enriched in heavy carbon during MIS 13 relative to
any other interglaciation of the last 800 kyr (Fig. 3). In contrast,
atmospheric CO2 concentrations were ∼ 240 ppm during
MIS 13, similar to other pre-MBT interglacial levels (Fig. 1). Ba / Fe
records of organic export productivity from the AZ that acts as a sink for
light carbon indicate no increase during this interglaciation, while Ca/Al
records from the SAZ indicate increased preservation and thus a deeper
lysocline and lower dissolved inorganic carbon (Jaccard et al., 2010). The
question thus becomes the following: if the ocean is heavily enriched in δ13C
during MIS 13 while CO2 and export productivity remained at low levels,
what reservoir contained the isotopically light carbon?
Paleoclimate records from the CLP indicate greater precipitation during MIS 13 relative to the other interglaciations (Liu, 1985; Yin and Guo,
2008). This greater precipitation has been attributed to increased monsoon
activity recognized throughout monsoonal areas of the Northern Hemisphere
and persisting through MIS 15, 14, and 13 (Yin and Guo, 2008;
Guo et al., 2009). Biogenic silica measurements from Lake Baikal exhibit
continuously high terrestrial productivity in central Asia throughout MIS 15
to MIS 13 (Prokopenko et al., 2002), whereas
sea-level reconstructions indicate that ice volume during MIS 14 was
considerably less relative to other glacial maxima of the last 800 kyr
(Fig. 11d) (Elderfield et al., 2012; Shakun et
al., 2015). Thus, the smaller ice sheets of MIS 14 would likely have had a
lesser effect on displacing forested areas of the Northern Hemisphere,
allowing greater terrestrial carbon storage to potentially persist through a
glacial cycle (Harden et al., 1992). We thus suggest that the increased
monsoonal precipitation and smaller ice volume during MIS 14 would have
combined to increase land biomass that continued into MIS 13. The Northern
Hemisphere thus had the potential to store light carbon in the terrestrial
reservoir resulting in the enriched δ13C MIS 13 signal seen in
the ocean basins (Yin and Guo, 2008).
Ocean circulation changes in the Atlantic basin
One explanation for the glacial–interglacial variations in atmospheric
CO2 invokes a dominant role by the Southern Ocean in storing and
releasing dissolved inorganic carbon (DIC) in the deep Southern Ocean, with
deep-ocean sequestration of atmospheric CO2 occurring through decreased
upwelling and vertical mixing of AABW (Sigman et al.,
2010). Expansion of Southern Ocean sea ice can also lower atmospheric
CO2 by insulating upwelled water from the atmosphere, thus reducing
outgassing, and by increasing the volume of AABW and its capacity to hold
DIC (Stephens and Keeling, 2000; Ferrari et al., 2014). According to this
framework, pre-MBT interglaciations with lower CO2 would be associated
with greater sea-ice extent and a larger volume of AABW, whereas post-MBT
interglaciations with higher CO2 suggest reduced sea-ice extent and
AABW volume. Glacial values of CO2 remain relatively constant
throughout the last 800 kyr (Fig. 1), suggesting that the change in
relative AABW volume before and after the MBT only occurred during
interglaciations.
Average interglacial ∂13C contours. Contour plots
of the average interglacial ∂13C values in the Atlantic
for (a) pre-MBT included MIS 13, (b) pre-MBT excluding
MIS 13 (enriched carbon isotope excursion), and (c) post-MBT. Red
colors indicate higher ∂13C values. Blue colors indicate
lower ∂13C values. The boundary between the two water masses
(NADW and AABW) is indicated at the 0.25 ‰ contour (Curry and Oppo,
2005).
This mechanism is consistent with ice-core evidence for greater sea-ice
extent during pre-MBT interglaciations (Wolff et al., 2006) and with
modeling results that show that interglacial AABW formation decreased after
the MBT through insolation-induced feedbacks on sea ice and surface water
density (Yin, 2013). Moreover, based on the Ba / Fe proxy of
organic matter flux to the deep ocean south of the Polar Front, Jaccard et al. (2013) argued that the deep
Southern Ocean reservoir was larger prior to the MBT.
Our analyses of changes in Atlantic δ13C over the last
800 kyr further support an important role of AABW in causing the post-MBT
increase in interglacial CO2. In particular, the steeper
latitudinal gradient between North and South Atlantic δ13C
records before the MBT reflects greater northward penetration of AABW,
whereas the post-MBT decrease in gradient suggests greater southward
penetration of NADW (Fig. 10b). These gradient changes are further
illustrated by contour plots of average interglacial δ13C
values in the Atlantic which show that prior to the MBT, AABW penetrated
north of the Equator, increasing the δ13C gradient
(Fig. 12a), in contrast to remaining south of the Equator after the MBT,
decreasing the gradient (Fig. 12c). Removal of MIS 13 and its associated
enriched carbon isotope excursion further highlights the greater volume of
AABW in the pre-MBT interglacial Atlantic (Fig. 12b). We note that a record
of the water mass tracer εNd from 6∘ N (Howe
and Piotrowski, 2017) is in good agreement with our North Atlantic regional
δ13C stack (Fig. 10a), with both records suggesting that
changes in volume of the interglacial AABW occurred south of the Equator.
This reorganization of the dominant interglacial water masses in the Atlantic
basin across the MBT, perhaps resulting from insolation-induced feedbacks
(Yin, 2013), would lead to a greater release of deep-ocean CO2
during the post-MBT interglaciations, with corresponding warmer
interglaciations (Fig. 5). An alternative explanation for the observed
decrease in latitudinal gradient could be changes in the isotopic composition
of AABW across this time period. However, modeling results of long-term
carbon fluctuations across this interval suggest that changes in the burial
rate of organic and inorganic carbon caused the δ13C
depletion – the opposite signal necessary to create the increased similarity
between northern- and southern-sourced waters (Hoogakker et al., 2006). Thus,
it is more likely explained by changes in AABW influence north of the
Equator.
Cross-spectral analysis of pre-MBT North and South Atlantic δ13C stacks indicates in-phase coherency between the records at the
eccentricity and obliquity frequencies. Similar tests for the post-MBT
δ13C stacks exhibit coherency at eccentricity, obliquity, and
precession frequencies, with the South Atlantic stack leading the North
Atlantic by ∼23∘ (7 kyr) in eccentricity,
∼18∘ (2 kyr) in obliquity, and
∼36∘ (2 kyr) in precession (Fig. S6). All
phase relationships overlap within uncertainty, suggesting that South
Atlantic δ13C leads North Atlantic δ13C by 2–7 kyr
following the MBT. This lead by the South Atlantic is most apparent during
terminations (Figs. 9, 12) and is most likely related to deglacial
mechanisms for ventilation of respired CO2 from the deep Southern Ocean
such as enhanced wind-driven upwelling or the melting of sea ice in response
to the bipolar seesaw (Cheng et al., 2009).
Conclusions
Using statistical analyses of multiple climate proxies, we have further
characterized the Mid-Brunhes Transition as an increase in interglacial
sea-surface and Antarctic temperatures, atmospheric CO2, and CH4
beginning with MIS 11. At the same time, our new analyses also document a
number of changes in other components of the climate system that began as
early as MIS 14 that suggest a more complex sequence of events prior to the
MBT, although their relationship to the MBT remains unclear. Figure 13
highlights key features in the sequence of events beginning with an increase
in Asian summer monsoon strength during MIS 15 that persisted through MIS 14
and into MIS 13. The strong monsoon strength during MIS 14 is associated
with a weak glaciation, which in combination would have been conducive to a
build-up of Northern Hemisphere land biomass. A continued strong Asian
summer monsoon during MIS 13 associated with greater precipitation would
have further sequestered land biomass and provided a reservoir for light
carbon, resulting in the oceans becoming unusually enriched in δ13C as recorded in the global benthic δ13C carbon isotope
excursion. MIS 12 was associated with the return of large ice sheets,
collapse of the Asian summer monsoon, and the first increase in amplitude of
Southern Hemisphere dust. A decrease in the latitudinal gradient of
interglacial Atlantic δ13C at the MBT suggests a reorganization
of the water masses in the basin and reduction in the size of interglacial
AABW, thus possibly explaining the increase in interglacial values of
atmospheric CO2 with corresponding increases in interglacial SSTs and
CH4. This evidence for a change in AABW is consistent with modeling
results that suggest that the MBT was forced by insolation (Yin, 2013).
Schematic representation of the sequence of events leading to the
Mid-Brunhes Transition. Corresponding Marine Isotope Stages are located on
the left side of each row. Boxes in a row indicate synchronous events.
All data related to this research, including input and analytical output
data, are available at the PANGAEA data publisher (Barth et al., 2018).
Additional data without persistent digital object identifiers used in this
study were retrieved from Bickert et al. (1997), Lisiecki et al. (2008), and
Pisias et al. (1997), and are available upon request.
The supplement related to this article is available online at: https://doi.org/10.5194/cp-14-2071-2018-supplement.
AMB, PUC, and NGP conceived the project. Statistical
tests were conducted by AMB and analyzed by AMB, FH, and NSB. The paper was
written by AMB, PUC, NSB, FH, and NGP.
The authors declare that they have no conflict of
interest.
This article is part of the special issue “Global Challenges
for our Common Future: a paleoscience perspective” – PAGES Young Scientists
Meeting 2017. It is a result of the 3rd Young Scientists Meeting (YSM),
Morillo de Tou, Spain, 7–9 May 2017.
Acknowledgements
We thank the two anonymous reviewers whose comments substantially improved this paper.
Edited by: Robert Barnett
Reviewed by: two anonymous referees
ReferencesBarth, A. M., Clark, P. U., Bill, N. S., He, F., and Pisias, N. G.: Climate
proxies of sea-surface temperature, carbon isotopes, and dust flux for the
last 800 kyr and the Mid-Brunhes Transition, PANGAEA,
10.1594/PANGAEA.894898, 2018.
Bickert, T. and Wefer, G.: Late Quaternary deep water circulation in the
South Atlantic: Reconstruction from carbonate dissolution and benthic stable
isotopes, in: The South Atlantic: Present and Past Circulation, edited by:
Wefer, G., Berger, W. H., Siedler, G., and Webb, D. J., 599–620,
Springer, New York, USA, 1996.
Bickert, T., Curry, W. B., and Wefer, G.: Late Pliocene to Holocene
(2.6–0 Ma) western equatorial Atlantic deep-water circulation: Inferences
from benthic stable isotopes, Proc. Ocean Drill. Program Sci. Results, 154,
239–254, 1997.
Cheng, H., Edwards, R. L., Broecker, W. S., Denton, G. H., Kong, X., Wang, Y.
J., Zhang, R., and Wang, X.: Ice Age Terminations, Science, 326, 248–252,
2009.
Cheng, X., Tian, J., and Wang, P.: Data report: Stable isotopes from Site
1143, in: Proc. ODP, Sci. Results, edited by: Prell, W. L., Wang, P., Blum,
P., Rea, D. K., and Clemens, S. C., 184 (Ocean Drilling Program), College
Station, TX, USA, 2004.Clark, P. U., Archer, D., Pollard, D., Blum, J. D., Rial, J. A., Brovkin, V.,
Mix, A. C., Pisias, N. G., and Roy, M.: The middle Pleistocene transition:
characteristics, mechanisms, and implications for long-term changes in
atmospheric pCO2, Quaternary Sci. Rev., 25, 3150–3184, 2006.Curry, W. B. and Oppo, D. W.: Glacial water mass geometry and the
distribution of δ13C of ΣCO2 in the
western Atlantic Ocean, Paleoceanography, 20, 1–12, 2005.
de Garidel-Thoron, T., Rosenthal, Y., Bassinot, F., and Beaufort, L.: Stable
sea surface temperatures in the western Pacific warm pool over the past 1.75
million years, Nature, 433, 294–298, 2005.deMenocal, P. B., Oppo, D. W., Fairbanks, R. G., and Prell, W. L.:
Pleistocene 13C variability of North Atlantic intermediate water.
Paleoceanography, 7, 229–250, 1992.deMenocal, P. B., Ruddiman, W. F., and Pokras, E. M.: Influences of high- and
low-latitude processes on African terrestrial climate: Pleistocene eolian
records from equatorial Atlantic Ocean Drilling Program Site 663,
Paleoceanography, 8, 209–242, 10.1029/93PA02688, 1993.
Elderfield, H., Ferretti, P., Greaves, M., Crowhurst, S., McCave, I. N.,
Hodell, D. A., and Piotrowski, A. M.: Evolution of Ocean Temperature and Ice
Volume Through the Mid-Pleistocene Climate Transition, Science, 337,
704–709, 2012.
EPICA-community-members: Eight glacial cycles from an Antarctic ice core,
Nature, 429, 623–628, 2004.
Etourneau, J., Martinez, P., Blanz, T., and Schneider, R.:
Pliocene-Pleistocene variability of upwelling activity, productivity, and
nutrient cycling in the Benguela region, Geology, 37, 871–874, 2009.
Ferrari, R., Jansen, M. F., Adkins, J. F., Burke, A., Stewart, A. L., and
Thompson, A. F.: Antarctic sea ice control on ocean circulation in present
and glacial climates, P. Natl. Acad. Sci. USA, 111, 8753–8758, 2014.Guo, Z. T., Berger, A., Yin, Q. Z., and Qin, L.: Strong asymmetry of
hemispheric climates during MIS-13 inferred from correlating China loess and
Antarctica ice records, Clim. Past, 5, 21–31,
10.5194/cp-5-21-2009, 2009.
Harden, J. W., Mark, R. K., Sundquist, E. T., and Stallard, R. F.: Dynamics
of soil carbon during deglaciation of the Laurentide Ice Sheet, Science, 258,
1921–1924, 1992.
Herbert, T. D., Peterson, L. C., Lawrence, K. T., and Liu, Z.: Tropical Ocean
Temperatures Over the Past 3.5 Million Years, Science, 328, 1530–1534, 2010.
Hodell, D., Charles, C., and Ninnemann, U.: Comparison of interglacial stages
in the South Atlantic sector of the southern ocean for the past 450 kyr:
implications for Marine Isotope Stage (MIS) 11, Global Planet. Change, 24,
7–26, 2000.
Hodell, D. A., Charles, C. D., Curtis, J. H., Mortyn, P. G., Ninnemann, U.
S., and Venz, K. A.: Data Report: Oxygen isotope stratigraphy of ODP Leg 177
Sites 1088, 1089, 1090, 1093, and 1094, Proc. Ocean Drill. Program Sci.
Results, 177, 1–26, 2003.Hodell, D. A., Channell, J. E. T., Curtis, J. H., Romero, O. E., and
Röhl, U.: Onset of “Hudson Strait” Heinrich events in the eastern North
Atlantic at the end of the middle Pleistocene transition (∼640 ka)?,
Paleoceanography, 23, PA4218, 10.1029/2008PA001591, 2008.Hoogakker, B. A. A., Rohling, E. J., Palmer, M. R., Tyrrell, T., and
Rothwell, R. G.: Underlying causes for long-term global ocean
δ13C fluctuations over the last 1.20 Myr, Earth Planet. Sc.
Lett., 248, 15–29, 2006.Howe, J. N. W. and Piotrowski, A. M.: Atlantic deep water provenance
decoupled from atmospheric CO2 concentration during the lukewarm
interglacials, Nat. Commun., 8, 1–7, 2017.
Howell, P., Pisias, N. G., Ballance, J., Baughman, J., and Ochs, L.: ARAND
Time-Series Analysis Software, Brown University, Providence, RI, USA, 2006.
Imbrie, J., Berger, A., Boyle, E. A., Clemens, S. C., Duffy, A., Howard, W.
R., Kukla, G., Kutzbach, J., Martinson, D. G., McIntyre, A., Mix, A. C.,
Molfino, B., Morley, J. J., Peterson, L. C., Pisias, N. G., Prell, W. L.,
Raymo, M. E., Shackleton, N. J., and Toggweiler, J. R.: On the structure and
origin of major glaciation cycles 2. The 100 000-year cycle,
Paleoceanography, 8, 699–735, 1993.
Jaccard, S. L., Galbraith, E. D., Sigman, D. M., and Haug, G. H.: A pervasive
link between Antarctic ice core and subarctic Pacific sediment records over
the past 800 kyrs, Quaternary Sci. Rev., 29, 206–212, 2010.
Jaccard, S. L., Hayes, C. T., Martinez-Garcia, A., Hodell, D. A., Anderson,
R. F., Sigman, D. M., and Haug, G. H.: Two modes of change in Southern Ocean
productivity over the past million years, Science, 339, 1419–1423, 2013.
Jansen, J. H. F., Kuijpers, A., and Troelstra, S. R.: A Mid-Brunhes Climatic
Event: Long-Term Changes in Global Atmosphere and Ocean Circulation, Science,
232, 619–622, 1986.
Jouzel, J., Masson-Delmotte, V., Cattani, O., Dreyfus, G., Falourd, S.,
Hoffman, G., Minster, B., Nouet, J., Barnola, J. M., Chappellaz, J., Fischer,
H., Gallet, J. C., Johnsen, S., Leuenberger, D., Loulergue, L., Luethi, D.,
Oerter, H., Parrenin, F., Raisbeck, G. M., Raynaud, D., Schilt, A.,
Schwander, J., Selmo, E., Souchez, R. A., Spahni, R., Stauffer, B.,
Steffensen, J. P., Stenni, B., Stocker, T. F., Tison, J. L., Werner, M., and
Wolff, E. W.: Orbital and Millennial Antarctic Climate Variability over the
Past 800 000 Years, Science, 317, 793–796, 2007.
Kemp, A. E. S., Grigorov, I., Pearce, R. B., and Naveira Garabato, A. C.:
Migration of the Antarctic Polar Front through the mid-Pleistocene
transition: evidence and climatic implications, Quaternary Sci. Rev., 29,
1993–2009, 2010.
Lambert, F., Delmonte, B., Petit, J. R., Bigler, M., Kaufmann, P. R.,
Hutterli, M. A., Stocker, T. F., Ruth, U., Steffensen, J. P., and Maggi, V.:
Dust-climate couplings over the past 800,000 years from the EPICA Dome C ice
core, Nature, 452, 616–619, 2008.
Lamy, F., Gersonde, R., Winckler, G., Esper, O., Jaeschke, A., Kuhn, G.,
Ullermann, J., Martinez-Garcia, A., Lambert, F., and Kilian, R.: Increased
Dust Deposition in the Pacific Southern Ocean During Glacial Periods,
Science, 343, 403–407, 2014.
Laskar, J., Robutel, P., Joutel, F., Gastineau, M., Correia, A. C. M., and
Levrard, B.: A long term numerical solution for the insolation quantities of
the Earth, Astron. Astrophys., 428, 261–285, 2004.Lawrence, K. T., Herbert, T. D., Brown, C. M., Raymo, M. E., and Haywood, A.
M.: High-amplitude variations in North Atlantic sea surface temperature
during the early Pliocene warm period, Paleoceanography, 24, PA2218,
10.1029/2008PA001669, 2009.
Li, L., Li, Q., Tian, J., Wang, P., Wang, H., and Liu, Z.: A 4-Ma record of
thermal evolution in the tropical western Pacific and its implications on
climate change, Earth Planet. Sc. Lett., 309, 10–20, 2011.
Lisiecki, L. E.: Atlantic overturning responses to obliquity and precession
over the last 3 Myr, Paleoceanography, 29, 71–86, 2014.Lisiecki, L. E. and Raymo, M. E.: A Pliocene-Pleistocene stack of 57 globally
distributed benthic δ18O records, Paleoceanography, 20,
1–17, 2005.
Lisiecki, L. E., Raymo, M. E., and Curry, W. B.: Atlantic overturning
responses to Late Pleistocene climate forcings, Nature, 456, 85–88, 2008.
Liu, T. S.: Loess and the Environment, China Ocean Press, Beijing, China,
1985.
Liu, Z. and Herbert, T. D.: High-latitude influence on the eastern equatorial
Pacific climate in the early Pleistocene eopch, Nature, 427, 720–723, 2004.Loulergue, L., Schilt, A., Spahni, R., Masson-Delmotte, V., Blunier, T.,
Lemieux, B., Barnola, J.-M., Raynaud, D., Stocker, T. F., and Chappellaz, J.:
Orbital and millennial-scale features of atmospheric CH4 over the
past 800 000 years, Nature, 453, 383–386, 2008.
Lüthi, D., Le Floch, M., Bereiter, B., Blunier, T., Barnola, J. M.,
Siegenthaler, U., Raynaud, D., Jouzel, J., Fischer, H., Kawamura, K., and
Stocker, T. F.: High-resolution carbon dioxide concentration record
650 000–800 000 years before present, Nature, 453, 379–382, 2008.Martínez-Garcia, A., Rosell-Melé, A., Geibert, W., Gersonde, R.,
Masqué, P., Gaspari, V., and Barbante, C.: Links between iron supply,
marine productivity, sea surface temperature, and CO2 over the last
1.1 Ma, Paleoceanography, 24, 1–14, 2009.
Martinez-Garcia, A., Rosell-Mele, A., Jaccard, S. L., Geibert, W., Sigman, D.
M., and Haug, G. H.: Southern Ocean dust-climate coupling over the past four
million years, Nature, 476, 312–315, 2011.
McIntyre, K., Ravelo, A. C., and Delaney, M. L.: North Atlantic Intermediate
Waters in the late Pliocene to early Pleistocene, Paleoceanography, 14,
324–335, 1999.
Medina-Elizalde, M. and Lea, D. W.: The mid-Pleistocene transition in the
tropical Pacific, Science, 310, 1009–1012, 2005.
Mix, A. C., Le, J., and Shackleton, N. J.: Benthic foraminiferal stable
isotope stratigraphy from Site 846: 0–1.8 Ma, Proc. Ocean Drill. Program
Sci. Results, 138, 839–847, 1995a.
Mix, A. C., Pisias, N. G., Rugh, W., Wilson, J., Morey, A., and Hagelberg, T.
K.: Benthic foraminifer stable isotope record from Site 849 (0–5 Ma): Local
and global climate changes, Proc. Ocean Drill. Program Sci. Results, 138,
371–412, 1995b.Mix, A. C., Morey, A. E., Pisias, N. G., and Hostetler, S. W.: Foraminiferal
faunal estimates of paleotemperature: Circumventing the no-analog problem
yields cool ice age tropics, Paleoceanography, 14, 350–359,
10.1029/1999PA900012, 1999.
Naafs, B. D. A., Hefter, J., Acton, G., Haug, G. H., Martínez-Garcia,
A., Pancost, R., and Stein, R.: Strengthening of North American dust sources
during the late Pliocene (2.7 Ma), Earth Planet. Sc. Lett., 317–318, 8–19,
2012.
Oppo, D. W., Fairbanks, R. G., Gordon, A. L., and Shackleton, N. J.: Late
Pleistocene Southern Ocean 13C variability, Paleoceanography, 5, 43–54,
1990.
Oppo, D. W., McManus, J. F., and Cullen, J. L.: Abrupt climate events
500 000 to 340 000 years ago; evidence from subpolar North Atlantic
sediments, Science, 279, 1335–1338, 1998.
Paillard, D., Labeyrie, L., and Yiou, P.: Macintosh program performs
time-series analysis, Eos Trans. AGU, 77, 379–379, 1996.
Pisias, N. G. and Mix, A. C.: Spatial and temporal oceanographic variability
of the eastern equatorial Pacific during the late Pleistocene: Evidence from
Radiolaria microfossils, Paleoceanography, 12, 381–393, 1997.
Pisias, N. G. and Moore, T. C.: The Evolution of Pleistocene Climate: A Time
Series Approach, Earth Planet. Sc. Lett., 52, 450–458, 1981.
Prokopenko, A. A., Williams, D. F., Kuzmin, D. V., Karabanov, E. B.,
Khursevich, G. K., and Peck, J. A.: Muted climate variations in continental
Siberia during the mid-Pleistocene epoch, Nature, 418, 65–68, 2002.
Prokopenko, A. A., Hinnov, L. A., Williams, D. F., and Kuzmin, M. I.: Orbital
forcing of continental climate during the Pleistocene: a complete
astronomically tuned climatic record from Lake Baikal, SE Siberia, Quaternary
Sci. Rev., 25, 3431–3457, 2006.
Raymo, M. E., Oppo, D., and Curry, W.: The mid-Pleistocene climate
transition: A deep sea carbon isotope perspective, Paleoceanography, 12,
546–559, 1997.Raymo, M. E., Oppo, D. W., Flower, B. P., Hodell, D. A., McManus, J. F.,
Venz, K. A., Kleiven, K. F., and McIntyre, K.: Stability of North Atlantic
water masses in face of pronounced climate variability during the
Pleistocene, Paleoceanography, 19, PA2008, 10.1029/2003PA000921, 2004.
Ruddiman, W. F., Raymo, M. E., Martinson, D. G., Clement, B. M., and Backman,
J.: Pleistocene evolution: Northern Hemisphere ice sheets and North Atlantic
Ocean, Paleoceanography, 4, 353–412, 1989.Russon, T., Elliot, M., Sadekov, A., Cabioch, G., Corrège, T., and De
Deckker, P.: Inter-hemispheric asymmetry in the early Pleistocene Pacific
warm pool, Geophys. Res. Lett., 37, L11601, 10.1029/2010GL043191, 2010.
Schaefer, G., Rodger, J. S., Hayward, B. W., Kennett, J. P., Sabaa, A. T.,
and Scott, G. H.: Planktic foraminiferal and sea surface temperature record
during the last 1 Myr across the Subtropical Front, Southwest Pacific, Mar.
Micropaleontol., 54, 191–212, 2005.
Shackleton, N. J., Berger, A., and Peltier, W. R.: An alternative
astronomical calibration of the Lower Pleistocene timescale based on ODP
Site 677, T. Roy. Soc. Edin. Earth, 81, 251–261, 1990.Shakun, J. D., Lea, D. W., Lisiecki, L. E., and Raymo, M. E.: An 800-kyr
record of global surface ocean δ18O and implications for ice
volume-temperature coupling, Earth Planet. Sc. Lett., 426, 58–68, 2015.Sigman, D. M., Hain, M. P., and Haug, G. H.: The polar ocean and glacial
cycles in atmospheric CO2 concentration, Nature, 466, 47–55, 2010.Stephens, B. B. and Keeling, R. F.: The influence of Antarctic sea ice on
glacial-interglacial CO2 variations, Nature, 404, 171–174, 2000.
Sun, Y. and An, Z.: Late Pliocene-Pleistocene changes in mass accumulation
rates of eolian deposits on the central Chinese Loess Plateau, J. Geophys.
Res., 110, 1–8, 2005.Tiedemann, R., Sarnthein, M., and Shackleton, N. J.: Astronomic timescale for
the Pliocene Atlantic δ18O and dust flux records of Ocean
Drilling Program site 659, Paleoceanography, 9, 619–638, 1994.Toggweiler, J. R., Russell, J. L., and Carson, S. R.: Midlatitude westerlies,
atmospheric CO2, and climate change during the ice ages,
Paleoceanography, 21, 1–15, 2006.
Venz, K. A. and Hodell, D. A.: New evidence for changes in Plio-Pleistocene
deep water circulation from Southern Ocean ODP Leg 177 Site 1090,
Palaeogeogr. Palaeocl., 182, 197–220, 2002.
Venz, K. A., Hodell, D. A., Stanton, C., and Warnke, D. A.: A 1.0 Myr Record
of Glacial North Atlantic Intermediate Water Variability from ODP Site 982 in
the Northeast Atlantic, Paleoceanography, 14, 42–52, 1999.Wefer, G., Berger, W. H., Bickert, T., Donner, B., Fischer, G., Kemle von
Mücke, S., Meinecke, G., Müller, P. J., Mulitza, S., Niebler, H.-S.,
Pätzold, J., Schmidt, H., Schneider, R. R., and Segl, M.: Late Quaternary
surface circulation of the South Atlantic: The stable isotope record and
implications for heat transport and productivity, in: The South Atlantic:
Present and Past Circulation, edited by: Wefer, G., Berger, W. H., Siedler,
G., and Webb, D., Springer, Berlin, Heidelberg, Germany, 461–502, 1996.
Wolff, E. W., Fischer, H., Fundel, F., Ruth, U., Twarloh, B., Littot, G. C.,
Mulvaney, R., Röthlisberger, R., de Angelis, M., Boutron, C. F., Hansson,
M., Jonsell, U., Hutterli, M. A., Lambert, F., Kaufmann, P., Stauffer, B.,
Stocker, T. F., Steffensen, J. P., Bigler, M., Siggaard-Andersen, M. L.,
Udisti, R., Becagli, S., Castellano, E., Severi, M., Wagenback, D., Barbante,
C., Gabrielli, P., and Gaspari, V.: Southern Ocean sea-ice extent,
productivity and iron flux over the past eight glacial cycles, Nature, 440,
491–496, 2006.
Yin, Q.: Insolation-induced mid-Brunhes transition in Southern Ocean
ventilation and deep-ocean temperature, Nature, 494, 222–225, 2013.Yin, Q. Z. and Guo, Z. T.: Strong summer monsoon during the cool MIS-13,
Clim. Past, 4, 29–34, 10.5194/cp-4-29-2008, 2008.