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.

SourceResearch connectionFrom yearThrough year

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

USGS2011

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

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

Journal of Geophysical Research B: Solid Earth↗

USGS2020

Generation of calc-alkaline magmas during crystallization at high oxygen fugacity: An experimental and petrologic study of tephras from Buldir Volcano, western Aleutian Arc, Alaska, USA

Despite agreement that calc-alkaline volcanism occurs at subduction zones and is responsible for the genesis of continental landmasses, there is no consensus on the source of the Fe-depleted signature hallmark to calc-alkaline volcanism. In this study, we utilize mafic tephras collected from Buldir Volcano to address the genesis of strongly calc-alkaline volcanic rocks (those with a low Tholeiitic Index; ≤0·7) in a segment of the western Aleutian Arc to determine if the eruptions are plausibly part of a liquid line of descent, if they are mixtures of crustal melts and parental magmas, or if they are mixtures of melts of the mantle and the subducting slab. We conducted a series of H 2 O-saturated phase equilibrium experiments (1175–1000°C; 100 MPa) in a rapid-quench cold-seal (MHC) apparatus on the most primitive natural lava from Buldir (9·34 wt % MgO) at oxidizing conditions near the Re–ReO 2 buffer. We confirmed that all experiments equilibrated 0·3 ± 0·23 log units above the Re–ReO 2 buffer (ΔQFM ∼ +2·8) using X-ray Absorption Near Edge Structure (XANES) spectroscopy. Chromite is the liquidus phase, followed by olivine, then plagioclase, then clinopyroxene, and finally hornblende. Once clinopyroxene saturates, spinel composition shifts to magnetite. We compared our experimental results to the major element geochemistry and petrology of six tephras (51·9–54·8 wt % SiO 2 ) from Buldir collected during the 2015 field season of the GeoPRISMS shared platform. Tephras contain olivine + plagioclase + clinopyroxene + spinel ± hornblende; plagioclase comprises most of the crystalline volume, followed by either olivine or hornblende. Spinel is ubiquitous; with Cr-rich spinel inclusions in olivine and hornblende, and magnetite in the groundmass. Variations in phenocryst assemblages and compositions between samples can be attributed to differences in pre-eruptive temperatures, where hotter samples are devoid of hornblende, and contain Fo-rich olivine and plagioclase with lower An-contents, owing to the position of the mineral-in curves at fluid-saturated conditions. Experimental glasses match the depletion in FeO T observed in the tephra whole rock compositions. The continuous depletion in FeO T is attributable to saturation of spinel as a liquidus phase (initially as chromite) and continuous crystallization through the experimental series (changing to magnetite at colder temperatures). In contrast to the natural samples, the experiments show enrichment in TiO 2 with decreasing MgO, suggesting that differentiation did not occur at 100 MPa on Buldir. The TiO 2 depletion in volcanic rocks from Buldir can be accounted for if hornblende crystallization occurs close to the liquidus of a parental magma; a condition that is met at higher pressures and hydrous conditions. The emerging picture for Buldir Island is that (1) oxidizing conditions are required to drive the observed depletions in FeO T via crystallization of spinel, and (2) elevated H 2 O contents and high pressures are required to saturate hornblende close to the liquidus to reproduce the entire suite of major elements. Our study provides a mechanism to generate the calc-alkaline trends observed at Buldir without requiring mixing of slab and mantle melts. We conclude that calc-alkaline volcanic rocks with extremely low Tholeiitic Indices (0·7), like those from Buldir, cannot be generated in absence of high oxygen fugacity, even at high pressure and/or elevated water pressures.

Alaska↗

USGS1998

Application of the top specified boundary layer (TSBL) approximation to initial characterization of an inland aquifer mineralization 1. Direct contact between fresh and saltwater

This paper presents a basic study in generalized terms that originates from two needs: (1) to understand the major mechanisms involved in the mineralization of groundwater of the Great Bend Prairie aquifer of Kansas by saltwater originating from a deeper Permian bedrock formation, and (2) to develop simple, robust tools that can readily be used for local assessment and management activities in the salt-affected region. A simplified basic conceptual model is adopted, incorporating two horizontal layers of porous medium which come into contact at a specific location within the model domain. The top layer is saturated with freshwater, and the bottom layer is saturated with saltwater. The paper considers various stages of approximation which can be useful for simplified simulation of the build-up of the transition zone (TZ) between the freshwater and the saltwater. The hierarchy of approximate approaches leads to the development of the top specified boundary layer (TSBL) method, which is the major tool used in this study for initial characterization of the development of the TZ. It is shown that the thickness of the TZ is mainly determined by the characteristic dispersivity. The build-up of the TZ is completed after a time period equal to the time needed to advect a fluid particle along the whole extent of the TZ. Potential applications and the effects of natural recharge and pumpage on salinity transport in the domain are discussed and evaluated in the context of demonstrating the practicality of the TSBL approach.

Journal of Contaminant Hydrology↗

USGS2001

Modeling consolidation and dewatering near the toe of the northern Barbados accretionary complex

At the toe of the northern Barbados accretionary complex, temperature and pore water chemistry data indicate that fluid flow is channeled along the décollement and other shallow thrust faults. We examine mechanisms that may prevent consolidation and maintain high permeability over large sections of the décollement. High-resolution bulk density data from five boreholes show that the décollement is well consolidated at some sites while other sites remain underconsolidated. Underconsolidated décollement behavior is associated with kilometer-scale negative-polarity seismic reflections from the décollement plane that have been interpreted to be fluid conduits. We use a coupled fluid flow/consolidation model to simulate the loading response of a 10-km-long by 680-m-thick slice of sediment as it enters the accretionary complex. The simulations capture 185 ka (5 km) of subduction, with a load function representing the estimated effective stress of the overriding accretionary prism (3.8° taper angle). Simulation results of bulk density in the décollement 3.2 km arcward of the deformation front are compared with observations. The results show that persistent high pore pressures at the arcward edge of the simulation domain can explain underconsolidated behavior. The scenario is consistent with previous modeling results showing that high pore pressures can propagate intermittently along the décollement from deeper in the complex. Simulated seaward fluxes in the décollement (1–14 cm yr −1 ) lie between previous estimates from modeling studies of steady state (<1 cm yr −1 ) and transient (>1 m yr −1 ) flow. Maximum simulated instantaneous fluid sources (2.5×10 −13 s −1 ) are comparable to previous estimates. The simulations show minor swelling of incoming sediments (fluid sources ∼−3×10 −15 s −1 ) up to 3 km before subduction that may help to explain small-scale shearing and normal faulting proximal to the protodécollement.

Atlantic Ocean, Northern Barbados Accretionary Com↗

USGS2003

Source mechanism of long-period events at Kusatsu-Shirane Volcano, Japan, inferred from waveform inversion of the effective excitation functions

We investigate the source mechanism of long-period (LP) events observed at Kusatsu-Shirane Volcano, Japan, based on waveform inversions of their effective excitation functions. The effective excitation function, which represents the apparent excitation observed at individual receivers, is estimated by applying an autoregressive filter to the LP waveform. Assuming a point source, we apply this method to seven LP events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1-3 Hz. The results of the waveform inversions show dominant volumetric change components accompanied by single force components, common to all the events analyzed, and suggesting a repeated activation of a sub-horizontal crack located 300 m beneath the summit crater lakes. Based on these results, we propose a model of the source process of LP seismicity, in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack. The rapid discharge of fluid causes the collapse of the fluid-filled crack and excites acoustic oscillations of the crack, which produce the characteristic waveforms observed in the LP events. The presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid-filled conduit linking the crack and the base of the crater lake. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

USGS2006

Gas slug ascent through changes in conduit diameter: Laboratory insights into a volcano-seismic source process in low-viscosity magmas

Seismic signals generated during the flow and degassing of low-viscosity magmas include long-period (LP) and very-long-period (VLP) events, whose sources are often attributed to dynamic fluid processes within the conduit. We present the results of laboratory experiments designed to investigate whether the passage of a gas slug through regions of changing conduit diameter could act as a suitable source mechanism. A vertical, liquid-filled glass tube featuring a concentric diameter change was used to provide canonical insights into potentially deep or shallow seismic sources. As gas slugs ascend the tube, we observe systematic pressure changes varying with slug size, liquid depth, tube diameter, and liquid viscosity. Gas slugs undergoing an abrupt flow pattern change upon entering a section of significantly increased tube diameter induce a transient pressure decrease in and above the flare and an associated pressure increase below it, which stimulates acoustic and inertial resonant oscillations. When the liquid flow is not dominantly controlled by viscosity, net vertical forces on the apparatus are also detected. The net force is a function of the magnitude of the pressure transients generated and the tube geometry, which dictates where, and hence when, the traveling pressure pulses can couple into the tube. In contrast to interpretations of related volcano-seismic data, where a single downward force is assumed to result from an upward acceleration of the center of mass in the conduit, our experiments suggest that significant downward forces can result from the rapid deceleration of relatively small volumes of downward-moving liquid. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

Detection of crystalline hematite mineralization on Mars by the Thermal Emission Spectrometer: evidence for near-surface water

The Thermal Emission Spectrometer (TES) instrument on the Mars Global Surveyor (MGS) mission has discovered a remarkable accumulation of crystalline hematite (α-Fe 2 O 3 ) that covers an area with very sharp boundaries approximately 350 by 350–750 km in size centered near 2°S latitude between 0° and 5°W longitude (Sinus Meridiani). Crystalline hematite is uniquely identified by the presence of fundamental vibrational absorption features centered near 300, 450, and >525 cm −1 and by the absence of silicate fundamentals in the 1000 cm −1 region. Spectral features resulting from atmospheric CO 2 , dust, and water ice were removed using a radiative transfer model. The spectral properties unique to Sinus Meridiani were emphasized by removing the average spectrum of the surrounding region. The depth and shape of the hematite fundamental bands show that the hematite is crystalline and relatively coarse grained (>5–10 μm). Diameters up to and greater than hundreds of micrometers are permitted within the instrumental noise and natural variability of hematite spectra. Hematite particles <5–10 μm in diameter (as either unpacked or hard-packed powders) fail to match the TES spectra. The spectrally derived areal abundance of hematite varies with particle size from ∼10% (>30 μm diameter) to 40–60% (10 μm diameter). The hematite in Sinus Meridiani is thus distinct from the fine-grained (diameter <5–10 μm), red, crystalline hematite considered, on the basis of visible, near-IR data, to be a minor spectral component in Martian bright regions like Olympus-Amazonis. Sinus Meridiani hematite is closely associated with a smooth, layered, friable surface that is interpreted to be sedimentary in origin. This material may be the uppermost surface in the region, indicating that it might be a late stage sedimentary unit or a layered portion of the heavily cratered plains units. We consider five possible mechanisms for the formation of coarse-grained, crystalline hematite. These processes fall into two classes depending on whether they require a significant amount of near-surface water: the first is chemical precipitation that includes origin by (1) precipitation from standing, oxygenated, Fe-rich water (oxide iron formations), (2) precipitation from Fe-rich hydrothermal fluids, (3) low-temperature dissolution and precipitation through mobile ground water leaching, and (4) formation of surface coatings, and the second is thermal oxidation of magnetite-rich lavas. Weathering and alteration processes, which produce nanophase and red hematite, are not consistent with the coarse, crystalline hematite observed in Sinus Meridiani. We prefer chemical precipitation models and favor precipitation from Fe-rich water on the basis of the probable association with sedimentary materials, large geographic size, distance from a regional heat source, and lack of evidence for extensive groundwater processes elsewhere on Mars. The TES results thus provide mineralogic evidence for probable large-scale water interactions. The Sinus Meridiani region may be an ideal candidate for future landed missions searching for biotic and prebiotic environments, and the physical characteristics of this site satisfy all of the engineering requirements for the missions currently planned.

Journal of Geophysical Research E: Planets↗

USGS2012

Compositions, growth mechanisms, and temporal relations of hydrothermal sulfide‐sulfate‐silica chimneys at the northern Cleft segment, Juan de Fuca Ridge

Three active hydrothermal vents forming sulfide mounds and chimneys (Monolith, Fountain, and Pipe Organ) and more widely distributed inactive chimneys are spatially related to a system of discontinuous fissures and young sheet flow lavas at the northern Cleft segment, Juan de Fuca Ridge. The formation of zoned tubular Curich chimneys (type I) on the Monolith sulfide mound is related to focused flow of high‐temperature (to 328°C) fluid. Bulbous chimneys (type II or “beehives”) at the Monolith and Fountain vents are products of diffuse high‐temperature (to 315°C) discharge. A broader zone of vigorous mixing between the hydrothermal fluid and seawater results in quench crystallization of anhydrite‐rich shells. Columnar Zn‐sulfide‐rich chimneys with narrow channelways (type III) are constructed where focused and relatively low‐temperature (261°C) fluid vents directly from the basalt substrate. The bulk chemistry (low Cu; high Pb, Ag, and SiO 2 contents), mineralogy (pyrite‐marcasite‐wurtzite‐amorphous silica‐anglesite), colloform and filamentous textures, and oxygen isotope characteristics of inactive (type IV) chimneys indicate a low‐temperature (<250°C) origin involving diffuse and sluggish flow patterns and conductive cooling. Seafloor observations and 210 Pb data indicate that (1) type IV chimneys are products of an earlier period of hydrothermal activity that ended no more than 60 years ago but prior to the sheet flow eruption, (2) the high‐temperature Monolith and Fountain vents are manifestations of the same heating event (shallow emplacement of magma) that led to the sheet flow eruption and recent megaplumes, and (3) the Pipe Organ Vent is in a very youthful stage of development and chimney deposition postdates the sheet flow eruption.

California, Oregon, Washington↗

USGS2018

The GFDL global atmosphere and land model AM4.0/LM4.0: 1. Simulation characteristics with prescribed SSTs

In this two‐part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100 km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part 1, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode—with prescribed sea surface temperatures (SSTs) and sea‐ice distribution—is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top‐of‐atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part 2, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.

Journal of Advances in Modeling Earth Systems↗

USGS1990

Plane-strain shear dislocations moving steadily in linear elastic diffusive solids

This paper derives the stress and pore pressure fields induced by a plane-strain shear (gliding edge) dislocation moving steadily at a constant speed V in a linear elastic, fluid-infiltrated (Biot) solid. Solutions are obtained for the limiting cases in which the plane containing the moving dislocation (y = 0) is permeable and impermeable to the diffusing species. Although the solutions for the permeable and impermeable planes are required to agree with each other and with the ordinary elastic solution in the limits of V = 0 (corresponding to drained response) and V = ∞ (corresponding to undrained response), the stress and pore pressure fields differ considerably for finite nonzero velocities. For the dislocation on the impermeable plane, the pore pressure is discontinuous on y = 0 and attains values which are equal in magnitude and opposite in sign as y = 0 is approached from above and below. The solution reveals the surprising result that the pore pressure on the impermeable plane is zero everywhere behind the moving dislocation (x < 0). For the dislocation on the permeable plane, the pore pressure is zero on y = 0 and attains its maximum at about (2c/V, 2c/V) where c is the diffusivity, and the origin of the coordinate system coincides with the dislocation. For the impermeable plane, the largest pore pressure change occurs at the origin.

Journal of Applied Mechanics↗

USGS2010

Supraslab earthquake clusters above the subduction plate boundary offshore Sanriku, northeastern Japan: Seismogenesis in a graveyard of detached seamounts?

Thousands of offshore repeating earthquakes with low‐angle thrust focal mechanisms occur along the subduction plate boundary of NE Japan. Double‐difference relocation methods using P ‐ and S ‐wave arrivals reveal clusters of events above these repeating events. To assure good depth control we restrict our study to events that are close to seismic stations. These “supraslab” earthquake clusters are regional features at depths of 25 to 50 km, and most of these clusters are below the depth of the forearc Moho, which we determined from converted waves. Seismicity over this depth range does not occur under the inland area of NE Japan except just below the vicinity of the arc volcanoes. Re‐entrants in the inner trench slope indicate that repeated collisions of seamounts have occurred in the past. Our preliminary interpretation of supraslab clusters is that they represent seismicity in seamounts detached from the Pacific plate during slab descent, driven by the resistance of seamounts to subduction. Detachment during slab descent probably occurs on the sedimented and hydrothermally altered seafloor on which seamounts were originally constructed since these are known as zones of weakness during active island growth. High fluid pressure produced during dehydration of clay minerals and other low‐temperature hydrous minerals could enable detachment at depths. Seamount crust is thus accreted to forearcs, possibly leading to a long‐term component of near‐coastal uplift. Supraslab earthquake clusters may be our most direct evidence of the fates of seamounts and suggest that tectonic underplating is actively occurring in this subduction system.

Northeastern Japan↗

USGS2003

Location of long-period events below Kilauea Volcano using seismic amplitudes and accurate relative relocation

We present methods for improving the location of long-period (LP) events, deep and shallow, recorded below Kilauea Volcano by the permanent seismic network. LP events might be of particular interest to understanding eruptive processes as their source mechanism is assumed to directly involve fluid transport. However, it is usually difficult or impossible to locate their source using traditional arrival time methods because of emergent wave arrivals. At Kilauea, similar LP waveform signatures suggest the existence of LP multiplets. The waveform similarity suggests spatially close sources, while catalog solutions using arrival time estimates are widely scattered beneath Kilauea's summit caldera. In order to improve estimates of absolute LP location, we use the distribution of seismic amplitudes corrected for station site effects. The decay of the amplitude as a function of hypocentral distance is used for inferring LP location. In a second stage, we use the similarity of the events to calculate their relative positions. The analysis of the entire LP seismicity recorded between January 1997 and December 1999 suggests that a very large part of the LP event population, both deep and shallow, is generated by a small number of compact sources. Deep events are systematically composed of a weak high-frequency onset followed by a low-frequency wave train. Aligning the low-frequency wave trains does not lead to aligning the onsets indicating the two parts of the signal are dissociated. This observation favors an interpretation in terms of triggering and resonance of a magmatic conduit. Instead of defining fault planes, the precise relocation of similar LP events, based on the alignment of the high-energy low-frequency wave trains, defines limited size volumes. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2007

Methane gas hydrate effect on sediment acoustic and strength properties

To improve our understanding of the interaction of methane gas hydrate with host sediment, we studied: (1) the effects of gas hydrate and ice on acoustic velocity in different sediment types, (2) effect of different hydrate formation mechanisms on measured acoustic properties (3) dependence of shear strength on pore space contents, and (4) pore pressure effects during undrained shear. A wide range in acoustic p-wave velocities (V p ) were measured in coarse-grained sediment for different pore space occupants. V p ranged from less than 1 km/s for gas-charged sediment to 1.77–1.94 km/s for water-saturated sediment, 2.91–4.00 km/s for sediment with varying degrees of hydrate saturation, and 3.88–4.33 km/s for frozen sediment. V p measured in fine-grained sediment containing gas hydrate was substantially lower (1.97 km/s). Acoustic models based on measured V p indicate that hydrate which formed in high gas flux environments can cement coarse-grained sediment, whereas hydrate formed from methane dissolved in the pore fluid may not. The presence of gas hydrate and other solid pore-filling material, such as ice, increased the sediment shear strength. The magnitude of that increase is related to the amount of hydrate in the pore space and cementation characteristics between the hydrate and sediment grains. We have found, that for consolidation stresses associated with the upper several hundred meters of sub-bottom depth, pore pressures decreased during shear in coarse-grained sediment containing gas hydrate, whereas pore pressure in fine-grained sediment typically increased during shear. The presence of free gas in pore spaces damped pore pressure response during shear and reduced the strengthening effect of gas hydrate in sands.

Journal of Petroleum Science and Engineering↗

USGS2016

Determining CO 2 storage potential during miscible CO 2 enhanced oil recovery: Noble gas and stable isotope tracers

Rising atmospheric carbon dioxide (CO 2 ) concentrations are fueling anthropogenic climate change. Geologic sequestration of anthropogenic CO 2 in depleted oil reservoirs is one option for reducing CO 2 emissions to the atmosphere while enhancing oil recovery. In order to evaluate the feasibility of using enhanced oil recovery (EOR) sites in the United States for permanent CO 2 storage, an active multi-stage miscible CO 2 flooding project in the Permian Basin (North Ward Estes Field, near Wickett, Texas) was investigated. In addition, two major natural CO 2 reservoirs in the southeastern Paradox Basin (McElmo Dome and Doe Canyon) were also investigated as they provide CO 2 for EOR operations in the Permian Basin. Produced gas and water were collected from three different CO 2 flooding phases (with different start dates) within the North Ward Estes Field to evaluate possible CO 2 storage mechanisms and amounts of total CO 2 retention. McElmo Dome and Doe Canyon were sampled for produced gas to determine the noble gas and stable isotope signature of the original injected EOR gas and to confirm the source of this naturally-occurring CO 2 . As expected, the natural CO 2 produced from McElmo Dome and Doe Canyon is a mix of mantle and crustal sources. When comparing CO 2 injection and production rates for the CO 2 floods in the North Ward Estes Field, it appears that CO 2 retention in the reservoir decreased over the course of the three injections, retaining 39%, 49% and 61% of the injected CO 2 for the 2008, 2010, and 2013 projects, respectively, characteristic of maturing CO 2 miscible flood projects. Noble gas isotopic composition of the injected and produced gas for the flood projects suggest no active fractionation, while δ 13 C CO 2 values suggest no active CO 2 dissolution into formation water, or mineralization. CO 2 volumes capable of dissolving in residual formation fluids were also estimated along with the potential to store pure-phase supercritical CO 2 . Using a combination of dissolution trapping and residual trapping, both volumes of CO 2 currently retained in the 2008 and 2013 projects could be justified, suggesting no major leakage is occurring. These subsurface reservoirs, jointly considered, have the capacity to store up to 9 years of CO 2 emissions from an average US powerplant.

International Journal of Greenhouse Gas Control↗

USGS2010

Seismicity and fluid geochemistry at Lassen Volcanic National Park, California: Evidence for two circulation cells in the hydrothermal system

Seismic analysis and geochemical interpretations provide evidence that two separate hydrothermal cells circulate within the greater Lassen hydrothermal system. One cell originates south to SW of Lassen Peak and within the Brokeoff Volcano depression where it forms a reservoir of hot fluid (235–270°C) that boils to feed steam to the high-temperature fumarolic areas, and has a plume of degassed reservoir liquid that flows southward to emerge at Growler and Morgan Hot Springs. The second cell originates SSE to SE of Lassen Peak and flows southeastward along inferred faults of the Walker Lane belt (WLB) where it forms a reservoir of hot fluid (220–240°C) that boils beneath Devils Kitchen and Boiling Springs Lake, and has an outflow plume of degassed liquid that boils again beneath Terminal Geyser. Three distinct seismogenic zones (identified as the West , Middle , and East seismic clusters) occur at shallow depths (<6 km) in Lassen Volcanic National Park, SW to SSE of Lassen Peak and adjacent to areas of high-temperature (≤161°C) fumarolic activity (Sulphur Works, Pilot Pinnacle, Little Hot Springs Valley, and Bumpass Hell) and an area of cold, weak gas emissions (Cold Boiling Lake). The three zones are located within the inferred Rockland caldera in response to interactions between deeply circulating meteoric water and hot brittle rock that overlies residual magma associated with the Lassen Volcanic Center. Earthquake focal mechanisms and stress inversions indicate primarily N–S oriented normal faulting and E–W extension, with some oblique faulting and right lateral shear in the East cluster. The different focal mechanisms as well as spatial and temporal earthquake patterns for the East cluster indicate a greater influence by regional tectonics and inferred faults within the WLB. A fourth, deeper (5–10 km) seismogenic zone (the Devils Kitchen seismic cluster) occurs SE of the East cluster and trends NNW from Sifford Mountain toward the Devils Kitchen thermal area where fumarolic temperatures are ≤123°C. Lassen fumaroles discharge geothermal gases that indicate mixing between a N 2 -rich, arc-type component and gases derived from air-saturated meteoric recharge water. Most gases have relatively weak isotopic indicators of upper mantle or volcanic components, except for gas from Sulphur Works where δ 13 C–CO 2 , δ 34 S–H 2 S, and δ 15 N–N 2 values indicate a contribution from the mantle and a subducted sediment source in an arc volcanic setting.

California↗

USGS1998

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

Journal of Geophysical Research B: Solid Earth↗

USGS1978

Pressure gradients and boiling as mechanisms for localizing ore in porphyry systems

Fluid inclusions in ore zones of porphyry systems indicate that extensive boiling of hydrothermal fluids accompanies deposition of ore and gangue minerals. The boiling commonly accompanied a change from a lithostatic to a hydrostatic environment during evolution of an epizonal stock. Pressure gradients near the margin of the stock can determine whether ore or only a diffuse zone of mineralization is formed. A sharp drop in pressure in an epizonal environment is more likely to cause extensive boiling than a comparable change in a deeper environment, as the slope of the boiling curve steepens with an increase 'in pressure. The drop in pressure causes the hydrothermal fluids to boil and creates a crackle (stockwork) breccia, which hosts the veinlets of gangue quartz and ore minerals. The boiling selectively partitions CO 2 , H 2 S, and HCl into the vapor phase, changing the pH, composition, ionic strength, and thus the solubility product of metal complexes in the remaining liquid and causing the ore and gangue to come out of solution. Fluid inclusions trapped from boiling solutions can exhibit several forms, depending on the physical and chemical conditions of the hydrothermal fluid from which they were trapped. In one case, inclusions when heated can homogenize to either liquid or vapor at the same temperature, which is the true boiling temperature. In another case, homogenization of various inclusions can occur through a range of temperatures. The latter case results from the trapping of mixture of liquid and vapor. Variations in salinity can result from boiling of the hydrothermal fluid, or intermittent incorporation of high-salinity fluids from the magma, or trapping of fluids of varying densities at pressure-temperature conditions above the critical point of the fluid. In places, paleopressure-temperature transition zones can be recognized by fluid-inclusion homogenization temperatures and phase relationships and by the presence of anhydrite daughter minerals. Boiling of a hydrothermal fluid in the porphyry environment affects light stable isotopes. Hydrogen is preferentially fractionated into the vapor phase from water boiling below 223°C; above this temperature deuterium is selectively enriched in the vapor phase. In certain environments boiling creates a vapor-dominated system in which the condensate is swept away by meteoric waters and the H/D in the residual fluids is progressively increased through time.

Utah↗