Intercomparison of methods to estimate gross primary production based on CO2 and COS flux measurements

Separating the components of ecosystem-scale carbon exchange is crucial in order to develop better models and future predictions of the terrestrial carbon cycle. However, there are several uncertainties and unknowns related to current photosynthesis estimates. In this study, we evaluate four different methods for estimating photosynthesis at a boreal forest at the ecosystem scale, of which two are based on carbon dioxide (CO2) flux measurements and two on carbonyl sulfide (COS) flux measurements. The CO2-based methods use traditional flux partitioning and artificial neural networks to separate the net CO2 flux into respiration and photosynthesis. The COS-based methods make use of a unique 5-year COS flux data set and involve two different approaches to determine the leaf-scale relative uptake ratio of COS and CO2 (LRU), of which one (LRUCAP) was developed in this study. LRUCAP was based on a previously tested stomatal optimization theory (CAP), while LRUPAR was based on an empirical relation to measured radiation. For the measurement period 2013-2017, the artificial neural network method gave a GPP estimate very close to that of traditional flux partitioning at all timescales. On average, the COS-based methods gave higher GPP estimates than the CO2-based estimates on daily (23 % and 7 % higher, using LRUPAR and LRUCAP, respectively) and monthly scales (20 % and 3 % higher), as well as a higher cumulative sum over 3 months in all years (on average 25 % and 3 % higher). LRUCAP was higher than LRU estimated from chamber measurements at high radiation, leading to underestimation of midday GPP relative to other GPP methods. In general, however, use of LRUCAP gave closer agreement with CO2-based estimates of GPP than use of LRUPAR. When extended to other sites, LRUCAP may be more robust than LRUPAR because it is based on a physiological model whose parameters can be estimated from simple measurements or obtained from the literature. In contrast, the empirical radiation relation in LRUPAR may be more site-specific. However, this requires further testing at other measurement sites.Peer reviewe

A carbon sink-driven approach to estimate gross primary production from microwave satellite observations

Global estimation of Gross Primary Production (GPP) - the uptake of atmospheric carbon dioxide by plants through photosynthesis - is commonly based on optical satellite remote sensing data. This presents a source-driven approach since it uses the amount of absorbed light, the main driver of photosynthesis, as a proxy for GPP. Vegetation Optical Depth (VOD) estimates obtained from microwave sensors provide an alternative and independent data source to estimate GPP on a global scale, which may complement existing GPP products. Recent studies have shown that VOD is related to aboveground biomass, and that both VOD and temporal changes in VOD relate to GPP. In this study, we build upon this concept and propose a model for estimating GPP from VOD. Since the model is driven by vegetation biomass, as observed through VOD, it presents a carbon sink-driven approach to quantify GPP and, therefore, is conceptually different from common source-driven approaches. The model developed in this study uses single frequencies from active or passive microwave VOD retrievals from C-, X- and Ku-band (Advanced Scatterometer (ASCAT) and Advanced Microwave Scanning Radiometer for Earth Observation (AMSR-E)) to estimate GPP at the global scale. We assessed the ability for temporal and spatial extrapolation of the model using global GPP from FLUXCOM and in situ GPP from FLUXNET. We further performed upscaling of in situ GPP based on different VOD data sets and compared these estimates with the FLUXCOM and MODerate-resolution Imaging Spectroradiometer (MODIS) GPP products. Our results show that the model developed for individual grid cells using VOD and change in VOD as input performs well in predicting temporal patterns in GPP for all VOD data sets. For spatial extrapolation of the model, however, additional input variables are needed to represent the spatial variability of the VOD-GPP relationship due to differences in vegetation type. As additional input variable, we included the grid cell median VOD (as a proxy for vegetation cover), which increased the model performance during cross validation. Mean annual GPP obtained for AMSR-E X-band data tends to overestimate mean annual GPP for FLUXCOM and MODIS but shows comparable latitudinal patterns. Overall, our findings demonstrate the potential of VOD for estimating GPP. The sink-driven approach provides additional information about GPP independent of optical data, which may contribute to our knowledge about the carbon source-sink balance in different ecosystems

Scaling carbon fluxes from eddy covariance sites to globe: Synthesis and evaluation of the FLUXCOM approach

FLUXNET assembles globally-distributed eddy covariance-based estimates of carbon fluxes between the biosphere and the atmosphere. Since eddy covariance flux towers have a relatively small footprint and are distributed unevenly across the world, upscaling the observations is necessary in order to obtain global-scale estimates of biosphere-atmosphere exchange from the flux tower network. Based on cross-consistency checks with atmospheric inversions, sun-induced fluorescence (SIF) and dynamic global vegetation models (DGVM), we provide here a systematic assessment of the latest upscaling efforts for gross primary production (GPP) and net ecosystem exchange (NEE) of the FLUXCOM initiative, where different machine learning methods, forcing datasets, and sets of predictor variables were employed. Spatial patterns of mean GPP are consistent among FLUXCOM and DGVM ensembles (R2 > 0.94 at 1° spatial resolution) while the majority of DGVMs are outside the FLUXCOM range for 70 % of the land surface. Global mean GPP magnitudes for 2008–2010 from FLUXCOM members vary within 106 and 130 PgC yr−1 with the largest uncertainty in the tropics. Seasonal variations of independent SIF estimates agree better with FLUXCOM GPP (mean global pixel-wise R2 ~ 0.75) than with GPP from DGVMs (mean global pixel wise R2 ~ 0.6). Seasonal variations of FLUXCOM NEE show good consistency with atmospheric inversion-based net land carbon fluxes, particularly for temperate and boreal regions (R2 > 0.92). Interannual variability of global NEE in FLUXCOM is underestimated compared to inversions and DGVMs. The FLUXCOM version which uses also meteorological inputs shows a strong co-variation of interannual patterns with inversions (R2 = 0.88 for 2001–2010). Mean regional NEE from FLUXCOM shows larger uptake than inversion and DGVM-based estimates, particularly in the tropics with discrepancies of up to several hundred gC m2 yr−1. These discrepancies can only partly be reconciled by carbon loss pathways that are implicit in inversions but not captured by the flux tower measurements such as carbon emissions from fires and water bodies. We hypothesize that a combination of systematic biases in the underlying eddy covariance data, in particular in tall tropical forests, and a lack of site-history effects on NEE in FLUXCOM are likely responsible for the too strong tropical carbon sink estimated by FLUXCOM. Furthermore, as FLUXCOM does not account for CO2 fertilization effects carbon flux trends are not realistic. Overall, current FLUXCOM estimates of mean annual and seasonal cycles of GPP as well as seasonal NEE variations provide useful constraints of global carbon cycling, while interannual variability patterns from FLUXCOM are valuable but require cautious interpretation. Exploring the diversity of Earth Observation data and of machine learning concepts along with improved quality and quantity of flux tower measurements will facilitate further improvements of the FLUXCOM approach overall

Vegetation optical depth at L-band and above ground biomass in the tropical range: Evaluating their relationships at continental and regional scales

Abstract The relationship between vegetation optical depth (VOD) retrieved by L-band SMOS radiometer and forest above ground biomass (AGB) was investigated in tropical areas of Africa and South America. VOD was retrieved from the latest version of level 2 SMOS algorithm, while reference AGB was obtained from a pantropical database, encompassing a large number of ground plot data derived from field surveys conducted on both continents. In Africa and South-America, VOD increased with AGB, reaching saturation at about 350 Mg ha−1. The strength of the relation was improved selecting VOD data in appropriate seasons, characterized by a higher dynamic range of values. The capability of VOD data to estimate AGB was further investigated using Random Forest decision trees, adding to VOD selected climate variables from the Climatic Research Unit (temperature, potential evapotranspiration, and precipitation) and water deficit data, and validating regression tests with ground data from the reference AGB database. The results for the five analyzed years indicate that the best estimates of AGB are obtained by the joined use of VOD and potential evapotranspiration input data, but all climate variables brought an improvement in AGB estimates. AGB estimates were relatively stable for the considered period, with limited variations possibly due to changes in biomass and to data quality of VOD and of climate variables. The VOD signal and estimated AGB were also analyzed according to ecological homogeneous units (ecoregions), evidencing data clusters, partially overlapped to each other, in the VOD - AGB plane

Assessing the relationship between microwave vegetation optical depth and gross primary production

At the global scale, the uptake of atmospheric carbon dioxide by terrestrial ecosystems through photosynthesis is commonly estimated through vegetation indices or biophysical properties derived from optical remote sensing data. Microwave observations of vegetated areas are sensitive to different components of the vegetation layer than observations in the optical domain and may therefore provide complementary information on the vegetation state, which may be used in the estimation of Gross Primary Production (GPP). However, the relation between GPP and Vegetation Optical Depth (VOD), a biophysical quantity derived from microwave observations, is not yet known. This study aims to explore the relationship between VOD and GPP. VOD data were taken from different frequencies (L-, C-, and X-band) and from both active and passive microwave sensors, including the Advanced Scatterometer (ASCAT), the Soil Moisture Ocean Salinity (SMOS) mission, the Advanced Microwave Scanning Radiometer for Earth Observation System (AMSR-E) and a merged VOD data set from various passive microwave sensors. VOD data were compared against FLUXCOM GPP and Solar-Induced chlorophyll Fluorescence (SIF) from the Global Ozone Monitoring Experiment-2 (GOME-2). FLUXCOM GPP estimates are based on the upscaling of flux tower GPP observations using optical satellite data, while SIF observations present a measure of photosynthetic activity and are often used as a proxy for GPP. For relating VOD to GPP, three variables were analyzed: original VOD time series, temporal changes in VOD (ΔVOD), and positive changes in VOD (ΔVOD≥0). Results show widespread positive correlations between VOD and GPP with some negative correlations mainly occurring in dry and wet regions for active and passive VOD, respectively. Correlations between VOD and GPP were similar or higher than between VOD and SIF. When comparing the three variables for relating VOD to GPP, correlations with GPP were higher for the original VOD time series than for ΔVOD or ΔVOD≥0 in case of sparsely to moderately vegetated areas and evergreen forests, while the opposite was true for deciduous forests. Results suggest that original VOD time series should be used jointly with changes in VOD for the estimation of GPP across biomes, which may further benefit from combining active and passive VOD data