SEARCH · Search USGS

Results for “Physics of Fluids”

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

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7Linked to original sources

USGS2010

Geophysical characterization of subsurface properties relevant to the hydrology of the Standard Mine in Elk Basin, Colorado

Geophysical data were collected at the Standard Mine in Elk Basin near Crested Butte, Colorado, to help improve the U.S. Environmental Protection Agency's understanding of the hydrogeologic controls in the basin and how they affect surface and groundwater interactions with nearby mine workings. These data are discussed in the context of geologic observations at the site, the details of which are provided in a separate report. This integrated approach uses the geologic observations to help constrain subsurface information obtained from the analysis of surface geophysical measurements, which is a critical step toward using the geophysical data in a meaningful hydrogeologic framework. This approach combines the benefit of many direct but sparse field observations with spatially continuous but indirect measurements of physical properties through the use of geophysics. Surface geophysical data include: (1) electrical resistivity profiles aimed at imaging variability in subsurface structures and fluid content; (2) self-potentials, which are sensitive to mineralized zones at this site and, to a lesser extent, shallow-flow patterns; and (3) magnetic measurements, which provide information on lateral variability in near-surface geologic features, although there are few magnetic minerals in the rocks at this site. Results from the resistivity data indicate a general two-layer model in which an upper highly resistive unit, 3 to 10 meters thick, overlies a less resistive unit that is imaged to depths of 20 to 25 meters. The high resistivity of the upper unit likely is attributed to unsaturated conditions, meaning that the contact between the upper and lower units may correspond to the water table. Significant lateral heterogeneity is observed because of the presence of major features such as the Standard and Elk fault veins, as well as highly heterogeneous joint distributions. Very high resistivities (greater than 10 kiloohmmeters) are observed in locations that may correspond to more silicified, lower porosity rock. Several thin (2 to 3 meters deep and up to tens of meters wide) low-resistivity features in the very near surface coincide with observed surface-water drainage features at the site. These are limited to depths less than 3 meters and may indicate surface and very shallow groundwater flowing downhill on top of less permeable bedrock. The data do not clearly point to discrete zones of high infiltration, but these cannot be ruled out given the heterogeneous nature of joints in the shallow subsurface. Disseminated and localized electrically conductive mineralization do not appear to play a strong role in controlling the resistivity values, which generally are high throughout the site. The self-potential analysis highlights the Standard fault vein, the northwest (NW) Elk vein near the Elk portal, and several polymetallic quartz veins. These features contain sulfide minerals in the subsurface that form an electrochemical cell that produces their distinct self-potential signal. A smaller component of the self-potential signal is attributed to relatively moderate topographically driven shallow groundwater flow, which is most prevalent in the vicinity of Elk Creek and to a lesser extent in the area of surface-water drainage below the Level 5 portal. Given the anomalies associated with the electrochemical weathering near the Standard fault vein, it is not possible to completely rule out downward infiltration of surface water and shallow groundwater intersected by the fault, though this is an unlikely scenario given the available data. Magnetic data show little variation, consistent with the mostly nonmagnetic host rocks and mineralization at the site, which is verified by magnetic susceptibility measurements and X-ray diffraction mineralogy data on local rock samples. The contact between the Ohio Creek Member of the Mesaverde Formation and Wasatch Formation coincides with a change in character of the magnetic signature, though

Colorado↗

USGS2011

Volcanic plume height measured by seismic waves based on a mechanical model

In August 2008 an unmonitored, largely unstudied Aleutian volcano, Kasatochi, erupted catastrophically. Here we use seismic data to infer the height of large eruptive columns such as those of Kasatochi based on a combination of existing fluid and solid mechanical models. In so doing, we propose a connection between a common, observable, short-period seismic wave amplitude to the physics of an eruptive column. To construct a combined model, we estimate the mass ejection rate of material from the vent on the basis of the plume height, assuming that the height is controlled by thermal buoyancy for a continuous plume. Using the estimated mass ejection rate, we then derive the equivalent vertical force on the Earth through a momentum balance. Finally, we calculate the far-field surface waves resulting from the vertical force. The model performs well for recent eruptions of Kasatochi and Augustine volcanoes if v , the velocity of material exiting the vent, is 120-230 m s -1 . The consistency between the seismically inferred and measured plume heights indicates that in these cases the far-field ~1 s seismic energy radiated by fluctuating flow in the volcanic jet during the eruption is a useful indicator of overall mass ejection rates. Thus, use of the model holds promise for characterizing eruptions and evaluating ash hazards to aircraft in real time on the basis of far-field short-period seismic data. This study emphasizes the need for better measurements of eruptive plume heights and a more detailed understanding of the full spectrum of seismic energy radiated coeruptively.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

Stable isotope evolution and paleolimnology of ancient Lake Creede

The lacustrine carbonate and travertine (tufa) deposits of ancient Lake Creede preserve a remarkable record of the isotopic evolution of the lake. That record indicates that the δ18O of the lake water, and by analogy its salinity, evolved through evaporation. Limited and less reliable data on hydrous minerals and fluid inclusions in early diagenetic carbonates indicate that the δD of the lake waters also evolved through evaporation. The isotope data place restrictions on models of the physical limnology of the lake and its evolution. The closed-basin Lake Creede formed shortly after collapse of the 26.9 Ma Creede caldera. Throughout most of its history it occupied the northern three quarters of the moat between the resurgent dome and wall of the caldera. The Creede Formation was deposited in the basin, dominantly as lacustrine sediments. Travertine mounds interfinger with Creede Formation sediments along the inner and outer margins of the lake basin. An estimated one-half of the original thickness of the Creede Formation has been lost mainly to erosion although scattered remnants of the upper portion remain on the caldera walls. Two diamond core holes (CCM-1 and CCM-2) sampled the uneroded portion of the Creede Formation as part of the U.S. Continental Drilling Program. Volcaniclastic material, including tuff units deposited directly into the lake and ash washed in from the watershed, compose the main lithologies of the Creede Formation. These volcaniclastic strata were produced by episodic ring-fracture volcanism. Lacustrine carbonates make up about 15% of the section sampled by drill core. They occur as 1 mm to 2 cm low-Mg calcite laminae alternating with siliciclastic laminae in scattered intervals throughout the preserved section. The carbonate laminae are accumulations of 5–20 µm crystallites (microsparites) and brine shrimp fecal pellets (peloids) composed mainly of microsparite particles. Low-Mg calcite also occurs as an early diagenetic replacement of gypsum or ikaite (CaCO3 ·6H2O) crystals grown displacively in the muds and silts near the water-sediment interface (rice grains). Other studies indicate that aragonite was the original CaCO3 precipitate forming the microsparite and peloidal laminae and that it converted to calcite during burial diagenesis. Samples from CCM-2 and nearby outcrop do not appear to have undergone significant isotope exchange during recrystallization. Samples from CCM-1 and nearby outcrop, however, appear to have undergone extensive oxygen isotope exchange with meteoric water-dominated fluids possibly during a local 17.6 Ma hydrothermal event. The δ18O-δ13C data set produced by microsampling of individual carbonate lamellae and rice grains is exceptional in several aspects and provides important clues concerning the evolution of limnologic structure of the lake and its chemical and isotopic composition. Travertine and ikaite pseudomorphs in travertine deposits extend the record an additional 330 m above the collar of CCM-2. The δ18O values on CCM-2 samples range from 10.4‰ to 37.3‰ and δ13C values range from –10.8‰ to 9.6‰. The data fall into two distinct groups, a covariant group and an invariant group. The covariant group shows a strong negative covariance and a large range of δ18O and δ13C values. The negative covariance is opposite that normally reported for lacustrine carbonates. The large range of δ18O and δ13C values requires that the carbonates precipitated from waters have a large range of temperature and carbon and oxygen isotopic composition. The invariant group has a narrow range of large δ18O values (35‰ ± 2‰) and a wide range of δ13C values (–10.8‰ to 9.6‰), indicating precipitation from waters with a narrow range of temperature and δ18O but a wide range in δ13C of aqueous carbon. The ranges of isotope values for microsparite and peloid samples are virtually identical; two-thirds are in the covariant group. By contrast, the values for almost all rice grain samples are in the invariant group. The range in δ18O for all samples reflects precipitation from waters having varying proportions of deep, cold evaporated lake water and shallow, warmer meteoric water. The range for δ13C reflects varying proportions of organic carbon and carbon of volcanic or atmospheric origin, probably dominantly volcanic, in the aqueous carbon. Changes in the detailed carbon-oxygen isotope systematics with stratigraphic position define three periods of isotopic evolution of Lake Creede. Period I is represented by the lowest ~200 m of Creede Formation core in CCM-2. Analyses of individual microsparite and peloidal carbonate laminae within single thin sections of samples from that interval are tightly grouped. The data set as a whole shows a negative covariance. Rice grains are not found in this interval. Period II is represented by the succeeding 120 m of core in CCM-2. In that interval, δ13C-δ18O values for individual microsparite and peloidal carbonate laminae within single thin sections show strong negative covariance, and the set of values for the entire interval also shows strong negative covariance. Rice grains occur near the top of the interval. Period III is represented by the upper 225 m of CCM-2 core. In this interval, rice grains are abundant and δ13C-δ18O values for microsparite and peloidal laminae as well as rice grains fall in the invariant group. During Period I the lake was well mixed and the oxygen isotopic composition of the lake in the productive zone was only slightly influenced by short-term (e.g., annual) variations in the water budget of the lake. In Period II the lake was stratified, possibly with annual overturn. The productive zone included the mixolimnion and the isotopic composition of the microsparites and peloids reflected mixtures of shallow surface (meteoric) water containing volcanic or atmospheric CO2 (epilimnion) and cold underlying waters, the oxygen isotopic compositions of which had evolved through evaporation and were dominated by CO2 produced by the oxidation of organic matter (hypolimnion). During Period III the lake remained stratified. The productive zone was in the hypolimnion, probably due to a thinning of the epilimnion resulting from an increase in the surface area of the lake or a decrease in input waters reflecting a climate change. An upsection increase in values of δ18O for the heaviest samples during Periods I and II indicates evaporative concentration of 18O and, by analogy, salinity in the hypolimnion. The δD-δ18O evolution of the lake is inferred on theoretical evaporation trends, comparison to Mono Lake, and measurement of the δD in fluid inclusions in a calcite pseudomorph after ikaite. The δD-δ18O composition of the lake water followed a curved path that eventually hooked over at a nearly constant δ18O value for the lake of 2‰ ± 2‰ Travertine (tufa) mounds formed along the inner and outer margins of the lake in a zone of mixing of warm, volcanic CO2-bearing, meteoric waters and lake water. Ikaite crystals formed on the mounds from unmixed saline lake water, probably below the thermocline. As the position of the thermocline deepened, likely following the spring runoff, the ikaite was replaced by calcite and the resulting “pearls” were covered with travertine deposited from mixed meteoric and lake waters. The upsection increase in δ18O values of the carbonates, the long period of invariance of large δ18OH2O values, the presence of brine shrimp fecal pellets, and the inferred hooked δD-δ18O path are consistent with evidence from other studies that Lake Creede obtained significant salinity rather early in its history and certainly by the time the lake became permanently stratified.

Colorado↗

USGS2014

A depth-averaged debris-flow model that includes the effects of evolving dilatancy. I. Physical basis

To simulate debris-flow behaviour from initiation to deposition, we derive a depth-averaged, two-phase model that combines concepts of critical-state soil mechanics, grain-flow mechanics and fluid mechanics. The model's balance equations describe coupled evolution of the solid volume fraction, m , basal pore-fluid pressure, flow thickness and two components of flow velocity. Basal friction is evaluated using a generalized Coulomb rule, and fluid motion is evaluated in a frame of reference that translates with the velocity of the granular phase, v s . Source terms in each of the depth-averaged balance equations account for the influence of the granular dilation rate, defined as the depth integral of ∇⋅ v s . Calculation of the dilation rate involves the effects of an elastic compressibility and an inelastic dilatancy angle proportional to m − m eq , where m eq is the value of m in equilibrium with the ambient stress state and flow rate. Normalization of the model equations shows that predicted debris-flow behaviour depends principally on the initial value of m − m eq and on the ratio of two fundamental timescales. One of these timescales governs downslope debris-flow motion, and the other governs pore-pressure relaxation that modifies Coulomb friction and regulates evolution of m . A companion paper presents a suite of model predictions and tests.

Proceedings of the Royal Society A↗

USGS2017

New methodology for computing tsunami generation by subaerial landslides: Application to the 2015 Tyndall Glacier landslide, Alaska

Landslide-generated tsunamis pose significant hazards and involve complex, multiphase physics that are challenging to model. We present a new methodology in which our depth-averaged two-phase model D-Claw is used to seamlessly simulate all stages of landslide dynamics as well as tsunami generation, propagation, and inundation. Because the model describes the evolution of solid and fluid volume fractions, it treats both landslides and tsunamis as special cases of a more general class of phenomena. Therefore, the landslide and tsunami can be efficiently simulated as a single-layer continuum with evolving solid-grain concentrations, and with wave generation via direct longitudinal momentum transfer—a dominant physical mechanism that has not been previously addressed in this manner. To test our methodology, we used D-Claw to model a large subaerial landslide and resulting tsunami that occurred on 17 October 2015, in Taan Fjord near the terminus of Tyndall Glacier, Alaska. Modeled shoreline inundation patterns compare well with those observed in satellite imagery.

Alaska↗

USGS2011

Recent faulting in western Nevada revealed by multi-scale seismic reflection

The main goal of this study is to compare different reflection methods used to image subsurface structure within different physical environments in western Nevada. With all the methods employed, the primary goal is fault imaging for structural information toward geothermal exploration and seismic hazard estimation. We use seismic CHIRP (a swept-frequency marine acquisition system), weight drop (an accelerated hammer source), and two different vibroseis systems to characterize fault structure. We focused our efforts in the Reno metropolitan area and the area within and surrounding Pyramid Lake in northern Nevada. These different methods have provided valuable constraints on the fault geometry and activity, as well as associated fluid movement. These are critical in evaluating the potential for large earthquakes in these areas, and geothermal exploration possibilities near these structures.

Nevada↗

USGS2014

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Colloidal and physical transport textures exhibited by electrum and naumannite in bonanza epithermal veins from western USA, and their significance

It is reasonably clear that disequilibrium or “far-from equilibrium” conditions lead to the formation of silica colloids and their deposition in many epithermal deposits. This implies ore-forming solutions had elevated concentrations of dissolved silica, well in excess of amorphous silica saturation. We have previously demonstrated that such colloidal silica particles were deposited in epithermal veins as silica gels and opal, which may later progress along a path to crystallize into more thermodynamically favored (less-soluble) silica phases such as quartz and chalcedony. Also, in some deposits, amorphous silica is co-deposited with precious-metal minerals, such as electrum in the banded super-bonanza ores of the Sleeper deposit (NV). Ore-mineral textures from some western USA bonanza epithermal ores indicate that two precious-metal phases (electrum and naumannite, Ag2Se) form colloidal particles that are transported by ore-forming fluids and are deposited either by aggregation (by sticking to other precious metal-particles) to make dendrites, or are deposited on the “lee” side of protrusion along vein walls (or perhaps by both processes). We can infer by analogy to silica that this also implies that ore-forming solutions contained elevated (supersaturated) dissolved concentrations of both gold and silver that formed colloidal particles under disequilibrium (often chaotic) conditions. Thus physical transport and deposition textures seem to indicate the presence of strongly precious-metal-enriched ore forming fluids, which led to (not surprisingly) the bonanza grades of these remarkable ores. What causes such a precious-metal-rich solution is debatable, but that is the subject of our continued investigations.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2021

Formation of miarolitic-class, segregation-type pegmatites in the Taishanmiao batholith, China: The role of pressure fluctuations and volatile exsolution during pegmatite formation in a closed, isochoric system

The Taishanmiao granitic batholith, located in the Eastern Qinling Orogen in Henan Province, China, contains numerous small (mostly tens of centimeters in maximum dimension) bodies exhibiting textures and mineralogy characteristics of simple quartz and alkali feldspar pegmatites. Analysis of melt inclusions (MI) and fluid inclusions (FI) in pegmatitic quartz, combined with Rhyolite-MELTS modeling of the crystallization of the granite, have been applied to develop a conceptual model of the physical and geochemical processes associated with the formation of the pegmatites. These results allow us to consider the formation of the Taishanmiao pegmatites within the context of varios models that have been proposed for pegmatite formation. Field observations and geochemical data indicate that the pegmatites represent the latest stage in the crystallization of the Taishanmiao granite and occupy ≤4 vol% of the syenogranite phase of the batholith. Results of Rhyolite-MELTS modeling suggest that the pegmatite-forming melts can be produced through continuous fractional crystallization of the Taishanmiao granitic magma, consistent with the designation of the pegmatites as a miarolitic class, segregation-type pegmatites rather than the more common intrusive-type of pegmatite. The mineral assemblage predicted by Rhyolite-MELTS after ~96% of the original granite-forming melt had crystallized consists of ~51 vol% alkali feldspar, 34 vol% quartz, 14 vol% plagioclase, 0.1 vol% biotite, and 1 vol% magnetite, similar to the alkali feldspar + quartz dominated mineralogy of the pegmatites. Moreover, the modeled residual melt composition following crystallization of ~96% of the original melt is similar to the composition of homogenized MI in quartz within the pegmatite. Rhyolite-MELTS predicts that the granite-forming melt remained volatile-undersaturated during crystallization of the batholith and contained ~6.3 wt% H 2 O and ~500 ppm CO 2 after ~96% crystallization when the pegmatites began to develop. The Rhyolite-MELTS prediction that the melt was volatile-undersaturated at the time the pegmatites began to form, but became volatile-saturated during the early stages of pegmatite formation, is consistent with the presence of some inclusion assemblages consisting of only MI, while others contain co-existing MI and FI. The relationship between halogen (F and Cl) and Na abundances in MI is also consistent with the interpretation that the very earliest stages of pegmatite formation occurred in the presence of a volatile-undersaturated melt and that the melt became volatile saturated as crystallization progressed. We propose a closed system, isochoric model for the formation of the pegmatites. Accordingly, the Taishanmiao granite crystallized isobarically at ~3.3 kbar, and the pegmatites began to form at ~734 °C and ~ 3.3 kbar, after ~96% of the original granitic melt had crystallized. During the final stages of crystallization of the granite, small pockets of the remaining residual melt became isolated within the enclosing granite and evolved as constant mass (closed), constant volume (isochoric) systems, similar to the manner in which volatile-rich melt inclusions in igneous phenocrysts evolve during post-entrapment crystallization under isochoric conditions. As a result of the negative volume change associated with crystallization, pressure in the pegmatite initially decreases as crystals form, and this leads to volatile exsolution from the melt phase. The changing PTX conditions produce a pressure-induced “liquidus deficit” that is analogous to liquidus undercooling and results in crystal growth as required to return the system to equilibrium PTX conditions. Owing to the complex closed system, isochoric PVTX evolution of the melt-crystal-volatile system, the pressure does not decrease rapidly or monotonically during pegmatite formation but, rather, gradually fluctuates such that at some stages in the evolution of the pegmatite the pressure is decreasing while at other times the pressure increases as the system cools to maintain mass and volume balance. This behavior, in turn, leads to alternating episodes of precipitation and dissolution that serve to coarsen (ripen) the crystals to produce the pegmatitic texture. The evolution of the pegmatitic melt described here is analogous to that which has been well-documented to occur in volatile-rich MI that undergo closed system, isochoric, post-entrapment crystallization.

American Mineralogist↗

USGS2018

Time-dependent pore filling

Capillarity traps fluids in porous media during immiscible fluid displacement. Most field situations involve relatively long time scales, such as hydrocarbon migration into reservoirs, resource recovery, nonaqueous phase liquid remediation, geological CO 2 storage, and sediment‐atmosphere interactions. Yet laboratory studies and numerical simulations of capillary phenomena rarely consider the impact of time on these processes. We use time‐lapse microphotography to record the evolution of saturation in air‐ or hydrocarbon‐filled capillary tubes submerged in water to investigate long‐term pore filling phenomena beyond imbibition. Microphotographic sequences capture a lively pore filling history where various concurrent physical phenomena coexist. Dissolution and diffusion play a central role. Observations indicate preferential transport of the wetting liquid along corners, vapor condensation, capillary flow induced by asymmetrical interfaces, and interface pinning that defines the diffusion length. Other processes include internal snap‐offs, fluid redistribution, and changes in wettability as fluids dissolve into each other. Overall, the rate of pore filling is diffusion‐controlled for a given interfacial configuration; diffusive transport takes place at a constant rate for pinned interfaces and is proportional to the square root of time for free interfaces where the diffusion length increases with time.

Water Resources Research↗

USGS1971

Problems on the origin of ore deposits in the lower Ordovician formations of east Tennessee

Most recent workers in the East Tennessee zinc and barite districts are in general agreement that the host breccias were formed by solution-collapse processes in early Middle Ordovician time, probably in an ancient carbonate aquifer system of regional extent. There is little agreement, however, on the source and nature of the ore fluids from which the epigenetic minerals were formed. In this symposium Hanshaw and his colleagues have presented evidence in support of the concept that much secondary dolomite may have been formed by ground water. On the other hand, Roedder has shown that fluid inclusion studies indicate hot saline brines (75 degrees -150 degrees C) as the source fluids for the ore minerals, and he suggests that they may have been deeply circulating connate brines. Other investigators in the past have proposed a variety of source fluids, including juvenile waters, meteoric waters that have become enriched by the leaching of particular elements from the surrounding rock, and sea water. In the zinc deposits, sphalerite is accompanied locally by minor amounts of galena and pyrite, and occasionally, traces of chalcopyrite. Dolomite is the chief gangue mineral and is associated with lesser amounts of quartz, calcite, barite, fluorite, anhydrite, gypsum, and bitumen. Some deposits containing barite as the dominant ore mineral have various amounts of fluorite and pyrite; sphalerite and galena may be present locally along with the other gangue minerals typical of the zinc deposits. The paragenetic sequences of the minerals have been determined for most of the mines and many of the prospects, and, except for minor refinements, will be little changed by additional studies. The real problem, though, is the correlation of these depositional sequences from mine to mine and from district to district. Perhaps the major difficulty in the correlation is the seemingly uncomplicated mineralogy that disguises a rather complex regional paragenesis. Questions regarding the routes the ore fluids travelled to reach the sites of deposition, the reasons why the ores are concentrated in some areas more than others, and the reasons why the mineralogy varies from area to area are questions that are unanswered or at best only partly answered. Studies to date have produced convincing evidence concerning the nature of the structures and the physical conditions responsible for the localization of ore in the Kingsport Formation (Lower Ordovician). On the other hand, the relation of the Kingsport ores to the smaller deposits in the underlying Lower Ordovician and Cambrian rocks is but little understood, and no final answer to the origin of these deposits can be reached until this relationship has been satisfactorily explained. Fluid inclusion studies, trace element analysis, and isotopic studies of all sorts are needed to further our understanding of these seemingly simple ores. Continued detailed studies of field relations, both regional and local, are of major importance, not only for what they will reveal in themselves, but also to perfect the frames of reference in which to place the growing mass of laboratory data.

Tennessee↗

USGS1984

Hydrogen speciation in synthetic quartz

The dominant hydrogen impurity in synthetic quartz is molecular H2O. H-OH groups also occur, but there is no direct evidence for the hydrolysis of Si-O-Si bonds to yield Si-OH HO-Si groups. Molecular H2O concentrations in the synthetic quartz crystals studied range from less than 10 to 3,300 ppm (H/Si), and decrease smoothly by up to an order of magnitude with distance away from the seed. OH- concentrations range from 96 to 715 ppm, and rise smoothly with distance away from the seed by up to a factor of three. The observed OH- is probably all associated with cationic impurities, as in natural quartz. Molecular H2O is the dominant initial hydrogen impurity in weak quartz. The hydrolytic weakening of quartz may be caused by the transformation H2O + Si-O-Si ??? 2SiOH, but this may be a transitory change with the SiOH groups recombining to form H2O, and the average SiOH concentration remaining very low. Synthetic quartz is strengthened when the H2O is accumulated into fluid inclusions and cannot react with the quartz framework. ?? 1984 Springer-Verlag.

Physics and Chemistry of Minerals↗

USGS2007

The role of fault zone drilling

The objective of fault-zone drilling projects is to directly study the physical and chemical processes that control deformation and earthquake generation within active fault zones. An enormous amount of field, laboratory, and theoretical work has been directed toward the mechanical and hydrological behavior of faults over the past several decades. Nonetheless, it is currently impossible to differentiate between – or even adequately constrain – the numerous conceptual models of active faults proposed over the years. For this reason, the Earth science community is left in the untenable position of having no generally accepted paradigm for the mechanical behavior of faults at depth. One of the primary causes for this dilemma is the difficulty of either directly observing or inferring physical properties and deformation mechanisms along faults at depth, as well as the need to observe directly key parameters such as the state of stress acting on faults at depth, pore fluid pressure (and its possible variation in space and time), and processes associated with earthquake nucleation and rupture. Today, we know very little about the composition of active faults at depth, their constitutive properties, the state of in situ stress or pore pressure within fault zones, the origin of fault-zone pore fluids, or the nature and significance of time-dependent fault-zone processes.

Book chapter↗

USGS2025

Potential for continental scientific drilling to inform fault mechanics and earthquake science

Our understanding of fault mechanics and earthquake processes remains limited, largely due to minimal direct observations near active faults at seismogenic depths. This lack of data restricts our ability to accurately assess and mitigate both natural and human-induced seismic hazards. However, recent advancements in drilling capabilities and downhole sensing technologies offer an opportunity: the ability to observe the physical conditions within a volume near active fault zones. In this contribution, we highlight how scientific drilling can provide access to the near-fault environment, enabling measurements of the stress, temperature, fluid pressure, and rock properties at depths where ruptures initiate, propagate, and arrest. These observations are essential to refine models of earthquake nucleation and dynamic rupture, bridging gaps between laboratory experiments, numerical simulations, and surface observations. These insights can advance fundamental understanding in earthquake science but also support the development of more effective seismic hazard assessments and risk mitigation strategies.

Seismica↗

USGS2020

Impacts of hydrothermal plume processes on oceanic metal cycles and transport

Chemical, physical and biological processes in hydrothermal plumes control the flux of elements from hydrothermal vents to the global oceans. The timescales of these processes range from less than a second, as the hydrothermal fluid mixes with seawater at the seafloor, to decades, as the plume disperses over thousands of kilometers. Integrating hydrothermal geochemistry throughout the lifetime of the plume reveals some well constrained processes, along with many surprises. For instance, contrary to the idea that metals are removed from the hydrothermal plume via oxidation, a survey of recent datasets reveals that oxidation of iron and manganese does not consistently result in their removal from the plume, and that manganese may be lost from the water column more rapidly than iron. These observations suggest that the understanding of element transport in hydrothermal plumes is incomplete, partly due to the change in removal processes as the plume disperses from less than 1 km from the vent to more than 4,000 km. We suggest that characterizing the plume based on regions that retain some reduced components versus those that are fully oxidized, in addition to buoyancy, will illuminate the nature of the dominant processes and allow a more complete understanding of the ultimate fate of hydrothermally derived metals.

Nature Geoscience↗

USGS2015

Hydrothermal response to a volcano-tectonic earthquake swarm, Lassen, California

The increasing capability of seismic, geodetic, and hydrothermal observation networks allows recognition of volcanic unrest that could previously have gone undetected, creating an imperative to diagnose and interpret unrest episodes. A November 2014 earthquake swarm near Lassen Volcanic National Park, California, which included the largest earthquake in the area in more than 60 years, was accompanied by a rarely observed outburst of hydrothermal fluids. Although the earthquake swarm likely reflects upward migration of endogenous H 2 O-CO 2 fluids in the source region, there is no evidence that such fluids emerged at the surface. Instead, shaking from the modest sized (moment magnitude 3.85) but proximal earthquake caused near-vent permeability increases that triggered increased outflow of hydrothermal fluids already present and equilibrated in a local hydrothermal aquifer. Long-term, multiparametric monitoring at Lassen and other well-instrumented volcanoes enhances interpretation of unrest and can provide a basis for detailed physical modeling.

California↗

USGS2001

Updip limit of the seismogenic zone beneath the accretionary prism of southwest Japan: An effect of diagenetic to low-grade metamorphic processes and increasing effective stress

Off southwest Japan the seaward limit of coseismic displacement (or updip limit of the seismogenic zone) of the 1946 M w 8.3 thrust earthquake reaches to 4 km depth and ∼40 km landward of the trench. This limit coincides with the estimated location of the 150 °C isotherm, and has been linked to changes in physical properties associated with the smectite to illite clay-mineral transition. Here we show that this limit correlates with a suite of diagenetic to low-grade metamorphic processes characterized by (1) declining fluid production and decreasing fluid pressure ratio (λ * ) and (2) active clay, carbonate, and zeolite cementation and the transition to pressure solution and quartz cementation. These diagenetic to low-grade metamorphic changes cause the onset of velocity weakening during thrust faulting, an increase in effective stress, and strengthening of the hanging wall, which together combine to produce recordable earthquakes.

Geology↗

USGS2022

Detection of aseismic slip and poroelastic reservoir deformation at the North Brawley Geothermal Field from 2009 to 2019

The North Brawley Geothermal Field, located within the Brawley Seismic Zone of Southern California, presents a case study for understanding seismic hazards linked to fluid injection and geothermal energy extraction. An earthquake swarm near the geothermal field in 2012 included two earthquakes with magnitudes greater than 5 and was potentially preceded by a years-long aseismic slip transient. To better understand ground deformation around the geothermal field, including its evolution with time and its physical mechanisms, we analyze deformation before, during, and after the swarm using ground- and satellite-based geodetic techniques between 2009 and 2019. We integrate observations from GNSS, Sentinel-1, TerraSAR-X, UAVSAR, and leveling surveys into a single deformation history. Modeling of this new collection of observations at the North Brawley Geothermal Field provides evidence for 80% more pre-swarm aseismic slip than previously recognized from 2009 to 2012. During the 2012 Brawley swarm, our geodetic slip inversions closely match the results of seismic waveform inversions from the swarm events. After the 2012 swarm, surface deformation is dominated by poroelastic deformation of a shallow fluid reservoir at <1 km depth rather than fault slip. The deformation history and seismicity catalogs at North Brawley suggest a cessation of fault-related slip during the ∼7 years after the 2012 earthquake swarm.

California↗