Search USGS⌕ Search

Geology topics

Clifford I. Voss

Publications and source records attributed to Clifford I. Voss.

At least 19 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↗
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↗
USGS2022

Hydrologic modeling of a perennial firn aquifer in southeast Greenland

A conceptual model, based on field observations and assumed physics of a perennial firn aquifer near Helheim Glacier (southeast Greenland), is evaluated via steady-state 2-D simulation of liquid water flow and energy transport with phase change. The simulation approach allows natural representation of flow and energy advection and conduction that occur in vertical meltwater recharge through the unsaturated zone and in lateral flow within the saturated aquifer. Agreement between measured and simulated aquifer geometry, temperature, and recharge and discharge rates confirms that the conceptual field-data-based description of the aquifer is consistent with the primary physical processes of groundwater flow, energy transport and phase change. Factors that are found to control simulated aquifer configuration include surface temperature, meltwater recharge rate, residual total-water saturation and capillary fringe thickness. Simulation analyses indicate that the size of perennial firn aquifers depends primarily on recharge rates from surface snowmelt. Results also imply that the recent aquifer expansion, likely due to a warming climate, may eventually produce lakes on the ice-sheet surface that would affect the surface energy balance.

Journal of Glaciology↗
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

Simulating water and heat transport with freezing and cryosuction in unsaturated soil: Comparing an empirical, semi-empirical and physically-based approach

Freezing of unsaturated soil is an important process that influences runoff and infiltration in cold-climate regions. We used a simple numerical model to simulate water and heat transport with phase change in unsaturated soil via three different approaches: empirical, semi-empirical and physically based. We compared the performance and parameterization of each approach through testing on three experimental datasets. All approaches reproduced the observed unsaturated freezing process satisfactorily. The empirical cryosuction equation used in this study managed to capture observed cryosuction with a fixed empirical parameter value. The semi-empirical version therefore does not require calibration of a specific frozen soil related parameter. In view of simplicity, small computational demand and accurate performance, all three approaches are suitable for implementation in land-use schemes, catchment scale hydrological models, or multi-dimensional thermo-hydrological models.

Advances in Water Resources↗
USGS2020

The firn meltwater Retention Model Intercomparison Project (RetMIP): Evaluation of nine firn models at four weather station sites on the Greenland ice sheet

Perennial snow, or firn, covers 80 % of the Greenland ice sheet and has the capacity to retain surface meltwater, influencing the ice sheet mass balance and contribution to sea-level rise. Multilayer firn models are traditionally used to simulate firn processes and estimate meltwater retention. We present, intercompare and evaluate outputs from nine firn models at four sites that represent the ice sheet's dry snow, percolation, ice slab and firn aquifer areas. The models are forced by mass and energy fluxes derived from automatic weather stations and compared to firn density, temperature and meltwater percolation depth observations. Models agree relatively well at the dry-snow site while elsewhere their meltwater infiltration schemes lead to marked differences in simulated firn characteristics. Models accounting for deep meltwater percolation overestimate percolation depth and firn temperature at the percolation and ice slab sites but accurately simulate recharge of the firn aquifer. Models using Darcy's law and bucket schemes compare favorably to observed firn temperature and meltwater percolation depth at the percolation site, but only the Darcy models accurately simulate firn temperature and percolation at the ice slab site. Despite good performance at certain locations, no single model currently simulates meltwater infiltration adequately at all sites. The model spread in estimated meltwater retention and runoff increases with increasing meltwater input. The highest runoff was calculated at the KAN_U site in 2012, when average total runoff across models ( ±2 σ ) was 353±610 mm w.e. (water equivalent), about 27±48 % of the surface meltwater input. We identify potential causes for the model spread and the mismatch with observations and provide recommendations for future model development and firn investigation.

The Cryosphere↗
USGS2020

Increasing threat of coastal groundwater hazards from sea-level rise in California

Projected sea-level rise will raise coastal water tables, resulting in groundwater hazards that threaten shallow infrastructure and coastal ecosystem resilience. Here we model a range of sea-level rise scenarios to assess the responses of water tables across the diverse topography and climates of the California coast. With 1 m of sea-level rise, areas flooded from below are predicted to expand ~50–130 m inland, and low-lying coastal communities such as those around San Francisco Bay are most at risk. Coastal topography is a controlling factor; long-term rising water tables will intercept low-elevation drainage features, allowing for groundwater discharge that damps the extent of shoaling in ~70% (68.9–82.2%) of California’s coastal water tables. Ignoring these topography-limited responses increases flooded-area forecasts by ~20% and substantially underestimates saltwater intrusion. All scenarios estimate that areas with shallow coastal water tables will shrink as they are inundated by overland flooding or are topographically limited from rising inland.

California↗
USGS2020

Wildfire-initiated talik development exceeds current thaw projections: Observations and models from Alaska's continuous permafrost zone

As the Arctic warms and wildfire occurrence increases, talik formation in permafrost regions is projected to expand and affect the cycling of water and carbon. Yet, few unified field and modeling studies have examined this process in detail, particularly in areas of continuous permafrost. We address this gap by presenting multimethod, multiseasonal geophysical measurements of permafrost and liquid‐water content that reveal substantial talik development in response to recent wildfire in continuous permafrost of boreal Alaska. Results from observation‐based cryohydrogeologic model simulations suggest that predisturbance subsurface conditions are key factors influencing thaw response to fire disturbance and air temperature warming. Our high‐resolution integrated study illustrates enhanced vulnerability of boreal continuous permafrost, with observed talik formation that exceeds coarse‐scale model projections by ~100 years even under the most extreme future emissions scenario. Results raise important scaling questions for representing extreme permafrost thaw phenomena of growing widespread importance in large‐scale predictive models.

Alaska↗
USGS2020

The influence of snow cover, air temperature, and groundwater flow on the active-layer thermal regime of Arctic hillslopes drained by water tracks

Permafrost in Arctic watersheds limits soil biological activity to a thin, seasonally thawed active layer that contributes water to streams. In many hillslopes, relatively wet drainage features called water tracks have distinct freeze-thaw patterns that affect groundwater flow and storage, and thus the export of heat and solutes to Arctic streams. This study uses groundwater flow and energy transport models to examine potential controls on the timing and duration of freeze–thaw conditions and the magnitude of temperature fluctuations within water tracks and their adjacent hillslopes. The simulated length of the active-layer thaw season varies by 1 month over the range of snow-cover and mean annual air-temperature scenarios simulated. The timing and duration of freezing is particularly sensitive to depth and duration of snow cover. Thus, the deeper snowpack covers that can accumulate in water tracks contribute to their more persistent thaw conditions and their ability to conduct groundwater downslope. A three-dimensional simulation shows that during the summer thaw season, the water track captures groundwater laterally from half way across the hillslope. The models presented here elucidate key mechanisms driving small-scale variation in the active-layer thermal regime of tundra hillslopes, which may be responsible for changes in drainage-network geometry and Arctic biogeochemical fluxes under a warming climate.

Alaska↗
USGS2020

Water tracks enhance water flow above permafrost in upland Arctic Alaska hillslopes

Upland permafrost regions occupy approximately one third of the Arctic landscape. In upland regions, hydrologic fluxes are influenced by water tracks, curvilinear features on hillslopes that preferentially fill with and route water in response to snowmelt and rainfall when the soil above continuous permafrost thaws in the summer. As continued warming of the Arctic may alter hydrologic cycling leading to increased frequency of extreme hydrologic events like drought and flooding as well as modification of biogeochemical cycling, it is imperative to untangle the interplay between precipitation, runoff, and subsurface flow as water is routed from upland Arctic regions to the Arctic Ocean. This study quantifies how ground surface temperatures affect groundwater discharge from hillslopes with water tracks in the upland Arctic by employing a three-dimensional, physically based subsurface flow model with variable saturation and freeze and thaw capabilities that is calibrated to field measurements from the Upper Kuparuk River watershed on the North Slope of Alaska, USA. Model analysis indicates that higher ground surface temperatures along water track hillslopes promote increases in groundwater discharge where water tracks act as conduits for large recharge events and continue to discharge groundwater into the autumn after the adjacent hillslope has frozen. Simulating the conditions that distinguish water tracks from their hillslope watersheds changes subsurface water storage and ground thermal responses but does not alter the total magnitude of groundwater discharge outside of parameter uncertainty. These findings suggest that water tracks play a complex and critical role in hydrologic cycles of the upland Arctic.

Alaska↗
USGS2019

SUTRA, a model for saturated-unsaturated, variable-density groundwater flow with solute or energy transport—Documentation of generalized boundary conditions, a modified implementation of specified pressures and concentrations or temperatures, and the lake capability

Version 3.0 of the SUTRA groundwater modeling program offers three new capabilities: generalized boundary conditions, a modified implementation of specified pressures and concentrations or temperatures, and lakes. Two new types of “generalized” boundary conditions facilitate simulation of a wide range of hydrologic processes that interact with the groundwater model, such as rivers, drains, and evapotranspiration. For generalized-flow boundary conditions, gain (inflow) or loss (outflow) of fluid mass varies linearly with pressure, subject to optional upper and lower limits on flow and (or) pressure. For generalized-transport boundary conditions, gain or loss of solute mass or energy varies linearly with concentration or temperature, respectively. Two of the original types of SUTRA boundary conditions—specified-pressure and specified-concentration or temperature—have been modified such that user-specified, conductance-like factors (known as GNUP and GNUU in previous versions of SUTRA) are no longer required. The new lake capability works with all types of SUTRA boundary conditions, including the new generalized boundary conditions, to enable simulation of the interaction of groundwater flow and transport with lake water “ponded” on the surface of a three-dimensional model. SUTRA uses the topography of the top surface of the model, or, optionally, user-specified lake-bottom elevations, to identify potential lakes automatically. Increases and decreases in lake stage can cause lakes to coalesce and divide, respectively. The lake capability may be used with saturated or unsaturated flow and solute or energy transport.

Techniques and Methods↗
USGS2019

Modeling δ18O as an early indicator of regime shift arising from salinity stress in coastal vegetation

In many important coastal habitats, a combination of increasing soil salinization due to sea level rise, reduced precipitation and storm surges may induce regime shift from salinity-intolerant glycophytic vegetation to salinity-tolerant halophytic species. Early detection of regime shift due to salinity stress in vegetation may facilitate conservation efforts. It has been shown that the 18 O value of water in the xylem of trees can be used as a surrogate for salinity in the rooting zone of plants. Coupling measured δ 18 O values in the tree xylem with simulated δ 18 O values in trees and salinity in the vadose zone can be used to investigate competitive responses of glycophytic versus halophytic trees. MANTRA-O18 simulations suggest that the impacts of salinization on diminishing the resilience of salinity-intolerant trees can be detected up to 25 years before the glycophytic trees are threatened with regime shift to halophytic species. This early detection provides critical lead time and valuable information and insights useful for planning adaptation strategy to mitigate against the adverse impacts of sea level rise and climate change.

Hydrogeology Journal↗
USGS2019

Development of perennial thaw zones in boreal hillslopes enhances potential mobilization of permafrost carbon

Permafrost thaw alters subsurface flow in boreal regions that in turn influences the magnitude, seasonality, and chemical composition of streamflow. Prediction of these changes is challenged by incomplete knowledge of timing, flowpath depth, and amount of groundwater discharge to streams in response to thaw. One important phenomenon that may affect flow and transport through boreal hillslopes is development of lateral perennial thaw zones (PTZs), the existence of which is here supported by geophysical observations and cryohydrogeologic modeling. Model results link thaw to enhanced and seasonally-extended baseflow, which have implications for mobilization of soluble constituents. Results demonstrate the sensitivity of PTZ development to organic layer thickness and near-surface factors that mediate heat exchange at the atmosphere/ground-surface interface. Study findings suggest that PTZs serve as a detectable precursor to accelerated permafrost degradation. This study provides important contextual insight on a fundamental thermo-hydrologic process that can enhance terrestrial-to-aquatic transfer of permafrost carbon, nitrogen, and mercury previously sequestered in thawing watersheds.

Alaska↗
USGS2018

Celebrating 50 years of SWIMs (Salt Water Intrusion Meetings)

The Salt Water Intrusion Meetings, or SWIMs, are a series of meetings that focus on seawater intrusion in coastal aquifers and other salinisation processes. 2018 marks the 50th year of the SWIM and the 25th biennial meeting. The SWIM proceedings record half a century of research progress on site characterisation, geophysical and geochemical techniques, variable-density flow, modelling, and water management. The SWIM is positioning itself to remain a viable platform for discussing the coastal aquifer management challenges of the next 50 years.

Hydrogeology Journal↗
USGS2018

The role of frozen soil in groundwater discharge predictions for warming alpine watersheds

Climate warming may alter the quantity and timing of groundwater discharge to streams in high alpine watersheds due to changes in the timing of the duration of seasonal freezing in the subsurface and snowmelt recharge. It is imperative to understand the effects of seasonal freezing and recharge on groundwater discharge to streams in warming alpine watersheds as streamflow originating from these watersheds is a critical water resource for downstream users. This study evaluates how climate warming may alter groundwater discharge due to changes in seasonally frozen ground and snowmelt using a 2‐D coupled flow and heat transport model with freeze and thaw capabilities for variably saturated media. The model is applied to a representative snowmelt‐dominated watershed in the Rocky Mountains of central Colorado, USA, with snowmelt time series reconstructed from a 12 year data set of hydrometeorological records and satellite‐derived snow covered area. Model analyses indicate that the duration of seasonal freezing in the subsurface controls groundwater discharge to streams, while snowmelt timing controls groundwater discharge to hillslope faces. Climate warming causes changes to subsurface ice content and duration, rerouting groundwater flow paths but not altering the total magnitude of future groundwater discharge outside of the bounds of hydrologic parameter uncertainties. These findings suggest that frozen soil routines play an important role for predicting the future location of groundwater discharge in watersheds underlain by seasonally frozen ground.

Colorado↗
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↗
USGS2018

Most atolls will be uninhabitable by the mid-21st century because of sea-level rise exacerbating wave-driven flooding

Sea levels are rising, with the highest rates in the tropics, where thousands of low-lying coral atoll islands are located. Most studies on the resilience of these islands to sea-level rise have projected that they will experience minimal inundation impacts until at least the end of the 21st century. However, these have not taken into account the additional hazard of wave-driven overwash or its impact on freshwater availability. We project the impact of sea-level rise and wave-driven flooding on atoll infrastructure and freshwater availability under a variety of climate change scenarios. We show that, on the basis of current greenhouse gas emission rates, the nonlinear interactions between sea-level rise and wave dynamics over reefs will lead to the annual wave-driven overwash of most atoll islands by the mid-21st century. This annual flooding will result in the islands becoming uninhabitable because of frequent damage to infrastructure and the inability of their freshwater aquifers to recover between overwash events. This study provides critical information for understanding the timing and magnitude of climate change impacts on atoll islands that will result in significant, unavoidable geopolitical issues if it becomes necessary to abandon and relocate low-lying island states.

Science Advances↗