The ice core record of atmospheric CO2 variability during the Last Glacial Period: new insights from timing and isotopes

Abstract

<p>Atmospheric carbon dioxide (CO<sub>2</sub>) concentrations during the last glacial period (70,000 – 23,000 years ago) fluctuated on millennial timescales closely following variations in Antarctic temperature. This close coupling has suggested that the sources and sinks driving millennial scale CO<sub>2</sub> changes are dominated by processes in the Southern Ocean. However, recent work revealed centennial-scale increases in CO<sub>2</sub> during abrupt climate events of the last deglaciation which may represent a second mechanism of carbon cycle variability. </p><p>Here we analyze a high resolution CO<sub>2</sub> record from the last glacial period from the West Antarctic Ice Sheet (WAIS Divide) that precisely defines the timing of CO<sub>2</sub> changes with respect to Antarctic ice core proxies for temperature, dust delivery, and sea-ice extent down to the centennial-timescale. Although CO<sub>2</sub> closely tracks all these proxies over millennia, peak CO<sub>2</sub> levels most often lag behind all proxies by a few hundred years. This decoupling from Antarctic climate variability is most prominent during the onset of DO interstadial events when CO<sub>2</sub>, CH<sub>4</sub> and Greenland temperature all increase simultaneously. Regression analysis suggests that the CO<sub>2</sub> variations can be explained by a combination of two mechanisms: one operating on the time scale of Antarctic climate variability, and a second responding on the Dansgaard-Oeschger time scale.  </p><p>Recent δ<sup>13</sup>C-CO<sub>2</sub> data from the last glacial period support our finding that CO<sub>2</sub> variability is the sum of multiple mechanisms.  The Antarctic climate variability is likely associated with the release of respired organic carbon from the deep ocean.  Superimposed on these oscillations are two types of centennial-scale changes: CO<sub>2</sub> increases and δ<sup>13</sup>C-CO<sub>2</sub> minima in the middle of Heinrich stadials and ii) CO<sub>2</sub> increases and small changes in δ<sup>13</sup>C-CO<sub>2 </sub>that at the onset of DO interstadial event.</p><p>To provide a comprehensive and quantitative constraint on the mechanisms of CO<sub>2</sub> variability during the last glacial period, we run a large suite of transient box model experiments (n = 500) forced with varying combinations of forcings based on proxy time-series (e.g. AABW formation, NADW formation, ocean temperature, dust delivery, and sea-ice extent).  Using data constraints from the ice core records of CO<sub>2</sub>, δ<sup>13</sup>C-CO<sub>2</sub> and mean ocean temperature, we arrive at an ensemble of scenarios that can explain a large amount of the centennial and millennial-scale variability observed in the ice core record. Parsing this into a series of factorial experiments we find that Southern Hemisphere processes can explain 80% of the observed variability and Northern Hemisphere processes account for the remaining 20%.  A further breakdown on the level of individual mechanisms is marred by the high degree of correlation between carbon cycle forcings likely operating in the Southern Hemisphere.  None-the-less, our results highlight how multiple mechanisms operating over multiple timescales may have interacted during the last glacial period to drive changes in atmospheric CO<sub>2</sub>. </p>