Modelling decadal trends and the impact of extreme events on carbon fluxes in a temperate deciduous forest using a terrestrial biosphere model

<p>Changing climatic conditions pose a challenge to accurately estimating the carbon sequestration potential of terrestrial vegetation, which is often mediated by nitrogen availability. The close coupling between the nitrogen and carbon cycles controls plant productivity and shapes the structu...

Full description

Saved in:
Bibliographic Details
Main Authors: T. Thum, T. Miinalainen, O. Seppälä, H. Croft, C. Rogers, R. Staebler, S. Caldararu, S. Zaehle
Format: Article
Language:English
Published: Copernicus Publications 2025-04-01
Series:Biogeosciences
Online Access:https://bg.copernicus.org/articles/22/1781/2025/bg-22-1781-2025.pdf
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:<p>Changing climatic conditions pose a challenge to accurately estimating the carbon sequestration potential of terrestrial vegetation, which is often mediated by nitrogen availability. The close coupling between the nitrogen and carbon cycles controls plant productivity and shapes the structure and functional dynamics of ecosystems. However, how carbon and nitrogen interactions affect both carbon fluxes and plant functional traits in dynamic ecotones, which are experiencing biotic and abiotic changes, remains unclear. In this work, we use in situ measurements of leaf chlorophyll content (Chl<span class="inline-formula"><sub>Leaf</sub></span>, 2013–2016) and the leaf area index (LAI, 1998–2018) to parameterize the seasonal dynamics of the QUINCY (QUantifying Interactions between terrestrial Nutrient CYcles and the climate system) terrestrial biosphere model (TBM) to simulate the carbon fluxes at the Borden Forest Research Station flux tower site, Ontario, Canada, over 22 years from 1996 to 2018. Our goals are to assess the additional value of using Chl<span class="inline-formula"><sub>Leaf</sub></span> in the model parameterization, to study how well QUINCY can capture observed trends related to the carbon cycle at the site, and to investigate how well the processes associated with a drought year and its legacy effects are captured by the model.</p> <p>QUINCY was able to simulate leaf-level maximum carboxylation capacity (<span class="inline-formula"><i>V</i><sub>c(max),25</sub></span>), Chl<span class="inline-formula"><sub>Leaf</sub></span> and leaf nitrogen in a manner quite consistent with observations. The model with the improved parameterization captured observed daily gross primary production (GPP) well (<span class="inline-formula"><i>r</i><sup>2</sup>=0.80</span>, root mean square error (RMSE) <span class="inline-formula">=</span> 2.2 <span class="inline-formula"><math xmlns="http://www.w3.org/1998/Math/MathML" id="M7" display="inline" overflow="scroll" dspmath="mathml"><mrow class="unit"><mi mathvariant="normal">µ</mi><mi mathvariant="normal">mol</mi><mspace linebreak="nobreak" width="0.125em"/><msup><mi mathvariant="normal">m</mi><mrow><mo>-</mo><mn mathvariant="normal">2</mn></mrow></msup><mspace linebreak="nobreak" width="0.125em"/><msup><mi mathvariant="normal">s</mi><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup></mrow></math><span><svg:svg xmlns:svg="http://www.w3.org/2000/svg" width="64pt" height="15pt" class="svg-formula" dspmath="mathimg" md5hash="784b0d66eae26c4b23e1e5ea7f3270b5"><svg:image xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="bg-22-1781-2025-ie00001.svg" width="64pt" height="15pt" src="bg-22-1781-2025-ie00001.png"/></svg:svg></span></span>). Nevertheless, we found that although observed GPP increased significantly during the study period (22.4 <span class="inline-formula"><math xmlns="http://www.w3.org/1998/Math/MathML" id="M8" display="inline" overflow="scroll" dspmath="mathml"><mrow class="unit"><mi mathvariant="normal">gC</mi><mspace linebreak="nobreak" width="0.125em"/><msup><mi mathvariant="normal">m</mi><mn mathvariant="normal">2</mn></msup><mspace linebreak="nobreak" width="0.125em"/><msup><mi mathvariant="normal">yr</mi><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mspace width="0.125em" linebreak="nobreak"/><msup><mi mathvariant="normal">yr</mi><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup></mrow></math><span><svg:svg xmlns:svg="http://www.w3.org/2000/svg" width="73pt" height="15pt" class="svg-formula" dspmath="mathimg" md5hash="ca9a96800dbe2e55c82d4df620d5784c"><svg:image xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="bg-22-1781-2025-ie00002.svg" width="73pt" height="15pt" src="bg-22-1781-2025-ie00002.png"/></svg:svg></span></span>) and net ecosystem exchange (NEE) shifted towards a stronger sink, these trends were not captured in the model. Instead, QUINCY showed a significant increasing trend for total ecosystem respiration (TER) that was not present in the observations. The severe drought in 2007 strongly affected observed carbon fluxes, lowering both GPP and TER in the following year as well. QUINCY was able to capture some of the decrease in GPP and TER in 2007. However, the legacy effect of the drought in 2008 was not captured by the model. These results call for further work on representing legacy effects in TBMs, as these can have long-lasting impacts on ecosystem functioning.</p>
ISSN:1726-4170
1726-4189