Mountain water resources are of particular importance for downstream populations but are threatened by decreasing water storage in snowpack and glaciers. Groundwater contribution to mountain ...streamflow, once assumed to be relatively small, is now understood to represent an important water source to streams. This review presents an overview of research on groundwater in high mountain environments (As classified by Meybeck et al. (2001) as very high, high, and mid‐altitude mountains). Coarse geomorphic units, like talus, alluvium, and moraines, are important stores and conduits for high mountain groundwater. Bedrock aquifers contribute to catchment streamflow through shallow, weathered bedrock but also to higher order streams and central valley aquifers through deep fracture flow and mountain‐block recharge. Tracer and water balance studies have shown that groundwater contributes substantially to streamflow in many high mountain catchments, particularly during low‐flow periods. The percentage of streamflow attributable to groundwater varies greatly through time and between watersheds depending on the geology, topography, climate, and spatial scale. Recharge to high mountain aquifers is spatially variable and comes from a combination of infiltration from rain, snowmelt, and glacier melt, as well as concentrated recharge beneath losing streams, or through fractures and swallow holes. Recent advances suggest that high mountain groundwater may provide some resilience—at least temporarily—to climate‐driven glacier and snowpack recession. A paucity of field data and the heterogeneity of alpine landscapes remain important challenges, but new data sources, tracers, and modeling methods continue to expand our understanding of high mountain groundwater flow.
This article is categorized under:
Science of Water > Hydrological Processes
Science of Water > Water and Environmental Change
Science of Water > Methods
Conceptual model of high mountain hydrogeological processes. Water is stored below the ground surface in talus slopes, moraines, alluvial deposits, wetlands, and bedrock aquifers. Groundwater flow is impacted by geologic heterogeneity and structures, and permafrost as it drains to springs, mountain rivers, wetlands, or contributes to mountain‐block recharge, helping to sustain mountain streamflow or basin water resources, particularly during dry periods.
Climate change is expected to increase regional and global air temperatures and significantly alter precipitation regimes. These projected changes in meteorological conditions will likely influence ...subsurface thermal regimes. Increases in groundwater and soil temperatures could impact groundwater quality, harm groundwater-sourced ecosystems, and contribute to the geotechnical failure of critical infrastructure. Furthermore, permafrost thaw induced by rising subsurface temperatures will likely alter surface and subsurface hydrology in high altitude and/or latitude regions and exacerbate the rate of anthropogenic climate change by releasing stored carbon into the atmosphere.
This contribution discusses the theory and development of subsurface heat transport equations for cold and temperate regions. Analytical solutions to transient forms of the conduction equation and the conduction–advection equation with and without freezing are detailed. In addition, recently developed groundwater flow and heat transport models that can accommodate freezing and thawing processes are briefly summarized. These models can be applied to simulate climate change-induced permafrost degradation and dormant aquifer activation in cold regions.
Several previous reviews have focused on the impact of climate change on subsurface hydraulic regimes and groundwater resources, but this is the first synthesis of studies considering the influence of future climate change on subsurface thermal regimes in cold and temperate regions. The current gaps in this body of knowledge are highlighted, and recommendations are made for improving future studies by linking atmospheric global climate models to subsurface heat transport models that consider heat advection via groundwater flow.
Permafrost thaw due to climate warming modifies hydrological processes by increasing hydrological connectivity between aquifers and surface water bodies and increasing groundwater storage. While ...previous studies have documented arctic river baseflow increases and changing wetland and lake distributions, the hydrogeological processes leading to these changes remain poorly understood. This study uses a coupled heat and groundwater flow numerical model with dynamic freezing and thawing processes and an improved set of boundary conditions to simulate the impacts of climate warming on permafrost distribution and groundwater discharge to surface water bodies. We show a spatial shift in groundwater discharge from upslope to downslope and a temporal shift with increasing groundwater discharge during the winter season due to the formation of a lateral supra-permafrost talik underlying the active layer. These insights into changing patterns of groundwater discharge help explain observed changes in arctic baseflow and wetland patterns and are important for northern water resources and ecosystem management.
Recent climate change has reduced the spatial extent and thickness of permafrost in many discontinuous permafrost regions. Rapid permafrost thaw is producing distinct landscape changes in the Taiga ...Plains of the Northwest Territories, Canada. As permafrost bodies underlying forested peat plateaus shrink, the landscape slowly transitions into unforested wetlands. The expansion of wetlands has enhanced the hydrologic connectivity of many watersheds via new surface and near‐surface flow paths, and increased streamflow has been observed. Furthermore, the decrease in forested peat plateaus results in a net loss of boreal forest and associated ecosystems. This study investigates fundamental processes that contribute to permafrost thaw by comparing observed and simulated thaw development and landscape transition of a peat plateau‐wetland complex in the Northwest Territories, Canada from 1970 to 2012. Measured climate data are first used to drive surface energy balance simulations for the wetland and peat plateau. Near‐surface soil temperatures simulated in the surface energy balance model are then applied as the upper boundary condition to a three‐dimensional model of subsurface water flow and coupled energy transport with freeze‐thaw. Simulation results demonstrate that lateral heat transfer, which is not considered in many permafrost models, can influence permafrost thaw rates. Furthermore, the simulations indicate that landscape evolution arising from permafrost thaw acts as a positive feedback mechanism that increases the energy absorbed at the land surface and produces additional permafrost thaw. The modeling results also demonstrate that flow rates in local groundwater flow systems may be enhanced by the degradation of isolated permafrost bodies.
Key Points:
Observed permafrost thaw rates are compared to results from a 3‐D groundwater flow and heat transfer model
Lateral heat flow can accelerate discontinuous permafrost thaw and land cover change in peatlands
Degradation of discontinuous permafrost enhances local groundwater flow
In northern regions, transportation infrastructure can experience severe structural damages due to permafrost degradation. Water infiltration and subsurface water flow under an embankment affect the ...energy balance of roadways and underlying permafrost. However, the quantification of the processes controlling these changes and a detailed investigation of their thermal impacts remain largely unknown due to a lack of available long‐term embankment temperature data in permafrost regions. Here, we report observations of heat advection linked to surface water infiltration and subsurface flow based on a 9‐year (from 2009 to 2017) thermal monitoring at an experimental road test site built on ice‐rich permafrost conditions in southwestern Yukon, Canada. Our results show that snowmelt water infiltration in the spring rapidly increases temperature in the upper portion of the embankment. The earlier disappearance of snow deposited at the embankment slope increases the thawing period and the temperature gradient in the embankment compared with the natural ground. Infiltrated summer rainfall water lowered the near‐surface temperatures and subsequently warmed embankment fill materials down to 3.6‐m depth. Heat advection caused by the flow of subsurface water produced warming rates at depth in the embankment subgrade up to two orders of magnitude faster than by atmospheric warming (heat conduction). Subsurface water flow promoted permafrost thawing under the road embankment and led to an increase in active layer thickness. We conclude that the thermal stability of roadways along the Alaska Highway corridor is not maintainable in situations where water is flowing under the infrastructure unless mitigation techniques are used. Severe structural damages to the highway embankment are expected to occur in the next decade.
Our manuscript presents three important novel findings:
The infiltration of spring snowmelt causes a very rapid temperature increase in the road embankment, within just a few days.
The infiltration of summer rainfall lowers near‐surface temperatures while simultaneously warming embankment fill materials at depth.
The advection of heat due to porewater flow through subsurface water flow leads to warming rates that can be up to two orders of magnitude faster than by just atmospheric warming.
► We present a method for calculating pore-water flux from temperature time series. ► The new method uses advancements in signal processing to analyze raw field data. ► Many temperature sensor pairs ...in profiles can be analyzed to find ideal spacing. ► Field examples include high-resolution DTS profile in shallow hyporheic flow path. ► Method is incorporated into freely-available computer program named VFLUX.
Heat is a useful tracer for quantifying groundwater–surface water interaction, but analyzing large amounts of raw thermal data has many challenges. We present a computer program named VFLUX, written in the MATLAB computing language, for processing raw temperature time series and calculating vertical water flux in shallow sub-surface-water systems. The step-by-step workflow synthesizes several recent advancements in signal processing, and adds new techniques for calculating flux rates with large numbers of temperature records from high-resolution sensor profiles. The program includes functions for quantitatively evaluating the ideal spacing between sensor pairs, and for performing error and sensitivity analyses for the heat transport model due to thermal parameter uncertainty. The program synchronizes and resamples temperature data from multiple sensors in a vertical profile, isolates the diurnal signal from each time series and extracts its amplitude and phase angle information using Dynamic Harmonic Regression (DHR), and calculates vertical water flux rates between multiple sensor pairs using heat transport models. Flux rates are calculated every 1-to-2
h using four similar analytical methods. One or more “sliding analysis windows” can be used to automatically identify any number of variably spaced sensor pairs for flux calculations, which is necessary when a single vertical profile contains many sensors, such as in a high-resolution fiber-optic distributed temperature sensing (DTS) profile. We demonstrate the new method by processing two field temperature time series datasets collected using discrete temperature sensors and a high-resolution DTS profile. The analyses of field data show vertical flux rates significantly decreasing with depth at high-spatial resolution as the sensor profiles penetrate shallow, curved hyporheic flow paths, patterns which may have been obscured without the unique analytical abilities of VFLUX.
The tropical glaciers of the Cordillera Blanca, Peru, are rapidly retreating, resulting in complex impacts on the hydrology of the upper Río Santa watershed. The effect of this retreat on water ...resources is evaluated by analyzing historical and recent time series of daily discharge at nine measurement points. Using the Mann-Kendall nonparametric statistical test, the significance of trends in three hydrograph parameters was studied. Results are interpreted using synthetic time series generated from a hydrologic model that calculates hydrographs based on glacier retreat sequences. The results suggest that seven of the nine study watersheds have probably crossed a critical transition point, and now exhibit decreasing dry-season discharge. Our results suggest also that once the glaciers completely melt, annual discharge will be lower than present by 2-30% depending on the watershed. The retreat influence on discharge will be more pronounced during the dry season than at other periods of the year. At La Balsa, which measures discharge from the upper Río Santa, the glacier retreat could lead to a decrease in dry-season average discharge of 30%.
Groundwater flow in cold regions containing permafrost accelerates climate-warming-driven thaw and changes thaw patterns. Simulation analyses of groundwater flow and heat transport with freeze/thaw ...in typical cold-regions terrain with nested flow indicate that early thaw rate is particularly enhanced by flow, the time when adverse environmental impacts of climate-warming-induced permafrost loss may be severest. For the slowest climate-warming rate predicted by the Intergovernmental Panel on Climate Change (IPCC), once significant groundwater flow begins, thick permafrost layers can vanish in several hundred years, but survive over 1,000 years where flow is minimal. Large-scale thaw depends mostly on the balance of heat advection and conduction in the supra-permafrost zone. Surface-water bodies underlain by open taliks allow slow sub-permafrost flow, with lesser influence on regional thaw. Advection dominance over conduction depends on permeability and topography. Groundwater flow around permafrost and flow through permafrost impact thaw differently; the latter enhances early thaw rate. Air-temperature seasonality also increases early thaw. Hydrogeologic heterogeneity and topography strongly affect thaw rates/patterns. Permafrost controls the groundwater/surface-water-geomorphology system; hence, prediction and mitigation of impacts of thaw on ecology, chemical exports and infrastructure require improved hydrogeology/permafrost characterization and understanding.
Surface energy balance (SEB) strongly influences the thermal state of permafrost, cryohydrological processes, and infrastructure stability. Road construction and snow accumulation affect the energy ...balance of underlying permafrost. Herein, we use an experimental road section of the Alaska Highway to develop a SEB model to quantify the surface energy components and ground surface temperature (GST) for different land cover types with varying snow regimes and properties. Simulated and measured ground temperatures are in good agreement, and our results show that the quantity of heat entering the embankment center and slope is mainly controlled by net radiation, and less by the sensible heat flux. In spring, lateral heat flux from the embankment center leads to earlier disappearance of snowpack on the embankment slope. In winter, the insulation created by the snow cover on the embankment slope reduces heat loss by a factor of three compared with the embankment center where the snow is plowed. The surface temperature offsets are 5.0°C and 7.8°C for the embankment center and slope, respectively. Furthermore, the heat flux released on the embankment slope exponentially decreases with increasing snow depth, and linearly decreases with earlier snow cover in fall and shorter snow‐covered period in spring.
Actively heated fiber-optic distributed temperature sensing (aFO-DTS) measures soil moisture content at sub-meter intervals across kilometres of fiber-optic cable. The technology has great potential ...for environmental monitoring but calibration at field scales with variable soil conditions is challenging. To better understand and quantify the errors associated with aFO-DTS soil moisture measurements, we use a parametric numerical modeling approach to evaluate different error factors for uniform soil. A thermo-hydrogeologic, unsaturated numerical model is used to simulate a 0.01 m by 0.01 m two-dimensional domain, including soil and a fiber-optic cable. Results from the model are compared to soil moisture values calculated using the commonly used Tcum calibration method for aFO-DTS. The model is found to have high accuracy between measured and observed saturations for static hydrologic conditions but shows discrepancies for more realistic settings with active recharge. We evaluate the performance of aFO-DTS soil moisture calculations for various scenarios, including varying recharge duration and heterogeneous soils. The aFO-DTS accuracy decreases as the variability in soil properties and intensity of recharge events increases. Further, we show that the burial of the fiber-optic cable within soil may adversely affect calculated results. The results demonstrate the need for careful selection of calibration data for this emerging method of measuring soil moisture content.