Modeling the thermal dynamics of the active layer at two contrasting permafrost sites on Svalbard and on the Tibetan Plateau

Abstract. Employing a one-dimensional, coupled thermal and hydraulic numerical model, we quantitatively analyze high-resolution, multi-year data from the active layers at two contrasting permafrost sites. The model implements heat conduction with the de Vries parameterization, heat convection with water and vapor flow, freeze-thaw transition parameterized with a heuristic soil-freezing characteristic, and liquid water flow with the Mualem-van Genuchten parameterization. The model is driven by measured temperatures and water contents at the upper and lower boundary with all required material properties deduced from the measured data. The aims are (i) to ascertain the applicability of such a rather simple model, (ii) to quantify the dominating processes, and (iii) to discuss possible causes of remaining deviations. Soil temperatures and water contents as well as characteristic quantities like thaw depth and duration of the isothermal plateau could be reproduced. Heat conduction is found to be the dominant process by far at both sites, with non-conductive transport contributing a maximum of some 3% to the mean heat flux at the Spitsbergen site, most of the time very much less, and practically negligible at the Tibetan site. Hypotheses discussed to explain the remaining deviations between measured and simulated state variables include, besides some technical issues, infiltration of snow melt, dry cracking with associated vapor condensation, and mechanical soil expansion in detail.

[1]  G. L. Guymon,et al.  A coupled heat and moisture transport model for Arctic soils , 1974 .

[2]  Howard E. Epstein,et al.  Active‐Layer Hydrology in Nonsorted Circle Ecosystems of the Arctic Tundra , 2007 .

[3]  Douglas L. Kane,et al.  Thermal response of the active layer to climatic warming in a permafrost environment , 1991 .

[4]  T. Koike,et al.  Frozen soil parameterization in a distributed biosphere hydrological model , 2009 .

[5]  R. Harlan,et al.  Analysis of coupled heat‐fluid transport in partially frozen soil , 1973 .

[6]  L. Hinzman,et al.  A distributed thermal model for calculating soil temperature profiles and depth of thaw in permafrost regions , 1998 .

[7]  P. Deuflhard.  A modified Newton method for the solution of ill-conditioned systems of nonlinear equations with application to multiple shooting , 1974 .

[8]  Arden L. Buck,et al.  New Equations for Computing Vapor Pressure and Enhancement Factor , 1981 .

[9]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[10]  J. C. Jaeger,et al.  Conduction of Heat in Solids , 1952 .

[11]  F. Vermolen,et al.  A comparison of numerical models for one-dimensional Stefan problems , 2006 .

[12]  Timothy A. Davis,et al.  Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method , 2004, TOMS.

[13]  Richard L. Snyder,et al.  Coupling of heat, water vapor, and liquid water fluxes to compute evaporation in bare soils , 2008 .

[14]  F. Gasse,et al.  A 13,000-year climate record from western Tibet , 1991, Nature.

[15]  Yoshisuke Nakano,et al.  Effect of a Freezing Zone of Finite Width on the Thermal Regime of Soils , 1971 .

[16]  J. Boike,et al.  Seasonal snow cover on frozen ground: Energy balance calculations of a permafrost site near Ny‐Ålesund, Spitsbergen , 2003 .

[17]  Bernd Etzelmüller,et al.  Recent advances in permafrost modelling , 2008 .

[18]  Paul R. Lowe,et al.  The Computation of Saturation Vapor Pressure , 1974 .

[19]  N. Afgan,et al.  Transfer processes in the plant environment , 1975 .

[20]  L. Hinzman,et al.  Non-conductive heat transfer associated with frozen soils , 2001 .

[21]  J. Philip.  EVAPORATION, AND MOISTURE AND HEAT FIELDS IN THE SOIL , 1957 .

[22]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[23]  Vladimir E. Romanovsky,et al.  Thawing of the Active Layer on the Coastal Plain of the Alaskan Arctic , 1997 .

[24]  J. Boike,et al.  Water, heat and solute dynamics of a mud boil, Spitsbergen , 2008 .

[25]  J. Boike,et al.  Quantifying the thermal dynamics of a permafrost site near Ny‐Ålesund, Svalbard , 2001 .

[26]  S. Tulaczyk,et al.  Response of subglacial sediments to basal freeze-on 1. Theory and comparison to observations from beneath the West Antarctic Ice Sheet , 2003 .

[27]  A. G. Agwu Nnanna,et al.  Phase-change phenomena in porous media – a non-local thermal equilibrium model , 2001 .

[28]  Helén Engelmark,et al.  Numerical modelling of phase change in freezing and thawing unsaturated soil , 1993 .

[29]  Masaru Mizoguchi,et al.  Water Flow and Heat Transport in Frozen Soil , 2004 .

[30]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..