Exploring Terrestrial CO₂, ¹³CO₂, and ¹⁴CO₂ Fluxes with the LPJ Dynamic Vegetation Model

B. Fischer-Femal¹, B. Poulter¹, S. Basu^2,3, J. Miller⁴, A. Kaushik^5,4 and S. Lehman⁶

¹NASA Goddard Space Flight Center (GSFC), Greenbelt, MD 20771; 541-285-8204, E-mail: brenden.j.fischer-femal@nasa.gov
²NASA Goddard Space Flight Center (GSFC), Global Monitoring and Assimilation Office, Greenbelt, MD 20771
³Earth System Science Interdisciplinary Center, University of Maryland, College Park, MD 20740
⁴NOAA Global Monitoring Laboratory (GML), Boulder, CO 80305
⁵Cooperative Institute for Research in Environmental Sciences (CIRES), University of Colorado, Boulder, CO 80309
⁶Institute of Arctic and Alpine Research (INSTAAR), University of Colorado, Boulder, CO 80309

Predicting future terrestrial carbon storage is essential to understanding feedbacks within the carbon-climate system and informing policies concerning CO₂ emissions mitigation and land management. To predict the effect of warmer temperatures and higher CO₂ on vegetation and soils, process-based models require accurate numerical representations of key processes, including photosynthesis, respiration, turnover rates of below-ground carbon, and plant community change. However, representations of these processes carry significant uncertainties and differ considerably between models. In this study, we add carbon isotopes (¹³C and ¹⁴C) to the Lund-Potsdam-Jena (LPJ) dynamic global vegetation model (DGVM), which also includes representations of fire, agricultural management, and vegetation population dynamics. Simulation of carbon isotope fluxes in LPJ provides for the application of additional observational constraints that can help separate the influences of gross primary productivity (GPP) and ecosystem respiration (R_ECO) on the seasonal dynamics of net ecosystem exchange (NEE) as well as providing additional information on plant water use efficiency and below-ground carbon turnover rates. Here we simulate land fluxes for CO₂, ¹³CO₂, and ¹⁴CO₂ from 1901-2022 using LPJ in a set of model runs changing parameters related to the temperature response of photosynthesis and respiration, isotope fractionation, age class, permafrost, and land use change. Modeled global land-atmosphere isotope disequilibrium fluxes of ¹³C and ¹⁴C from 1901-2022 are compared to previous estimates of global isotope disequilibrium fluxes. Simulated fluxes of bulk CO₂, ¹³C, and ¹⁴C associated with unique combinations of model parameters and assumptions are then transported using the TM5 atmospheric transport model and evaluated by comparison to observations of all three quantities (CO₂, d¹³CO₂, and D¹⁴CO₂) in measurements obtained from NOAA global monitoring sites around the world.

Figure 1. Modelled δ¹³C values of global carbon pools (vegetation, soil, and litter) and fluxes (gross primary productivity, ecosystem respiration, heterotrophic respiration, fire, and landuse change) for 3 LPJ simulations which test the effects of permafrost (pf), ageclass (ac), and Vcmax temperature scaling (KK) and all using the TRENDY S3 protocol.