Importance of fossil fuel emission uncertainties over Europe for CO2 modeling: model intercomparison

Inverse modeling techniques used to quantify surface carbon fluxes commonly assume that the uncertainty of fossil fuel CO2 (FFCO2) emissions is negligible and that intra-annual variations can be neglected. To investigate these assumptions, we analyzed the differences between four fossil fuel emission inventories with spatial and temporal differences over Europe and their impact on the model simulated CO2 concentration. Large temporal flux variations characterize the hourly fields (~40 % and ~80 % for the seasonal and diurnal cycles, peak-to-peak) and annual country totals differ by 10 % on average and up to 40 % for some countries (i.e., the Netherlands). These emissions have been prescribed to seven different transport models, resulting in 28 different FFCO2 concentrations fields. The modeled FFCO2 concentration time series at surface sites using time-varying emissions show larger seasonal cycles (+2 ppm at the Hungarian tall tower (HUN)) and smaller diurnal cycles in summer (-1 ppm at HUN) than when using constant emissions. The concentration range spanned by all simulations varies between stations, and is generally larger in winter (up to ~10 ppm peak-to-peak at HUN) than in summer (~5 ppm). The contribution of transport model differences to the simulated concentration std-dev is 2–3 times larger than the contribution of emission differences only, at typical European sites used in global inversions. These contributions to the hourly (monthly) std-dev's amount to ~1.2 (0.8) ppm and ~0.4 (0.3) ppm for transport and emissions, respectively. First comparisons of the modeled concentrations with 14C-based fossil fuel CO2 observations show that the large transport differences still hamper a quantitative evaluation/validation of the emission inventories. Changes in the estimated monthly biosphere flux (Fbio) over Europe, using two inverse modeling approaches, are relatively small (less that 5 %) while changes in annual Fbio (up to ~0.15 % GtC yr-1) are only slightly smaller than the differences in annual emission totals and around 30 % of the mean European ecosystem carbon sink. These results point to an urgent need to improve not only the transport models but also the assumed spatial and temporal distribution of fossil fuel emission inventories

Non-steady-state, non-uniform transpiration rate and leaf anatomy effects on the progressive stable isotope enrichment of leaf water along monocot leaves

This study focuses on the spatial patterns of transpiration-driven water isotope enrichment (Δlw) along monocot leaves. It has been suggested that these spatial patterns are the result of competing effects of advection and (back-)diffusion of water isot

Land surface model parameter optimisation using in situ flux data : Comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2)

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 and USCCC. We acknowledge the financial support to the eddy covariance data harmonisation provided by CarboEuropeIP, FAO-GTOS-TCO, iLEAPS, Max Planck Institute for Biogeochemistry, National Science Foundation, University of Tuscia, Universiteì Laval, 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 and the University of Virginia.Peer reviewedPublisher PD

Energy and CO2 fluxes in temperate deciduous broadleaf forests: using local add-covariance data for a global optimization of ORCHIDEE model

Non peer reviewedPublisher PD

Towards a representation of priming on soil carbon decomposition in the global land biosphere model ORCHIDEE (version 1.9.5.2)

Priming of soil carbon decomposition encompasses different processes through which the decomposition of native (already present) soil organic matter is amplified through the addition of new organic matter, with new inputs typically being more labile than the native soil organic matter. Evidence for priming comes from laboratory and field experiments, but to date there is no estimate of its impact at global scale and under the current anthropogenic perturbation of the carbon cycle. Current soil carbon decomposition models do not include priming mechanisms, thereby introducing uncertainty when extrapolating short-term local observations to ecosystem and regional to global scale. In this study we present a simple conceptual model of decomposition priming, called PRIM, able to reproduce laboratory (incubation) and field (litter manipulation) priming experiments. Parameters for this model were first optimized against data from 20 soil incubation experiments using a Bayesian framework. The optimized parameter values were evaluated against another set of soil incubation data independent from the ones used for calibration and the PRIM model reproduced the soil incubations data better than the original, CENTURY-type soil decomposition model, whose decomposition equations are based only on first-order kinetics. We then compared the PRIM model and the standard first-order decay model incorporated into the global land biosphere model ORCHIDEE (Organising Carbon and Hydrology In Dynamic Ecosystems). A test of both models was performed at ecosystem scale using litter manipulation experiments from five sites. Although both versions were equally able to reproduce observed decay rates of litter, only ORCHIDEE-PRIM could simulate the observed priming (R² = 0.54)in cases where litter was added or removed. This result suggests that a conceptually simple and numerically tractable representation of priming adapted to global models is able to capture the sign and magnitude of the priming of litter and soil organic matter

Constraining a global ecosystem model with multi-site eddy-covariance data

Peer reviewedPublisher PD

Modeling of biospheric CO2 gross fluxes via oxygen isotopes in a spruce forest canopy: a 222Rn calibrated box model approach

One-dimensional box model estimates of biospheric CO2 gross fluxes are presented. The results are based on measurements performed during the EUROSIBERIAN CARBONFLUX intensive campaign between July 27 to August 1, 1999 in a natural Picea abies forest in Russia. CO2 mixing ratios and stable isotope ratios of CO2 were measured on flask samples taken in two heights within the canopy. Simultaneously, soil and leaf samples were collected and analysed to derive the 18O/16O ratio of the respective water reservoirs and the 13C/12C ratio of the leaf tissue. The main objective of this project was to investigate biospheric gas exchange with soil and vegetation and, thereby, take advantage of the potential of the 18O/16O ratio in atmospheric CO2. Via exchange of oxygen isotopes with associated liquid water reservoirs, leaf CO2 assimilation fluxes generally enrich while soil CO2 respiration fluxes generally deplete the 18O/16O ratio of atmospheric CO2. In the model, we parameterised intra canopy transport by exploiting soil-borne 222Rn as a tracer for turbulent transport. Our approach showed that, in principle, a net ecosystem CO2 flux can be separated into assimilation and respiration fluxes using oxygen isotopes. However, quantitative partitioning is highly sensitive to the respective discrimination factors, and, therefore, also on the parameterisation of internal leaf CO2 concentrations and gradients

Data-based estimates of the ocean carbon sink variability – First results of the Surface Ocean pCO2 Mapping intercomparison (SOCOM)

Using measurements of the surface-ocean CO2 partial pressure (pCO2) and 14 different pCO2 mapping methods recently collated by the Surface Ocean pCO2 Mapping intercomparison (SOCOM) initiative, variations in regional and global sea–air CO2 fluxes are investigated. Though the available mapping methods use widely different approaches, we find relatively consistent estimates of regional pCO2 seasonality, in line with previous estimates. In terms of interannual variability (IAV), all mapping methods estimate the largest variations to occur in the eastern equatorial Pacific. Despite considerable spread in the detailed variations, mapping methods that fit the data more closely also tend to agree more closely with each other in regional averages. Encouragingly, this includes mapping methods belonging to complementary types – taking variability either directly from the pCO2 data or indirectly from driver data via regression. From a weighted ensemble average, we find an IAV amplitude of the global sea–air CO2 flux of 0.31 PgC yr−1 (standard deviation over 1992–2009), which is larger than simulated by biogeochemical process models. From a decadal perspective, the global ocean CO2 uptake is estimated to have gradually increased since about 2000, with little decadal change prior to that. The weighted mean net global ocean CO2 sink estimated by the SOCOM ensemble is −1.75 PgC yr−1 (1992–2009), consistent within uncertainties with estimates from ocean-interior carbon data or atmospheric oxygen trend

Assessing a New Clue to How Much Carbon Plants Take Up

Current climate models disagree on how much carbon dioxide land ecosystems take up for photosynthesis. Tracking the stronger carbonyl sulfide signal could help