Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Fluid Mechanics”

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

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

At least 19 recordsLinked to original sources
USGS2017

Turbulence, entrainment and low-order description of a transitional variable-density jet

Geophysical flows occur over a large range of scales, with Reynolds numbers and Richardson numbers varying over several orders of magnitude. For this study, jets of different densities were ejected vertically into a large ambient region, considering conditions relevant to some geophysical phenomena. Using particle image velocimetry, the velocity fields were measured for three different gases exhausting into air – specifically helium, air and argon. Measurements focused on both the jet core and the entrained ambient. Experiments considered relatively low Reynolds numbers from approximately 1500 to 10 000 with Richardson numbers near 0.001 in magnitude. These included a variety of flow responses, notably a nearly laminar jet, turbulent jets and a transitioning jet in between. Several features were studied, including the jet development, the local entrainment ratio, the turbulent Reynolds stresses and the eddy strength. Compared to a fully turbulent jet, the transitioning jet showed up to 50 % higher local entrainment and more significant turbulent fluctuations. For this condition, the eddies were non-axisymmetric and larger than the exit radius. For turbulent jets, the eddies were initially smaller and axisymmetric while growing with the shear layer. At lower turbulent Reynolds number, the turbulent stresses were more than 50 % higher than at higher turbulent Reynolds number. In either case, the low-density jet developed faster than a comparable non-buoyant jet. Quadrant analysis and proper orthogonal decomposition were also utilized for insight into the entrainment of the jet, as well as to assess the energy distribution with respect to the number of eigenmodes. Reynolds shear stresses were dominant in Q1 and Q3 and exhibited negligible contributions from the remaining two quadrants. Both analysis techniques showed that the development of stresses downstream was dependent on the Reynolds number while the spanwise location of the stresses depended on the Richardson number.

Journal of Fluid Mechanics↗
USGS2024

Pulsing in the Ahu‘ailaʻau pond-spillway system during the 2018 Kilauea Eruption: A dynamical systems perspective

During the 2018 K ī "> ī ī lauea lower East Rift Zone eruption, lava from 24 fissures inundated more than 8000 acres of land, destroying more than 700 structures over three months. Eruptive activity eventually focused at a single vent characterized by a continuously fed lava pond that was drained by a narrow spillway into a much wider, slower channelized flow. The spillway exhibited intervals of ‘pulsing’ behaviour in which the lava depth and velocity were observed to oscillate on time scales of several minutes. At the time, this was attributed to variations in vesiculation originating at depth. Here, we construct a toy fluid dynamical model of the pond–spillway system, and present an alternative hypothesis in which pulsing is generated at the surface, within this system. We posit that the appearance of pulsing is due to a supercritical Hopf bifurcation driven by an increase in the Reynolds number. Asymptotics for the limit cycle near the bifurcation point are derived with averaging methods and compare favourably with the cycle periodicity. Because oscillations in the pond were not observable directly due to the elevation of the cone rim and an obscuring volcanic plume, we model the observations using a spatially averaged Saint-Venant model of the spillway forced by the pond oscillator. The predicted spillway cycle periodicity and waveforms compare favourably with observations made during the eruption. The unusually well-documented nature of this eruption enables estimation of the viscosity of the erupting lava.

Hawaii↗
USGS2006

On the interpretation of energy and energy fluxes of nonlinear internal waves: An example from Massachusetts Bay

A self-consistent formalism to estimate baroclinic energy densities and fluxes resulting from the propagation of internal waves of arbitrary amplitude is derived using the concept of available potential energy. The method can be applied to numerical, laboratory or field data. The total energy flux is shown to be the sum of the linear energy flux ??? u??? p??? dz (primes denote baroclinic quantities), plus contributions from the non-hydrostatic pressure anomaly and the self-advection of kinetic and available potential energy. Using highly resolved observations in Massachusetts Bay, it is shown that due to the presence of nonlinear internal waves periodically propagating in the area, ??? u??? p??? dz accounts for only half of the total flux. The same data show that equipartition of available potential and kinetic energy can be violated, especially when the nonlinear waves begin to interact with the bottom. ?? 2006 Cambridge University Press.

Massachusetts Bay↗
USGS1977

Laboratory hydraulic fracturing experiments in intact and pre-fractured rock

Laboratory hydraulic fracturing experiments were conducted to investigate two factors which could influence the use of the hydrofrac technique for in-situ stress determinations; the possible dependence of the breakdown pressure upon the rate of borehole pressurization, and the influence of pre-existing cracks on the orientation of generated fractures. The experiments have shown that while the rate of borehole pressurization has a marked effect on breakdown pressures, the pressure at which hydraulic fractures initiate (and thus tensile strength) is independent of the rate of borehole pressurization when the effect of fluid penetration is negligible. Thus, the experiments indicate that use of breakdown pressures rather than fracture initiation pressures may lead to an erroneous estimate of tectonic stresses. A conceptual model is proposed to explain anomalously high breakdown pressures observed when fracturing with high viscosity fluids. In this model, initial fracture propagation is presumed to be stable due to large differences between the borehole pressure and that within the fracture. In samples which contained pre-existing fractures which were 'leaky' to water, we found it possible to generate hydraulic fractures oriented parallel to the direction of maximum compression if high viscosity drilling mud was used as the fracturing fluid. ?? 1977.

International Journal of Rock Mechanics and Mining↗
USGS1998

Effects of slip, slip rate, and shear heating on the friction of granite

The stability of fault slip is sensitive to the way in which frictional strength responds to changes in slip rate and in particular to the effective velocity dependence of steady state friction Δμ ss /Δ ln V . This quantity can vary substantially with displacement, temperature and slip rate. To investigate the physical basis for this behavior and the possible influence of shear heating, we slid initially bare granite surfaces in unconfined rotary shear to displacements of hundreds of millimeters at normal stresses, σ n of 10 and 25 MPa and at room temperature. We imposed step changes in slip rate within the range 10 −2 to 10 3.5 μm/s and also monitored frictional heating with thermistors embedded in the granite. The transient response of μ to slip rate steps was fit to a rate- and state-dependent friction law using two state variables to estimate the values of several parameters in the constitutive law. The first 20 mm of slip shows rising friction and falling Δμ ss /Δ ln V ; further slip shows roughly constant friction, Δμ ss /Δ ln V and parameter values, suggesting that a steady state condition is reached on the fault surface. At V ≤ 10 μm/s, Δμ ss /Δ ln V = −0.004 ± 0.001. At higher rates the response is sensitive to normal stress: At σ n = 25 MPa granite shows a transition to effective velocity strengthening (Δμ ss /Δ ln V = 0.008 ± 0.004) at the highest slip rates tested. At 10 MPa granite shows a less dramatic change to Δμ ss /Δ ln V ≈ 0 at the highest rates. The maximum temperature measured in the granite is ∼60°C at 25 MPa and 10 3.5 μm/s. Temperatures are in general agreement with a numerical model of heat conduction which assumes spatially homogeneous frictional heating over the sliding surface. The simplest interpretation of our measurements of Δμ ss /Δ ln V is that the granite is inherently velocity weakening (∂μ ss /∂ ln V < 0) and temperature strengthening (∂μ ss /∂ T −1 < 0) at all velocities. At high slip rates the response of μ to changes in temperature from shear heating may outweigh the response to changing velocity, such that the net effect Δμ ss /Δ ln V > 0 mimics velocity strengthening. These results have implications for the frictional behavior of faults during earthquakes. High slip rates may cause a switch to effective velocity strengthening which could limit peak coseismic slip rate and stress drop. For fluid-saturated faults, strengthening by this mechanism may be partly or fully offset by weakening due to thermal pressurization of a poorly drained pore fluid.

Journal of Geophysical Research B: Solid Earth↗
USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗
USGS2006

Constraints on the mechanism of long-term, steady subsidence at Medicine Lake volcano, northern California, from GPS, leveling, and InSAR

Leveling surveys across Medicine Lake volcano (MLV) have documented subsidence that is centered on the summit caldera and decays symmetrically on the flanks of the edifice. Possible mechanisms for this deformation include fluid withdrawal from a subsurface reservoir, cooling/crystallization of subsurface magma, loading by the volcano and dense intrusions, and crustal thinning due to tectonic extension (Dzurisin et al., 1991 [Dzurisin, D., Donnelly-Nolan, J.M., Evans, J.R., Walter, S.R., 1991. Crustal subsidence, seismicity, and structure near Medicine Lake Volcano, California. Journal of Geophysical Research 96, 16, 319-16, 333.]; Dzurisin et al., 2002 [Dzurisin, D., Poland, M.P., Bürgmann, R., 2002. Steady subsidence of Medicine Lake Volcano, Northern California, revealed by repeated leveling surveys. Journal of Geophysical Research 107, 2372, doi:10.1029/2001JB000893.]). InSAR data that approximate vertical displacements are similar to the leveling results; however, vertical deformation data alone are not sufficient to distinguish between source mechanisms. Horizontal displacements from GPS were collected in the Mt. Shasta/MLV region in 1996, 1999, 2000, 2003, and 2004. These results suggest that the region is part of the western Oregon block that is rotating about an Euler pole in eastern Oregon. With this rotation removed, most sites in the network have negligible velocities except for those near MLV caldera. There, measured horizontal velocities are less than predicted from ∼10 km deep point and dislocation sources of volume loss based on the leveling data; therefore volumetric losses simulated by these sources are probably not causing the observed subsidence at MLV. This result demonstrates that elastic models of subsurface volume change can provide misleading results where additional geophysical and geological constraints are unavailable, or if only vertical deformation is known. The deformation source must be capable of causing broad vertical deformation with comparatively smaller horizontal displacements. Thermoelastic contraction of a column of hot rock beneath the volcano cannot reproduce the observed ratio of vertical to horizontal surface displacements. Models that determine deformation due to loading by the volcano and dense intrusions can be made to fit the pattern of vertical displacements by assuming a weak upper crust beneath MLV, though the subsidence rates due to surface loading must be lower than the observed displacements. Tectonic extension is almost certainly occurring based on fault orientations and focal mechanisms, but does not appear to be a major contributor to the observed deformation. We favor a model that includes a combination of sources, including extension and loading of a hot weak crust with thermal contraction of a cooling mass of rock beneath MLV, which are processes that are probably occurring at MLV. Future microgravity surveys and the planned deployment of an array of continuous GPS stations as part of a Plate Boundary Observatory volcano cluster will help to refine this model.

California↗
USGS2016

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface↗
USGS1980

Hydraulic fracturing in granite under geothermal conditions

The experimental hydraulic fracturing of granite under geothermal conditions produces tensile fracture at rapid fluid injection rates and shear fracture at slow injection rates and elevated differential stress levels. A sudden burst of acoustic emission activity accompanies tensile fracture formation whereas the acoustic emission rate increases exponentially prior to shear fracture. Temperature does not significantly affect the failure mechanism, and the experimental results have not demonstrated the occurrence of thermal fracturing. A critical result of these experiments is that fluid injection at intermediate rates and elevated differential stress levels increases permeability by more than an order of magnitude without producing macroscopic fractures, and low-level acoustic emission activity occurs simultaneously near the borehole and propagates outward into the specimen with time. Permeability measurements conducted at atmospheric pressure both before and after these experiments show that increased permeability is produced by permanent structural changes in the rock. Although results of this study have not demonstrated the occurrence of thermal fracturing, they suggest that fluid injection at certain rates in situ may markedly increase local permeability. This could prove critical to increasing the efficiency of heat exchange for geothermal energy extraction from hot dry rock. ?? 1980.

International Journal of Rock Mechanics and Mining↗
USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗
USGS2007

Scales of equilibrium and disequilibrium during cleavage formation in chlorite and biotite-grade phyllites, SE Vermont

Detailed electron microprobe analyses of phyllosilicates in crenulated phyllites from south-eastern Vermont show that grain-scale zoning is common, and sympathetic zoning in adjacent minerals is nearly universal. We interpret this to reflect a pressure-solution mechanism for cleavage development, where precipitation from a very small fluid reservoir fractionated that fluid. Multiple analyses along single muscovite, biotite and chlorite grains (30–200 μ m in length) show zoning patterns indicating Tschermakitic substitutions in muscovite and both Tschermakitic and di/trioctahedral substitutions in biotite and chlorite. Using cross-cutting relationships and mineral chemistry it is shown that these patterns persist in cleavages produced at metamorphic conditions of chlorite-grade, chlorite-grade overprinted by biotite-grade and biotite-grade. Zoning patterns are comparable in all three settings, requiring a similar cleavage-forming mechanism independent of metamorphic grade. Moreover, the use of 40 Ar/ 39 Ar geochronology demonstrates this is true regardless of age. Furthermore, samples with chlorite-grade cleavages overprinted by biotite porphyroblasts suggest the closure temperatures for the diffusion of Al, Si, Mg and Fe ions are greater than the temperature of the biotite isograd (>∼400 °C). Parallel and smoothly fanning tie lines produced by coexisting muscovite–chlorite, and muscovite–biotite pairs on compositional diagrams demonstrate effectively instantaneous chemical equilibrium and probably indicate simultaneous crystallization. These results do not support theories suggesting cleavages form in fluid-dominated systems. If crenulation cleavages formed in systems in which the chemical potentials of all major components are fixed by an external reservoir, then the compositions of individual grains defining these cleavages would be uniform. On the contrary, the fine-scale chemical zoning observed probably reflects a grain-scale process consistent with a pressure-solution mechanism in which the aqueous activities of major components are defined by local dissolution and precipitation. Thus the role of fluids was probably limited to one of catalysing pressure-solution and fluids apparently did not drive cleavage development.

Journal of Metamorphic Geology↗
USGS2012

Seismic constraints on the nature of lower crustal reflectors beneath the extending Southern Transition Zone of the Colorado Plateau, Arizona

We determine the reflection polarity and exploit variations in P and S wave reflectivity and P wave amplitude versus offset (AVO) to constrain the origin of lower crustal reflectivity observed on new three-component seismic data recorded across the structural transition of the Colorado Plateau. The near vertical incidence reflection data were collected by Stanford University in 1989 as part of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that traversed the Arizona Transition Zone of the Colorado Plateau. The results of independent waveform modeling methods are consistent with much of the lower crustal reflectivity resulting from thin, high-impedance layers. The reflection polarity of the cleanest lower crustal events is positive, which implies that these reflections result from high-velocity contrasts, and the waveform character indicates that the reflectors are probably layers less than or approximately equal to 200 m thick. The lower crustal events are generally less reflective to incident S waves than to P waves, which agrees with the predicted behavior of high-velocity mafic layering. Analysis of the P wave AVO character of lower crustal reflections demonstrates that the events maintain a constant amplitude with offset, which is most consistent with a mafic-layering model. One exception is a high-amplitude (10 dB above background) event near the base of lower crustal reflectivity which abruptly decreases in amplitude at increasing offsets. The event has a pronounced S wave response, which along with its negative AVO trend is a possible indication of the presence of fluids in the lower crust. The Arizona Transition Zone is an active but weakly extended province, which causes us to discard models of lower crustal layering resulting from shearing because of the high degree of strain required to create such layers. Instead, we favor horizontal basaltic intrusions as the primary origin of high-impedance reflectors based on (1) The fact that most xenoliths in eruptive basalts of the Transition Zone are of mafic igneous composition, (2) indications that a pulse of magmatic activity crossed the Transition Zone in the late Tertiary period, and (3) the high regional heat flow observed in the Transition Zone. The apparent presence of fluids near the base of the reflective zone may indicate a partially molten intrusion. We present a mechanism by which magma can be trapped and be induced to intrude horizontally at rheologic contrasts in extending crust.

Arizona↗
USGS2004

Non-double-couple microearthquakes at Long Valley caldera, California, provide evidence for hydraulic fracturing

Most of 26 small (0.4≲ M ≲3.1) microearthquakes at Long Valley caldera in mid-1997, analyzed using data from a dense temporary network of 69 digital three-component seismometers, have significantly non-double-couple focal mechanisms, inconsistent with simple shear faulting. We determined their mechanisms by inverting P - and S -wave polarities and amplitude ratios using linear-programming methods, and tracing rays through a three-dimensional Earth model derived using tomography. More than 80% of the mechanisms have positive (volume increase) isotropic components and most have compensated linear-vector dipole components with outward-directed major dipoles. The simplest interpretation of these mechanisms is combined shear and extensional faulting with a volume-compensating process, such as rapid flow of water, steam, or CO 2 into opening tensile cracks. Source orientations of earthquakes in the south moat suggest extensional faulting on ESE-striking subvertical planes, an orientation consistent with planes defined by earthquake hypocenters. The focal mechanisms show that clearly defined hypocentral planes in different locations result from different source processes. One such plane in the eastern south moat is consistent with extensional faulting, while one near Casa Diablo Hot Springs reflects en echelon right-lateral shear faulting. Source orientations at Mammoth Mountain vary systematically with location, indicating that the volcano influences the local stress field. Events in a ‘spasmodic burst’ at Mammoth Mountain have practically identical mechanisms that indicate nearly pure compensated tensile failure and high fluid mobility. Five earthquakes had mechanisms involving small volume decreases, but these may not be significant. No mechanisms have volumetric moment fractions larger than that of a force dipole, but the reason for this fact is unknown.

California↗
USGS2012

Geochemical evidence for water‐rock interaction along the San Andreas and Garlock Faults of California

Mesozoic granitoid rocks adjacent to the San Andreas fault in central California have retained their radiogenic Ar for the last 70 m.y. but have, generally, the highest 18 O and H 2 O + contents and the lowest D contents of all the granitoid rocks in California. The geographical coincidence of the D, 18 O, and H 2 O + patterns with the present trace of the San Andreas fault leave little doubt that some kind of groundwater circulation system has operated in the vicinity of the fault in central California. Similar isotopic patterns exist in rocks along the Garlock fault. These water‐rock interactions probably took place at temperatures <200°C, although depth, extent, and timing are not resolved. Stable isotope compositions of rocks from many localities of the earth have provided unambiguous proof that meteoric waters have descended to depths at least as great as 8 to 10 km when a heat source such as a cooling pluton was present. Massive circulation of groundwater is common in the upper crust of tectonically active areas and may affect the frictional stress state along major faults. The geochemical evidence for extensive deep circulation of groundwater and the relatively high permeabilities for most rocks of the upper crust recently compiled by Brace (1980) argues against nearly lithostatic fluid pressures as an ambient condition of the upper crust. Groundwater circulation along the San Andreas provides an efficient mechanism for diffusing the heat flow anomaly that arises in heat transport calculations that are based on thermal conduction alone.

California↗
USGS1997

Source and path effects in the wave fields of tremor and explosions at Stromboli Volcano, Italy

The wave fields generated by Strombolian activity are investigated using data from small-aperture seismic arrays deployed on the north flank of Stromboli and data from seismic and pressure transducers set up near the summit crater. Measurements of slowness and azimuth as a function of time clearly indicate that the sources of tremor and explosions are located beneath the summit crater at depths shallower than 200 m with occasional bursts of energy originating from sources extending to a depth of 3 km. Slowness, azimuth, and particle motion measurements reveal a complex composition of body and surface waves associated with topography, structure, and source properties. Body waves originating at depths shallower than 200 m dominate the wave field at frequencies of 0.5–2.5 Hz, and surface waves generated by the surficial part of the source and by scattering sources distributed around the island dominate at frequencies above 2.5 Hz. The records of tremor and explosions are both dominated by SH motion. Far-field records from explosions start with radial motion, and near-field records from those events show dominantly horizontal motion and often start with a low-frequency (1–2 Hz) precursor characterized by elliptical particle motion, followed within a few seconds by a high-frequency radial phase (1–10 Hz) accompanying the eruption of pyroclastics. The dominant component of the near- and far-field particle motions from explosions, and the timing of air and body wave phases observed in the near field, are consistent with a gas-piston mechanism operating on a shallow (<200 m deep), vertical crack-like conduit. Models of a degassing fluid column suggest that noise emissions originating in the collective oscillations of bubbles ascending in the magma conduit may provide an adequate self-excitation mechanism for sustained tremor generation at Stromboli.

Journal of Geophysical Research B: Solid Earth↗
USGS2018

Investigating the mixing efficiencies of liquid-to-liquid chemical injection manifolds for aquatic invasive species management

Aquatic invasive species (AIS) have spread throughout the United States via major rivers and tributaries. Locks and dams positioned along affected waterways, specifically lock chambers, are being evaluated as potential management sites to prevent further expansion into new areas. Recent research has shown that infusion of chemicals (e.g., carbon dioxide) into water can block or kill several invasive organisms and could be a viable option at navigational structures such as lock chambers because chemical infusion would not interfere with vessel passage or lock operation. Chemical treatments near lock structures will require large-scale fluid-mechanic systems and significant energy. Mixing must extend to all stagnation regions within a lock structure to prevent the passage of an invasive fish. This work describes the performance of both wall- and floor-based CO 2 -infused-water to water injection manifolds targeted for lock structures in terms of mixing time, mixing homogeneity, injection efficiency, and operational power requirements. Both systems have strengths and weaknesses so selection recommendations are given for applications such as open systems and closed systems.

Journal of Fluids Engineering↗
USGS1994

Precursory swarms of long-period events at Redoubt Volcano (1989-1990), Alaska: Their origin and use as a forecasting tool

During the eruption of Redoubt Volcano from December 1989 through April 1990, the Alaska Volcano Observatory issued advance warnings of several tephra eruptions based on changes in seismic activity related to the occurrence of precursory swarms of long-period (LP) seismic events (dominant period of about 0.5 s). The initial eruption on December 14 occurred after 23 years of quiescence and was heralded by a 23-hour swarm of LP events that ended abruptly with the eruption. After a series of vent-clearing explosions over the next few days, dome growth began on December 21. Another swarm, with LP events similar to those of the first, began on the 26th and ended in a major tephra eruption on January 2. Eruptions continued over the next two weeks and then ceased until February 15, when a large eruption initiated a long phase of repetitive dome-building and dome-destroying episodes that continued into April. Warnings were issued before the major events on December 14 and January 2, but as the eruptive sequence continued after January 2, the energy of the swarms decreased and forecasting became more difficult. A significant but less intense swarm preceded the February 15 eruption, which was not forecast. This eruption destroyed the only seismograph on the volcanic edifice and stymied forecasting until March 4, when the first of three new stations was installed within 3 km of the active vent. From March 4 to the end of the sequence on April 21, there were eight eruptions, six of which were preceded by detectable swarms of LP events. Although weak, these swarms provided the basis for warnings issued before the eruptions on March 23 and April 6. The initial swarm on December 13 had the following features: (1) short duration (23 hours); (2) a rapidly accelerating rate of seismic energy release over the first 18 hours of the swarm, followed by a decline of activity during the 5 hours preceding the eruption; (3) a magnitude range from −0.4 to 1.6; (4) nearly identical LP signatures with a dominant period near 0.5 s; (5) dilatational first motions everywhere; and (6) a stationary source location at a depth of 1.4 km beneath the crater. This occurrence of long-period events suggests a model involving the interaction of magma with groundwater in which magmatic gases, steam and water drive a fixed conduit at a stationary point throughout the swarm. The initiation of that sequence of events is analogous to the failure of a pressure-relief valve connecting a lower, supercharged magma-dominated reservoir to a shallow hydrothermal system. A three-dimensional model of a vibrating fluid-filled crack recently developed by Chouet is found to be compatible with the seismic data and yields the following parameters for the LP source: crack length, 280–380 m; crack width, 140–190 m; crack thickness, 0.05–0.20 m; crack stiffness, 100–200; sound speed of fluid, 0.8–1.3 km/s; compressional-wave speed of rock, 5.1 km/s; density ratio of fluid to rock, ≈0.4; and ratio of bulk modulus of fluid to rigidity of rock, 0.03–0.07. The fluid-filled crack is excited intermittently by an impulsive pressure drop that varies in magnitude within the range of 0.4 to 40 bar. Such disturbance appears to be consistent with a triggering mechanism associated with choked flow conditions in the crack.

Alaska↗
USGS1982

Laboratory measurements of reservoir rock from the Geysers geothermal field, California

Rock samples taken from two outcrops, as well as rare cores from three well bores at the Geysers geothermal field, California, were tested at temperatures and pressures similar to those found in the geothermal field. Both intact and 30?? sawcut cylinders were deformed at confining pressures of 200-1000 bars, pore pressure of 30 bars and temperatures of 150?? and 240??C. Thin-section and X-ray analysis revealed that some borehole samples had undergone extensive alteration and recrystallization. Constant strain rate tests of 10-4 and 10-6 per sec gave a coefficient of friction of 0.68. Due to the highly fractured nature of the rocks taken from the production zone, intact samples were rarely 50% stronger than the frictional strength. This result suggests that the Geysers reservoir can support shear stresses only as large as its frictional shear strength. Velocity of p-waves (6.2 km/sec) was measured on one sample. Acoustic emission and sliding on a sawcut were related to changes in pore pressure. b-values computed from the acoustic emissions generated during fluid injection were typically about 0.55. An unusually high b-value (approximately 1.3) observed during sudden injection of water into the sample may have been related to thermal cracking. ?? 1982.

International Journal of Rock Mechanics and Mining↗