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 145 records · Page 8Linked to original sources

USGS1998

Reply to comment by E. Holzbecher on "Constant-concentration boundary condition: Lessons from the HYDROCOIN variable-density groundwater benchmark problem"

We appreciate the opportunity to continue a dialogue on solute-transport modeling issues and perhaps to clarify some of our original concerns. It appears to us that Holzbecher [this issue] fails to consider the physical implications of groundwater flow and solute transport along no-flow boundaries. A no-flow boundary coincides with a streamline along which the fluid velocity is generally nonzero. Because the velocity is nonzero, both advective and dispersive transport processes can be operative along a no-flow boundary.

Water Resources Research↗

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↗

USGS1956

Diatremes and craters attributed to natural explosions

Diatremes - volcanic pipes attributed to explosion - and craters have been studied to infer the ultimate causes and physical conditions attending natural explosive processes. Initial piercement of diatremes on the Navajo reservation, Arizona was probably along a fracture propagated by a high-pressure aqueous fluid. Gas rising at high velocity along the fracture would become converted to a gas-solid fluidized system by entrainment of wall- rock fragments. The first stages of widening of the vent are probably accomplished mainly by simple abrasion of the high-velocity fluidized system on the walls of the fracture. As the vent widens, its enlargement may be accelerated by inward spalling of the walls. The inferred mechanics of the Navajo-Hopi diatremes is used to illustrate the possibility of diatreme formation over a molten salt mass.

Trace Elements Memorandum↗

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↗

USGS2005

Generation and verification of theoretical rating curves in the Whitewater River basin, Kansas

[1] A new method for generating stage-discharge relations (rating curves) for geomorphically stable channels is presented and applied to two streams in the Whitewater River basin, Kansas. The approach converts measurements of stage into discharge using a fluid mechanically based model. The model does not use empirical roughness coefficients, such as Manning coefficients, but rather determines channel roughness from field measurements of the (1) channel geometry, (2) the physical roughness of the bed, banks, and floodplain, and (3) the vegetation density on the banks and floodplain. These measurements are used to calculate explicitly the drag on the small-scale topographic features on the boundary, the drag on the stems and branches of woody vegetation, and the friction on the bed, banks, and floodplain. The theoretical rating curves produced by the model for two study reaches, which are near U.S. Geological Survey (USGS) streamflow-gauging stations, are in good agreement with direct measurements of discharge made by the USGS. Our method has the potential of providing accurate estimates of stream flows less expensively than conventional gauging methods. In addition, the method can be used to obtain more accurate discharge estimates than conventional indirect methods for determining discharge, which are based on estimates of Manning's coefficient.

Journal of Geophysical Research F: Earth Surface↗

USGS1982

Field, model, and computer simulation study of some aspects of the origin and distribution of Colorado Plateau-type uranium deposits

Numerous hypotheses have been proposed to account for the nature and distribution of tabular uranium and vanadium-uranium deposits of the Colorado Plateau. In one of these hypotheses it is suggested that the deposits resulted from geochemical reactions at the interface between a relatively stagnant groundwater solution and a dynamic, ore-carrying groundwater solution which permeated the host sandstones (Shawe, 1956; Granger, et al., 1961; Granger, 1968, 1976; and Granger and Warren, 1979). The study described here was designed to investigate some aspects of this hypothesis, particularly the nature of fluid flow in sands and sandstones, the nature and distribution of deposits, and the relations between the deposits and the host sandstones. The investigation, which was divided into three phases, involved physical model, field, and computer simulation studies. During the initial phase of the investigation, physical model studies were conducted in porous-media flumes. These studies verified the fact that humic acid precipitates could form at the interface between a humic acid solution and a potassium aluminum sulfate solution and that the nature and distribution of these precipitates were related to flow phenomena and to the nature and distribution of the host porous-media. During the second phase of the investigation field studies of permeability and porosity patterns in Holocene stream deposits were investigated and the data obtained were used to design more realistic porous media models. These model studies, which simulated actual stream deposits, demonstrated that precipitates possess many characteristics, in terms of their nature and relation to host sandstones, that are similar to ore deposits of the Colorado Plateau. The final phase of the investigation involved field studies of actual deposits, additional model studies in a large indoor flume, and computer simulation studies. The field investigations provided an up-to-date interpretation of the depositional environments of the host sandstones in the Slick Rock District and data on the nature and distribution of the ore deposits which are found to be directly related to the architecture of the host sandstones which acted as conduits for the transport of mineralized groundwaters. Large-scale model studies, designed to simulate Grants Mineral Belt deposits, demonstrated that precipitates had characteristics similar to those of actual uranium deposits and data obtained from these studies strongly supported the hypothesis that the ores formed soon after deposition of the host sandstones and that their distribution was largely controlled by permeability and porosity patterns established at the time of deposition of the host sandstones. A numerical model was developed during the second and third stages of the investigation that can predict favorable locations for mineralization given sufficient data on porosity, hydraulic conductivity, the distribution and thickness of sandstone hosts, and an estimate of the initial hydrologic conditions. The model was successfully tested using data from the Slick Rock District.

Open-File Report↗

USGS2009

Incorporating both physical and kinetic limitations in quantifying dissolved oxygen flux to aquatic sediments

Traditionally, dissolved oxygen (DO) fluxes have been calculated using the thin-film theory with DO microstructure data in systems characterized by fine sediments and low velocities. However, recent experimental evidence of fluctuating DO concentrations near the sediment-water interface suggests that turbulence and coherent motions control the mass transfer, and the surface renewal theory gives a more mechanistic model for quantifying fluxes. Both models involve quantifying the mass transfer coefficient (k) and the relevant concentration difference (?? C ). This study compared several empirical models for quantifying k based on both thin-film and surface renewal theories, as well as presents a new method for quantifying ?? C (dynamic approach) that is consistent with the observed DO concentration fluctuations near the interface. Data were used from a series of flume experiments that includes both physical and kinetic uptake limitations of the flux. Results indicated that methods for quantifying k and ?? C using the surface renewal theory better estimated the DO flux across a range of fluid-flow conditions. ?? 2009 ASCE.

Journal of Environmental Engineering↗

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↗