Carbon cycle variability during the last millennium and last deglaciation

The exchange of carbon on earth is one of the fundamental processes that sustains life and regulates climate. Since the onset of the Industrial Revolution, the burning of fossil fuels and anthropogenic land conversion have altered the carbon cycle, increasing carbon dioxide in the atmosphere to levels that are unprecedented in the last 800,000 years. This rapid rise in atmospheric carbon dioxide is driving current climate change and further increases are projected to dominate future climate change. However, the fate of the carbon cycle in response to climate change remains uncertain. Insight into how the carbon cycle may change in the future can come from an understanding how it has changed in the past. Key constraints on past carbon cycle variability come from the concentration and stable isotopic composition of atmospheric carbon dioxide recorded in polar ice cores, but reconstructing these histories has been a significant analytical challenge. This thesis presents a new, more precise method for measuring the stable isotopic composition of carbon in carbon dioxide (δ¹³C of CO₂) from polar ice. The new method is then used to reconstruct the atmospheric history of δ¹³C of CO₂ during the last millennium (~770-1900 C.E.) and last deglaciation (~20,000-10,000 years before present). Previously, methods for measuring the δ¹³C of CO₂ had been limited to precision of greater than ±0.05‰. The method presented here combines an ice grater air extraction method and micro-volume equipped dual-inlet mass spectrometer to make high-precision measurements on very small samples of fossil CO₂. The precision as determined by replicate analysis is ±0.018‰. The method also provides high-precision measurements of the CO₂ (±2 ppm) and N2O (±4 ppb). A new high-resolution (~20 year spacing) record of the δ¹³C of CO₂ from 770-1900 C.E is presented that suggests land carbon controlled atmospheric CO₂ variability prior to the Industrial Revolution. A deconvolution of the CO₂ fluxes to the atmosphere provides a well-constrained estimate of the evolution of land carbon stocks. The relationship between climate and land carbon for this time period constrains future climate-carbon cycle sensitivity, but an additional process affecting land carbon is required to explain the data. This missing process may be related to early anthropogenic land cover change or patterns of drought. A long-standing problem in the field of paleoclimatology is a complete mechanistic understanding of the 80 ppm increase in atmospheric CO₂ during the last deglaciation. A horizontal ice core on the Taylor Glacier in Antarctica allowed for the recovery of well-dated, large ice samples spanning the last deglaciation. From this unique archive, a new δ¹³C of CO₂ of very high resolution (50-150 year spacing) is reconstructed. A box model of the carbon cycle is used to construct a framework of the evolution of the carbon cycle during deglaciation. During the Last Glacial Maximum, the lower CO₂ concentration accompanied by only a minor shift in δ¹³C of CO₂ relative to the early Holocene is consistent with a more efficient biological pump in the Southern ocean, limited air-sea gas exchange around Antarctica, and colder ocean temperatures. The temporal evolution of these factors, as informed by timing of proxy data, reconciles the non-linear relationship between CO₂ and δ¹³C of CO₂ from the Last Glacial Maximum to the pre-Industrial. However, the data also reveal very fast changes in δ¹³C of CO₂ that suggest a rapid emission of depleted carbon to the atmosphere on the centennial timescale that is not captured in current models

Author Correction: The potential of gypsum speleothems for paleoclimatology: application to the Iberian Roman Humid Period.

An amendment to this paper has been published and can be accessed via a link at the top of the paper

Quantification of drought during the collapse of the classic Maya civilization.

The demise of Lowland Classic Maya civilization during the Terminal Classic Period (~800 to 1000 CE) is a well-cited example of how past climate may have affected ancient societies. Attempts to estimate the magnitude of hydrologic change, however, have met with equivocal success because of the qualitative and indirect nature of available climate proxy data. We reconstructed the past isotopic composition (δ18O, δD, 17O-excess, and d-excess) of water in Lake Chichancanab, Mexico, using a technique that involves isotopic analysis of the structurally bound water in sedimentary gypsum, which was deposited under drought conditions. The triple oxygen and hydrogen isotope data provide a direct measure of past changes in lake hydrology. We modeled the data and conclude that annual precipitation decreased between 41 and 54% (with intervals of up to 70% rainfall reduction during peak drought conditions) and that relative humidity declined by 2 to 7% compared to present-day conditions.ERC 339694 (Water Isotopes of Hydrated Minerals

Abrupt Holocene ice loss due to thinning and ungrounding in the Weddell Sea Embayment

The extent of grounded ice and buttressing by the Ronne Ice Shelf, which provides resistance to the outflow of ice streams, moderate West Antarctic Ice Sheet stability. During the Last Glacial Maximum, the ice sheet advanced and was grounded near the Weddell Sea continental shelf break. The timing of subsequent ice sheet retreat and the relative roles of ice shelf buttressing and grounding line changes remain unresolved. Here we use an ice core record from grounded ice at Skytrain Ice Rise to constrain the timing and speed of early Holocene ice sheet retreat. Measured δ18O and total air content suggest that the surface elevation of Skytrain Ice Rise decreased by about 450 m between 8.2 and 8.0 kyr before 1950 CE (±0.13 kyr). We attribute this elevation change to dynamic thinning due to flow changes induced by the ungrounding of ice in the area. Ice core sodium concentrations suggest that the ice front of this ungrounded ice shelf then retreated about 270 km (±30 km) from 7.7 to 7.3 kyr before 1950 CE. These centennial-scale changes demonstrate how quickly ice mass can be lost from the West Antarctic Ice Sheet due to changes in grounded ice without extensive ice shelf calving. Our findings both support and temporally constrain ice sheet models that exhibit rapid ice loss in the Weddell Sea sector in the early Holocene

Multiple carbon cycle mechanisms associated with the glaciation of Marine Isotope Stage 4

Here we use high-precision carbon isotope data (δ13C-CO2) to show atmospheric CO2 during Marine Isotope Stage 4 (MIS 4, ~70.5-59 ka) was controlled by a succession of millennial-scale processes. Enriched δ13C-CO2 during peak glaciation suggests increased ocean carbon storage. Variations in δ13C-CO2 in early MIS 4 suggest multiple processes were active during CO2 drawdown, potentially including decreased land carbon and decreased Southern Ocean air-sea gas exchange superposed on increased ocean carbon storage. CO2 remained low during MIS 4 while δ13C-CO2 fluctuations suggest changes in Southern Ocean and North Atlantic air-sea gas exchange. A 7 ppm increase in CO2 at the onset of Dansgaard-Oeschger event 19 (72.1 ka) and 27 ppm increase in CO2 during late MIS 4 (Heinrich Stadial 6, ~63.5-60 ka) involved additions of isotopically light carbon to the atmosphere. The terrestrial biosphere and Southern Ocean air-sea gas exchange are possible sources, with the latter event also involving decreased ocean carbon storage

Spatial pattern of accumulation at Taylor Dome during the last glacial inception: stratigraphic constraints from Taylor Glacier

A new ice core retrieved from the Taylor Glacier blue ice area contains ice and air spanning the Marine Isotope Stage (MIS) 5/4 transition (74 to 65 ka), a period of global cooling and glacial inception. Dating the ice and air bubbles in the new ice core reveals an ice age-gas age difference (Δage) approaching 10 ka during MIS 4, implying very low accumulation at the Taylor Glacier accumulation zone on the northern flank of Taylor Dome. A revised chronology for the Taylor Dome ice core (80 to 55 ka), situated to the south of the Taylor Glacier accumulation zone, shows that Δage did not exceed 2.5 ka at that location. The difference in Δage between the new Taylor Glacier ice core and the Taylor Dome ice core implies a spatial gradient in snow accumulation across Taylor Dome that intensified during the last glacial inception and through MIS 4

Spatial pattern of accumulation at Taylor Dome during Marine Isotope Stage 4: stratigraphic constraints from Taylor Glacier

New ice cores retrieved from the Taylor Glacier (Antarctica) blue ice area contain ice and air spanning the Marine Isotope Stage (MIS) 5–4 transition, a period of global cooling and ice sheet expansion. We determine chronologies for the ice and air bubbles in the new ice cores by visually matching variations in gas- and ice-phase tracers to preexisting ice core records. The chronologies reveal an ice age–gas age difference (Δage) approaching 10 ka during MIS 4, implying very low snow accumulation in the Taylor Glacier accumulation zone. A revised chronology for the analogous section of the Taylor Dome ice core (84 to 55 ka), located to the south of the Taylor Glacier accumulation zone, shows that Δage did not exceed 3 ka. The difference in Δage between the two records during MIS 4 is similar in magnitude but opposite in direction to what is observed at the Last Glacial Maximum. This relationship implies that a spatial gradient in snow accumulation existed across the Taylor Dome region during MIS 4 that was oriented in the opposite direction of the accumulation gradient during the Last Glacial Maximum

Precise interpolar phasing of abrupt climate change during the last ice age

The last glacial period exhibited abrupt Dansgaard–Oeschger climatic oscillations, evidence of which is preserved in a variety of Northern Hemisphere palaeoclimate archives¹. Ice cores show that Antarctica cooled during the warm phases of the Greenland Dansgaard–Oeschger cycle and vice versa[superscript 2,3], suggesting an interhemispheric redistribution of heat through a mechanism called the bipolar seesaw[superscript 4–6]. Variations in the Atlantic meridional overturning circulation (AMOC) strength are thought to have been important, but much uncertainty remains regarding the dynamics and trigger of these abrupt events[superscript 7–9]. Key information is contained in the relative phasing of hemispheric climate variations, yet the large, poorly constrained difference between gas age and ice age and the relatively low resolution of methane records from Antarctic ice cores have so far precluded methane-based synchronization at the required sub-centennial precision[superscript 2,3,10]. Here we use a recently drilled high-accumulation Antarctic ice core to show that, on average, abrupt Greenland warming leads the corresponding Antarctic cooling onset by 218 ± 92 years (2σ) for Dansgaard–Oeschger events, including the Bølling event; Greenland cooling leads the corresponding onset of Antarctic warming by 208 ± 96 years. Our results demonstrate a north-to-south directionality of the abrupt climatic signal, which is propagated to the Southern Hemisphere high latitudes by oceanic rather than atmospheric processes. The similar interpolar phasing of warming and cooling transitions suggests that the transfer time of the climatic signal is independent of the AMOC background state. Our findings confirm a central role for ocean circulation in the bipolar seesaw and provide clear criteria for assessing hypotheses and model simulations of Dansgaard–Oeschger dynamics