Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Fluid Mechanics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources

Three-dimensional seismic structure and moment tensors of non-double-couple earthquakes at the Hengill-Grensdalur volcanic complex, Iceland

The volcanic and geothermal areas of Iceland are rich sources of non-double-couple (non-DC) earthquakes. A state-of-the-art digital seismometer network deployed at the Hengill–Grensdalur volcanic complex in 1991 recorded 4000 small earthquakes. We used the best recorded of these to determine 3-D V P and V P /V S structure tomographically and accurate earthquake moment tensors. The V P field is dominated by high seismic wave speed bodies interpreted as solidified intrusions. A widespread negative (−4 per cent) V P /V S anomaly in the upper 4km correlates with the geothermal field, but is too strong to be caused solely by the effect of temperature upon liquid water or the presence of vapour, and requires in addition mineralogical or lithological differences between the geothermal reservoir and its surroundings. These may be caused by geothermal alteration. Well-constrained moment tensors were obtained for 70 of the best-recorded events by applying linear programming methods to P - and S -wave polarities and amplitude ratios. About 25 per cent of the mechanisms are, within observational error, consistent with DC mechanisms consistent with shear faulting. The other 75 per cent have significantly non-DC mechanisms. Many have substantial explosive components, one has a substantial implosive component, and the deviatoric component of many is strongly non-DC. Many of the non-DC mechanisms are consistent, within observational error, with simultaneous tensile and shear faulting. However, the mechanisms occupy a continuum in source-type parameter space and probably at least one additional source process is occurring. This may be fluid flow into newly formed cracks, causing partial compensation of the volumetric component. Studying non-shear earthquakes such as these has great potential for improving our understanding of geothermal processes and earthquake source processes in general.

Geophysical Journal International

Impact of pore fluid chemistry on fine-grained sediment fabric and compressibility

Fines, defined here as grains or particles, less than 75 μm in diameter, exist nearly ubiquitously in natural sediment, even those classified as coarse. Macroscopic sediment properties, such as compressibility, which relates applied effective stress to the resulting sediment deformation, depend on the fabric of fines. Unlike coarse grains, fines have sizes and masses small enough to be more strongly influenced by electrical interparticle forces than by gravity. These electrical forces acting through pore fluids are influenced by pore fluid chemistry changes. Macroscopic property dependence on pore fluid chemistry must be accounted for in sediment studies involving subsurface flow and sediment stability analyses, as well as in engineered flow situations such as groundwater pollutant remediation, hydrocarbon migration, or other energy resource extraction applications. This study demonstrates how the liquid limit‐based electrical sensitivity index can be used to predict sediment compressibility changes due to pore fluid chemistry changes. Laboratory tests of electrical sensitivity, sedimentation, and compressibility illustrate mechanisms linking microscale and macroscale processes for selected pure, end‐member fines. A specific application considered here is methane extraction via depressurization of gas hydrate‐bearing sediment, which causes a dramatic pore water salinity drop concurrent with sediment being compressed by the imposed effective stress increase.

Journal of Geophysical Research

Observations of static Coulomb stress triggering of the November 2011 M5.7 Oklahoma earthquake sequence

In November 2011, a M 5.0 earthquake occurred less than a day before a M 5.7 earthquake near Prague, Oklahoma, which may have promoted failure of the mainshock and thousands of aftershocks along the Wilzetta fault, including a M 5.0 aftershock. The M 5.0 foreshock occurred in close proximity to active fluid injection wells; fluid injection can cause a buildup of pore fluid pressure, decrease the fault strength, and may induce earthquakes. Keranen et al . [ 2013 ] links the M 5.0 foreshock with fluid injection, but the relationship between the foreshock and successive events has not been investigated. Here we examine the role of coseismic Coulomb stress transfer on earthquakes that follow the M 5.0 foreshock, including the M 5.7 mainshock. We resolve the static Coulomb stress change onto the focal mechanism nodal plane that is most consistent with the rupture geometry of the three M ≥ 5.0 earthquakes, as well as specified receiver fault planes that reflect the regional stress orientation. We find that Coulomb stress is increased, e.g., fault failure is promoted, on the nodal planes of ~60% of the events that have focal mechanism solutions, and more specifically, that the M 5.0 foreshock promoted failure on the rupture plane of the M 5.7 mainshock. We test our results over a range of effective coefficient of friction values. Hence, we argue that the M 5.0 foreshock, induced by fluid injection, potentially triggered a cascading failure of earthquakes along the complex Wilzetta fault system.

Oklahoma

Episodic fluid flow in the Nankai accretionary complex: Timescale, geochemistry, flow rates, and fluid budget

Down-hole geochemical anomalies encountered in active accretionary systems can be used to constrain the timing, rates, and localization of fluid flow. Here we combine a coupled flow and solute transport model with a kinetic model for smectite dehydration to better understand and quantify fluid flow in the Nankai accretionary complex offshore of Japan. Compaction of sediments and clay dehydration provide fluid sources which drive the model flow system. We explicitly include the consolidation rate of underthrust sediments in our calculations to evaluate the impact that variations in this unknown quantity have on pressure and chloride distribution. Sensitivity analysis of steady state pressure solutions constrains bulk and flow conduit permeabilities. Steady state simulations with 30% smectite in the incoming sedimentary sequence result in minimum chloride concentrations at site 808 of 550 m M , but measured chlorinity is as low as 447 m M . We simulate the transient effects of hydrofracture or a strain event by assuming an instantaneous permeability increase of 3–4 orders of magnitude along a flow conduit (in this case the décollement), using steady state results as initial conditions. Transient results with an increase in décollement permeability from 10 −16 m 2 to 10 −13 m 2 and 20% smectite reproduce the observed chloride profile at site 808 after 80–160 kyr. Modeled chloride concentrations are highly sensitive to the consolidation rate of underthrust sediments, such that rapid compaction of underthrust material leads to increased freshening. Pressures within the décollement during transient simulations rise rapidly to a significant fraction of lithostatic and remain high for at least 160 kyr, providing a mechanism for maintaining high permeability. Flow rates at the deformation front for transient simulations are in good agreement with direct measurements, but steady state flow rates are 2–3 orders of magnitude smaller than observed. Fluid budget calculations indicate that nearly 71% of the incoming water in the sediments leaves the accretionary wedge via diffuse flow out the seafloor, 0–5% escapes by focused flow along the décollement, and roughly 1% is subducted.

Journal of Geophysical Research B: Solid Earth

Repeating coupled earthquakes at Shishaldin Volcano, Alaska

Since it last erupted in 1999, Shishaldin Volcano, Aleutian Islands, Alaska, has produced hundreds to thousands of long-period (1-2 Hz; LP) earthquakes every day with no other sign of volcanic unrest. In 2002, the earthquakes also exhibited a short-period (4-7 Hz; SP) signal occurring between 3 and 15 s before the LP phase. Although the SP phase contains higher frequencies than the LP phase, its spectral content is still well below that expected of brittle failure events. The SP phase was never observed without the LP phase, although LP events continued to occur in the absence of the precursory signal. The two-phased events are termed "coupled events", reflecting a triggered relationship between two discrete event types. Both phases are highly repetitive in time series, suggestive of stable, non-destructive sources. Waveform cross-correlation and spectral coherence are used to extract waveforms from the continuous record and determine precise P-wave arrivals for the SP phase. Although depths are poorly constrained, the SP phase is believed to lie at shallow (<4 km) depths just west of Shishaldin's summit. The variable timing between the SP and LP arrivals indicates that the trigger mechanism between the phases itself moves at variable speeds. A model is proposed in which the SP phase results from fluid moving within the conduit, possibly around an obstruction and the LP phase results from the coalescence of a shallow gas bubble. The variable timing is attributed to changes in gas content within the conduit. The destruction of the conduit obstacle on November 21, 2002 resulted in the abrupt disappearance of the SP phase.

Alaska

Constraints on the composition of ore fluids and implications for mineralising events at the Cleo gold deposit, Eastern Goldfields Province, Western Australia

The Cleo gold deposit, 55 km south of Laverton in the Eastern Goldfields Province of Western Australia, is characterised by banded iron-formation (BIF)-hosted ore zones in the gently dipping Sunrise Shear Zone and high-grade vein-hosted ore in the Western Lodes. There is evidence that gold mineralisation in the Western Lodes (which occurred at ca 2655 Ma) post-dates the majority of displacement along the Sunrise Shear Zone, but it remains uncertain if the ore in both structures formed simultaneously or separately. Overall, the Pb, Nd, Sr, C. O and S isotopic compositions of ore-related minerals from both the Western Lodes and ore zones in the Sunrise Shear Zone are similar. Early low-salinity aqueous-carbonic fluids and late high-salinity fluids with similar characteristics are trapped in inclusions in quartz veins from both the Sunrise Shear Zone and the Western Lodes. The early CO2, CO2-H2O, and H2O- dominant inclusions are interpreted as being related to ore formation, and to have formed from a single low-salinity aqueous-carbonic fluid as a result of intermittent fluid immiscibility. Homogenisation temperatures indicate that these inclusions were trapped at approximately 280??C and at approximately 4 km depth, in the deeper epizonal range. Differences between the ore zones are detected in the trace-element composition of gold samples, with gold from the Sunrise Shear Zone enriched in Ni, Pb, Sn, Te and Zn, and depleted In As, Bi, Cd, Cu and Sb, relative to gold from the Western Lodes. Although there are differences in gold composition between the Sunrise Shear Zone and Western Lodes, and hence the metal content of ore fluids may have varied slightly between the different ore zones, no other systematic fluid or solute differences are detected between the ore zones. Given the fact that the ore fluids in each zone have very similar bulk properties, the considerable differences in gold grade, sulfide mineral abundance, and ore textures between the two ore zones most likely result from different gold-deposition mechanisms. The association of ore zones in the Sunrise Shear Zone with pyrite-replaced BIF suggests that wall-rock sulfidation was the most significant mechanism of gold precipitation, through the destabilisation of gold-bisulfide complexes. The Western Lodes, however, do not exhibit any host-rock preference and multistage veins commonly contain coarse-grained gold. Fluid-inclusion characteristics and breccia textures in veins in the Western Lodes suggest that rapid pressure changes, brought about by intermittent release of overpressured fluids and concomitant phase separation, are likely to have caused the destabilisation of gold-thiocomplexes, leading to formation of higher-grade gold ore zones.

Australian Journal of Earth Sciences

Mechanical and chemical compaction in fine-grained shallow-water limestones

Artificial compaction of in-situ cores of sediments resulted in: 1) reduction of sediment thickness by 50 percent and more; 2) reduction of initial porosities of 65 to 75 percent to 35 to 45 percent; 3) creation of megascopic textures almost identical to many ancient lime mud- and wackestone; 4) creation of organic, wispy "stylolite-like" layers; 5) chemical compaction, evidenced by thin sections showing quartz grains piercing mollusc shells without causing fractures and SEM evidence of solutional interfitting of 1 to 4-mu m-size aragonitic carbonate grains; 6) obliteration of pellets and birdseye or fenestral voids in those sediments where early cementation was lacking; obliteration of identifiable marine grasses and vertical "root" tube voids; 8) mashing of sediment-filled circular burrows to produce ellipsoidal structures. Significant mechanical compaction resulted from pressures simulating less than 1,000 ft of burial. Increasing loads to more than 10,000 ft did not significantly increase compaction. Chemical compaction was detected only in cores compacted to pressures greater than 10,000 ft of burial. These experiments suggest that chemical compaction would begin at much shallower depths given geologic time. Experiments that caused chemical compaction lend support to the hypothesis that cement required to produce a low-porosity/low-permeability fine-grained limestone is derived internally. Dissolution, ion diffusion, and reprecipitation are the most likely processes for creating significant thicknesses of dense limestones. Continuation of chemical compaction after significant porosity reduction necessitates expulsion of connate fluids, possibly including hydrocarbons.--Modified journal abstract.

Journal of Sedimentary Petrology

Geologic context of recurring slope lineae in Melas and Coprates Chasmata, Mars

One of the major Mars discoveries of recent years is the existence of recurring slope lineae (RSL), which suggests that liquid water occurs on or near the surface of Mars today. These dark and narrow features emerge from steep, rocky exposures and incrementally grow, fade, and reform on a seasonal basis and are detected in images from the High Resolution Imaging Science Experiment camera. RSL are known to occur at scattered midlatitude and equatorial sites with little spatial connection to one another. One major exception is the steep, low-albedo slopes of Melas and Coprates Chasmata, in Valles Marineris where RSL are detected among diverse geologic surfaces (e.g., bedrock and talus) and landforms (e.g., inselbergs and landslides). New images show topographic changes including sediment deposition on active RSL slopes. Midwall locations in Coprates and Melas appear to have more areally extensively abundant RSL and related fans as compared with other RSL sites found on Mars. Water budget estimates for regional RSL are on the order of 10 5 to 10 6 m 3 of fluid, for depths of 10 to 100mm, and suggest that a significant amount of near-surface watermight be present. Many RSL are concentrated near local topographic highs, such as ridge crests or peaks, which is challenging to explain via groundwater or ice without a recharge mechanism. Collectively, results provide additional support for the notion that significant amounts of near-surface water can be found on Mars today and suggest that a widespread mechanism, possibly related to the atmosphere, is recharging RSL sources.

Journal of Geophysical Research E: Planets

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

Hydraulic sand-model studies of miscible-fluid flow

Hydraulic sand models are useful physical tools in the investigation of the transition zone that occurs between salt and fresh ground water in coastal aquifers. Such models are used to demonstrate the effects of transport mechanisms that influence the dynamic behavior and the shape of the transition zone. The techniques employed in obtaining in-place measurements of solute concentrations are generally the stumbling block in generating data for two-dimensional dispersion systems. Two in-place measurement techniques were used in the studies described: (1) conductivity probes when salt was used as a tracer; and (2) photoelectric cells when organic dye was used as a tracer. Results indicate that conductivity methods are more reliable; however, care must be exercised inasmuch as the probes tend to disturb the fluid flow.

Journal of Research of the U.S. Geological Survey

Thermal conductivity of hydrate-bearing sediments

A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate–saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate–bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

Journal of Geophysical Research B: Solid Earth

Granular-flow rheology: Role of shear-rate number in transition regime

This paper examines the rationale behind the semiempirical formulation of a generalized viscoplastic fluid (GVF) model in the light of the Reiner-Rivlin constitutive theory and the viscoplastic theory, thereby identifying the parameters that control the rheology of granular flow. The shear-rate number ( N ) proves to be among the most significant parameters identified from the GVF model. As N → 0 and N →∞, the GVF model can reduce asymptotically to the theoretical stress versus shear-rate relations in the macroviscous and grain-inertia regimes, respectively, where the grain concentration ( C ) also plays a major role in the rheology of granular flow. Using available data obtained from the rotating-cylinder experiments of neutrally buoyant solid spheres dispersing in an interstitial fluid, the shear stress for granular flow in transition between the two regimes proves dependent on N and C in addition to some material constants, such as the coefficient of restitution. The insufficiency of data on rotating-cylinder experiments cannot presently allow the GVF model to predict how a granular flow may behave in the entire range of N ; however, the analyzed data provide an insight on the interrelation among the relevant dimensionless parameters.

Journal of Engineering Mechanics

P wave crustal velocity structure in the greater Mount Rainier area from local earthquake tomography

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate ( M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

Journal of Geophysical Research B: Solid Earth

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth

Generalized viscoplastic modeling of debris flow

Various concepts have been proposed or used in the development of Theological models for debris flow. The earliest model developed by Bagnold was based on the concept of the “dispersive” pressure generated by grain collisions. Bagnold's concept appears to be theoretically sound, but his empirical model has been found to be inconsistent with most theoretical models developed from non‐Newtonian fluid mechanics. Although the generality of Bagnold's model is still at issue, debris‐flow modelers in Japan have generally accepted Takahashi's formulas derived from Bagnold's model. Some efforts have recently been made by theoreticians in non‐Newtonian fluid mechanics to modify or improve Bagnold's concept or model. A viable rheological model should consist both of a rate‐independent part and a ratedependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i.e., the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical, for general use in debris‐flow modeling. In fact, Bagnold's model is found to be only a particular case of the GVF model. Analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold's simplified assumption of constant grain concentration.

Journal of Hydraulic Engineering

Salinity variations in submarine hydrothermal systems by layered double-diffusive convection

Various mechanisms have been proposed to explain the salinity variations in vent fluids of seafloor geothermal systems. New experiments reacting diabase and evolved seawater were carried out to reproduce earlier published observations of Cl depletions attributed to formation of an ephemeral Cl-bearing mineral. The absence of any Cl depletions in the present study suggests that the formation of Cl-bearing minerals is not sufficiently widespread to account for the observed salinity variations in the vent fluids. A re-evaluation of both field and laboratory evidence has led to a new model for subseafloor circulation that accounts for salinity variations as well as other chemical and mineralogic observations. In place of a simple single-pass convection system, we propose that the seafloor systems consist of two vertically nested convection cells in which a brine layer at depth heats and drives an overlying seawater cell. Such layering of salinities, a process known in fluid mechanics as double-diffusive convection, is an expected result when convection is induced in saline fluids. The process provides for stable high-temperature heat transfer upward from the cracking front adjacent to the magma, and for limited chemical exchange of the brine with the overlying seawater to explain salinity variations and high metal contents in the vent fluids. The brine also provides an effective medium to produce the secondary mineral assemblages observed in rocks from the mid-ocean ridges and ophiolites unsuccessfully produced in laboratory studies using seawater. The brine originates from the two-phase separation of seawater during magmatic/tectonic events and accumulates and remains relatively stable in the region immediately above the magma chamber.

Journal of Geology

Review: The hydrogeology of critical mineral resources relevant to the energy transition

Attaining the goals of the international treaty on climate change (the Paris Agreement) will greatly increase the demand for the critical minerals required to implement clean-energy technologies. This poses both challenges and opportunities to the hydrogeologic community from several perspectives. Here, important insights that the hydrogeological sciences have to offer for mineral exploration, mineral production, and addressing environmental issues related to mining and mine decommissioning are summarized. This study focuses on copper, cobalt, lithium, and rare earths to represent the broad spectrum of critical minerals and illustrate their relevance by referring to projected demands and production rates. The current understanding of the hydrogeologic processes that form major deposits of these minerals are then summarized. Ore is defined as the naturally occurring material from which minerals of economic value can be extracted, where most ore deposits are the products of complex hydrogeologic couplings between fluid flow, heat transport, solute transport, chemical reactions, and mechanical deformation. Exploration models for the discovery of deeper, hidden deposits are potentially informed by hydrogeologic theory and hydrogeochemical processes. Hydrogeologic understanding and methods are also essential to production and recovery. Longstanding challenges are mine dewatering and (conversely) mine water supply, as well as mineral-extraction practices such as spoil heap leaching and in situ mining. New challenges arise from element extraction from subsurface brines. Finally, the quantity of water use and potential environmental impacts of mining on water quality are at the core of ‘social license’: the approval and acceptance of society to mining activities.

Hydrogeology Journal

Hydrological response to earthquakes in the Haibara well, central Japan - II. Possible mechanism inferred from time-varying hydraulic properties

28 coseismic groundwater level decreases have been observed at the Haibara well, Shizuoka prefecture, central Japan, from 1981 to 1997. These groundwater level changes cannot be explained as the poroelastic response to coseismic static strain. We use the atmospheric pressure and tidal responses of the well, rock properties measured on core samples from the same formation and pumping test results to characterize the hydraulic and mechanical properties of the aquifer. The responses of the Haibara well to the M2 Earth tide constituent and to atmospheric pressure have varied over time. In particular, increasing amplitude and decreasing phase lags were observed after the 1993 pumping test, as well as after earthquakes that caused coseismic water level changes. The tidal response, together with the surface load efficiency derived from the atmospheric pressure response, is used to estimate the mechanical properties of the aquifer. The largest amplitude of the M2 constituent, 2.2 mm, is small enough to imply that pore fluid in this system is approximately twice as compressible as water, possibly due to the presence of a small amount of exsolved gas. Diffusion of a coseismic pressure drop near the well could account for the observed time histories of the water level changes. The time histories of the water level drops are well matched by the decay of a coseismic pressure drop at least 80 m away from the well. Removal of a small amount of gas from the formation in that location might in turn explain the coseismic pressure drops.

Geophysical Journal International