Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

Light limitation regulates the response of autumn terrestrial carbon uptake to warming

Abstract

Global warming is projected to shift the phenology and increase the productivity of northern ecosystems1,2,3,4,5,6. Both changes will further feed back to climate through biophysical and biogeochemical processes and are critical for future prediction7,8. However, it remains unclear whether warming and the extended growing season, especially in autumn, would lead to increased net ecosystem carbon uptake9,10. Here we analyse satellite observations, field measurements and model simulations and show a prevalent radiation limitation on carbon uptake in northern ecosystems, especially in autumn. By comparing the start and end of the growing season estimated from vegetation indices and from solar-induced chlorophyll fluorescence (a proxy for gross primary production11,12 (GPP)), we find a greater change in greenness-based start and end of season than that from GPP, mostly caused by the radiation limitation on photosynthesis. This radiation limitation explains the contrasting responses of autumn net carbon exchanges to warming, using both eddy covariance measurements and model simulations from Coupled Model Intercomparison Project Phase 5. Regions with weak radiation limitation benefit more from warming and enhanced vegetation greenness in autumn, where GPP increases can outweigh the warming-induced respiration carbon losses. With continued warming, radiation limitation will increase and exert a strong upper bound on northern ecosystems to act as carbon sinks.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Trends of spring growth onset and autumn dormancy onset in 2001–2017.
Fig. 2: Phenological sensitivity (γ) between GPP and greenness for spring and autumn.
Fig. 3: Relationship between temperature anomaly and NEE for spring and autumn at 28 eddy covariance flux tower sites.
Fig. 4: Latitudinal comparison of radiation limitations derived from CSIF and CMIP5 models.
Fig. 5: Future radiation limitation simulated by CMIP5 models under the RCP8.5 scenario.

Similar content being viewed by others

Data availability

The MODIS VI and land surface temperature data are from https://lpdaac.usgs.gov/node/841 and https://lpdaac.usgs.gov/node/834, the CSIF data are from https://doi.org/10.17605/OSF.IO/8XQY6, the FLUXNET2015 dataset is from http://fluxnet.fluxdata.org/data/fluxnet2015-dataset/ and the CMIP5 model simulations are from https://esgf-node.llnl.gov/search/cmip5/. Raw data for Figs. 1–5 are available via FigShare at https://doi.org/10.6084/m9.figshare.11986764.

Code availability

The code used for the phenology retrieval and data analysis are available from Y.Z. upon request.

References

  1. Nemani, R. R. Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science 300, 1560–1563 (2003).

    CAS  Google Scholar 

  2. Zhu, Z. et al. Greening of the Earth and its drivers. Nat. Clim. Change 6, 791–795 (2016).

    CAS  Google Scholar 

  3. Jeong, S.-J., Ho, C.-H., Gim, H.-J. & Brown, M. E. Phenology shifts at start vs. end of growing season in temperate vegetation over the Northern Hemisphere for the period 1982–2008. Glob. Change Biol. 17, 2385–2399 (2011).

    Google Scholar 

  4. Buitenwerf, R., Rose, L. & Higgins, S. I. Three decades of multi-dimensional change in global leaf phenology. Nat. Clim. Change 5, 364–368 (2015).

    Google Scholar 

  5. Fridley, J. D., Lynn, J. S., Grime, J. P. & Askew, A. P. Longer growing seasons shift grassland vegetation towards more-productive species. Nat. Clim. Change 6, 865–868 (2016).

    Google Scholar 

  6. Schwartz, M. D., Ahas, R. & Aasa, A. Onset of spring starting earlier across the Northern Hemisphere. Glob. Change Biol. 12, 343–351 (2006).

    Google Scholar 

  7. Richardson, A. D. et al. Climate change, phenology, and phenological control of vegetation feedbacks to the climate system. Agric. For. Meteorol. 169, 156–173 (2013).

    Google Scholar 

  8. Peñuelas, J., Rutishauser, T. & Filella, I. Phenology feedbacks on climate change. Science 324, 887–888 (2009).

    Google Scholar 

  9. Piao, S. et al. Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 451, 49–52 (2008).

    CAS  Google Scholar 

  10. Keenan, T. F. et al. Net carbon uptake has increased through warming-induced changes in temperate forest phenology. Nat. Clim. Change 4, 598–604 (2014).

    CAS  Google Scholar 

  11. Sun, Y. et al. OCO-2 advances photosynthesis observation from space via solar-induced chlorophyll fluorescence. Science 358, eaam5747 (2017).

    Google Scholar 

  12. Guanter, L. et al. Global and time-resolved monitoring of crop photosynthesis with chlorophyll fluorescence. Proc. Natl Acad. Sci. USA 111, E1327–E1333 (2014).

    CAS  Google Scholar 

  13. Richardson, A. D. et al. Ecosystem warming extends vegetation activity but heightens vulnerability to cold temperatures. Nature 560, 368–371 (2018).

    CAS  Google Scholar 

  14. Ueyama, M., Iwata, H. & Harazono, Y. Autumn warming reduces the CO2 sink of a black spruce forest in interior Alaska based on a nine-year eddy covariance measurement. Glob. Change Biol. 20, 1161–1173 (2014).

    Google Scholar 

  15. Barichivich, J. et al. Large-scale variations in the vegetation growing season and annual cycle of atmospheric CO2 at high northern latitudes from 1950 to 2011. Glob. Change Biol. 19, 3167–3183 (2013).

    Google Scholar 

  16. Zhang, Y., Joiner, J., Alemohammad, S. H., Zhou, S. & Gentine, P. A global spatially contiguous solar-induced fluorescence (CSIF) dataset using neural networks. Biogeosciences 15, 5779–5800 (2018).

    CAS  Google Scholar 

  17. Porcar-Castell, A. et al. Linking chlorophyll a fluorescence to photosynthesis for remote sensing applications: mechanisms and challenges. J. Exp. Bot. 65, 4065–4095 (2014).

    CAS  Google Scholar 

  18. Krause, G. H. & Weis, E. Chlorophyll fluorescence and photosynthesis: the basics. Annu. Rev. Plant Physiol. Plant Mol. Biol. 42, 313–349 (1991).

    CAS  Google Scholar 

  19. Jeong, S. J. et al. Application of satellite solar-induced chlorophyll fluorescence to understanding large-scale variations in vegetation phenology and function over northern high latitude forests. Remote Sens. Environ. 190, 178–187 (2017).

    Google Scholar 

  20. Monteith, J. L. Solar radiation and productivity in tropical ecosystems. J. Appl. Ecol. 9, 747–766 (1972).

    Google Scholar 

  21. Wohlfahrt, G. et al. Trade-offs between global warming and day length on the start of the carbon uptake period in seasonally cold ecosystems. Geophys. Res. Lett. 40, 6136–6142 (2013).

    Google Scholar 

  22. Anav, A. et al. Evaluation of land surface models in reproducing satellite derived leaf area index over the high-latitude northern hemisphere. Part II: Earth system models. Remote Sens. 5, 3637–3661 (2013).

    Google Scholar 

  23. Hu, J., Moore, D. J. P., Burns, S. P. & Monson, R. K. Longer growing seasons lead to less carbon sequestration by a subalpine forest. Glob. Change Biol. 16, 771–783 (2010).

    Google Scholar 

  24. Schuur, E. A. G. et al. The effect of permafrost thaw on old carbon release and net carbon exchange from tundra. Nature 459, 556–559 (2009).

    CAS  Google Scholar 

  25. Davidson, E. A. & Janssens, I. A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440, 165–173 (2006).

    CAS  Google Scholar 

  26. Karl, T. R. & Trenberth, K. E. Modern global climate change. Science 302, 1719–1723 (2003).

    CAS  Google Scholar 

  27. Huang, M. et al. Velocity of change in vegetation productivity over northern high latitudes. Nat. Ecol. Evol. 1, 1649–1654 (2017).

    Google Scholar 

  28. Fu, Y. H. et al. Declining global warming effects on the phenology of spring leaf unfolding. Nature 526, 104–107 (2015).

    CAS  Google Scholar 

  29. Piao, S. et al. Weakening temperature control on the interannual variations of spring carbon uptake across northern lands. Nat. Clim. Change 7, 359–363 (2017).

    CAS  Google Scholar 

  30. Keenan, T. F. & Riley, W. J. Greening of the land surface in the world’s cold regions consistent with recent warming. Nat. Clim. Change 8, 825–828 (2018).

    CAS  Google Scholar 

  31. Shen, M. et al. No evidence of continuously advanced green-up dates in the Tibetan Plateau over the last decade. Proc. Natl Acad. Sci. USA 110, E2329–E2329 (2013).

    CAS  Google Scholar 

  32. Zhang, G. et al. Mapping paddy rice planting areas through time series analysis of MODIS land surface temperature and vegetation index data. ISPRS J. Photogramm. Remote Sens. 106, 157–171 (2015).

    Google Scholar 

  33. Zhang, Y. et al. Consistency between sun-induced chlorophyll fluorescence and gross primary production of vegetation in North America. Remote Sens. Environ. 183, 154–169 (2016).

    Google Scholar 

  34. Frankenberg, C. et al. New global observations of the terrestrial carbon cycle from GOSAT: patterns of plant fluorescence with gross primary productivity. Geophys. Res. Lett. 38, 1–6 (2011).

    Google Scholar 

  35. Zhang, Y. et al. On the relationship between sub-daily instantaneous and daily total gross primary production: implications for interpreting satellite-based SIF retrievals. Remote Sens. Environ. 205, 276–289 (2018).

    Google Scholar 

  36. Joiner, J. et al. Global monitoring of terrestrial chlorophyll fluorescence from moderate-spectral-resolution near-infrared satellite measurements: methodology, simulations, and application to GOME-2. Atmos. Meas. Tech. 6, 2803–2823 (2013).

    Google Scholar 

  37. Köhler, P. et al. Global retrievals of solar-induced chlorophyll fluorescence with TROPOMI: first results and intersensor comparison to OCO-2. Geophys. Res. Lett. 45, 10456–10463 (2018).

    Google Scholar 

  38. Zhang, Y., Joiner, J., Gentine, P. & Zhou, S. Reduced solar-induced chlorophyll fluorescence from GOME-2 during Amazon drought caused by dataset artifacts. Glob. Change Biol. 24, 2229–2230 (2018).

    Google Scholar 

  39. Reichstein, M. et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm. Glob. Change Biol. 11, 1424–1439 (2005).

    Google Scholar 

  40. Dušek, J., Čížková, H., Stellner, S., Czerný, R. & Květ, J. Fluctuating water table affects gross ecosystem production and gross radiation use efficiency in a sedge-grass marsh. Hydrobiologia 692, 57–66 (2012).

    Google Scholar 

  41. Pugh, C. A., Reed, D. E., Desai, A. R. & Sulman, B. N. Wetland flux controls: how does interacting water table levels and temperature influence carbon dioxide and methane fluxes in northern Wisconsin? Biogeochemistry 137, 15–25 (2018).

    CAS  Google Scholar 

  42. Migliavacca, M. et al. Using digital repeat photography and eddy covariance data to model grassland phenology and photosynthetic CO2 uptake. Agric. For. Meteorol. 151, 1325–1337 (2011).

    Google Scholar 

  43. Roerink, G. J., Menenti, M. & Verhoef, W. Reconstructing cloudfree NDVI composites using Fourier analysis of time series. Int. J. Remote Sens. 21, 1911–1917 (2000).

    Google Scholar 

  44. Piao, S., Fang, J., Zhou, L., Ciais, P. & Zhu, B. Variations in satellite-derived phenology in China’s temperate vegetation. Glob. Change Biol. 12, 672–685 (2006).

    Google Scholar 

  45. Gonsamo, A., Chen, J. M. & D’Odorico, P. Deriving land surface phenology indicators from CO2 eddy covariance measurements. Ecol. Indic. 29, 203–207 (2013).

    CAS  Google Scholar 

  46. Gu, L. Response of a deciduous forest to the Mount Pinatubo eruption: enhanced photosynthesis. Science 299, 2035–2038 (2003).

    Google Scholar 

  47. Sims, D. A. et al. Parallel adjustments in vegetation greenness and ecosystem CO2 exchange in response to drought in a Southern California chaparral ecosystem. Remote Sens. Environ. 103, 289–303 (2006).

    Google Scholar 

Download references

Acknowledgements

This work used eddy covariance data acquired and shared by the FLUXNET community, including these networks: AmeriFlux, AfriFlux, AsiaFlux, CarboAfrica, CarboEuropeIP, CarboItaly, CarboMont, ChinaFlux, Fluxnet-Canada, GreenGrass, ICOS, KoFlux, LBA, NECC, OzFlux-TERN, TCOS-Siberia and USCCC. The ERA-Interim reanalysis data are provided by ECMWF and processed by LSCE. The FLUXNET eddy covariance data processing and harmonization were carried out by the European Fluxes Database Cluster, AmeriFlux Management Project and Fluxdata project of FLUXNET, with the support of CDIAC and ICOS Ecosystem Thematic Center, and the OzFlux, ChinaFlux and AsiaFlux offices. Y.Z. acknowledges funding from NASA ROSES Advanced Information System Technology (grant no. 80NSSC17K0285). P.G. acknowledges support from NASA ROSES Terrestrial Hydrology (grant no. 80NSSC18K0998), NOAA MAPP (grant no. NA17OAR4310127) and European Research Council (grant no. CU18-3746). A.P.W. acknowledges support from the NASA Modeling, Analysis, and Prediction (MAP) programme (grant no. NASA 80NSSC17K0265).

Author information

Authors and Affiliations

Authors

Contributions

Y.Z. and P.G. conceived and designed the study. Y.Z. performed the data analysis. Y.Z., S.Z., R.C., A.P.W. and P.G. contributed to the interpretation of the results. Y.Z. drafted the manuscript. All authors participated in discussions and the editing of the manuscript.

Corresponding author

Correspondence to Yao Zhang.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Climate Change thanks Su-Jong Jeong and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Comparison between CSIF OCO-2 SIF and GPP at multiple scales.

a Comparison between OCO-2 SIF (SIFOCO-2), CSIF, NDVI and EVI for the northern ecosystems (>30°N). The error bands associated with SIFOCO-2 and CSIF represent the standard deviation of all observations in northern ecosystems, with the solid lines for the mean. For readability, error bands for NDVI and EVI are not shown. b-g Comparison between the spring growth onset and autumn dormancy onset derived from CSIF and flux tower estimated GPP at 40 sites from the FLUXNET2015 tier 1 dataset (Supplementary Table 3). b Locations of the 40 sites. c Comparison between eddy-covariance (EC) estimated GPP and CSIF at 4-day temporal resolution. d, e Comparison between SOS and EOS derived from EC estimated GPP and CSIF for all site-years. f, g Comparison between the interannual anomaly of SOS and EOS (ΔSOS and ΔEOS) for each site. The anomalies were calculated as the phenological dates for each year minus the multi-year average. Only sites with at least 5 years of observations were used for this interannual comparison.

Extended Data Fig. 2 Average growth onset and dormancy onset for 2001–2017 for four vegetation indicators.

a, b NDVI, c, d EVI and e, f CSIF, g, h GOME-2 SIF. The median value from four methods (wSpline-Thr, wHANTS-MR, wPolyfit-MR, wDLogistic) were shown. The growth onset and dormancy onset for GOME-2 SIF were derived from a mean seasonal cycle between 2007-2017, while others were derived for each year and averaged over the entire period.

Extended Data Fig. 3 MODIS VI processing schemes for phenology retrieval.

The nighttime LST (1:30 am local time) from Aqua satellite were used as it is close to the minimum LST.

Extended Data Fig. 4 A schematic diagram showing the extension of fPAR based growing season and its impact on GPP based growing season.

Left panels show how γ can be calculated. b shows a zoom in of autumn dormancy from schematic seasonal cycles a. Solid blue line in b represents the original fPAR×LUE and dashed-lines shows that with one day delay (corresponding to a one-day delay of greenness-based EOS). As a result, GPP based EOS delayed by γ day, from EOS0 to EOS1. Open circles and dots in b indicate the fPAR×LUE and PAR values for EOS0 and EOS1, respectively. At this new EOS date, fPAR×LUE increases by Δ(fPAR×LUE) and PAR decreases by ΔPAR. Note that Δ(fPAR×LUE) does not equal to ΔPAR, but rather \(\left[ {1 + \frac{{\Delta \left( {{\mathrm{fPAR}} \times {\mathrm{LUE}}} \right)}}{{{\mathrm{fPAR}} \times {\mathrm{LUE}}}}} \right]\left[ {1 - \frac{{\Delta {\mathrm{PAR}}}}{{{\mathrm{PAR}}}}} \right] = 1\) (see Methods). The right panels show two examples of how latitude and plant species (evergreen or deciduous) may affect γ values. Both PAR and fPAR×LUE in c and d were normalized by the PAR and fPAR×LUE values at EOS0, respectively. With increases in latitude, γ decreases, suggesting an increase in radiation limitation. Evergreen forests also exhibit relatively smaller γ values compared to deciduous species because of the relatively smaller percentage decrease of fPAR×LUE.

Extended Data Fig. 5 Comparison between radiation limitation derived from CSIF and that derived from GOME-2 SIF.

Left column is for spring and right column is for autumn.

Extended Data Fig. 6 Comparison between the radiation limitation using different threshold values to retrieve SOS and EOS.

Three thresholds for wSplineThr method to retrieve SOS and EOS were compared here. a,b threshold is minimum +30% of the seasonal magnitude, c,d threshold is minimum +25%, of the seasonal magnitude, e,f threshold is minimum +35% of the seasonal magintude. The difference in radiation limitation between the 35% and 25% thresholds is 0.9% (95% confidence interval: -2.1% to 14.0%).

Extended Data Fig. 7 Trend in the change of radiation limitation (1-γ) during 2001-2017.

Pie charts shows the areal percentages with decreasing trend (‘–‘) or increasing trend (‘+’).

Extended Data Fig. 8 Change of normalized derivative of PAR (parameter b in Equations. 6 or 9) along latitude and season.

Thick solid lines represent the median value of multi-year average SOS and EOS along latitude. Thin dashed and dotted lines represent the timings when b reaches maximum or minimum for each latitude bins, respectively. In the spring, advancing of the SOS would lead to increase of b and thus a greater radiation limitation. In the autumn, delaying of the EOS would lead to decrease of b until the minimum, also leading to a greater radiation limitation.

Extended Data Fig. 9 Correlation between temperature sensitivity of autumn NEE (\({\upgamma}_{\mathrm{T}}^{{\mathrm{NEE}}}\)) and radiation limitation.

Positive \({\upgamma}_{\mathrm{T}}^{{\mathrm{NEE}}}\) values indicate that ecosystem will release carbon (postive NEE) when temperature increases. Shaded area indicate 95% confidence interval of the regression.

Extended Data Fig. 10 Radiation limitation and temperature responses of autumn NEE for CMIP5 models.

Net ecosystem exchange (NEE) is calculated as the summation of autotropic and heterotopic respiration minus gross primary production, with a negative value indicates a carbon sink. Similar with inset in Fig. 3b, autumn NEE sensitivity is categorized into increase (blue) or decrease (red) carbon uptake with autumn warming. Only the gridcells that are limited by temperature is used for analysis (GPP based EOS delays when autumn temperature increases). Autumn is defined as 60 days interval around the GPP determined multi-year average EOS date. GPP, ecosystem respiration (Ra+Rh) and air temperature were interpolated from monthly to daily values to get the autumn average. Upper and lower boundary of the boxes represent the 25 and 75 percentiles, with the white lines in the middle the median values. The horizontal dashed green dashed line indicates the threshold to separate the negative and positive temperature sensitivity to ecosystem autumn NEE based on flux tower analysis. The grey shaded area represents the standard deviation of the threshold estimated by bootstrapping (n=2000). *** represent significance level at 0.001.

Supplementary information

Supplementary Information

Supplementary Figs. 1–6, Tables 1–4 and references.

Reporting Summary

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Commane, R., Zhou, S. et al. Light limitation regulates the response of autumn terrestrial carbon uptake to warming. Nat. Clim. Chang. 10, 739–743 (2020). https://doi.org/10.1038/s41558-020-0806-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41558-020-0806-0

Search

Quick links

Nature Briefing Microbiology

Sign up for the Nature Briefing: Microbiology newsletter — what matters in microbiology research, free to your inbox weekly.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing: Microbiology