Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change-carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables. We showed a case study on modeled soil carbon storage under two forcing datasets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin-up.
Xingjie Lu, Zhenggang Du, Yuanyuan Huang, David Lawrence, Erik Kluzek, Nathan Collier, Danica Lombardozzi, Negin Sobhani, Edward A. G. Schuur, Yiqi Luo. 2020. Full implementation of matrix approach to biogeochemistry module of Community Land Model version 5 (CLM5). JAMES, under review
Code and Data (Download Files)