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

Acknowledgements. This work was mainly funded by the EU FP7 CARBONES project (contracts FP7-SPACE-2009-1-242316), with also a small contribution from GEOCARBON project (ENV.2011.4.1.1-1-283080). 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, USCCC. We acknowledge the financial support to the eddy covariance data harmonization provided by CarboEuropeIP, FAO-GTOS-TCO, iLEAPS, Max Planck Institute for Biogeochemistry, National Science Foundation, University of Tuscia, Université Laval and 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, University of Virginia. Philippe Ciais acknowledges support from the European Research Council through Synergy grant ERC-2013-SyG-610028 “IMBALANCE-P”. The authors wish to thank M. Jung for providing access to the GPP MTE data, which were downloaded from the GEOCARBON data portal (https://www.bgc-jena.mpg.de/geodb/projects/Data.php). The authors are also grateful to computing support and resources provided at LSCE and to the overall ORCHIDEE project that coordinate the development of the code (http://labex.ipsl.fr/orchidee/index.php/about-the-team).Peer reviewedPublisher PD

Disentangling the impacts of anthropogenic aerosols on terrestrial carbon cycle during 1850-2014

Aerosols have a dimming and cooling effect and change hydrological regimes, thus affecting carbon fluxes, which are sensitive to climate. Aerosols also scatter sunlight, which increases the fraction of diffuse radiation, increasing photosynthesis. There remains no clear conclusion whether the impact of aerosols on land carbon fluxes is larger through diffuse radiation change than through changes in other climate variables. In this study, we quantified the overall physical impacts of anthropogenic aerosols on land C fluxes and explored the contribution from each factor using a set of factorial simulations driven by climate and aerosol data from the IPSL-CM6A-LR experiments during 1850-2014. A newly-developed land surface model which distinguishes diffuse and direct radiation in canopy radiation transmission, ORCHIDEE_DF, was used. Specifically, a sub-grid scheme was developed to distinguish the cloudy and clear sky conditions. We found that anthropogenic aerosol emissions since 1850 cumulatively enhanced the land C sink by 22.6 PgC. 78% of this C sink enhancement is contributed by aerosol-induced increase in the diffuse radiation fraction, much larger than the effect of the aerosol-induced dimming. The cooling of anthropogenic aerosols has different impacts in different latitudes but overall increases the global land C sink. The dominant role of diffuse radiation changes found in this study implies that future aerosol emissions may have a much stronger impacts on the C cycle through changing radiation quality than through changing climate alone. Earth system models need to consider the diffuse radiation fertilization effect to better evaluate the impacts of climate change mitigation scenarios

Climatic and biotic factors influencing regional declines and recovery of tropical forest biomass from the 2015/16 El Niño

International audienceThe 2015/16 El Niño brought severe drought and record-breaking temperatures in the tropics. Here, using satellite-based L-band microwave vegetation optical depth, we mapped changes of above-ground biomass (AGB) during the drought and in subsequent years up to 2019. Over more than 60% of drought-affected intact forests, AGB reduced during the drought, except in the wettest part of the central Amazon, where it declined 1 y later. By the end of 2019, only 40% of AGB reduced intact forests had fully recovered to the predrought level. Using random-forest models, we found that the magnitude of AGB losses during the drought was mainly associated with regionally distinct patterns of soil water deficits and soil clay content. For the AGB recovery, we found strong influences of AGB losses during the drought and of γ . γ is a parameter related to canopy structure and is defined as the ratio of two relative height (RH) metrics of Geoscience Laser Altimeter System (GLAS) waveform data—RH25 (25% energy return height) and RH100 (100% energy return height; i.e., top canopy height). A high γ may reflect forests with a tall understory, thick and closed canopy, and/or without degradation. Such forests with a high γ ( γ ≥ 0.3) appear to have a stronger capacity to recover than low- γ ones. Our results highlight the importance of forest structure when predicting the consequences of future drought stress in the tropics

Carbonyl sulfide : comparing a mechanistic representation of the vegetation uptake in a land surface model and the leaf relative uptake approach

Land surface modellers need measurable proxies to constrain the quantity of carbon dioxide (CO2) assimilated by continental plants through photosynthesis, known as gross primary production (GPP). Carbonyl sulfide (COS), which is taken up by leaves through their stomates and then hydrolysed by photosynthetic enzymes, is a candidate GPP proxy. A former study with the ORCHIDEE land surface model used a fixed ratio of COS uptake to CO2 uptake normalised to respective ambient concentrations for each vegetation type (leaf relative uptake, LRU) to compute vegetation COS fluxes from GPP. The LRU approach is known to have limited accuracy since the LRU ratio changes with variables such as photosynthetically active radiation (PAR): while CO2 uptake slows under low light, COS uptake is not light limited. However, the LRU approach has been popular for COS-GPP proxy studies because of its ease of application and apparent low contribution to uncertainty for regional-scale applications. In this study we refined the COS-GPP relationship and implemented in ORCHIDEE a mechanistic model that describes COS uptake by continental vegetation. We compared the simulated COS fluxes against measured hourly COS fluxes at two sites and studied the model behaviour and links with environmental drivers. We performed simulations at a global scale, and we estimated the global COS uptake by vegetation to be -756 Gg S yr(-1) , in the middle range of former studies (-490 to -1335 Gg S yr(-1)). Based on monthly mean fluxes simulated by the mechanistic approach in ORCHIDEE, we derived new LRU values for the different vegetation types, ranging between 0.92 and 1.72, close to recently published averages for observed values of 1.21 for C-4 and 1.68 for C-3 plants. We transported the COS using the monthly vegetation COS fluxes derived from both the mechanistic and the LRU approaches, and we evaluated the simulated COS concentrations at NOAA sites. Although the mechanistic approach was more appropriate when comparing to high-temporal-resolution COS flux measurements, both approaches gave similar results when transporting with monthly COS fluxes and evaluating COS concentrations at stations. In our study, uncertainties between these two approaches are of secondary importance compared to the uncertainties in the COS global budget, which are currently a limiting factor to the potential of COS concentrations to constrain GPP simulated by land surface models on the global scale.Peer reviewe

Chlorophyll a fluorescence illuminates a path connecting plant molecular biology to Earth-system science

Remote sensing methods enable detection of solar-induced chlorophyll a fluorescence. However, to unleash the full potential of this signal, intensive cross-disciplinary work is required to harmonize biophysical and ecophysiological studies. For decades, the dynamic nature of chlorophyll a fluorescence (ChlaF) has provided insight into the biophysics and ecophysiology of the light reactions of photosynthesis from the subcellular to leaf scales. Recent advances in remote sensing methods enable detection of ChlaF induced by sunlight across a range of larger scales, from using instruments mounted on towers above plant canopies to Earth-orbiting satellites. This signal is referred to as solar-induced fluorescence (SIF) and its application promises to overcome spatial constraints on studies of photosynthesis, opening new research directions and opportunities in ecology, ecophysiology, biogeochemistry, agriculture and forestry. However, to unleash the full potential of SIF, intensive cross-disciplinary work is required to harmonize these new advances with the rich history of biophysical and ecophysiological studies of ChlaF, fostering the development of next-generation plant physiological and Earth-system models. Here, we introduce the scale-dependent link between SIF and photosynthesis, with an emphasis on seven remaining scientific challenges, and present a roadmap to facilitate future collaborative research towards new applications of SIF.Peer reviewe

Contrasting effects of CO₂ fertilization, land-use change and warming on seasonal amplitude of Northern Hemisphere CO₂ exchange

Continuous atmospheric CO₂ monitoring data indicate an increase in the amplitude of seasonal CO₂-cycle exchange (SCA_(NBP)) in northern high latitudes. The major drivers of enhanced SCA_(NBP) remain unclear and intensely debated, with land-use change, CO₂ fertilization and warming being identified as likely contributors. We integrated CO₂-flux data from two atmospheric inversions (consistent with atmospheric records) and from 11 state-of-the-art land-surface models (LSMs) to evaluate the relative importance of individual contributors to trends and drivers of the SCA_(NBP) of CO₂ fluxes for 1980–2015. The LSMs generally reproduce the latitudinal increase in SCA_(NBP) trends within the inversions range. Inversions and LSMs attribute SCA_(NBP) increase to boreal Asia and Europe due to enhanced vegetation productivity (in LSMs) and point to contrasting effects of CO₂ fertilization (positive) and warming (negative) on SCA_(NBP). Our results do not support land-use change as a key contributor to the increase in SCA_(NBP). The sensitivity of simulated microbial respiration to temperature in LSMs explained biases in SCA_(NBP) trends, which suggests that SCA_(NBP) could help to constrain model turnover times

ORCHIDEE-MICT (revision 4126), a land surface model for the high-latitudes: model description and validation

Abstract. The high-latitude regions of the northern hemisphere are a nexus for the interaction between land surface physical properties and their exchange of carbon and energy with the atmosphere. At these latitudes, two carbon pools of planetary significance – those of the permanently frozen soils (permafrost), and of the great expanse of boreal forest – are vulnerable to destabilization in the face of currently observed climatic warming, the speed and intensity of which are expected to increase with time. Improved projections of future Arctic and boreal ecosystem transformation require improved land surface models that integrate processes specific to these cold biomes. To this end, this study lays out relevant new parameterizations in the ORCHIDEE-MICT land surface model. These describe the interactions between soil carbon, soil temperature and hydrology, and their resulting feedbacks on water and CO2 fluxes, in addition to a recently-developed fire module. Outputs from ORCHIDEE-MICT, when forced by two climate input data sets, are extensively evaluated against: (i) temperature gradients between the atmosphere and deep soils; (ii) the hydrological components comprising the water balance of the largest high-latitude basins, and (iii) CO2 flux and carbon stock observations. The model performance is good with respect to empirical data, despite a simulated excessive plant water stress and a positive land surface temperature bias. In addition, acute model sensitivity to the choice of input forcing data suggests that the calibration of model parameters is strongly forcing-dependent. Overall, we suggest that this new model design is at the forefront of current efforts to reliably estimate future perturbations to the high-latitude terrestrial environment. </jats:p

What eddy-covariance measurements tell us about prior land flux errors in CO₂-flux inversion schemes

To guide the future development of CO₂-atmospheric inversion modeling systems, we analyzed the errors arising from prior information about terrestrial ecosystem fluxes. We compared the surface fluxes calculated by a process-based terrestrial ecosystem model with daily averages of CO₂ flux measurements at 156 sites across the world in the FLUXNET network. At the daily scale, the standard deviation of the model-data fit was 2.5 gC·m⁻²·d⁻¹; temporal autocorrelations were significant at the weekly scale (>0.3 for lags less than four weeks), while spatial correlations were confined to within the first few hundred kilometers (<0.2 after 200 km). Separating out the plant functional types did not increase the spatial correlations, except for the deciduous broad-leaved forests. Using the statistics of the flux measurements as a proxy for the statistics of the prior flux errors was shown not to be a viable approach. A statistical model allowed us to upscale the site-level flux error statistics to the coarser spatial and temporal resolutions used in regional or global models. This approach allowed us to quantify how aggregation reduces error variances, while increasing correlations. As an example, for a typical inversion of grid point (300 km × 300 km) monthly fluxes, we found that the prior flux error follows an approximate e-folding correlation length of 500 km only, with correlations from one month to the next as large as 0.6

What eddy-covariance measurements tell us about prior land flux errors in CO2-flux inversion schemes

0.2 after 200 km). Separating out the plant functional types did not increase the spatial correlations, except for the deciduous broad-leaved forests. Using the statistics of the flux measurements as a proxy for the statistics of the prior flux errors was shown not to be a viable approach. A statistical model allowed us to upscale the site-level flux error statistics to the coarser spatial and temporal resolutions used in regional or global models. This approach allowed us to quantify how aggregation reduces error variances, while increasing correlations. As an example, for a typical inversion of grid point (300 km × 300 km) monthly fluxes, we found that the prior flux error follows an approximate e-folding correlation length of 500 km only, with correlations from one month to the next as large as 0.6

Intercomparison of atmospheric Carbonyl Sulfide (TransCom-COS; Part one): Evaluating the impact of transport and emissions on tropospheric variability using ground-based and aircraft data

We present a comparison of atmospheric transport model simulations for carbonyl sulfide (COS), within the framework of the ongoing atmospheric tracer transport model intercomparison project “TransCom”. Seven atmospheric transport models participated in the inter-comparison experiment and provided simulations of COS mixing ratios in the troposphere over a 9-year period (2010–2018), using prescribed state-of-the-art surface fluxes for various components of the atmospheric COS budget: biospheric sink, oceanic source, sources from fire and industry. Since the biosphere is the largest sink of COS, we tested sink estimates produced by two different biosphere models. The main goals of TransCom-COS are (a) to investigate the impact of the transport uncertainty and emission distribution in simulating the spatio-temporal variability of COS mixing ratios in the troposphere, and (b) to assess the sensitivity of simulated tropospheric COS mixing ratios to the seasonal and diurnal variability of the COS biosphere fluxes. To this end, a control case with state-of-the-art seasonal fluxes of COS was constructed. Models were run with the same fluxes and without chemistry to isolate transport differences. Further, two COS flux scenarios were compared: one using a biosphere flux with a monthly time resolution and the other using a biosphere flux with a three-hourly time resolution. In addition, we investigated the sensitivity of the simulated concentrations to different biosphere fluxes and to indirect oceanic emissions through dimethylsulfide (DMS) and carbon disulfide (CS2). The modelled COS mixing ratios were assessed against in-situ observations from surface stations and aircraft