Abstract. 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. With the eddy covariance (EC) technique, net fluxes of carbon dioxide (CO2) and other trace gases as well as water and energy fluxes can be measured at the ecosystem level. These flux measurements are a main source for understanding biosphere–atmosphere interactions and feedbacks through cross-site analysis, model–data integration, and upscaling. The raw fluxes measured with the EC technique require extensive and laborious data processing. While there are standard tools1 available in an open-source environment for processing high-frequency (10 or 20 Hz) data into half-hourly quality-checked fluxes, there is a need for more usable and extensible tools for the subsequent post-processing steps. We tackled this need by developing the REddyProc package in the cross-platform language R that provides standard CO2-focused post-processing routines for reading (half-)hourly data from different formats, estimating the u* threshold, as well as gap-filling, flux-partitioning, and visualizing the results. In addition to basic processing, the functions are extensible and allow easier integration in extended analysis than current tools. New features include cross-year processing and a better treatment of uncertainties. A comparison of REddyProc routines with other state-of-the-art tools resulted in no significant differences in monthly and annual fluxes across sites. Lower uncertainty estimates of both u* and resulting gap-filled fluxes by 50 % with the presented tool were achieved by an improved treatment of seasons during the bootstrap analysis. Higher estimates of uncertainty in daytime partitioning (about twice as high) resulted from a better accounting for the uncertainty in estimates of temperature sensitivity of respiration. The provided routines can be easily installed, configured, and used. Hence, the eddy covariance community will benefit from the REddyProc package, allowing easier integration of standard post-processing with extended analysis. 1http://fluxnet.fluxdata.org/2017/10/10/toolbox-a-rolling-list-of-softwarepackages-for-flux-related-data-processing/, last access: 17 August 2018

The flux of carbon dioxide from the soil to the atmosphere (soil respiration) is one of the major fluxes in the global carbon cycle. At present, the accumulated field observation data cover a wide range of geographical locations and climate conditions. However, there are still large uncertainties in the magnitude and spatiotemporal variation of global soil respiration. Using a global soil respiration dataset, we developed a climate-driven model of soil respiration by modifying and updating Raich's model, and the global spatiotemporal distribution of soil respiration was examined using this model. The model was applied at a spatial resolution of 0.5° and a monthly time step. Soil respiration was divided into the heterotrophic and autotrophic components of respiration using an empirical model. The estimated mean annual global soil respiration was 91 Pg C yr-1 (between 1965 and 2012; Monte Carlo 95% confidence interval: 87–95 Pg C yr-1) and increased at the rate of 0.09 Pg C yr-2. The contribution of soil respiration from boreal regions to the total increase in global soil respiration was on the same order of magnitude as that of tropical and temperate regions, despite a lower absolute magnitude of soil respiration in boreal regions. The estimated annual global heterotrophic respiration and global autotrophic respiration were 51 and 40 Pg C yr-1, respectively. The global soil respiration responded to the increase in air temperature at the rate of 3.3 Pg C yr-1 °C−1, and Q10 = 1.4. Our study scaled up observed soil respiration values from field measurements to estimate global soil respiration and provide a data-oriented estimate of global soil respiration. Our results, including the modeled spatiotemporal distribution of global soil respiration, are based on a semi-empirical model parameterized with over one thousand data points. We expect that these spatiotemporal estimates will provide a benchmark for future studies and also help to constrain process-oriented models

Plant phenological development is orchestrated through subtle changes in photoperiod, temperature, soil moisture and nutrient availability. Presently, the exact timing of plant development stages and their response to climate and management practices 5 are crudely represented in land surface models. As visual observations of phenology are laborious, there is a need to supplement long-term observations with automated techniques such as those provided by digital repeat photography at high temporal and spatial resolution. We present the first synthesis from a growing observational network of digital cameras installed on towers across Europe above deciduous and evergreen 10 forests, grasslands and croplands, where vegetation and atmosphere CO2 fluxes are measured continuously. Using colour indices from digital images and using piecewise regression analysis of time-series, we explored whether key changes in canopy phenology could be detected automatically across different land use types in the network. The piecewise regression approach could capture the start and end of the growing 15 season, in addition to identifying striking changes in colour signals caused by flowering and management practices such as mowing. Exploring the dates of green up and senescence of deciduous forests extracted by the piecewise regression approach against dates estimated from visual observations we found that these phenological events could be detected adequately (RMSE< 8 and 11 days for leaf out and leaf fall 20 respectively). We also investigated whether the seasonal patterns of red, green and blue colour fractions derived from digital images could be modelled mechanistically using the PROSAIL model parameterised with information of seasonal changes in canopy leaf area and leaf chlorophyll and carotenoid concentrations. From a model sensitivity analysis we found that variations in colour fractions, and in particular the late spring 25 “green hump” observed repeatedly in deciduous broadleaf canopies across the network, are essentially dominated by changes in the respective pigment concentrations. Using the model we were able to explain why this spring maximum in green signal is often observed out of phase with the maximum period of canopy photosynthesis in ecosystems across Europe. Coupling such quasi-continuous digital records of canopy colours with co-located CO2 flux measurements will improve our understanding of how changes in growing season length are likely to shape the capacity of European ecosystems to sequester CO2 in the future

Abstract. Existing dynamic global vegetation models (DGVMs) have a limited ability in reproducing phenology and decadal dynamics of vegetation greenness as observed by satellites. These limitations in reproducing observations reflect a poor understanding and description of the environmental controls on phenology, which strongly influence the ability to simulate longer-term vegetation dynamics, e.g. carbon allocation. Combining DGVMs with observational data sets can potentially help to revise current modelling approaches and thus enhance the understanding of processes that control seasonal to long-term vegetation greenness dynamics. Here we implemented a new phenology model within the LPJmL (Lund Potsdam Jena managed lands) DGVM and integrated several observational data sets to improve the ability of the model in reproducing satellite-derived time series of vegetation greenness. Specifically, we optimized LPJmL parameters against observational time series of the fraction of absorbed photosynthetic active radiation (FAPAR), albedo and gross primary production to identify the main environmental controls for seasonal vegetation greenness dynamics. We demonstrated that LPJmL with new phenology and optimized parameters better reproduces seasonality, inter-annual variability and trends of vegetation greenness. Our results indicate that soil water availability is an important control on vegetation phenology not only in water-limited biomes but also in boreal forests and the Arctic tundra. Whereas water availability controls phenology in water-limited ecosystems during the entire growing season, water availability co-modulates jointly with temperature the beginning of the growing season in boreal and Arctic regions. Additionally, water availability contributes to better explain decadal greening trends in the Sahel and browning trends in boreal forests. These results emphasize the importance of considering water availability in a new generation of phenology modules in DGVMs in order to correctly reproduce observed seasonal-to-decadal dynamics of vegetation greenness.