Engineering Papers⌕ Search

Engineering topics

Peylin, Philippe

Publications and source records attributed to Peylin, Philippe.

Trends and Drivers of Terrestrial Sources and Sinks of Carbon Dioxide: An Overview of the TRENDY Project

The terrestrial biosphere plays a major role in the global carbon cycle, and there is a recognized need for regularly updated estimates of land-atmosphere exchange at regional and global scales. An international ensemble of Dynamic Global Vegetation Models (DGVMs), known as the “Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide” (TRENDY) project, quantifies land biophysical exchange processes and biogeochemistry cycles in support of the annual Global Carbon Budget assessments and the REgional Carbon Cycle Assessment and Processes, phase 2 project. DGVMs use a common protocol and set of driving data sets. A set of factorial simulations allows attribution of spatio-temporal changes in land surface processes to three primary global change drivers: changes in atmospheric CO 2 , climate change and variability, and Land Use and Land Cover Changes (LULCC). Here, we describe the TRENDY project, benchmark DGVM performance using remote-sensing and other observational data, and present results for the contemporary period. Simulation results show a large global carbon sink in natural vegetation over 2012–2021, attributed to the CO 2 fertilization effect (3.8 ± 0.8 PgC/yr) and climate (–0.58 ± 0.54 PgC/yr). Forests and semi-arid ecosystems contribute approximately equally to the mean and trend in the natural land sink, and semi-arid ecosystems continue to dominate interannual variability. The natural sink is offset by net emissions from LULCC (–1.6 ± 0.5 PgC/yr), with a net land sink of 1.7 ± 0.6 PgC/yr. Despite the largest gross fluxes being in the tropics, the largest net land-atmosphere exchange is simulated in the extratropical regions.

58 GEOSCIENCES↗

Carbon and Water Fluxes of the Boreal Evergreen Needleleaf Forest Biome Constrained by Assimilating Ecosystem Carbonyl Sulfide Flux Observations

Gross primary production (GPP) by boreal forests is highly sensitive to environmental changes. However, GPP simulated by land surface models (LSMs) remains highly uncertain due to the lack of direct photosynthesis observations at large scales. Carbonyl sulfide (COS) has emerged as a promising proxy to improve the representation of GPP in LSMs. Because COS is absorbed by vegetation following the same diffusion pathway as CO2 during photosynthesis and not emitted back to the atmosphere, incorporating a mechanistic representation of vegetation COS uptake in LSMs allows using COS observations to refine GPP representation. Here, we perform ecosystem COS flux and GPP data assimilations to constrain the COS- and GPP-related parameters in the ORCHIDEE LSM for boreal evergreen needleleaf forests (BorENF). Assimilating ecosystem COS fluxes at Hyytiälä forest increases the simulated net ecosystem COS uptake by 14%. This increase largely results from changes in the internal conductance to COS, highlighting the need to improve the representation of COS internal diffusion and consumption. Moreover, joint assimilation of ecosystem COS flux and GPP at Hyytiälä improves the simulated latent heat flux, contrary to the GPP-only data assimilation, which fails to do so. Finally, we scaled this assimilation framework up to the boreal region and find that the joint assimilation of COS at Hyytiälä and GPP fluxes at 10 BorENF sites increases the modeled vegetation COS uptake up to 18%, but not GPP. Therefore, this study encourages the use of COS flux observations to inform GPP and latent heat flux representations in LSMs.

54 ENVIRONMENTAL SCIENCES↗

Assimilation of multiple datasets results in large differences in regional- to global-scale NEE and GPP budgets simulated by a terrestrial biosphere model

In spite of the importance of land ecosystems in offsetting carbon dioxide emissions released by anthropogenic activities into the atmosphere, the spatiotemporal dynamics of terrestrial carbon fluxes remain largely uncertain at regional to global scales. Over the past decade, data assimilation (DA) techniques have grown in importance for improving these fluxes simulated by terrestrial biosphere models (TBMs), by optimizing model parameter values while also pinpointing possible parameterization deficiencies. Although the joint assimilation of multiple data streams is expected to constrain a wider range of model processes, their actual benefits in terms of reduction in model uncertainty are still under-researched, also given the technical challenges. In this study, we investigated with a consistent DA framework and the ORCHIDEE-LMDz TBM–atmosphere model how the assimilation of different combinations of data streams may result in different regional to global carbon budgets. To do so, we performed comprehensive DA experiments where three datasets (in situ measurements of net carbon exchange and latent heat fluxes, spaceborne estimates of the normalized difference vegetation index, and atmospheric CO 2 concentration data measured at stations) were assimilated alone or simultaneously. We thus evaluated their complementarity and usefulness to constrain net and gross C land fluxes. We found that a major challenge in improving the spatial distribution of the land C sinks and sources with atmospheric CO 2 data relates to the correction of the soil carbon imbalance.

54 ENVIRONMENTAL SCIENCES↗

Are Land‐Use Change Emissions in Southeast Asia Decreasing or Increasing?

Abstract Southeast Asia is a region known for active land‐use changes (LUC) over the past 60 years; yet, how trends in net CO 2 uptake and release resulting from LUC activities (net LUC flux) have changed through past decades remains uncertain. The level of uncertainty in net LUC flux from process‐based models is so high that it cannot be concluded that newer estimates are necessarily more reliable than older ones. Here, we examined net LUC flux estimates of Southeast Asia for the 1980s−2010s from older and newer sets of Dynamic Global Vegetation Model simulations (TRENDY v2 and v7, respectively), and forcing data used for running those simulations, along with two book‐keeping estimates (H&N and BLUE). These estimates yielded two contrasting historical LUC transitions, such that TRENDY v2 and H&N showed a transition from increased emissions from the 1980s to 1990s to declining emissions in the 2000s, while TRENDY v7 and BLUE showed the opposite transition. We found that these contrasting transitions originated in the update of LUC forcing data, which reduced the loss of forest area during the 1990s. Further evaluation of remote sensing studies, atmospheric inversions, and the history of forestry and environmental policies in Southeast Asia supported the occurrence of peak emissions in the 1990s and declining thereafter. However, whether LUC emissions continue to decline in Southeast Asia remains uncertain as key processes in recent years, such as conversion of peat forest to oil‐palm plantation, are yet to be represented in the forcing data, suggesting a need for further revision.

54 ENVIRONMENTAL SCIENCES↗

Dynamic global vegetation models underestimate net CO 2 flux mean and inter-annual variability in dryland ecosystems

Despite their sparse vegetation, dryland regions exert a huge influence over global biogeochemical cycles because they cover more than 40% of the world surface (Schimel 2010 Science 327 418–9). It is thought that drylands dominate the inter-annual variability (IAV) and long-term trend in the global carbon (C) cycle (Poulter et al 2014 Nature 509 600–3, Ahlstrom et al 2015 Science 348 895–9, Zhang et al 2018 Glob. Change Biol. 24 3954–68). Projections of the global land C sink therefore rely on accurate representation of dryland C cycle processes; however, the dynamic global vegetation models (DGVMs) used in future projections have rarely been evaluated against dryland C flux data. Here, we carried out an evaluation of 14 DGVMs (TRENDY v7) against net ecosystem exchange (NEE) data from 12 dryland flux sites in the southwestern US encompassing a range of ecosystem types (forests, shrub- and grasslands). We find that all the models underestimate both mean annual C uptake/release as well as the magnitude of NEE IAV, suggesting that improvements in representing dryland regions may improve global C cycle projections. Across all models, the sensitivity and timing of ecosystem C uptake to plant available moisture was at fault. Spring biases in gross primary production (GPP) dominate the underestimate of mean annual NEE, whereas models' lack of GPP response to water availability in both spring and summer monsoon are responsible for inability to capture NEE IAV. Errors in GPP moisture sensitivity at high elevation forested sites were more prominent during the spring, while errors at the low elevation shrub and grass-dominated sites were more important during the monsoon. We propose a range of hypotheses for why model GPP does not respond sufficiently to changing water availability that can serve as a guide for future dryland DGVM developments. Our analysis suggests that improvements in modeling C cycle processes across more than a quarter of the Earth's land surface could be achieved by addressing the moisture sensitivity of dryland C uptake.

54 ENVIRONMENTAL SCIENCES↗

Global Carbon Budget 2018

Accurate assessment of anthropogenic carbon dioxide ( CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – 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 methodology to quantify the five major components ofthe global carbon budget and their uncertainties. Fossil CO 2 emissions ( E FF ) are based on energy statistics and cement production data, while emissions from land use and land-use change ( E LUC ), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO 2 concentration is measured directly and its growth rate ( G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink ( S OCEAN ) and terrestrial CO 2 sink ( S LAND ) are estimated with global process models constrained by observations. The resulting carbon budget imbalance ( B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1 σ . For the last decade available (2008–2017), E FF was 9.4±0.5 GtC yr -1 , E LUC 1.5±0.7 GtC yr -1 , G ATM 4.7±0.02 GtC yr -1 , S OCEAN 2.4±0.5 GtC yr -1 , and S LAND 3.2±0.8 GtC yr -1 , with a budget imbalance B IM of0.5 GtC yr -1 indicating overestimated emissions and/or underestimated sinks. For the year 2017 alone, the growth in E FF was about 1.6 %and emissions increased to 9.9±0.5 GtC yr -1 . Also for 2017, E LUC was 1.4±0.7 GtC yr -1 , G ATM was 4.6±0.2 GtC yr -1 , S OCEAN was 2.5±0.5 GtC yr -1 , and S LAND was 3.8±0.8 GtC yr -1 ,with a B IM of 0.3 GtC. The global atmospheric CO 2 concentration reached 405.0±0.1 ppm averaged over 2017.For 2018, preliminary data for the first 6–9 months indicate a renewed growth in E FF of + 2.7 % (range of 1.8 % to 3.7 %) based on national emission projections for China, the US, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. The analysis presented here shows that the mean and trend in the five components of the global carbon budget are consistently estimated over the period of 1959–2017,but discrepancies of up to 1 GtC yr -1 persist for the representation of semi-decadal variability in CO 2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations show (1) no consensus in the mean and trend in land-use change emissions, (2) a persistent low agreement among the different methods on the magnitude of the land CO 2 flux in the northern extra-tropics,and (3) an apparent underestimation of the CO 2 variability by ocean models, originating outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018, 2016,2015a, b, 2014, 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2018 .

54 ENVIRONMENTAL SCIENCES↗

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.

58 GEOSCIENCES↗

Reviews and syntheses: Carbonyl sulfide as a multi-scale tracer for carbon and water cycles

Abstract. For the past decade, observations of carbonyl sulfide (OCS or COS)have been investigated as a proxy for carbon uptake by plants. OCS is destroyed by enzymes that interact with CO 2 during photosynthesis, namely carbonic anhydrase (CA) and RuBisCO, where CA is the more important one. The majority of sources of OCS to the atmosphere are geographically separated from this large plant sink,whereas the sources and sinks of CO 2 are co-located in ecosystems. The draw down of OCS can therefore be related to the uptake of CO 2 without the added complication of co-located emissions comparable in magnitude. Here we review the state of our understanding of the global OCS cycle and its applications to ecosystem carbon cycle science. OCS uptake is correlated well toplant carbon uptake, especially at the regional scale. OCS can be used in conjunction with other independent measures of ecosystem function, like solar-induced fluorescence and carbon and water isotope studies. More work needs to be done to generate global coverage for OCS observations and to link this powerful atmospheric tracer to systems where fundamental questions concerning the carbon and water cycle remain.

54 ENVIRONMENTAL SCIENCES↗

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

Here, large uncertainties in land surface models (LSMs) simulations still arise from inaccurate forcing, poor description of land surface heterogeneity (soil and vegetation properties), incorrect model parameter values and incomplete representation of biogeochemical processes. The recent increase in the number and type of carbon cycle-related observations, including both in situ and remote sensing measurements, has opened a new road to optimize model parameters via robust statistical model–data integration techniques, in order to reduce the uncertainties of simulated carbon fluxes and stocks. In this study we present a carbon cycle data assimilation system that assimilates three major data streams, namely the Moderate Resolution Imaging Spectroradiometer (MODIS)-Normalized Difference Vegetation Index (NDVI) observations of vegetation activity, net ecosystem exchange (NEE) and latent heat (LE) flux measurements at more than 70 sites (FLUXNET), as well as atmospheric CO 2 concentrations at 53 surface stations, in order to optimize the main parameters (around 180 parameters in total) of the Organizing Carbon and Hydrology in Dynamics Ecosystems (ORCHIDEE) LSM (version 1.9.5 used for the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations). The system relies on a stepwise approach that assimilates each data stream in turn, propagating the information gained on the parameters from one step to the next. Overall, the ORCHIDEE model is able to achieve a consistent fit to all three data streams, which suggests that current LSMs have reached the level of development to assimilate these observations. The assimilation of MODIS-NDVI (step 1) reduced the growing season length in ORCHIDEE for temperate and boreal ecosystems, thus decreasing the global mean annual gross primary production (GPP). Using FLUXNET data (step 2) led to large improvements in the seasonal cycle of the NEE and LE fluxes for all ecosystems (i.e., increased amplitude for temperate ecosystems). The assimilation of atmospheric CO 2 , using the general circulation model (GCM) of the Laboratoire de Météorologie Dynamique (LMDz; step 3), provides an overall constraint (i.e., constraint on large-scale net CO 2 fluxes), resulting in an improvement of the fit to the observed atmospheric CO 2 growth rate. Thus, the optimized model predicts a land C (carbon) sink of around 2.2 PgC yr -1 (for the 2000–2009 period), which is more compatible with current estimates from the Global Carbon Project (GCP) than the prior value. The consistency of the stepwise approach is evaluated with back-compatibility checks. The final optimized model (after step 3) does not significantly degrade the fit to MODIS-NDVI and FLUXNET data that were assimilated in the first two steps, suggesting that a stepwise approach can be used instead of the more “challenging” implementation of a simultaneous optimization in which all data streams are assimilated together. Most parameters, including the scalar of the initial soil carbon pool size, changed during the optimization with a large error reduction. This work opens new perspectives for better predictions of the land carbon budgets.

58 GEOSCIENCES↗