Gross Primary Production Measurements Research Articles

Using automated machine learning for the upscaling of gross primary productivity

Abstract. Estimating gross primary productivity (GPP) over space and time is fundamental for understanding the response of the terrestrial biosphere to climate change. Eddy covariance flux towers provide in situ estimates of GPP at the ecosystem scale, but their sparse geographical distribution limits larger-scale inference. Machine learning (ML) techniques have been used to address this problem by extrapolating local GPP measurements over space using satellite remote sensing data. However, the accuracy of the regression model can be affected by uncertainties introduced by model selection, parameterization, and choice of explanatory features, among others. Recent advances in automated ML (AutoML) provide a novel automated way to select and synthesize different ML models. In this work, we explore the potential of AutoML by training three major AutoML frameworks on eddy covariance measurements of GPP at 243 globally distributed sites. We compared their ability to predict GPP and its spatial and temporal variability based on different sets of remote sensing explanatory variables. Explanatory variables from only Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data and photosynthetically active radiation explained over 70 % of the monthly variability in GPP, while satellite-derived proxies for canopy structure, photosynthetic activity, environmental stressors, and meteorological variables from reanalysis (ERA5-Land) further improved the frameworks' predictive ability. We found that the AutoML framework Auto-sklearn consistently outperformed other AutoML frameworks as well as a classical random forest regressor in predicting GPP but with small performance differences, reaching an r2 of up to 0.75. We deployed the best-performing framework to generate global wall-to-wall maps highlighting GPP patterns in good agreement with satellite-derived reference data. This research benchmarks the application of AutoML in GPP estimation and assesses its potential and limitations in quantifying global photosynthetic activity.

Using SAR-based products to calculate potato carbon uptake in a tropical Andean region

ABSTRACT Gross primary productivity (GPP) is an essential parameter to estimate the efficiency of carbon transfer in terrestrial ecosystems. The daily GPP has been monitored in the past using mainly optical satellite imagery. However, GPP has never been monitored using satellite Synthetic Aperture Radar (SAR) imagery. We evaluate the possibility of using Sentinel-1 SAR data to estimate GPP and compare with field GPP measurements using eddy covariance (EC) systems in three commercial potato fields, located in the Andean region of Colombia, using three different water irrigation levels. Raw and processed radar data from Sentinel-1 were compared with the daily and accumulated GPP from the EC. Results indicate that SAR data has a high correlation with the accumulated GPP measured by the EC in the field. The normalized radar backscatter and radar brilliance coefficients with VH polarization show a good correlation with the EC accumulated GPP in irrigated potato fields (R2: 0.77–0.81), while radar vegetation indices show a good correlation with the EC accumulated GPP in potato fields with no irrigation (R2≈0.82). In particular, the accumulated GPP during the whole potato crop cycle was estimated with good accuracy (~5% error with irrigation and ~10% error with no irrigation).

A novel approach for retrieving GPP of evergreen forest regions of India using random forest regression

Gross Primary Productivity (GPP) is a crucial variable of global carbon cycle for determining the ecosystem’s health. Various methods are devised to quantify GPP and upscale it in both time and space. The most common methods are physical model and eddy covariance-based estimation, which are very restricted to surrounding area of study, only. The alternative methods are empirical (e.g., LUE, CASA, SCARF, and MODIS) and Machine learning (ML) models that employ remote sensing satellite data and geographical factors. However, for using ML models, ground-based measurements of GPP is a very important factor, which is not available in most places. We propose an alternative and effective way of estimating the GPP using the ML model and data from various flux sites around the globe for a particular plant functional type (PFT). In the present study, RF is used as ML model, which is trained on global GPP data from evergreen forest and implemented in Indian region. The key findings indicated that ML-based GPP is highly accurate and hence, we generated 20 years of time series GPP dataset (2001–2020). We validated with ground-based flux tower observations during 2016–2018 for three sites (very limited datasets) in India and compared them with MODIS GPP. The coefficient of determination (R2) value of the ML-based model was 0.84 with root mean square error (RMSE) of 1.45 gC m−2 Day−1 and mean absolute error (MAE) of 0.838 gC m−2 Day−1. The proposed approach is highly accurate and far better than the MODIS-based GPP. Therefore, it can be further extended to other forest types for a holistic assessment of the carbon cycle of a region.

Satellite-based assessment of water use and leaf area efficiencies of dryland conifer forests along an aridity gradient

Dryland forests worldwide are increasingly threatened by drought stress due to climate change. Understanding the relationships between forest structure and function is essential for managing dryland forests to adapt to these changes. We investigated the structure-function relationships in four dryland conifer forests distributed along a semiarid to subhumid climatic aridity gradient. Forest structure was represented by leaf area index (LAI) and function by gross primary productivity (GPP), evapotranspiration (ET), and the derived efficiencies of water use (WUE = GPP/ET) and leaf area (LAE = GPP/LAI). Estimates of GPP and ET were based on the observed relationships between high-resolution vegetation indices from VENμS and Sentinel-2A satellites and flux data from three eddy covariance towers in the study regions between November 2015 to October 2018. The red-edge-based MERIS Terrestrial Chlorophyll Index (MTCI) from VENμS and Sentinel-2A showed strong correlations to flux tower GPP and ET measurements for the three sites (R2cal > 0.91, R2val > 0.84). Using our approach, we showed that as LAI decreased with decreasing aridity index (AI) (i.e., dryer conditions), estimated GPP and ET decreased (R2 > 0.8 to LAI), while WUE (R2 = 0.68 to LAI) and LAE increased. The observed global-scale patterns are associated with a variety of forest vegetation characteristics, at the local scale, such as tree species composition and density. However, our results point towards a canopy-level mechanism, where the ecosystem-LAI and resultant proportion of sun-exposed vs. shaded leaves are primary determinants of WUE and LAE along the studied climatic aridity gradient. This work demonstrates the importance of high-resolution (spatially and spectrally) remote sensing data conjugated with flux tower data for monitoring dryland forests and understanding the intricate structure-function interactions in their response to drying conditions.

Accelerated vegetative growth measured by gross primary productivity in China from 1980 to 2018

Vegetation response to climate change can be characterized as long-term trends, but with periodic oscillations. Revealing vegetation response trends, with oscillations, is essential for formulating climate change mitigation and adaptation planning, as well as adopting ecological protection and restoration measures. Here, the long-term trend and periodic oscillations in vegetation changes were explored through a measure of gross primary production (GPP), an indicator of vegetation photosynthesis, using data of both high spatial resolution (1 km2) and a long-term of nearly 4 decades, from 1980 to 2018. The GPP data were estimated for China’s terrestrial ecosystem and validated against the eddy covariance-based observations of 70 site-year in China and compared to other GPP data products. The results showed the estimated GPP can explain 76% of spatial–temporal variance with the same data accuracy as other products, but with the advantage of both high spatial resolution and long time periods that other products do not have. From the estimated GPP data, Chinese terrestrial vegetation growth was found to have a periodic oscillation of close to 2.79 years and a long-term increasing trend of 59.8 Mt C per year over 92% of all vegetated land in China from 1980 to 2018. Moreover, it was found that vegetation growth accelerated significantly in the two recent decades over the previous two decades. The growth rate was 83.5 Mt C a−1 from 2000 to 2018, which was 2.31 times faster than the rate from 1980 to 2000. The underlying mechanism of the vegetation responses was analyzed by ridge regression, a method to decrease the effects of multicollinearity. Long-term variability and trends in vegetation could be attributed to a warming minimum air temperature when examined with precipitation, maximum air temperature, and solar shortwave radiation. Though more influences, such as those from freezing and thawing and human activities, should be explored in the future, this study offers insight that warming, with increasing daily minimum temperatures, is likely to be a vital contributor to ecosystem changes, and will be essential knowledge for mitigation of effects and adaptation to future climate.

Remote sensing of daily evapotranspiration and gross primary productivity of four forest ecosystems in East Asia using satellite multi-channel passive microwave measurements

Satellite multi-channel passive microwave observations provide valuable information of canopy vertical structure and vegetation water content, while their potential in the simultaneous estimation of forest evapotranspiration (ET) and gross primary productivity (GPP) has not been explored. This study combined satellite multi-channel microwave Emissivity Difference Vegetation Indices (EDVIs) from 6.9 to 36.5 GHz with the Penman-Monteith (PM), Priestley-Taylor (PT) and light use efficiency (LUE) models to estimate daily ET and GPP over four forest ecosystems in East Asia. In-situ measurements of daily ET, GPP and meteorological variables from ChinaFLUX network (2003 to 2010) were used to calculate and validate the method. An exponential relationship was found between EDVIs and canopy conductance across the forests. By inverting EDVIs-based LUE model coupled with PM and PT models when canopy carbon uptake and water content reached a maximum in growing seasons, the maximum LUE (EDVI-LUEmax) was derived and found to agree with other independent studies. Site validations showed that the EDVI-based method performed a good simultaneous estimation of daily ET and GPP at a tropical forest and three sub-tropical and temperate forests (R2 of 0.64 to 0.94 and RMSE of −11.7% to 19.6% for ET, R2 of 0.39 to 0.73 and RMSE of 31.2% to 50.4% for GPP). The modeled relationship between GPP and ET also agreed that from in-situ measurements. The results were also found to correlate with microwave frequencies. In-situ ET and GPP showed a stronger correlation with EDVI at 18.7 and 36.5 GHz than that at 10.7 and 18.7 GHz across the forests. EDVI-LUEmax of a temperate mixed forest increased as EDVI incorporated lower frequency microwave emissions from deeper canopy layers, which reduced the negative bias in GPP estimation. Overall, this study demonstrated the potential of multi-channel microwave measurements in estimating ET and GPP with the potential effects of vegetation vertical structure.

The magnitude and pace of photosynthetic recovery after wildfire in California ecosystems

Wildfire modifies the short- and long-term exchange of carbon between terrestrial ecosystems and the atmosphere, with impacts on ecosystem services such as carbon uptake. Dry western US forests historically experienced low-intensity, frequent fires, with patches across the landscape occupying different points in the fire-recovery trajectory. Contemporary perturbations, such as recent severe fires in California, could shift the historic stand-age distribution and impact the legacy of carbon uptake on the landscape. Here, we combine flux measurements of gross primary production (GPP) and chronosequence analysis using satellite remote sensing to investigate how the last century of fires in California impacted the dynamics of ecosystem carbon uptake on the fire-affected landscape. A GPP recovery trajectory curve of more than five thousand fires in forest ecosystems since 1919 indicated that fire reduced GPP by [Formula: see text] g C m[Formula: see text] y[Formula: see text]([Formula: see text]) in the first year after fire, with average recovery to prefire conditions after [Formula: see text] y. The largest fires in forested ecosystems reduced GPP by [Formula: see text] g C m[Formula: see text] y[Formula: see text] (n = 401) and took more than two decades to recover. Recent increases in fire severity and recovery time have led to nearly [Formula: see text] MMT CO[Formula: see text] (3-y rolling mean) in cumulative forgone carbon uptake due to the legacy of fires on the landscape, complicating the challenge of maintaining California's natural and working lands as a net carbon sink. Understanding these changes is paramount to weighing the costs and benefits associated with fuels management and ecosystem management for climate change mitigation.

Capturing the spatiotemporal variations in the gross primary productivity in coastal wetlands by integrating eddy covariance, Landsat, and MODIS satellite data: A case study in the Yangtze Estuary, China

Accurate monitoring of the spatiotemporal variations in the gross primary productivity (GPP) in coastal wetlands is essential for blue carbon quantification. However, the currently available moderate-spatial-resolution GPP algorithms and products may contain large uncertainties and may rarely meet the demands of tidal wetland research. In this paper, we examine two statistical methods, namely, simple linear regression and random forest regression, to capture the high spatial resolution GPP variations in the salt marshes in the Yangtze Estuary by integrating multi-source data, including eddy covariance, Landsat, and Moderate-resolution Imaging Spectroradiometer (MODIS) satellite data. The satellite-derived photosynthetically active radiation (PAR), enhanced vegetation index (EVI), normalized difference vegetation index (NDVI), and normalized difference water index (NDWI) are used individually and in different combinations to drive the two statistical methods and investigate their performances in estimating the GPP. The results show that together with the PAR, the EVI generally has the greatest potential for GPP estimation (R2 > 0.75, RMSE < 6.80 μmol m−2 s−1) based on simple linear regression. However, the NDVI outperforms the EVI in wetlands with little tidal flooding and plant litter. The random forest method improves the performance in terms of upscaling the GPP measurements to large regions, and the combination of the EVI, NDVI, NDWI, and PAR performs the best, indicating that the use of the more advanced algorithm and the incorporation of a remotely sensed index that reflects the tidal influence are beneficial to GPP estimation in coastal salt marshes. Moreover, using Landsat data with a high spatial resolution (30 m) yields a much more accurate GPP estimation than using MODIS data with a 500 m resolution. Our results demonstrate that the use of high spatial resolution data, the proper use of remotely sensed indices, and the incorporation of a good combination of these indices through an advanced algorithm such as a machine learning algorithm are vital for capturing the nearly instantaneous GPP variations in coastal wetlands at large spatial scales.

Observations of Satellite Land Surface Phenology Indicate That Maximum Leaf Greenness Is More Associated With Global Vegetation Productivity Than Growing Season Length

AbstractVegetation green leaf phenology directly impacts gross primary productivity (GPP) of terrestrial ecosystems. Satellite observations of land surface phenology (LSP) provide an important means to monitor the key timing of vegetation green leaf development. However, differences between satellite‐derived LSP proxies and in situ measurements of GPP make it difficult to quantify the impact of climate‐induced changes in green leaf phenology on annual GPP. Here, we used 1,110 site‐years of GPP measurements from eddy‐covariance towers in association with time series of satellite LSP observations from 2000 to 2014 to show that while satellite LSP explains a large proportion of variation in annual GPP, changes in green‐leaf‐based growing season length (GSL, leaf development period from spring to autumn) had less impact on annual GPP by ∼30% than GSL changes in GPP‐based photosynthetic duration. Further, maximum leaf greenness explained substantially more variance in annual GPP than green leaf GSL, highlighting the role of future vegetation greening trends on large‐scale carbon budgets. Site‐level variability contributes a substantial proportion of annual GPP variance in the model based on LSP metrics, suggesting the importance of local environmental factors altering regional GPP. We conclude that satellite LSP‐based inferences regarding large‐scale dynamics in GPP need to consider changes in both green leaf GSL and maximum greenness.

Land surface phenology indicators retrieved across diverse ecosystems using a modified threshold algorithm

Land surface phenology (LSP), the study of the seasonal vegetation dynamics from remote sensing imagery, provides crucial information for plant monitoring and reflects the responses of ecosystems to climate change. The Moderate Resolution Imaging Spectroradiometer (MODIS) phenology product (MCD12Q2) provides global LSP information, but it has large spatial gaps in many regions, especially in ecosystems where rainfall influences phenology more than temperature. This study aimed to improve spatial coverage of LSP retrieval in these ecosystems. To do so, we used a regionally modified threshold algorithm for LSP retrievals, which were tested over continental Australia as it includes diverse landscapes of arid, mesic, and forest environments. We generated LSP metrics annually from 2003 to 2018 using satellite Enhanced Vegetation Index (EVI) time series at 500 m resolution, including the start, peak, end, and length of growing seasons, the minimum EVI value prior to and after the peak date, the seasonal maximum EVI value, the integral EVI value during the growing season (an approximation of productivity), and seasonal amplitude (maximum EVI value minus minimum EVI). Our regionally optimised algorithm improved the spatial coverage of LSP information in Australia from only 26 % of the continent to 70 % averaged across 16 years. Our results showed that the growing season amplitude was low (EVI < 0.1) over arid/semi-arid shrublands and savannas, tropical and subtropical savannas, and temperate evergreen forests, whose LSP metrics were captured by our regional algorithm and not by the global product. Some ecosystems, such as arid/semi-arid shrublands and savannas, showed more irregular phenology with low seasonal dynamics, and the growing seasons could skip a year or occur more than once in a year depending on climate conditions. Our algorithm was more sensitive to ecosystems with low seasonal amplitudes. We found that the detectability of LSP increases as the growing season amplitude increases, regardless of vegetation cover. Evaluation of the LSP metrics using eddy covariance flux tower measurements of gross primary productivity (GPP) demonstrated the reliability and accuracy of the algorithm. These improved LSP retrievals provide a greater understanding of the vegetation phenology across diverse ecosystems, especially savanna, shrubland, and evergreen forest ecosystems that cover more than 30 % of the land globally. The LSP provides essential information for ecological and agricultural studies such as quantifying bushfire fuel accumulation and forest carbon cycling, whilst enhancing our capacity for quantifying ecological responses to climate change.

Assessing the Use of Sentinel-2 Data for Spatio-Temporal Upscaling of Flux Tower Gross Primary Productivity Measurements

The conservation, restoration and sustainable use of wetlands is the target of several international agreements, among which are the Sustainable Development Goals (SDGs). Earth Observation (EO) technologies can assist national authorities in monitoring activities and the environmental status of wetlands to achieve these targets. In this study, we assess the capabilities of the Sentinel-2 instrument to model Gross Primary Productivity (GPP) as a proxy for the monitoring of ecosystem health. To estimate the spatial and temporal variation of GPP, we develop an empirical model correlating in situ measurements of GPP, eight Sentinel-2 derived vegetation indexes (VIs), and different environmental drivers of GPP. The model automatically performs an interdependency analysis and selects the model with the highest accuracy and statistical significance. Additionally, the model is upscaled across larger areas and monthly maps of GPP are produced. The study methodology is applied in a marsh ecosystem located in Doñana National Park, Spain. In this application, a combination of the red-edge chlorophyll index (CLr) and rainfall data results in the highest correlation with in situ measurements of GPP and is used for the model formulation. This yields a coefficient of determination (R2) of 0.93, Mean Absolute Error (MAE) equal to 0.52 gC m−2 day−1, Root Mean Squared Error (RMSE) equal to 0.63 gC m−2 day−1, and significance level p &lt; 0.05. The model outputs are compared with the MODIS GPP global product (MOD17) for reference; an enhancement of the estimation of GPP is found in the applied methodology.

A Low-Cost Method of Measuring the In Situ Primary Productivity of Periphyton Communities of Lentic Waters.

Measuring the in situ primary productivity of periphyton during the growing season gradient can elucidate the quantitative effect of environmental drivers (mainly phosphorus concentration and light intensity) and species composition on primary productivity. Primary productivity is mainly driven by light intensity, temperature, availability of nutrients, and distribution of the ionic species of the carbonate system in the respective depths of the euphotic zone. It is a complex system that is very difficult to simulate in the laboratory. This cheap, transportable, and easy-to-build floating barge allows measuring the primary productivity accurately-directly under the actual natural conditions. The methodology is based on measuring the primary productivity in real time using noninvasive oxygen sensors integrated into tightly sealed glass jars, enabling online oxygen flux monitoring and providing new insights into metabolic activities. Detailed seasonal in situ measurements of gross primary productivity of microbial mats (or other benthic organisms) can improve current knowledge of the processes controlling primary productivity dynamics in lentic waters.

How well do recently reconstructed solar-induced fluorescence datasets model gross primary productivity?

The collection of various long-term reconstructed solar-induced fluorescence (SIF) datasets derived at a range of spatio-temporal scales provides new opportunities for modelling vegetation dynamics, in particular, gross primary productivity (GPP). Information about the proximity of the reconstructed SIF (SIFr) datasets to GPP across land cover types and climatic conditions provides important support for a better application of these products for modelling applications. We conducted a multiscale analysis of four different long-term (12 years, 2007–2018) high-resolution global SIFr datasets (0.05° × 0.05°), namely – CSIF (Contiguous SIF), GOSIF (Global OCO-2 SIF), LUE-SIF (Light Use Efficiency SIF), and HSIF (Harmonized SIF) - at 4-day, 8-day, and monthly time scales and found that for the majority of sites, the SIFr is linearly related to ground-based GPP measurements with the eddy covariance method. While the relationship between SIFr and GPP (i.e., the slope - GPP/SIFr) varied significantly across the SIFr datasets, sites, and land cover types, all four SIFr datasets were unequivocally a better predictor of GPP compared to remotely sensed vegetation indices – NDVI (normalized difference vegetation index) and EVI (enhanced vegetation index), sensed by the MODIS satellite. Furthermore, we also analyzed SIF-GPP relationships during drought vs non-drought conditions and found that for about 30% of the sites, comprising mostly non-forests site, the SIF-GPP relationship became weaker (decreased R2) with a lower slope during drought conditions compared to non-drought conditions. Among the four different products, the CSIF (at 4-day timescale) and GOSIF (at 8-day timescale) predicted GPP better compared to LUE-SIF and HSIF across all land cover types. Owing to their long-term availability (since 2000 for CSIF and GOSIF), these SIFr datasets combined with proxies of ecosystem properties can be used to appropriately capture vegetation dynamics and the interannual variabilities across a wide range of climatic conditions.

Developing a parsimonious canopy model (PCM v1.0) to predict forest gross primary productivity and leaf area index of deciduous broad-leaved forest

Abstract. Temperate forest ecosystems play a crucial role in governing global carbon and water cycles. However, unprecedented global warming presents fundamental alterations to the ecological functions (e.g., carbon uptake) and biophysical variables (e.g., leaf area index) of forests. The quantification of forest carbon uptake, gross primary productivity (GPP), as the largest carbon flux has a direct consequence on carbon budget estimations. Part of this assimilated carbon stored in leaf biomass is related to the leaf area index (LAI), which is closely linked to and is of critical significance in the water cycle. There already exist a number of models to simulate dynamics of LAI and GPP; however, the level of complexity, demanding data, and poorly known parameters often prohibit the model applicability over data-sparse and large domains. In addition, the complex mechanisms associated with coupling the terrestrial carbon and water cycles poses a major challenge for integrated assessments of interlinked processes (e.g., accounting for the temporal dynamics of LAI for improving water balance estimations and soil moisture availability for enhancing carbon balance estimations). In this study, we propose a parsimonious forest canopy model (PCM) to predict the daily dynamics of LAI and GPP with few required inputs, which would also be suitable for integration into state-of-the-art hydrologic models. The light use efficiency (LUE) concept, coupled with a phenology submodel, is central to PCM (v1.0). PCM estimates total assimilated carbon based on the efficiency of the conversion of absorbed photosynthetically active radiation into biomass. Equipped with the coupled phenology submodel, the total assimilated carbon partly converts to leaf biomass, from which prognostic and temperature-driven LAI is simulated. The model combines modules for the estimation of soil hydraulic parameters based on pedotransfer functions and vertically weighted soil moisture, considering the underground root distribution, when soil moisture data are available. We test the model on deciduous broad-leaved forest sites in Europe and North America, as selected from the FLUXNET network. We analyze the model's parameter sensitivity on the resulting GPP and LAI and identified, on average, 10 common sensitive parameters at each study site (e.g., LUE and SLA). The model's performance is evaluated in a validation period, using in situ measurements of GPP and LAI (when available) at eddy covariance flux towers. The model adequately captures the daily dynamics of observed GPP and LAI at each study site (Kling–Gupta efficiency, KGE, varies between 0.79 and 0.92). Finally, we investigate the cross-location transferability of model parameters and derive a compromise parameter set to be used across different sites. The model also showed robustness with the compromise single set of parameters, applicable to different sites, with an acceptable loss in model skill (on average ±8 %). Overall, in addition to the satisfactory performance of the PCM as a stand-alone canopy model, the parsimonious and modular structure of the developed PCM allows for a smooth incorporation of carbon modules to existing hydrologic models, thereby facilitating the seamless representation of coupled water and carbon cycle components, i.e., prognostic simulated vegetation leaf area index (LAI) would improve the representation of the water cycle components (i.e., evapotranspiration), while GPP predictions would benefit from the simulated soil water storage from a hydrologic model.

Attributing differences of solar-induced chlorophyll fluorescence (SIF)-gross primary production (GPP) relationships between two C4 crops: corn and miscanthus

There remains limited information to characterize the solar-induced chlorophyll fluorescence (SIF)-gross primary production (GPP) relationship in C4 cropping systems. The annual C4 crop corn and perennial C4 crop miscanthus differ in phenology, canopy structure and leaf physiology. Investigating the SIF-GPP relationships in these species could deepen our understanding of SIF-GPP relationships within C4 crops. Using in situ canopy SIF and GPP measurements for both species along with leaf-level measurements, we found considerable differences in the SIF-GPP relationships between corn and miscanthus, with a stronger SIF-GPP relationship and higher slope of SIF-GPP observed in corn compared to miscanthus. These differences were mainly caused by leaf physiology. For miscanthus, high non-photochemical quenching (NPQ) under high light, temperature and water vapor deficit (VPD) conditions caused a large decline of fluorescence yield (Φ F ), which further led to a SIF midday depression and weakened the SIF-GPP relationship. The larger slope in corn than miscanthus was mainly due to its higher GPP in mid-summer, largely attributed to the higher leaf photosynthesis and less NPQ. Our results demonstrated variation of the SIF-GPP relationship within C4 crops and highlighted the importance of leaf physiology in determining canopy SIF behaviors and SIF-GPP relationships.

Different Responses of Solar-Induced Chlorophyll Fluorescence at the Red and Far-Red Bands and Gross Primary Productivity to Air Temperature for Winter Wheat

Solar-induced chlorophyll fluorescence (SIF) is closely related to the light-reaction process and has been recognized as a good indicator for tracking gross primary productivity (GPP). Nevertheless, it has not been widely examined how SIF and GPP respond to temperature. Here, we explored the linkage mechanisms between SIF and GPP in winter wheat based on continuous measurements of canopy SIF (cSIF), GPP, and meteorological data. To separately explore the structural and physiological mechanisms underlying the SIF–GPP relationship, we studied the temperature responses of the estimated light use efficiency (LUEp), canopy-level chlorophyll fluorescence yield (cSIFyield) and photosystem-level chlorophyll fluorescence yield (ΦF) estimated using canopy-scale remote sensing measurements. We found that GPP, red canopy SIF (cSIF688) and far-red canopy SIF (cSIF760) all exhibited a decreasing trend during overwintering periods. However, GPP and cSIF688 showed relatively more obvious changes in response to air temperature (Ta) than cSIF760 did. In addition, the LUEp responded sensitively to Ta (the correlation coefficient, r = 0.83, p-value &lt; 0.01). The cSIFyield_688 and ΦF_688 (ΦF at 688 nm) also exhibited significantly positive correlations with Ta (r &gt; 0.7, p-value &lt; 0.05), while cSIFyield_760 and ΦF_760 (ΦF at 760 nm) were weakly correlated with Ta (r &lt; 0.3, p-value &gt; 0.05) during overwintering periods. The results also show that LUEp was more sensitive to Ta than ΦF, which caused changes in the LUEp/ΦF ratio in response to Ta. By considering the influence of Ta, the GPP estimation based on the total SIF emitted at the photosystem level (tSIF) was improved (with R2 increased by more than 0.12 for tSIF760 and more than 0.05 for tSIF688). Therefore, our results indicate that the LUEp/ΦF ratio is affected by temperature conditions and highlights that the SIF–GPP model should consider the influence of temperature.

Estimation of Global Cropland Gross Primary Production from Satellite Observations by Integrating Water Availability Variable in Light-Use-Efficiency Model

Satellite-based models have been widely used to estimate gross primary production (GPP) of terrestrial ecosystems. Although they have many advantages for mapping spatiotemporal variations of regional or global GPP, the performance in agroecosystems is relatively poor. In this study, a light-use-efficiency model for cropland GPP estimation, named EF-LUE, driven by remote sensing data, was developed by integrating evaporative fraction (EF) as limiting factor accounting for soil water availability. Model parameters were optimized first using CO2 flux measurements by eddy covariance system from flux tower sites, and the optimized parameters were further spatially extrapolated according to climate zones for global cropland GPP estimation in 2001–2019. The major forcing datasets include the fraction of absorbed photosynthetically active radiation (FAPAR) data from the Copernicus Global Land Service System (CGLS) GEOV2 dataset, EF from the ETMonitor model, and meteorological forcing variables from ERA5 data. The EF-LUE model was first evaluated at flux tower site-level, and the results suggested that the proposed EF-LUE model and the LUE model without using water availability limiting factor, both driven by flux tower meteorology data, explained 82% and 74% of the temporal variations of GPP across crop sites, respectively. The overall KGE increased from 0.73 to 0.83, NSE increased from 0.73 to 0.81, and RMSE decreased from 2.87 to 2.39 g C m−2 d−1 in the estimated GPP after integrating EF in the LUE model. These improvements may be largely attributed to parameters optimized for different climatic zones and incorporating water availability limiting factor expressed by EF into the light-use-efficiency model. At global scale, the verification by GPP measurements from cropland flux tower sites showed that GPP estimated by the EF-LUE model driven by ERA5 reanalysis meteorological data and EF from ETMonitor had overall the highest R2, KGE, and NSE and the smallest RMSE over the four existing GPP datasets (MOD17 GPP, revised EC-LUE GPP, GOSIF GPP and PML-V2 GPP). The global GPP from the EF-LUE model could capture the significant negative GPP anomalies during drought or heat-wave events, indicating its ability to express the impacts of the water stress on cropland GPP.

Remotely Monitoring Vegetation Productivity in Two Contrasting Subtropical Forest Ecosystems Using Solar-Induced Chlorophyll Fluorescence

Subtropical forests can sequester a larger amount of atmospheric carbon dioxide (CO2) relative to other terrestrial ecosystems through photosynthetic activity and act as an important role in mitigating global climate warming. Compared with the model-based gross primary production (GPP) products, satellite-derived solar-induced fluorescence (SIF) opens a new window for quantification. Here, we used the remotely sensed SIF retrievals, two satellite-driven GPP products including MODIS (GPPMOD) and BESS (GPPBESS), and tower-based GPP measurements at two contrasting subtropical forests to provide a systematic analysis. Our results revealed that GPP and the associated environmental factors exhibited distinct seasonal patterns. However, the peak GPP values had large differences, with stronger GPP in the evergreen needleleaf forest site (8.76 ± 0.71 g C m−2 d−1) than that in the evergreen broadleaf forest site (5.71 ± 0.31 g C m−2 d−1). The satellite-derived SIF retrievals showed great potential in quantifying the variability in GPP, especially for the evergreen needleleaf forest with r reaching up to 0.909 (p &lt; 0.01). GPPMOD and GPPBESS showed distinctly different performances for the two subtropical forests, whereas the GPP estimates by exclusive use of satellite-based SIF data promised well to the tower-based GPP observations. Multi-year evaluation again confirmed the good performance of the SIF-based GPP estimates. These findings will provide an alternative framework for quantifying the magnitude of forest GPP and advance our understanding of the carbon sequestration capacity of subtropical forest ecosystems.

An ecological and methodological assessment of benthic gross primary production in northern lakes

AbstractBenthic gross primary production (GPP) is often the most important part of aquatic food webs in northern lakes, which are gradually warming and receiving increased terrestrial colored dissolved organic carbon loadings due to global change. Yet, measurements of benthic GPP are fairly uncommon, and methods and unit dimensions of benthic GPP are unstandardized and rarely compared. In this study, we measured benthic GPP in 27 headwater lakes from three regions in northern Sweden and analyzed potential constraining drivers of benthic GPPz rates at discrete depths and estimates of benthic GPP averages across the whole lake, as well as across the littoral zone. We also compared in situ measurements of benthic GPP averages across the whole lake with modeled values using the “autotrophic structuring model.” We found that benthic GPPz rates were best explained by, and positively related to, available light (i.e., a function of depth and water color) and temperature. Benthic GPP averages across the whole lake, on the contrary, were best explained by the relative size of the littoral area, which is a measure that combines lake bathymetry and water color. The comparison between in situ measured and modeled estimates of benthic GPP averages across the whole lake revealed that (1) the autotrophic structuring model underestimates GPP at low values and overestimates GPP at high values compared with measured data, and that (2) measured values were related to temperature, which is not included as a variable in the autotrophic structuring model. Considering future predicted changes impacting northern latitude lakes, our results suggest that increased lake water temperatures can to some extent mitigate the negative impacts of reduced light availability from lake browning on benthic GPPz rates. The combined impact of these changes on benthic GPP averages across the whole lake will depend on, and be moderated by, lake bathymetry determining the relative size of the littoral area.

Diurnal and Seasonal Dynamics of Solar‐Induced Chlorophyll Fluorescence, Vegetation Indices, and Gross Primary Productivity in the Boreal Forest

AbstractRemote sensing of solar‐induced chlorophyll fluorescence (SIF) provides a powerful proxy for gross primary productivity (GPP). It is particularly promising in boreal ecosystems where seasonal downregulation of photosynthesis occurs without significant changes in canopy structure or chlorophyll content. The use of SIF as a proxy for GPP is complicated by inherent non‐linearities due to both physical (illumination effects) and ecophysiological (light use efficiencies) controls at fine spatial (tower/leaf) and temporal (half‐hourly) scales. To study the SIF‐GPP relationship, we investigated the diurnal and seasonal dynamics of continuous tower‐based measurements of SIF, GPP, and common vegetation indices at the Southern Old Black Spruce Site (SOBS) in Saskatchewan, CA over the course of two years. We find that SIF outperforms other vegetation indices as a proxy for GPP at all temporal scales but shows a non‐linear relationship with GPP at a half‐hourly resolution. At small temporal scales, SIF and GPP are predominantly driven by light and non‐linearity between SIF and GPP is due to the light saturation of GPP. Averaged over daily and monthly scales, the relationship between SIF and GPP is linear due to a reduction in the observed PAR range. Seasonal changes in the light responses of SIF and GPP are driven by changes in light use efficiency which co‐vary with changes in temperature, while illumination and canopy structure partially linearize the SIF‐GPP relationship. Additionally, we find that the SIF‐GPP relationship has a seasonal dependency. Our results help clarify the utility of SIF for estimating carbon assimilation in boreal forests.