skip to main content


Title: Solar‐induced chlorophyll fluorescence is strongly correlated with terrestrial photosynthesis for a wide variety of biomes: First global analysis based on OCO‐2 and flux tower observations
Abstract

Solar‐induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite‐observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3 km × 2.25 km) measured from the Orbiting Carbon Observatory‐2 (OCO‐2) provides the first opportunity to examine the SIF–GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO‐2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO‐2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757 nm (R2 = 0.72,p < 0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes (R2 = 0.57–0.79,p < 0.0001) except evergreen broadleaf forests (R2 = 0.16,p < 0.05) at the daily timescale. A higher slope was found for C4grasslands and croplands than for C3ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome‐specific SIF–GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO‐2 SIF generally had a better performance for predicting GPP than satellite‐derived vegetation indices and a light use efficiency model. The universal SIF–GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO‐2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies.

 
more » « less
Award ID(s):
1637459
NSF-PAR ID:
10059121
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Global Change Biology
Volume:
24
Issue:
9
ISSN:
1354-1013
Page Range / eLocation ID:
p. 3990-4008
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Accurate estimation of terrestrial gross primary productivity (GPP) remains a challenge despite its importance in the global carbon cycle. Chlorophyll fluorescence (ChlF) has been recently adopted to understand photosynthesis and its response to the environment, particularly with remote sensing data. However, it remains unclear how ChlF and photosynthesis are linked at different spatial scales across the growing season. We examined seasonal relationships between ChlF and photosynthesis at the leaf, canopy, and ecosystem scales and explored how leaf‐level ChlF was linked with canopy‐scale solar‐induced chlorophyll fluorescence (SIF) in a temperate deciduous forest at Harvard Forest, Massachusetts,USA. Our results show that ChlF captured the seasonal variations of photosynthesis with significant linear relationships between ChlF and photosynthesis across the growing season over different spatial scales (R= 0.73, 0.77, and 0.86 at leaf, canopy, and satellite scales, respectively;P < 0.0001). We developed a model to estimateGPPfrom the tower‐based measurement ofSIFand leaf‐level ChlF parameters. The estimation ofGPPfrom this model agreed well with flux tower observations ofGPP(R= 0.68;P < 0.0001), demonstrating the potential ofSIFfor modelingGPP. At the leaf scale, we found that leafFq/Fm, the fraction of absorbed photons that are used for photochemistry for a light‐adapted measurement from a pulse amplitude modulation fluorometer, was the best leaf fluorescence parameter to correlate with canopySIFyield (SIF/APAR,R= 0.79;P < 0.0001). We also found that canopySIFandSIF‐derivedGPP(GPPSIF) were strongly correlated to leaf‐level biochemistry and canopy structure, including chlorophyll content (R= 0.65 for canopyGPPSIFand chlorophyll content;P < 0.0001), leaf area index (LAI) (R= 0.35 for canopyGPPSIFandLAI;P < 0.0001), and normalized difference vegetation index (NDVI) (R= 0.36 for canopyGPPSIFandNDVI;P < 0.0001). Our results suggest that ChlF can be a powerful tool to track photosynthetic rates at leaf, canopy, and ecosystem scales.

     
    more » « less
  2. Abstract

    Arctic‐boreal landscapes are experiencing profound warming, along with changes in ecosystem moisture status and disturbance from fire. This region is of global importance in terms of carbon feedbacks to climate, yet the sign (sink or source) and magnitude of the Arctic‐boreal carbon budget within recent years remains highly uncertain. Here, we provide new estimates of recent (2003–2015) vegetation gross primary productivity (GPP), ecosystem respiration (Reco), net ecosystem CO2exchange (NEE;Reco − GPP), and terrestrial methane (CH4) emissions for the Arctic‐boreal zone using a satellite data‐driven process‐model for northern ecosystems (TCFM‐Arctic), calibrated and evaluated using measurements from >60 tower eddy covariance (EC) sites. We used TCFM‐Arctic to obtain daily 1‐km2flux estimates and annual carbon budgets for the pan‐Arctic‐boreal region. Across the domain, the model indicated an overall average NEE sink of −850 Tg CO2‐C year−1. Eurasian boreal zones, especially those in Siberia, contributed to a majority of the net sink. In contrast, the tundra biome was relatively carbon neutral (ranging from small sink to source). Regional CH4emissions from tundra and boreal wetlands (not accounting for aquatic CH4) were estimated at 35 Tg CH4‐C year−1. Accounting for additional emissions from open water aquatic bodies and from fire, using available estimates from the literature, reduced the total regional NEE sink by 21% and shifted many far northern tundra landscapes, and some boreal forests, to a net carbon source. This assessment, based on in situ observations and models, improves our understanding of the high‐latitude carbon status and also indicates a continued need for integrated site‐to‐regional assessments to monitor the vulnerability of these ecosystems to climate change.

     
    more » « less
  3. Abstract

    Recent advances in satellite observations of solar‐induced chlorophyll fluorescence (SIF) provide a new opportunity to constrain the simulation of terrestrial gross primary productivity (GPP). Accurate representation of the processes driving SIF emission and its radiative transfer to remote sensing sensors is an essential prerequisite for data assimilation. Recently, SIF simulations have been incorporated into several land surface models, but the scaling of SIF from leaf‐level to canopy‐level is usually not well‐represented. Here, we incorporate the simulation of far‐red SIF observed at nadir into the Community Land Model version 5 (CLM5). Leaf‐level fluorescence yield was simulated by a parametric simplification of the Soil Canopy‐Observation of Photosynthesis and Energy fluxes model (SCOPE). And an efficient and accurate method based on escape probability is developed to scale SIF from leaf‐level to top‐of‐canopy while taking clumping and the radiative transfer processes into account. SIF simulated by CLM5 and SCOPE agreed well at sites except one in needleleaf forest (R2 > 0.91, root‐mean‐square error <0.19 W⋅m−2⋅sr−1⋅μm−1), and captured the day‐to‐day variation of tower‐measured SIF at temperate forest sites (R2 > 0.68). At the global scale, simulated SIF generally captured the spatial and seasonal patterns of satellite‐observed SIF. Factors including the fluorescence emission model, clumping, bidirectional effect, and leaf optical properties had considerable impacts on SIF simulation, and the discrepancies between simulate d and observed SIF varied with plant functional type. By improving the representation of radiative transfer for SIF simulation, our model allows better comparisons between simulated and observed SIF toward constraining GPP simulations.

     
    more » « less
  4. Abstract

    Climate change is rapidly altering composition, structure, and functioning of the boreal biome, across North America often broadly categorized into ecoregions. The resulting complex changes in different ecoregions present a challenge for efforts to accurately simulate carbon dioxide (CO2) and energy exchanges between boreal forests and the atmosphere with terrestrial ecosystem models (TEMs). Eddy covariance measurements provide valuable information for evaluating the performance of TEMs and guiding their development. Here, we compiled a boreal forest model benchmarking dataset for North America by harmonizing eddy covariance and supporting measurements from eight black spruce (Picea mariana)-dominated, mature forest stands. The eight forest stands, located in six boreal ecoregions of North America, differ in stand characteristics, disturbance history, climate, permafrost conditions and soil properties. By compiling various data streams, the benchmarking dataset comprises data to parameterize, force, and evaluate TEMs. Specifically, it includes half-hourly, gap-filled meteorological forcing data, ancillary data essential for model parameterization, and half-hourly, gap-filled or partitioned component flux data on CO2(net ecosystem production, gross primary production [GPP], and ecosystem respiration [ER]) and energy (latent [LE] and sensible heat [H]) and their daily aggregates screened based on half-hourly gap-filling quality criteria. We present a case study with the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) to: (1) demonstrate the utility of our dataset to benchmark TEMs and (2) provide guidance for model development and refinement. Model skill was evaluated using several statistical metrics and further examined through the flux responses to their environmental controls. Our results suggest that CLASSIC tended to overestimate GPP and ER among all stands. Model performance regarding the energy fluxes (i.e., LE and H) varied greatly among the stands and exhibited a moderate correlation with latitude. We identified strong relationships between simulated fluxes and their environmental controls except for H, thus highlighting current strengths and limitations of CLASSIC.

     
    more » « less
  5. Abstract

    We mapped tidal wetland gross primary production (GPP) with unprecedented detail for multiple wetland types across the continental United States (CONUS) at 16‐day intervals for the years 2000–2019. To accomplish this task, we developed the spatially explicit Blue Carbon (BC) model, which combined tidal wetland cover and field‐based eddy covariance tower data into a single Bayesian framework, and used a super computer network and remote sensing imagery (Moderate Resolution Imaging Spectroradiometer Enhanced Vegetation Index). We found a strong fit between the BC model and eddy covariance data from 10 different towers (r2= 0.83,p< 0.001, root‐mean‐square error = 1.22 g C/m2/day, average error was 7% with a mean bias of nearly zero). When compared with NASA's MOD17 GPP product, which uses a generalized terrestrial algorithm, the BC model reduced error by approximately half (MOD17 hadr2= 0.45,p< 0.001, root‐mean‐square error of 3.38 g C/m2/day, average error of 15%). The BC model also included mixed pixels in areas not covered by MOD17, which comprised approximately 16.8% of CONUS tidal wetland GPP. Results showed that across CONUS between 2000 and 2019, the average daily GPP per m2was 4.32 ± 2.45 g C/m2/day. The total annual GPP for the CONUS was 39.65 ± 0.89 Tg C/year. GPP for the Gulf Coast was nearly double that of the Atlantic and Pacific Coasts combined. Louisiana alone accounted for 15.78 ± 0.75 Tg C/year, with its Atchafalaya/Vermillion Bay basin at 4.72 ± 0.14 Tg C/year. The BC model provides a robust platform for integrating data from disparate sources and exploring regional trends in GPP across tidal wetlands.

     
    more » « less