Geology topics

A.M. Provost

Publications and source records attributed to A.M. Provost.

3 recordsLinked to original sources

USGS2012

Glaciation and regional groundwater flow in the Fennoscandian shield

Regional-scale groundwater flow modeling of the Fennoscandian shield suggests that groundwater flow can be strongly affected by future climate change and glaciation. We considered variable-density groundwater flow in a 1500-km-long and approximately 10-km-deep cross-section through southern Sweden. Groundwater flow and shield brine transport in the cross-sectional model were analyzed under projected surface conditions for the next 140 ka. Simulations suggest that blockage of recharge and discharge by low-permeability permafrost or cold-based ice causes sinking of brine and consequent freshening of near-surface water in areas of natural discharge. Although recharge of basal meltwater is limited by the requirement that water pressure at the base of the ice sheet not exceed the pressure exerted by the weight of the ice, warm-based ice with basal melting creates a potential for groundwater recharge rates much larger than those of present, ice-free conditions. In the simulations, regional-scale redistribution of recharged water by subsurface flow is minor over the duration of a glacial advance (approximately 10 ka). During glacial retreat, significant upward flow of groundwater may occur below the ice sheet owing to pressure release. If the mechanical loading efficiency of the rocks is high, both subsurface penetration of meltwater during glacial advance and up-flow during glacial retreat are reduced because of loading-induced pressure changes. The maximum rate of groundwater discharge in the simulations occurs at the receding ice margin, and some discharge occurs below incursive postglacial seas. Recharge of basal meltwater could decrease the concentration of dissolved solids significantly below present-day levels at depths of up to several kilometers and may bring oxygenated conditions to an otherwise reducing chemical environment for periods exceeding 10 ka.

Fennoscandian Shield↗

USGS2009

Recent experimental data may point to a greater role for osmotic pressures in the subsurface

Uncertainty about the origin of anomalous fluid pressures in certain geologic settings has caused researchers to take a second look at osmosis, or flow driven by chemical potential differences, as a pressure‐generating process in the subsurface. Interest in geological osmosis has also increased because of an in situ experiment by Neuzil (2000) suggesting that Pierre Shale could generate large osmotic pressures when highly compacted. In the last few years, additional laboratory and in situ experiments have greatly increased the number of data on osmotic properties of argillaceous formations, but they have not been systematically examined. In this paper we compile these data and explore their implications for osmotic pressure generation in subsurface systems. Rather than base our analysis on osmotic efficiencies, which depend strongly on concentration, we calculated values of a quantity we term osmotic specific surface area ( A so ) that, in principle, is a property of the porous medium only. The A so values are consistent with a surprisingly broad spectrum of osmotic behavior in argillaceous formations, and all the formations tested exhibited at least a modest ability to generate osmotic pressure. It appears possible that under appropriate conditions some formations can be highly effective osmotic membranes able to generate osmotic pressures exceeding 30 MPa (3 km of head) at porosities as high as ∼0.1 and pressures exceeding 10 MPa at porosities as high as ∼0.2. These findings are difficult to reconcile with the lack of compelling field evidence for osmotic pressures, and we propose three explanations for the disparity: (1) Our analysis is flawed and argillaceous formations are less effective osmotic membranes than it suggests; (2) the necessary subsurface conditions, significant salinity differences within intact argillaceous formations, are rare; or (3) osmotic pressures are unlikely to be detected and are not recognized when encountered. The last possibility, that osmotic pressures routinely escape detection or are attributed to other mechanisms, has important implications for understanding subsurface flow regimes.

Water Resources Research↗

USGS2002

SUTRA: A model for 2D or 3D saturated-unsaturated, variable-density ground-water flow with solute or energy transport

SUTRA (Saturated-Unsaturated Transport) is a computer program that simulates fluid movement and the transport of either energy or dissolved substances in a subsurface environment. This upgraded version of SUTRA adds the capability for three-dimensional simulation to the former code (Voss, 1984), which allowed only two-dimensional simulation. The code employs a two- or three-dimensional finite-element and finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated: 1) fluid density-dependent saturated or unsaturated ground-water flow; and 2) either (a) transport of a solute in the ground water, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay; or (b) transport of thermal energy in the ground water and solid matrix of the aquifer. SUTRA may also be used to simulate simpler subsets of the above processes. A flow-direction-dependent dispersion process for anisotropic media is also provided by the code and is introduced in this report. As the primary calculated result, SUTRA provides fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA flow simulation may be employed for two-dimensional (2D) areal, cross sectional and three-dimensional (3D) modeling of saturated ground-water flow systems, and for cross sectional and 3D modeling of unsaturated zone flow. Solute-transport simulation using SUTRA may be employed to model natural or man-induced chemical-species transport including processes of solute sorption, production, and decay. For example, it may be applied to analyze ground-water contaminant transport problems and aquifer restoration designs. In addition, solute-transport simulation with SUTRA may be used for modeling of variable-density leachate movement, and for cross sectional modeling of saltwater intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between freshwater and saltwater. SUTRA energy-transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal-energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. Mesh construction, which is quite flexible for arbitrary geometries, employs quadrilateral finite elements in 2D Cartesian or radial-cylindrical coordinate systems, and hexahedral finite elements in 3D systems. 3D meshes are currently restricted to be logically rectangular; in other words, they are similar to deformable finite-difference-style grids. Permeabilities may be anisotropic and may vary in both direction and magnitude throughout the system, as may most other aquifer and fluid properties. Boundary conditions, sources and sinks may be time dependent. A number of input data checks are made to verify the input data set. An option is available for storing intermediate results and restarting a simulation at the intermediate time. Output options include fluid velocities, fluid mass and solute mass or energy budgets, and time-varying observations at points in the system. Both the mathematical basis for SUTRA and the program structure are highly general, and are modularized to allow for straightforward addition of new methods or processes to the simulation. The FORTRAN-90 coding stresses clarity and modularity rather than efficiency, providing easy access for later modifications.

Water-Resources Investigations Report↗