Abstract. Arctic and sub-arctic regions are amongst the most susceptible regions on Earth to global warming and climate change. Understanding and predicting the impact of climate change in these regions require a proper process representation of the interactions between climate, the carbon cycle, and hydrology in Earth system models. This study focuses on Land Surface Models (LSMs) that represent the lower boundary condition of General Circulation Models (GCMs) and Regional Climate Models (RCMs), which simulate climate change evolution at the global and regional scales, respectively. LSMs 15 typically utilize a standard soil configuration with a depth of no more than 4 meters, whereas for cold, permafrost regions, field experiments show that attention to deep soil profiles is needed to understand and close the water and energy balances, which are tightly coupled through the phase change. To address this, we design and run a series of model experiments with a one-dimensional LSM, called CLASS (Canadian Land Surface Scheme), as embedded in the MESH (Modélisation Environmentale Communautaire -Surface and Hydrology) modelling system, to (1) characterize the effect of soil profile 20 depth under different climate conditions and in the presence of parameter uncertainty, and (2) develop a methodology for temperature profile initialization in permafrost regions, where the system has an extended memory, by the use of paleorecords and bootstrapping. Our study area is in Norman Wells, Northwest Territories of Canada, where measurements of soil temperature profiles and historical reconstructed climate data are available. Our results demonstrate that the adequate depth of soil profile in an LSM varies for warmer and colder conditions and is sensitive to model parameters and the uncertainty 25 around them. In general, however, we show that a minimum of 20 meters of soil profile is essential to adequately represent the temperature dynamics. Our results also indicate the significance of model initialization in permafrost regions and our proposed spin-up method requires running the LSM over more than 300 years of reconstructed climate time series.Hydrol. Earth Syst. Sci. Discuss., https://doi