Search

Carbon–nitrogen (C–N) interactions regulate N availability for plant growth and for emissions of nitrous oxide (N2O) and the uptake of carbon dioxide. Future projections of these terrestrial greenhouse gas fluxes are strikingly divergent, leading to major uncertainties in projected global warming. Here we analyse the large increase in terrestrial N2O emissions over the past 21 000 years as reconstructed from ice-core isotopic data and presented in part 1 of this study. Remarkably, the increase occurred in two steps, each realized over decades and within a maximum of 2 centuries, at the onsets of the major deglacial Northern Hemisphere warming events. The data suggest a highly dynamic and responsive global N cycle. The increase may be explained by an increase in the flux of reactive N entering and leaving ecosystems or by an increase in N2O yield per unit N converted. We applied the LPX-Bern dynamic global vegetation model in deglacial simulations forced with Earth system model climate data to investigate N2O emission patterns, mechanisms, and C–N coupling. The N2O emission changes are mainly attributed to changes in temperature and precipitation and the loss of land due to sea-level rise. LPX-Bern simulates a deglacial increase in N2O emissions but underestimates the reconstructed increase by 47 %. Assuming time-independent N sources in the model to mimic progressive N limitation of plant growth results in a decrease in N2O emissions in contrast to the reconstruction. Our results appear consistent with suggestions of (a) biological controls on ecosystem N acquisition and (b) flexibility in the coupling of the C and N cycles during periods of rapid environmental change. A dominant uncertainty in the explanation of the reconstructed N2O emissions is the poorly known N2O yield per N lost through gaseous pathways and its sensitivity to soil conditions. The deglacial N2O record provides a constraint for future studies. more...

Using high-precision and centennial-resolution ice core information on atmospheric nitrous oxide concentrations and its stable nitrogen and oxygen isotopic composition, we quantitatively reconstruct changes in the terrestrial and marine N2O emissions over the last 21 000 years. Our reconstruction indicates that N2O emissions from land and ocean increased over the deglaciation largely in parallel by 1.7±0.3 and 0.7±0.3 TgN yr−1, respectively, relative to the Last Glacial Maximum level. However, during the abrupt Northern Hemisphere warmings at the onset of the Bølling–Allerød warming and the end of the Younger Dryas, terrestrial emissions respond more rapidly to the northward shift in the Intertropical Convergence Zone connected to the resumption of the Atlantic Meridional Overturning Circulation. About 90 % of these large step increases were realized within 2 centuries at maximum. In contrast, marine emissions start to slowly increase already many centuries before the rapid warmings, possibly connected to a re-equilibration of subsurface oxygen in response to previous changes. Marine emissions decreased, concomitantly with changes in atmospheric CO2 and δ13C(CO2), at the onset of the termination and remained minimal during the early phase of Heinrich Stadial 1. During the early Holocene a slow decline in marine N2O emission of 0.4 TgN yr−1 is reconstructed, which suggests an improvement of subsurface water ventilation in line with slowly increasing Atlantic overturning circulation. In the second half of the Holocene total emissions remain on a relatively constant level, but with significant millennial variability. The latter is still difficult to attribute to marine or terrestrial sources. Our N2O emission records provide important quantitative benchmarks for ocean and terrestrial nitrogen cycle models to study the influence of climate on nitrogen turnover on timescales from several decades to glacial–interglacial changes. more...

Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH4) budget over 2000–2012 (Saunois et al., 2016), we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000–2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000–2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008–2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16–32] Tg CH4 yr−1 higher methane emissions over the period 2008–2012 compared to 2002–2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002–2006 and 2008–2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13CH4. However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13CH4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper. more...

We present a global distribution of surface methane (CH4) emission estimates for 2000–2012 derived using the CarbonTracker Europe-CH4 (CTE-CH4) data assimilation system. In CTE-CH4, anthropogenic and biospheric CH4 emissions are simultaneously estimated based on constraints of global atmospheric in situ CH4 observations. The system was configured to either estimate only anthropogenic or biospheric sources per region, or to estimate both categories simultaneously. The latter increased the number of optimizable parameters from 62 to 78. In addition, the differences between two numerical schemes available to perform turbulent vertical mixing in the atmospheric transport model TM5 were examined. Together, the system configurations encompass important axes of uncertainty in inversions and allow us to examine the robustness of the flux estimates. The posterior emission estimates are further evaluated by comparing simulated atmospheric CH4 to surface in situ observations, vertical profiles of CH4 made by aircraft, remotely sensed dry-air total column-averaged mole fraction (XCH4) from the Total Carbon Column Observing Network (TCCON), and XCH4 from the Greenhouse gases Observing Satellite (GOSAT). The evaluation with non-assimilated observations shows that posterior XCH4 is better matched with the retrievals when the vertical mixing scheme with faster interhemispheric exchange is used. Estimated posterior mean total global emissions during 2000–2012 are 516 ± 51 Tg CH4 yr−1, with an increase of 18 Tg CH4 yr−1 from 2000–2006 to 2007–2012. The increase is mainly driven by an increase in emissions from South American temperate, Asian temperate and Asian tropical TransCom regions. In addition, the increase is hardly sensitive to different model configurations ( 4 yr−1 difference), and much smaller than suggested by EDGAR v4.2 FT2010 inventory (33 Tg CH4 yr−1), which was used for prior anthropogenic emission estimates. The result is in good agreement with other published estimates from inverse modelling studies (16–20 Tg CH4 yr−1). However, this study could not conclusively separate a small trend in biospheric emissions (−5 to +6.9 Tg CH4 yr−1) from the much larger trend in anthropogenic emissions (15–27 Tg CH4 yr−1). Finally, we find that the global and North American CH4 balance could be closed over this time period without the previously suggested need to strongly increase anthropogenic CH4 emissions in the United States. With further developments, especially on the treatment of the atmospheric CH4 sink, we expect the data assimilation system presented here will be able to contribute to the ongoing interpretation of changes in this important greenhouse gas budget. more...

The global methane (CH4) budget is becoming an increasingly important component for managing realistic pathways to mitigate climate change. This relevance, due to a shorter atmospheric lifetime and a stronger warming potential than carbon dioxide, is challenged by the still unexplained changes of atmospheric CH4 over the past decade. Emissions and concentrations of CH4 are continuing to increase, making CH4 the second most important human-induced greenhouse gas after carbon dioxide. Two major difficulties in reducing uncertainties come from the large variety of diffusive CH4 sources that overlap geographically, and from the destruction of CH4 by the very short-lived hydroxyl radical (OH). To address these difficulties, we have established a consortium of multi-disciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate research on the methane cycle, and producing regular (∼ biennial) updates of the global methane budget. This consortium includes atmospheric physicists and chemists, biogeochemists of surface and marine emissions, and socio-economists who study anthropogenic emissions. Following Kirschke et al. (2013), we propose here the first version of a living review paper that integrates results of top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models, inventories and data-driven approaches (including process-based models for estimating land surface emissions and atmospheric chemistry, and inventories for anthropogenic emissions, data-driven extrapolations). For the 2003–2012 decade, global methane emissions are estimated by top-down inversions at 558 Tg CH4 yr−1, range 540–568. About 60 % of global emissions are anthropogenic (range 50–65 %). Since 2010, the bottom-up global emission inventories have been closer to methane emissions in the most carbon-intensive Representative Concentrations Pathway (RCP8.5) and higher than all other RCP scenarios. Bottom-up approaches suggest larger global emissions (736 Tg CH4 yr−1, range 596–884) mostly because of larger natural emissions from individual sources such as inland waters, natural wetlands and geological sources. Considering the atmospheric constraints on the top-down budget, it is likely that some of the individual emissions reported by the bottom-up approaches are overestimated, leading to too large global emissions. Latitudinal data from top-down emissions indicate a predominance of tropical emissions (∼ 64 % of the global budget, 4 yr−1, range 51–72, −14 %) and higher emissions in Africa (86 Tg CH4 yr−1, range 73–108, +19 %) than bottom-up values used as prior estimates. Overall, uncertainties for anthropogenic emissions appear smaller than those from natural sources, and the uncertainties on source categories appear larger for top-down inversions than for bottom-up inventories and models. The most important source of uncertainty on the methane budget is attributable to emissions from wetland and other inland waters. We show that the wetland extent could contribute 30–40 % on the estimated range for wetland emissions. Other priorities for improving the methane budget include the following: (i) the development of process-based models for inland-water emissions, (ii) the intensification of methane observations at local scale (flux measurements) to constrain bottom-up land surface models, and at regional scale (surface networks and satellites) to constrain top-down inversions, (iii) improvements in the estimation of atmospheric loss by OH, and (iv) improvements of the transport models integrated in top-down inversions. The data presented here can be downloaded from the Carbon Dioxide Information Analysis Center (http://doi.org/10.3334/CDIAC/GLOBAL_METHANE_BUDGET_2016_V1.1) and the Global Carbon Project. more...

Records of stable oxygen isotope ratios in tree rings are valuable tools to reconstruct past climatic conditions and investigate the response of trees to those conditions. So far the use of stable oxygen isotope signatures of tree rings has not been systematically evaluated in dynamic global vegetation models (DGVMs). DGVMs integrate many hydrological and physiological processes and their application could improve proxy-model comparisons and the interpretation of oxygen isotope records. Here we present an approach to simulate leaf water and stem cellulose δ18O of trees using the LPX-Bern DGVM (LPX-Bern). Our results lie within a few per mil of measured tree ring δ18O of 31 different forest stands mainly located in Europe. Temporal means over the last 5 decades as well as interannual variations for a subset of sites in Switzerland are captured. A sensitivity analysis reveals that relative humidity, temperature, and the water isotope boundary conditions have the largest influence on simulated stem cellulose δ18O, followed by all climatic factors combined, whereas increasing atmospheric CO2 and nitrogen deposition exert no impact. We conclude that simulations with LPX-Bern are useful for investigating large-scale oxygen isotope patterns of tree ring cellulose to elucidate the importance of different environmental factors on isotope variations and therefore help to reduce uncertainties in the interpretation of δ18O of tree rings. more...

Eight surface observation sites providing quasi-continuous measurements of atmospheric methane mixing ratios have been operated since the mid-2000's in Siberia. For the first time in a single work, we assimilate 1 year of these in situ observations in an atmospheric inversion. Our objective is to quantify methane surface fluxes from anthropogenic and wetland sources at the mesoscale in the Siberian lowlands for the year 2010. To do so, we first inquire about the way the inversion uses the observations and the way the fluxes are constrained by the observation sites. As atmospheric inversions at the mesoscale suffer from mis-quantified sources of uncertainties, we follow recent innovations in inversion techniques and use a new inversion approach which quantifies the uncertainties more objectively than the previous inversion systems. We find that, due to errors in the representation of the atmospheric transport and redundant pieces of information, only one observation every few days is found valuable by the inversion. The remaining high-resolution quasi-continuous signal is representative of very local emission patterns difficult to analyse with a mesoscale system. An analysis of the use of information by the inversion also reveals that the observation sites constrain methane emissions within a radius of 500 km. More observation sites than the ones currently in operation are then necessary to constrain the whole Siberian lowlands. Still, the fluxes within the constrained areas are quantified with objectified uncertainties. Finally, the tolerance intervals for posterior methane fluxes are of roughly 20 % (resp. 50 %) of the fluxes for anthropogenic (resp. wetland) sources. About 50–70 % of Siberian lowlands emissions are constrained by the inversion on average on an annual basis. Extrapolating the figures on the constrained areas to the whole Siberian lowlands, we find a regional methane budget of 5–28 TgCH4 for the year 2010, i.e. 1–5 % of the global methane emissions. As very few in situ observations are available in the region of interest, observations of methane total columns from the Greenhouse Gas Observing SATellite (GOSAT) are tentatively used for the evaluation of the inversion results, but they exhibit only a marginal signal from the fluxes within the region of interest. more...

An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (δD-δ18Oice, δ18Oatm, total air content, CO2, CH4, N2O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the "basal clean ice facies". Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies. It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate "from within", and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice–bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future "oldest ice" drilling location in Antarctica. more...

Wetlands are the world's largest natural source of methane, a powerful greenhouse gas. The strong sensitivity of methane emissions to environmental factors such as soil temperature and moisture has led to concerns about potential positive feedbacks to climate change. This risk is particularly relevant at high latitudes, which have experienced pronounced warming and where thawing permafrost could potentially liberate large amounts of labile carbon over the next 100 years. However, global models disagree as to the magnitude and spatial distribution of emissions, due to uncertainties in wetland area and emissions per unit area and a scarcity of in situ observations. Recent intensive field campaigns across the West Siberian Lowland (WSL) make this an ideal region over which to assess the performance of large-scale process-based wetland models in a high-latitude environment. Here we present the results of a follow-up to the Wetland and Wetland CH4 Intercomparison of Models Project (WETCHIMP), focused on the West Siberian Lowland (WETCHIMP-WSL). We assessed 21 models and 5 inversions over this domain in terms of total CH4 emissions, simulated wetland areas, and CH4 fluxes per unit wetland area and compared these results to an intensive in situ CH4 flux data set, several wetland maps, and two satellite surface water products. We found that (a) despite the large scatter of individual estimates, 12-year mean estimates of annual total emissions over the WSL from forward models (5.34 ± 0.54 Tg CH4 yr−1), inversions (6.06 ± 1.22 Tg CH4 yr−1), and in situ observations (3.91 ± 1.29 Tg CH4 yr−1) largely agreed; (b) forward models using surface water products alone to estimate wetland areas suffered from severe biases in CH4 emissions; (c) the interannual time series of models that lacked either soil thermal physics appropriate to the high latitudes or realistic emissions from unsaturated peatlands tended to be dominated by a single environmental driver (inundation or air temperature), unlike those of inversions and more sophisticated forward models; (d) differences in biogeochemical schemes across models had relatively smaller influence over performance; and (e) multiyear or multidecade observational records are crucial for evaluating models' responses to long-term climate change. more...

Simulating the spatio-temporal dynamics of inundation is key to understanding the role of wetlands under past and future climate change. Earlier modelling studies have mostly relied on fixed prescribed peatland maps and inundation time series of limited temporal coverage. Here, we describe and assess the the Dynamical Peatland Model Based on TOPMODEL (DYPTOP), which predicts the extent of inundation based on a computationally efficient TOPMODEL implementation. This approach rests on an empirical, grid-cell-specific relationship between the mean soil water balance and the flooded area. DYPTOP combines the simulated inundation extent and its temporal persistency with criteria for the ecosystem water balance and the modelled peatland-specific soil carbon balance to predict the global distribution of peatlands. We apply DYPTOP in combination with the LPX-Bern DGVM and benchmark the global-scale distribution, extent, and seasonality of inundation against satellite data. DYPTOP successfully predicts the spatial distribution and extent of wetlands and major boreal and tropical peatland complexes and reveals the governing limitations to peatland occurrence across the globe. Peatlands covering large boreal lowlands are reproduced only when accounting for a positive feedback induced by the enhanced mean soil water holding capacity in peatland-dominated regions. DYPTOP is designed to minimize input data requirements, optimizes computational efficiency and allows for a modular adoption in Earth system models. more...

During the last glacial cycle, Greenland temperature showed many rapid temperature variations, the so-called Dansgaard–Oeschger (DO) events. The past atmospheric methane concentration closely followed these temperature variations, which implies that the warmings recorded in Greenland were probably hemispheric in extent. Here we substantially extend and complete the North Greenland Ice Core Project (NGRIP) methane record from the Preboreal Holocene (PB) back to the end of the last interglacial period with a mean time resolution of 54 yr. We relate the amplitudes of the methane increases associated with DO events to the amplitudes of the local Greenland NGRIP temperature increases derived from stable nitrogen isotope (δ15N) measurements, which have been performed along the same ice core (Kindler et al., 2014). We find the ratio to oscillate between 5 parts per billion (ppb) per °C and 18 ppb °C−1 with the approximate frequency of the precessional cycle. A remarkably high ratio of 25.5 ppb °C−1 is reached during the transition from the Younger Dryas (YD) to the PB. Analysis of the timing of the fast methane and temperature increases reveals significant lags of the methane increases relative to NGRIP temperature for DO events 5, 9, 10, 11, 13, 15, 19, and 20. These events generally have small methane increase rates and we hypothesize that the lag is caused by pronounced northward displacement of the source regions from stadial to interstadial. We further show that the relative interpolar concentration difference (rIPD) of methane is about 4.5% for the stadials between DO events 18 and 20, which is in the same order as in the stadials before and after DO event 2 around the Last Glacial Maximum. The rIPD of methane remains relatively stable throughout the full last glacial, with a tendency for elevated values during interstadial compared to stadial periods. more...

Tropical wetlands are estimated to represent about 50% of the natural wetland methane (CH4) emissions and explain a large fraction of the observed CH4 variability on timescales ranging from glacial–interglacial cycles to the currently observed year-to-year variability. Despite their importance, however, tropical wetlands are poorly represented in global models aiming to predict global CH4 emissions. This publication documents a first step in the development of a process-based model of CH4 emissions from tropical floodplains for global applications. For this purpose, the LPX-Bern Dynamic Global Vegetation Model (LPX hereafter) was slightly modified to represent floodplain hydrology, vegetation and associated CH4 emissions. The extent of tropical floodplains was prescribed using output from the spatially explicit hydrology model PCR-GLOBWB. We introduced new plant functional types (PFTs) that explicitly represent floodplain vegetation. The PFT parameterizations were evaluated against available remote-sensing data sets (GLC2000 land cover and MODIS Net Primary Productivity). Simulated CH4 flux densities were evaluated against field observations and regional flux inventories. Simulated CH4 emissions at Amazon Basin scale were compared to model simulations performed in the WETCHIMP intercomparison project. We found that LPX reproduces the average magnitude of observed net CH4 flux densities for the Amazon Basin. However, the model does not reproduce the variability between sites or between years within a site. Unfortunately, site information is too limited to attest or disprove some model features. At the Amazon Basin scale, our results underline the large uncertainty in the magnitude of wetland CH4 emissions. Sensitivity analyses gave insights into the main drivers of floodplain CH4 emission and their associated uncertainties. In particular, uncertainties in floodplain extent (i.e., difference between GLC2000 and PCR-GLOBWB output) modulate the simulated emissions by a factor of about 2. Our best estimates, using PCR-GLOBWB in combination with GLC2000, lead to simulated Amazon-integrated emissions of 44.4 ± 4.8 Tg yr−1. Additionally, the LPX emissions are highly sensitive to vegetation distribution. Two simulations with the same mean PFT cover, but different spatial distributions of grasslands within the basin, modulated emissions by about 20%. Correcting the LPX-simulated NPP using MODIS reduces the Amazon emissions by 11.3%. Finally, due to an intrinsic limitation of LPX to account for seasonality in floodplain extent, the model failed to reproduce the full dynamics in CH4 emissions but we proposed solutions to this issue. The interannual variability (IAV) of the emissions increases by 90% if the IAV in floodplain extent is accounted for, but still remains lower than in most of the WETCHIMP models. While our model includes more mechanisms specific to tropical floodplains, we were unable to reduce the uncertainty in the magnitude of wetland CH4 emissions of the Amazon Basin. Our results helped identify and prioritize directions towards more accurate estimates of tropical CH4 emissions, and they stress the need for more research to constrain floodplain CH4 emissions and their temporal variability, even before including other fundamental mechanisms such as floating macrophytes or lateral water fluxes. more...

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways. more...

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions. more...

The Wetland and Wetland CH4 Intercomparison of Models Project (WETCHIMP) was created to evaluate our present ability to simulate large-scale wetland characteristics and corresponding methane (CH4) emissions. A multi-model comparison is essential to evaluate the key uncertainties in the mechanisms and parameters leading to methane emissions. Ten modelling groups joined WETCHIMP to run eight global and two regional models with a common experimental protocol using the same climate and atmospheric carbon dioxide (CO2) forcing datasets. We reported the main conclusions from the intercomparison effort in a companion paper (Melton et al., 2013). Here we provide technical details for the six experiments, which included an equilibrium, a transient, and an optimized run plus three sensitivity experiments (temperature, precipitation, and atmospheric CO2 concentration). The diversity of approaches used by the models is summarized through a series of conceptual figures, and is used to evaluate the wide range of wetland extent and CH4 fluxes predicted by the models in the equilibrium run. We discuss relationships among the various approaches and patterns in consistencies of these model predictions. Within this group of models, there are three broad classes of methods used to estimate wetland extent: prescribed based on wetland distribution maps, prognostic relationships between hydrological states based on satellite observations, and explicit hydrological mass balances. A larger variety of approaches was used to estimate the net CH4 fluxes from wetland systems. Even though modelling of wetland extent and CH4 emissions has progressed significantly over recent decades, large uncertainties still exist when estimating CH4 emissions: there is little consensus on model structure or complexity due to knowledge gaps, different aims of the models, and the range of temporal and spatial resolutions of the models. more...

Rapid changes in atmospheric methane (CH4), temperature and precipitation are documented by Greenland ice core data both for glacial times (the so called Dansgaard-Oeschger (D-O) events) as well as for a cooling event in the early Holocene (the 8.2 kyr event). The onsets of D-O warm events are paralleled by abrupt increases in CH4 by up to 250 ppb in a few decades. Vice versa, the 8.2 kyr event is accompanied by an intermittent decrease in CH4 of about 80 ppb over 150 yr. The abrupt CH4 changes are thought to mainly originate from source emission variations in tropical and boreal wet ecosystems, but complex process oriented bottom-up model estimates of the changes in these ecosystems during rapid climate changes are still missing. Here we present simulations of CH4 emissions from northern peatlands with the LPJ-Bern dynamic global vegetation model. The model represents CH4 production and oxidation in soils and transport by ebullition, through plant aerenchyma, and by diffusion. Parameters are tuned to represent site emission data as well as inversion-based estimates of northern wetland emissions. The model is forced with climate input data from freshwater hosing experiments using the NCAR CSM1.4 climate model to simulate an abrupt cooling event. A concentration reduction of ~10 ppb is simulated per degree K change of mean northern hemispheric surface temperature in peatlands. Peatland emissions are equally sensitive to both changes in temperature and in precipitation. If simulated changes are taken as an analogy to the 8.2 kyr event, boreal peatland emissions alone could only explain 23% of the 80 ppb decline in atmospheric methane concentration. This points to a significant contribution to source changes from low latitude and tropical wetlands to this event. more...

Global wetlands are believed to be climate sensitive, and are the largest natural emitters of methane (CH4). Increased wetland CH4 emissions could act as a positive feedback to future warming. The Wetland and Wetland CH4 Inter-comparison of Models Project (WETCHIMP) investigated our present ability to simulate large-scale wetland characteristics and corresponding CH4 emissions. To ensure inter-comparability, we used a common experimental protocol driving all models with the same climate and carbon dioxide (CO2) forcing datasets. The WETCHIMP experiments were conducted for model equilibrium states as well as transient simulations covering the last century. Sensitivity experiments investigated model response to changes in selected forcing inputs (precipitation, temperature, and atmospheric CO2 concentration). Ten models participated, covering the spectrum from simple to relatively complex, including models tailored either for regional or global simulations. The models also varied in methods to calculate wetland size and location, with some models simulating wetland area prognostically, while other models relied on remotely sensed inundation datasets, or an approach intermediate between the two. Four major conclusions emerged from the project. First, the suite of models demonstrate extensive disagreement in their simulations of wetland areal extent and CH4 emissions, in both space and time. Simple metrics of wetland area, such as the latitudinal gradient, show large variability, principally between models that use inundation dataset information and those that independently determine wetland area. Agreement between the models improves for zonally summed CH4 emissions, but large variation between the models remains. For annual global CH4 emissions, the models vary by ±40% of the all-model mean (190 Tg CH4 yr−1). Second, all models show a strong positive response to increased atmospheric CO2 concentrations (857 ppm) in both CH4 emissions and wetland area. In response to increasing global temperatures (+3.4 °C globally spatially uniform), on average, the models decreased wetland area and CH4 fluxes, primarily in the tropics, but the magnitude and sign of the response varied greatly. Models were least sensitive to increased global precipitation (+3.9 % globally spatially uniform) with a consistent small positive response in CH4 fluxes and wetland area. Results from the 20th century transient simulation show that interactions between climate forcings could have strong non-linear effects. Third, we presently do not have sufficient wetland methane observation datasets adequate to evaluate model fluxes at a spatial scale comparable to model grid cells (commonly 0.5°). This limitation severely restricts our ability to model global wetland CH4 emissions with confidence. Our simulated wetland extents are also difficult to evaluate due to extensive disagreements between wetland mapping and remotely sensed inundation datasets. Fourth, the large range in predicted CH4 emission rates leads to the conclusion that there is both substantial parameter and structural uncertainty in large-scale CH4 emission models, even after uncertainties in wetland areas are accounted for. more...

Reconstructions of past atmospheric methane concentrations are available from ice cores from both Greenland and Antarctica. The difference observed between the two polar methane concentration levels represents a valuable constraint on the geographical location of the methane sources. Here we present new high-resolution methane records from the North Greenland Ice Core Project (NGRIP) and the European Project for Ice Coring in Antarctica (EPICA) Dronning Maud Land (EDML) ice cores covering Termination 1, the Last Glacial Maximum, and parts of the last glacial back to 32 000 years before present. Due to the high resolution of the records, the synchronisation between the ice cores from NGRIP and EDML is considerably improved, and the interpolar concentration difference of methane is determined with unprecedented precision and temporal resolution. Relative to the mean methane concentration, we find a rather stable positive relative interpolar difference throughout the record with its minimum value of 3.7 ± 0.7 % between 21 900–21 200 years before present, which is higher than previously estimated in this interval close to the Last Glacial Maximum. This implies that Northern Hemisphere boreal wetland sources were never completely shut off during the peak glacial, as suggested from previous bipolar methane concentration records. Starting at 21 000 years before present, i.e. several millennia prior to the transition into the Holocene, the relative interpolar difference becomes even more positive and stays at a fairly stable level of 6.5 ± 0.8 % during Termination 1. We thus find that the boreal and tropical methane sources increased by approximately the same factor during Termination 1. We hypothesise that latitudinal shifts in the Intertropical Convergence Zone (ITCZ) and the monsoon system contribute, either by dislocation of the methane source regions or, in case of the ITCZ, also by changing the relative atmospheric volumes of the Northern and Southern Hemispheres, to the subtle variations in the relative interpolar concentration difference on glacial/interglacial as well as on millennial time scales. more...

Natural methane (CH4) emissions from wet ecosystems are an important part of today's global CH4 budget. Climate affects the exchange of CH4 between ecosystems and the atmosphere by influencing CH4 production, oxidation, and transport in the soil. The net CH4 exchange depends on ecosystem hydrology, soil and vegetation characteristics. Here, the LPJ-WHyMe global dynamical vegetation model is used to simulate global net CH4 emissions for different ecosystems: northern peatlands (45°–90° N), naturally inundated wetlands (60° S–45° N), rice agriculture and wet mineral soils. Mineral soils are a potential CH4 sink, but can also be a source with the direction of the net exchange depending on soil moisture content. The geographical and seasonal distributions are evaluated against multi-dimensional atmospheric inversions for 2003–2005, using two independent four-dimensional variational assimilation systems. The atmospheric inversions are constrained by the atmospheric CH4 observations of the SCIAMACHY satellite instrument and global surface networks. Compared to LPJ-WHyMe the inversions result in a~significant reduction in the emissions from northern peatlands and suggest that LPJ-WHyMe maximum annual emissions peak about one month late. The inversions do not put strong constraints on the division of sources between inundated wetlands and wet mineral soils in the tropics. Based on the inversion results we diagnose model parameters in LPJ-WHyMe and simulate the surface exchange of CH4 over the period 1990–2008. Over the whole period we infer an increase of global ecosystem CH4 emissions of +1.11 Tg CH4 yr−1, not considering potential additional changes in wetland extent. The increase in simulated CH4 emissions is attributed to enhanced soil respiration resulting from the observed rise in land temperature and in atmospheric carbon dioxide that were used as input. The long-term decline of the atmospheric CH4 growth rate from 1990 to 2006 cannot be fully explained with the simulated ecosystem emissions. However, these emissions show an increasing trend of +3.62 Tg CH4 yr−1 over 2005–2008 which can partly explain the renewed increase in atmospheric CH4 concentration during recent years. more...

Millennial variability is a robust feature of many paleoclimate records, at least throughout the last several glacial cycles. Here we use the mean signal from Antarctic climate events 1 to 4 to probe the EPICA Dome C temperature proxy reconstruction through the last 500 ka for similar millennial-scale events. We find that clusters of millennial events occurred in a regular fashion over half of the time during this with a mean recurrence interval of 21 kyr. We find that there is no consistent link between ice-rafted debris deposition and millennial variability. Instead we speculate that changes in the zonality of atmospheric circulation over the North Atlantic form a viable alternative to freshwater release from icebergs as a trigger for millennial variability. We suggest that millennial changes in the zonality of atmospheric circulation over the North Atlantic are linked to precession via sea-ice feedbacks and that this relationship is modified by the presence of the large, Northern Hemisphere ice sheets during glacial periods. more...

Gas is trapped in polar ice sheets at ~50–120 m below the surface and is therefore younger than the surrounding ice. Firn densification models are used to evaluate this ice age-gas age difference (Δage) in the past. However, such models need to be validated by data, in particular for periods colder than present day on the East Antarctic plateau. Here we bring new constraints to test a firn densification model applied to the EPICA Dome C (EDC) site for the last 50 kyr, by linking the EDC ice core to the EPICA Dronning Maud Land (EDML) ice core, both in the ice phase (using volcanic horizons) and in the gas phase (using rapid methane variations). We also use the structured 10Be peak, occurring 41 kyr before present (BP) and due to the low geomagnetic field associated with the Laschamp event, to experimentally estimate the Δage during this event. Our results seem to reveal an overestimate of the Δage by the firn densification model during the last glacial period at EDC. Tests with different accumulation rates and temperature scenarios do not entirely resolve this discrepancy. Although the exact reasons for the Δage overestimate at the two EPICA sites remain unknown at this stage, we conclude that current densification model simulations have deficits under glacial climatic conditions. Whatever the cause of the Δage overestimate, our finding suggests that the phase relationship between CO2 and EDC temperature previously inferred for the start of the last deglaciation (lag of CO2 by 800±600 yr) seems to be overestimated. more...

The EPICA (European Project for Ice Coring in Antarctica) Dome C drilling in East Antarctica has now been completed to a depth of 3260 m, at only a few meters above bedrock. Here we present the new EDC3 chronology, which is based on the use of 1) a snow accumulation and mechanical flow model, and 2) a set of independent age markers along the core. These are obtained by pattern matching of recorded parameters to either absolutely dated paleoclimatic records, or to insolation variations. We show that this new time scale is in excellent agreement with the Dome Fuji and Vostok ice core time scales back to 100 kyr within 1 kyr. Discrepancies larger than 3 kyr arise during MIS 5.4, 5.5 and 6, which points to anomalies in either snow accumulation or mechanical flow during these time periods. We estimate that EDC3 gives accurate event durations within 20% (2σ) back to MIS11 and accurate absolute ages with a maximum uncertainty of 6 kyr back to 800 kyr. more...

To interpret new high resolution climate records it becomes more and more important to know about the succession of climate events. Such knowledge is hard to get especially when dealing with different types of climate archives. Even for ice cores a direct synchronization between ice cores from Greenland and Antarctica has not been possible so far due to the lack of time markers occurring in both hemispheres. Fortunately, variations in the time series of global gas records can be used as indirect time markers. Here we discuss in detail the steps that are necessary to synchronize ice cores via global gas records exemplified on the synchronization of the EPICA ice core from Dronning Maud Land to a Greenland record from North GRIP. more...

While there are no indications of mixing back to 800 000 years in the EPICA Dome C ice core record, comparison with marine sediment records shows significant differences in the timing and duration of events prior to stage 11 (~430 ka, thousands of years before 1950). A relationship between the isotopic composition of atmospheric oxygen (δ18O of O2, noted δ18Oatm) and daily northern hemisphere summer insolation has been observed for the youngest four climate cycles. Here we use this relationship with new δ18O of O2 measurements to show that anomalous flow in the bottom 500 m of the core distorts the duration of events by up to a factor of 2. By tuning δ18Oatm to orbital precession we derive a corrected thinning function and present a revised age scale for the interval corresponding to Marine Isotope Stages 11–20 in the EPICA Dome C ice core. Uncertainty in the phasing of δ18Oatm with respect to insolation variations in the precession band limits the accuracy of this new agescale to ±6 kyr (thousand of years). The previously reported ~30 kyr duration of interglacial stage 11 is unchanged. In contrast, the duration of stage 15.1 is reduced by a factor of 2, from 31 to 16 kyr. more...

No abstract available.

Information on how species distributions and ecosystem services are impacted by anthropogenic climate change is important for adaptation planning. Palaeo data suggest that Abies alba formed forests under significantly warmer-than-present conditions in Europe and might be a native substitute for widespread drought-sensitive temperate and boreal tree species such as beech (Fagus sylvatica) and spruce (Picea abies) under future global warming conditions. Here, we combine pollen and macrofossil data, modern observations, and results from transient simulations with the LPX-Bern dynamic global vegetation model to assess past and future distributions of A. alba in Europe. LPX-Bern is forced with climate anomalies from a run over the past 21 000 years with the Community Earth System Model, modern climatology, and with 21st-century multimodel ensemble results for the high-emission RCP8.5 and the stringent mitigation RCP2.6 pathway. The simulated distribution for present climate encompasses the modern range of A. alba, with the model exceeding the present distribution in north-western and southern Europe. Mid-Holocene pollen data and model results agree for southern Europe, suggesting that at present, human impacts suppress the distribution in southern Europe. Pollen and model results both show range expansion starting during the Bølling-Allerød warm period, interrupted by the Younger Dryas cold, and resuming during the Holocene. The distribution of A. alba expands to the north-east in all future scenarios, whereas the potential (currently unrealized) range would be substantially reduced in southern Europe under RCP8.5. A. alba maintains its current range in central Europe despite competition by other thermophilous tree species. Our combined palaeoecological and model evidence suggest that A. alba may ensure important ecosystem services including stand and slope stability, infrastructure protection, and carbon sequestration under significantly warmer-than-present conditions in central Europe., (© 2015 John Wiley & Sons Ltd.) more...

Nitrous oxide (N2O) is an important greenhouse gas and ozone-depleting substance that has anthropogenic as well as natural marine and terrestrial sources. The tropospheric N2O concentrations have varied substantially in the past in concert with changing climate on glacial-interglacial and millennial timescales. It is not well understood, however, how N2O emissions from marine and terrestrial sources change in response to varying environmental conditions. The distinct isotopic compositions of marine and terrestrial N2O sources can help disentangle the relative changes in marine and terrestrial N2O emissions during past climate variations. Here we present N2O concentration and isotopic data for the last deglaciation, from 16,000 to 10,000 years before present, retrieved from air bubbles trapped in polar ice at Taylor Glacier, Antarctica. With the help of our data and a box model of the N2O cycle, we find a 30 per cent increase in total N2O emissions from the late glacial to the interglacial, with terrestrial and marine emissions contributing equally to the overall increase and generally evolving in parallel over the last deglaciation, even though there is no a priori connection between the drivers of the two sources. However, we find that terrestrial emissions dominated on centennial timescales, consistent with a state-of-the-art dynamic global vegetation and land surface process model that suggests that during the last deglaciation emission changes were strongly influenced by temperature and precipitation patterns over land surfaces. The results improve our understanding of the drivers of natural N2O emissions and are consistent with the idea that natural N2O emissions will probably increase in response to anthropogenic warming. more...

The increasing carbon dioxide (CO2 ) concentration in the atmosphere in combination with climatic changes throughout the last century are likely to have had a profound effect on the physiology of trees: altering the carbon and water fluxes passing through the stomatal pores. However, the magnitude and spatial patterns of such changes in natural forests remain highly uncertain. Here, stable carbon isotope ratios from a network of 35 tree-ring sites located across Europe are investigated to determine the intrinsic water-use efficiency (iWUE), the ratio of photosynthesis to stomatal conductance from 1901 to 2000. The results were compared with simulations of a dynamic vegetation model (LPX-Bern 1.0) that integrates numerous ecosystem and land-atmosphere exchange processes in a theoretical framework. The spatial pattern of tree-ring derived iWUE of the investigated coniferous and deciduous species and the model results agreed significantly with a clear south-to-north gradient, as well as a general increase in iWUE over the 20th century. The magnitude of the iWUE increase was not spatially uniform, with the strongest increase observed and modelled for temperate forests in Central Europe, a region where summer soil-water availability decreased over the last century. We were able to demonstrate that the combined effects of increasing CO2 and climate change leading to soil drying have resulted in an accelerated increase in iWUE. These findings will help to reduce uncertainties in the land surface schemes of global climate models, where vegetation-climate feedbacks are currently still poorly constrained by observational data., (© 2014 John Wiley & Sons Ltd.) more...

Ecosystem nitrous oxide (N2O) emissions respond to changes in climate and CO2 concentration as well as anthropogenic nitrogen (N) enhancements. Here, we aimed to quantify the responses of natural ecosystem N2O emissions to multiple environmental drivers using a process-based global vegetation model (DyN-LPJ). We checked that modelled annual N2O emissions from nonagricultural ecosystems could reproduce field measurements worldwide, and experimentally observed responses to step changes in environmental factors. We then simulated global N2O emissions throughout the 20th century and analysed the effects of environmental changes. The model reproduced well the global pattern of N2O emissions and the observed responses of N cycle components to changes in environmental factors. Simulated 20th century global decadal-average soil emissions were c. 8.2-9.5 Tg N yr(-1) (or 8.3-10.3 Tg N yr(-1) with N deposition). Warming and N deposition contributed 0.85±0.41 and 0.80±0.14 Tg N yr(-1), respectively, to an overall upward trend. Rising CO2 also contributed, in part, through a positive interaction with warming. The modelled temperature dependence of N2O emission (c. 1 Tg N yr(-1) K(-1)) implies a positive climate feedback which, over the lifetime of N2O (114 yr), could become as important as the climate-carbon cycle feedback caused by soil CO2 release., (© 2012 The Authors. New Phytologist © 2012 New Phytologist Trust.) more...

The causes of past changes in the global methane cycle and especially the role of marine methane hydrate (clathrate) destabilization events are a matter of debate. Here we present evidence from the North Greenland Ice Core Project ice core based on the hydrogen isotopic composition of methane [deltaD(CH4)] that clathrates did not cause atmospheric methane concentration to rise at the onset of Dansgaard-Oeschger (DO) events 7 and 8. Box modeling supports boreal wetland emissions as the most likely explanation for the interstadial increase. Moreover, our data show that deltaD(CH4) dropped 500 years before the onset of DO 8, with CH4 concentration rising only slightly. This can be explained by an early climate response of boreal wetlands, which carry the strongly depleted isotopic signature of high-latitude precipitation at that time. more...

Atmospheric methane is an important greenhouse gas and a sensitive indicator of climate change and millennial-scale temperature variability. Its concentrations over the past 650,000 years have varied between approximately 350 and approximately 800 parts per 10(9) by volume (p.p.b.v.) during glacial and interglacial periods, respectively. In comparison, present-day methane levels of approximately 1,770 p.p.b.v. have been reported. Insights into the external forcing factors and internal feedbacks controlling atmospheric methane are essential for predicting the methane budget in a warmer world. Here we present a detailed atmospheric methane record from the EPICA Dome C ice core that extends the history of this greenhouse gas to 800,000 yr before present. The average time resolution of the new data is approximately 380 yr and permits the identification of orbital and millennial-scale features. Spectral analyses indicate that the long-term variability in atmospheric methane levels is dominated by approximately 100,000 yr glacial-interglacial cycles up to approximately 400,000 yr ago with an increasing contribution of the precessional component during the four more recent climatic cycles. We suggest that changes in the strength of tropical methane sources and sinks (wetlands, atmospheric oxidation), possibly influenced by changes in monsoon systems and the position of the intertropical convergence zone, controlled the atmospheric methane budget, with an additional source input during major terminations as the retreat of the northern ice sheet allowed higher methane emissions from extending periglacial wetlands. Millennial-scale changes in methane levels identified in our record as being associated with Antarctic isotope maxima events are indicative of ubiquitous millennial-scale temperature variability during the past eight glacial cycles. more...

Past atmospheric methane concentrations show strong fluctuations in parallel to rapid glacial climate changes in the Northern Hemisphere superimposed on a glacial-interglacial doubling of methane concentrations. The processes driving the observed fluctuations remain uncertain but can be constrained using methane isotopic information from ice cores. Here we present an ice core record of carbon isotopic ratios in methane over the entire last glacial-interglacial transition. Our data show that the carbon in atmospheric methane was isotopically much heavier in cold climate periods. With the help of a box model constrained by the present data and previously published results, we are able to estimate the magnitude of past individual methane emission sources and the atmospheric lifetime of methane. We find that methane emissions due to biomass burning were about 45 Tg methane per year, and that these remained roughly constant throughout the glacial termination. The atmospheric lifetime of methane is reduced during cold climate periods. We also show that boreal wetlands are an important source of methane during warm events, but their methane emissions are essentially shut down during cold climate conditions. more...

The rate of change of climate codetermines the global warming impacts on natural and socioeconomic systems and their capabilities to adapt. Establishing past rates of climate change from temperature proxy data remains difficult given their limited spatiotemporal resolution. In contrast, past greenhouse gas radiative forcing, causing climate to change, is well known from ice cores. We compare rates of change of anthropogenic forcing with rates of natural greenhouse gas forcing since the Last Glacial Maximum and of solar and volcanic forcing of the last millennium. The smoothing of atmospheric variations by the enclosure process of air into ice is computed with a firn diffusion and enclosure model. The 20th century increase in CO(2) and its radiative forcing occurred more than an order of magnitude faster than any sustained change during the past 22,000 years. The average rate of increase in the radiative forcing not just from CO(2) but from the combination of CO(2), CH(4), and N(2)O is larger during the Industrial Era than during any comparable period of at least the past 16,000 years. In addition, the decadal-to-century scale rate of change in anthropogenic forcing is unusually high in the context of the natural forcing variations (solar and volcanoes) of the past millennium. Our analysis implies that global climate change, which is anthropogenic in origin, is progressing at a speed that is unprecedented at least during the last 22,000 years. more...

A high-resolution deuterium profile is now available along the entire European Project for Ice Coring in Antarctica Dome C ice core, extending this climate record back to marine isotope stage 20.2, approximately 800,000 years ago. Experiments performed with an atmospheric general circulation model including water isotopes support its temperature interpretation. We assessed the general correspondence between Dansgaard-Oeschger events and their smoothed Antarctic counterparts for this Dome C record, which reveals the presence of such features with similar amplitudes during previous glacial periods. We suggest that the interplay between obliquity and precession accounts for the variable intensity of interglacial periods in ice core records. more...

Polar ice cores contain, in trapped air bubbles, an archive of the concentrations of stable atmospheric gases. Of the major non-CO2 greenhouse gases, methane is measured quite routinely, while nitrous oxide is more challenging, with some artefacts occurring in the ice and so far limited interpretation. In the recent past, the ice cores provide the only direct measure of the changes that have occurred during the industrial period; they show that the current concentration of methane in the atmosphere is far outside the range experienced in the last 650,000 years; nitrous oxide is also elevated above its natural levels. There is controversy about whether changes in the pre-industrial Holocene are natural or anthropogenic in origin. Changes in wetland emissions are generally cited as the main cause of the large glacial-interglacial change in methane. However, changing sinks must also be considered, and the impact of possible newly described sources evaluated. Recent isotopic data appear to finally rule out any major impact of clathrate releases on methane at these time-scales. Any explanation must take into account that, at the rapid Dansgaard-Oeschger warmings of the last glacial period, methane rose by around half its glacial-interglacial range in only a few decades. The recent EPICA Dome C (Antarctica) record shows that methane tracked climate over the last 650,000 years, with lower methane concentrations in glacials than interglacials, and lower concentrations in cooler interglacials than in warmer ones. Nitrous oxide also shows Dansgaard-Oeschger and glacial-interglacial periodicity, but the pattern is less clear. more...

The European Project for Ice Coring in Antarctica Dome C ice core enables us to extend existing records of atmospheric methane (CH4) and nitrous oxide (N2O) back to 650,000 years before the present. A combined record of CH4 measured along the Dome C and the Vostok ice cores demonstrates, within the resolution of our measurements, that preindustrial concentrations over Antarctica have not exceeded 773 +/- 15 ppbv (parts per billion by volume) during the past 650,000 years. Before 420,000 years ago, when interglacials were cooler, maximum CH4 concentrations were only about 600 ppbv, similar to lower Holocene values. In contrast, the N2O record shows maximum concentrations of 278 +/- 7 ppbv, slightly higher than early Holocene values. more...