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

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

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

Peer reviewedPublisher PD

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

Improving Estimates of Gross Primary Productivity by Assimilating Solar-Induced Fluorescence Satellite Retrievals in a Terrestrial Biosphere Model Using a Process-Based SIF Model

Abstract Over the last few years, solar-induced chlorophyll fluorescence (SIF) observations from space have emerged as a promising resource for evaluating the spatio-temporal distribution of gross primary productivity (GPP) simulated by global terrestrial biosphere models. SIF can be used to improve GPP simulations by optimizing critical model parameters through statistical Bayesian data assimilation techniques. A prerequisite is the availability of a functional link between GPP and SIF in terrestrial biosphere models. Here we present the development of a mechanistic SIF observation operator in the ORCHIDEE (Organizing Carbon and Hydrology In Dynamic Ecosystems) terrestrial biosphere model. It simulates the regulation of photosystem II fluorescence quantum yield at the leaf level thanks to a novel parameterization of non-photochemical quenching as a function of temperature, photosynthetically active radiation, and normalized quantum yield of photochemistry. It emulates the radiative transfer of chlorophyll fluorescence to the top of the canopy using a parametric simplification of the SCOPE (Soil Canopy Observation Photosynthesis Energy) model. We assimilate two years of monthly OCO-2 (Orbiting Carbon Observatory-2) SIF product at 0.5° (2015?2016) to optimize ORCHIDEE photosynthesis and phenological parameters over an ensemble of grid points for all plant functional types. The impact on the simulated GPP is considerable with a large decrease of the global scale budget by 28 GtC/year over the period 1990?2009. The optimized GPP budget (134/136 GtC/year over 1990?2009/2001?2009) remarkably agrees with independent GPP estimates, FLUXSAT (137 GtC/year over 2001?2009) in particular and FLUXCOM (121 GtC/year over 1990?2009). Our results also suggest a biome dependency of the SIF-GPP relationship that needs to be improved for some plant functional types.Peer reviewe

Spatiotemporal patterns of terrestrial gross primary production: A review

This is the final version of the article. Available from American Geophysical Union via the DOI in this record.There is another record for this publication in ORE at http://hdl.handle.net/10871/30934Great advances have been made in the last decade in quantifying and understanding the spatiotemporal patterns of terrestrial gross primary production (GPP) with ground, atmospheric, and space observations. However, although global GPP estimates exist, each data set relies upon assumptions and none of the available data are based only on measurements. Consequently, there is no consensus on the global total GPP and large uncertainties exist in its benchmarking. The objective of this review is to assess how the different available data sets predict the spatiotemporal patterns of GPP, identify the differences among data sets, and highlight the main advantages/disadvantages of each data set. We compare GPP estimates for the historical period (1990-2009) from two observation-based data sets (Model Tree Ensemble and Moderate Resolution Imaging Spectroradiometer) to coupled carbon-climate models and terrestrial carbon cycle models from the Fifth Climate Model Intercomparison Project and TRENDY projects and to a new hybrid data set (CARBONES). Results show a large range in the mean global GPP estimates. The different data sets broadly agree on GPP seasonal cycle in terms of phasing, while there is still discrepancy on the amplitude. For interannual variability (IAV) and trends, there is a clear separation between the observation-based data that show little IAV and trend, while the process-based models have large GPP variability and significant trends. These results suggest that there is an urgent need to improve observation-based data sets and develop carbon cycle modeling with processes that are currently treated either very simplistically to correctly estimate present GPP and better quantify the future uptake of carbon dioxide by the world's vegetation.European Commission's Seventh Framework Programme. Grant Numbers: 238366, 28267

An assessment of the carbon balance of Arctic tundra:Comparisons among observations, process models, and atmospheric inversions

Although Arctic tundra has been estimated to cover only 8% of the global land surface, the large and potentially labile carbon pools currently stored in tundra soils have the potential for large emissions of carbon (C) under a warming climate. These emissions as radiatively active greenhouse gases in the form of both CO<sub>2</sub> and CH<sub>4</sub> could amplify global warming. Given the potential sensitivity of these ecosystems to climate change and the expectation that the Arctic will experience appreciable warming over the next century, it is important to assess whether responses of C exchange in tundra regions are likely to enhance or mitigate warming. In this study we compared analyses of C exchange of Arctic tundra between 1990 and 2006 among observations, regional and global applications of process-based terrestrial biosphere models, and atmospheric inversion models. Syntheses of flux observations and inversion models indicate that the annual exchange of CO<sub>2</sub> between Arctic tundra and the atmosphere has large uncertainties that cannot be distinguished from neutral balance. The mean estimate from an ensemble of process-based model simulations suggests that Arctic tundra has acted as a sink for atmospheric CO<sub>2</sub> in recent decades, but based on the uncertainty estimates it cannot be determined with confidence whether these ecosystems represent a weak or a strong sink. Tundra was 0.6 °C warmer in the 2000s compared to the 1990s. The central estimates of the observations, process-based models, and inversion models each identify stronger sinks in the 2000s compared with the 1990s. Some of the process models indicate that this occurred because net primary production increased more in response to warming than heterotrophic respiration. Similarly, the observations and the applications of regional process-based models suggest that CH<sub>4</sub> emissions from Arctic tundra have increased from the 1990s to 2000s because of the sensitivity of CH<sub>4</sub> emissions to warmer temperatures. Based on our analyses of the estimates from observations, process-based models, and inversion models, we estimate that Arctic tundra was a sink for atmospheric CO<sub>2</sub> of 110 Tg C yr<sup>−1</sup> (uncertainty between a sink of 291 Tg C yr<sup>−1</sup> and a source of 80 Tg C yr<sup>−1</sup>) and a source of CH<sub>4</sub> to the atmosphere of 19 Tg C yr<sup>−1</sup> (uncertainty between sources of 8 and 29 Tg C yr<sup>−1</sup>). The suite of analyses conducted in this study indicate that it is important to reduce uncertainties in the observations, process-based models, and inversions in order to better understand the degree to which Arctic tundra is influencing atmospheric CO<sub>2</sub> and CH<sub>4</sub> concentrations. The reduction of uncertainties can be accomplished through (1) the strategic placement of more CO<sub>2</sub> and CH<sub>4</sub> monitoring stations to reduce uncertainties in inversions, (2) improved observation networks of ground-based measurements of CO<sub>2</sub> and CH<sub>4</sub> exchange to understand exchange in response to disturbance and across gradients of climatic and hydrological variability, and (3) the effective transfer of information from enhanced observation networks into process-based models to improve the simulation of CO<sub>2</sub> and CH<sub>4</sub> exchange from Arctic tundra to the atmosphere

Spatiotemporal patterns of terrestrial gross primary production: A review

This is the final version of the article. Available from American Geophysical Union via the DOI in this record.There is another record for this publication in ORE at http://hdl.handle.net/10871/21007Great advances have been made in the last decade in quantifying and understanding the spatiotemporal patterns of terrestrial gross primary production (GPP) with ground, atmospheric, and space observations. However, although global GPP estimates exist, each data set relies upon assumptions and none of the available data are based only on measurements. Consequently, there is no consensus on the global total GPP and large uncertainties exist in its benchmarking. The objective of this review is to assess how the different available data sets predict the spatiotemporal patterns of GPP, identify the differences among data sets, and highlight the main advantages/disadvantages of each data set. We compare GPP estimates for the historical period (1990-2009) from two observation-based data sets (Model Tree Ensemble and Moderate Resolution Imaging Spectroradiometer) to coupled carbon-climate models and terrestrial carbon cycle models from the Fifth Climate Model Intercomparison Project and TRENDY projects and to a new hybrid data set (CARBONES). Results show a large range in the mean global GPP estimates. The different data sets broadly agree on GPP seasonal cycle in terms of phasing, while there is still discrepancy on the amplitude. For interannual variability (IAV) and trends, there is a clear separation between the observation-based data that show little IAV and trend, while the process-based models have large GPP variability and significant trends. These results suggest that there is an urgent need to improve observation-based data sets and develop carbon cycle modeling with processes that are currently treated either very simplistically to correctly estimate present GPP and better quantify the future uptake of carbon dioxide by the world's vegetation.European Commission's Seventh Framework Programme. Grant Numbers: 238366, 28267

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)

Land surface models (LSMs), which form the land component of earth system models, rely on numerous processes for describing carbon, water and energy budgets, often associated with highly uncertain parameters. Data assimilation (DA) is a useful approach for optimising the most critical parameters in order to improve model accuracy and refine future climate predictions. In this study, we compare two different DA methods for optimising the parameters of seven plant functional types (PFTs) of the ORCHIDEE LSM using daily averaged eddy-covariance observations of net ecosystem exchange and latent heat flux at 78 sites across the globe. We perform a technical investigation of two classes of minimisation methods – local gradient-based (the L-BFGS-B algorithm, limited memory Broyden–Fletcher–Goldfarb–Shanno algorithm with bound constraints) and global random search (the genetic algorithm) – by evaluating their relative performance in terms of the model–data fit and the difference in retrieved parameter values. We examine the performance of each method for two cases: when optimising parameters at each site independently (“single-site” approach) and when simultaneously optimising the model at all sites for a given PFT using a common set of parameters (“multi-site” approach). We find that for the single site case the random search algorithm results in lower values of the cost function (i.e. lower model–data root mean square differences) than the gradient-based method; the difference between the two methods is smaller for the multi-site optimisation due to a smoothing of the cost function shape with a greater number of observations. The spread of the cost function, when performing the same tests with 16 random first-guess parameters, is much larger with the gradient-based method, due to the higher likelihood of being trapped in local minima. When using pseudo-observation tests, the genetic algorithm results in a closer approximation of the true posterior parameter value in the L-BFGS-B algorithm. We demonstrate the advantages and challenges of different DA techniques and provide some advice on using it for the LSM parameter optimisation.</p

Multicriteria evaluation of discharge simulation in Dynamic Global Vegetation Models

PublishedJournal Article© 2015. American Geophysical Union. All Rights Reserved. In this study, we assessed the performance of discharge simulations by coupling the runoff from seven Dynamic Global Vegetation Models (DGVMs; LPJ, ORCHIDEE, Sheffield-DGVM, TRIFFID, LPJ-GUESS, CLM4CN, and OCN) to one river routing model for 16 large river basins. The results show that the seasonal cycle of river discharge is generally modeled well in the low and middle latitudes but not in the high latitudes, where the peak discharge (due to snow and ice melting) is underestimated. For the annual mean discharge, the DGVMs chained with the routing model show an underestimation. Furthermore, the 30 year trend of discharge is also underestimated. For the interannual variability of discharge, a skill score based on overlapping of probability density functions (PDFs) suggests that most models correctly reproduce the observed variability (correlation coefficient higher than 0.5; i.e., models account for 50% of observed interannual variability) except for the Lena, Yenisei, Yukon, and the Congo river basins. In addition, we compared the simulated runoff from different simulations where models were forced with either fixed or varying land use. This suggests that both seasonal and annual mean runoff has been little affected by land use change but that the trend itself of runoff is sensitive to land use change. None of the models when considered individually show significantly better performances than any other and in all basins. This suggests that based on current modeling capability, a regional-weighted average of multimodel ensemble projections might be appropriate to reduce the bias in future projection of global river discharge.National Natural Science Foundation of China. Grant Numbers: 41125004, 31321061, Chinese Ministry of Environmental Protection. Grant Number: 201209031, 111 Project. Grant Number: B14001, National Youth Top-notch Talent Support Program in China, Imbalance-P ERC-synergy, TRENDY, Global River Discharge Cente