A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Climate extremes and the carbon cycle

The terrestrial biosphere is a key component of the global carbon cycle and its carbon balance is strongly influenced by climate. Continuing environmental changes are thought to increase global terrestrial carbon uptake. But evidence is mounting that climate extremes such as droughts or storms can lead to a decrease in regional ecosystem carbon stocks and therefore have the potential to negate an expected increase in terrestrial carbon uptake. Here we explore the mechanisms and impacts of climate extremes on the terrestrial carbon cycle, and propose a pathway to improve our understanding of present and future impacts of climate extremes on the terrestrial carbon budget.

Links between atmospheric carbon dioxide, the land carbon reservoir and climate over the past millennium

The stability of terrestrial carbon reservoirs is thought to be closely linked to variations in climate 1, but the magnitude of carbon–climate feedbacks has proved dificult to constrain for both modern 2–4 and millennial 5–13 timescales. Reconstructions of atmospheric CO2 concentrations for the past thousand years have shown fluctuations on multidecadal to centennial timescales 5–7, but the causes of these fluctuations are unclear. Here we report high-resolution carbon isotope measurements of CO2 trapped within the ice of the West Antarctic Ice Sheet Divide ice core for the past 1,000 years. We use a deconvolution approach 14 to show that changes in terrestrial organic carbon stores best explain the observed multidecadal variations in the 13 C of CO2 and in CO2 concentrations from 755 to 1850 CE. If significant long-term carbon emissions came from pre-industrial anthropogenic land-use changes over this interval, the emissions must have been offset by a natural terrestrial sink for 13 C-depleted carbon, such as peatlands. We find that on multidecadal timescales, carbon cycle changes seem to vary with reconstructed regional climate changes. We conclude that climate variability could be an important control of fluctuations in land carbon storage on these timescales.

Comparative carbon cycle dynamics of the present and last interglacial

Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO₂ remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO₂ fertilization, land use, wild fire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO₂ dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO₂ and δ¹³ CO₂ changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO₂ dynamics from 8ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO₂ changes after 122 ka BP. This failure to simulate late-Eemian CO₂ dynamics could be a result of the imposed forcings such as prescribed CaCO₃ accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO₂ dynamics eshallow water CaCO₃ accumulation, peat and permafrost carbon dynamics are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO₂ dynamics.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (d13Catm ), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of d13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 years BP). The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4 permil shift to heavier values between the mean d13Catm level in the Penultimate (~ 140 000 years BP) and Last Glacial Maximum (~ 22 000 years BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 years, but with different phasing and magnitudes. Furthermore, a 5000 years lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS 5.5 (120 000 years BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Stable carbon and oxygen isotope record of planktonic foraminifera across the Paleocene-Eocene boundary

The early Cenozoic marine carbon isotopic record is marked by a long-term shift from high d13C values in the late Paleocene to values that are 2 to 3 lower in the early Eocene. The shift is recorded in fossil carbonates from each ocean basin and represents a large change in the distribution of 12C between the ocean and other carbon reservoirs. Superimposed upon this long-term shift are several distinct carbon isotopic negative excursions that are also recorded globally. These carbon isotopic 'events' near the Paleocene-Eocene boundary provide strati-graphic information that can facilitate intersite correlations between marine and non-marine sequences. Here we present a detailed marine carbon isotopic stratigraphy across the Paleocene-Eocene boundary that is constrained by calcareous nannofossil and planktonic foraminifera bio-stratigraphy and magnetostratigraphy. We show that several distinct carbon isotopic changes are recorded in uppermost Paleocene and lowermost Eocene marine biogenic carbonate sediments. At least one of these isotopic changes in the ocean's carbon isotopic composition was transmitted to terrestrial carbon reservoirs, including plant biomass via atmospheric CO2. As a consequence of this exchange of 12C between the ocean and terrestrial carbon reservoirs, it is possible to use carbon isotope stratigraphy to correlate the uppermost Paleocene and lowermost Eocene non-fossiliferous terrestrial sediments of the Paris Basin with marine sequences.

Seawater carbonate chemistry and carbon and oxygen isotopes during experiments with planktonic foraminifera Orbulina universa and Globigerina bulloides, 1997

Stable oxygen and carbon isotope measurements on biogenic calcite and aragonite have become standard tools for reconstructing past oceanographic and climatic change. In aquatic organisms, 18O/16O ratios in the shell carbonate are a function of the ratio in the sea water and the calcification temperature (Epstein et al., 1953). In contrast, 13C/12C ratios are controlled by the ratio of dissolved inorganic carbon in sea water and physiological processes such as respiration and symbiont photosynthesis (Spero et al., 1991, doi:10.1029/91PA02022). These geochemical proxies have been used with analyses of foraminifera shells to reconstruct global ice volumes (Shackleton and Opdyke, 1973, doi:10.1016/0033-5894(73)90052-5), surface and deep ocean temperatures (Broecker, 1986, doi:10.1016/0033-5894(86)90087-6; Labeyrie et al., 1987, doi:10.1038/327477a0), ocean circulation changes (Duplessy et al., 1988, doi:10.1029/PA003i003p00343) and glacial-interglacial exchange between the terrestrial and oceanic carbon pools (Sackleton, 1977). Here, we report experimental measurements on living symbiotic and non-symbiotic plankton foraminifera (Orbulina universa and Globigerina bulloides respectively) showing that the 13C/12C and 18O/16O ratios of the calcite shells decrease with increasing seawater [CO3 2-]. Because glacial-period oceans had higher pH and [CO3 2-] than today (Sanyal et al., 1995, doi:10.1038/373234a0), these new relationships confound the standard interpretation of glacial foraminiferal stable-isotope data. In particular, the hypothesis that the glacial-interglacial shift in the 13C/12C ratio was due to a transfer of terrestrial carbon into the ocean(Shackleton ,1977) can be explained alternatively by an increase in ocean alkalinity (Lea et al., 1996). A carbonate-concentration effect could also help explain some of the extreme stable-isotope variations during the Proterozoic and Phanerozoic aeons (Kaufman et al., 1993, doi:10.1016/0012-821X(93)90254-7).

Potential responses of soil organic carbon to global environmental change

Recent improvements in our understanding of the dynamics of soil carbon have shown that 20–40% of the approximately 1,500 Pg of C stored as organic matter in the upper meter of soils has turnover times of centuries or less. This fast-cycling organic matter is largely comprised of undecomposed plant material and hydrolyzable components associated with mineral surfaces. Turnover times of fast-cycling carbon vary with climate and vegetation, and range from <20 years at low latitudes to >60 years at high latitudes. The amount and turnover time of C in passive soil carbon pools (organic matter strongly stabilized on mineral surfaces with turnover times of millennia and longer) depend on factors like soil maturity and mineralogy, which, in turn, reflect long-term climate conditions. Transient sources or sinks in terrestrial carbon pools result from the time lag between photosynthetic uptake of CO2 by plants and the subsequent return of C to the atmosphere through plant, heterotrophic, and microbial respiration. Differential responses of primary production and respiration to climate change or ecosystem fertilization have the potential to cause significant interrannual to decadal imbalances in terrestrial C storage and release. Rates of carbon storage and release in recently disturbed ecosystems can be much larger than rates in more mature ecosystems. Changes in disturbance frequency and regime resulting from future climate change may be more important than equilibrium responses in determining the carbon balance of terrestrial ecosystems.

Evolution of the stable carbon isotope composition of atmospheric CO₂ over the last glacial cycle

We present new δ¹³C measurements of atmospheric CO₂ covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in δ¹³C(atm) of 0.5‰ occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in δ¹³C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO₂]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of δ¹³C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.

A new stepwise carbon cycle data assimilation system using multiple data streams to constrain the simulated land surface carbon cycle

Acknowledgements. This work was mainly funded by the EU FP7 CARBONES project (contracts FP7-SPACE-2009-1-242316), with also a small contribution from GEOCARBON project (ENV.2011.4.1.1-1-283080). This work used eddy covariance data acquired by the FLUXNET community and in particular by the following networks: AmeriFlux (U.S. Department of Energy, Biological and Environmental Research, Terrestrial Carbon Program; DE-FG02-04ER63917 and DE-FG02-04ER63911), AfriFlux, AsiaFlux, CarboAfrica, CarboEuropeIP, CarboItaly, CarboMont, ChinaFlux, Fluxnet-Canada (supported by CFCAS, NSERC, BIOCAP, Environment Canada, and NRCan), GreenGrass, KoFlux, LBA, NECC, OzFlux, TCOS-Siberia, USCCC. We acknowledge the financial support to the eddy covariance data harmonization provided by CarboEuropeIP, FAO-GTOS-TCO, iLEAPS, Max Planck Institute for Biogeochemistry, National Science Foundation, University of Tuscia, Université Laval and Environment Canada and US Department of Energy and the database development and technical support from Berkeley Water Center, Lawrence Berkeley National Laboratory, Microsoft Research eScience, Oak Ridge National Laboratory, University of California-Berkeley, University of Virginia. Philippe Ciais acknowledges support from the European Research Council through Synergy grant ERC-2013-SyG-610028 “IMBALANCE-P”. The authors wish to thank M. Jung for providing access to the GPP MTE data, which were downloaded from the GEOCARBON data portal (https://www.bgc-jena.mpg.de/geodb/projects/Data.php). The authors are also grateful to computing support and resources provided at LSCE and to the overall ORCHIDEE project that coordinate the development of the code (http://labex.ipsl.fr/orchidee/index.php/about-the-team).

(Table 1) Oxygen and carbon isotopes for foraminiferal taxa at DSDP Hole 76-540

We analyzed the oxygen and carbon isotopic composition of planktonic and benthic foraminifers picked from 13 late Eocene to late Oligocene samples from DSDP Site 540 (23°49.73'N, 84°22.25'W, 2926 m water depth) from the Gulf of Mexico. An enrichment occurs in 18O of about 0.5 to 0.8 per mil in both benthic foraminifers and surface-dwelling planktonic foraminifers between the latest Eocene and early Oligocene. This early Oligocene maximum is followed by lower 18O values. A 1.2 per mil d13C decrease in both benthic and planktonic foraminiferal data occurs from the late Eocene to the late Oligocene. There is a correspondence of the 13C signal to deep-sea records; however, the amplitude of this change is greater than previously seen in deep-sea cores, possibly as a result of proximity to terrestrial sources of carbon. The covarying isotopic changes in both benthic and planktonic foraminifers suggest global causes, such as ice volume increases and increased terrestrial carbon input to the ocean. However, during the latter part of the record (early-late Oligocene), the increases in the benthic 18O without accompanying increases observed with planktonic foraminifers suggest that changes in only one part of the system occurred; one potential explanation being a decrease in bottom-water temperatures without concomitant changes in the surface waters. The 18O differences between species of planktonic foraminifers and the difference between planktonic and benthic 18O data indicate that diagenesis problems are minimal. These preliminary results are encouraging given that these cores are partially lithified.

Estimates of global sources and sinks of carbon from land cover and land use changes

The distribution of sources and sinks of carbon over the land surface is dominated by changes in land use such as deforestation, reforestation, and agricultural management. Despite, the importance of land-use change in dominating long-term net terrestrial fluxes of carbon, estimates of the annual flux are uncertain relative to other terms in the global carbon budget. The interaction of the nitrogen cycle via atmospheric N inputs and N limitation with the carbon cycle contributes to the uncertain effect of land use change on terrestrial carbon uptake. This study uses two different land use datasets to force the geographically explicit terrestrial carbon-nitrogen coupled component of the Integrated Science Assessment Model (ISAM) to examine the response of terrestrial carbon stocks to historical LCLUC (cropland, pastureland and wood harvest) while accounting for changes in N deposition, atmospheric CO2 and climate. One of the land use datasets is based on satellite data (SAGE) while the other uses population density maps (HYDE), which allows this study to investigate how global LCLUC data construction can affect model estimated emissions. The timeline chosen for this study starts before the Industrial Revolution in 1765 to the year 2000 because of the influence of rising population and economic development on regional LCLUC. Additionally, this study evaluates the impact that resulting secondary forests may have on terrestrial carbon uptake. The ISAM model simulations indicate that uncertainties in net terrestrial carbon fluxes during the 1990s are largely due to uncertainties in regional LCLUC data. Also results show that secondary forests increase the terrestrial carbon sink but secondary tropical forests carbon uptake are constrained due to nutrient limitation.

Marine and estuarine reservoir effects in central Queensland, Australia: Determination of Delta R values

As a component of archaeological investigations on the central Queensland coast, a series of five marine shell specimens live-collected between A.D. 1904 and A.D. 1929 and 11 shell/ charcoal paired samples from archaeological contexts were radiocarbon dated to determine local DeltaR values. The object of the study was to assess the potential influence of localized variation in marine reservoir effect in accurately determining the age of marine and estuarine shell from archaeological deposits in the area. Results indicate that the routinely applied DeltaR value of -5 +/- 35 for northeast Australia is erroneously calculated. The determined values suggest a minor revision to Reimer and Reimer's (2000) recommended value for northeast Australia from DeltaR = +11 +/- 5 to + 12 +/- 7, and specifically for central Queensland to DeltaR = +10 +/- 7, for near-shore open marine environments. In contrast, data obtained from estuarine shell/charcoal pairs demonstrate a general lack of consistency, suggesting estuary-specific patterns of variation in terrestrial carbon input and exchange with the open ocean. Preliminary data indicate that in some estuaries, at some time periods, a DeltaR value of more than - 155 +/- 55 may be appropriate, In estuarine contexts in central Queensland, a localized estuary-specific correction factor is recommended to account for geographical and temporal variation in C-14 activity. (C) 2002 Wiley Periodicals.

Modelling Amazonian forest eddy covariance data: a comparison of big leaf versus sun/shade models for the C-14 tower at Manaus I. Canopy photosynthesis

In this study, we concentrate on modelling gross primary productivity using two simple approaches to simulate canopy photosynthesis: "big leaf" and "sun/shade" models. Two approaches for calibration are used: scaling up of canopy photosynthetic parameters from the leaf to the canopy level and fitting canopy biochemistry to eddy covariance fluxes. Validation of the models is achieved by using eddy covariance data from the LBA site C14. Comparing the performance of both models we conclude that numerically (in terms of goodness of fit) and qualitatively, (in terms of residual response to different environmental variables) sun/shade does a better job. Compared to the sun/shade model, the big leaf model shows a lower goodness of fit and fails to respond to variations in the diffuse fraction, also having skewed responses to temperature and VPD. The separate treatment of sun and shade leaves in combination with the separation of the incoming light into direct beam and diffuse make sun/shade a strong modelling tool that catches more of the observed variability in canopy fluxes as measured by eddy covariance. In conclusion, the sun/shade approach is a relatively simple and effective tool for modelling photosynthetic carbon uptake that could be easily included in many terrestrial carbon models.

A geochemical and paleoceanographic investigation using marine molluscs of the late Quaternary marine submergences Quebec, Ontario, British Columbia

Sedjrrlents deposited in the Late Quaternary marine sUbrnergences that follov'ted the deglaciation of Ontario} Quebec., and 6ritlst-1 Columbia often contaln an abundant nlarlne invertebrate macrofauna. The rnacrofauna~ dotYllnated by aragonitic pelecypods} is fully preserved In their original mineralogy and cherrlistry 8S deternl1ned by x-ray dlffractlon., scannlng electron tl-,lcroscoDY., trace and r1l1 nor elet11ent analyses and stable isotopes. Ttle trace elernent and stable isotope geochen-Ilstry of chernlcal1y unaltered aragorlitlc molluscs can be used to determine paleoter1-lperatures and paleosallnltles." HO\Never} corrections need to be tllade \fvtlen deterrTIlnlng oxygen-isotope paleotenlperi:ttures due to the lnfluence of isotopically 11gtlt glaciol rneltv-laters and reduced sal1nltles. Ttle eastern Laurentide Ice Sheet probably had an o:~ygen lS0tOP1C composition as low as -8e) 0/00 (Sr1[IW). In additl0fl} corrections need to be rnade to the carbonlsotope values, before salinity deterrnlnatlons are t11ade., due to the reJjuctlon of the terrestrial carbon bl0rnass during glac1al maxlrna. Using geochernlcal data frot11 537 marlne n-'8crolnvertebrates frorTI 72 localities in soutt-,easter Ontarl0 and southern Quebec, it tras been deterrnined that the Late Quaternary Char1lplaln Sea \N6S density stratified along salinity and temperatlJre gradients. The deep-\h/aters of tt-,e Charnplaln Sea tlad salinities that ranged frorn 31 to 36 ppt} and terrlperatures of 00 to 5°C. Conversely.. the st1alloy./-\f*later regirrle of ttle Ctlarnplaln Sea tlad sal1nltles that ranged fron-, 24 to 33 ppt} Y.tltt1 terrlperatures ranglng from 5° to 15°C. Tr,8 rrlajorl rnlnor1 and trace e1et1-,ent geochernlcal analysls of 155 marine lnvertebrates frorn 4 10C611t1es of tt-,e Late Quaternary Ft. Langley Forrnatlon and Capl1ano Sedlments;. souttl\Nestern Brltlsh Columblal suggest l t~lat the 'waters of the o-,arlne lnundation that fol1o....ved the retreating Cordl11eran Ice Sheet had sal1nltles ranglng frorn 32 to 3f. DPt.