Geology topics

Jeffrey M. McKenzie

Publications and source records attributed to Jeffrey M. McKenzie.

11 recordsLinked to original sources

USGS2024

SUTRA— A code for simulation of saturated-unsaturated, variable-density groundwater flow with solute or energy transport—Documentation of the version 4.0 enhancements—Freeze-thaw capability, saturation and relative-permeability relations, spatially varying properties, and enhanced budget and velocity outputs

Version 4.0 of the Saturated-Unsaturated Transport (SUTRA) software code provides the capability to simulate the freezing and thawing of groundwater during energy transport simulations under saturated and unsaturated conditions. In addition to the types of hydrogeologic processes that SUTRA has been able to simulate in the past, this version can be used to study the effects of the freeze-thaw process on the flow and energy dynamics of hydrogeologic systems. The freeze-thaw simulation capability accounts for the latent heat of fusion and allows thermal property values to vary with changing total-water saturation, liquid-water saturation, and ice saturation. It allows the effective permeability of the porous medium to change as a result of freezing and thawing. This version also provides several user-selectable relations for the dependence of total-water saturation on fluid pressure, the dependence of liquid-water saturation on temperature during freezing and thawing, and the dependence of relative permeability on liquid saturation, as well as three user-selectable formulae for defining the bulk thermal conductivity of a mixture of solid grains, liquid water, ice, and air. For unsaturated simulations without freezing, the selectable total-water saturation relations eliminate the need for the user to program these and their associated relative-permeability functions, as had been required in previous SUTRA versions. Optional nonlinear dependence of fluid density on temperature, which covers the range from supercooled (about −50 degrees Celsius) to superheated (about 400 degrees Celsius), is also provided. Additionally, this version makes it possible to spatially vary parameters that, in previous versions of SUTRA, were required to be spatially uniform: solid-matrix properties, adsorption parameters, and parameters for production of solute mass or energy. Spatial variation is also allowed for the newly included freeze-thaw process parameters. Additional enhancements provide (1) output of water-mass and energy budgets that include values of all component terms in the governing balance equations, and (2) output of Darcy velocities (fluid fluxes), in addition to the velocity output provided by previous SUTRA versions. These enhanced outputs allow fuller interpretation of simulation results, especially for freeze-thaw phenomena. The set of processes simulated by this version of SUTRA are useful for studying a wide range of hydrogeologic system types, conditions, and questions. For cryohydrogeologic simulations, however, this version of the code is limited in that (1) it does not simulate thermomechanical effects of freeze-thaw, (2) pressure changes due to water density change during freezing are neglected, (3) ice saturation cannot exceed the initial porosity of the simulated medium, and (4) cryosuction, the migration of liquid water toward freezing fronts, is neglected. Furthermore, this version does not account for air flow or for water vaporization and sublimation under unsaturated conditions.

Techniques and Methods↗

USGS2024

A history of cryohydrogeology modeling and recent advancements through the integration of solute transport

Groundwater flow systems and permafrost are interrelated because permafrost thaw enhances permeability, while groundwater flow can advect heat and accelerate permafrost thaw (McKenzie et al. 2021). Given amplified climate change in cold regions, there is renewed interest in ‘cryohydrogeology’, the study of groundwater in cold regions. Many data-driven studies have shown that permafrost thaw is leading to activated aquifers and increased baseflow across the pan-Arctic region (e.g. Walvoord and Striegl 2007, Evans et al. 2020). Empirical evidence of a subsurface ‘replumbing’ (Walvoord and Kurylyk 2016) in permafrost regions raises questions about the fate of sequestered contaminants in the North (Langer et al. 2023). We will discuss the history of and emerging opportunities in cryohydrogeological modeling, with a focus on recent contaminant transport modeling.

Conference Paper↗

USGS2023

Subsurface porewater flow accelerates talik development under the Alaska Highway, Yukon: A prelude to road collapse and permafrost thaw?

The presence of taliks (perennially unfrozen zones in permafrost areas) adversely affects the thermal stability of infrastructure in cold regions, including roads. The role of heat advection on talik development and feedback on permafrost degradation has not been quantified methodically in this context. We incorporate a surface energy balance model into a coupled groundwater flow and energy transport numerical model (SUTRA-ice). The model, calibrated with long-term observations (1997–2018 on the Alaska Highway), is used to investigate and quantify the role of heat advection on talik initiation and development under a road embankment. Over the 25-year simulation period, the new model is driven by reconstructed meteorological data and has a good agreement with near surface soil temperatures. The model successfully reproduces the increasing depth to the permafrost table (mean absolute error <0.2 m), and talik development. The results demonstrate that heat advection provides an additional energy source that expedites the rate of permafrost thaw and roughly doubles the rate of permafrost table deepening, compared to purely conductive thawing. Talik initially formed and grew over time under the combined effect of water flow, snow insulation, road construction and climate warming. Talik formation creates a new thermal state under the road embankment, resulting in acceleration of underlying permafrost degradation, due to the positive feedback of heat accumulation created by trapped unfrozen water. In a changing climate, mobile water flow will play a more important role in permafrost thaw and talik development under road embankments, and is likely to significantly increase maintenance costs and reduce the long-term stability of the infrastructure.

Yukon↗

USGS2021

Surface energy balance of sub-Arctic roads with varying snow regimes and properties in permafrost regions

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.

Yukon↗

USGS2021

Invited perspective: What lies beneath a changing Arctic?

As permafrost thaws in the Arctic, new subsurface pathways open for the transport of groundwater, energy, and solutes. We identify different ways that these subsurface changes are driving observed surface consequences, including the potential for increased contaminant transport, modification to water resources, and enhanced rates of infrastructure (e.g. buildings and roads) damage. Further, as permafrost thaws it allows groundwater to transport carbon, nutrients, and other dissolved constituents from terrestrial to aquatic environments via progressively deeper subsurface flow paths. Cryohydrogeology, the study of groundwater in cold regions, should be included in northern research initiatives to account for this hidden catalyst of environmental and societal change.

The Cryosphere↗

USGS2019

Hydrogeology of an alpine talus aquifer: Cordillera Blanca, Peru

The dramatic loss of glacial mass in low latitudes is causing shifts in downstream water availability and use during the driest months of the year. The world’s largest concentration of tropical glaciers lies in the Cordillera Blanca range of Peru, where glacial runoff is declining and regional stresses are emerging over water resources. Throughout the Cordillera Blanca, groundwater inputs from alpine meadow–talus complexes, locally known as pampas, supply proglacial streams with up to 80% of their flow during the region’s dry season. Structural knowledge of the pampa aquifers is needed to estimate their drainable groundwater storage capacity and residence time, to elucidate the role and importance of alpine groundwater storage in the regional water budget of the Cordillera Blanca. To understand the structure of these proglacial aquifers, multiple near-surface geophysical methods were implemented in a proglacial valley near dense networks of spring-fed tributaries. Geophysical results and borehole logs suggest groundwater is stored in a confined aquifer composed of buried talus deposits overlain by lacustrine clay, while deeper portions of the unit, 10–15 m in depth, are relatively clay-free and more hydraulically conductive. Based on these findings and assumptions of aquifer porosity, the pampas of the Callejon de Huaylas may store from 0.006 to 0.02 km 3 of groundwater. Furthermore, these findings suggest that the talus aquifers of the Cordillera Blanca were formed in proglacial lakes, followed by infilling with fine lacustrine sediments that confine lower units and allow for groundwater discharge to springs via macropores and preferential flow.

Hydrogeology Journal↗

USGS2018

GSFLOW-GRASS v1.0.0: GIS-enabled hydrologic modeling of coupled groundwater–surface-water systems

The importance of water moving between the atmosphere and aquifers has led to efforts to develop and maintain coupled models of surface water and groundwater. However, developing inputs to these models is usually time-consuming and requires extensive knowledge of software engineering, often prohibiting their use by many researchers and water managers, thus reducing these models' potential to promote science-driven decision-making in an era of global change and increasing water resource stress. In response to this need, we have developed GSFLOW–GRASS, a bundled set of open-source tools that develops inputs for, executes, and graphically displays the results of GSFLOW, the U.S. Geological Survey's coupled groundwater and surface-water flow model. In order to create a robust tool that can be widely implemented over diverse hydro(geo)logic settings, we built a series of GRASS GIS extensions that automatically discretizes a topological surface-water flow network that is linked with an underlying gridded groundwater domain. As inputs, GSFLOW–GRASS requires at a minimum a digital elevation model, a precipitation and temperature record, and estimates of channel parameters and hydraulic conductivity. We demonstrate the broad applicability of the toolbox by successfully testing it in environments with varying degrees of drainage integration, landscape relief, and grid resolution, as well as the presence of irregular coastal boundaries. These examples also show how GSFLOW–GRASS can be implemented to examine the role of groundwater–surface-water interactions in a diverse range of water resource and land management applications.

Geoscientific Model Development↗

USGS2018

Groundwater flow and heat transport for systems undergoing freeze-thaw: Intercomparison of numerical simulators for 2D test cases

In high-elevation, boreal and arctic regions, hydrological processes and associated water bodies can be strongly influenced by the distribution of permafrost. Recent field and modeling studies indicate that a fully-coupled multidimensional thermo-hydraulic approach is required to accurately model the evolution of these permafrost-impacted landscapes and groundwater systems. However, the relatively new and complex numerical codes being developed for coupled non-linear freeze-thaw systems require verification. This issue is addressed by means of an intercomparison of thirteen numerical codes for two-dimensional test cases with several performance metrics (PMs). These codes comprise a wide range of numerical approaches, spatial and temporal discretization strategies, and computational efficiencies. Results suggest that the codes provide robust results for the test cases considered and that minor discrepancies are explained by computational precision. However, larger discrepancies are observed for some PMs resulting from differences in the governing equations, discretization issues, or in the freezing curve used by some codes.

Advances in Water Resources↗

USGS2016

Using diurnal temperature signals to infer vertical groundwater-surface water exchange

Heat is a powerful tracer to quantify fluid exchange between surface water and groundwater. Temperature time series can be used to estimate pore water fluid flux, and techniques can be employed to extend these estimates to produce detailed plan-view flux maps. Key advantages of heat tracing include cost-effective sensors and ease of data collection and interpretation, without the need for expensive and time-consuming laboratory analyses or induced tracers. While the collection of temperature data in saturated sediments is relatively straightforward, several factors influence the reliability of flux estimates that are based on time series analysis (diurnal signals) of recorded temperatures. Sensor resolution and deployment are particularly important in obtaining robust flux estimates in upwelling conditions. Also, processing temperature time series data involves a sequence of complex steps, including filtering temperature signals, selection of appropriate thermal parameters, and selection of the optimal analytical solution for modeling. This review provides a synthesis of heat tracing using diurnal temperature oscillations, including details on optimal sensor selection and deployment, data processing, model parameterization, and an overview of computing tools available. Recent advances in diurnal temperature methods also provide the opportunity to determine local saturated thermal diffusivity, which can improve the accuracy of fluid flux modeling and sensor spacing, which is related to streambed scour and deposition. These parameters can also be used to determine the reliability of flux estimates from the use of heat as a tracer.

Groundwater↗

USGS2016

Influence of vertical and lateral heat transfer on permafrost thaw, peatland landscape transition, and groundwater flow

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.

Northwest Territories↗

USGS2014

New permafrost is forming around shrinking Arctic lakes, but will it last?

Widespread lake shrinkage in cold regions has been linked to climate warming and permafrost thaw. Permafrost aggradation, however, has been observed within the margins of recently receded lakes, in seeming contradiction of climate warming. Here permafrost aggradation dynamics are examined at Twelvemile Lake, a retreating lake in interior Alaska. Observations reveal patches of recently formed permafrost within the dried lake margin, colocated with discrete bands of willow shrub. We test ecological succession, which alters shading, infiltration, and heat transport, as the driver of aggradation using numerical simulation of variably saturated groundwater flow and heat transport with phase change (i.e., freeze-thaw). Simulations support permafrost development under current climatic conditions, but only when net effects of vegetation on soil conditions are incorporated, thus pointing to the role of ecological succession. Furthermore, model results indicate that permafrost aggradation is transitory with further climate warming, as new permafrost thaws within seven decades.

Alaska↗