A model parameter sensitivity approach to understand soil thermal and hydrological linkages and their influence on ecosystem carbon dynamics

Accelerated warming of the Arctic can affect the global climate system by thawing permafrost and exposing organic carbon in soils to decompose and release greenhouse gases into the atmosphere. We used a process-based biosphere model (DVM-DOS-TEM) designed to simulate biophysical and biogeochemical interactions between the soil, vegetation, and atmosphere. We varied soil and environmental parameters to assess the impact on cryohydrological and biogeochemical outputs in the model. We analyzed the responses of ecosystem carbon balances to permafrost thaw by running site-level simulations at two long-term tundra ecological monitoring sites in Alaska: Eight Mile Lake (EML) and Imnavait Creek Watershed (IMN), which are characterized by similar tussock tundra vegetation but differing soil drainage conditions and climate. Model outputs showed agreement with field observations at both sites for soil physical properties and ecosystem CO₂ fluxes. Model simulations of Net Ecosystem Exchange (NEE) showed an overestimation during the frozen season (higher CO₂ emissions) at EML with a mean NEE of 26.98 ± 4.83 gC/m²/month compared to observational mean of 22.01 ± 5.67 gC/m²/month, and during the fall months at IMN, with a modeled mean of 19.21 ± 7.49 gC/m²/month compared to observation mean of 11.9 ± 4.45 gC/m²/month. Our results underscore the importance of representing the impact of soil drainage conditions on the thawing of permafrost soils, particularly poorly drained soils, which will drive the magnitude of carbon released at sites across the high-latitude tundra. These findings can help improve predictions of net carbon releases from thawing permafrost, ultimately contributing to a better understanding of the impact of Arctic warming on the global climate system.