Coupled multilayer canopy-permafrost model (CryoGrid) for the simulation of forest trajectories in permafrost underlain boreal forests.

• Main Authors: Stuenzi, Simone Maria, Boike, Julia, Gädeke, Anne, Herzschuh, Ulrike, Kruse, Stefan, Pestryakova, Luidmila A., Westermann, Sebastian, Langer, Moritz
• Format: Software
• Language: English
• Published: 2021
• Subjects: Permafrost; Land surface model; CryoGrid; Boreal forest; Zakharov; Spasskaya; permafrost; Siberia
• Online Access: https://zenodo.org/record/4603668; https://doi.org/10.5281/zenodo.4603668

CryoGrid is a land-surface scheme dedicated to modeling of ground temperatures in permafrost environments. Here, the one-dimensional land surface model (CryoGrid) is adapted for application in vegetated areas by coupling a multilayer canopy model (CLM-ml v0). This model setup is used to reproduce the energy transfer and thermal regime at a study site in mixed boreal forest in Eastern Siberia. The vegetation module forms the upper boundary layer of the coupled vegetation-permafrost model and replaces the surface energy balance equation used for common CryoGrid representations. The coupled model was first described in the following article which has been published in Biogeosciences: Stuenzi, S. M., Boike, J., Cable, W., Herzschuh, U., Kruse, S., Pestryakova, L. A., Schneider von Deimling, T., Westermann, S., Zakharov, E. S., and Langer, M.: Variability of the surface energy balance in permafrost-underlain boreal forest, Biogeosciences, 18, 343–365, https://doi.org/10.5194/bg-18-343-2021, 2021. The model code for this publication can be found here: https://doi.org/10.5281/zenodo.4317106 Here, we add a paramterization for deciduous forest to simulate the leafless state of deciduous-dominated regions outside of the short vegetative period in summer. This is achieved by allowing for a separate leaf area index defined by a rough parameterization of a leaf-on and a leaf-off season (10. October - 10. April) based on literature values from Spasskaya Pad. Further, more realistic mixed canopy compositions can now be simulated by allowing for a certain percentage of deciduous taxa within the simulated forest stand. In addition, we add a parameterization for coupling forest density (LAI) to fine root biomass. Further, we have implemented a new relationship for phase partitioning of water in frozen soil (freeze curve). The parameters are set to the default values that were used for the simulations in the article. Parameters different from the default values can be specified in the main script run_CG_RUN_1D_STANDARD.m (general ...