Merging a mechanistic enzymatic model of soil heterotrophic respiration into an ecosystem model in two AmeriFlux sites of northeastern USA.
The objective of this work was to insert into an ecosystem model a more mechanistic, but still parsimonious, model of environmental factors controlling heterotrohic respiration (Rh) and evaluate the model performance in terms of soil and ecosystem respiration. The Dual Arrhenius and Michaelis-Menten (DAMM) model simulates Rh using Michaelis-Menten, Arrhenius, and diffusion functions.
Heterotrophic respiration (Rh), microbial processing of soil organic matter to carbon dioxide (CO2), is a major, yet highly uncertain, carbon (C) flux from terrestrial systems to the atmosphere. Temperature sensitivity of Rh is often represented with a simple Q10 function in ecosystem models and earth system models (ESMs), sometimes accompanied by an empirical soil moisture modifier. More explicit representation of the effects of soil moisture, substrate supply, and their interactions with temperature has been proposed as a way to disentangle the confounding factors of apparent temperature sensitivity of Rh and improve the performance of ecosystem models and ESMs.
Here, we merged the DAMM soil flux model with a parsimonious ecosystem flux model, FöBAAR (Forest Biomass, Assimilation, Allocation and Respiration). We used high-frequency soil flux data from automated soil chambers and landscape-scale ecosystem fluxes from eddy covariance towers at two AmeriFlux sites (Harvard Forest, MA and Howland Forest, ME) in the northeastern USA to estimate parameters, validate the merged model, and to quantify the uncertainties in a multiple constraints approach. The optimized DAMM-FöBAAR model better captured the seasonal and inter-annual dynamics of soil respiration compared to the FöBAAR-only model for the Harvard Forest. However, DAMM-FöBAAR showed improvement over FöBAAR-only at the boreal transition Howland Forest only in unusually dry years. The frequency of synoptic-scale dry periods is lower at Howland. At both sites, the declining trend of Rh during drying events was captured by the DAMM-FöBAAR model. While the DAMM functions require a few more parameters than a simple Q10 function, we demonstrate that they can be included in an ecosystem model and reduce model-data mismatches. Moreover, the mechanistic structure of the soil moisture effects using DAMM functions should be more generalizable than the wide variety of empirical functions that are commonly used, and these DAMM functions could be readily incorporated into other ecosystem models and ESMs.