The closing door of climate targets

Robust evidence from a range of climate–carbon cycle models shows that the maximum warming relative to pre-industrial times caused by the emissions of carbon dioxide is nearly proportional to the total amount of emitted anthropogenic carbon (1, 2). This proportionality is a reasonable approximation for simulations covering many emissions scenarios for the time frame 1750 to 2500 (1). This linear relationship is remarkable given the different complexities of the models and the wide range of emissions scenarios considered. It has direct implications for the possibility of achieving internationally agreed climate targets such as those mentioned in the Copenhagen Accord and the Cancun Agreements (3, 4). Here I explain some of the implications of the linear relationship between peak warming and total cumulative carbon emissions.

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.

On the evolution of the oceanic component of the IPSL climate models from CMIP3 to CMIP5: A mean state comparison

This study analyses the impact on the oceanic mean state of the evolution of the oceanic component (NEMO) of the climate model developed at Institut Pierre Simon Laplace (IPSL-CM), from the version IPSL-CM4, used for third phase of the Coupled Model Intercomparison Project (CMIP3), to IPSL-CM5A, used for CMIP5. Several modifications have been implemented between these two versions, in particular an interactive coupling with a biogeochemical module, a 3-band model for the penetration of the solar radiation, partial steps at the bottom of the ocean and a set of physical parameterisations to improve the representation of the impact of turbulent and tidal mixing. A set of forced and coupled experiments is used to single out the effect of each of these modifications and more generally the evolution of the oceanic component on the IPSL coupled models family. Major improvements are located in the Southern Ocean, where physical parameterisations such as partial steps and tidal mixing reinforce the barotropic transport of water mass, in particular in the Antarctic Circumpolar Current) and ensure a better representation of Antarctic bottom water masses. However, our analysis highlights that modifications, which substantially improve ocean dynamics in forced configuration, can yield or amplify biases in coupled configuration. In particular, the activation of radiative biophysical coupling between biogeochemical cycle and ocean dynamics results in a cooling of the ocean mean state. This illustrates the difficulty to improve and tune coupled climate models, given the large number of degrees of freedom and the potential compensating effects masking some biases.

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.

Fruit production in three masting tree species does not rely on stored carbon reserves

Fruiting is typically considered to massively burden the seasonal carbon budget of trees. The cost of reproduction has therefore been suggested as a proximate factor explaining observed mast-fruiting patterns. Here, we used a large-scale, continuous 13C labeling of mature, deciduous trees in a temperate Swiss forest to investigate to what extent fruit formation in three species with masting reproduction behavior (Carpinus betulus, Fagus sylvatica, Quercus petraea) relies on the import of stored carbon reserves. Using a free-air CO2 enrichment system, we exposed trees to 13C-depleted CO2 during 8 consecutive years. By the end of this experiment, carbon reserve pools had significantly lower δ13C values compared to control trees. δ13C analysis of new biomass during the first season after termination of the CO2 enrichment allowed us to distinguish the sources of built-in carbon (old carbon reserves vs. current assimilates). Flowers and expanding leaves carried a significant 13C label from old carbon stores. In contrast, fruits and vegetative infructescence tissues were exclusively produced from current, unlabeled photoassimilates in all three species, including F. sylvatica, which had a strong masting season. Analyses of δ13C in purified starch from xylem of fruit-bearing shoots revealed a complete turn-over of starch during the season, likely due to its usage for bud break. This study is the first to directly demonstrate that fruiting is independent from old carbon reserves in masting trees, with significant implications for mechanistic models that explain mast seeding.

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.

Anthropogenic perturbation of the carbon fluxes from land to ocean

A substantial amount of the atmospheric carbon taken up on land through photosynthesis and chemical weathering is transported laterally along the aquatic continuum from upland terrestrial ecosystems to the ocean. So far, global carbon budget estimates have implicitly assumed that the transformation and lateral transport of carbon along this aquatic continuum has remained unchanged since pre-industrial times. A synthesis of published work reveals the magnitude of present-day lateral carbon fluxes from land to ocean, and the extent to which human activities have altered these fluxes. We show that anthropogenic perturbation may have increased the flux of carbon to inland waters by as much as 1.0 Pg C yr(-1) since pre-industrial times, mainly owing to enhanced carbon export from soils. Most of this additional carbon input to upstream rivers is either emitted back to the atmosphere as carbon dioxide (similar to 0.4 Pg C yr(-1)) or sequestered in sediments (similar to 0.5 Pg C yr(-1)) along the continuum of freshwater bodies, estuaries and coastal waters, leaving only a perturbation carbon input of similar to 0.1 Pg C yr(-1) to the open ocean. According to our analysis, terrestrial ecosystems store similar to 0.9 Pg C yr(-1) at present, which is in agreement with results from forest inventories but significantly differs from the figure of 1.5 Pg C yr(-1) previously estimated when ignoring changes in lateral carbon fluxes. We suggest that carbon fluxes along the land-ocean aquatic continuum need to be included in global carbon dioxide budgets.

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

Simulating Carbon Stocks and Fluxes of an African Tropical Montane Forest with an Individual-Based Forest Model

Tropical forests are carbon-dense and highly productive ecosystems. Consequently, they play an important role in the global carbon cycle. In the present study we used an individual-based forest model (FORMIND) to analyze the carbon balances of a tropical forest. The main processes of this model are tree growth, mortality, regeneration, and competition. Model parameters were calibrated using forest inventory data from a tropical forest at Mt. Kilimanjaro. The simulation results showed that the model successfully reproduces important characteristics of tropical forests (aboveground biomass, stem size distribution and leaf area index). The estimated aboveground biomass (385 t/ha) is comparable to biomass values in the Amazon and other tropical forests in Africa. The simulated forest reveals a gross primary production of 24 tcha-1yr-1. Modeling above- and belowground carbon stocks, we analyzed the carbon balance of the investigated tropical forest. The simulated carbon balance of this old-growth forest is zero on average. This study provides an example of how forest models can be used in combination with forest inventory data to investigate forest structure and local carbon balances.

Water-use efficiency and transpiration across European forests during the Anthropocene

The Earth’s carbon and hydrologic cycles are intimately coupled by gas exchange through plant stomata1, 2, 3. However, uncertainties in the magnitude4, 5, 6 and consequences7, 8 of the physiological responses9, 10 of plants to elevated CO2 in natural environments hinders modelling of terrestrial water cycling and carbon storage11. Here we use annually resolved long-term δ13C tree-ring measurements across a European forest network to reconstruct the physiologically driven response of intercellular CO2 (Ci) caused by atmospheric CO2 (Ca) trends. When removing meteorological signals from the δ13C measurements, we find that trees across Europe regulated gas exchange so that for one ppmv atmospheric CO2 increase, Ci increased by ~0.76 ppmv, most consistent with moderate control towards a constant Ci/Ca ratio. This response corresponds to twentieth-century intrinsic water-use efficiency (iWUE) increases of 14 ± 10 and 22 ± 6% at broadleaf and coniferous sites, respectively. An ensemble of process-based global vegetation models shows similar CO2 effects on iWUE trends. Yet, when operating these models with climate drivers reintroduced, despite decreased stomatal opening, 5% increases in European forest transpiration are calculated over the twentieth century. This counterintuitive result arises from lengthened growing seasons, enhanced evaporative demand in a warming climate, and increased leaf area, which together oppose effects of CO2-induced stomatal closure. Our study questions changes to the hydrological cycle, such as reductions in transpiration and air humidity, hypothesized to result from plant responses to anthropogenic emissions.

Global carbon budget 2013

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil-fuel combustion and cement production (EFF) are based on energy statistics, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2 and land cover change (some including nitrogen–carbon interactions). All uncertainties are reported as ± 1 σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2003–2012), EFF was 8.6 ± 0.4 GtC yr − 1, ELUC 0.9 ± 0.5 GtC yr − 1, GATM 4.3 ± 0.1 GtC yr − 1, S OCEAN 2.5 ± 0.5 GtC yr − 1, and S LAND 2.8 ± 0.8 GtC yr − 1. For year 2012 alone, EFF grew to 9.7 ± 0.5 GtC yr − 1, 2.2 % above 2011, reflecting a continued growing trend in these emissions, GATM was 5.1 ± 0.2 GtC yr − 1, SOCEANwas 2.9 ± 0.5 GtC yr −1, and assuming an ELU Cof 1.0 ± 0.5 GtC yr − 1 (based on the 2001–2010 average), SLAND was 2.7 ± 0.9 GtC yr − 1. GATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high EFF. The global atmospheric CO2 con- centration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that EFF will increase by 2.1 % (1.1–3.1 %) to 9.9 ± 0.5 GtC in 2013, 61 % above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO2 will reach about 535 ± 55 GtC for 1870–2013, about 70 % from EFF (390 ± 20 GtC) and 30 % from ELUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center.

Interaction of corroding iron with bentonite in the ABM1 experiment at Äspö, Sweden: A microscopic approach

Bentonite and iron metals are common materials proposed for use in deep-seated geological repositories for radioactive waste. The inevitable corrosion of iron leads to interaction processes with the clay which may affect the sealing properties of the bentonite backfill. The objective of the present study was to improve our understanding of this process by studying the interface between iron and compacted bentonite in a geological repository-type setting. Samples of MX-80 bentonite samples which had been exposed to an iron source and elevated temperatures (up to 115ºC) for 2.5 y in an in situ experiment (termed ABM1) at the Äspö Hard Rock Laboratory, Sweden, were investigated by microscopic means, including scanning electron microscopy, μ-Raman spectroscopy, spatially resolved X-ray diffraction, and X-ray fluorescence. The corrosion process led to the formation of a ~100 mm thick corrosion layer containing siderite, magnetite, some goethite, and lepidocrocite mixed with the montmorillonitic clay. Most of the corroded Fe occurred within a 10 mm-thick clay layer adjacent to the corrosion layer. An average corrosion depth of the steel of 22–35 μm and an average Fe2+ diffusivity of 1–26×10–13 m2/s were estimated based on the properties of the Fe-enriched clay layer. In that layer, the corrosion-derived Fe occurred predominantly in the clay matrix. The nature of this Fe could not be identified. No indications of clay transformation or newly formed clay phases were found. A slight enrichment of Mg close to the Fe–clay contact was observed. The formation of anhydrite and gypsum, and the dissolution of some SiO2 resulting from the temperature gradient in the in situ test, were also identified. © 2014, Clay Minerals Society. All right reserved.

Holocene biomass burning and global dynamics of the carbon cycle

Fire regimes have changed during the Holocene due to changes in climate, vegetation, and in human practices. Here, we hypothesise that changes in fire regime may have affected the global CO2 concentration in the atmosphere through the Holocene. Our data are based on quantitative reconstructions of biomass burning deduced from stratified charcoal records from Europe, and South-, Central- and North America, and Oceania to test the fire-carbon release hypothesis. In Europe the significant increase of fire activity is dated ≈6000 cal. yr ago. In north-eastern North America burning activity was greatest before 7500 years ago, very low between 7500–3000 years, and has been increasing since 3000 years ago. In tropical America, the pattern is more complex and apparently latitudinally zonal. Maximum burning occurred in the southern Amazon basin and in Central America during the middle Holocene, and during the last 2000 years in the northern Amazon basin. In Oceania, biomass burning has decreased since a maximum 5000 years ago. Biomass burning has broadly increased in the Northern and Southern hemispheres throughout the second half of the Holocene associated with changes in climate and human practices. Global fire indices parallel the increase of atmospheric CO2 concentration recorded in Antarctic ice cores. Future issues on carbon dynamics relatively to biomass burning are discussed to improve the quantitative reconstructions.

Radiocarbon dating of modern peat profiles: Pre- and post-bomb C-14 variations in the construction of age-depth models

We present studies of 9 modern (up to 400-yr-old) peat sections from Slovenia, Switzerland, Austria, Italy, and Finland. Precise radiocarbon dating of modern samples is possible due to the large bomb peak of atmospheric 14C concentration in 1963 and the following rapid decline in the 14C level. All the analyzed 14C profiles appeared concordant with the shape of the bomb peak of atmospheric 14C concentration, integrated over some time interval with a length specific to the peat section. In the peat layers covered by the bomb peak, calendar ages of individual peat samples could be determined almost immediately, with an accuracy of 23 yr. In the pre-bomb sections, the calendar ages of individual dated samples are determined in the form of multi-modal probability distributions of about 300 yr wide (about AD 16501950). However, simultaneous use of the post-bomb and pre-bomb 14C dates, and lithological information, enabled the rejection of most modes of probability distributions in the pre-bomb section. In effect, precise age-depth models of the post-bomb sections have been extended back in time, into the wiggly part of the 14C calibration curve.

Comparison of 3D kinetic and hydrodynamic models to ROSINA-COPS measurements of the neutral coma of 67P/Churyumov-Gerasimenko

67P/Churyumov-Gerasimenko (67P) is a Jupiter-family comet and the object of investigation of the European Space Agency mission Rosetta. This report presents the first full 3D simulation results of 67P’s neutral gas coma. In this study we include results from a direct simulation Monte Carlo method, a hydrodynamic code, and a purely geometric calculation which computes the total illuminated surface area on the nucleus. All models include the triangulated 3D shape model of 67P as well as realistic illumination and shadowing conditions. The basic concept is the assumption that these illumination conditions on the nucleus are the main driver for the gas activity of the comet. As a consequence, the total production rate of 67P varies as a function of solar insolation. The best agreement between the model and the data is achieved when gas fluxes on the night side are in the range of 7% to 10% of the maximum flux, accounting for contributions from the most volatile components. To validate the output of our numerical simulations we compare the results of all three models to in situ gas number density measurements from the ROSINA COPS instrument. We are able to reproduce the overall features of these local neutral number density measurements of ROSINA COPS for the time period between early August 2014 and January 1 2015 with all three models. Some details in the measurements are not reproduced and warrant further investigation and refinement of the models. However, the overall assumption that illumination conditions on the nucleus are at least an important driver of the gas activity is validated by the models. According to our simulation results we find the total production rate of 67P to be constant between August and November 2014 with a value of about 1 × 10²⁶ molecules s⁻¹.