Global Gross Primary Production Research Articles

Why models underestimate West African tropical forest primary productivity

Tropical forests dominate terrestrial photosynthesis, yet there are major contradictions in our understanding due to a lack of field studies, especially outside the tropical Americas. A recent field study indicated that West African forests have among the highest forests gross primary productivity (GPP) yet observed, contradicting models that rank them lower than Amazonian forests. Here, we show possible reasons for this data-model mismatch. We found that biometric GPP measurements are on average 56.3% higher than multiple global GPP products at the study sites. The underestimation of GPP largely disappears when a standard photosynthesis model is informed by local field-measured values of (a) fractional absorbed photosynthetic radiation (fAPAR), and (b) photosynthetic traits. Remote sensing products systematically underestimate fAPAR (33.9% on average at study sites) due to cloud contamination issues. The study highlights the potential widespread underestimation of tropical forests GPP and carbon cycling and hints at the ways forward for model and input data improvement.

Terrestrial photosynthesis inferred from plant carbonyl sulfide uptake.

Terrestrial photosynthesis, or gross primary production (GPP), is the largest carbon flux in the biosphere, but its global magnitude and spatiotemporal dynamics remain uncertain1. The global annual mean GPP is historically thought to be around 120 PgC yr-1 (refs. 2-6), which is about 30-50 PgC yr-1 lower than GPP inferred from the oxygen-18 (18O) isotope7 and soil respiration8. This disparity is a source of uncertainty in predicting climate-carbon cycle feedbacks9,10. Here we infer GPP from carbonyl sulfide, an innovative tracer for CO2 diffusion from ambient air to leaf chloroplasts through stomata and mesophyll layers. We demonstrate that explicitly representing mesophyll diffusion is important for accurately quantifying the spatiotemporal dynamics of carbonyl sulfide uptake by plants. From the estimate of carbonyl sulfide uptake by plants, we infer a global contemporary GPP of 157 (±8.5) PgC yr-1, which is consistent with estimates from 18O (150-175 PgC yr-1) and soil respiration ( PgC yr-1), but with an improved confidence level. Our global GPP is higher than satellite optical observation-driven estimates (120-140 PgC yr-1) that are used for Earth system model benchmarking. This difference predominantly occurs in the pan-tropical rainforests and is corroborated by ground measurements11, suggesting a more productive tropics than satellite-based GPP products indicated. As GPP is a primary determinant of terrestrial carbon sinks and may shape climate trajectories9,10, our findings lay a physiological foundation on which the understanding and prediction of carbon-climate feedbacks can be advanced.

Estimating Global Gross Primary Production Using an Improved MODIS Leaf Area Index Dataset

Remote sensing and process-coupled ecological models are widely used for the simulation of GPP, which plays a key role in estimating and monitoring terrestrial ecosystem productivity. However, most such models do not differentiate the C3 and C4 photosynthetic pathways and neglect the effect of nitrogen content on Vmax and Jmax, leading to considerable bias in the estimation of gross primary productivity (GPP). Here, we developed a model driven by the leaf area index, climate, and atmospheric CO2 concentration to estimate global GPP with a spatial resolution of 0.1° and a temporal interval of 1 day from 2000 to 2022. We validated our model with ground-based GPP measurements at 128 flux tower sites, which yielded an accuracy of 72.3%. We found that the global GPP ranged from 116.4 PgCyear−1 to 133.94 PgCyear−1 from 2000 to 2022, with an average of 125.93 PgCyear−1. We also found that the global GPP showed an increasing trend of 0.548 PgCyear−1 during the study period. Further analyses using the structure equation model showed that atmospheric CO2 concentration and air temperature were the main drivers of the global GPP changes, total associations of 0.853 and 0.75, respectively, while precipitation represented a minor but negative contribution to global GPP.

Improving global gross primary productivity estimation using two-leaf light use efficiency model by considering various environmental factors via machine learning

Distinguishing gross primary productivity (GPP) into sunlit (GPPsu) and shaded (GPPsh) components is critical for understanding the carbon exchange between the atmosphere and terrestrial ecosystems under climate change. Recently, the two-leaf light use efficiency (TL-LUE) model has proven effective for simulating global GPPsu and GPPsh. However, no known physical method has focused on integrating the overall constraint of intricate environmental factors on photosynthetic capability, and seasonal differences in the foliage clumping index (CI), which most likely influences GPP estimation in LUE models. Here, we propose the TL-CRF model, which uses the random forest technique to integrate various environmental variables, particularly for terrestrial water storage (TWS), into the TL-LUE model. Moreover, we consider seasonal differences in CI at a global scale. Based on 267 global eddy covariance flux sites, we explored the functional response of vegetation photosynthesis to key environmental factors, and trained and evaluated the TL-CRF model. The TL-CRF model was then used to simulate global eight-day GPP, GPPsu, and GPPsh from 2002 to 2020. The results show that the relative prediction error of environmental stress factors on the maximum LUE is reduced by approximately 52 % when these factors are integrated via the RF model. Thus the accuracy of global GPP estimation (R2 = 0.87, RMSE = 0.94 g C m−2 d−1, MAE = 0.61 g C m−2 d−1) in the TL-CRF model is greater than that (R2 = 0.76, RMSE = 2.18 g C m−2 d−1, MAE = 1.50 g C m−2 d−1) in the TL-LUE model, although this accuracy awaits further investigation among the released GPP products. TWS exerts the greatest control over ecosystem photosynthesis intensity, making it a suitable water indicator. Furthermore, the results confirm an optimal minimum air temperature for photosynthesis. Overall, these findings indicate a promising method for producing a new global GPP dataset, advancing our understanding of the dynamics and interactions between photosynthesis and environmental factors.

A 2001–2022 global gross primary productivity dataset using an ensemble model based on the random forest method

Abstract. Advancements in remote sensing technology have significantly contributed to the improvement of models for estimating terrestrial gross primary productivity (GPP). However, discrepancies in the spatial distribution and interannual variability within GPP datasets pose challenges to a comprehensive understanding of the terrestrial carbon cycle. In contrast to previous models that rely on remote sensing and environmental variables, we developed an ensemble model based on the random forest method (denoted ERF model). This model used GPP outputs from established models: Eddy Covariance Light Use Efficiency (EC-LUE), GPP estimate model based on Kernel Normalized Difference Vegetation Index (GPP-kNDVI), GPP estimate model based on Near-Infrared Reflectance of Vegetation (GPP-NIRv), Revised-EC-LUE, Vegetation Photosynthesis Model (VPM), and GPP estimate model based on the Moderate Resolution Imaging Spectroradiometer (MODIS). These outputs were used as inputs to estimate GPP. The ERF model demonstrated superior performance, explaining 85.1 % of the monthly GPP variations at 170 sites and surpassing the performance of selected GPP estimate models (67.7 %–77.5 %) and an independent random forest model using remote sensing and environmental variables (81.5 %). Additionally, the ERF model improved accuracy across each month and with various subranges, mitigating the issue of “high-value underestimation and low-value overestimation” in GPP estimates. Over the period from 2001 to 2022, the global GPP estimated by the ERF model was 132.7 PgC yr−1, with an increasing trend of 0.42 PgC yr−2, which is comparable to or slightly better than the accuracy of other mainstream GPP datasets in terms of validation results of GPP observations independent of FLUXNET (i.e., ChinaFLUX). Importantly, for a growing number of GPP datasets, our study provides a way to integrate these GPP datasets, which may lead to a more reliable estimate of global GPP.

Increasing exposure of global croplands productivity to growing season heatwaves under climate warming

Abstract Growing season heatwaves that occur simultaneously over global croplands can negatively impact global food baskets. The long-term changes of growing season heatwaves, as well as their impacts on croplands productivity, are crucial to food security, but remain unclear. Here, we investigated changes in the frequency, intensity and magnitude of growing season heatwaves from the past to the future over the global croplands, based on observations and Coupled Model Intercomparison Project Phase 6 models. We introduced an index, gross primary productivity (GPP) exposure, as a proxy of the overall impact of heatwaves on cropland productivity. The results show that the frequency and intensity of growing season heatwaves have increased since 1950 and will continue throughout the 21st century. The increase of the annual accumulated magnitude of growing season heatwaves in the future is mainly contributed by the increase of heatwave frequency. This leads to a global-scale increase in the GPP exposure to growing season heatwaves, with Asia, North America, and Europe being the most affected. The continued increase in GPP exposure is dominated by increases in heatwaves rather than GPP itself. Under the lower emission scenario SSP1-2.6, the global cropland GPP exposure will reduce by 86.11% and 330.47% relative to that under SSP2-4.5 and SSP5-8.5 scenarios, respectively, by the end of 21st century. Our results provide crucial insights into potential impacts of heatwaves on cropland productivity and hence food security.

Emergent constraints on historical and future global gross primary productivity.

Terrestrial gross primary productivity (GPP) is the largest carbon flux in the global carbon cycle and plays a crucial role in terrestrial carbon sequestration. However, historical and future global GPP estimates still vary markedly. In this study, we reduced uncertainties in global GPP estimates by employing an innovative emergent constraint method on remote sensing-based GPP datasets (RS-GPP), using ground-based estimates of GPP from flux towers as the observational constraint. Using this approach, the global GPP in 2001-2014 was estimated to be 126.8 ± 6.4 PgC year-1, compared to the original RS-GPP ensemble mean of 120.9 ± 10.6 PgC year-1, which reduced the uncertainty range by 39.6%. Independent space- and time-based (different latitudinal zones, different vegetation types, and individual year) constraints further confirmed the robustness of the global GPP estimate. Building on these insights, we extended our constraints to project global GPP estimates in 2081-2100 under various Shared Socioeconomic Pathway (SSP) scenarios: SSP126 (140.6 ± 9.3 PgC year-1), SSP245 (153.5 ± 13.4 PgC year-1), SSP370 (170.7 ± 16.9 PgC year-1), and SSP585 (194.1 ± 23.2 PgC year-1). These findings have important implications for understanding and projecting climate change, helping to develop more effective climate policies and carbon reduction strategies.

Development and evaluation of the interactive Model for Air Pollution and Land Ecosystems (iMAPLE) version 1.0

Abstract. Land ecosystems are important sources and sinks of atmospheric components. In turn, air pollutants affect the exchange rates of carbon and water fluxes between ecosystems and the atmosphere. However, these biogeochemical processes are usually not well presented in Earth system models, limiting the explorations of interactions between land ecosystems and air pollutants from regional to global scales. Here, we develop and validate the interactive Model for Air Pollution and Land Ecosystems (iMAPLE) by upgrading the Yale Interactive Terrestrial Biosphere Model with process-based water cycles, fire emissions, wetland methane (CH4) emissions, and trait-based ozone (O3) damage. Within iMAPLE, soil moisture and temperature are dynamically calculated based on the water and energy balance in soil layers. Fire emissions are dependent on dryness, lightning, population, and fuel load. Wetland CH4 is produced but consumed through oxidation, ebullition, diffusion, and plant-mediated transport. The trait-based scheme unifies O3 sensitivity of different plant functional types (PFTs) with the leaf mass per area. Validations show correlation coefficients (R) of 0.59–0.86 for gross primary productivity (GPP) and 0.57–0.84 for evapotranspiration (ET) across the six PFTs at 201 flux tower sites and yield an average R of 0.68 for CH4 emissions at 44 sites. Simulated soil moisture and temperature match reanalysis data with high R above 0.86 and low normalized mean biases (NMBs) within 7 %, leading to reasonable simulations of global GPP (R=0.92, NMB=1.3 %) and ET (R=0.93, NMB=-10.4 %) against satellite-based observations for 2001–2013. The model predicts an annual global area burned of 507.1 Mha, close to the observations of 475.4 Mha with a spatial R of 0.66 for 1997–2016. The wetland CH4 emissions are estimated to be 153.45 Tg [CH4] yr−1 during 2000–2014, close to the multi-model mean of 148 Tg [CH4] yr−1. The model also shows reasonable responses of GPP and ET to the changes in diffuse radiation and yields mean O3 damage of 2.9 % to global GPP. iMAPLE provides an advanced tool for studying the interactions between land ecosystems and air pollutants.

Impact of seasonal global land surface temperature (LST) change on gross primary production (GPP) in the early 21st century

Understanding the impact of global land surface temperature (LST) changes on gross primary production (GPP) is crucial for addressing global sustainability challenges effectively. This article explores the effects of winter and summer temperature variations on GPP from 2001 to 2020 on a global scale, employing a range of modeling techniques to investigate the complex relationship between temperature changes and GPP. The study employs various modeling approaches, including linear regression models, artificial neural networks, Random Forest (RF) models, and the XGBoost algorithm to examine both linear and non-linear relationships between global temperature changes and GPP. The results reveal that the RF and XGBoost models effectively capture the non-linear relationship between LST and GPP during both summer and winter, demonstrating a high level of statistical significance (P < 0.01) and achieving R2 values of 0.598. Furthermore, the study applies a Sen+MK trend analysis model to identify six distinct trend patterns in LST and GPP during both summer and winter seasons. In summer, the area exhibiting a non-significant increasing trend (NSI) for GPP covers 65,066,274.77 km2, whereas in winter, a strong significant decreasing trend (SSD) for GPP spans 64,537,108.31 km2. Notably, GPP patterns closely mirror those of LST, suggesting that rising high-temperature conditions during summer lead to reduced GPP, while increased low temperatures during winter promote GPP growth. Robust correlations between LST and GPP are observed under various trend conditions in both summer and winter, displaying strong statistical significance (P < 0.01), with SSD achieving a maximum R2 of 0.594. These findings contribute significantly to our comprehensive understanding of the dynamic relationship between LST and GPP and offer valuable insights for addressing the sustainable development of global climate change challenges.

Estimating the CO2 Fertilization Effect on Extratropical Forest Productivity From Flux‐Tower Observations

AbstractThe land sink of anthropogenic carbon emissions, a crucial component of mitigating climate change, is primarily attributed to the CO2 fertilization effect on global gross primary productivity (GPP). However, direct observational evidence of this effect remains scarce, hampered by challenges in disentangling the CO2 fertilization effect from other long‐term confounding drivers, particularly climatic changes. Here, we introduce a novel statistical approach to separate the CO2 fertilization effect on photosynthetic carbon uptake using eddy covariance (EC) records across 38 extratropical forest sites. We find the median stimulation rate of GPP to be 3.2 ± 0.9 gC m−2 yr−1 ppm−1 (or 16.4 ± 4.2% per 100 ppm) under increasing atmospheric CO2 across these sites, respectively. To validate the robustness of our findings, we test our statistical method using factorial simulations of an ensemble of process‐based land surface models. We address additional factors, including nitrogen deposition and land management, that may impact plant productivity, potentially confounding the attribution to the CO2 fertilization effect. Assuming these site‐specific effects offset to some extent across sites as random factors, the estimated median value still reflects the strength of the CO2 fertilization effect. However, disentanglement of these long‐term effects, often inseparable by timescale, requires further causal research. Our study provides direct evidence that the photosynthetic stimulation is maintained under long‐term CO2 fertilization across multiple EC sites. Such observation‐based quantification is key to constraining the long‐standing uncertainties in the land carbon cycle under rising CO2 concentrations.

Mutual inhibition effects of elevated CO2 and climate change on global forest GPP

The impact of CO2 fertilization on enhancing global forest gross primary productivity (GPP) is acknowledged, but its interaction with climate factors—air temperature (Tem), precipitation (Pre), vapor pressure deficit (VPD), and radiation (Rad)—remains unclear. In this study, global forest GPP trends from 1982 to 2018 were examined using BEPS, NIRv, FLUXCOM, and revised EC-LUE datasets, with interannual trends of 5.618 (p < 0.01), 5.831 (p < 0.01), 0.227, and 6.566 g C m−2 yr−1 (p < 0.01), respectively. Elevated CO2 was identified as the primary driver of GPP trends, with the dominant area ranging from 51.11% to 90.37% across different GPP datasets. In the NIRv and revised EC-LUE datasets, the positive impact of CO2 on GPP showed a decrease of 0.222 g C m−2 yr−1, while the negative impact of Rad increased by 0.007 g C m−2 yr−1. An inhibitory relationship was found between the actual effects of elevated CO2 and climate change on GPP in most forest types. At lower latitudes, Tem primarily constrained CO2 fertilization, while at higher latitudes, VPD emerged as the key limiting factor. This was mainly attributed to the potential trade-off or competition between elevated CO2 and climate change in influencing GPP, with strategic resource allocation varying across different forest ecosystems. This study highlights the significant inhibitory effects of elevated CO2 and climate change on global forest GPP, providing insights into the dynamic responses of forest ecosystems to changing environments.

Estimating global annual gross primary production based on satellite-derived phenology and maximal carbon uptake capacity

The high uncertainty regarding global gross primary production (GPP) remains unresolved. This study explored the relationships between phenology, physiology, and annual GPP to provide viable alternatives for accurate estimation. A statistical model of integrated phenology and physiology (SMIPP) was developed using GPP data from 145 FLUXNET sites to estimate the annual GPP for various vegetation types. By employing the SMIPP model driven by satellite-derived datasets of the global carbon uptake period (CUP) and maximal carbon uptake capacity (GPPmax), the global annual GPP was estimated for the period from 2001 to 2018. The results demonstrated that the SMIPP model accurately predicted annual GPP, with relative root mean square error values ranging from 11.20 to 19.29% for forest types and 20.49–35.71% for non-forest types. However, wetlands, shrublands, and evergreen forests exhibited relatively low accuracies. The average, trend, and interannual variation of global GPP during 2001–2018 were 132.6 Pg C yr−1, 0.25 Pg C yr−2, and 1.57 Pg C yr−1, respectively. They were within the ranges estimated in other global GPP products. Sensitivity analysis revealed that GPPmax had comparable effects to CUP in high-latitude regions but significantly greater impacts at the global scale, with sensitivity coefficients of 0.85 ± 0.23 for GPPmax and 0.46 ± 0.28 for CUP. This study provides a simple and practical method for estimating global annual GPP and highlights the influence of GPPmax and CUP on global-scale annual GPP.

Terrestrial Ecosystem Model in R (TEMIR) version 1.0: simulating ecophysiological responses of vegetation to atmospheric chemical and meteorological changes

Abstract. The newly developed offline land ecosystem model Terrestrial Ecosystem Model in R (TEMIR) version 1.0 is described here. This version of the model simulates plant ecophysiological (e.g., photosynthetic and stomatal) responses to varying meteorological conditions and concentrations of CO2 and ground-level ozone (O3) based on prescribed meteorological and atmospheric chemical inputs from various sources. Driven by the same meteorological data used in the GEOS-Chem chemical transport model, this allows asynchronously coupled experiments with GEOS-Chem simulations with unique coherency for investigating biosphere–atmosphere chemical interactions. TEMIR agrees well with FLUXNET site-level gross primary productivity (GPP) in terms of both the diurnal and monthly cycles (correlation coefficients R2&gt;0.85 and R2&gt;0.8, respectively) for most plant functional types (PFTs). Grass and shrub PFTs have larger biases due to generic model representations. The model performs best when driven by local site-level meteorology rather than reanalyzed gridded meteorology. Simulation using gridded meteorology agrees well for annual GPP in seasonality and spatial distribution with a global average of 134 Pg C yr−1. Application of Monin–Obukhov similarity theory to infer canopy conditions from gridded meteorology does not improve model performance, predicting an increase of +7 % in global GPP. Present-day O3 concentrations simulated by GEOS-Chem and an O3 damage scheme at high sensitivity show a 2 % reduction in global GPP with prominent reductions of up to 15 % in eastern China and the eastern USA. Regional correlations are generally unchanged when O3 is present and biases are reduced, especially for regions with high O3 damage. An increase in atmospheric CO2 concentration of 20 ppmv from the level in 2000 to the level in 2010 modestly decreases O3 damage due to reduced stomatal uptake, consistent with ecophysiological understanding. Our work showcases the utility of this version of TEMIR for evaluating biogeophysical responses of vegetation to changes in atmospheric composition and meteorological conditions.

Progress and challenges in remotely sensed terrestrial carbon fluxes

ABSTRACT Accurate evaluation of terrestrial carbon balance is essential for designing climate change mitigation policies, and capabilities of remote sensing techniques in monitoring carbon fluxes are widely recognized for their great contributions to regional and global carbon budget accounting. In this review, we synthesized satellite-based data and methodologies to estimate the main flux components of terrestrial carbon balance and their uncertainties over the past two decades. The global gross primary production (GPP) during the period 2001–2022 is 134 ± 14 PgC yr−1, and nearly half of them occurs in tropical forest regions such as South America and Africa. Less than 2% of global GPP is converted into a net carbon sink of 2.28 ± 1.12 PgC yr−1 using satellite-based atmospheric inversion during 2015–2020, and this sink is comparable to the stock change-based estimate (2.49 PgC yr−1) but twice as large as model-based estimate (1.08 ± 0.78 PgC yr−1). By decomposing satellite-derived net carbon balance into different terms including satellite-derived carbon emissions from land-use change and wildfires (3.55 PgC yr−1), we inferred that ~ 43% of global GPP would be respired through soil microbes (57.1 PgC yr−1), but which is higher than the previous bottom-up estimate (39–46 PgC yr−1). We then propose that an accurate remote sensing of terrestrial carbon balance requires to enhance representations of photosynthetic responses to rising CO2 and disturbances, develop satellite-constrained belowground carbon dynamics and separate natural fluxes from anthropogenic CO2 emissions, by integrating multi-source satellite sensors in orbit, revolutionized remote sensing capabilities with focused field campaigns in data-scarce regions.

Global datasets of hourly carbon and water fluxes simulated using a satellite-based process model with dynamic parameterizations

Abstract. Diagnostic terrestrial biosphere models (TBMs) forced by remote sensing observations have been a principal tool for providing benchmarks on global gross primary productivity (GPP) and evapotranspiration (ET). However, these models often estimate GPP and ET at coarse daily or monthly steps, hindering analysis of ecosystem dynamics at the diurnal (hourly) scales, and prescribe some essential parameters (i.e., the Ball–Berry slope (m) and the maximum carboxylation rate at 25 °C (Vcmax25)) as constant, inducing uncertainties in the estimates of GPP and ET. In this study, we present hourly estimations of global GPP and ET datasets at a 0.25° resolution from 2001 to 2020 simulated with a widely used diagnostic TBM – the Biosphere–atmosphere Exchange Process Simulator (BEPS). We employed eddy covariance observations and machine learning approaches to derive and upscale the seasonally varied m and Vcmax25 for carbon and water fluxes. The estimated hourly GPP and ET are validated against flux observations, remote sensing, and machine learning-based estimates across multiple spatial and temporal scales. The correlation coefficients (R2) and slopes between hourly tower-measured and modeled fluxes are R2=0.83, regression slope =0.92 for GPP, and R2=0.72, regression slope =1.04 for ET. At the global scale, we estimated a global mean GPP of 137.78±3.22 Pg C yr−1 (mean ± 1 SD) with a positive trend of 0.53 Pg C yr−2 (p&lt;0.001), and an ET of 89.03±0.82×103 km3 yr−1 with a slight positive trend of 0.10×103 km3 yr−2 (p&lt;0.001) from 2001 to 2020. The spatial pattern of our estimates agrees well with other products, with R2=0.77–0.85 and R2=0.74–0.90 for GPP and ET, respectively. Overall, this new global hourly dataset serves as a “handshake” among process-based models, remote sensing, and the eddy covariance flux network, providing a reliable long-term estimate of global GPP and ET with diurnal patterns and facilitating studies related to ecosystem functional properties, global carbon, and water cycles. The hourly GPP and ET estimates are available at https://doi.org/10.57760/sciencedb.ecodb.00163 (Leng et al., 2023a) and the accumulated daily datasets are available at https://doi.org/10.57760/sciencedb.ecodb.00165 (Leng et al., 2023b).

Comparative assessment of leaf photosynthetic capacity datasets for estimating terrestrial gross primary productivity

The maximum Rubisco carboxylation rate normalized to 25 °C (Vcmax25) is a key parameter in terrestrial biosphere models for simulating carbon cycling. Recently, global distributions of Vcmax25 have been derived through various methods and different data, including field measurements, ecological optimality theory (EOT), leaf chlorophyll content (LCC), and solar-induced chlorophyll fluorescence (SIF). However, direct validation poses challenges due to high uncertainty arising from limited ground-based observations. This study conducted an indirect evaluation of four Vcmax25 datasets by assessing the accuracy of gross primary productivity (GPP) simulated using the Biosphere-atmosphere Exchange Process Simulator (BEPS) at both site and global scales. Results indicate that, compared to utilizing Vcmax25 fixed by plant functional types (PFT) derived from field measurements, incorporating Vcmax25 derived from SIF and LCC (SIF + LCC), or solely LCC, into BEPS significantly reduces simulated errors in the annual total GPP, with a 23.2 %–25.1 % decrease in the average absolute bias across 196 FLUXNET2015 sites. Daily GPP for evergreen needleleaf forests, deciduous broadleaf forests, shrublands, grasslands, and croplands shows a 7.8 %–27.6 % decrease in absolute bias, primarily attributed to reduced simulation errors during off-peak seasons of vegetation growth. Conversely, the annual total GPP error simulated using EOT-derived Vcmax25 increases slightly (2.2 %) compared to that simulated using PFT-fixed Vcmax25. This is primarily due to a significant overestimation in evergreen broadleaf forests and underestimation in croplands, despite slight increased accuracy for other PFTs. The global annual GPP simulated using Vcmax25 with seasonal variations (i.e., LCC Vcmax25 and SIF + LCC Vcmax25) yields a 4.3 %–7.3 % decrease compared to that simulated using PFT-fixed Vcmax25. Compared to FLUXCOM and GOSIF GPP products, the GPP simulated based on SIF + LCC Vcmax25 and LCC Vcmax25 demonstrates better consistency (R2 = 0.91–0.93, RMSE = 314.2–376.6 g C m−2 yr−1). This study underscores the importance of accurately characterizing the spatiotemporal variations in Vcmax25 for the accurate simulation of global vegetation productivity.

A novel hybrid modelling framework for GPP estimation: Integrating a multispectral surface reflectance based [formula omitted] simulator into the process-based model

Terrestrial gross primary productivity (GPP) is the key element in the carbon cycle process. Accurate GPP estimation hinges on the maximum carboxylation rate (Vcmax,025). The high uncertainty in deriving ecosystem-level Vcmax,025 has long hampered efforts toward the performance of the GPP model. Recently studies suggest the strong relationship between spectral reflectance and Vcmax,025. We proposed the multispectral surface reflectance-driven Vcmax,025 simulator using the fully connected deep neural network and built the hybrid modelling framework for GPP estimation by integrating the data-driven Vcmax,025 simulator in the process-based model. The performance of hybrid GPP model was evaluated at 95 flux sites. The result shows that the multispectral surface reflectance-driven Vcmax,025 simulator acquires the satisfactory estimation, with correlation coefficient (R), root mean square error (RMSE) and median absolute percentage error (MdAPE) ranging from 0.34 to 0.80, 14 to 43 μmol m−2 s−1 and 21 % to 66 % across different land cover types, respectively. The hybrid framework generates good GPP estimates with R, RMSE and MdAPE varying from 0.76 to 0.89, 1.79 to 6.16 μmol m−2 s−1 and 27 % to 90 %, respectively. Compared with EVI-driven method, the multispectral surface reflectance significantly improves the Vcmax,025 and GPP estimates, with MdAPE declining by 0.6 %–18 % and 1 % to 21 %, respectively. The Shapley value analysis reveals that red (620–670 nm), near-infrared (841–876 nm) and shortwave infrared (1628–1652 nm and 2105–2155 nm) are the key bands for Vcmax,025 estimation. This study highlights the potential of multispectral surface reflectance for quantifying ecosystem-level Vcmax,025. The new hybrid framework fully extracts the information of all available spectral bands using deep learning to reduce parameter uncertainty while maintains the description of photosynthetic process to ensure its physical reasonability. It can serve as a powerful tool for accurate global GPP estimation.

Divergent trajectories of future global gross primary productivity and evapotranspiration of terrestrial vegetation in Shared Socioeconomic Pathways

Understanding the future trends of carbon and water fluxes between terrestrial ecosystems and the atmosphere is crucial for predicting Earth's climate dynamics. This study employs an advanced numerical approach to project global gross primary productivity (GPP) and evapotranspiration (ET) from 2001 to 2100 under various climate scenarios based on Shared Socioeconomic Pathways (SSPs). To improve predictions of vegetation dynamics, we introduce a novel model (CoLM-PVPM), an enhancement of the Common Land Model version 2014 (CoLM2014), incorporating a prognostic vegetation phenology model (PVPM). Compared to CoLM2014 that relies on satellite-based leaf area index (LAI) inputs, CoLM-PVPM predicts LAI time series using climate variables. Model validation using historical data from 2001 to 2010 demonstrates PVPM in capturing spatiotemporal variations in satellite LAI. Our modeling results indicate that annual averaged LAI and total GPP increase under SSP1–2.6 but decrease under SSP2–4.5, SSP3–7.0, and SSP5–8.5 by 2100. By comparison, annual total ET consistently increases under all SSP scenarios by 2100. Global annual averaged LAI is highly correlated with annual total GPP in all scenarios, while its correlation with annual total ET weakens in SSP2–4.5, SSP3–7.0, and SSP5–8.5. Global annual total vapor pressure deficit (VPD) and precipitation are highly correlated with annual total ET in all scenarios. As emission levels increase, the negative correlation between annual total VPD and GPP strengthens, while the correlation between annual total precipitation and GPP weakens. This research presents an improved model for predicting terrestrial vegetation processes and underscores the importance of low carbon emission scenarios in maintaining carbon-water balances in specific regions.

Measuring stomatal and guard cell metrics for plant physiology and growth using StoManager1.

Automated guard cell detection and measurement are vital for understanding plant physiological performance and ecological functioning in global water and carbon cycles. Most current methods for measuring guard cells and stomata are laborious, time-consuming, prone to bias, and limited in scale. We developed StoManager1, a high-throughput tool utilizing geometrical, mathematical algorithms, and convolutional neural networks to automatically detect, count, and measure over 30 guard cell and stomatal metrics, including guard cell and stomatal area, length, width, stomatal aperture area/guard cell area, orientation, stomatal evenness, divergence, and aggregation index. Combined with leaf functional traits, some of these StoManager1-measured guard cell and stomatal metrics explained 90% and 82% of tree biomass and intrinsic water use efficiency (iWUE) variances in hardwoods, making them substantial factors in leaf physiology and tree growth. StoManager1 demonstrated exceptional precision and recall (mAP@0.5 over 0.96), effectively capturing diverse stomatal properties across over 100 species. StoManager1 facilitates the automation of measuring leaf stomatal and guard cells, enabling broader exploration of stomatal control in plant growth and adaptation to environmental stress and climate change. This has implications for global gross primary productivity (GPP) modeling and estimation, as integrating stomatal metrics can enhance predictions of plant growth and resource usage worldwide. Easily accessible open-source code and standalone Windows executable applications are available on a GitHub repository (https://github.com/JiaxinWang123/StoManager1) and Zenodo (https://doi.org/10.5281/zenodo.7686022).

Gross primary production of Mediterranean watersheds: Using isotope mass balance approach to improve estimations

AbstractGlobal‐scale estimates of carbon fluxes from satellite data‐driven models are constrained by considerable uncertainties regarding Gross Primary Production (GPP) and the lack of the watershed‐scale measurements required for model calibration. Recently conducted global modelling efforts indicate that semiarid ecosystems dominate the increasing trends and inter‐annual variation of net CO2 exchange with the atmosphere, but semi‐arid regions have received little attention with regard to GPP estimation. In this study, we used the distinct isotope effect of transpiration and evaporation to calculate transpiration losses and subsequently CO2 uptake by terrestrial vegetation through the water and carbon cycle using the water use efficiency of plants. By studying two Mediterranean watersheds with contrasted environmental conditions over several hydrological years, we found a strong dependence of GPP on annual and seasonal water availability. The results demonstrated that when compared to GPP values obtained in worldwide biomes using biological methods, our isotope approach was validated, highlighting the limitations of satellite‐data‐driven models like MODIS in capturing the impact of water stress on photosynthesis and GPP estimates. These results encourage investigation of GPP by the isotope mass balance approach where direct carbon flux measurements are rare or absent in order to help to substantiate, modify or shed doubt on interpolated GPP for those regions and achieve consensus on global GPP estimates. Given the relevant role of semi‐arid ecosystems in the global carbon balance as well as the limitation of existing data sets, our improved method based on the isotope mass balance approach helps to obtain rapid and affordable estimates of GPP for semi‐arid ecosystems.