TY - JOUR
T1 - Across-model spread and shrinking in predicting peatland carbon dynamics under global change
AU - Hou, Enqing
AU - Ma, Shuang
AU - Huang, Yuanyuan
AU - Zhou, Yu
AU - Kim, Hyung Sub
AU - López-Blanco, Efrén
AU - Jiang, Lifen
AU - Xia, Jianyang
AU - Tao, Feng
AU - Williams, Christopher
AU - Williams, Mathew
AU - Ricciuto, Daniel
AU - Hanson, Paul J.
AU - Luo, Yiqi
N1 - Funding Information:
The authors thank Ning Wei for discussion on the work, and Professor Nigel Roulet and two anonymous reviewers for their constructive comments on this manuscript. This study was supported by the National Natural Science Foundation of China (32271644, 31870464), Guangdong Basic and Applied Basic Research Foundation (2022B1515020014), National Science Foundation Grant DEB (1655499; 2017884), US Department of Energy (DOE), Terrestrial Ecosystem Sciences Grant (DE‐SC0020227), and the subcontract 4000158404 from Oak Ridge National Laboratory (ORNL) to Northern Arizona University. ORNL's work was supported by DOE, Office of Science, Office of Biological and Environmental Research. ORNL is managed by UT‐Battelle, LLC, for DOE under contract DE‐AC05‐00OR22725. Hyung‐Sub Kim was supported by the National Research Foundation of Korea (NRF‐2018R1A2B6001012). Efrén López‐Blanco and Mathew Williams thank Luke Smallman for support in calibrating the DALEC2 model for the site. Efrén López‐Blanco has received funding from GreenFeedBack (Greenhouse gas fluxes and earth system feedbacks) funded by the European Union's HORIZON research and innovation program under grant agreement No 101056921. Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.
Funding Information:
The authors thank Ning Wei for discussion on the work, and Professor Nigel Roulet and two anonymous reviewers for their constructive comments on this manuscript. This study was supported by the National Natural Science Foundation of China (32271644, 31870464), Guangdong Basic and Applied Basic Research Foundation (2022B1515020014), National Science Foundation Grant DEB (1655499; 2017884), US Department of Energy (DOE), Terrestrial Ecosystem Sciences Grant (DE-SC0020227), and the subcontract 4000158404 from Oak Ridge National Laboratory (ORNL) to Northern Arizona University. ORNL's work was supported by DOE, Office of Science, Office of Biological and Environmental Research. ORNL is managed by UT-Battelle, LLC, for DOE under contract DE-AC05-00OR22725. Hyung-Sub Kim was supported by the National Research Foundation of Korea (NRF-2018R1A2B6001012). Efrén López-Blanco and Mathew Williams thank Luke Smallman for support in calibrating the DALEC2 model for the site. Efrén López-Blanco has received funding from GreenFeedBack (Greenhouse gas fluxes and earth system feedbacks) funded by the European Union's HORIZON research and innovation program under grant agreement No 101056921. Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.
Publisher Copyright:
© 2023 John Wiley & Sons Ltd.
PY - 2023/5
Y1 - 2023/5
N2 - Large across-model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across-model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO2 levels and five temperature levels. We developed a matrix-based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter-model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter-model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO2 (eCO2) concentrations in the measurements, simulated NEP responded positively to eCO2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across-model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.
AB - Large across-model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across-model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO2 levels and five temperature levels. We developed a matrix-based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter-model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter-model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO2 (eCO2) concentrations in the measurements, simulated NEP responded positively to eCO2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across-model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.
KW - SPRUCE experiment
KW - across-model spread
KW - carbon residence time
KW - environmental scalar
KW - land carbon dynamics
KW - matrix model
KW - traceability analysis
UR - http://www.scopus.com/inward/record.url?scp=85149743648&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85149743648&partnerID=8YFLogxK
U2 - 10.1111/gcb.16643
DO - 10.1111/gcb.16643
M3 - Article
C2 - 36799318
AN - SCOPUS:85149743648
SN - 1354-1013
VL - 29
SP - 2759
EP - 2775
JO - Global Change Biology
JF - Global Change Biology
IS - 10
ER -