Climate warming has caused increased air temperature as well as increased subsurface temperature. Many previous studies on subsurface warming simplified heat advection by neglecting the horizontal component of regional groundwater flow or even neglected heat advection accompanying groundwater flow. In this study, the simultaneous control of heat advection and conduction on subsurface warming is numerically investigated in a 2D hypothetical basin cross section. By calculating the increment of subsurface temperature, we find that heat advection could accelerate subsurface warming. In a given basin, subsurface warming in the recharge area with downward groundwater flow is more significant than that in the discharge area with upward groundwater flow. By using a 1D model with vertical groundwater flow only for comparison, we find that in any part of a basin, even if the horizontal flux is very small compared with the vertical flux, neglecting the horizontal flux would underestimate the propagation depth of climate warming. This implies that when the propagation depth of climate warming is a priori known variable, existing 1D models would overestimate downward groundwater flux or underestimate upward groundwater flux. Moreover, we find pumping would cause deeper propagation depths of climate warming by accelerating groundwater circulation, whereas basin-scale heterogeneity and anisotropy of hydraulic conductivity would cause shallower propagation depths of climate warming because of the relative dominance of horizontal flow. By demonstrating the importance of 2D groundwater flow on subsurface warming, the results provide new insight into understanding the regulation of temperature among the atmosphere, the hydrosphere and the lithosphere.

Rising temperatures in the Arctic and subarctic are driving the rapid thaw of permafrost by reducing permafrost cooling, increasing active layer thickness, and promoting talik formation. In this study, the cyrohydrogeology of a permafrost mound located within the discontinuous permafrost zone near Umiujaq (Nunavik, Quebec, Canada) is characterized through the analysis of a dataset covering more than two decades of monitoring. This dataset captures a high degree of interannual variability in air temperature and ground thermal conditions, as well as the formation and closure of a supra-permafrost talik. Data indicate that variable saturation and advective heat transport directly contribute to the expansion and contraction of the talik. Data further indicate the presence of two distinct thermo-hydrologic settings resulting from differences in surface conditions, as well as subsurface thermal and flow regimes. The first, found at the top of the mound feature, is characterized by very low moisture contents (& lt;0.05 m(3)/m(3)), while the second, found at the side of the mound feature, shows higher annual moisture contents that strongly influence the dynamics of heat and groundwater flow. The data were synthesized into a detailed conceptual model of the cyrohydrogeological dynamics that highlights the important role of hydrogeological characterization and long-term data sets in understanding the effects of groundwater flow on seasonal frost and permafrost dynamics. Specifically, the results presented here show that in the absence of long-term data sets, longer-period transient phenomena such as talik opening and closure may be misrepresented as uni-directional feedback loops, as opposed to highly dynamic temporary phenomena.

Advective heat transported by water percolating into discontinuities in frozen ground can rapidly increase temperatures at depth because it provides a thermal shortcut between the atmosphere and the subsurface. Here, we develop a conceptual model that incorporates the main heat-exchange processes in a rock cleft. Laboratory experiments and numerical simulations based on the model indicate that latent heat release due to initial ice aggradation can rapidly warm cold bedrock and precondition it for later thermal erosion of cleft ice by advected sensible heat. The timing and duration of water percolation both affect the ice-level change if initial aggradation and subsequent erosion are of the same order of magnitude. The surplus advected heat is absorbed by cleft ice loss and runoff from the cleft so that this energy is not directly detectable in ground temperature records. Our findings suggest that thawing-related rockfall is possible even in cold permafrost if meltwater production and flow characteristics change significantly. Advective warming could rapidly affect failure planes beneath large rock masses and failure events could therefore differ greatly from common magnitude reaction-time relations. Copyright (C) 2011 John Wiley & Sons, Ltd.