Search

Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20∘ S to 20∘ N) the spring onset of elevated CH4 emissions starts 3 d earlier, and the CH4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

While wetlands are the largest natural source of methane (CH4 ) to the atmosphere, they represent a large source of uncertainty in the global CH4 budget due to the complex biogeochemical controls on CH4 dynamics. Here we present, to our knowledge, the first multi-site synthesis of how predictors of CH4 fluxes (FCH4) in freshwater wetlands vary across wetland types at diel, multiday (synoptic), and seasonal time scales. We used several statistical approaches (correlation analysis, generalized additive modeling, mutual information, and random forests) in a wavelet-based multi-resolution framework to assess the importance of environmental predictors, nonlinearities and lags on FCH4 across 23 eddy covariance sites. Seasonally, soil and air temperature were dominant predictors of FCH4 at sites with smaller seasonal variation in water table depth (WTD). In contrast, WTD was the dominant predictor for wetlands with smaller variations in temperature (e.g., seasonal tropical/subtropical wetlands). Changes in seasonal FCH4 lagged fluctuations in WTD by ~17±11days, and lagged air and soil temperature by median values of 8±16 and 5±15days, respectively. Temperature and WTD were also dominant predictors at the multiday scale. Atmospheric pressure (PA) was another important multiday scale predictor for peat-dominated sites, with drops in PA coinciding with synchronous releases of CH4 . At the diel scale, synchronous relationships with latent heat flux and vapor pressure deficit suggest that physical processes controlling evaporation and boundary layer mixing exert similar controls on CH4 volatilization, and suggest the influence of pressurized ventilation in aerenchymatous vegetation. In addition, 1- to 4-h lagged relationships with ecosystem photosynthesis indicate recent carbon substrates, such as root exudates, may also control FCH4. By addressing issues of scale, asynchrony, and nonlinearity, this work improves understanding of the predictors and timing of wetland FCH4 that can inform future studies and models, and help constrain wetland CH4 emissions.

Ground reference data are an essential prerequisite for supervised crop mapping. The lack of a low-cost and efficient ground referencing method results in pervasively limited reference data and hinders crop classification. In this study, we apply a convolutional neural network (CNN) model to explore the efficacy of automatic ground truthing via Google Street View (GSV) images in two distinct farming regions: Illinois and the Central Valley in California. We demonstrate the feasibility and reliability of our new ground referencing technique by performing pixel-based crop mapping at the state level using the cloud-based Google Earth Engine platform. The mapping results are evaluated using the United States Department of Agriculture (USDA) crop data layer (CDL) products. From ~ 130,000 GSV images, the CNN model identified ~ 9,400 target crop images. These images are well classified into crop types, including alfalfa, almond, corn, cotton, grape, rice, soybean, and pistachio. The overall GSV image classification accuracy is 92% for the Central Valley and 97% for Illinois. Subsequently, we shifted the image geographical coordinates 2–3 times in a certain direction to produce 31,829 crop reference points: 17,358 in Illinois, and 14,471 in the Central Valley. Evaluation of the mapping results with CDL products revealed satisfactory coherence. GSV-derived mapping results capture the general pattern of crop type distributions for 2011–2019. The overall agreement between CDL products and our mapping results is indicated by R2 values of 0.44–0.99 for the Central Valley and 0.81–0.98 for Illinois. To show the applicational value of the proposed method in other countries, we further mapped rice paddy (2014–2018) in South Korea which yielded fairly well outcomes (R2 = 0.91). These results indicate that GSV images used with a deep learning model offer an efficient and cost-effective alternative method for ground referencing, in many regions of the world.

Quantifying global terrestrial photosynthesis is essential to understanding the global carbon cycle and the climate system. Remote sensing has played a pivotal role in advancing our understanding of photosynthesis from leaf to global scale; however, substantial uncertainties still exist. In this review, we provide a historical overview of theory, modeling, and observations of photosynthesis across space and time for decadal intervals beginning in the 1950s. Then we identify the key uncertainties in global photosynthesis estimates, including evaluating light intercepted by canopies, biophysical forcings, the structure of light use efficiency models and their parameters, like photosynthetic capacity, and relationships between sun-induced chlorophyll fluorescence and canopy photosynthesis. Finally, we review new opportunities with big data and recently launched or planned satellite missions.

International audience; Global vegetation and land-surface models embody interdisciplinary scientific understanding of the behaviour of plants and ecosystems, and are indispensable to project the impacts of environmental change on vegetation and the interactions between vegetation and climate. However, systematic errors and persistently large differences among carbon and water cycle projections by different models highlight the limitations of current process formulations. In this review, focusing on core plant functions in the terrestrial carbon and water cycles, we show how unifying hypotheses derived from eco-evolutionary optimality (EEO) principles can provide novel, parameter-sparse representations of plant and vegetation processes. We present case studiesthat demonstrate how EEO generate parsimonious representations of core, leaf-level processes that are individually testable and supported by evidence. EEO approaches to photosynthesis and primary production, dark respiration, and stomatal behaviour are ripe for implementation in global models. EEO approaches to other important traits, including the leaf economics spectrum and applications of EEO at the community level are active research areas. Independently tested modules emerging from EEO studies could profitably be integrated into modelling frameworks that account for the multiple time scales on which plants and plant communities adjust toenvironmental change.

Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions, including their seasonality, due to quasi-continuous and high temporal resolution of flux measurements, coincident measurements of carbon, water, and energy fluxes, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we 1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4- community-product/). FLUXNET-CH4 includes half-hourly and daily gap-filled and non gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we 2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally, because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands and because freshwater wetlands are a substantial source of total atmospheric CH4 emissions; and 3) provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions, but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20° S to 20° N) the spring onset of elevated CH4 emissions starts three days earlier, and the CH4 emission season lasts 4 days longer, for each degree C increase in mean annual air temperature. On average, the onset of increasing CH4 emissions lags soil warming by one month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling, and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). The FLUXNET-CH4 dataset provides an open-access resource for CH4 flux synthesis, has a range of applications, and is unique in that it includes coupled measurements of important CH4 drivers such as GPP and temperature. Although FLUXNET-CH4 could certainly be improved by adding more sites in tropical ecosystems and by increasing the number of site-years at existing sites, it is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4408468. Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/, and a complete list of the 79 individual site data DOIs is provided in Table 2 in the Data Availability section of this document., Earth System Science Data Discussions, ISSN:1866-3591

Leaf-level gas exchange data support the mechanistic understanding of plant fluxes of carbon and water. These fluxes inform our understanding of ecosystem function, are an important constraint on parameterization of terrestrial biosphere models, are necessary to understand the response of plants to global environmental change, and are integral to efforts to improve crop production. Collection of these data using gas analyzers can be both technically challenging and time consuming, and individual studies generally focus on a small range of species, restricted time periods, or limited geographic regions. The high value of these data is exemplified by the many publications that reuse and synthesize gas exchange data, however the lack of metadata and data reporting conventions make full and efficient use of these data difficult. Here we propose a reporting format for leaf-level gas exchange data and metadata to provide guidance to data contributors on how to store data in repositories to maximize their discoverability, facilitate their efficient reuse, and add value to individual datasets. For data users, the reporting format will better allow data repositories to optimize data search and extraction, and more readily integrate similar data into harmonized synthesis products. The reporting format specifies data table variable naming and unit conventions, as well as metadata characterizing experimental conditions and protocols. For common data types that were the focus of this initial version of the reporting format, i.e., survey measurements, dark respiration, carbon dioxide and light response curves, and parameters derived from those measurements, we took a further step of defining required additional data and metadata that would maximize the potential reuse of those data types. To aid data contributors and the development of data ingest tools by data repositories we provided a translation table comparing the outputs of common gas exchange instruments. Extensive consultation with data collectors, data users, instrument manufacturers, and data scientists was undertaken in order to ensure that the reporting format met community needs. The reporting format presented here is intended to form a foundation for future development that will incorporate additional data types and variables as gas exchange systems and measurement approaches advance in the future. The reporting format is published in the U.S. Department of Energy's ESS-DIVE data repository, with documentation and future development efforts being maintained in a version control system. publishedVersion

Continuous monitoring of vegetation indices (VIs) the fraction of absorbed photosynthetically active radiation (fPAR) and leaf area index (LAI) through satellite remote sensing has advanced our understanding of biosphere–atmosphere interactions. Substantial efforts have been put into monitoring individual variables in the field, but options to concurrently monitor VIs, fPAR, and LAI in-situ have been lacking. In this paper, we present the Smart Surface Sensing System (4S), which automatically collects, transfers and processes VIs, fPAR and LAI data streams. The 4S consists of a microcomputer, controller and camera, a multi-spectral spectrometer built in with a light-emitting diode (LED) and an internet connection. Lab testing and field observations in a rice paddy site that experiences wet summer monsoon seasons confirmed the linear response of 4S to light intensities in the blue, green, red and near-infrared spectral channels, with wide ranging temperatures and humidity having only a minor impact on 4S throughout the growing season. Applied over an entire rice growing season (day of year [DOY] 120 - 248), VIs and fPAR from 4S were linearly related to corresponding VIs from a reference spectrometer (R2 = 0.98; NDVI, R2 = 0.96; EVI) and the LAI-2200 instrument (R2 = 0.76), respectively. Integration of gap fraction-based LAI from LED sensors and a green index from the micro-camera allowed tracking of the seasonality of green LAI. The continuous and diverse nature of 4S observations highlights its potential for evaluating satellite remote sensing products. We believe that 4S will be useful for the expansion of ecological sensing networks across multiple spatial and temporal scales.

The ECOsystem Space Thermal Radiometer Experiment on the Space Station (ECOSTRESS) was launched to the International Space Station on June 29, 2018 by the National Aeronautics and Space Administration (NASA). The primary scientific focus of ECOSTRESS is on evapotranspiration (ET), which occurs as Level 3 (L3) latent heat flux (LE) data products. This data is generated from the Level 2 Earth Surface Temperature and Emissivity Product (L2_LSTE), along with the atmospheric and ancillary surface data. Here, we provide the first validation (Stage 1, preliminary) of the global product ECOSTRESS clear - sky ET (L3_ET_PT - JPL, Version 6.0) against LE measurements at 82 eddy covariance sites around the world. Overall, the ECOSTRESS ET product performs well against site measurements (instant clear sky / flyover time: r 2 = 0.88; overall bias = 8%; normalized root mean square error, RMSE = 6%) . ET uncertainty was generally constant across climate zones, biome types, and times of day (ECOSTRESS samples the diurnal cycle), although temperate sites are over-represented. ECOSTRESS 70 m high spatial resolution improved correlations by 85% and RMSE by 62%, relative to 1 km pixels. This document serves as a reference for the accuracy of ECOSTRESS L3 ET and the Stage 1 validation status for subsequent science that continues to use this data. Fil: Fisher, Joshua B. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Lee, Brian. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Purdy, Adam J. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Halverson, Gregory H. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Dohlen, Matthew B. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Cawse - Nicholson, Kerry. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Wang, Audrey. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Anderson, Ray G. U.S. Department of Agriculture; Estados Unidos Fil: Aragon, Bruno. King Abdullah University of Science and Technology. Division of Biological and Environmental Science and Engineering. Hydrology, Agriculture and Land Observation Group (HALO); Arabia Saudita Fil: Arain, M. Alfat. McMaster University. School of Geography and Earth Sciences; Canadá Fil: Baldocchi, Dennis D. University of California, Berkeley. Department of Environmental Science, Policy and Management; Estados Unidos Fil: Baker, John M. U.S. Department of Agriculture; Estados Unidos Fil: Barral, Helene. University Montpellier. HydroSciences Montpellier (HSM); Francia Fil: Bernacchi, Carl J. University of Illinois Urbana‐Champaign. Institute for Sustainability, Energy, and Environment and Carl R. Woese Institute for Genomic Biology. Center for Advanced Bioenergy and Bioproducts Innovation; Estados Unidos. U.S. Department of Agriculture, Global Change and Photosynthesis Research Unit; Estados Unidos Fil: Bernhofer, Christian. Technische Universität Dresden, Institute of Hydrology and Meteorology; Alemania Fil: Biraud, Sebastien C. Lawrence Berkeley National Laboratory; Estados Unidos Fil: Bohrer, Gil. Ohio State University. Department of Civil, Environmental and Geodetic Engineering; Estados Unidos Fil: Brunsell, Nathaniel. University of Kansas. Department of Geography and Atmospheric Science; Estados Unidos Fil: Cappelaere, Bernard. University Montpellier. HydroSciences Montpellier (HSM); Francia Fil: Castro - Contreras, Saulo. University of Alberta. Earth and Atmospheric Sciences Department. Centre for Earth Observation Sciences (CEOS); Canadá Fil: Chum, Junghwa. National Institute of Forest Science; Corea del Sur Fil: Conrad, Bryan. University of Kansas. Department of Geography and Atmospheric Science; Estados Unidos Fil: Cremonese, Edoardo. Environmental Protection Agency of Aosta Valley; Italia Fil: Demarty, Jerome. University Montpellier. HydroSciences Montpellier (HSM); Francia Fil: Desai, Ankur R. University of Wisconsin‐Madison. Department of Atmospheric and Oceanic Sciences; Estados Unidos Fil: De Ligne, Anne. University of Liege. TERRA Teaching and Research Centre. Gembloux Agro‐Bio Tech; Bélgica Fil: Foltynova, Lenka. Czech Academy of Sciences. Global Change Research Institute ; República Checa Fil: Goulden, Michael L. University of California Irvine. Department of Earth System Science; Estados Unidos Fil: Griffis, Timothy J. University of Minnesota. Department of Soil, Water, and Climate; Estados Unidos Fil: Grunwald, Thomas. Technische Universität Dresden, Institute of Hydrology and Meteorology; Alemania Fil: Johnson, Mark S. University of British Columbia. Department of Earth, Ocean, and Atmospheric Sciences; Canadá Fil: Kang, Minseok. National Center for AgroMeteorology; Corea del Sur Fil: Kelbe, Dave. Xerra Earth Observation Institute; Nueva Zelanda Fil: Kowalska, Natalia. Czech Academy of Sciences. Global Change Research Institute ; República Checa Fil: Jong - Hwan, Lim. National Institute of Forest Science; Corea del Sur Fil: Mainassara, Ibrahim. University Montpellier. HydroSciences Montpellier (HSM); Francia Fil: McCabe, Matthew F. King Abdullah University of Science and Technology. Division of Biological and Environmental Science and Engineering. Hydrology, Agriculture and Land Observation Group (HALO); Arabia Saudita Fil: Missik, Justine E.C. Washington State University. Department of Civil and Environmental Engineering. Laboratory for Atmospheric Research; Estados Unidos Fil: Mohanty, Binayak P. Texas A&M University. Texas Water Observatory; Estados Unidos Fil: Moore, Caitlin E. University of Illinois Urbana‐Champaign. Institute for Sustainability, Energy, and Environment and Carl R. Woese Institute for Genomic Biology. Center for Advanced Bioenergy and Bioproducts Innovation; Estados Unidos Fil: Morillas, Laura. University of British Columbia. Department of Earth, Ocean, and Atmospheric Sciences; Canadá Fil: Morrison, Ross. Centre for Ecology and Hydrology; Reino Unido Fil: Munger, J. Willians. Harvard University. School of Engineering and Applied Sciences; Estados Unidos Fil: Posse, Gabriela. Instituto Nacional de Tecnología Agropecuaria (INTA). Instituto de Clima y Agua; Argentina Fil: Richardson, Andrew D. Northern Arizona University Center for Ecosystem Science and Society; Estados Unidos Fil: Russell, Eric. S. Washington State University. Department of Civil and Environmental Engineering. Laboratory for Atmospheric Research; Estados Unidos Fil: Ryu, Youngryel. Seoul National University . Department of Landscape Architecture and Rural Systems Engineering; Corea del Sur Fil: Sanchez - Azofeifa, Arturo. University of Alberta. Earth and Atmospheric Sciences Department. Centre for Earth Observation Sciences (CEOS); Canadá Fil: Schmidt, Marius. IInstitute of Bio‐ and Geosciences: Agrosphere (IBG‐3) Forschungszentrum Jülich IBG‐3; Alemania Fil: Schwartz, Efrat. Weizmann Institute of Science. Department of Earth and Planetary Sciences, ; Israel Fil: Sharp, Iain. University of Alberta. Earth and Atmospheric Sciences Department. Centre for Earth Observation Sciences (CEOS); Canadá Fil: Sigut, Ladislav. Czech Academy of Sciences. Global Change Research Institute ; República Checa Fil: Tang, Yao. Georgia Institute of Technology. School of Civil and Environmental Engineering; Estados Unidos Fil: Hulley, Glynn. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos Fil: Anderson, Martha. U.S. Department of Agriculture‐Agricultural Research Service; Estados Unidos Fil: Hain, Christopher. NASA Marshall Space Flight Center; Estados Unidos Fil: French, Andrew. U.S. Department of Agriculture; Estados Unidos Fil: Wood, Eric. Princeton University. Department of Civil and Environmental Engineering; Estados Unidos Fil: Hook, Simón. California Institute of Technology. Jet Propulsion Laboratory; Estados Unidos

Changes in vegetation structure are expected to influence the redistribution of heat and moisture; however, how variations in the leaf area index (LAI) affect this global energy partitioning is not yet quantified. Here, we estimate that a unit change in LAI leads to 3.66 ± 0.45 and −3.26 ± 0.41 W m−2 in latent (LE) and sensible (H) fluxes, respectively, over the 1982–2016 period. Analysis of an ensemble of data-driven products shows that these sensitivities increase by about 20% over the observational period, prominently in regions with a limited water supply, probably because of an increased transpiration/evaporation ratio. Global greening has caused a decrease in the Bowen ratio (B = H/LE) of −0.010 ± 0.002 per decade, which is attributable to the increased evaporative surface. Such a direct LAI effect on energy fluxes is largely modulated by plant functional types (PFTs) and background climate conditions. Land surface models (LSMs) misrepresent this vegetation control, possibly due to underestimation of the biophysical responses to changes in the water availability and poor representation of LAI dynamics. Changes in the leaf area index alter the distribution of heat and moisture. The change in energy partitioning related to leaf area, increasing latent and decreasing sensible fluxes over the observational period 1982–2016, is moderated by plant functional type and background climate.

The maximum rate of carboxylation (Vcmax ) is an essential leaf trait determining the photosynthetic capacity of plants. Existing approaches for estimating Vcmax at large scale mainly rely on empirical relationships with proxies such as leaf nitrogen/chlorophyll content or hyperspectral reflectance, or on complicated inverse models from gross primary production or solar-induced fluorescence. A novel mechanistic approach based on the assumption that plants optimize resource investment coordinating with environment and growth has been shown to accurately predict C3 plant Vcmax based on mean growing season environmental conditions. However, the ability of optimality theory to explain seasonal variation in Vcmax has not been fully investigated. Here, we adapt an optimality-based model to simulate daily Vcmax,25C (Vcmax at a standardized temperature of 25°C) by incorporating the effects of antecedent environment, which affects current plant functioning, and dynamic light absorption, which coordinates with plant functioning. We then use seasonal Vcmax,25C field measurements from 10 sites across diverse ecosystems to evaluate model performance. Overall, the model explains about 83% of the seasonal variation in C3 plant Vcmax,25C across the 10 sites, with a medium root mean square error of 12.3 μmol m-2 s-1 , which suggests that seasonal changes in Vcmax,25C are consistent with optimal plant function. We show that failing to account for acclimation to antecedent environment or coordination with dynamic light absorption dramatically decreases estimation accuracy. Our results show that optimality-based approach can accurately reproduce seasonal variation in canopy photosynthetic potential, and suggest that incorporating such theory into next-generation trait-based terrestrial biosphere models would improve predictions of global photosynthesis.

Sun-induced chlorophyll fluorescence (SiF) is increasingly used as a proxy for vegetation canopy photosynthesis. While ground-based, airborne, and satellite observations have demonstrated a strong linear relationship between SiF and gross primary production (GPP) at seasonal scales, their relationships at high temporal resolution across diurnal to seasonal scales remain unclear. In this study, far-red canopy SiF, GPP, and absorbed photosynthetically active radiation (APAR) were continuously monitored using automated spectral systems and an eddy flux tower over an entire growing season in a rice paddy. At half-hourly resolution, strong linear relationships between SiF and GPP (R2 = 0.76) and APAR and GPP (R2 = 0.76) for the whole growing season were observed. We found that relative humidity, diffuse PAR fraction, and growth stage influenced the relationships between SiF and GPP, and APAR and GPP, and incorporating those factors into multiple regression analysis led to improvements up to R2 = 0.83 and R2 = 0.88, respectively. Relationships between LUEp (=GPP/APAR) and LUEf (=SiF/APAR) were inconsistent at half-hourly and weak at daily resolutions (R2 = 0.24). Both at diurnal and seasonal time scales with half-hourly resolution, we found considerably stronger linear relationships between SiF and APAR than between either SiF and GPP or APAR and GPP. Overall, our results indicate that for subdiurnal temporal resolution, canopy SiF in the rice paddy is above all a very good proxy for APAR at diurnal and seasonal time scales and that therefore SiF-based GPP estimation needs to take into account relevant environmental information to model LUEp. These findings can help develop mechanistic links between canopy SiF and GPP across multiple temporal scales.

Conventional process-based crop simulation models and agro-land surface models require numerous forcing variables and input parameters. The regional application of these crop simulation models is complicated by factors concerning input data requirements and parameter uncertainty. In addition, the empirical remotely sensed regional scale crop yield estimation method does not enable growth process modeling. In this study, we developed a process-based rice yield estimation model by integrating an assimilate allocation module into the satellite remote sensing-derived and biophysical process-based Breathing Earth System Simulator (BESS). Normalized accumulated gross primary productivity ( G P P n o r m - a c c u ) was used as a scaler for growth development, and the relationships between G P P n o r m - a c c u and dry matter partitioning coefficients were determined from the eddy covariance and biometric measurements at the Cheorwon Rice paddy KoFlux site. Over 95% of the variation in the dry matter allocation coefficients of rice grain could be explained by G P P n o r m - a c c u . The dynamics of dry matter distribution among different rice components were simulated, and the annual grain yields were estimated. BESS-Rice simulated GPP and dry matter partitioning dynamics, and rice yields were evaluated against in-situ measurements at three paddy rice sites registered in KoFlux. The results showed that BESS-Rice performed well in terms of rice productivity estimation, with average root mean square error (RMSE) value of 2.2 g C m−2 d−1 (29.5%) and bias of –0.5 g C m−2 d−1 (–7.1%) for daily GPP, and an average RMSE value of 534.8 kg ha−1 (7.7%) and bias of 242.1 kg ha−1 (3.5%) for the annual yield, respectively. BESS-Rice is much simpler than conventional crop models and this helps to reduce the uncertainty related to the forcing variables and input parameters and can result in improved regional yield estimation. The process-based mechanism of BESS-Rice also enables an agronomic diagnosis to be made and the potential impacts of climate change on rice productivity to be investigated.

Time series of vegetation indices (e.g. normalized difference vegetation index [NDVI]) and color indices (e.g. green chromatic coordinate [GCC]) based on radiometric measurements are now available at different spatial and temporal scales ranging from weekly satellite observations to sub-hourly in situ measurements by means of near-surface remote sensing (e.g. spectral sensors or digital cameras). In situ measurements are essential for providing validation data for satellite-derived vegetation indices. In this study we used a recently developed method to calculate NDVI from near-infrared (NIR) enabled digital cameras (NDVIC) at 17 sites (for a total of 74 year-sites) encompassing six plant functional types (PFT) from the PhenoCam network. The seasonality of NDVIC was comparable to both NDVI measured by ground spectral sensors and by the moderate resolution imaging spectroradiometer (MODIS). We calculated site- and PFT-specific scaling factors to correct NDVIC values and recommend the use of site-specific NDVI from MODIS in order to scale NDVIC. We also compared GCC extracted from red-green-blue images to NDVIC and found PFT-dependent systematic differences in their seasonalities. During senescence, NDVIC lags behind GCC in deciduous broad-leaf forests and grasslands, suggesting that GCC is more sensitive to changes in leaf color and NDVIC is more sensitive to changes in leaf area. In evergreen forests, NDVIC peaks later than GCC in spring, probably tracking the processes of shoot elongation and new needle formation. Both GCC and NDVIC can be used as validation tools for the MODIS Land Cover Dynamics Product (MCD12Q2) for deciduous broad-leaf spring phenology, whereas NDVIC is more comparable than GCC with autumn phenology derived from MODIS. For evergreen forests, we found a poor relationship between MCD12Q2 and camera-derived phenology, highlighting the need for more work to better characterize the seasonality of both canopy structure and leaf biochemistry in those ecosystems. Our results demonstrate that NDVIC is in excellent agreement with NDVI obtained from spectral measurements, and that NDVIC and GCC can complement each other in describing ecosystem phenology. Additionally, NDVIC allows the detection of structural changes in the canopy that cannot be detected by visible-wavelength imagery.

International audience; Disentangling the individual contributions from vegetation and soil in measured canopy reflectance is a grand challenge to the remote sensing and ecophysiology communities. Since Solar Induced chlorophyll Fluorescence (SIF) is uniquely emitted from vegetation, it can be used to evaluate how well reflectance-based vegetation indices (VIs) can separate the vegetation and soil components. Due to the residual soil background contributions, Near-infrared (NIR) reflectance of vegetation (NIRv) and Difference Vegetation index (DVI) present offsets when compared to SIF (i.e., the value of NIRv or DVI is non-zero when SIF is zero). In this study, we proposed a simple framework for estimating the true NIR reflectance of vegetation from Hyperspectral measurements (NIRvH) with minimal soil impacts. NIRvH takes advantage of the spectral shape variations in the red-edge region to minimize the soil effects. We evaluated the capability of NIRvH, NIRv and DVI in isolating the true NIR reflectance of vegetation using the data from both the model-based simulations and Hyperspectral Plant imaging spectrometer (HyPlant) measurements. Benchmarked by simultaneously measured SIF, NIRvH has the smallest offset (0–0.037), as compared to an intermediate offset of 0.047–0.062 from NIRv, and the largest offset of 0.089–0.112 from DVI. The magnitude of the offset can vary with different soil reflectance spectra across spatio-temporal scales, which may lead to bias in the downstream NIRv-based photosynthesis estimates. NIRvH and SIF measurements from the same sensor platform avoided complications due to different geometry, footprint and time of observation across sensors when studying the radiative transfer of reflected photons and SIF. In addition, NIRvH was primarily determined by canopy structure rather than chlorophyll content and soil brightness. Our work showcases that NIRvH is promising for retrieving canopy structure parameters such as leaf area index and leaf inclination angle, and for estimating fluorescence yield with current and forthcoming hyperspectral satellite measurements.

Soil respiration (Rs) determines land surface carbon balance; however, there have been few studies that measured Rs in heterogeneous urban landscapes. Here, we investigated the spatial and temporal variations in Rs in six land cover types (mixed forest, deciduous broadleaf forest, evergreen needleleaf forest, lawn, wetland, and bare land) in Seoul Forest Park, Republic of Korea, between March 2013 and September 2014, which included a wet (2013) and an extremely dry (2014) summer. Spatially, there was a three-fold difference (0.48–1.45 kgC m−2) in annual Rs among the six land cover types. The soil organic carbon stock at a depth of 0.1 m explained 72% of the spatial variation in the annual Rs across the land cover types. During the entire study period, the soil temperature explained 82–97% of the temporal variation in Rs among different land cover types. Comparing the two summers, the 2014 drought only resulted in a decrease in Rs in the lawn plots (25%), which was driven by a reduction in the leaf area index and the fine root density. The temperature sensitivity of Rs in 2014 (dry summer) compared to 2013 (wet summer) was significantly lower in mixed forest, deciduous broadleaf forest, and lawn, and did not change in evergreen needleleaf forest, wetland, or bare land. The differences in Rs in these drought responses highlight the importance of the careful selection of land cover type during park planning to better manage carbon cycles.

The savanna complex is a highly diverse global biome that occurs within the seasonally dry tropical to sub-tropical equatorial latitudes and are structurally and functionally distinct from grasslands and forests. Savannas are open-canopy environments that encompass a broad demographic continuum, often characterised by a changing dominance between C3-tree and C4-grass vegetation, where frequent environmental disturbances such as fire modulates the balance between ephemeral and perennial life forms. Climate change is projected to result in significant changes to the savanna floristic structure, with increases to woody biomass expected through CO2 fertilisation in mesic savannas and increased tree mortality expected through increased rainfall interannual variability in xeric savannas. The complex interaction between vegetation and climate that occurs in savannas has traditionally challenged terrestrial biosphere models (TBMs), which aim to simulate the interaction between the atmosphere and the land surface to predict responses of vegetation to changing in environmental forcing. In this review, we examine whether TBMs are able to adequately represent savanna fluxes and what implications potential deficiencies may have for climate change projection scenarios that rely on these models. We start by highlighting the defining characteristic traits and behaviours of savannas, how these differ across continents and how this information is (or is not) represented in the structural framework of many TBMs. We highlight three dynamic processes that we believe directly affect the water use and productivity of the savanna system: phenology, root-water access and fire dynamics. Following this, we discuss how these processes are represented in many current-generation TBMs and whether they are suitable for simulating savanna fluxes.Finally, we give an overview of how eddy-covariance observations in combination with other data sources can be used in model benchmarking and intercomparison frameworks to diagnose the performance of TBMs in this environment and formulate road maps for future development. Our investigation reveals that many TBMs systematically misrepresent phenology, the effects of fire and root-water access (if they are considered at all) and that these should be critical areas for future development. Furthermore, such processes must not be static (i.e. prescribed behaviour) but be capable of responding to the changing environmental conditions in order to emulate the dynamic behaviour of savannas. Without such developments, however, TBMs will have limited predictive capability in making the critical projections needed to understand how savannas will respond to future global change.

Understanding the long-term performance of global satellite leaf area index (LAI) products is important for global change research. However, few effort has been devoted to evaluating the long-term time-series consistencies of LAI products. This study compared four long-term LAI products (GLASS, GLOBMAP, LAI3g, and TCDR) in terms of trends, interannual variabilities, and uncertainty variations from 1982 through 2011. This study also used four ancillary LAI products (GEOV1, MERIS, MODIS C5, and MODIS C6) from 2003 through 2011 to help clarify the performances of the four long-term LAI products. In general, there were marked discrepancies between the four long-term LAI products. During the pre-MODIS period (1982-1999), both linear trends and interannual variabilities of global mean LAI followed the order GLASS>LAI3g>TCDR>GLOBMAP. The GLASS linear trend and interannual variability were almost 4.5 times those of GLOBMAP. During the overlap period (2003-2011), GLASS and GLOBMAP exhibited a decreasing trend, TCDR no trend, and LAI3g an increasing trend. GEOV1, MERIS, and MODIS C6 also exhibited an increasing trend, but to a much smaller extent than that from LAI3g. During both periods, the R2 of detrended anomalies between the four long-term LAI products was smaller than 0.4 for most regions. Interannual variabilities of the four long-term LAI products were considerably different over the two periods, and the differences followed the order GLASS>LAI3g>TCDR>GLOBMAP. Uncertainty variations quantified by a collocation error model followed the same order. Our results indicate that the four long-term LAI products were neither intraconsistent over time nor interconsistent with each other. These inconsistencies may be due to NOAA satellite orbit changes and MODIS sensor degradation. Caution should be used in the interpretation of global changes derived from the four long-term LAI products.

Remote sensing of far-red sun-induced chlorophyll fluorescence (SIF) has emerged as an important tool for studying gross primary productivity (GPP) at the global scale. However, the relationship between SIF and GPP at the canopy scale lacks a clear mechanistic explanation. This is largely due to the poorly characterized role of the relative contributions from canopy structure and leaf physiology to the variability of the top-of-canopy, observed SIF signal. In particular, the effect of the canopy structure beyond light absorption is that only a fraction (fesc) of the SIF emitted from all leaves in the canopy can escape from the canopy due to the strong scattering of near-infrared radiation. We combined rice, wheat and corn canopy-level in-situ datasets to study how the physiological and structural components of SIF individually relate to measures of photosynthesis. At seasonal time scales, we found a considerably strong positive correlation (R2 = 0.4–0.6) of fesc to the seasonal dynamics of the photosynthetic light use efficiency (LUEP), while the estimated physiological SIF yield was almost entirely uncorrelated to LUEP both at seasonal and diurnal time scales, with the partial exception of wheat. Consistent with these findings, the canopy structure and radiation component of SIF, defined as the product of APAR and fesc, explained the relationship of observed SIF to GPP and even outperformed GPP estimation based on observed SIF at two of the three sites investigated. These results held for both half-hourly and daily mean values. In contrast, the total emitted SIF, obtained by normalizing observed SIF for fesc, improved only the relationship to APAR but considerably decreased the correlation to GPP for all three crops. Our findings demonstrate the dominant role of canopy structure in the SIF-GPP relationship and establish a strong, mechanistic link between the near-infrared reflectance of vegetation (NIRV) and the relevant canopy structure information contained in the SIF signal. These insights are expected to be useful in improving remote sensing based GPP estimates.

Global vegetation greening trend stalled after the late 1990s due to increased atmospheric water demand., Atmospheric vapor pressure deficit (VPD) is a critical variable in determining plant photosynthesis. Synthesis of four global climate datasets reveals a sharp increase of VPD after the late 1990s. In response, the vegetation greening trend indicated by a satellite-derived vegetation index (GIMMS3g), which was evident before the late 1990s, was subsequently stalled or reversed. Terrestrial gross primary production derived from two satellite-based models (revised EC-LUE and MODIS) exhibits persistent and widespread decreases after the late 1990s due to increased VPD, which offset the positive CO2 fertilization effect. Six Earth system models have consistently projected continuous increases of VPD throughout the current century. Our results highlight that the impacts of VPD on vegetation growth should be adequately considered to assess ecosystem responses to future climate conditions.

Solar-induced chlorophyll fluorescence (SIF) provides us with new opportunities to understand the physiological and structural dynamics of vegetation from leaf to global scales. However, the relationships between SIF and gross primary productivity (GPP) are not fully understood, which is mainly due to the challenges of decoupling structural and physiological factors that control the relationships. Here, we report the results of continuous observations of canopy-level SIF, GPP, absorbed photosynthetically active radiation (APAR), and chlorophyll: carotenoid index (CCI) in a temperate evergreen needleleaf forest. To understand the mechanisms underlying the relationship between GPP and SIF, we investigated the relationships of light use efficiency (LUEp), chlorophyll fluorescence yield (ΦF), and the fraction of emitted SIF photons escaping from the canopy (fesc) separately. We found a strongly non-linear relationship between GPP and SIF at diurnal and seasonal time scales (R2 = 0.91 with a hyperbolic regression function, daily). GPP saturated with APAR, while SIF did not. Also, there were differential responses of LUEp and ΦF to air temperature. While LUEp reached saturation at high air temperatures, ΦF did not saturate. We found that the canopy-level chlorophyll: carotenoid index was strongly correlated to canopy-level ΦF (R2 = 0.84) implying that ΦF could be more closely related to pigment pool changes rather than LUEp. In addition, we found that the fesc contributed to a stronger SIF-GPP relationship by partially capturing the response of LUEp to diffuse light. These findings can help refine physiological and structural links between canopy-level SIF and GPP in evergreen needleleaf forest.

Satellite image fusion methods that improve spatial and temporal resolution have significant potential to advance understanding of ecosystem dynamics in space and time. However, systematic evaluations of image fusion methods against in situ spectral data are lacking. Here, we used a suite of in situ spectral data collected at 60 elementary sampling units (10 × 10 m) covering 15 Landsat pixel (30 × 30 m) plots and one Moderate Resolution Imaging Spectroradiometer (MODIS) pixel (250 × 250 m) throughout the entire growing season in a heterogeneous rice paddy landscape to evaluate four state-of-the-art image fusion NDVI products. They include the Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM), Flexible Spatiotemporal DAta Fusion (FSDAF), SaTellite dAta IntegRation (STAIR), and the CubeSat Enabled Spatio-Temporal Enhancement Method (CESTEM); the former three blended Landsat and MODIS data, whereas the latter combined CubeSats, Landsat, and MODIS observations. All fusion products showed strong linear relationships against in situ data when combining all spatial and temporal observations (R2: 0.73 to 0.93) although there were partly negative biases (–1% to –9%). These biases resulted from forcing data to image fusion algorithms, such as Landsat (–4%) and MODIS (–7%). Performance difference between fusion methods were considerably larger for spatial than for temporal variation. Furthermore, Landsat NDVI explained only 17–22% of spatial variation against in situ spectral data, which can be translated into weak performance of image fusion products to predict spatial variability in NDVI. Image fusion products that relied on spatial interpolation showed large biases (–15% to –30%) for a vegetation plot surrounded by mixed land cover plots. Our results highlight key sources of uncertainty and will be instrumental in improving satellite image fusion methods to monitor land surface phenology in space and time.

There is an increasing need for approaches to determine reference emission levels and implement policies to address the objectives of Reducing Emissions from Deforestation and Forest Degradation, plus improving forest management, carbon stock enhancement and conservation (REDD+). Important aspects of approaching emissions reductions include coordination and sharing of technology, data, protocols and experiences within and among countries to maximize resources and apply knowledge to build robust monitoring, reporting and verification (MRV) systems. We propose that enhancing the multiple facets of interoperability could facilitate implementation of REDD+ programs and actions. For this case, interoperability is a collective effort with the ultimate goal of sharing and using information to produce knowledge and apply knowledge gained, by removing conceptual, technological, organizational and cultural barriers. These efforts must come from various actors and institutions, including government ministrie...

It is assumed that the ratio of actual evapotranspiration (AET) to potential evapotranspiration (PET) is mostly controlled by the soil water content available for ET. This control is formulated using the soil moisture stress function (SSF), where the evaporative fraction (EF) or the fraction of the AET to PET (fPET) is assumed to be either a linear or a non-linear function of soil moisture. We examine the effectiveness of the soil moisture stress function to quantify soil moisture control on EF or fPET over a dryland wheat field in Victoria, Australia. Micrometeorological observations from two cropping seasons were used for the analysis. The efficacy of a root-density-weighted soil moisture estimate in predicting EF and fPET was investigated as against the commonly assumed fixed-depth root zone soil moisture. However, results indicate a strong relationship between EF and available soil water fraction (AWF) in the root zone only when solar radiation is higher than 5 MJ/m2/day. As the rooting depth increases with vegetation growth, SSF exhibits the strongest correlation with AWF for increasing soil profile depth. In the early and harvesting crop growth stages, ET is constrained mostly by surface soil moisture (0–5 cm). In the mid-growth stages, ET is strongly influenced by soil moisture in the root zone (0–60 cm). The shape of SSF, however, changes significantly between the two years (2012 and 2013). It is inferred that different temporal rainfall patterns between the years caused wheat’s different response to water stress.

Several global gross primary production (GPP) and evapotranspiration (ET) remote sensing products exist, mainly provided by machine-learning (e.g. MPI-BGC) and semi-empirical (e.g. MODIS) approaches. Process-based approaches have the advantage of representing the atmosphere-vegetation-soil system and associated fluxes as an organic integration, but their sophistication results in a lack of high spatiotemporal resolution continuous products. Targeting this gap, we reported a new set of global 8-day composite 1-km resolution GPP and ET products from 2000 to 2015, using a simplified process-based model, the Breathing Earth System Simulator (BESS). BESS couples atmosphere and canopy radiative transfer, photosynthesis and evapotranspiration, and uses MODIS atmosphere and land data and other satellite data sources as inputs. We evaluated BESS products against FLUXNET observations at site scale (total of 113 sites, 742 site years), and against MPI-BGC products at global scale. At site scale, BESS 8-day products agreed with FLUXNET observations with R 2 = 0.67 and RMSE = 2.58 gC m − 2 d − 1 for GPP, and R 2 = 0.62 and RMSE = 0.78 mm d − 1 for ET, respectively, and they captured a majority of seasonal variability, about half of spatial variability, and a minority of interannual variability in FLUXNET observations. At global scale, BESS mean annual sum GPP and ET maps agreed with MPI-BGC products with R 2 = 0.93 and RMSE = 229 gC m − 2 y − 1 for GPP, and R 2 = 0.90 and RMSE = 118 mm y − 1 for ET, respectively. Over the period of 2001–2011, BESS quantified the mean global GPP and ET as 122 ± 25 PgC y − 1 and 65 × 10 3 ± 11 × 10 3 km 3 y − 1 , respectively, with a significant ascending GPP trend by 0.27 PgC y − 2 ( p

Digital cover photography (DCP) has emerged as an indirect method to measure gap fraction of vegetation canopies. However, as with other photographic methods, determining camera relative exposure value (REV) and threshold for pixel classification, cause substantial uncertainties in gap fraction estimates. Here we propose a new method to improve the measurement of gap fraction under various solar zenith angles (SZAs), sky conditions, and canopy structures. This method computes gap fractions of ambiguous vegetation or sky pixels using an unsaturated raw image from DCP and a reconstructed sky image from the raw image, thus taking full advantage of the potential of raw image processing. This is combined with pre-classification of pixels that are unambiguously canopy and sky to greatly reduce light scattering effects that are likely to be present within the canopy. To test the sensitivity of the new method, we acquired images at one-hour intervals between 20 and 85° of SZAs under closed, half-closed, and open canopies with REV settings from 0 to −5. The new method showed little variation in gap fractions across the diverse SZAs in closed, half-closed, and open canopies. A perforated panel experiment, which was used to test the accuracy of the estimated gap fractions, confirmed that the new method accurately estimated gap fractions across a range of hole size, gap fractions and SZAs. We conclude that the new method opens new opportunities to estimate gap fractions accurately regardless of solar positions from open to closed canopies, and is a significant advance for accurate and precise monitoring of canopy cover and leaf area index (LAI), and for calibrating and evaluating satellite remote sensing LAI products.

Clumping index quantifies the level of foliage aggregation, relative to a random distribution, and is a key structural parameter of plant canopies and is widely used in ecological and meteorological models. In this study, the inter- and intra-annual variations in clumping index values, derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) BRDF product, are investigated at six forest sites, including conifer forests, a mixed deciduous forest and an oak-savanna system. We find that the clumping index displays large seasonal variation, particularly for the deciduous sites, with the magnitude in clumping index values at each site comparable on an intra-annual basis, and the seasonality of clumping index well captured after noise removal. For broadleaved and mixed forest sites, minimum clumping index values are usually found during the season when leaf area index is at its maximum. The magnitude of MODIS clumping index is validated by ground data collected from 17 sites. Validation shows that the MODIS clumping index can explain 75% of variance in measured values (bias = 0.03 and rmse = 0.08), although with a narrower amplitude in variation. This study suggests that the MODIS BRDF product has the potential to produce good seasonal trajectories of clumping index values, but with an improved estimation of background reflectance.

Vegetation monitoring can be used to detect CO2 leakage in carbon capture and storage (CCS) sites because it can monitor a large area at a relatively low cost. However, a rapidly responsive, sensitive, and cost-effective plant parameters must be suggested for vegetation monitoring to be practically utilized as a CCS management strategy. To screen the proper plant parameters for leakage monitoring, a greenhouse experiment was conducted by exposing kale (Brassica oleracea var. viridis), a sensitive plant, to 10%, 20%, and 40% soil CO2 concentrations. Water and water with CO2 stress treatments were also introduced to examine the parameters differentiating CO2 stress from water stresses. We tested the hypothesis that chlorophyl fluorescence parameters would be early and sensitive indicator to detect CO2 leakage. The results showed that the fluorescence parameters of effective quantum yield of photosystem II (Y(II)), detected the difference between CO2 treatments and control earlier than any other parameters, such as chlorophyl content, hyperspectral vegetation indices, and biomass. For systematic comparison among many parameters, we proposed an indicator evaluation score (IES) method based on four categories: CO2 specificity, early detection, field applicability, and cost. The IES results showed that fluorescence parameters (Y(II)) had the highest IES scores, and the parameters from spectral sensors (380–800 nm wavelength) had the second highest values. We suggest the IES system as a useful tool for evaluating new parameters in vegetation monitoring.

Rapid urbanization has created large areas of impervious surface areas globally. As there is little carbon input by plants into soils under impervious surfaces, soil organic carbon (SOC) stocks under impervious surfaces generally have been overlooked in the urban carbon budgets. Here we investigate SOC stocks under impervious concrete surfaces and vegetative surfaces across soil profiles to a depth of 5 m in urban housing complexes in Seoul, Republic of Korea. In the top 1 m of the profile, SOC stocks under vegetative surfaces were three times greater than those under impervious surfaces. However, we discovered that unexpectedly high SOC stocks appeared in deeper soil layers under both surface types, which led to comparable SOC stocks at a depth of 5 m beneath the impervious surface (16.9 ± 1.9 kgC m−2) and at the vegetative surface (22.3 ± 2.2 kgC m−2). Consequently, the ratio of SOC stocks at depths of 1 m to 5 m were 16% in impervious surfaces and 34% in vegetative surfaces, suggesting conventional soil sampling at 1 m depth could miss large SOC. Stable isotope data ( δ 13C and δ 15N) combined with historical aerial photographs revealed that cropland that existed until the 1970s formed the high SOC cultural layer in deeper soils. Our results highlight that deep soils under impervious surfaces could be overlooked carbon hotspots in urban ecosystems. We believe this finding could help city planners and policy makers to assess regional carbon budgets and to reduce carbon footprint by recycling the deep SOC excavated from various construction projects towards sustainable urban development.

Addressing the impact of vegetation architecture on shortwave radiation transfer in land surface models is important for accurate weather forecasting, carbon budget estimates, and climate predictions. This paper investigates to what extent it is possible to retrieve structural parameters of two different parameterization schemes from direct transmittance derived from digital hemispherical photography and 3D radiative transfer modeling for two study sites with different vegetation canopy architectures. Neglecting the representation of 3D canopy structure in radiative transfer schemes leads to significant errors in shortwave radiation partitioning (up to 3.5 times more direct transmittance in the 3D model). Structural parameters, referred to as whole canopy ‘clumping indices’, were obtained in order to evaluate the impact of angular variation in clumping on shortwave radiation transfer. Impacts on photosynthesis were evaluated at site level with the UKESM land surface model, JULES. A comparison between flux tower derived and modeled photosynthesis indicates that considering zenith angular variations of structural parameters in the radiative transfer scheme of the UKESM land surface model significantly improves photosynthesis prediction in light limited ecosystems (from RMSE = 2.91 μmol CO2.m − 2.s − 1 to RMSE = 1.51 μmol CO2.m − 2.s − 1, 48% smaller), typically with enhanced photosynthesis from bottom layers.

Rice is a major crop that feeds more than half of the global population. Rice paddies play a complex role in the carbon cycle by emitting CH4 to the atmosphere while either sequestering or releasing CO2 to the atmosphere. Here, we present 3.5 years of eddy covariance measurements of CO2 and CH4 fluxes over an intermittently irrigated, single crop rice paddy in South Korea in tandem with carbon stock measurements for leaves, stems, grains, and roots. The rice paddy acted as a slight sink or neutral in CO2 (−47 ± 51 g C m−2 y−1; mean ± 95% confidence interval) with 941 ± 130 g C m−2 y−1 and 894 ± 108 g C m−2 y−1 as the gross primary production (GPP) and ecosystem respiration (Reco), respectively. The mean annual net primary productivity (NPP) was 666 ± 31 g C m−2 y−1, which was partitioned to leaves (9%), stems (30%), grains (51%), and roots (10%), for a harvest index of 57%. The paddy emitted 20.6 ± 1.5 g C m−2 y−1 in the form of CH4, of which between 82 and 96% was emitted during the cultivation period (~132 days) . The CH4 flux showed a mid-season decline caused by mid-season drainage, which produced aerobic soil conditions that inhibited CH4 emissions. The proportion of estimated ebullitive flux to total CH4 flux ranged from 10 to 17% during the cultivation period. GPP, Reco, and CH4 fluxes varied by 12%, 11%, and 7%, respectively, over the 3-year period. When considering the harvest and the global warming potential of CH4, the paddy shifted from a greenhouse gas sink to a strong source (2,741 ± 468 g CO2 eq m−2 y−1). These findings will improve our understanding of carbon budgets and cycles, and help develop greenhouse gas and water management strategies for rice paddies.

International audience; Terrestrial gross primary productivity (GPP) plays an essential role in the global carbon cycle, but the quantification of the spatial and temporal variations in photosynthesis is still largely uncertain. Our work aimed to investigate the potential of remote sensing to provide new insights into plant photosynthesis at a fine spatial resolution. This goal was achieved by exploiting high-resolution images acquired with the FLuorescence EXplorer (FLEX) airborne demonstrator HyPlant. The sensor was flown over a mixed forest, and the images collected were elaborated to obtain two independent indicators of plant photosynthesis. First, maps of sun-induced chlorophyll fluorescence (F), a novel indicator of plant photosynthetic activity, were successfully obtained at both the red and far-red peaks (r2 = 0.89 and p 0.05). The spatial relationships found at high resolution provide valuable insight into the critical role of spatial heterogeneity in controlling the relationship between the far-red F and the GPP, indicating the need to consider this heterogeneity at a coarser resolution.

Statistical gross primary productivity (GPP) estimation from remote sensing observations has mostly been attempted on the basis of multispectral observations. To make full use of the information contained in vegetation spectra, however, hyperspectral observations should be used in combination with appropriate multivariate methods. Nevertheless, only very few previous studies attempted to estimate GPP directly from hyperspectral observations and did so on the basis of reflectance, with only a limited number of temporally discontinuous observations. In this study, we used long-term, continuous, half-hourly hyperspectral observations covering the visible and near-infrared spectral range to estimate GPP directly from upwelling irradiance using partial least square (PLS) regression in a rice paddy. To gain a better understanding of processes underlying the PLS estimation, we used extensive complementary field observations to run process-based simulations using the SCOPE model. We then applied PLS regression to the simulated hyperspectral data in the same way as for the observations and disentangled contributions related to relevant physiological processes, namely sun-induced chlorophyll fluorescence (SIF) and xanthophyll cycle-related spectral changes (XC). We found that upwelling hyperspectral irradiance in the visible and near-infrared spectral range predicted GPP better than reflectance. Furthermore, PLS-based GPP estimates outperformed both far-red SIF and widely used vegetation index-based methods. However, the most relevant information for the observation-based PLS-models was not clearly related to XC or SIF as the near-infrared spectral range showed comparable performance. Also, the simple average of upwelling irradiance over the 850–900 nm range outperformed the other non-multivariate approaches, including far-red SIF. These results held for the evaluation in terms of the seasonal variation of GPP, while there was apparently a small contribution of SIF and XC for the diurnal variation. The simulation-based analysis showed that SIF and XC contributed useful information to both GPP and photosynthetic light use efficiency (LUE) estimates at both seasonal and diurnal time scales. The strongest unique contribution from either SIF or XC, however, was to the diurnal variation of GPP and XC showed considerably better performance than SIF. We did not find improvements when combining the spectral regions of XC (500–570 nm) and SIF (650–800 nm) to estimate GPP. While SIF showed improvements when combined with the remaining spectral information excluding XC, this was not the case for XC. Our approach combines the strengths of process-based modeling with multivariate statistical analysis to improve our understanding of the usable information content in vegetation spectra and is highly relevant for further developing suitable methods for GPP estimation at large scales.

Solar-induced chlorophyll fluorescence (SIF) has emerged as a leading approach for remote sensing of gross primary productivity (GPP). While SIF has an intrinsic, underlying relationship with canopy light capture and light use efficiency, these physiological relationships are obscured by the fact that satellites observe a small and variable fraction of total emitted canopy SIF. Upon emission, most SIF photons are reabsorbed or scattered within the canopy, preventing their observation remotely. The complexities of the radiative transfer process, which vary across time and space, limit our ability to reliably infer physiological processes from SIF observations. Here, we propose an approach for estimating the fraction of total emitted near-infrared SIF (760 nm) photons that escape the canopy by combining the near-infrared reflectance of vegetation (NIRV) and the fraction of absorbed photosynthetically active radiation (fPAR), two widely available remote sensing products. Our approach relies on the fact that NIRV is resilient against soil background contamination, allowing us to reliably calculate the bidirectional reflectance factor of vegetation, which in turn conveys information about the escape ratio of SIF photons. Our NIRV-based approach explains variations in the escape ratio with an R2 of 0.91 and an RMSE of 1.48% across a series of simulations where canopy structure, soil brightness, and sun-sensor-canopy geometry are varied. The approach is applicable to conditions of low leaf area index and fractional vegetation cover. We show that correcting for the escape ratio of SIF using NIRV provides robust estimates of total emitted SIF, providing for the possibility of studying physiological variations of fluorescence yield at the global scale.

Clumping index, the measure of foliage grouping relative to a random distribution of leaves in space, is a key structural parameter of plant canopies that influences canopy radiation regimes and controls canopy photosynthesis and other land–atmosphere interactions. In this study, we retrieve the clumping index using the original 275 m resolution data of the Multi-angle Imaging SpectroRadiometer (MISR) instrument over a set of sites representing diverse biomes and different canopy structures. Also for the first time, the MISR derived clumping index values are directly validated with both in-situ vertical profiles and seasonal trajectories of clumping index. Our results illustrate that MISR data with 275 m allow clumping index estimates at much more pertinent scales (both spatial and temporal) than previous maps from Polarization and Directionality of Earth Reflectances (POLDER) and Moderate Resolution Imaging Spectroradiometer (MODIS) for modeling local carbon and energy fluxes.

The world is experiencing remarkable climate change, which alters vegetation structural and functional status, and the biosphere feedback to the climate might amplify or dampen regional and global climate change. The exchange of energy and mass across the land-atmosphere interface is an essential indicator of the interaction between ecosystem and environment, and it has been widely measured at landscape scale since the establishment of FLUXNET, a global network of micrometeorological flux measurement sites (Baldocchi, 2008; Baldocchi et al., 2001). Furthermore, these site observations have been upscaled to global scale in conjunction with satellite remote sensing data, providing opportunities to investigate the spatial and temporal variations of carbon and water cycles from a macro perspective (Beer et al., 2010; Martin Jung et al., 2010). However, such upscaling approaches are based on data-driven models, which need to be well calibrated using site observations (M. Jung, Reichstein, & Bondeau, 2009; Papale & Valentini, 2003; Xiao et al., 2010; Yang et al., 2007). Consequently, those models are limited by the representativeness, quantity and quality of the training datasets (Sundareshwar et al., 2006), as well as the lack of independent reference datasets for validation. Furthermore, by using the upscaling approach predictors are directly linked with target fluxes through machine learning, while the intermediate variables along with intrinsic mechanisms are hidden.

Savanna ecosystems are one of the most dominant and complex terrestrial biomes that derives from a distinct vegetative surface comprised of co-dominant tree and grass populations. While these two vegetation types co-exist functionally, demographically they are not static, but are dynamically changing in response to environmental forces such as annual fire events and rainfall variability. Modelling savanna environments with the current generation of terrestrial biosphere models (TBMs) has presented many problems, particularly describing fire frequency and intensity, phenology, leaf biochemistry of C3 and C4 photosynthesis vegetation, and root water uptake. In order to better understand why TBMs perform so poorly in savannas, we conducted a model inter-comparison of 6 TBMs and assessed their performance at simulating latent energy (LE) and gross primary productivity (GPP) for five savanna sites along a rainfall gradient in northern Australia. Performance in predicting LE and GPP was measured using an empirical benchmarking system, which ranks models by their ability to utilise meteorological driving information to predict the fluxes. On average, the TBMs performed as well as a multi-linear regression of the fluxes against solar radiation, temperature and vapour pressure deficit, but were outperformed by a more complicated nonlinear response model that also included the leaf area index (LAI). This identified that the TBMs are not fully utilising their input information effectively in determining savanna LE and GPP, and highlights that savanna dynamics cannot be calibrated into models and that there are problems in underlying model processes. We identified key weaknesses in a model's ability to simulate savanna fluxes and their seasonal variation, related to the representation of vegetation by the models and root water uptake. We underline these weaknesses in terms of three critical areas for development. First, prescribed tree-rooting depths must be deep enough, enabling the extraction of deep soil water stores to maintain photosynthesis and transpiration during the dry season. Second, models must treat grasses as a co-dominant interface for water and carbon exchange, rather than a secondary one to trees. Third, models need a dynamic representation of LAI that encompasses the dynamic phenology of savanna vegetation and its response to rainfall interannual variability. We believe this study is the first to assess how well TBMs simulate savanna ecosystems, and that these results will be used to improve the representation of savannas ecosystems in future global climate model studies.

The savanna complex is a highly diverse global biome that occurs within the seasonally dry tropical to sub-tropical equatorial latitudes. Savannas are open-canopy environments that encompass a broad demographic continuum, often characterised by a dynamically changing dominance between C3-tree and C4-grass vegetation, where frequent environmental disturbances such as fire modulates the balance between ephemeral and perennial life forms. Climate change is projected to result in significant changes to the savanna floristic structure, with increases to woody biomass expected through CO2 fertilisation in mesic savannas and increased tree mortality expected through increased rainfall interannual variability in xeric savannas. The complex interaction between vegetation and climate that occurs in savannas has traditionally challenged current-generation terrestrial biosphere models (TBMs), which aim to simulate the interaction between the atmosphere and the land-surface to predict responses of vegetation to changing in environmental forcing. In this review, we examine whether TBMs are able to adequately represent savanna dynamics and what implications potential deficiencies may have for climate change projection scenarios that rely on these models. We start by highlighting the defining characteristic traits and behaviours of savanna, how these differ across continents, and how this information is (or is not) represented in the structural framework of many TBMs. We highlight three dynamic processes that we believe directly affect the water-use and productivity of the savanna system, namely: phenology; root-water access; and fire dynamics. Following this, we discuss how these processes are represented in many current generation TBMs and whether they are suitable for simulating savanna dynamics. Finally, we give an overview of how eddy-covariance observations in combination with other data sources, can be used in model benchmarking and inter-comparison frameworks to diagnose the performance of TBMs in this environment and formulate roadmaps for future development. Our investigation reveals that many TBMs systematically misrepresent phenology, effects of fire and root-water access (if they are considered at all) and that these should be critical areas for future development. Furthermore, such processes must not be static (i.e. prescribed behaviour), but be capable of responding to the changing environmental conditions in order to emulate the dynamic behaviour of savannas. Without such developments, however, TBMs will have limited predictive capability in making the critical projections needed to understand how savannas will respond to future global change.

The effectiveness of using leveled digital camera for measuring leaf inclination angles was investigated in this study as an inexpensive and convenient alternative to existing approaches. The new method is validated with manual leaf angle measurements for various broadleaf tree species common to hemi-boreal region of Estonia and the tropical forests of Hawai’i Islands. The acquired leaf angle distributions suggest that planophile case might be more appropriate than the commonly assumed spherical as the general approximation of leaf orientation while modeling the radiation transmission through the canopies of (hemi)-boreal broadleaf stands. However, direct leaf inclination measurements should be obtained whenever possible, as there will always exist a large variety of leaf orientation, both among different species and in the space–time domain within a single species. The camera method tested in this study provides a new robust and affordable tool to obtain this information.

The El Niño-Southern Oscillation exerts a large influence on global climate regimes and on the global carbon cycle. Although El Niño is known to be associated with a reduction of the global total land carbon sink, results based on prognostic models or measurements disagree over the relative contribution of photosynthesis to the reduced sink. Here, we provide an independent remote sensing-based analysis on the impact of the 2015–2016 El Niño on global photosynthesis using six global satellite-based photosynthesis products and a global solar-induced fluorescence (SIF) dataset. An ensemble of satellite-based photosynthesis products showed a negative anomaly of −0.7 ± 1.2 PgC in 2015, but a slight positive anomaly of 0.05 ± 0.89 PgC in 2016, which when combined with observations of the growth rate of atmospheric carbon dioxide concentrations suggests that the reduction of the land residual sink was likely dominated by photosynthesis in 2015 but by respiration in 2016. The six satellite-based products unanimously identified a major photosynthesis reduction of −1.1 ± 0.52 PgC from savannahs in 2015 and 2016, followed by a highly uncertain reduction of −0.22 ± 0.98 PgC from rainforests. Vegetation in the Northern Hemisphere enhanced photosynthesis before and after the peak El Niño, especially in grasslands (0.33 ± 0.13 PgC). The patterns of satellite-based photosynthesis ensemble mean were corroborated by SIF, except in rainforests and South America, where the anomalies of satellite-based photosynthesis products also diverged the most. We found the inter-model variation of photosynthesis estimates was strongly related to the discrepancy between moisture forcings for models. These results highlight the importance of considering multiple photosynthesis proxies when assessing responses to climatic anomalies. This article is part of a discussion meeting issue ‘The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'.

A growing literature is reporting on how the terrestrial carbon cycle is experiencing year-to-year variability because of climate anomalies and trends caused by global change. As CO2 concentration records in the atmosphere exceed 50 years and as satellite records reach over 30 years in length, we are becoming better able to address carbon cycle variability and trends. Here we review how variable the carbon cycle is, how large the trends in its gross and net fluxes are, and how well the signal can be separated from noise. We explore mechanisms that explain year-to-year variability and trends by deconstructing the global carbon budget. The CO2 concentration record is detecting a significant increase in the seasonal amplitude between 1958 and now. Inferential methods provide a variety of explanations for this result, but a conclusive attribution remains elusive. Scientists have reported that this trend is a consequence of the greening of the biosphere, stronger northern latitude photosynthesis, more photosynthesis by semi-arid ecosystems, agriculture and the green revolution, tropical temperature anomalies, or increased winter respiration. At the global scale, variability in the terrestrial carbon cycle can be due to changes in constituent fluxes, gross primary productivity, plant respiration and heterotrophic (microbial) respiration, and losses due to fire, land use change, soil erosion, or harvesting. It remains controversial whether or not there is a significant trend in global primary productivity (due to rising CO2, temperature, nitrogen deposition, changing land use, and preponderance of wet and dry regions). The degree to which year-to-year variability in temperature and precipitation anomalies affect global primary productivity also remains uncertain. For perspective, interannual variability in global gross primary productivity is relatively small (on the order of 2 Pg-C y-1) with respect to a large and uncertain background (123 +/- 4 Pg-C y-1), and detected trends in global primary productivity are even smaller (33 Tg-C y-2). Yet residual carbon balance methods infer that the terrestrial biosphere is experiencing a significant and growing carbon sink. Possible explanations for this large and growing net land sink include roles of land use change and greening of the land, regional enhancement of photosynthesis, and down regulation of plant and soil respiration with warming temperatures. Longer time series of variables needed to provide top-down and bottom-up assessments of the carbon cycle are needed to resolve these pressing and unresolved issues regarding how, why, and at what rates gross and net carbon fluxes are changing.

Reliable and fine resolution estimates of surface net-radiation are required for estimating latent and sensible heat fluxes between the land surface and the atmosphere. However, currently, fine resolution estimates of net-radiation are not available and consequently it is challenging to develop multi-year estimates of evapotranspiration at scales that can capture land surface heterogeneity and are relevant for policy and decision-making. We developed and evaluated a global net-radiation product at 5 km and 8-day resolution by combining mutually consistent atmosphere and land data from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra. Comparison with net-radiation measurements from 154 globally distributed sites (414 site-years) from the FLUXNET and Surface Radiation budget network (SURFRAD) showed that the net-radiation product agreed well with measurements across seasons and climate types in the extratropics (Wilmott's index ranged from 0.74 for boreal to 0.63 for Mediterranean sites). Mean absolute deviation between the MODIS and measured net-radiation ranged from 38.0 +/- 1.8 W.m(-2) in boreal to 72.0 +/- 4.1 W.m(-2) in the tropical climates. The mean bias was small and constituted only 11%, 0.7%, 8.4%, 4.2%, 13.3%, and 5.4% of the mean absolute error in daytime net-radiation in boreal, Mediterranean, temperate-continental, temperate, semi-arid, and tropical climate, respectively. To assess the accuracy of the broader spatiotemporal patterns, we upscaled error-quantified MODIS net-radiation and compared it with the net-radiation estimates from the coarse spatial (1 degrees x 1 degrees) but high temporal resolution gridded net-radiation product from the Clouds and Earth's Radiant Energy System (CERES). Our estimates agreed closely with the net-radiation estimates from the CERES. Difference between the two was less than 10 Wm(-2) in 94% of the total land area. MODIS net-radiation product will be a valuable resource for the science community studying turbulent fluxes and energy budget at the Earth's surface.

International audience; Interannual variability in biosphere-atmosphere exchange of CO2 is driven by a diverse range of biotic and abiotic factors. Replicating this variability thus represents the ‘acid test’ for terrestrial biosphere models. Although such models are commonly used to project responses to both normal and anomalous variability in climate, they are rarely tested explicitly against inter-annual variability in observations. Herein, using standardized data from the North American Carbon Program, we assess the performance of 16 terrestrial biosphere models and 3 remote sensing products against long-term measurements of biosphere-atmosphere CO2 exchange made with eddy-covariance flux towers at 11 forested sites in North America. Instead of focusing on model-data agreement we take a systematic, variability-oriented approach and show that although the models tend to reproduce the mean magnitude of the observed annual flux variability, they fail to reproduce the timing. Large biases in modeled annual means are evident for all models. Observed interannual variability is found to commonly be on the order of magnitude of the mean fluxes. None of the models consistently reproduce observed interannual variability within measurement uncertainty. Underrepresentation of variability in spring phenology, soil thaw and snowpack melting, and difficulties in reproducing the lagged response to extreme climatic events are identified as systematic errors, common to all models included in this study.

a b s t r a c t The regular monitoring of evapotranspiration from satellites has been limited because of discontinu- ous temporal coverage, resulting in snapshots at a particular point in space and time. We developed a temporal upscaling scheme using satellite-derived instantaneous estimates of evapotranspiration to produce a daily-sum evapotranspiration averaged over an 8-day interval. We tested this scheme against measured evapotranspiration data from 34 eddy covariance flux towers covering seven plant functional types from boreal to tropical climatic zones. We found that the ratio of a half-hourly-sum of potential solar radiation (extraterrestrial solar irradiance on a plane parallel to the Earth's surface) between 10:00 hh and 14:00 hh to a daily-sum of potential solar radiation provides a robust scaling factor to convert a half-hourly measured evapotranspiration to an estimate of a daily-sum; the estimated and measured daily sum evapotranspiration showed strong linear relation (r 2 = 0.92) and small bias (−2.7%). By com- parison, assuming a constant evaporative fraction (the ratio of evapotranspiration to available energy) during the daytime, although commonly used for temporal upscaling, caused 13% underestimation of evapotranspiration on an annual scale. The proposed temporal upscaling scheme requires only latitude, longitude and time as input. Thus it will be useful for developing continuous evapotranspiration estimates in space and time, which will improve continuous monitoring of hydrological cycle from local to global scales.

International audience; We propose the Breathing Earth System Simulator (BESS), an upscaling approach to quantify global gross primary productivity and evapotranspiration using MODIS with a spatial resolution of 1-5 km and a temporal resolution of 8 days. This effort is novel because it is the first system that harmonizes and utilizes MODIS Atmosphere and Land products on the same projection and spatial resolution over the global land. This enabled us to use the MODIS Atmosphere products to calculate atmospheric radiative transfer for visual and near infrared radiation wave bands. Then we coupled atmospheric and canopy radiative transfer processes, with models that computed leaf photosynthesis, stomatal conductance and transpiration on the sunlit and shaded portions of the vegetation and soil. At the annual time step, the mass and energy fluxes derived from BESS showed strong linear relations with measurements of solar irradiance (r(2) = 0.95, relative bias: 8%), gross primary productivity (r(2) = 0.86, relative bias: 5%) and evapotranspiration (r(2) = 0.86, relative bias: 15%) in data from 33 flux towers that cover seven plant functional types across arctic to tropical climatic zones. A sensitivity analysis revealed that the gross primary productivity and evapotranspiration computed in BESS were most sensitive to leaf area index and solar irradiance, respectively. We quantified the mean global terrestrial estimates of gross primary productivity and evapotranpiration between 2001 and 2003 as 118 +/- 26 PgC yr(-1) and 500 +/- 104 mm yr(-1) (equivalent to 63,000 +/- 13,100 km(3) yr(-1)), respectively. BESS-derived gross primary productivity and evapotranspiration estimates were consistent with the estimates from independent machine-learning, data-driven products, but the process-oriented structure has the advantage of diagnosing sensitivity of mechanisms. The process-based BESS is able to offer gridded biophysical variables everywhere from local to the total global land scales with an 8-day interval over multiple years.