Search NASA⌕ Search

SEARCH · Search NASA

Results for “GPP”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Data for Attributing Differences of Solar-Induced Chlorophyll Fluorescence (SIF)-Gross Primary Production (GPP) Relationships between Two C4 Crops, Corn and Miscanthus

Information to characterize the solar-induced chlorophyll fluorescence (SIF)-gross primary production (GPP) relationship in C4 cropping systems remains limited. 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.

Feedstock Production↗

Can upscaling ground nadir SIF to eddy covariance footprint improve the relationship between SIF and GPP in croplands?

Ground solar-induced chlorophyll fluorescence (SIF) is important for the mechanistic understanding of the dynamics of vegetation gross primary production (GPP) at fine spatiotemporal scales. However, eddy covariance (EC) observations generally cover larger footprint areas than ground SIF observations (a bare fiber with nadir), and this footprint mismatch between nadir SIF and GPP could complicate the canopy SIF-GPP relationships. Here, we upscaled nadir SIF observations to EC footprint and investigated the change in SIF-GPP relationships after the upscaling in cropland. We included 13 site-years data in our study, with seven site-years corn, four site-years soybeans, and two site-years miscanthus, all located in the US Corn Belt. All sites’ crop nadir SIF observations collected from the automated FluoSpec2 system (a hemispheric-nadir system) were upscaled to the GPP footprint-based SIF using vegetation indices (VIs) calculated from high spatiotemporal satellite reflectance data. We found that SIF-GPP relationships were not substantially changed after upscaling nadir SIF to GPP footprint at our crop sites planted with corn, soybean, and miscanthus, with R 2 change after the upscaling ranging from -0.007 to 0.051 and root mean square error (RMSE) difference from -0.658 to 0.095 umol m –2 s –1 relative to original nadir SIF-GPP relationships across all the site-years. The variation of the SIF-GPP relationship within each species across different site-years was similar between the original nadir SIF and upscaled SIF. Different VIs, EC footprint models, and satellite data led to marginal differences in the SIF-GPP relationships when upscaling nadir SIF to EC footprint. Furthermore, our study provided a methodological framework to correct this spatial mismatch between ground nadir SIF and GPP observations for croplands and potentially for other ecosystems. Our results also demonstrated that the spatial mismatch between ground nadir SIF and GPP might not significantly affect the SIF-GPP relationship in cropland that are largely homogeneous.

54 ENVIRONMENTAL SCIENCES↗

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. Furthermore, 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.

54 ENVIRONMENTAL SCIENCES↗

Adjusting solar-induced fluorescence to nadir-viewing provides a better proxy for GPP

Solar-induced fluorescence (SIF) provides key information for inferring terrestrial gross primary productivity (GPP). However, there is little research on analyzing the performance of nadir-adjusted SIF (SIF nadir ) on GPP, compared to total SIF emitted by all leaves (SIF total ) for reducing the viewing angle effects and estimating GPP. Besides, there have been controversial opinions on whether SIF total is better correlated to GPP than raw SIF observations (SIF obs ). Here we systematically analyzed the relationship between raw/nadir/total SIF (i.e., SIF obs , SIF nadir and SIF total ) and GPP and investigated the underlying mechanism, using multi-angular field measurements and eddy covariance data in wheat and corn crops at sub-daily scale. We further compared the performance of SIF obs , SIF nadir and SIF total in estimating GPP using the TROPOspheric Monitoring Instrument (TROPOMI) SIF and concurrent AmeriFlux measurements at daily scale. Furthermore, the results indicate that diurnal SIF nadir has stronger correlations to GPP than SIF obs for field measurements, with the increase of mean coefficient of determination (R 2 ) by 0.05–0.07 for far-red band and 0.11–0.20 for red band. SIF nadir shows comparable performance with SIF total for both far-red and red bands. Although the viewing angle effects dominate the difference in estimating GPP between SIF nadir and SIF obs , the correlation between light use efficiency (LUE) and fesc further determines the different performance of SIF nadir and SIF total in estimating GPP. TROPOMI-based analysis further confirms that SIF nadir overall has higher correlations to AmeriFlux GPP than SIF obs for different plant functional types and shows similar performance with SIF total . Compared to SIF total , the estimation of SIF nadir independent of viewing angle effects does not require any canopy structure parameters, and thus offers promising potential for reliably estimating regional and global GPP.

54 ENVIRONMENTAL SCIENCES↗

Gross primary production (GPP) and red solar induced fluorescence (SIF) respond differently to light and seasonal environmental conditions in a subalpine conifer forest

The phenology of montane conifer forests is likely to shift in response to climate change and altered seasonal dynamics of light, temperature, and moisture. Solar-induced fluorescence (SIF) is expected to provide substantial improvement for mapping temporal changes in evergreen gross primary production (GPP) over greenness-based remote sensing indices. The utility of SIF to monitor seasonal changes in the phenology of conifer photosynthesis depends on the degree to which GPP and SIF respond in synchrony to key environmental drivers. However, to what extent SIF and GPP become decoupled by responding differently to the combined effects of light and other environmental conditions remains unknown. The goal of this study was to characterize the responses of GPP and SIFred to a suite of environmental drivers at the half-hour time scale and determine how these relationships change across seasons. We analyzed one year of tower-based SIFred and eddy covariance-derived GPP data from a conifer forest at Niwot Ridge, Colorado. We compared the light responses of GPP and SIFred across the year, finding that SIFred increased in response to light earlier in the year than did GPP. The light response of GPP had a positive temperature dependence in spring, and this dependency reversed in summer due to increased evaporative demand, while the light response of SIFred was less temperature dependent. Using artificial neural network ensemble analysis, we found that from spring to summer, SIFred did not exhibit a parallel response to the seasonally dynamic temperature and moisture controls on GPP. In summer SIFred was not correlated with canopy conductance, suggesting that SIF is less sensitive to stomatal control than GPP. Furthermore, our results suggest that, in conifers, photosystems begin to activate in spring prior to when water becomes available for photosynthesis, presenting a challenge for the use of SIF as a phenological indicator in conifer forests.

54 ENVIRONMENTAL SCIENCES↗

SIF-based GPP modeling for evergreen forests considering the seasonal variation in maximum photochemical efficiency

Solar-induced chlorophyll fluorescence (SIF) has shown great potential in estimating gross primary production (GPP). However, their quantitative relationship is not invariant, which undermines the reliability of empirical SIF-based GPP estimation at fine spatiotemporal scales, especially under extreme conditions. In this study, we developed a parsimonious mechanistic model for SIF-based GPP estimation in evergreen needle forests (ENF) by employing the Mechanistic Light Response framework and Eco-Evolutionary theory to describe the light and dark reactions during photosynthesis, respectively. Specifically, we found that considering the seasonal variation in a key parameter of the MLR framework, the maximum photochemical efficiency of photosystem II (Φ PSIImax ), can avoid the GPP overestimation in winter and early spring due to the relatively low environmental sensitivity of SIF. Compared to the estimates from other benchmark models, our GPP estimates were closer to the 1: 1 line and had higher accuracy (average R 2 = 0.86, RMSE=1.99 μmol m -2 s -1 ) across sites. Furthermore, the changes in the relationship between SIF and J (refers to the electron transport rate) contribute a lot to the dynamic SIF–GPP relationship in this study, while the J–GPP relationship is less variant when the temperature drops. Further, the seasonal variation in the SIF–J relationship, especially the reduction in its slope at low temperatures, is found largely explained by the Φ PSIImax . These results indicate the importance of the uncertainty caused by the variation in the SIF–J relationship for SIF-based GPP estimation, and the consideration of changes in Φ PSIImax under extreme conditions (such as severe winter in this study) is crucial for the improvement of GPP estimation via SIF.

59 BASIC BIOLOGICAL SCIENCES↗

Quantifying Northern High Latitude Gross Primary Productivity (GPP) Using Carbonyl Sulfide (OCS)

The northern high latitude (NHL, 40°N to 90°N) is where the second peak region of gross primary productivity (GPP) other than the tropics. The summer NHL GPP is about 80% of the tropical peak, but both regions are still highly uncertain. Carbonyl sulfide (OCS) provides an important proxy for photosynthetic carbon uptake. Here we optimize the OCS plant uptake fluxes across the NHL by fitting atmospheric concentration simulation with the GEOS-CHEM global transport model to the aircraft profiles acquired over Alaska during NASA's Carbon in Arctic Reservoirs Vulnerability Experiment (2012–2015). We use the empirical biome-specific linear relationship between OCS plant uptake flux and GPP to derive the six plant uptake OCS fluxes from different GPP data. Such GPP-based fluxes are used to drive the concentration simulations. We evaluate the simulations against the independent observations at two ground sites of Alaska. The optimized OCS fluxes suggest the NHL plant uptake OCS flux of -247 Gg S year -1 , about 25% stronger than the ensemble mean of the six GPP-based OCS fluxes. GPP-based OCS fluxes systematically underestimate the peak growing season across the NHL, while a subset of models predict early start of season in Alaska, consistent with previous studies of net ecosystem exchange. The OCS optimized GPP of 34 PgC yr -1 for NHL is also about 25% more than the ensembles mean from six GPP data. Further work is needed to fully understand the environmental and biotic drivers and quantify their rate of photosynthetic carbon uptake in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Radiance-based NIR v as a proxy for GPP of corn and soybean

Substantial uncertainty exists in daily and sub-daily gross primary production (GPP) estimation, which dampens accurate monitoring of the global carbon cycle. Here we find that near-infrared radiance of vegetation (NIR v,Rad ), defined as the product of observed NIR radiance and normalized difference vegetation index (NDVI), can accurately estimate corn and soybean GPP at daily and half-hourly time scales, benchmarked with multi-year tower-based GPP at three sites with different environmental and irrigation conditions. Overall, NIR v,Rad explains 84% and 78% variations of half-hourly GPP for corn and soybean, respectively, outperforming NIR reflectance of vegetation (NIR v,Ref ), enhanced vegetation index (EVI), and far-red solar-induced fluorescence (SIF760). The strong linear relationship between NIR v,Rad and absorbed photosynthetically active radiation by green leaves (APAR green ), and that between APAR green and GPP, explain the good NIR v,Rad -GPP relationship. The NIR v,Rad -GPP relationship is robust and consistent across sites. The scalability and simplicity of NIR v,Rad indicate a great potential to estimate daily or sub-daily GPP from high-resolution and/or long-term satellite remote sensing data.

54 ENVIRONMENTAL SCIENCES↗

SLOPE daily and 250 m gross primary productivity (GPP) for the CONUS, 2000-2019, Carbon Monitoring System (CMS)

The SatelLite Only Photosynthesis Estimation (SLOPE) GPP product is a daily, 250 m GPP dataset covering the Contiguous United States (CONUS) from 2000 to present with 1 day latency. Gross primary productivity (GPP) quantifies the amount of carbon dioxide (CO2) fixed by plants through photosynthesis. Although as a key quantity of terrestrial ecosystems, there is a lack of high-spatial-and-temporal-resolution, real-time, and observation-based GPP products. This product has been developed to address this critical gap, leveraging a number of MODIS land and atmosphere products. There are three distinct features of the SLOPE GPP production algorithm. (1) SLOPE couples machine learning models with MODIS atmosphere and land products to accurately estimate PAR. (2) SLOPE couples highly efficient and pragmatic gap-filling and filtering algorithms with surface reflectance acquired by both Terra and Aqua MODIS satellites to derive a soil-adjusted NIRv (SANIRv) dataset. (3) SLOPE couples a temporal pattern recognition approach with a long-term Crop Data Layer (CDL) product to predict dynamic C4 crop fraction. PAR, SANIRv and C4 fraction are used to drive a parsimonious model with only two slope parameters to estimate GPP along with a quantitative uncertainty on a per-pixel and daily basis. The slope GPP product has a R2 = 0.84 and a root-mean-square error (RMSE) of 1.65 gC m-2 d-1, evaluated against from 50 AmeriFlux eddy covariance sites (332 site-years). Archived data from 2000 to 2019 are publicly available at the NASA's Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC). Data from 2020 are available from the authors upon request. All data are projected in the standard MODIS Land Integerized Sinusoidal tile map projection. Each processing tile is in size of 4800 pixels by 4800 pixels, representing approximately 1200 km by 1200 km land region. In addition to the GPP product, SLOPE PAR, SANIRV, and C4 fraction, along with their uncertainties, are also released.

Jiang, Chongya↗

Escape Ratio Contributes More Than Fluorescence Yield to SIF-GPP Relationship Over Crops and Rainforest

Solar-induced chlorophyll fluorescence (SIF) is an effective indicator to track the gross primary productivity (GPP). However, there is still a lack of a clear understanding for the contribution of physiological and structural factors to the SIF–GPP relationship at the canopy scale. To quantify the influence of different SIF components, particularly the photon escape ratio ($f$ esc ) and fluorescence yield (Φ F ), on the SIF–GPP relationship, this study evaluated the performance of various vegetation indices (VIs), SIF components, light use efficiency (LUE), and GPP over two typical biomes (crops and rainforest), using a range of satellite remote sensing products. In August of each year from 2018 to 2020, both SIF and GPP over United States (U.S.) Corn Belt are higher than those over the Amazon rainforest, attributed to the consistent pattern of higher $f$ esc and LUE over crops than over rainforest, as well as Φ F . Furthermore, the structural signals represented by $f$ esc (R = 0.41–0.64) can better capture the LUE variations than Φ F (R = 0.10–0.30) for each biome. In conclusion, this study highlights that $f$ esc , determined by canopy structure, has great potential to capture LUE and GPP changes within and across biomes.

60 APPLIED LIFE SCIENCES↗

Comparison of Phenology Estimated from Reflectance-Based Indices and Solar-Induced Chlorophyll Fluorescence (SIF) Observations in a Temperate Forest Using GPP-Based Phenology as the Standard

We assessed the performance of reflectance-based vegetation indices and solar-induced chlorophyll fluorescence (SIF) datasets with various spatial and temporal resolutions in monitoring the Gross Primary Production (GPP)-based phenology in a temperate deciduous forest. The reflectance-based indices include the green chromatic coordinate (GCC), field measured and satellite remotely sensed Normalized Difference Vegetation Index (NDVI); and the SIF datasets include ground-based measurement and satellite-based products. We found that, if negative impacts due to coarse spatial and temporal resolutions are effectively reduced, all these data can serve as good indicators of phenological metrics for spring. However, the autumn phenological metrics derived from all reflectance-based datasets are later than the those derived from ground-based GPP estimates (flux sites). This is because the reflectance-based observations estimate phenology by tracking physiological properties including leaf area index (LAI) and leaf chlorophyll content (Chl), which does not reflect instantaneous changes in phenophase transitions, and thus the estimated fall phenological events may be later than GPP-based phenology. In contrast, we found that SIF has a good potential to track seasonal transition of photosynthetic activities in both spring and fall seasons. The advantage of SIF in estimating the GPP-based phenology lies in its inherent link to photosynthesis activities such that SIF can respond quickly to all factors regulating phenological events. Despite uncertainties in phenological metrics estimated from current spaceborne SIF observations due to their coarse spatial and temporal resolutions, dates in middle spring and autumn—the two most important metrics—can still be reasonably estimated from satellite SIF. Our study reveals that SIF provides a better way to monitor GPP-based phenological metrics.

59 BASIC BIOLOGICAL SCIENCES↗

Accounting for Changes in Radiation Improves the Ability of SIF to Track Water Stress-Induced Losses in Summer GPP in a Temperate Deciduous Forest

Global observations of solar-induced chlorophyll fluorescence (SIF) are available from multiple satellite platforms, and SIF is increasingly used as a proxy for photosynthetic activity and ecosystem productivity. Because the relationship between SIF and gross primary productivity (GPP) depends on a variety of factors including ecosystem type and environmental conditions, it is necessary to study SIF observations across various spatiotemporal scales and ecosystems. To explore how SIF signals relate to productivity over a temperate deciduous forest, we deployed a PhotoSpec spectrometer system at the University of Michigan Biological Station AmeriFlux site (US-UMB) in the northern Lower Peninsula of Michigan during the 2018 and 2019 growing seasons. We found that SIF correlated with GPP across diurnal and seasonal cycles (R 2 = 0.61 and 0.64 for 90-min- and daily-averaged data), but that SIF signals were more strongly related to downwelling radiation than GPP (R 2 = 0.91 for daily-averaged data). The dependence of SIF on radiation obscured the impact of intraseasonal drought in the SIF timeseries, but drought stress was apparent as a decrease in relative SIF, which exhibited a stronger correlation with GPP (R 2 = 0.56) than other remotely sensed data over the drought period. These results highlight the potential of SIF for detecting stress-induced losses in forest productivity. Additionally, we found that the red:far-red SIF ratio did not exhibit a response to water stress-induced losses in productivity, but was largely driven by seasonal and interannual changes in canopy structure, as well as by synoptic changes in downwelling radiation.

54 ENVIRONMENTAL SCIENCES↗

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗

Assimilation of multiple datasets results in large differences in regional- to global-scale NEE and GPP budgets simulated by a terrestrial biosphere model

In spite of the importance of land ecosystems in offsetting carbon dioxide emissions released by anthropogenic activities into the atmosphere, the spatiotemporal dynamics of terrestrial carbon fluxes remain largely uncertain at regional to global scales. Over the past decade, data assimilation (DA) techniques have grown in importance for improving these fluxes simulated by terrestrial biosphere models (TBMs), by optimizing model parameter values while also pinpointing possible parameterization deficiencies. Although the joint assimilation of multiple data streams is expected to constrain a wider range of model processes, their actual benefits in terms of reduction in model uncertainty are still under-researched, also given the technical challenges. In this study, we investigated with a consistent DA framework and the ORCHIDEE-LMDz TBM–atmosphere model how the assimilation of different combinations of data streams may result in different regional to global carbon budgets. To do so, we performed comprehensive DA experiments where three datasets (in situ measurements of net carbon exchange and latent heat fluxes, spaceborne estimates of the normalized difference vegetation index, and atmospheric CO 2 concentration data measured at stations) were assimilated alone or simultaneously. We thus evaluated their complementarity and usefulness to constrain net and gross C land fluxes. We found that a major challenge in improving the spatial distribution of the land C sinks and sources with atmospheric CO 2 data relates to the correction of the soil carbon imbalance.

54 ENVIRONMENTAL SCIENCES↗

Underestimates of Grassland Gross Primary Production in MODIS Standard Products

As the biggest carbon flux of terrestrial ecosystems from photosynthesis, gross primary productivity (GPP) is an important indicator in understanding the carbon cycle and biogeochemical process of terrestrial ecosystems. Despite advances in remote sensing-based GPP modeling, spatial and temporal variations of GPP are still uncertain especially under extreme climate conditions such as droughts. As the only official products of global spatially explicit GPP, MOD17A2H (GPP MOD ) has been widely used to assess the variations of carbon uptake of terrestrial ecosystems. However, systematic assessment of its performance has rarely been conducted especially for the grassland ecosystems where inter-annual variability is high. Based on a collection of GPP datasets (GPP EC ) from a global network of eddy covariance towers (FluxNet), we compared GPP MOD and GPP EC at all FluxNet grassland sites with more than five years of observations. We evaluated the performance and robustness of GPP MOD in different grassland biomes (tropical, temperate, and alpine) by using a bootstrapping method for calculating 95% confident intervals (CI) for the linear regression slope, coefficients of determination (R 2 ), and root mean square errors (RMSE). We found that GPP MOD generally underestimated GPP by about 34% across all biomes despite a significant relationship (R 2 = 0.66 (CI, 0.63–0.69), RMSE = 2.46 (2.33–2.58) g Cm -2 day -1 ) for the three grassland biomes. GPP MOD had varied performances with R2 values of 0.72 (0.68–0.75) (temperate), 0.64 (0.59–0.68) (alpine), and 0.40 (0.27–0.52) (tropical). Thus, GPP MOD performed better in low GPP situations (e.g., temperate grassland type), which further indicated that GPP MOD underestimated GPP. The underestimation of GPP could be partly attributed to the biased maximum light use efficiency (ε max ) values of different grassland biomes. The uncertainty of the fraction of absorbed photosynthetically active radiation (FPAR) and the water scalar based on the vapor pressure deficit (VPD) could have other reasons for the underestimation. Therefore, more accurate estimates of GPP for different grassland biomes should consider improvements in ε max , FPAR, and the VPD scalar. Our results suggest that the community should be cautious when using MODIS GPP products to examine spatial and temporal variations of carbon fluxes.

54 ENVIRONMENTAL SCIENCES↗

Multi‐model analysis of climate impacts on plant photosynthesis in China during 2000–2015

Abstract Differences, arising from differences in gross primary production (GPP) model structures and driving forces, have fuelled arguments concerning interannual changes of GPP in China since 2000. To better investigate the interannual variability of GPP and its covariance with climate factors in China, this study adopted a multi‐model analysis based on three GPP models (i.e., Terrestrial Ecosystem Carbon flux model [TEC], Breathing Earth System Simulator model [BESS], and MOD17 GPP model). The results show that annual GPP in China increased by 0.021–0.057 Pg C year −1 from 2000 to 2015 attributable to atmospheric‐CO 2 fertilization effects and favourable climate change, that is, increasing precipitation ( P r ) and temperature ( T a ). However, northern China and southern China had a large difference in the amplitude of these GPP changes; annual GPP increased by 0.017–0.039 Pg C year −1 in northern China but only 0.001–0.018 Pg C year −1 in southern China. Northern China and southern China occupy contrasting climate zones and this contrast produced different interannual variability of GPP through different mechanisms. Northern China has a dry climate with GPP changes sensitive to P r . As a result, more P r along with higher T a in northern China produced the strong uptrend of GPP from 2000 to 2015. In contrast, southern China has a wet climate with its GPP sensitive to solar radiation and T a . For the interval of 2000–2015, decreasing radiation plus drought exerted a negative influence on GPP in southern China. This study highlights the diverse mechanisms in which climate change affects GPP in dry and wet climate zones. A robust multi‐model analysis is preferred to reduce uncertainties arising from a single GPP model and its driving data.

Yan, Hao↗