Search USGS⌕ Search

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.

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

At least 91 records · Page 5Linked to original sources
USGS2009

Water, something peculiar

Some chemical and physical properties of water are discussed and compared with those of other fluids. For instance, the boiling point is much higher than one would expect considering the molecular weight of water. The heat capacity is also much higher but the viscosity is not. The dielectric constant is exceptionally high. These and other properties of water can be explained by the geometry of the water molecule and the structure of water or ice.

Hydrological Sciences Bulletin↗
USGS2021

Hydrothermal alteration can result in pore pressurization and volcano instability

The collapse of a volcanic flank can be destructive and deadly. Hydrothermal alteration is common to volcanoes worldwide and is thought to promote volcano instability by decreasing rock strength. However, some laboratory studies have shown that not all alteration reduces rock strength. Our new laboratory data for altered rhyodacites from Chaos Crags (Lassen volcanic center, California, USA) show that pore- and crack-filling mineral precipitation can reduce porosity and permeability and increase strength, Young's modulus, and cohesion. A significant reduction in permeability, by as much as four orders of magnitude, will inhibit fluid circulation and create zones of high pore fluid pressure. We explored the consequences of pore fluid pressurization on volcano stability using large-scale numerical modeling. Upscaled physical and mechanical properties for hydrothermally altered rocks were used as input parameters in our modeling. Results show that a high-pore-pressure zone within a volcano increases volcano deformation and that increasing the size of this zone increases the observed deformation. Hydrothermal alteration associated with mineral precipitation, and increases to rock strength, can therefore promote pore pressurization and volcano deformation, increasing the likelihood of volcano spreading, flank collapses, and phreatic/phreatomagmatic explosions. We conclude that porosity-decreasing alteration, explored here, and porosity-increasing alteration can both promote volcano instability and collapse, but by different mechanisms. Hydrothermal alteration should therefore be monitored at volcanoes worldwide and incorporated into hazard assessments.

Geology↗
USGS1988

Borehole geophysics applied to ground-water investigations

The purpose of this manual is to provide hydrologists, geologists, and others who have the necessary training with the basic information needed to apply the most useful borehole-geophysical-logging techniques to the solution of problems in ground-water hydrology. Geophysical logs can provide information on the construction of wells and on the character of the rocks and fluids penetrated by those wells, in addition to changes in the character of these factors with time. The response of well logs is caused by: petrophysical factors; the quality; temperature, and pressure of interstitial fluids; and ground-water flow. Qualitative and quantitative analysis of the analog records and computer analysis of digitized logs are used to derive geohydrologic information. This information can then be extrapolated vertically within a well and laterally to other wells using logs. The physical principles by which the mechanical and electronic components of a logging system measure properties of rocks, fluids and wells, and the principles of measurement need to be understood to correctly interpret geophysical logs. Planning the logging operation involves selecting the equipment and the logs most likely to provide the needed information. Information on well construction and geohydrology are needed to guide this selection. Quality control of logs is an important responsibility of both the equipment operator and log analyst and requires both calibration and well-site standardization of equipment. Logging techniques that are widely used in ground-water hydrology or that have significant potential for application to this field include: spontaneous potential, resistance, resistivity, gamma, gamma spectrometry, gamma-gamma, neutron, acoustic velocity, acoustic televiewer, caliper, and fluid temperature, conductivity, and flow. The following topics are discussed for each of these techniques: principles and instrumentation, calibration and standardization, volume of investigation, extraneous effects, and interpretation and applications.

Open-File Report↗
USGS1990

Borehole geophysics applied to ground-water investigations

The purpose of this manual is to provide hydrologists, geologists, and others who have the necessary background in hydrogeology with the basic information needed to apply the most useful borehole-geophysical-logging techniques to the solution of problems in ground-water hydrology. Geophysical logs can provide information on the construction of wells and on the character of the rocks and fluids penetrated by those wells, as well as on changes in the character of these factors over time. The response of well logs is caused by petrophysical factors, by the quality, temperature, and pressure of interstitial fluids, and by ground-water flow. Qualitative and quantitative analysis of analog records and computer analysis of digitized logs are used to derive geohydrologic information. This information can then be extrapolated vertically within a well and laterally to other wells using logs. The physical principles by which the mechanical and electronic components of a logging system measure properties of rocks, fluids, and wells, as well as the principles of measurement, must be understood if geophysical logs are to be interpreted correctly. Plating a logging operation involves selecting the equipment and the logs most likely to provide the needed information. Information on well construction and geohydrology is needed to guide this selection. Quality control of logs is an important responsibility of both the equipment operator and the log analyst and requires both calibration and well-site standardization of equipment. Logging techniques that are widely used in ground-water hydrology or that have significant potential for application to this field include spontaneous potential, resistance, resistivity, gamma, gamma spectrometry, gamma-gamma, neutron, acoustic velocity, acoustic televiewer, caliper, and fluid temperature, conductivity, and flow. The following topics are discussed for each of these techniques: principles and instrumentation, calibration and standardization, volume of investigation, extraneous effects, and interpretation and applications.

Techniques of Water-Resources Investigations↗
USGS2011

An approach to modeling coupled thermal-hydraulic-chemical processes in geothermal systems

Interactions between hydrothermal fluids and rock alter mineralogy, leading to the formation of secondary minerals and potentially significant physical and chemical property changes. Reactive transport simulations are essential for evaluating the coupled processes controlling the geochemical, thermal and hydrological evolution of geothermal systems. The objective of this preliminary investigation is to successfully replicate observations from a series of hydrothermal laboratory experiments [Morrow et al., 2001] using the code TOUGHREACT. The laboratory experiments carried out by Morrow et al. [2001] measure permeability reduction in fractured and intact Westerly granite due to high-temperature fluid flow through core samples. Initial permeability and temperature values used in our simulations reflect these experimental conditions and range from 6.13 × 10−20 to 1.5 × 10−17 m2 and 150 to 300 °C, respectively. The primary mineralogy of the model rock is plagioclase (40 vol.%), K-feldspar (20 vol.%), quartz (30 vol.%), and biotite (10 vol.%). The simulations are constrained by the requirement that permeability, relative mineral abundances, and fluid chemistry agree with experimental observations. In the models, the granite core samples are represented as one-dimensional reaction domains. We find that the mineral abundances, solute concentrations, and permeability evolutions predicted by the models are consistent with those observed in the experiments carried out by Morrow et al. [2001] only if the mineral reactive surface areas decrease with increasing clay mineral abundance. This modeling approach suggests the importance of explicitly incorporating changing mineral surface areas into reactive transport models.

Conference Paper↗
USGS2012

GFDL's ESM2 global coupled climate-carbon Earth System Models. Part I: physical formulation and baseline simulation characteristics

We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory's previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.

Journal of Climate↗
USGS2018

Sulfur cycle in the Valles Caldera volcanic complex, New Mexico – Letter 1: Sulfate sources in aqueous system, and implications for S isotope record in Gale Crater on Mars

Initial in situ sulfur (S) isotope measurements of the Martian bedrock in Gale Crater have revealed an unexpectedly wide range of δ 34 S values (−47 to +28%). Generally, it is unclear what processes could have contributed to these large isotope fractionations. Therefore, we studied S sources and aqueous SO 2− 4 cycling in the Valles Caldera volcanic complex, New Mexico to better understand S isotope fractionations related to S degassing, hydrothermal activity, and low-temperature processes in aqueous environment. Overall, our study demonstrates that volcanic systems show large spatial heterogeneity in δ 34 S. Magmatic S sources are obvious in steam-dominated H 2 S degassing and precipitation of secondary minerals from hydrothermal fluids with low δ 34 S values of +0.9 ±3%. Locally, however, hydrothermal processes have resulted in more negative δ 34 S values in sulfide minerals (−18 to −4%) and more positive δ 34 S values in sulfate minerals (−1 to +3%). Major aqueous SO 2− 4 sources are oxidation of H 2 S from modern hydrothermal gas emission, and oxidation and dissolution of sulfide and sulfate minerals present in the hydrothermally altered bedrock and crater-lake sediments. The δ 34 S of aqueous SO 2− 4 in surface water and groundwater varies widely (−8 to +5%) and is similar to major S endmembers that undergo oxidation and/or dissolution by active hydrological system. Minor SO 2− 4 contributions with more positive δ 34 S values (+9 to +14%) come from deeply circulating geothermal fluids and negligible amounts from atmospheric deposition (+5 to +7% in snow). Elevated SO 2− 4 contents are mainly associated with modern and past H 2 S emissions and oxidations near the surface. On regional scale, however, most of the intracaldera bedrock is S-depleted, thus the SO 2− 4 contents are usually low in the surface aquatic system and younger sedimentary lake deposits formed at times of negligible near surface hydrothermal activity. In general, magmatic-hydrothermal processes apparently cause the largest δ 34 S variation in S-bearing minerals on volcanic terrains. Therefore, we infer that the measured wide range of δ 34 S values in the Gale sediments by the Curiosity rover on Mars can be explained by S isotope composition of magmatic-hydrothermal sulfide and sulfate minerals that were present in the initial igneous/volcanic rocks prior to crater formation. Later aqueous processes involved oxidation and dissolution of S minerals initially present in these rocks and led to subsequent formation of diagenetic fluids and alteration products enriched in SO 2− 4 with relatively large δ 34 S variation. Additionally, physical erosion, transport and deposition of detrital hydrothermal S minerals from igneous/volcanic rocks might be in part responsible for the measured wide range of δ 34 S in Gale Crater. These unique S isotope results, measured in situ on another planet for the first time, imply the importance of magmatic-hydrothermal fluids in S transport on early Mars and their subsequent alteration in low-temperature aqueous environments.

New Mexico↗
USGS2023

Multiphysics modelling in PyLith: Poroelasticity

PyLith, a community, open-source code for modelling quasi-static and dynamic crustal deformation with an emphasis on earthquake faulting, has recently been updated with a flexible multiphysics implementation. We demonstrate the versatility of the multiphysics implementation by extending the code to model fully coupled continuum poromechanics. We verify the newly incorporated physics using standard benchmarks for a porous medium saturated with a slightly compressible fluid. The benchmarks include the 1-D consolidation problem as outlined by Terzaghi, Mandel’s problem for the 2-D case, and Cryer’s problem for the 3-D case. All three benchmarks have been added to the PyLith continuous integration test suite. We compare the closed form analytical solution for each benchmark against solutions generated by our updated code, and lastly, demonstrate that the poroelastic material formulation may be used alongside the existing fault implementation in PyLith.

Geophysical Journal International↗
USGS1978

A theoretical basis for exploration for native copper in northern Wisconsin

Exploration for native copper in the Keweenawan lavas of northern Wisconsin has been concentrated in areas of relatively shallow overburden that have sparse to numerous outcrops. Lack of success of this exploration suggests that if large deposits, comparable to those of northern Michigan, are present, they are more likely to be found in one or more of the large tracts that have few, if any, exposures, away from the ' copper ranges. ' A hydrologic model that may explain the mineralization of the classic native-copper district of Michigan could be helpful in suggesting that certain covered tracts are more favorable than others, and in narrowing the targets for physical exploration within these tracts. This model involves updip migration of a hydrothermal fluid of metamorphic origin, formed when ground water contained in the interstices of lava flows and conglomerate beds was carried to great depth along the axis of the Lake Superior syncline. The largest reservoirs of such buried ground water would be expected in the peripheral rather than central parts of the basins in which the lavas originally accumulated because the ratio of porous fragmental flow tops and conglomerate beds to massive basalt should increase towards the margins of a basin. Thus, the areas most promising for copper deposits should be those updip from places where the trough of the later formed Lake Superior syncline crosses the marginal parts of a lava basin. (Woodard-USGS)

Circular↗
USGS2006

Is amplitude loss of sonic waveforms due to intrinsic attenuation or source coupling to the medium?

Sonic waveforms acquired in gas-hydrate-bearing sediments indicate strong amplitude loss associated with an increase in sonic velocity. Because the gas hydrate increases sonic velocities, the amplitude loss has been interpreted as due to intrinsic attenuation caused by the gas hydrate in the pore space, which apparently contradicts conventional wave propagation theory. For a sonic source in a fluid-filled borehole, the signal amplitude transmitted into the formation depends on the physical properties of the formation, including any pore contents, in the immediate vicinity of the source. A signal in acoustically fast material, such as gas-hydrate-bearing sediments, has a smaller amplitude than a signal in acoustically slower material. Therefore, it is reasonable to interpret the amplitude loss in the gas-hydrate-bearing sediments in terms of source coupling to the surrounding medium as well as intrinsic attenuation. An analysis of sonic waveforms measured at the Mallik 5L-38 well, Northwest Territories, Canada, indicates that a significant part of the sonic waveform's amplitude loss is due to a source-coupling effect. All amplitude analyses of sonic waveforms should include the effect of source coupling in order to accurately characterize the formation's intrinsic attenuation.

Scientific Investigations Report↗
USGS2012

Book review: Earthquakes and water

It is really nice to see assembled in one place a discussion of the documented and hypothesized hydrologic effects of earthquakes. The book is divided into chapters focusing on particular hydrologic phenomena including liquefaction, mud volcanism, stream discharge increases, groundwater level, temperature and chemical changes, and geyser period changes. These hydrologic effects are inherently fascinating, and the large number of relevant publications in the past decade makes this summary a useful milepost. The book also covers hydrologic precursors and earthquake triggering by pore pressure. A natural need to limit the topics covered resulted in the omission of tsunamis and the vast literature on the role of fluids and pore pressure in frictional strength of faults. Regardless of whether research on earthquake-triggered hydrologic effects ultimately provides insight into the physics of earthquakes, the text provides welcome common ground for interdisciplinary collaborations between hydrologists and seismologists. Such collaborations continue to be crucial for investigating hypotheses about the role of fluids in earthquakes and slow slip. Review info: Earthquakes and Water. By Wang, C.-Y. and Manga, M., 2010. ISBN: 9783642008092, 218 pp.

Geofluids↗
USGS1974

Alteration and fluid inclusion studies of the porphyry copper ore body at Bingham, Utah

Distribution patterns for biotitic alteration, sericitic alteration, and distinctive fluid-inclusion types in igneous host rocks of the porphyry copper ore body at Bingham, Utah, have been determined by petrographic examination of about 300 samples. These patterns are related to differences in original rock composition, variations in physical-chemical conditions during periods of intrusion and mineralization, and spatial position within the ore body. The distribution of biotitic (potassium-silicate) alteration assemblages and high-salinity fluid inclusions generally follows the crudely triangular form of the disseminated copper ore zone. Variations in abundance of hydrothermal biotite are attributed to differences in original mafic mineral content of the igneous host rocks. Biotitic alteration and initial copper mineralization were accomplished by high-salinity fluids concentrated during final crystallization of the monzonitic parent magma; genetic continuity between magmatic and hydrothermal stages is indicated . Pervasive sericitic alteration of plagioclase is confined to rocks in the northern one-third of the Bingham stock; a subzone of argillic alteration in the north-central part of the ore body occurs within the broader area of sericitic alteration. Fluids responsible for sericitic and argillic alteration were channeled by a broad zone of northeast-trending fractures. Hydrothermal minerals and high-salinity fluid inclusions occur within a large volume of shattered rock. Boiling of fluids during crystallization of the aplitic porphyry may account for the shattering. Sericitic (and argillic) alteration were apparently super-imposed on the earlier biotitic assemblage as the hydrothermal system cooled. Cooling and hydrolytic alteration were promoted by progressive introduction of meteoric waters. The many generations of inclusions trapped from boiling fluids in the temperature range 400 degrees to 600 degrees C suggest that the system was recharged repeatedly during the period of mineralization. Estimated fluid pressures of about 800 bars in the early stages of mineralization correspond to a lithostatic load of about 3 km; pressures were even lower (less than 200 bars) in the later stages and were probably controlled by hydrostatic conditions. Distribution patterns for biotitic alteration, sericitic alteration, and distinctive fluid-inclusion types in igneous host rocks of the porphyry copper ore body at Bingham, Utah, have been determined by petrographic examination of about 300 samples. These patterns are related to differences in original rock composition, variations in physical-chemical conditions during periods of intrusion and mineralization, and spatial position within the ore body.The distribution of biotitic (potassium-silicate) alteration assemblages and high-salinity fluid inclusions generally follows the crudely triangular form of the disseminated copper ore zone. Variations in abundance of hydrothermal biotite are attributed to differences in original mafic mineral content of the igneous host rocks. Biotitic alteration and initial copper mineralization were accomplished by high-salinity fluids concentrated during final crystallization of the monzonitic parent magma; genetic continuity between magmatic and hydrothermal stages is indicated.Pervasive sericitic alteration of plagioclase is confined to rocks in the northern one-third of the Bingham stock; a subzone of argillic alteration in the north-central part of the ore body occurs within the broader area of sericitic alteration. Fluids responsible for sericitic and argillic alteration were channeled by a broad zone of northeast-trending fractures.Hydrothermal minerals and high-salinity fluid inclusions occur within a large volume of shattered rock. Boiling of fluids during crystallization of the aplitic porphyry may account for the shattering. Sericitic (and argillic) alteration were apparently super-imposed on the earlier biotitic assemblage as the hydrothermal system cooled. Cooling and hydrolytic alteration were promoted by progressive introduction of meteoric waters. The many generations of inclusions trapped from boiling fluids in the temperature range 400 degrees to 600 degrees C suggest that the system was recharged repeatedly during the period of mineralization. Estimated fluid pressures of about 800 bars in the early stages of mineralization correspond to a lithostatic load of about 3 km; pressures were even lower (less than 200 bars) in the later stages and were probably controlled by hydrostatic conditions.

Utah↗
USGS1991

Numerical simulations of hydrothermal circulation resulting from basalt intrusions in a buried spreading center

A two-dimensional, one by two-kilometer section through the seafloor was simulated with a numerical model to investigate coupled fluid and heat flow resulting from basalt intrusions in a buried spreading center. Boundary and initial conditions and physical properties of both sediments and basalt were constrained by field surveys and drilling in the Guaymas Basin, central Gulf of California. Parametric variations in these studies included sediment and basalt permeability, anisotropy in sediment permeability, and the size of heat sources. Faults were introduced through new intrusions both before and after cooling. Background heat input caused fluid convection at velocities ≤ 3 cm a −1 through shallow sediments. Eighty to ninety percent of the heat introduced at the base of the simulations exited through the upper, horizontal surface, even when the vertical boundaries were made permeable to fluid flow. The simulated injection of a 25–50 m thick basalt intrusion at a depth of 250 m resulted in about 10 yr of pore-fluid expulsion through the sea-floor in all cases, leaving the sediments above the intrusions strongly underpressured. A longer period of fluid recharge followed, sometimes accompanied by reductions in total seafloor heat output of 10% in comparison to pre-intrusion values. Additional discharge-recharge events were dispersed chaotically through the duration of the cooling period. These cycles in heat and fluid flow resulted from the response of the simulated system to a thermodynamic shock, the sudden emplacement of a large heat source, and not from mechanical displacement of sediments and pore fluids, which was not simulated. Water/rock mass ratios calculated from numerical simulations are in good agreement with geochemical estimates from materials recovered from the Guaymas Basin, assuming a bulk basalt permeability value of at least 10 −17 m 2 /(10 −2 mD). The addition of faults through intrusions and sediments in these simulations did not facilitate continuous, rapid venting. Increased heat input at the base of the faults resulted in temporarily greater fluid discharge, but the flow could not be sustained because the modeled system could not recharge cold fluid quickly enough to remove sufficient heat through the vents.

Earth and Planetary Science Letters↗
USGS2014

Practical limitations on the use of diurnal temperature signals to quantify groundwater upwelling

Groundwater upwelling to streams creates unique habitat by influencing stream water quality and temperature; upwelling zones also serve as vectors for contamination when groundwater is degraded. Temperature time series data acquired along vertical profiles in the streambed have been applied to simple analytical models to determine rates of vertical fluid flux. These models are based on the downward propagation characteristics (amplitude attenuation and phase-lag) of the surface diurnal signal. Despite the popularity of these models, there are few published characterizations of moderate-to-strong upwelling. We attribute this limitation to the thermodynamics of upwelling, under which the downward conductive signal transport from the streambed interface occurs opposite the upward advective fluid flux. Governing equations describing the advection–diffusion of heat within the streambed predict that under upwelling conditions, signal amplitude attenuation will increase, but, counterintuitively, phase-lag will decrease. Therefore the extinction (measurable) depth of the diurnal signal is very shallow, but phase lag is also short, yielding low signal to noise ratio and poor model sensitivity. Conversely, amplitude attenuation over similar sensor spacing is strong, yielding greater potential model sensitivity. Here we present streambed thermal time series over a range of moderate to strong upwelling sites in the Quashnet River, Cape Cod, Massachusetts. The predicted inverse relationship between phase-lag and rate of upwelling was observed in the field data over a range of conditions, but the observed phase-lags were consistently shorter than predicted. Analytical solutions for fluid flux based on signal amplitude attenuation return results consistent with numerical models and physical seepage meters, but the phase-lag analytical model results are generally unreasonable. Through numerical modeling we explore reasons why phase-lag may have been over-predicted by the analytical models, and develop guiding relations of diurnal temperature signal extinction depth based on stream diurnal signal amplitude, upwelling magnitude, and streambed thermal properties that will be useful in designing future experiments.

Journal of Hydrology↗
USGS2019

Pressure core analysis of geomechanical and fluid flow properties of seals associated with gas hydrate-bearing reservoirs in the Krishna-Godavari Basin, offshore India

Physical properties of the sediment directly overlying a gas hydrate reservoir provide important controls on the effectiveness of depressurizing that reservoir to extract methane from gas hydrate as an energy resource. The permeability of overlying sediment determines if a gas hydrate reservoir’s upper contact will provide an effective seal that enables efficient reservoir depressurization. Compressibility, stiffness and strength indicate how overlying sediment will deform as the in situ stress changes during production, providing engineering data for well designs. Assessing these properties requires minimally-disturbed sediment. India’s National Gas Hydrates Program Expedition 2 (NGHP-02) provided an opportunity to study these seal sediment properties, reducing disturbance from gas exsolution and bubble growth by collecting a pressure core from the seal sediment just above the primary gas hydrate reservoir at Site NGHP-02-08 in Area C of the Krishna-Godavari Basin. The effective stress chamber (ESC) and the direct shear chamber (DSC) devices in the suite of Pressure Core Characterization Tools (PCCTs) were used to measure permeability, compressibility, stiffness and shear strength at the in situ vertical stress. Geotechnical properties of the predominantly fine-grained seal layer at in situ vertical stress are in typical clay sediment ranges, with low measured permeability (0.02 mD), high compressibility (Cc = 0.26 – 0.33) and low shear strength (404 kPa). Though pressure and temperature were maintained throughout the collection and measurement process to stabilize gas hydrate, the lack of effective stress in the pressure core storage chamber and the chamber pressurization with methane-free water caused core expansion and gas hydrate in a thin coarser-grained layer to dissolve. The PCCTs can reapply in situ stress with incremental loading steps during a consolidation test to account for sediment compaction. Gas hydrate dissolution can be limited by storing cores just above freezing temperatures, and by using solid spacers to reduce the storage chamber’s free volume.

Marine and Petroleum Geology↗
USGS2008

Comment on “Flow and tracer transport in a single fracture: A stochastic model and its relation to some field observations” by L. Moreno et al.

Moreno et al. [1988] (hereinafter referred to as MT) used a particle-tracking scheme to investigate the physics of solute movement in a variable-aperture planar fracture. The spatially heterogeneous fluid velocity was assumed to be the only mechanism of solute movement; local or pore scale dispersion and molecular diffusion were assumed to be negligible. The particle-tracking scheme used by MT consisted of routing particles from node to node in a finite difference grid. In this scheme, the direction of an individual particle is randomly selected and the probability associated with the particle movement in a given direction is proportional to the fluid flux in that direction. The same method was used by Desbarats [1990] to investigate advective transport in aquifers composed of two porous media of different hydraulic conductivities.

Water Resources Research↗
USGS1997

The physics of debris flows

Recent advances in theory and experimentation motivate a thorough reassessment of the physics of debris flows. Analyses of flows of dry, granular solids and solid-fluid mixtures provide a foundation for a comprehensive debris flow theory, and experiments provide data that reveal the strengths and limitations of theoretical models. Both debris flow materials and dry granular materials can sustain shear stresses while remaining static; both can deform in a slow, tranquil mode characterized by enduring, frictional grain contacts; and both can flow in a more rapid, agitated mode characterized by brief, inelastic grain collisions. In debris flows, however, pore fluid that is highly viscous and nearly incompressible, composed of water with suspended silt and clay, can strongly mediate intergranular friction and collisions. Grain friction, grain collisions, and viscous fluid flow may transfer significant momentum simultaneously. Both the vibrational kinetic energy of solid grains (measured by a quantity termed the granular temperature) and the pressure of the intervening pore fluid facilitate motion of grains past one another, thereby enhancing debris flow mobility. Granular temperature arises from conversion of flow translational energy to grain vibrational energy, a process that depends on shear rates, grain properties, boundary conditions, and the ambient fluid viscosity and pressure. Pore fluid pressures that exceed static equilibrium pressures result from local or global debris contraction. Like larger, natural debris flows, experimental debris flows of ???10 m3 of poorly sorted, water-saturated sediment invariably move as an unsteady surge or series of surges. Measurements at the base of experimental flows show that coarse-grained surge fronts have little or no pore fluid pressure. In contrast, finer-grained, thoroughly saturated debris behind surge fronts is nearly liquefied by high pore pressure, which persists owing to the great compressibility and moderate permeability of the debris. Realistic models of debris flows therefore require equations that simulate inertial motion of surges in which high-resistance fronts dominated by solid forces impede the motion of low-resistance tails more strongly influenced by fluid forces. Furthermore, because debris flows characteristically originate as nearly rigid sediment masses, transform at least partly to liquefied flows, and then transform again to nearly rigid deposits, acceptable models must simulate an evolution of material behavior without invoking preternatural changes in material properties. A simple model that satisfies most of these criteria uses depth-averaged equations of motion patterned after those of the Savage-Hutter theory for gravity-driven flow of dry granular masses but generalized to include the effects of viscous pore fluid with varying pressure. These equations can describe a spectrum of debris flow behaviors intermediate between those of wet rock avalanches and sediment-laden water floods. With appropriate pore pressure distributions the equations yield numerical solutions that successfully predict unsteady, nonuniform motion of experimental debris flows.

Reviews of Geophysics↗
USGS2015

Tremor-genic slow slip regions may be deeper and warmer and may slip slower than non-tremor-genic regions

Slow slip events (SSEs) are observed worldwide and often coincide with tectonic tremor. Notable examples of SSEs lacking observed tectonic tremor, however, occur beneath Kīlauea Volcano, Hawaii, the Boso Peninsula, Japan, near San Juan Bautista on the San Andreas Fault, California, and recently in Central Ecuador. These SSEs are similar to other worldwide SSEs in many ways (e.g., size or duration), but lack the concurrent tectonic tremor observed elsewhere; instead, they trigger swarms of regular earthquakes. We investigate the physical conditions that may distinguish these non-tremor-genic SSEs from those associated with tectonic tremor, including slip velocity, pressure, temperature, fluids, and fault asperities, although we cannot eliminate the possibility that tectonic tremor may be obscured in highly attenuating regions. Slip velocities of SSEs at Kīlauea Volcano (∼10 −6 m/s) and Boso Peninsula (∼10 −7 m/s) are among the fastest SSEs worldwide. Kīlauea Volcano, the Boso Peninsula, and Central Ecuador are also among the shallowest SSEs worldwide, and thus have lower confining pressures and cooler temperatures in their respective slow slip zones. Fluids also likely contribute to tremor generation, and no corresponding zone of high v p /v s has been noted at Kīlauea or Boso. We suggest that the relatively faster slip velocities at Kīlauea Volcano and the Boso Peninsula result from specific physical conditions that may also be responsible for triggering swarms of regular earthquakes adjacent to the slow slip, while different conditions produce slower SSE velocities elsewhere and trigger tectonic tremor.

Alaska, California, Hawaii↗