Remote sensing, regional ground temperature and ground ice observations, and numerical simulation were used to investigate the size, distribution, and activity of ice wedges in fine-grained mineral and organic soils across the forest-tundra transition in uplands east of the Mackenzie Delta. In the northernmost dwarf-shrub tundra, ice wedge polygons cover up to 40% of the ground surface, with the wedges commonly exceeding 3 m in width. The largest ice wedges are in peatlands where thermal contraction cracking occurs more frequently than in nearby hummocky terrain with fine-grained soils. There are fewer ice wedges, rarely exceeding 2 m in width, in uplands to the south and none have been found in mineral soils of the tall-shrub tundra, although active ice wedges are found there throughout peatlands. In the spruce forest zone, small, relict ice wedges are restricted to peatlands. At tundra sites, winter temperatures at the top of permafrost are lower in organic than mineral soils because of the shallow permafrost table, occurrence of phase change at 0°C, and the relatively high thermal conductivity of icy peat. Due to these factors and the high coefficient of thermal contraction of frozen saturated peat, ice wedge cracking and growth is more common in peatlands than in mineral soil. However, the high latent heat content of saturated organic active layer soils may inhibit freezeback, particularly where thick snow accumulates, making the permafrost and the ice wedges in spruce forest polygonal peatlands susceptible to degradation following alteration of drainage or climate warming.

No abstract

Owing to the mathematical similarities between heat and mass transport, the multi-species transport model MT3DMS should be able to simulate heat transport if the effects of buoyancy and changes in viscosity are small. Although in several studies solute models have been successfully applied to simulate heat transport, these studies failed to provide any rigorous test of this approach. In the current study, we carefully evaluate simulations of a single borehole ground source heat pump (GSHP) system in three scenarios: a pure conduction situation, an intermediate case, and a convection-dominated case. Two evaluation approaches are employed: first, MT3DMS heat transport results are compared with analytical solutions. Second, simulations by MT3DMS, which is finite difference, are compared with those by the finite element code FEFLOW and the finite difference code SEAWAT. Both FEFLOW and SEAWAT are designed to simulate heat flow. For each comparison, the computed results are examined based on residual errors. MT3DMS and the analytical solutions compare satisfactorily. MT3DMS and SEAWAT results show very good agreement for all cases. MT3DMS and FEFLOW two-dimensional (2D) and three-dimensional (3D) results show good to very good agreement, except that in 3D there is somewhat deteriorated agreement close to the heat source where the difference in numerical methods is thought to influence the solution. The results suggest that MT3DMS can be successfully applied to simulate GSHP systems, and likely other systems with similar temperature ranges and gradients in saturated porous media.

No abstract

Infiltrating river water carries the temperature signal of the river into the adjacent aquifer. While the diurnal temperature fluctuations are strongly dampened, the seasonal fluctuations are much less attenuated and can be followed into the aquifer over longer distances. In one-dimensional model with uniform properties, this signal is propagated with a retarded velocity, and its amplitude decreases exponentially with distance. Therefore, time shifts in seasonal temperature signals between rivers and groundwater observation points may be used to estimate infiltration rates and near-river groundwater velocities. As demonstrated in this study, however, the interpretation is nonunique under realistic conditions. We analyze a synthetic test case of a two-dimensional cross section perpendicular to a losing stream, accounting for multi-dimensional flow due to a partially penetrating channel, convective-conductive heat transport within the aquifer, and heat exchange with the underlying aquitard and the land surface. We compare different conceptual simplifications of the domain in order to elaborate on the importance of different system elements. We find that temperature propagation within the shallow aquifer can be highly influenced by conduction through the unsaturated zone and into the underlying aquitard. In contrast, regional groundwater recharge has no major effect on the simulated results. In our setup, multi-dimensionality of the flow field is important only close to the river. We conclude that over-simplistic analytical models can introduce substantial errors if vertical heat exchange at the aquifer boundaries is not accounted for. This has to be considered when using seasonal temperature fluctuations as a natural tracer for bank infiltration.