Geology topics

Clifford I. Voss

Publications and source records attributed to Clifford I. Voss.

At least 37 records · Page 2Linked to original sources

USGS2017

The role of uplift and erosion in the persistence of saline groundwater in the shallow subsurface

In many regions of the world, the shallow (<300 m) subsurface is replenished with meteoric recharge within a few centuries or millennia, but in some regions saline groundwater persists despite abundant rainfall. Analyses of the flushing rate of shallow groundwater usually consider the permeability and recharge rate and a static landscape. The influence of landscape evolution can become important over millions of years, however. Here we present numerical simulations of fluid flow and transport in the top 1 km of a sedimentary foreland basin dominated by aquitards, where the rate of uplift and erosion (20 m Ma −1 ) balances that of meteoric flushing. Paleozoic age saline groundwater and brine persist at shallow depths that might otherwise have contained potable water. Similar hydrogeologic conditions, and uplift and erosion rates, likely exist in many other regions of the world, where a moving landscape has probably never been considered as an important contributor to groundwater quality.

Geophysical Research Letters↗

USGS2017

Seawater-flooding events and impact on freshwater lenses of low-lying islands: Controlling factors, basic management and mitigation

An unprecedented set of hydrologic observations was collected after the Dec 2008 seawater-flooding event on Roi-Namur, Kwajalein Atoll, Republic of the Marshall Islands. By two days after the seawater flooding that occurred at the beginning of dry season, the observed salinity of water withdrawn by the island’s main skimming well increased to 100% seawater concentration, but by ten days later already decreased to only 10–20% of seawater fraction. However, the damaging impact on the potability of the groundwater supply (when pumped water had concentrations above 1% seawater fraction) lasted 22 months longer. The data collected make possible analyses of the hydrologic factors that control recovery and management of the groundwater-supply quality on Roi-Namur and on similar low-lying islands. With the observed data as a guide, three-dimensional numerical-model simulation analyses reveal how recovery is controlled by the island’s hydrology. These also allow evaluation of the efficacy of basic water-quality management/mitigation alternatives and elucidate how groundwater withdrawal and timing of the seawater-flooding event affect the length of recovery. Simulations show that, as might be expected, by adding surplus captured rainwater as artificial recharge, the freshwater-lens recovery period (after which potable groundwater may again be produced) can be shortened, with groundwater salinity remaining lower even during the dry season, a period during which no artificial recharge is applied. Simulations also show that the recovery period is not lengthened appreciably by groundwater withdrawals during recovery. Simulations further show that had the flooding event occurred at the start of the wet season, the recovery period would have been about 25% (5.5 months) shorter than actually occurred during the monitored flood that occurred at the dry-season start. Finally, analyses show that artificial recharge improves freshwater-lens water quality, making possible longer use of groundwater as a water supply throughout each year, even when no seawater flooding has occurred.

Journal of Hydrology↗

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↗

USGS2015

Effect of permafrost thaw on the dynamics of lakes recharged by ice-jam floods: case study in Yukon Flats, Alaska

Large river floods are a key water source for many lakes in fluvial periglacial settings. Where permeable sediments occur, the distribution of permafrost may play an important role in the routing of floodwaters across a floodplain. This relationship is explored for lakes in the discontinuous permafrost of Yukon Flats, interior Alaska, using an analysis that integrates satellite-derived gradients in water surface elevation, knowledge of hydrogeology, and hydrologic modeling. We observed gradients in water surface elevation between neighboring lakes ranging from 0.001 to 0.004. These high gradients, despite a ubiquitous layer of continuous shallow gravel across the flats, are consistent with limited groundwater flow across lake basins resulting from the presence of permafrost. Permafrost impedes the propagation of floodwaters in the shallow subsurface and constrains transmission to “fill-and-spill” over topographic depressions (surface sills), as we observed for the Twelvemile-Buddy Lake pair following a May 2013 ice-jam flood on the Yukon River. Model results indicate that permafrost table deepening of 1–11 m in gravel, depending on watershed geometry and subsurface properties, could shift important routing of floodwater to lakes from overland flow (fill-and-spill) to shallow groundwater flow (“fill-and-seep”). Such a shift is possible in the next several hundred years of ground surface warming, and may bring about more synchronous water level changes between neighboring lakes following large flood events. This relationship offers a potentially useful tool, well-suited to remote sensing, for identifying long-term changes in shallow groundwater flow resulting from thawing of permafrost.

Alaska↗

USGS2014

Understanding heat and groundwater flow through continental flood basalt provinces: insights gained from alternative models of permeability/depth relationships for the Columbia Plateau, USA

Heat-flow mapping of the western USA has identified an apparent low-heat-flow anomaly coincident with the Columbia Plateau Regional Aquifer System, a thick sequence of basalt aquifers within the Columbia River Basalt Group (CRBG). A heat and mass transport model (SUTRA) was used to evaluate the potential impact of groundwater flow on heat flow along two different regional groundwater flow paths. Limited in situ permeability (k) data from the CRBG are compatible with a steep permeability decrease (approximately 3.5 orders of magnitude) at 600–900 m depth and approximately 40°C. Numerical simulations incorporating this permeability decrease demonstrate that regional groundwater flow can explain lower-than-expected heat flow in these highly anisotropic (k x /k z ~ 10 4 ) continental flood basalts. Simulation results indicate that the abrupt reduction in permeability at approximately 600 m depth results in an equivalently abrupt transition from a shallow region where heat flow is affected by groundwater flow to a deeper region of conduction-dominated heat flow. Most existing heat-flow measurements within the CRBG are from shallower than 600 m depth or near regional groundwater discharge zones, so that heat-flow maps generated using these data are likely influenced by groundwater flow. Substantial k decreases at similar temperatures have also been observed in the volcanic rocks of the adjacent Cascade Range volcanic arc and at Kilauea Volcano, Hawaii, where they result from low-temperature hydrothermal alteration.

Idaho;Oregon;Washington↗

USGS2014

Fluid pressure responses for a Devil's Slide-like system: problem formulation and simulation

This study employs a hydrogeologic simulation approach to investigate subsurface fluid pressures for a landslide-prone section of the central California, USA, coast known as Devil's Slide. Understanding the relative changes in subsurface fluid pressures is important for systems, such as Devil's Slide, where slope creep can be interrupted by episodic slip events. Surface mapping, exploratory core, tunnel excavation records, and dip meter data were leveraged to conceptualize the parameter space for three-dimensional (3D) Devil's Slide-like simulations. Field observations (i.e. seepage meter, water retention, and infiltration experiments; well records; and piezometric data) and groundwater flow simulation (i.e. one-dimensional vertical, transient, and variably saturated) were used to design the boundary conditions for 3D Devil's Slide-like problems. Twenty-four simulations of steady-state saturated subsurface flow were conducted in a concept-development mode. Recharge, heterogeneity, and anisotropy are shown to increase fluid pressures for failure-prone locations by up to 18.1, 4.5, and 1.8% respectively. Previous estimates of slope stability, driven by simple water balances, are significantly improved upon with the fluid pressures reported here. The results, for a Devil's Slide-like system, provide a foundation for future investigations

California↗

USGS2014

Climate change impacts on the temperature and magnitude of groundwater discharge from shallow, unconfined aquifers

Cold groundwater discharge to streams and rivers can provide critical thermal refuge for threatened salmonids and other aquatic species during warm summer periods. Climate change may influence groundwater temperature and flow rates, which may in turn impact riverine ecosystems. This study evaluates the potential impact of climate change on the timing, magnitude, and temperature of groundwater discharge from small, unconfined aquifers that undergo seasonal freezing and thawing. Seven downscaled climate scenarios for 2046–2065 were utilized to drive surficial water and energy balance models (HELP3 and ForHyM2) to obtain future projections for daily ground surface temperature and groundwater recharge. These future surface conditions were then applied as boundary conditions to drive subsurface simulations of variably saturated groundwater flow and energy transport. The subsurface simulations were performed with the U.S. Geological Survey finite element model SUTRA that was recently modified to include the dynamic freeze-thaw process. The SUTRA simulations indicate a potential rise in the magnitude (up to 34%) and temperature (up to 3.6°C) of groundwater discharge to the adjacent river during the summer months due to projected increases in air temperature and precipitation. The thermal response of groundwater to climate change is shown to be strongly dependent on the aquifer dimensions. Thus, the simulations demonstrate that the thermal sensitivity of aquifers and baseflow-dominated streams to decadal climate change may be more complex than previously thought. Furthermore, the results indicate that the probability of exceeding critical temperature thresholds within groundwater-sourced thermal refugia may significantly increase under the most extreme climate scenarios.

Water Resources Research↗

USGS2014

Water resources management in the Ganges Basin: a comparison of three strategies for conjunctive use of groundwater and surface water

The most difficult water resources management challenge in the Ganges Basin is the imbalance between water demand and seasonal availability. More than 80 % of the annual flow in the Ganges River occurs during the 4-month monsoon, resulting in widespread flooding. During the rest of the year, irrigation, navigation, and ecosystems suffer because of water scarcity. Storage of monsoonal flow for utilization during the dry season is one approach to mitigating these problems. Three conjunctive use management strategies involving subsurface water storage are evaluated in this study: Ganges Water Machine (GWM), Pumping Along Canals (PAC), and Distributed Pumping and Recharge (DPR). Numerical models are used to determine the efficacy of these strategies. Results for the Indian State of Uttar Pradesh (UP) indicate that these strategies create seasonal subsurface storage from 6 to 37 % of the yearly average monsoonal flow in the Ganges exiting UP over the considered range of conditions. This has clear implications for flood reduction, and each strategy has the potential to provide irrigation water and to reduce soil waterlogging. However, GWM and PAC require significant public investment in infrastructure and management, as well as major shifts in existing water use practices; these also involve spatially-concentrated pumping, which may induce land subsidence. DPR also requires investment and management, but the distributed pumping is less costly and can be more easily implemented via adaptation of existing water use practices in the basin.

Ganges Basin↗

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↗

USGS2014

Analytical solutions for benchmarking cold regions subsurface water flow and energy transport models: one-dimensional soil thaw with conduction and advection

Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.

Advances in Water Resources↗

USGS2013

The transboundary non-renewable Nubian Aquifer System of Chad, Egypt, Libya and Sudan: Classical groundwater questions and parsimonious hydrogeologic analysis and modeling

Parsimonious groundwater modeling provides insight into hydrogeologic functioning of the Nubian Aquifer System (NAS), the world’s largest non-renewable groundwater system (belonging to Chad, Egypt, Libya, and Sudan). Classical groundwater-resource issues exist (magnitude and lateral extent of drawdown near pumping centers) with joint international management questions regarding transboundary drawdown. Much of NAS is thick, containing a large volume of high-quality groundwater, but receives insignificant recharge, so water-resource availability is time-limited. Informative aquifer data are lacking regarding large-scale response, providing only local-scale information near pumps. Proxy data provide primary underpinning for understanding regional response: Holocene water-table decline from the previous pluvial period, after thousands of years, results in current oasis/sabkha locations where the water table still intersects the ground. Depletion is found to be controlled by two regional parameters, hydraulic diffusivity and vertical anisotropy of permeability. Secondary data that provide insight are drawdowns near pumps and isotope-groundwater ages (million-year-old groundwaters in Egypt). The resultant strong simply structured three-dimensional model representation captures the essence of NAS regional groundwater-flow behavior. Model forecasts inform resource management that transboundary drawdown will likely be minimal—a nonissue—whereas drawdown within pumping centers may become excessive, requiring alternative extraction schemes; correspondingly, significant water-table drawdown may occur in pumping centers co-located with oases, causing oasis loss and environmental impacts.

Hydrogeology Journal↗

USGS2013

Impacts on groundwater recharge areas of megacity pumping: analysis of potential contamination of Kolkata, India, water supply

Water supply to the world's megacities is a problem of quantity and quality that will be a priority in the coming decades. Heavy pumping of groundwater beneath these urban centres, particularly in regions with low natural topographic gradients, such as deltas and floodplains, can fundamentally alter the hydrological system. These changes affect recharge area locations, which may shift closer to the city centre than before development, thereby increasing the potential for contamination. Hydrogeological simulation analysis allows evaluation of the impact on past, present and future pumping for the region of Kolkata, India, on recharge area locations in an aquifer that supplies water to over 13 million people. Relocated recharge areas are compared with known surface contamination sources, with a focus on sustainable management of this urban groundwater resource. The study highlights the impacts of pumping on water sources for long-term development of stressed city aquifers and for future water supply in deltaic and floodplain regions of the world.

West Bengal↗

USGS2013

Linkages between lake shrinkage/expansion and sublacustrine permafrost distribution determined from remote sensing of interior Alaska, USA

[1] Linkages between permafrost distribution and lake surface-area changes in cold regions have not been previously examined over a large scale because of the paucity of subsurface permafrost information. Here, a first large-scale examination of these linkages is made over a 5150 km 2 area of Yukon Flats, Alaska, USA, by evaluating the relationship between lake surface-area changes during 1979–2009, derived from Landsat satellite data, and sublacustrine groundwater flow-path connectivity inferred from a pioneering, airborne geophysical survey of permafrost. The results suggest that the shallow (few tens of meters) thaw state of permafrost has more influence than deeper permafrost conditions on the evolving water budgets of lakes on a multidecadal time scale. In the region studied, these key shallow aquifers have high hydraulic conductivity and great spatial variability in thaw state, making groundwater flow and associated lake level evolution particularly sensitive to climate change owing to the close proximity of these aquifers to the atmosphere.

Alaska↗

USGS2013

Impacts of climate, lake size, and supra- and sub-permafrost groundwater flow on lake-talik evolution, Yukon Flats, Alaska (USA)

In cold regions, hydrologic systems possess seasonal and perennial ice-free zones (taliks) within areas of permafrost that control and are enhanced by groundwater flow. Simulation of talik development that follows lake formation in watersheds modeled after those in the Yukon Flats of interior Alaska (USA) provides insight on the coupled interaction between groundwater flow and ice distribution. The SUTRA groundwater simulator with freeze–thaw physics is used to examine the effect of climate, lake size, and lake–groundwater relations on talik formation. Considering a range of these factors, simulated times for a through-going sub-lake talik to form through 90 m of permafrost range from ∼200 to > 1,000 years (vertical thaw rates < 0.1–0.5 m yr −1 ). Seasonal temperature cycles along lake margins impact supra-permafrost flow and late-stage cryologic processes. Warmer climate accelerates complete permafrost thaw and enhances seasonal flow within the supra-permafrost layer. Prior to open talik formation, sub-lake permafrost thaw is dominated by heat conduction. When hydraulic conditions induce upward or downward flow between the lake and sub-permafrost aquifer, thaw rates are greatly increased. The complexity of ground-ice and water-flow interplay, together with anticipated warming in the arctic, underscores the utility of coupled groundwater-energy transport models in evaluating hydrologic systems impacted by permafrost.

Alaska↗

USGS2013

Permafrost thaw in a nested groundwater-flow system

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

Hydrogeology Journal↗

USGS2013

Brine intrusion by upconing for a high-level nuclear waste repository at Forsmark: Scoping calculations

SSM currently reviews a license application for a spent nuclear fuel repository that is proposed to be located at Forsmark, Sweden. The repository is to be situated at 500 m depth in the rock and copper canisters are deposited in holes excavated from the tunnel system. To protect the canisters they are surrounded by a bentonite clay buffer, which is to swell when getting in contact with water. The swelling properties are dependent on the salt content of the water and excessively high salt contents may inhibit the swelling. Thus it is important to ensure that the bentonite is not subjected to water with too high salt contents. The salt content of the groundwater increases with depth and is expected to reach levels that may affect buffer performance at large depths. When excavating the repository very high hydraulic gradients are established and water and salt movement from the depth to the repository, so-called ‘upconing’, could possibly occur. The objective of this study is to evaluate the possibility of salt-water migration to the repository. This objective is motivated by the adverse impacts of water with too high salinity entering the repository and by the uncertainty of the relevant hydraulic and hydrogeochemical conditions at the Forsmark site at great depths. To analyse density dependent flow and salt transport at the Forsmark site the USGS’ SUTRA code is used. This study proceeds by finding critical model cases for which upconing does or does not occur, while assessing whether the parameterizations of these cases are realistic for the Forsmark site. In addition, the fall of the upconed salt mound (i.e. downconing) following closure of the repository is also evaluated. In particular the objectives are (1) to determine the factors that control saltwater upconing in a hydrogeological setting representative of Forsmark; (2) to relate these factors to the plausible conditions prevailing at the repository site; (3) to investigate whether the proposed repository is likely to generate saltwater upconing, given the range of uncertainty in hydrogeologic structure and parameter values; and (4) to evaluate the timing of upconing (salinization) and the timing of downconing (freshening) following repository closure for cases where upconing occurs. The results of this simulation analysis show that upconing behavior is strongly affected by the ratio of permeability to porosity in any zone in which upconing might occur. Within the full range of parameters that are likely to occur at the Forsmark site, the model yields either no significant upconing at all during the operational period of the repository or intrusion of brine-type waters after only one to a few decades.

Forsmark↗

USGS2012

Influence of permafrost distribution on groundwater flow in the context of climate-driven permafrost thaw: Example from Yukon Flats Basin, Alaska, United States

Understanding the role of permafrost in controlling groundwater flow paths and fluxes is central in studies aimed at assessing potential climate change impacts on vegetation, species habitat, biogeochemical cycling, and biodiversity. Recent field studies in interior Alaska show evidence of hydrologic changes hypothesized to result from permafrost degradation. This study assesses the hydrologic control exerted by permafrost, elucidates modes of regional groundwater flow for various spatial permafrost patterns, and evaluates potential hydrologic consequences of permafrost degradation. The Yukon Flats Basin (YFB), a large (118,340 km 2 ) subbasin within the Yukon River Basin, provides the basis for this investigation. Model simulations that represent an assumed permafrost thaw sequence reveal the following trends with decreasing permafrost coverage: (1) increased groundwater discharge to rivers, consistent with historical trends in base flow observations in the Yukon River Basin, (2) potential for increased overall groundwater flux, (3) increased spatial extent of groundwater discharge in lowlands, and (4) decreased proportion of suprapermafrost (shallow) groundwater contribution to total base flow. These trends directly affect the chemical composition and residence time of riverine exports, the state of groundwater-influenced lakes and wetlands, seasonal river-ice thickness, and stream temperatures. Presently, the YFB is coarsely mapped as spanning the continuous-discontinuous permafrost transition that model analysis shows to be a critical threshold; thus, the YFB may be on the verge of major hydrologic change should the current permafrost extent decrease. This possibility underscores the need for improved characterization of permafrost and other hydrogeologic information in the region via geophysical techniques, remote sensing, and ground-based observations.

Alaska↗

USGS2012

Modeling radium distribution in coastal aquifers during sea level changes: The Dead Sea case

We present a new approach to studying the behavior of radium isotopes in a coastal aquifer. In order to simulate radium isotope distributions in the dynamic flow field of the Dead Sea aquifer, a multi-species density dependent flow model (SUTRA-MS) was used. Field data show that the activity of 226 Ra decreases from 140 to 60 dpm/L upon entering the aquifer from the Dead Sea, and then further decreases linearly due to mixing with Ra-poor fresh water. On the other hand, an increase is observed in the activity of the shorter-lived isotopes (up to 52 dpm/L 224 Ra and 31 dpm/L 223 Ra), which are relatively low in Dead Sea water (up to 2.5 dpm/L 224 Ra and 0.5 dpm/L 223 Ra). The activities of the short lived radium isotopes also decrease with decreasing salinity, which is due to the effect of salinity on the adsorption of radium. The relationship between 224 Ra and salinity suggests that the adsorption partition coefficient ( K ) is linearly related to salinity. Simulations of the steady-state conditions, show that the distance where equilibrium activity is attained for each radium isotope is affected by the isotope half-life, K and the groundwater velocity, resulting in a longer distance for the long-lived radium isotopes. K affects the radium distribution in transient conditions, especially that of the long-lived radium isotopes. The transient conditions in the Dead Sea system, with a 1 m/yr lake level drop, together with the radium field data, constrains K to be relatively low (<10). Thus, the sharp decrease in 226 Ra cannot be explained by adsorption, and it is better explained by removal via coprecipitation, probably with barite or celestine.

Dead Sea↗