Transient hysteresis of near-surface permafrost response to external forcing

Estimates of changes in near-surface permafrost (NSP) area S (p) relative to change in globally averaged surface air temperature T (g) are made by using the global climate model developed at the A.M. Obukhov Institute of Atmospheric Physics RAS (IAP RAS CM). For ensemble of runs forced by scenarios constructed as return-to-preindustrial continuations of the RCP (Representative Concentration Pathways) scenarios family, a possibility of transient hysteresis in dependence of S (p) versus T (g) is exhibited: in some temperature range which depends on imposed scenario of external forcing, NSP area is larger, at the same global mean surface air temperature, in a warming climate than in a cooling climate. This hysteresis is visible more clearly for scenarios with higher concentration of greenhouse gases in the atmosphere in comparison to those in which this concentration is lower. Hysteresis details are not sensitive to the type of the prescribed continuation path which is used to return the climate to the preindustrial state. The multiple-valued dependence of S (p) on T (g) arises due to dependence of soil state in the regions of extra-tropical wetlands and near the contemporary NSP boundaries on sign of external climatic forcing. To study the dependence of permafrost hysteresis on amplitude and temporal scale of external forcing, additional model runs are performed. These runs are forced by idealised scenarios of atmospheric CO2 content varying, depending on run, with periods from 100 to 1,000 year and with different amplitudes. It is shown that the above-mentioned hysteresis is related to the impact of phase transitions of soil water on apparent inertia of the system as well as to the impact of soil state on atmospheric hydrological cycle and radiation transfer in the atmosphere.