Search

Wood density is a crucial variable linked to mechanical, physiological, and ecological properties. In this study, we analyzed an extensive dataset of over 48,000 wood density samples collected from 2,920 trees. Our aim was to explore variations in wood density, at both inter-tree and intra-tree levels, along with the factors contributing to these variations. Inter-tree variations reveal significant differences in wood density among eight dominant species, highlighting their role in shaping wood density. As tree species exhibit specific spatial distributions associated with microhabitats, we anticipated a link between wood density distribution and microhabitat. Using a feature selection approach and random forest model, we identified six predictors, including satellite-based vegetation indexes, topographic variables, and soil sand content, capable of predicting 91 % of spatial wood density variations. The Normalized Difference Vegetation Index (NDVI) positively representing the amount of carbon within trees correlated with wood density, while the Normalized Difference Water Index (NDWI), reflecting water content, and soil sand content showed negative associations. Geomorphons and soil sand context provided insights into wood density variations and specific landforms. Lower wood density values were linked to landforms with low geomorphons (summit, ridge, or shoulder), whereas higher wood density was found in landforms with high geomorphons (valley, depression, or hollow areas). Furthermore, our study highlighted the importance of considering intra-tree variation, a facet often overlooked in previous research. Interestingly, the magnitude of intra-tree variation is comparable to, and in some species even exceeds, that of inter-tree variations. The intra-tree wood density samples display significant differences both vertically along the height and radially from the center to the bark zones of trees. These variations are influenced by tree growing strategy, living conditions, and physiological structure. In summary, our research delved into the multifaceted features of wood density, shedding light on critical aspects of this fundamental variable. [ABSTRACT FROM AUTHOR]

Wood density is a crucial variable linked to mechanical, physiological, and ecological properties. In this study, we analyzed an extensive dataset of over 48,000 wood density samples collected from 2,920 trees. Our aim was to explore variations in wood density, at both inter-tree and intra-tree levels, along with the factors contributing to these variations. Intertree variations reveal significant differences in wood density among eight dominant species, highlighting their role in shaping wood density. As tree species exhibit specific spatial distributions associated with microhabitats, we anticipated a link between wood density distribution and microhabitat. Using a feature selection approach and random forest model, we identified six predictors, including satellite-based vegetation indexes, topographic variables, and soil sand content, capable of predicting 91% of spatial wood density variations. The Normalized Difference Vegetation Index (NDVI) positively representing the amount of carbon within trees correlated with wood density, while the Normalized Difference Water Index (NDWI), reflecting water content, and soil sand content showed negative associations. Geomorphons and soil sand context provided insights into wood density variations and specific landforms. Lower wood density values were linked to landforms with low geomorphons (summit, ridge, or shoulder), whereas higher wood density was found in landforms with high geomorphons (valley, depression, or hollow areas). Furthermore, our study highlighted the importance of considering intra-tree variation, a facet often overlooked in previous research. Interestingly, the magnitude of intra-tree variation is comparable to, and in some species even exceeds, that of inter-tree variations. The intra-tree wood density samples display significant differences both vertically along the height and radially from the center to the bark zones of trees. These variations are influenced by tree growing strategy, living conditions, and physiological structure. In summary, our research delved into the multifaceted features of wood density, shedding light on critical aspects of this fundamental variable. [ABSTRACT FROM AUTHOR]

Terrestrial water storage (TWS) is an integrative hydrological state that is key for our understanding of the global water cycle. The TWS observation from the GRACE missions has, therefore, been instrumental in the calibration and validation of hydrological models and understanding the variations in the hydrological storage. The models, however, still show significant uncertainties in reproducing observed TWS variations, especially for the interannual variability (IAV) at the global scale. Here, we diagnose the regions dominating the variance in globally integrated TWS IAV and the sources of the errors in two data-driven hydrological models that were calibrated against global TWS, snow water equivalent, evapotranspiration, and runoff data. We used (1) a parsimonious process-based hydrological model, the Strategies to INtegrate Data and BiogeochemicAl moDels (SINDBAD) framework and (2) a machine learning, physically based hybrid hydrological model (H2M) that combines a dynamic neural network with a water balance concept. While both models agree with the Gravity Recovery and Climate Experiment (GRACE) that global TWS IAV is largely driven by the semi-arid regions of southern Africa, the Indian subcontinent and northern Australia, and the humid regions of northern South America and the Mekong River basin, the models still show errors such as the overestimation of the observed magnitude of TWS IAV at the global scale. Our analysis identifies modeling error hotspots of the global TWS IAV, mostly in the tropical regions including the Amazon, sub-Saharan regions, and Southeast Asia, indicating that the regions that dominate global TWS IAV are not necessarily the same as those that dominate the error in global TWS IAV. Excluding those error hotspot regions in the global integration yields large improvements in the simulated global TWS IAV, which implies that model improvements can focus on improving processes in these hotspot regions. Further analysis indicates that error hotspot regions are associated with lateral flow dynamics, including both sub-pixel moisture convergence and across-pixel lateral river flow, or with interactions between surface processes and groundwater. The association of model deficiencies with land processes that delay the TWS variation could, in part, explain why the models cannot represent the observed lagged response of TWS IAV to precipitation IAV in hotspot regions that manifest as errors in global TWS IAV. Our approach presents a general avenue to better diagnose model simulation errors for global data streams to guide efficient and focused model development for regions and processes that matter the most. [ABSTRACT FROM AUTHOR]

Terrestrial water storage (TWS) is an integrative hydrological state that is key for our understanding of the global water cycle. The TWS observation from GRACE missions has, therefore, been instrumental in calibration and validation of hydrological models and understanding the variations of the hydrological storages. The models, however, still show significant uncertainties in reproducing observed TWS variations, especially for the interannual variability (IAV) at the global scale. Here, we diagnose the regions dominating the variance of globally integrated TWS IAV, and sources of the errors in two data-driven hydrological models that were calibrated against global TWS, snow water equivalent, evapotranspiration, and runoff data: 1) a parsimonious process-based hydrological model, the Strategies to INtegrate Data and BiogeochemicAl moDels (SINDBAD) framework, and 2) a machine learning-physically based hybrid hydrological model (H2M) that combines a dynamic neural network with a water balance concept. While both models agree with GRACE that global TWS IAV is largely driven by the semi-arid regions of southern Africa, Indian subcontinent and northern Australia, and the humid regions of northern South America and Mekong River Basin, the models still show errors such as overestimation of the observed magnitude of TWS IAV at the global scale. Our analysis identifies modeling error hotspots of the global TWS IAV mostly in the tropical regions including Amazon, sub-Saharan regions, and Southeast Asia, indicating that the regions that dominate global TWS IAV are not necessarily the same as those that dom15 inate the error in global TWS IAV. Excluding those error hotspot regions in the global integration yields large improvements of simulated global TWS IAV, which implies that model improvements can focus on improving processes in these hotspot regions. Further analysis indicates that error hotspot regions are associated with lateral flow dynamics, including both sub-pixel moisture convergence and across pixel lateral river flow, or with interactions between surface processes and groundwater. The association of model deficiencies with land processes that delay the TWS variation could, in part, explain why the models cannot represent the observed lagged response of TWS IAV to precipitation IAV in hotspot regions that manifest to errors in global TWS IAV. Our approach presents a general avenue to better diagnose model simulation errors for global data streams to guide efficient and focused model development for regions and processes that matter the most. [ABSTRACT FROM AUTHOR]

The eddy-covariance technique measures carbon, water, and energy fluxes between the land surface and the atmosphere at hundreds of sites globally. Collections of standardised and homogenised flux estimates such as the LaThuile , Fluxnet2015 , National Ecological Observatory Network (NEON), Integrated Carbon Observation System (ICOS), AsiaFlux, AmeriFlux, and Terrestrial Ecosystem Research Network (TERN)/OzFlux data sets are invaluable to study land surface processes and vegetation functioning at the ecosystem scale. Space-borne measurements give complementary information on the state of the land surface in the surroundings of the towers. They aid the interpretation of the fluxes and support the benchmarking of terrestrial biosphere models. However, insufficient quality and frequent and/or long gaps are recurrent problems in applying the remotely sensed data and may considerably affect the scientific conclusions. Here, we describe a standardised procedure to extract, quality filter, and gap-fill Earth observation data from the MODIS instruments and the Landsat satellites. The methods consistently process surface reflectance in individual spectral bands, derived vegetation indices, and land surface temperature. A geometrical correction estimates the magnitude of land surface temperature as if seen from nadir or 40 ∘ off-nadir. Finally, we offer the community living data sets of pre-processed Earth observation data, where version 1.0 features the MCD43A4/A2 and MxD11A1 MODIS products and Landsat Collection 1 Tier 1 and Tier 2 products in a radius of 2 km around 338 flux sites. The data sets we provide can widely facilitate the integration of activities in the eddy-covariance, remote sensing, and modelling fields. [ABSTRACT FROM AUTHOR]

So far, various studies have aimed at decomposing the integrated terrestrial water storage variations observed by satellite gravimetry (GRACE, GRACE-FO) with the help of large-scale hydrological models. While the results of the storage decomposition depend on model structure, little attention has been given to the impact of the way that vegetation is represented in these models. Although vegetation structure and activity represent the crucial link between water, carbon, and energy cycles, their representation in large-scale hydrological models remains a major source of uncertainty. At the same time, the increasing availability and quality of Earth-observation-based vegetation data provide valuable information with good prospects for improving model simulations and gaining better insights into the role of vegetation within the global water cycle. In this study, we use observation-based vegetation information such as vegetation indices and rooting depths for spatializing the parameters of a simple global hydrological model to define infiltration, root water uptake, and transpiration processes. The parameters are further constrained by considering observations of terrestrial water storage anomalies (TWS), soil moisture, evapotranspiration (ET) and gridded runoff (Q) estimates in a multi-criteria calibration approach. We assess the implications of including varying vegetation characteristics on the simulation results, with a particular focus on the partitioning between water storage components. To isolate the effect of vegetation, we compare a model experiment in which vegetation parameters vary in space and time to a baseline experiment in which all parameters are calibrated as static, globally uniform values. Both experiments show good overall performance, but explicitly including varying vegetation data leads to even better performance and more physically plausible parameter values. The largest improvements regarding TWS and ET are seen in supply-limited (semi-arid) regions and in the tropics, whereas Q simulations improve mainly in northern latitudes. While the total fluxes and storages are similar, accounting for vegetation substantially changes the contributions of different soil water storage components to the TWS variations. This suggests an important role of the representation of vegetation in hydrological models for interpreting TWS variations. Our simulations further indicate a major effect of deeper moisture storages and groundwater–soil moisture–vegetation interactions as a key to understanding TWS variations. We highlight the need for further observations to identify the adequate model structure rather than only model parameters for a reasonable representation and interpretation of vegetation–water interactions. [ABSTRACT FROM AUTHOR]

Forest age can determine the capacity of a forest to uptake carbon from the atmosphere. However, a lack of global diagnostics that reflect the forest stage and associated disturbance regimes hampers the quantification of age-related differences in forest carbon dynamics. This study provides a new global distribution of forest age circa 2010, estimated using a machine learning approach trained with more than 40 000 plots using forest inventory, biomass and climate data. First, an evaluation against the plot-level measurements of forest age reveals that the data-driven method has a relatively good predictive capacity of classifying old-growth vs. non-old-growth (precision = 0.81 and 0.99 for old-growth and non-old-growth, respectively) forests and estimating corresponding forest age estimates (NSE = 0.6 – Nash–Sutcliffe efficiency – and RMSE = 50 years – root-mean-square error). However, there are systematic biases of overestimation in young- and underestimation in old-forest stands, respectively. Globally, we find a large variability in forest age with the old-growth forests in the tropical regions of Amazon and Congo, young forests in China, and intermediate stands in Europe. Furthermore, we find that the regions with high rates of deforestation or forest degradation (e.g. the arc of deforestation in the Amazon) are composed mainly of younger stands. Assessment of forest age in the climate space shows that the old forests are either in cold and dry regions or warm and wet regions, while young–intermediate forests span a large climatic gradient. Finally, comparing the presented forest age estimates with a series of regional products reveals differences rooted in different approaches and different in situ observations and global-scale products. Despite showing robustness in cross-validation results, additional methodological insights on further developments should as much as possible harmonize data across the different approaches. The forest age dataset presented here provides additional insights into the global distribution of forest age to better understand the global dynamics in the forest water and carbon cycles. The forest age datasets are openly available at 10.17871/ForestAgeBGI.2021 (Besnard et al., 2021). [ABSTRACT FROM AUTHOR]

The terrestrial forest carbon pool is poorly quantified, in particular in regions with low forest inventory capacity. By combining multiple satellite observations of synthetic aperture radar (SAR) backscatter around the year 2010, we generated a global, spatially explicit dataset of above-ground live biomass (AGB; dry mass) stored in forests with a spatial resolution of 1 ha. Using an extensive database of 110 897 AGB measurements from field inventory plots, we show that the spatial patterns and magnitude of AGB are well captured in our map with the exception of regional uncertainties in high-carbon-stock forests with AGB >250 Mgha-1 , where the retrieval was effectively based on a single radar observation. With a total global AGB of 522 Pg , our estimate of the terrestrial biomass pool in forests is lower than most estimates published in the literature (426–571 Pg). Nonetheless, our dataset increases knowledge on the spatial distribution of AGB compared to the Global Forest Resources Assessment (FRA) by the Food and Agriculture Organization (FAO) and highlights the impact of a country's national inventory capacity on the accuracy of the biomass statistics reported to the FRA. We also reassessed previous remote sensing AGB maps and identified major biases compared to inventory data, up to 120 % of the inventory value in dry tropical forests, in the subtropics and temperate zone. Because of the high level of detail and the overall reliability of the AGB spatial patterns, our global dataset of AGB is likely to have significant impacts on climate, carbon, and socio-economic modelling schemes and provides a crucial baseline in future carbon stock change estimates. The dataset is available at 10.1594/PANGAEA.894711 (Santoro, 2018). [ABSTRACT FROM AUTHOR]

So far, various studies aimed at decomposing the integrated terrestrial water storage variations observed by satellite gravimetry (GRACE, GRACE-FO) with the help of large-scale hydrological models. While the results of the storage decomposition depend on model structure, little attention has been given to the impact of the way how vegetation is represented in these models. Although vegetation structure and activity represent the crucial link between water, carbon and energy cycles, their representation in large-scale hydrological models remains a major source of uncertainty. At the same time, the increasing availability and quality of Earth observation-based vegetation data provide valuable information with good prospects for improving model simulations and gaining better insights into the role of vegetation within the global water cycle. In this study, we use observation-based vegetation information such as vegetation indices and rooting depths for spatializing the parameters of a simple global hydrological model to define infiltration, root water uptake and transpiration processes. The parameters are further constrained by considering observations of terrestrial water storage anomalies (TWS), soil moisture, evapotranspiration (ET) and gridded runoff (Q) estimates in a multi-criteria calibration approach. We assess the implications of including vegetation on the simulation results, with a particular focus on the partitioning between water storage components. To isolate the effect of vegetation, we compare a model experiment with vegetation parameters varying in space and time to a baseline experiment in which all parameters are calibrated as static, globally uniform values. Both experiments show good overall performance, but including vegetation data led to even better performance and more physically plausible parameter values. Largest improvements regarding TWS and ET were seen in supply-limited (semi-arid) regions and in the tropics, whereas Q simulations improve mainly in northern latitudes. While the total fluxes and storages are similar, accounting for vegetation substantially changes the contributions of snow and different soil water storage components to the TWS variations, with the dominance of an intermediate water pool that interacts with the fast plant accessible soil moisture and the delayed water storage. The findings indicate the important role of deeper moisture storages as well as groundwater-soil moisture-vegetation interactions as a key to understanding TWS variations. We highlight the need for further observations to identify the adequate model structure rather than only model parameters for a reasonable representation and interpretation of vegetation-water interactions. [ABSTRACT FROM AUTHOR]

The turnover time of terrestrial ecosystem carbon is an emergent ecosystem property that quantifies the strength of land surface on the global carbon cycle–climate feedback. However, observation- and modeling-based estimates of carbon turnover and its response to climate are still characterized by large uncertainties. In this study, by assessing the apparent whole ecosystem carbon turnover times (τ) as the ratio between carbon stocks and fluxes, we provide an update of this ecosystem level diagnostic and its associated uncertainties in high spatial resolution (0.083 ∘) using multiple, state-of-the-art, observation-based datasets of soil organic carbon stock (Csoil), vegetation biomass (Cveg) and gross primary productivity (GPP). Using this new ensemble of data, we estimated the global median τ to be 43-7+7 yr (median-difference to percentile 25+difference to percentile 75) when the full soil is considered, in contrast to limiting it to 1 m depth. Only considering the top 1 m of soil carbon in circumpolar regions (assuming maximum active layer depth is up to 1 m) yields a global median τ of 37-6+3 yr , which is longer than the previous estimates of 23-4+7 yr (Carvalhais et al., 2014). We show that the difference is mostly attributed to changes in global Csoil estimates. Csoil accounts for approximately 84 % of the total uncertainty in global τ estimates; GPP also contributes significantly (15 %), whereas Cveg contributes only marginally (less than 1 %) to the total uncertainty. The high uncertainty in Csoil is reflected in the large range across state-of-the-art data products, in which full-depth Csoil spans between 3362 and 4792 PgC. The uncertainty is especially high in circumpolar regions with an uncertainty of 50 % and a low spatial correlation between the different datasets (0.2

The terrestrial forest carbon pool is poorly quantified, in particular in regions with low forest inventory capacity. By combining multiple satellite observations of synthetic aperture radar (SAR) backscatter around the year 2010, we generated a global, spatially explicit dataset of above-ground forest biomass (dry mass, AGB) with a spatial resolution of 1 ha. Using an extensive database of 110,897 AGB measurements from field inventory plots, we show that the spatial patterns and magnitude of AGB are well captured in our map with the exception of regional uncertainties in high carbon stock forests with AGB > 250 Mg ha−1 where the retrieval was effectively based on a single radar observation. With a total global AGB of 522 Pg, our estimate of the terrestrial biomass pool in forests is lower than most estimates published in literature (426-571 Pg). Nonetheless, our dataset increases knowledge on the spatial distribution of AGB compared to the global Forest Resources Assessment (FRA) by the Food and Agriculture Organization (FAO) and highlights the impact of a country's national inventory capacity on the accuracy of the biomass statistics reported to the FRA. We also reassessed previous remote sensing AGB maps, and identify major biases compared to inventory data, up to 120 % of the inventory value in dry tropical forests, in the sub-tropics and temperate zone. Because of the high level of detail and the overall reliability of the AGB spatial patterns, our global dataset of AGB is likely to have significant impacts on climate, carbon and socio-economic modelling schemes, and provides a crucial baseline in future carbon stock changes estimates. The dataset is available at: https://doi.pangaea.de/10.1594/PANGAEA.894711 (Santoro, 2018). [ABSTRACT FROM AUTHOR]

FLUXNET comprises 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 to obtain global-scale estimates of biosphere–atmosphere exchange. Based on cross-consistency checks with atmospheric inversions, sun-induced fluorescence (SIF) and dynamic global vegetation models (DGVMs), here we provide 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 data sets and sets of predictor variables were employed. Spatial patterns of mean GPP are consistent across FLUXCOM and DGVM ensembles (R2>0.94 at 1 ∘ spatial resolution) while the majority of DGVMs show, for 70 % of the land surface, values outside the FLUXCOM range. 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 in 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 in 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 also uses meteorological inputs shows a strong co-variation in interannual patterns with inversions (R2=0.87 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 grammes of carbon per square metre per year. 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. [ABSTRACT FROM AUTHOR]

Climate variables carry signatures of variability at multiple timescales. How these modes of variability are reflected in the state of the terrestrial biosphere is still not quantified or discussed at the global scale. Here, we set out to gain a global understanding of the relevance of different modes of variability in vegetation greenness and its covariability with climate. We used >30 years of remote sensing records of the normalized difference vegetation index (NDVI) to characterize biosphere variability across timescales from submonthly oscillations to decadal trends using discrete Fourier decomposition. Climate data of air temperature (Tair) and precipitation (Prec) were used to characterize atmosphere–biosphere covariability at each timescale. Our results show that short-term (intra-annual) and longer-term (interannual and longer) modes of variability make regionally highly important contributions to NDVI variability: short-term oscillations focus in the tropics where they shape 27 % of NDVI variability. Longer-term oscillations shape 9 % of NDVI variability, dominantly in semiarid shrublands. Assessing dominant timescales of vegetation–climate covariation, a natural surface classification emerges which captures patterns not represented by conventional classifications, especially in the tropics. Finally, we find that correlations between variables can differ and even invert signs across timescales. For southern Africa for example, correlation between NDVI and Tair is positive for the seasonal signal but negative for short-term and longer-term oscillations, indicating that both short- and long-term temperature anomalies can induce stress on vegetation dynamics. Such contrasting correlations between timescales exist for 15 % of vegetated areas for NDVI with Tair and 27 % with Prec, indicating global relevance of scale-specific climate sensitivities. Our analysis provides a detailed picture of vegetation–climate covariability globally, characterizing ecosystems by their intrinsic modes of temporal variability. We find that (i) correlations of NDVI with climate can differ between scales, (ii) nondominant subsignals in climate variables may dominate the biospheric response, and (iii) possible links may exist between short-term and longer-term scales. These heterogeneous ecosystem responses on different timescales may depend on climate zone and vegetation type, and they are to date not well understood and do not always correspond to transitions in dominant vegetation types. These scale dependencies can be a benchmark for vegetation model evaluation and for comparing remote sensing products. [ABSTRACT FROM AUTHOR]

The Earth System Model Evaluation Tool (ESMValTool) is a community diagnostics and performance metrics tool designed to improve comprehensive and routine evaluation of Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project (CMIP). It has undergone rapid development since the first release in 2016 and is now a well-tested tool that provides end-to-end provenance tracking to ensure reproducibility. It consists of an easy-to-install, well documented Python package providing the core functionalities (ESMValCore) that performs common pre-processing operations and a diagnostic part that includes tailored diagnostics and performance metrics for specific scientific applications. Here we describe large-scale diagnostics of the second major release of the tool that supports the evaluation of ESMs participating in CMIP Phase 6 (CMIP6). ESMValTool v2.0 includes a large collection of diagnostics and performance metrics for atmospheric, oceanic, and terrestrial variables for the mean state, trends, and variability. ESMValTool v2.0 also successfully reproduces figures from the evaluation and projections chapters of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5) and incorporates updates from targeted analysis packages, such as the NCAR Climate Variability Diagnostics Package for the evaluation of modes of variability the Thermodynamic Diagnostic Tool (TheDiaTo) to evaluate the energetics of the climate system, as well as parts of AutoAssess that contains a mix of top-down performance metrics. The tool has been fully integrated into the Earth System Grid Federation (ESGF) infrastructure at the Deutsches Klima Rechenzentrum (DKRZ) to provide evaluation results from CMIP6 model simulations shortly after the output is published to the CMIP archive. A result browser has been implemented that enables advanced monitoring of the evaluation results by a broad user community at much faster timescales than what was possible in CMIP5. [ABSTRACT FROM AUTHOR]

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 (R² > 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 R² ~ 0.75) than with GPP from DGVMs (mean global pixel wise R² ~ 0.6). Seasonal variations of FLUXCOM NEE show good consistency with atmospheric inversion-based net land carbon fluxes, particularly for temperate and boreal regions (R² > 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 (R² = 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 m² 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. [ABSTRACT FROM AUTHOR]

Climate variables carry signatures of variability at multiple time scales. How these modes of variability are reflected in the state of the terrestrial biosphere is still not quantified, nor discussed at the global scale. Here, we set out to gain a global understanding of the relevance of different modes of variability in vegetation greenness and its co-variability with climate. We used > 30 years of remote sensing records of Normalized Difference Vegetation Index (NDVI) to characterize biosphere variability across time scales from sub-monthly oscillations to decadal trends using discrete Fourier decomposition. Climate data of air temperature (Tair) and precipitation (Prec) were used to characterize atmosphere-biosphere co-variability at each time scale. Our results show that short-term (intra-annual) and longer-term (inter-annual and longer) modes of variability make regionally highly important contributions to NDVI variability: Short-term oscillations focus in the tropics where they shape 27 % of NDVI variability. Longer-term oscillations shape 9 % of NDVI variability, dominantly in semi-arid shrublands. Assessing dominant time scales of vegetation-climate co-variation, a natural surface classification emerges which captures patterns not represented by conventional classifications, especially in the tropics. Finally, we find that correlations between variables can differ and even invert signs across time scales. For southern Africa for example, correlation between NDVI and Tair is positive for the seasonal signal, but negative for short-term and longer-term oscillations, indicating that both short and long-term temperature anomalies can induce stress on vegetation dynamics. Such contrasting correlations between time scales exist for 15 % of vegetated area for NDVI with Tair, and 27 % with Prec, indicating global relevance of scale-specific climate sensitivities. Our analysis provides a detailed picture of vegetation-climate co-variability globally, characterizing ecosystems by their intrinsic modes of temporal variability. We find that (i) correlations of NDVI with climate can differ between scales, (ii) non-dominant sub-signals in climate variables may dominate the biospheric response, and (iii) possible links may exist between short-term and longer-term scales. These heterogeneous ecosystem responses on different time scales may depend on climate zone and vegetation type, and are to date not well understood, nor always correspond to transitions in dominant vegetation types. These scale dependencies can be a benchmark for vegetation model evaluation and for comparing remote sensing products. [ABSTRACT FROM AUTHOR]

Water-use efficiency (WUE), defined as the ratio of carbon assimilation over evapotranspiration (ET), is a key metric to assess ecosystem functioning in response to environmental conditions. It remains unclear which factors control this ratio during periods of extended water limitation. Here, we used dry-down events occurring at eddy-covariance flux tower sites in the FLUXNET database as natural experiments to assess if and how decreasing soil-water availability modifies WUE at ecosystem scale. WUE models were evaluated by their performance to predict ET from both the gross primary productivity (GPP), which characterizes carbon assimilation at ecosystem scale, and environmental variables. We first compared two water-use efficiency models: the first was based on the concept of a constant underlying water-use efficiency, and the second augmented the first with a previously detected direct influence of radiation on transpiration. Both models predicting ET strictly from atmospheric covariates failed to reproduce observed ET dynamics for these periods, as they did not explicitly account for the effect of soil-water limitation. We demonstrate that an ET-attenuating soil-water-availability factor in junction with the additional radiation term was necessary to accurately predict ET flux magnitudes and dry-down lengths of these water-limited periods. In an analysis of the attenuation of ET for the 31 included FLUXNET sites, up to 50 % of the observed decline in ET was due to the soil-water-availability effect we identified in this study. We conclude by noting that the rates of ET decline differ significantly between sites with different vegetation and climate types and discuss the dependency of this rate on the variability of seasonal dryness. [ABSTRACT FROM AUTHOR]

Vegetation carbon use efficiency (CUE) is a key measure of carbon (C) transfer from the atmosphere to terrestrial biomass, and indirectly reflects how much C is released through autotrophic respiration from the vegetation to the atmosphere. Diagnosing the variability of CUE with climate and other environmental factors is fundamental to understand its driving factors, and to further fill the current gaps in knowledge about the environmental controls on CUE. Thus, to study CUE variability and its driving factors, this study established a global database of site-year CUE based on observations from 188 field measurement sites for five ecosystem types - forest, grass, wetland, crop and tundra. The spatial pattern of CUE was predicted from global climate and soil variables using Random Forest, and compared with estimates from Dynamic Global Vegetation Models (DGVMs) from the TRENDY model ensemble. Globally, we found two prominent CUE gradients in ecosystem types and latitude, that is, CUE varied with ecosystem types, being the highest in wetlands and lowest in grassland, and CUE decreased with latitude with the lowest CUE in tropics, and the highest CUE in higher latitude regions. CUE varied greatly between data-derived CUE and TRENDY-CUE, but also among TRENDY models. Both data-derived and TRENDY-CUE challenged the constant value of 0.5 for CUE, independent of environmental controls. However, given the role of CUE in controlling the spatial and temporal variability of the terrestrial biosphere C cycle, these results emphasize the need to better understand the biotic and abiotic controls on CUE to reduce the uncertainties in prognostic land-process model simulations. Finally, this study proposed a new estimate of net primary production based on CUE and gross primary production, offering another benchmark for net primary production comparison for global carbon modelling. [ABSTRACT FROM AUTHOR]

Water-use efficiency, defined as the ratio of carbon assimilation over evapotranspiration (ET), is a key metric to assess ecosystem functioning in response to environmental conditions. It remains unclear which factors control this ratio during periods of extended water-limitation, and current semi-empirical water-use efficiency models fail to reproduce observed ET dynamics for these periods. Here, we use dry-down events occurring at eddy-covariance flux tower sites in the FLUXNET database as natural experiments to assess if and how decreasing soil-water availability modifies water-use efficiency on ecosystem scale. We demonstrate that an attenuating soil-water availability factor in junction with a previously discovered additive radiation term is necessary to accurately predict ET flux magnitudes and dry-down lengths of these water-limited periods. In an analysis of the attenuation, 20–33% of the observed decline in ET was due to the previously unconsidered soil-water availability effect. We conclude by noting the rates of ET decline differ significantly between FLUXNET sites with tall and short vegetation types and discuss the dependency of this rate on the variability of seasonal dryness. [ABSTRACT FROM AUTHOR]

The GRACE satellites provide signals of total terrestrial water storage (TWS) variations over large spatial domains at seasonal to inter-annual timescales. While the GRACE data have been extensively and successfully used to assess spatio-temporal changes in TWS, little effort has been made to quantify the relative contributions of snowpacks, soil moisture, and other components to the integrated TWS signal across northern latitudes, which is essential to gain a better insight into the underlying hydrological processes. Therefore, this study aims to assess which storage component dominates the spatio-temporal patterns of TWS variations in the humid regions of northern mid- to high latitudes. To do so, we constrained a rather parsimonious hydrological model with multiple state-of-the-art Earth observation products including GRACE TWS anomalies, estimates of snow water equivalent, evapotranspiration fluxes, and gridded runoff estimates. The optimized model demonstrates good agreement with observed hydrological spatio-temporal patterns and was used to assess the relative contributions of solid (snowpack) versus liquid (soil moisture, retained water) storage components to total TWS variations. In particular, we analysed whether the same storage component dominates TWS variations at seasonal and inter-annual temporal scales, and whether the dominating component is consistent across small to large spatial scales. Consistent with previous studies, we show that snow dynamics control seasonal TWS variations across all spatial scales in the northern mid- to high latitudes. In contrast, we find that inter-annual variations of TWS are dominated by liquid water storages at all spatial scales. The relative contribution of snow to inter-annual TWS variations, though, increases when the spatial domain over which the storages are averaged becomes larger. This is due to a stronger spatial coherence of snow dynamics that are mainly driven by temperature, as opposed to spatially more heterogeneous liquid water anomalies, that cancel out when averaged over a larger spatial domain. The findings first highlight the effectiveness of our model--data fusion approach that jointly interprets multiple Earth observation data streams with a simple model. Secondly, they reveal that the determinants of TWS variations in snow-affected northern latitudes are scale-dependent. In particular, they seem to be not merely driven by snow variability, but rather are determined by liquid water storages on interannual timescales. We conclude that inferred driving mechanisms of TWS cannot simply be transferred from one scale to another, which is of particular relevance for understanding the short- and long-term variability of water resources. [ABSTRACT FROM AUTHOR]

Understanding of terrestrial carbon and water cycles is currently hampered by an uncertainty in how to capture the large variety of plant responses to drought. In FLUXNET, the global network of CO2 and H2O flux observations, many sites do not uniformly report the ancillary variables needed to study drought response physiology. To this end, we outline two data-driven indicators based on diurnal energy, water, and carbon flux patterns derived directly from the eddy covariance data and based on theorized physiological responses to hydraulic and non-stomatal limitations. Hydraulic limitations (i.e. intra-plant limitations on water movement) are proxied using the relative diurnal centroid (C* ET), which measures the degree to which the flux of evapotranspiration (ET) is shifted toward the morning. Non-stomatal limitations (e.g. inhibitions of biochemical reactions, RuBisCO activity, and/or mesophyll conductance) are characterized by the Diurnal Water-Carbon Index (DWCI), which measures the degree of coupling between ET and gross primary productivity (GPP) within each day. As a proof of concept we show the response of the metrics at six European sites during the 2003 heat wave event, showing a varied response of morning shifts and decoupling. Globally, we found indications of hydraulic limitations in the form of significantly high frequencies of morning-shifted days in dry/Mediterranean climates and savanna/ evergreen plant functional types (PFTs), whereas high frequencies of decoupling were dominated by dry climates and grassland/savanna PFTs indicating a prevalence of nonstomatal limitations in these ecosystems. Overall, both the diurnal centroid and DWCI were associated with high net radiation and low latent energy typical of drought. Using three water use efficiency (WUE) models, we found the mean differences between expected and observed WUE to be -0.09 to 0.44 μmol mmol-1 and -0.29 to -0.40 μmol mmol-1 for decoupled and morning-shifted days, respectively, compared to mean differences -1.41 to -1.42 μmol mmol-1 in dry conditions, suggesting that morning shifts/hydraulic responses are associated with an increase in WUE, whereas decoupling/non-stomatal limitations are not. [ABSTRACT FROM AUTHOR]

The use of dynamic global vegetation models (DGVMs) to estimate CO2 emissions from land-use and land-cover change (LULCC) offers a new window to account for spatial and temporal details of emissions and for ecosystem processes affected by LULCC. One drawback of LULCC emissions from DGVMs, however, is lack of observation constraint. Here, we propose a new method of using satellite- and inventory-based biomass observations to constrain historical cumulative LULCC emissions (ELUCc) from an ensemble of nine DGVMs based on emerging relationships between simulated vegetation biomass and ELUCc. This method is applicable on the global and regional scale. The original DGVM estimates of ELUCc range from 94 to 273 PgC during 1901-2012. After constraining by current biomass observations, we derive a best estimate of 155±50 PgC (1σ Gaussian error). The constrained LULCC emissions are higher than prior DGVM values in tropical regions but significantly lower in North America. Our emergent constraint approach independently verifies the median model estimate by biomass observations, giving support to the use of this estimate in carbon budget assessments. The uncertainty in the constrained ELUCc is still relatively large because of the uncertainty in the biomass observations, and thus reduced uncertainty in addition to increased accuracy in biomass observations in the future will help improve the constraint. This constraint method can also be applied to evaluate the impact of land-based mitigation activities. [ABSTRACT FROM AUTHOR]

The GRACE satellites provide signals of total terrestrial water storage (TWS) variations over large spatial domains at seasonal to inter-annual time scales. While the GRACE data have been extensively and successfully used to assess spatio-temporal changes in TWS, little effort has been made to quantify the relative contributions of snow pack, soil moisture and other components to the integrated TWS signal across northern latitudes, which is essential to gain a better insight into the underlying hydrological processes. Therefore, this study aims to assess which storage component dominates the spatio-temporal patterns of TWS variations in the humid regions of northern mid-to-high latitudes. To do so, we constrained a rather parsimonious hydrological model with multiple state-of-the-art Earth observation products including GRACE TWS anomalies, estimates of snow water equivalent, evapotranspiration fluxes, and gridded runoff estimates. The optimized model demonstrates good agreement with observed hydrological patterns, and was used to assess the relative contributions of solid (snow pack) versus liquid (soil moisture, retained water) storage components to total TWS variations. In particular, we analysed whether the same storage component dominates TWS variations at seasonal and inter-annual temporal scales, and whether the dominating component is consistent across small to large spatial scales. Consistent with previous studies, we show that snow dynamics control seasonal TWS variations across spatial scales in the northern mid-to-high latitudes. In contrast, we find that inter-annual variations of TWS are dominated by liquid water storages, comprising mainly of soil moisture. However, as the spatial domain over which the storages are averaged becomes larger, the relative contribution of snow to inter-annual TWS variations increases. This is due to a stronger spatial coherence of snow anomalies as opposed to spatially more heterogeneous liquid water anomalies that cancel out over large spatial domains. The findings first highlight the effectiveness of our model-data fusion approach that jointly interprets multiple Earth observation data streams with a simple model. Secondly, they reveal that the determinants of TWS variations in snow-affected northern latitudes are scale dependent. We conclude that inferred driving mechanisms of TWS cannot simply be transferred from one scale to another, which is of particular relevance for understanding the short and long-term variability of water resources. [ABSTRACT FROM AUTHOR]

Understanding of terrestrial carbon and water cycles is currently hampered by an uncertainty in how to capture the large variety of plant responses to drought across climates, ecological strategies, and environments. In FLUXNET, the global network of CO2 and H2O flux observations, many sites do not uniformly report the ancillary variables needed to study drought response physiology such as soil moisture, sap flux, or species composition. In this sense, the use of diurnal energy, water, and carbon flux patterns to derive clues on ecosystem water limitation responses at a daily resolution could prove valuable, if nothing less than a benchmark to test current hypotheses. To this end, we propose two data-driven indicators derived directly from the eddy covariance data and based on theorized physiological responses to hydraulic and non-stomatal limitations. Hydraulic limitations (i.e. intra-plant limitations to water movement) are proxied using the relative diurnal centroid (C*ET), which measures the degree to which the flux of evapotranspiration (ET) is shifted toward the morning. Non-stomatal limitations (e.g. inhibitions of biochemical reactions, Rubisco activity, and/or mesophyll conductance) are characterized by the Diurnal Water:Carbon Index (DWCI), which measures the degree of coupling between ET and gross primary productivity (GPP) within each day. Globally, we found indications of hydraulic limitations in the form of significantly high frequencies of morning shifted days in dry/Mediterranean climates and savanna/evergreen plant functional types (PFT), whereas high frequencies of decoupling were dominated by dry climates and grassland/savanna PFTs indicating a prevalence of non-stomatal limitations in these ecosystems. Overall, both the diurnal centroid and DWCI were associated with high net radiation and low latent energy typical of drought. Using three water use efficiency (WUE) models, we found the mean differences between expected and observed WUE to be -0.14 to 0.56 µmol/mmol and -0.52 to -0.64 µmol/mmol for decoupled and morning shifted days respectively compared to mean differences -1.41 to -1.43 µmol/mmol in dry conditions. These results suggest that morning shifts/hydraulic responses are associated with an increase in WUE whereas decoupling/non-stomatal limitations are not. [ABSTRACT FROM AUTHOR]

Accurate model representation of land- atmosphere carbon fluxes is essential for climate projections. However, the exact responses of carbon cycle processes to climatic drivers often remain uncertain. Presently, knowledge derived from experiments, complemented by a steadily evolving body of mechanistic theory, provides the main basis for developing such models. The strongly increasing availability of measurements may facilitate new ways of identifying suitable model structures using machine learning. Here, we explore the potential of gene expression programming (GEP) to derive relevant model formulations based solely on the signals present in data by automatically applying various mathematical transformations to potential predictors and repeatedly evolving the resulting model structures. In contrast to most other machine learning regression techniques, the GEP approach generates "readable" models that allow for prediction and possibly for interpretation. Our study is based on two cases: artificially generated data and real observations. Simulations based on artificial data show that GEP is successful in identifying prescribed functions, with the prediction capacity of the models comparable to four state-of-the-art machine learning methods (random forests, support vector machines, artificial neural networks, and kernel ridge regressions). Based on real observations we explore the responses of the different components of terrestrial respiration at an oak forest in south-eastern England. We find that the GEP-retrieved models are often better in prediction than some established respiration models. Based on their structures, we find previously unconsidered exponential dependencies of respiration on seasonal ecosystem carbon assimilation and water dynamics. We noticed that the GEP models are only partly portable across respiration components, the identification of a "general" terrestrial respiration model possibly prevented by equifinality issues. Overall, GEP is a promising tool for uncovering new model structures for terrestrial ecology in the data-rich era, complementing more traditional modelling approaches. [ABSTRACT FROM AUTHOR]

Water-use efficiency (WUE) is a fundamental property for the coupling of carbon and water cycles in plants and ecosystems. Existing model formulations predicting this variable differ in the type of response of WUE to the atmospheric vapor pressure deficit of water (VPD). We tested a representative WUE model on the ecosystem scale at 110 eddy covariance sites of the FLUXNET initiative by predicting evapotranspiration (ET) based on gross primary productivity (GPP) and VPD. We found that introducing an intercept term in the formulation increases model performance considerably, indicating that an additional factor needs to be considered. We demonstrate that this intercept term varies seasonally and we subsequently associate it with radiation. Replacing the constant intercept term with a linear function of global radiation was found to further improve model predictions of ET. Our new semiempirical ecosystem WUE formulation indicates that, averaged over all sites, this radiation term accounts for up to half (39-47 %) of transpiration. These empirical findings challenge the current understanding of water-use efficiency on the ecosystem scale. [ABSTRACT FROM AUTHOR]

Water-use efficiency (WUE) is a fundamental property for the coupling of carbon and water cycles in plants and ecosystems. Existing model formulations predicting this variable differ in the type of response of WUE to the atmospheric vapor pressure deficit of water (VPD). We tested a representative WUE model on ecosystem scale at 110 eddy-covariance sites of the FLUXNET initiative by predicting evapotranspiration (ET) based on gross primary productivity (GPP) and VPD. We found that introducing an intercept term in the formulation increases model performance considerably, indicating that an additional factor needs to be considered. We demonstrate that this intercept term varies seasonally and we subsequently associate it with radiation. Replacing the constant intercept term with a linear function of global radiation was found to further improve model predictions of ET. Our new semi-empirical ecosystem WUE formulation indicates that, averaged over all sites, this radiation term accounts for up to half (40-49 %) of transpiration. These empirical findings challenge the current understanding of water-use efficiency on ecosystem-scale. [ABSTRACT FROM AUTHOR]

Accurate modelling of land-atmosphere carbon fluxes is essential for future climate projections. However, the exact responses of carbon cycle processes to climatic drivers often remain uncertain. Presently, knowledge derived from experiments complemented with a steadily evolving body of mechanistic theory provides the main basis for developing the respective models. The strongly increasing availability of measurements may complicate the traditional hypothesis driven path to developing mechanistic models, but it may facilitate new ways of identifying suitable model structures using machine learning as well. Here we explore the potential to derive model formulations automatically from data based on gene expression programming (GEP). GEP automatically (re)combines various mathematical operators to model formulations that are further evolved, eventually identifying the most suitable structures. In contrast to most other machine learning regression techniques, the GEP approach generates models that allow for prediction and possibly for interpretation. Our study is based on two cases: artificially generated data and real observations. Simulations based on artificial data show that GEP is successful in identifying prescribed functions with the prediction capacity of the models comparable to four state-of-the-art machine learning methods (Random Forests, Support Vector Machines, Artificial Neural Networks, and Kernel Ridge Regressions). The case of real observations explores different components of terrestrial respiration at an oak forest in south-east England. We find that GEP retrieved models are often better in prediction than established respiration models. Furthermore, the structure of the GEP models offers new insights to driver selection and interactions. We find previously unconsidered exponential dependencies of respiration on seasonal ecosystem carbon assimilation and water dynamics. However, we also noticed that the GEP models are only partly portable across respiration components; equifinality issues possibly preventing the identification of a "general" terrestrial respiration model. Overall, GEP is a promising tool to uncover new model structures for terrestrial ecology in the data rich era, complementing the traditional approach of model building. [ABSTRACT FROM AUTHOR]

Understanding the global carbon (C) cycle is of crucial importance to map current and future climate dynamics relative to global environmental change. A full characterization of C cycling requires detailed information on spatiotemporal patterns of surface-atmosphere fluxes. However, relevant C cycle observations are highly variable in their coverage and reporting standards. Especially problematic is the lack of integration of vertical oceanic, inland freshwaters and terrestrial carbon dioxide (CO2) exchange. Here we adopt a data-driven approach to synthesize a wide range of observation-based spatially explicit surface-atmosphere CO2 fluxes from 2001 and 2010, to identify the state of today's observational opportunities and data limitation. The considered fluxes include vertical net exchange of open oceans, continental shelves, estuaries, rivers, and lakes, as well as CO2 fluxes related to gross primary productivity, terrestrial ecosystem respiration, fire emissions, loss of tropical aboveground C, harvested wood and crops, as well as fossil fuel and cement emissions. Spatially explicit CO2 fluxes are obtained through geostatistical and/or remote sensing-based upscaling; minimizing biophysical or biogeochemical assumptions encoded in process-based models. We estimate a global bottom-up net C exchange (NCE) between the surface (land, ocean, and coastal areas) and the atmosphere. Uncertainties for NCE and its components are derived using resampling. In most continental regions our NCE estimates agree well with independent estimates from other sources. This holds for Europe (mean ±1 SD: 0.80 ± 0.16 PgC/yr, positive numbers are sources to the atmosphere), Russia (-0.02 ± 0.49 PgC/yr), East Asia (1.76 ± 0.38 PgC/yr), South Asia (0.25 ± 0.16 PgC/yr), and Australia (0.22 ± 0.47 PgC/yr). Our NCE estimates also suggest large C sink in tropical areas. The global NCE estimate is -6.07 ± 3.38 PgC/yr. This global bottom-up value is the opposite direction of what is expected from the atmospheric growth rate of CO2, and would require an offsetting surface C source of 4.27±0.10 PgC/yr. This mismatch highlights large knowledge and observational gaps in tropical areas, particularly in South America, Africa, and Southeast Asia, but also in North America. Our uncertainty assessment provides the basis for designing new observation campaigns. In particular, we lack seasonal monitoring of shelf, estuary and inland water-atmosphere C exchange. Also, extensive pCO2 measurements are missing in the Southern Ocean. Most importantly, tropical land C fluxes suffer from a lack of in-situ observations. The consistent derivation of data uncertainties could serve as prior knowledge in multi-criteria optimization such as the Carbon Cycle Data Assimilation System (CCDAS) without overstating data credibility. Furthermore, the spatially explicit flux estimates may be used as a starting point to assess the validity of countries' claims of reducing net C emissions in climate change negotiations. [ABSTRACT FROM AUTHOR]

We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the JSBACH land-surface scheme, which is part of the MPI-Earth System Model v1. The simulated phenology and net land carbon balance were constrained by globally distributed observations of the fraction of absorbed photosynthetically active radiation (FAPAR, using the TIP-FAPAR product) and atmospheric CO2 at a global set of monitoring stations for the years 2005 to 2009. When constrained by FAPAR observations alone, the system successfully, and computationally efficiently, improved simulated growing-season average FAPAR, as well as its seasonality in the northern extra-tropics. When constrained by atmospheric CO2 observations alone, global net and gross carbon fluxes were improved, despite a tendency of the system to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, thereby increasing the overall appropriateness of the simulated biosphere dynamics and underlying parameter values. Our study thus demonstrates the value of multiple-data-stream assimilation for the simulation of terrestrial biosphere dynamics. It further highlights the potential role of remote sensing data, here the TIP-FAPAR product, in stabilising the strongly underdetermined atmospheric inversion problem posed by atmospheric transport and CO2 observations alone. Notwithstanding these advances, the constraint of the observations on regional gross and net CO2 flux patterns on the MPI-CCDAS is limited through the coarse-scale parametrisation of the biosphere model. We expect improvement through a refined initialisation strategy and inclusion of further biosphere observations as constraints. [ABSTRACT FROM AUTHOR]