Search USGS⌕ Search

Geology topics

C. E. Neuzil

Publications and source records attributed to C. E. Neuzil.

14 recordsLinked to original sources
USGS2013

Can shale safely host US nuclear waste?

"Even as cleanup efforts after Japan’s Fukushima disaster offer a stark reminder of the spent nuclear fuel (SNF) stored at nuclear plants worldwide, the decision in 2009 to scrap Yucca Mountain as a permanent disposal site has dimmed hope for a repository for SNF and other high-level nuclear waste (HLW) in the United States anytime soon. About 70,000 metric tons of SNF are now in pool or dry cask storage at 75 sites across the United States [Government Accountability Office, 2012], and uncertainty about its fate is hobbling future development of nuclear power, increasing costs for utilities, and creating a liability for American taxpayers [Blue Ribbon Commission on America’s Nuclear Future, 2012].However, abandoning Yucca Mountain could also result in broadening geologic options for hosting America’s nuclear waste. Shales and other argillaceous formations (mudrocks, clays, and similar clay-rich media) have been absent from the U.S. repository program. In contrast, France, Switzerland, and Belgium are now planning repositories in argillaceous formations after extensive research in underground laboratories on the safety and feasibility of such an approach [Blue Ribbon Commission on America’s Nuclear Future, 2012; Nationale Genossenschaft für die Lagerung radioaktiver Abfälle (NAGRA), 2010; Organisme national des déchets radioactifs et des matières fissiles enrichies, 2011]. Other nations, notably Japan, Canada, and the United Kingdom, are studying argillaceous formations or may consider them in their siting programs [Japan Atomic Energy Agency, 2012; Nuclear Waste Management Organization (NWMO), (2011a); Powell et al., 2010]."

Eos, Transactions, American Geophysical Union↗
USGS2012

Geologic isolation of nuclear waste at high latitudes: the role of ice sheets

Geologic isolation of high-level nuclear waste from the biosphere requires special consideration in countries at high latitudes (>40°N) owing to the possibility of future episodes of continental glaciation (Talbot 1999). It is now widely recognized that Pleistocene continental glaciations have had a profound effect on rates of sediment erosion (Cuffey & Paterson 2010) and deformation including tectonic thrusting (Pedersen 2005) as well as groundwater flow (Person et al. 2007; Lemieux et al. 2008a,b,c). In addition, glacial mechanical loads may have generated anomalous, or fossil, pore pressures within certain clay-rich confining units (e.g. Vinard et al. 2001). Because high-level nuclear wastes must be isolated from the biosphere as long as 1 million years (McMurry et al. 2003), the likelihood of one or more continental ice sheets overrunning high-latitude sites must be considered.

Geofluids↗
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↗
USGS2011

Hydromechanical effects of continental glaciation on groundwater systems

Hydromechanical effects of continental ice sheets may involve considerably more than the widely recognized direct compression of overridden terrains by ice load. Lithospheric flexure, which lags ice advance and retreat, appears capable of causing comparable or greater stress changes. Together, direct and flexural loading may increase fluid pressures by tens of MPa in geologic units unable to drain. If so, fluid pressures in low-permeability formations subject to glaciation may have increased and decreased repeatedly during cycles of Pleistocene glaciation and can again in the future. Being asynchronous and normally oriented, direct and flexural loading presumably cause normal and shear stresses to evolve in a complex fashion through much or all of a glacial cycle. Simulations of fractured rock predict permeability might vary by two to three orders of magnitude under similar stress changes as fractures at different orientations are subjected to changing normal and shear stresses and some become critically stressed. Uncertainties surrounding these processes and their interactions, and the confounding influences of surface hydrologic changes, make it challenging to delineate their effects on groundwater flow and pressure regimes with any specificity. To date, evidence for hydromechanical changes caused by the last glaciation is sparse and inconclusive, comprising a few pressure anomalies attributed to the removal of direct ice load. This may change as more data are gathered, and understanding of relevant processes is refined.

Geofluids↗
USGS2010

How permeable are clays and shales?

The permeability of argillaceous formations, although rarely measured and poorly understood, is commonly a critical parameter in analyses of subsurface flow. Data now available suggest a regular relation between permeability and porosity in clays and shales and permeabilities that, even at large scales, are significantly lower than usually assumed. Permeabilities between 10 −23 and 10 −17 m 2 have been obtained at porosities between 0.1 and 0.4 in both laboratory and regional studies. Although it is clear that transmissive fractures or other heterogeneities control the large-scale hydraulic behavior of certain argillaceous units, the permeability of many others is apparently scale independent. These results have significant implications for understanding fluid transport rates and abnormal pressure generation in basins, and could prove important for waste isolation efforts.

Water Resources Research↗
USGS2010

Flow through fractures

Flow through fractures is often idealized as flow between two parallel plates (plane Poiseuille flow). The opening or aperture between parallel plates is unambiguous and its relation to flowrate is well known. However, fractures in rock have uneven walls and a variable aperture. A model for flow in a fracture is proposed wherein the fracture is represented by a set of parallel plate openings with different apertures. The model leads to a modified Poiseuille equation for flow which includes an aperture frequency distribution for the fracture. Any arbitrary aperture distribution can be used; in order to simplify computation and demonstrate the properties of the model a log normal form of distribution is assumed. Even when an analytical form of the distribution is assumed, two parameters, rather than a single value representing ‘aperture size’ are required to determine flowrate. Models of aperture change for a fracture undergoing compression (fracture walls deforming) and extension (fracture walls separating) are developed which constrain the additional parameter and allow calculation of flowrate as a function of mean aperture. The theoretical relationships developed between mean aperture and flowrate can be used to interpret published laboratory data for single fractures.

Water Resources Research↗
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↗
USGS2006

Numerical modeling of a long-term in situ chemical osmosis experiment in the Pierre Shale, South Dakota

We have numerically modeled evolving fluid pressures and concentrations from a nine-year in situ osmosis experiment in the Pierre Shale, South Dakota. These data were obtained and recently interpreted by one of us (C.E.N.) as indicating a potentially significant role for chemical osmosis in media like the Pierre Shale. That analysis considered only the final pressure differentials among boreholes that were assumed to represent osmotic equilibrium. For this study, the system evolution was modeled using a recently developed transient model for membrane transport. The model simulates hydraulically and chemically driven fluid and solute transport. The results yield an estimate of the thickness of the water film between the clay platelets b of 40 A??, which corresponds to an osmotic efficiency ?? of 0.21 for the ambient pore water salinity of 3.5 g/l TDS. These values largely confirm the results of the earlier equilibrium analysis. However, the new model analysis provides additional constraints suggesting that intrinsic permeability k = 1.4 ?? 10-19 m2, specific storage Ss = 1.7 ?? 10-5 m-1, and diffusion coefficient D* = 6 ?? 10-11 m2/s. The k value is larger than certain independent estimates which range from 10-21 to 10-20; it may indicate opening of microcracks during the experiments. The fact that the complex transient pressure and concentration behavior for the individual wells could be reproduced quite accurately, and the inferred parameter values appear to be realistic for the Pierre Shale, suggests that the new model is a useful tool for modeling transient coupled flows in groundwater systems. ?? 2005 Elsevier Ltd. All rights reserved.

Advances in Water Resources↗
USGS2003

Hydromechanical coupling in geologic processes

Earth's porous crust and the fluids within it are intimately linked through their mechanical effects on each other. This paper presents an overview of such "hydromechanical" coupling and examines current understanding of its role in geologic processes. An outline of the theory of hydromechanics and rheological models for geologic deformation is included to place various analytical approaches in proper context and to provide an introduction to this broad topic for nonspecialists. Effects of hydromechanical coupling are ubiquitous in geology, and can be local and short-lived or regional and very long-lived. Phenomena such as deposition and erosion, tectonism, seismicity, earth tides, and barometric loading produce strains that tend to alter fluid pressure. Resulting pressure perturbations can be dramatic, and many so-called "anomalous" pressures appear to have been created in this manner. The effects of fluid pressure on crustal mechanics are also profound. Geologic media deform and fail largely in response to effective stress, or total stress minus fluid pressure. As a result, fluid pressures control compaction, decompaction, and other types of deformation, as well as jointing, shear failure, and shear slippage, including events that generate earthquakes. By controlling deformation and failure, fluid pressures also regulate states of stress in the upper crust. Advances in the last 80 years, including theories of consolidation, transient groundwater flow, and poroelasticity, have been synthesized into a reasonably complete conceptual framework for understanding and describing hydromechanical coupling. Full coupling in two or three dimensions is described using force balance equations for deformation coupled with a mass conservation equation for fluid flow. Fully coupled analyses allow hypothesis testing and conceptual model development. However, rigorous application of full coupling is often difficult because (1) the rheological behavior of geologic media is complex and poorly understood and (2) the architecture, mechanical properties and boundary conditions, and deformation history of most geologic systems are not well known. Much of what is known about hydromechanical processes in geologic systems is derived from simpler analyses that ignore certain aspects of solid-fluid coupling. The simplifications introduce error, but more complete analyses usually are not warranted. Hydromechanical analyses should thus be interpreted judiciously, with an appreciation for their limitations. Innovative approaches to hydromechanical modeling and obtaining critical data may circumvent some current limitations and provide answers to remaining questions about crustal processes and fluid behavior in the crust.

Hydrogeology Journal↗
USGS1998

Glaciation and regional ground-water flow in the Fennoscandian Shield: Site 94

Results from a regional-scale ground-water flow model of the Fennoscandian shield suggest that ground-water flow is strongly affected by surface conditions associated with climatic change and glaciation. The model was used to run a series of numerical simulations of variable-density ground-water flow in a 1500-km-long and approximately 10-km-deep cross-section that passes through southern Sweden. Ground-water flow and shield brine transport in the cross-sectional model are controlled by an assumed time evolution of surface conditions over the next 140 ka. Simulations show that, under periglacial conditions, permafrost may locally or extensively impede the free recharge or discharge of ground water. Below cold-based glacial ice, no recharge or discharge of ground water occurs. Both of these conditions result in the settling of shield brine and consequent freshening of near-surface water in areas of natural discharge blocked by permafrost. The presence of warm-based ice with basal melting creates a potential for ground-water recharge rates much larger than under present, ice-free conditions. Recharging basal meltwater can reach depths of a few kilometers in a few thousand years. The vast majority of recharged water is accommodated through storage in the volume of bedrock below the local area of recharge; regional (lateral) redistribution of recharged water by subsurface flow is minor over the duration of a glacial advance (~10 ka). During glacial retreat, the weight of the ice overlying a given surface location decreases, and significant upward flow of ground water may occur below the ice sheet due to pressure release, despite the continued potential for recharge of basal meltwater. Excess meltwater must exit from below the glacier through subglacial cavities and channels. Subsurface penetration of meltwater during glacial advance and up-flow during glacial retreat are greatest if the loading efficiency of the shield rock is low. The maximum rate of ground-water discharge occurs at the receding ice margin, and some discharge occurs below incursive post-glacial seas. The simulation results suggest that vertical movement of deep shield brines induced by the next few glacial cycles should not increase the concentration of dissolved solids significantly above present-day levels. However, the concentration of dissolved solids should decrease significantly at depths of up to several kilometers during periods of glacial meltwater recharge. The meltwater may reside in the subsurface for periods exceeding 10 ka and may bring oxygenated conditions to an otherwise reducing chemical environment.

SKI Report↗
USGS1983

Regional flow in the Dakota aquifer: A study of the role of confining layers

The Dakota Sandstone in South Dakota is one of the classic artesian aquifers; it was first studied by N. H. Darton at the turn of the century. Since then, hydrogeologists have debated the source of the large quantities of water which have been discharged by artesian flow from the Dakota. Among suggestions for the source of this water are (1) recharge of the aquifer at outcrops in the Black Hills, (2) compaction and compressive storage within the aquifer, (3)leakage through confining layers, and (4) upward flow from the underlying Madison Group limestones. A series of increasingly refined models of the aquifer system in South Dakota have been developed and used for numerical simulations of the ground-water flow. The simulations have provided estimates of leakage through the confining layers. The results indicate that, before development, most of the flow into and out of the Dakota Sandstone occurred as leakage through confining layers and, since development, most of the water released from storage has come from the confining layers. In situ and laboratory hydraulic conductivity measurements have been made for the Cretaceous shale confining layer which overlies the Dakota. These data indicate hydraulic conductivities which are one to three orders of magnitude lower than the conductivities indicated by the numerical analyses; this suggests that the leakage through the confining layer is largely through fractures. The fractures apparently did not influence the laboratory and in situ measurements. To test the conception of flow in the aquifer-confining layer system derived from our analyses, the transport of sulfate in the system was simulated. Simulations using a numerical ground-water transport model were reasonably successful in explaining the present distribution of sulfate in the system. This result increases confidence in the flow system implied by the flow simulations in which leakage through confining layers is dominant.

Nebraska, South Dakota↗
USGS1981

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-I. Theory

Transient pulse testing has been employed increasingly in the laboratory to measure the hydraulic properties of rock samples with low permeability. Several investigators have proposed a mathematical model in terms of an initial-boundary value problem to describe fluid flow in a transient pulse test. However, the solution of this problem has not been available. In analyzing data from the transient pulse test, previous investigators have either employed analytical solutions that are derived with the use of additional, restrictive assumptions, or have resorted to numerical methods. In Part I of this paper, a general, analytical solution for the transient pulse test is presented. This solution is graphically illustrated by plots of dimensionless variables for several cases of interest. The solution is shown to contain, as limiting cases, the more restrictive analytical solutions that the previous investigators have derived. A method of computing both the permeability and specific storage of the test sample from experimental data will be presented in Part II.

International Journal of Rock Mechanics and Mining↗
USGS1981

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-II. Application

In Part I a general analytical solution for the transient pulse test was presented. Part II presents a graphical method for analyzing data from a test to obtain the hydraulic properties of the sample. The general solution depends on both hydraulic conductivity and specific storage and, in theory, analysis of the data can provide values for both of these hydraulic properties. However, in practice, one of two limiting cases may apply in which case it is possible to calculate only hydraulic conductivity or the product of hydraulic conductivity times specific storage. In this paper we examine the conditions when both hydraulic parameters can be calculated. The analyses of data from two tests are presented. In Appendix I the general solution presented in Part I is compared with an earlier analysis, in which compressive storage in the sample is assumed negligible, and the error in calculated hydraulic conductivity due to this simplifying assumption is examined.

International Journal of Rock Mechanics and Mining↗