Carbon-water coupling across scales

Summary of our research

In our new paper published in Nature Geoscience (view it here), we identified a common bias in different state-of-the-art satellite-based estimates of terrestrial photosynthesis, and found that this bias is closely related to droughts and their apparent impact on the functioning of vegetation. In other words, common satellite-based vegetation monitoring methods underestimate the impact of droughts.

What are the implications of this bias for our understanding of the terrestrial C cycle and its year-to-year variability? We investigated this question using a new global implementation of the P-model (open-access code available here). Our results suggest that accurately capturing soil moisture stress in photosynthesis modelling is crucial for estimating the year-to-year variability in vegetation productivity and the impacts of drought extreme events. Not capturing this effect accurately leads to a substantial underestimation of this variability and to a systematic underestimation of drought extreme event impacts on vegetation productivity.

However, in our simulations, the locally strong soil moisture effect on vegetation productivity and its variability does not translate into a similarly strong water-control on the global-scale variability in terrestrial gross primary productivity. This is due to the compensating effects of different continents that are not synchronously drought-stressed and tend to compensate each other.

In view of a recent study by Humphrey et al. (2018) who found a clear relationship between terrestrial water storage and the atmospheric CO2 growth rate, this raises a fundamental question about the evident water-carbon coupling at the global scale: Is it an apparent phenomenon where terrestrial water storage is a proxy for some other factor that controls vegetation functioning? Or are we missing some element of plant access to water that has a longer memory and stronger year-to-year variability, so that the simulated compensating effects from different continents that we found would no longer neatly compensate each other at the global scale? This is an exciting lead for future research…

Data and code

  • Model output data from site-scale and global P-model simulations are available on Zenodo.
  • Code for the entire analysis and producing figures is available on Zenodo and on github. Compile (knit) RMarkdown file knit_soilm_global.Rmd for a full reproduction of results or si_soilm_global.Rmd for re-creating the Supplementary Information to our paper.
Benjamin Stocker
Benjamin Stocker
Group leader, Asst. Prof.

Heliocentrist human being.