Abstract
The mid-latitude westerly winds of the Southern Hemisphere play a central role in the global climate system via Southern Ocean upwelling1, carbon exchange with the deep ocean2, Agulhas leakage (transport of Indian Ocean waters into the Atlantic)3 and possibly Antarctic ice-sheet stability4. Meridional shifts of the Southern Hemisphere westerly winds have been hypothesized to occur5,6 in parallel with the well-documented shifts of the intertropical convergence zone7 in response to Dansgaard–Oeschger (DO) events— abrupt North Atlantic climate change events of the last ice age. Shifting moisture pathways to West Antarctica8 are consistent with this view but may represent a Pacific teleconnection pattern forced from the tropics9. The full response of the Southern Hemisphere atmospheric circulation to the DO cycle and its impact on Antarctic temperature remain unclear10. Here we use five ice cores synchronized via volcanic markers to show that the Antarctic temperature response to the DO cycle can be understood as the superposition of two modes: a spatially homogeneous oceanic ‘bipolar seesaw’ mode that lags behind Northern Hemisphere climate by about 200 years, and a spatially heterogeneous atmospheric mode that is synchronous with abrupt events in the Northern Hemisphere. Temperature anomalies of the atmospheric mode are similar to those associated with present-day Southern Annular Mode variability, rather than the Pacific–South American pattern. Moreover, deuterium-excess records suggest a zonally coherent migration of the Southern Hemisphere westerly winds over all ocean basins in phase with Northern Hemisphere climate. Our work provides a simple conceptual framework for understanding circum-Antarctic temperature variations forced by abrupt Northern Hemisphere climate change. We provide observational evidence of abrupt shifts in the Southern Hemisphere westerly winds, which have previously documented1,2,3 ramifications for global ocean circulation and atmospheric carbon dioxide. These coupled changes highlight the necessity of a global, rather than a purely North Atlantic, perspective on the DO cycle.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
£14.99 / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
£199.00 per year
only £3.90 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
Source Data (WDC sulfur data, volcanic tie points and water isotope data on synchronized chronologies) and derived products (stacks, PCA results, etc.) are available in the online version of the paper and in the NOAA palaeoclimate data archive.
References
Marshall, J. & Speer, K. Closure of the meridional overturning circulation through Southern Ocean upwelling. Nat. Geosci. 5, 171–180 (2012).
Toggweiler, J. R., Russell, J. L. & Carson, S. R. Midlatitude westerlies, atmospheric CO2, and climate change during the ice ages. Paleoceanography 21, PA2005 (2006).
Biastoch, A., Boning, C. W., Schwarzkopf, F. U. & Lutjeharms, J. R. E. Increase in Agulhas leakage due to poleward shift of Southern Hemisphere westerlies. Nature 462, 495–498 (2009).
Pritchard, H. D. et al. Antarctic ice-sheet loss driven by basal melting of ice shelves. Nature 484, 502–505 (2012).
Lee, S. Y., Chiang, J. C., Matsumoto, K. & Tokos, K. S. Southern Ocean wind response to North Atlantic cooling and the rise in atmospheric CO2: modeling perspective and paleoceanographic implications. Paleoceanography 26, PA1214 (2011).
Pedro, J. B. et al. Beyond the bipolar seesaw: toward a process understanding of interhemispheric coupling. Quat. Sci. Rev. 192, 27–46 (2018).
Schneider, T., Bischoff, T. & Haug, G. H. Migrations and dynamics of the intertropical convergence zone. Nature 513, 45–53 (2014).
Markle, B. R. et al. Global atmospheric teleconnections during Dansgaard–Oeschger events. Nat. Geosci. 10, 36–40 (2017).
Ding, Q., Steig, E. J., Battisti, D. S. & Kuttel, M. Winter warming in West Antarctica caused by central tropical Pacific warming. Nat. Geosci. 4, 398–403 (2011).
Buiron, D. et al. Regional imprints of millennial variability during the MIS 3 period around Antarctica. Quat. Sci. Rev. 48, 99–112 (2012).
Stocker, T. F. & Johnsen, S. J. A minimum thermodynamic model for the bipolar seesaw. Paleoceanography 18, 1087 (2003); correction 20, PA1002 (2005).
Rind, D. et al. Effects of glacial meltwater in the GISS coupled atmosphere–ocean model. 2. A bipolar seesaw in Atlantic Deep Water production. J. Geophys. Res. Atmos. 106, 27355–27365 (2001).
Kostov, Y. et al. Fast and slow responses of Southern Ocean sea surface temperature to SAM in coupled climate models. Clim. Dyn. 48, 1595–1609 (2017).
Jouzel, J. et al. Validity of the temperature reconstruction from water isotopes in ice cores. J. Geophys. Res. 102, 26471–26487 (1997).
WAIS Divide Project Members. Precise interpolar phasing of abrupt climate change during the last ice age. Nature 520, 661–665 (2015).
Merlivat, L. & Jouzel, J. Global climatic interpretation of the deuterium–oxygen 18 relationship for precipitation. J. Geophys. Res. 84, 5029–5033 (1979).
Masson-Delmotte, V. et al. GRIP deuterium excess reveals rapid and orbital-scale changes in Greenland moisture origin. Science 309, 118–121 (2005).
Masson-Delmotte, V. et al. Abrupt change of Antarctic moisture origin at the end of Termination II. Proc. Natl Acad. Sci. USA 107, 12091–12094 (2010).
Schmidt, G. A., LeGrande, A. N. & Hoffmann, G. Water isotope expressions of intrinsic and forced variability in a coupled ocean-atmosphere model. J. Geophys. Res. 112, D10103 (2007).
Uemura, R. et al. Ranges of moisture-source temperature estimated from Antarctic ice cores stable isotope records over glacial–interglacial cycles. Clim. Past 8, 1109–1125 (2012).
Sodemann, H. & Stohl, A. Asymmetries in the moisture origin of Antarctic precipitation. Geophys. Res. Lett. 36, L22803 (2009).
Lefebvre, W., Goosse, H., Timmermann, R. & Fichefet, T. Influence of the Southern Annular Mode on the sea ice–ocean system. J. Geophys. Res. 109, C09005 (2004).
Thompson, D. W. J. & Wallace, J. M. Annular modes in the extratropical circulation. Part I: month-to-month variability. J. Clim. 13, 1000–1016 (2000).
Sen Gupta, A. & England, M. H. Coupled ocean-atmosphere-ice response to variations in the Southern Annular Mode. J. Clim. 19, 4457–4486 (2006).
Liu, Z. et al. Transient simulation of last deglaciation with a new mechanism for Bolling–Allerod warming. Science 325, 310–314 (2009).
Ferreira, D., Marshall, J., Bitz, C. M., Solomon, S. & Plumb, A. Antarctic Ocean and sea ice response to ozone depletion: a two-time-scale problem. J. Clim. 28, 1206–1226 (2015).
Raisbeck, G. M. et al. An improved north–south synchronization of ice core records around the 41 kyr 10Be peak. Clim. Past 13, 217–229 (2017).
Grootes, P. M., Stuiver, M., White, J. W. C., Johnsen, S. & Jouzel, J. Comparison of oxygen isotope records from the GISP2 and GRIP Greenland ice cores. Nature 366, 552–554 (1993).
Rhodes, R. H. et al. Enhanced tropical methane production in response to iceberg discharge in the North Atlantic. Science 348, 1016–1019 (2015).
Dee, D. P. et al. The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 137, 553–597 (2011).
Severi, M. et al. Synchronisation of the EDML and EDC ice cores for the last 52 kyr by volcanic signature matching. Clim. Past 3, 367–374 (2007).
Severi, M., Udisti, R., Becagli, S., Stenni, B. & Traversi, R. Volcanic synchronisation of the EPICA-DC and TALDICE ice cores for the last 42 kyr BP. Clim. Past 8, 509–517 (2012).
Fujita, S., Parrenin, F., Severi, M., Motoyama, H. & Wolff, E. W. Volcanic synchronization of Dome Fuji and Dome C Antarctic deep ice cores over the past 216 kyr. Clim. Past 11, 1395–1416 (2015).
Sigl, M. et al. Insights from Antarctica on volcanic forcing during the Common Era. Nat. Clim. Chang. 4, 693–697 (2014).
Parrenin, F. et al. Volcanic synchronisation between the EPICA Dome C and Vostok ice cores (Antarctica) 0–145 kyr BP. Clim. Past 8, 1031–1045 (2012).
Veres, D. et al. The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years. Clim. Past 9, 1733–1748 (2013).
Buizert, C. et al. The WAIS Divide deep ice core WD2014 chronology – part 1: methane synchronization (68–31 ka BP) and the gas age–ice age difference. Clim. Past 11, 153–173 (2015).
Sigl, M. et al. The WAIS Divide deep ice core WD2014 chronology – part 2: annual-layer counting (0–31 ka BP). Clim. Past 12, 769–786 (2016).
Kawamura, K. et al. Northern Hemisphere forcing of climatic cycles in Antarctica over the past 360,000 years. Nature 448, 912–916 (2007).
Bazin, L. et al. An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120–800 ka. Clim. Past 9, 1715–1731 (2013).
Parrenin, F. et al. IceChrono1: a probabilistic model to compute a common and optimal chronology for several ice cores. Geosci. Model Dev. 8, 1473–1492 (2015).
Lemieux-Dudon, B. et al. Consistent dating for Antarctic and Greenland ice cores. Quat. Sci. Rev. 29, 8–20 (2010).
Fudge, T. J. et al. Variable relationship between accumulation and temperature in West Antarctica for the past 31,000 years. Geophys. Res. Lett. 43, 3795–3803 (2016).
Fortuin, J. & Oerlemans, J. Parameterization of the annual surface temperature and mass balance of Antarctica. Ann. Glaciol. 14, 78–84 (1990).
WAIS Divide Project Members. Onset of deglacial warming in West Antarctica driven by local orbital forcing. Nature 500, 440–444 (2013).
Jones, T. R. et al. Water isotope diffusion in the WAIS Divide ice core during the Holocene and last glacial. J. Geophys. Res. Earth 122, 290–309 (2017).
EPICA Community Members. One-to-one coupling of glacial climate variability in Greenland and Antarctica. Nature 444, 195–198 (2006).
Stenni, B. et al. The deuterium excess records of EPICA Dome C and Dronning Maud Land ice cores (East Antarctica). Quat. Sci. Rev. 29, 146–159 (2010).
Jouzel, J. et al. Orbital and millennial Antarctic climate variability over the past 800,000 years. Science 317, 793–796 (2007).
Kawamura, K. et al. State dependence of climatic instability over the past 720,000 years from Antarctic ice cores and climate modeling. Sci. Adv. 3, (2017).
Uemura, R. et al. Asynchrony between Antarctic temperature and CO2 associated with obliquity over the past 720,000 years. Nat. Commun. 9, 961 (2018).
Stenni, B. et al. Expression of the bipolar see-saw in Antarctic climate records during the last deglaciation. Nat. Geosci. 4, 46–49 (2011).
Landais, A. et al. A review of the bipolar see–saw from synchronized and high resolution ice core water stable isotope records from Greenland and East Antarctica. Quat. Sci. Rev. 114, 18–32 (2015).
Baumgartner, M. et al. NGRIP CH4 concentration from 120 to 10 kyr before present and its relation to a δ15N temperature reconstruction from the same ice core. Clim. Past 10, 903–920 (2014).
Mudelsee, M. Break function regression. Eur. Phys. J. Spec. Top. 174, 49–63 (2009).
Mudelsee, M. Ramp function regression: a tool for quantifying climate transitions. Comput. Geosci. 26, 293–307 (2000).
Clark, P. U., Pisias, N. G., Stocker, T. F. & Weaver, A. J. The role of the thermohaline circulation in abrupt climate change. Nature 415, 863–869 (2002).
Grootes, P. M. & Stuiver, M. Oxygen 18/16 variability in Greenland snow and ice with 10−3-to-105-year time resolution. J. Geophys. Res. 102, 26455–26470 (1997).
Alley, R. B., Anandakrishnan, S. & Jung, P. Stochastic resonance in the North Atlantic. Paleoceanography 16, 190–198 (2001).
Buizert, C. & Schmittner, A. Southern Ocean control of glacial AMOC stability and Dansgaard–Oeschger interstadial duration. Paleoceanography 30, 1595–1612 (2015).
Ditlevsen, P. D., Kristensen, M. S. & Andersen, K. K. The recurrence time of Dansgaard–Oeschger events and limits on the possible periodic component. J. Clim. 18, 2594–2603 (2005).
He, F. et al. Northern Hemisphere forcing of Southern Hemisphere climate during the last deglaciation. Nature 494, 81–85 (2013).
Liu, Z. et al. Younger Dryas cooling and the Greenland climate response to CO2. Proc. Natl Acad. Sci. 109 11101–11104 (2012).
He, F. Simulating Transient Climate Evolution of the Last Deglaciation with CCSM3. PhD thesis, Univ. Wisconsin-Madison (2010).
Pedro, J. B. et al. The spatial extent and dynamics of the Antarctic Cold Reversal. Nat. Geosci. 9, 51–55 (2016).
Sodemann, H., Schwierz, C. & Wernli, H. Interannual variability of Greenland winter precipitation sources: Lagrangian moisture diagnostic and North Atlantic Oscillation influence. J. Geophys. Res. 113, D03107 (2008).
Läderach, A. & Sodemann, H. A revised picture of the atmospheric moisture residence time. Geophys. Res. Lett. 43, 924–933 (2016).
Stohl, A., Forster, C., Frank, A., Seibert, P. & Wotawa, G. Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2. Atmos. Chem. Phys. 5, 2461–2474 (2005).
Locarnini, R. A. et al. World Ocean Atlas 2013, Volume 1: Temperature (eds Levitus, S. & Mishonov, A.) (NOAA/NESDIS, Silver Spring, 2013).
Kalnay, E. et al. The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc. 77, 437–471 (1996).
Acknowledgements
This work is funded by the US National Science Foundation (NSF), through grants ANT-1643394 (to C.B. and J.J.W.), ANT-1643355 (to T.J.F. and E.J.S) and AGS-1502990 (to F.H.); the Swiss National Science Foundation through grant 200021_143436 (to H.S.); the CNRS/INSU/LEFE projects IceChrono and CO2Role (to F.P.); JSPS KAKENHI grants 15H01731 (to K.G.-A., H.M. and M.H.), 15KK0027 (to K.K.) and 26241011 (to K.K., S.F. and H.M.); MEXT KAKENHI grant 17H06320 (to K.K., H.M. and R.U.); the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007–2013)/ERC grant agreement 610055 (to J.B.P.); and the NOAA Climate and Global Change Postdoctoral Fellowship programme, administered by the University Corporation for Atmospheric Research (F.H.). We acknowledge high-performance computing support from Yellowstone (ark:/85065/d7wd3xhc) provided by NCAR's Computational and Information Systems Laboratory, sponsored by the NSF. This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the US Department of Energy under contract number DE-AC05-00OR22725. This is TALDICE publication number 52.
Reviewer information
Nature thanks N. Abram, S. Davies and T. Stocker for their contribution to the peer review of this work.
Author information
Authors and Affiliations
Contributions
Data analysis by C.B., M. Severi, M. Sigl and J.J.W.; manuscript preparation by C.B.; volcanic ice-core synchronization by M. Sigl, M. Severi, C.B., J.R.M., F.P., S.F. and T.J.F.; GCM simulations and interpretation by F.H., J.B.P. and J.J.W.; moisture tagging/tracing experiments by B.R.M. and H.S.; ice-core water isotope analysis by K.G.-A., K.K., H.M., M.H., R.U., B.S. and E.J.S.; all authors discussed the results and contributed towards improving the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Volcanic synchronization of Antarctic ice cores.
a, Age offset between the WD201437,38 (WDC) and AICC201236 (TAL, EDML, EDC) age scales, with each dot representing a volcanic tie point. Yellow and blue triangles denote the timing of TAL–EDC and EMDL–EDC volcanic ties32,36, respectively. b, Overview of synchronizations between the ice cores used in this study. Grey arrows indicate previously published synchronizations and coloured arrows denote synchronizations performed here. Synchronizations within Antarctica are based purely on volcanic links; synchronization between WDC and NGRIP (Greenland) are based on atmospheric CH4 (green arrow).
Extended Data Fig. 2 Age uncertainty due to volcanic synchronization.
a, Interpolation uncertainty (1σ) for four different values of L (the spacing between two adjacent volcanic tie points), based on the layer-counted WD2014 age scale38. b, Interpolation uncertainty in synchronizing the four ice cores to WD2014 chronology. Grey vertical lines give the timing of DO events.
Extended Data Fig. 3 Site-specific stacks of δ18O and dln.
a, Stack of NGRIP δ18O (teal; left axis) and WDC CH4 (green; right axis) during DO warming. b, As in a, but for DO cooling. c, Stack of Antarctic δ18O at the indicated locations (see key) during DO warming. d, As in c, but for DO cooling. e, Stack of Antarctic dln at the indicated locations during DO warming. f, As in e, but for DO cooling. Isotope ratios are on the VSMOW scale.
Extended Data Fig. 4 Proportionality of the atmospheric response.
a–f, Comparison of stacks of just the major DO events (those following Heinrich events, namely, DO events 0, 1, 4, 8, 12 and 14; a, c, e) and just the minor DO events (the remainder; b, d, f). a, b, Stacks of NGRIP δ18O (left axis) and CH4 (right axis). c, d, Stacks of Antarctic δ18O at the indicated locations. e, f, As in c and d, but with the Antarctic mean subtracted. g, Proportionality of the atmospheric response for individual events (numbered). The NGRIP event size is found via regression of individual NGRIP events to the multi-event NGRIP δ18O stack normalized to unit variance (Fig. 2a). The Antarctic atmospheric response is found via multiple linear regression of single-site individual events to the atmospheric and oceanic modes (Fig. 2e). Shown are the average (dots) and standard deviation (error bars) of the response at EDC, DF and EDML (the cores with the strongest atmospheric response); the EDML response is multiplied by −1 because it has the opposite sign of the response at DF and EDC. Red and blue dots denote the major and minor DO events, respectively. Isotope ratios are on the VSMOW scale.
Extended Data Fig. 5 Reduction of the number of events in the δ18O stacks.
a, PC1, when stacking n = 2 and n = 8 randomly selected events. The thick line and shaded area represent the mean and ±1σ, respectively, obtained from 50,000 runs. The vertical yellow bands denote the 200-yr period after the abrupt DO event at t = 0. b, As in a, but for PC2. c, Fraction of signal variance explained by PC1 and PC2, when stacking 2 or 8 randomly selected events. Colour coding as in a and b.
Extended Data Fig. 6 Robustness of the atmospherically forced warming pattern.
a, PCA as a function of window length, with the fraction of variance explained by PC1 and PC2. b, Comparison of PC1 at a 400-yr window length to PC2 at a 2,000-yr window length, showing the crossover of the atmospheric response (that is, from PC1 to PC2) as a function of window length. c, Lag time of the Antarctic PC1 and PC2 response as a function of window length, assessed using the BREAKFIT55 and RAMPFIT56 routines (see Methods for the selection of routine). d, EOF1 at a 400-yr window length, expressed in units of δ18O (‰). e, EOF2 at a 2,000-yr window length, expressed in units of δ18O (‰). f, Slope of linear fit to δ18O stacks in the interval t = 0 to t = +200 yr, shown as the change (in ‰) during these 200 years. Isotope ratios are on the VSMOW scale.
Extended Data Fig. 7 Antarctic climate response to DO cooling.
a, Stack of NGRIP δ18O. b, Stack of Antarctic δ18O at indicated locations. c, As in b, but with the Antarctic mean subtracted. d, First two principal components of the Antarctic δ18O stacks, with the fraction of variance explained (offset for clarity). The lines show the BREAKFIT (PC1) and RAMPFIT (PC2) fits. e, f, Empirical orthogonal functions EOF1 and EOF2 associated with PC1 and PC2 in d, scaled to show the magnitude in units of δ18O (‰). Isotope ratios are on the VSMOW scale.
Extended Data Fig. 8 Moisture sources of Antarctic ice cores and the SAM.
a, Mass-weighted probability distribution functions of Antarctic moisture sources for the five ice cores of interest (5 × 10−5 deg−2 contour lines; the area-integrated probability distribution function equals 1). Distributions are calculated from reanalysis data30,67 using a Lagrangian source diagnostic21,66 (Methods). Parallels are plotted in 15° increments of latitude and meridians in 45° increments of longitude. b, SST (black)69 and relative humidity (RH; grey)70 as a function of latitude. The coloured curves give the latitudinal source distribution during a negative SAM phase (solid curves; SAM index <0) and a positive SAM phase (dashed curves; SAM index >0). The solid and open dots show the first moment of the source distribution during negative and positive SAM phases, respectively. We note that during a positive SAM phase, moisture sources for all core locations are located closer to the Antarctic continent. Source distribution data were obtained using water tagging experiments8 in the CAM (Methods).
Extended Data Fig. 9 SAM-like variability in zonal near-surface winds.
a, ERA-Interim reanalysis (1979–2016 annual means) zonal wind speed at a height of 10 m, regressed onto the SAM index (here the first principal component of sea-level pressure variability south of 20° S), expressed in metres per second per standard deviation in the index. b, As in a, but for internal variability in the CCSM3 TraCE model simulation25,62 during glacial climate before freshwater forcing of Heinrich stadial 1 (19.5 kyr bp to 19.01 kyr bp, decadal means). c, As in a, but for the response to North Atlantic freshwater forcing in the CCSM3 TraCE model (19.1 kyr bp to 18.9 kyr bp, decadal means, with the freshwater forcing applied at 19 kyr bp).
Supplementary information
Supplementary Data
This file contains Supplementary Data Tables 1-11. Full descriptions are provided within the data sheets
Rights and permissions
About this article
Cite this article
Buizert, C., Sigl, M., Severi, M. et al. Abrupt ice-age shifts in southern westerly winds and Antarctic climate forced from the north. Nature 563, 681–685 (2018). https://doi.org/10.1038/s41586-018-0727-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-018-0727-5
Keywords
This article is cited by
-
Two decades of deep ice cores from Antarctica
Nature (2024)
-
Near-synchronous Northern Hemisphere and Patagonian Ice Sheet variation over the last glacial cycle
Nature Geoscience (2024)
-
Antarctic evidence for an abrupt northward shift of the Southern Hemisphere westerlies at 32 ka BP
Nature Communications (2023)
-
Bipolar impact and phasing of Heinrich-type climate variability
Nature (2023)
-
The PaleoJump database for abrupt transitions in past climates
Scientific Reports (2023)