SEARCH · USGS Science

Results for “Journal of Geophysical Research B: Solid Earth”

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 181 records · Page 10Linked to original sources

USGS2000

Style and age of late Oligocene-early Miocene deformation in the southern Stillwater Range, west central Nevada: Paleomagnetism, geochronology, and field relations

Paleomagnetic and geochronologic data combined with geologic mapping tightly restrict the timing and character of a late Oligocene to early Miocene episode of large magnitude extension in the southern Stillwater Range and adjacent regions of west central Nevada. The southern Stillwater Range was the site of an Oligocene to early Miocene volcanic center comprising (1) 28.3 to 24.3 Ma intracaldera ash flow tuffs, lava flows, and subjacent plutons associated with three calderas, (2) 24.8 to 20.7 Ma postcaldera silicic dikes and domes, and (3) unconformably overlying 15.3 to 13.0 Ma dacite to basalt lava flows, plugs, and dikes. The caldera-related tuffs, lava flows, and plutons were tilted 60°-70° either west or east during the initial period of Cenozoic deformation that accommodated over 100% extension. Directions of remanent magnetization obtained from these extrusive and intrusive, caldera-related rocks are strongly deflected from an expected Miocene direction in senses appropriate for their tilt. A mean direction for these rocks after tilt correction, however, suggests that they were also affected by a moderate (33.4° ± 11.8°) component of counterclockwise vertical axis rotation. Paleomagnetic data indicate that the episode of large tilting occurred during emplacement of 24.8 to 20.7 Ma postcaldera dikes and domes. In detail, an apparent decrease in rotation with decreasing age of individual, isotopically dated bodies of the postcaldera group indicates that most tilting occurred between 24.4 and 24.2 Ma. The onset of tilting immediately following after the final caldera eruptions suggests that the magmatism and deformation were linked. Deformation was not driven by magma buoyancy, however, because tilting equally affected the caldera systems of different ages, including their plutonic roots. It is more likely that regional extension was focused in the southern Stillwater Range due to magmatic warming and reduction of tensile strength of the brittle crust. Faults that accommodated deformation in the southern Stillwater Range initially dipped steeply and cut deeply to expose more than 9 km of crustal section. The exposed crustal sections are probably rotated blocks above an unexposed basal detachment that lay near the early Miocene brittle-ductile transition.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

On rate-state and Coulomb failure models

We examine the predictions of Coulomb failure stress and rate-state frictional models. We study the change in failure time (clock advance) Δt due to stress step perturbations (i.e., coseismic static stress increases) added to "background" stressing at a constant rate (i.e., tectonic loading) at time t 0 . The predictability of Δt implies a predictable change in seismicity rate r(t)/r 0 , testable using earthquake catalogs, where r 0 is the constant rate resulting from tectonic stressing. Models of r(t)/r 0 , consistent with general properties of aftershock sequences, must predict an Omori law seismicity decay rate, a sequence duration that is less than a few percent of the mainshock cycle time and a return directly to the background rate. A Coulomb model requires that a fault remains locked during loading, that failure occur instantaneously, and that Δt is independent of t 0 . These characteristics imply an instantaneous infinite seismicity rate increase of zero duration. Numerical calculations of r(t)/r 0 for different state evolution laws show that aftershocks occur on faults extremely close to failure at the mainshock origin time, that these faults must be "Coulomb-like," and that the slip evolution law can be precluded. Real aftershock population characteristics also may constrain rate-state constitutive parameters; a may be lower than laboratory values, the stiffness may be high, and/or normal stress may be lower than lithostatic. We also compare Coulomb and rate-state models theoretically. Rate-state model fault behavior becomes more Coulomb-like as constitutive parameter a decreases relative to parameter b. This is because the slip initially decelerates, representing an initial healing of fault contacts. The deceleration is more pronounced for smaller a, more closely simulating a locked fault. Even when the rate-state Δt has Coulomb characteristics, its magnitude may differ by some constant dependent on b. In this case, a rate-state model behaves like a modified Coulomb failure model in which the failure stress threshold is lowered due to weakening, increasing the clock advance. The deviation from a non-Coulomb response also depends on the loading rate, elastic stiffness, initial conditions, and assumptions about how state evolves.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

The failure of earthquake failure models

In this study I show that simple heuristic models and numerical calculations suggest that an entire class of commonly invoked models of earthquake failure processes cannot explain triggering of seismicity by transient or "dynamic" stress changes, such as stress changes associated with passing seismic waves. The models of this class have the common feature that the physical property characterizing failure increases at an accelerating rate when a fault is loaded (stressed) at a constant rate. Examples include models that invoke rate state friction or subcritical crack growth, in which the properties characterizing failure are slip or crack length, respectively. Failure occurs when the rate at which these grow accelerates to values exceeding some critical threshold. These accelerating failure models do not predict the finite durations of dynamically triggered earthquake sequences (e.g., at aftershock or remote distances). Some of the failure models belonging to this class have been used to explain static stress triggering of aftershocks. This may imply that the physical processes underlying dynamic triggering differs or that currently applied models of static triggering require modification. If the former is the case, we might appeal to physical mechanisms relying on oscillatory deformations such as compaction of saturated fault gouge leading to pore pressure increase, or cyclic fatigue. However, if dynamic and static triggering mechanisms differ, one still needs to ask why static triggering models that neglect these dynamic mechanisms appear to explain many observations. If the static and dynamic triggering mechanisms are the same, perhaps assumptions about accelerating failure and/or that triggering advances the failure times of a population of inevitable earthquakes are incorrect.

Journal of Geophysical Research B: Solid Earth↗

USGS2002

Fault structure and mechanics of the Hayward Fault, California from double-difference earthquake locations

The relationship between small-magnitude seismicity and large-scale crustal faulting along the Hayward Fault, California, is investigated using a double-difference (DD) earthquake location algorithm. We used the DD method to determine high-resolution hypocenter locations of the seismicity that occurred between 1967 and 1998. The DD technique incorporates catalog travel time data and relative P and S wave arrival time measurements from waveform cross correlation to solve for the hypocentral separation between events. The relocated seismicity reveals a narrow, near-vertical fault zone at most locations. This zone follows the Hayward Fault along its northern half and then diverges from it to the east near San Leandro, forming the Mission trend. The relocated seismicity is consistent with the idea that slip from the Calaveras Fault is transferred over the Mission trend onto the northern Hayward Fault. The Mission trend is not clearly associated with any mapped active fault as it continues to the south and joins the Calaveras Fault at Calaveras Reservoir. In some locations, discrete structures adjacent to the main trace are seen, features that were previously hidden in the uncertainty of the network locations. The fine structure of the seismicity suggests that the fault surface on the northern Hayward Fault is curved or that the events occur on several substructures. Near San Leandro, where the more westerly striking trend of the Mission seismicity intersects with the surface trace of the (aseismic) southern Hayward Fault, the seismicity remains diffuse after relocation, with strong variation in focal mechanisms between adjacent events indicating a highly fractured zone of deformation. The seismicity is highly organized in space, especially on the northern Hayward Fault, where it forms horizontal, slip-parallel streaks of hypocenters of only a few tens of meters width, bounded by areas almost absent of seismic activity. During the interval from 1984 to 1998, when digital waveforms are available, we find that fewer than 6.5% of the earthquakes can be classified as repeating earthquakes, events that rupture the same fault patch more than one time. These most commonly are located in the shallow creeping part of the fault, or within the streaks at greater depth. The slow repeat rate of 2–3 times within the 15-year observation period for events with magnitudes around M = 1.5 is indicative of a low slip rate or a high stress drop. The absence of microearthquakes over large, contiguous areas of the northern Hayward Fault plane in the depth interval from ∼5 to 10 km and the concentrations of seismicity at these depths suggest that the aseismic regions are either locked or retarded and are storing strain energy for release in future large-magnitude earthquakes.

California↗

USGS2001

Strain accumulation and rotation in the Eastern California Shear Zone

Although the Eastern California Shear Zone (ECSZ) (strike ∼N25°W) does not quite coincide with a small circle drawn about the Pacific-North America pole of rotation, trilateration and GPS measurements demonstrate that the motion within the zone corresponds to right-lateral simple shear across a vertical plane (strike N33°W±5°) roughly parallel to the tangent to that local small circle (strike ∼N40°W). If the simple shear is released by slip on faults subparallel to the shear zone, the accumulated rotation is also released, leaving no secular rotation. South of the Garlock fault the principal faults (e.g., Calico-Blackwater fault) strike ∼N40°W, close enough to the strike of the vertical plane across which maximum right-lateral shear accumulates to almost wholly accommodate that accumulation of both strain and rotation by right-lateral slip. North of the Garlock fault dip slip as well as strike slip on the principal faults (strike ∼N20°W) is required to accommodate the simple shear accumulation. In both cases the accumulated rotation is released with the shear strain. The Garlock fault, which transects the ECSZ, is not offset by north-northwest striking faults nor, despite geological evidence for long-term left-lateral slip, does it appear at the present time to be accumulating left-lateral simple shear strain across the fault due to slip at depth. Rather the motion is explained by right-lateral simple shear across the orthogonal ECSZ. Left-lateral slip on the Garlock fault will release the shear strain accumulating there but would augment the accumulating rotation, resulting in a secular clockwise rotation rate ∼80 nrad yr −1 (4.6° Myr −1 ).

California↗

USGS2002

Complex earthquake rupture and local tsunamis

In contrast to far-field tsunami amplitudes that are fairly well predicted by the seismic moment of subduction zone earthquakes, there exists significant variation in the scaling of local tsunami amplitude with respect to seismic moment. From a global catalog of tsunami runup observations this variability is greatest for the most frequently occurring tsunamigenic subduction zone earthquakes in the magnitude range of 7 < M w < 8.5. Variability in local tsunami runup scaling can be ascribed to tsunami source parameters that are independent of seismic moment: variations in the water depth in the source region, the combination of higher slip and lower shear modulus at shallow depth, and rupture complexity in the form of heterogeneous slip distribution patterns. The focus of this study is on the effect that rupture complexity has on the local tsunami wave field. A wide range of slip distribution patterns are generated using a stochastic, self-affine source model that is consistent with the falloff of far-field seismic displacement spectra at high frequencies. The synthetic slip distributions generated by the stochastic source model are discretized and the vertical displacement fields from point source elastic dislocation expressions are superimposed to compute the coseismic vertical displacement field. For shallow subduction zone earthquakes it is demonstrated that self-affine irregularities of the slip distribution result in significant variations in local tsunami amplitude. The effects of rupture complexity are less pronounced for earthquakes at greater depth or along faults with steep dip angles. For a test region along the Pacific coast of central Mexico, peak nearshore tsunami amplitude is calculated for a large number ( N = 100) of synthetic slip distribution patterns, all with identical seismic moment ( M w = 8.1). Analysis of the results indicates that for earthquakes of a fixed location, geometry, and seismic moment, peak nearshore tsunami amplitude can vary by a factor of 3 or more. These results indicate that there is substantially more variation in the local tsunami wave field derived from the inherent complexity subduction zone earthquakes than predicted by a simple elastic dislocation model. Probabilistic methods that take into account variability in earthquake rupture processes are likely to yield more accurate assessments of tsunami hazards.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗

USGS2000

Synthetic aperture radar interferometry of Okmok volcano, Alaska: radar observations

ERS-1/ERS-2 synthetic aperture radar interferometry was used to study the 1997 eruption of Okmok volcano in Alaska. First, we derived an accurate digital elevation model (DEM) using a tandem ERS-1/ERS-2 image pair and the preexisting DEM. Second, by studying changes in interferometric coherence we found that the newly erupted lava lost radar coherence for 5-17 months after the eruption. This suggests changes in the surface backscattering characteristics and was probably related to cooling and compaction processes. Third, the atmospheric delay anomalies in the deformation interferograms were quantitatively assessed. Atmospheric delay anomalies in some of the interferograms were significant and consistently smaller than one to two fringes in magnitude. For this reason, repeat observations are important to confidently interpret small geophysical signals related to volcanic activities. Finally, using two-pass differential interferometry, we analyzed the preemptive inflation, coeruptive deflation, and posteruptive inflation and confirmed the observations using independent image pairs. We observed more than 140 cm of subsidence associated with the 1997 eruption. This subsidence occurred between 16 months before the eruption and 5 months after the eruption, was preceded by ∼18 cm of uplift between 1992 and 1995 centered in the same location, and was followed by ∼10 cm of uplift between September 1997 and 1998. The best fitting model suggests the magma reservoir resided at 2.7 km depth beneath the center of the caldera, which was ∼5 km from the eruptive vent. We estimated the volume of the erupted material to be 0.055 km 3 and the average thickness of the erupted lava to be ∼7.4 m. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

Dating of major normal fault systems using thermochronology: An example from the Raft River detachment, Basin and Range, western United States

Application of thermochronological techniques to major normal fault systems can resolve the timing of initiation and duration of extension, rates of motion on detachment faults, timing of ductile mylonite formation and passage of rocks through the crystal-plastic to brittle transition, and multiple events of extensional unroofing. Here we determine the above for the top-to-the-east Raft River detachment fault and shear zone by study of spatial gradients in 40 Ar/ 39 Ar and fission track cooling ages of footwall rocks and cooling histories and by comparison of cooling histories with deformation temperatures. Mica 40 Ar/ 39 Ar cooling ages indicate that extension-related cooling began at ∼25–20 Ma, and apatite fission track ages show that motion on the Raft River detachment proceeded until ∼7.4 Ma. Collective cooling curves show acceleration of cooling rates during extension, from 5–10°C/m.y. to rates in excess of 70–100°C/m.y. The apparent slip rate along the Raft River detachment, recorded in spatial gradients of apatite fission track ages, is 7 mm/yr between 13.5 and 7.4 Ma and is interpreted to record the rate of migration of a rolling hinge. Microstructural study of footwall mylonite indicates that deformation conditions were no higher than middle greenschist facies and that deformation occurred during cooling to cataclastic conditions. These data show that the shear zone and detachment fault represent a continuum produced by progressive exhumation and shearing during Miocene extension and preclude the possibility of a Mesozoic age for the ductile shear zone. Moderately rapid cooling in middle Eocene time likely records exhumation resulting from an older, oppositely rooted, extensional shear zone along the west side of the Grouse Creek, Raft River, and Albion Mountains.

Idaho, Utah↗

USGS2000

Strain accumulation across the Eastern California Shear Zone at latitude 36°30'N

The motion of a linear array of monuments extending across the Eastern California Shear Zone (ECSZ) has been measured from 1994 to 1999 with the Global Positioning System. The linear array is oriented N54°E, perpendicular to the tangent to the local small circle drawn about the Pacific-North America pole of rotation, and the observed motion across the ECSZ is approximated by differential rotation about that pole. The observations suggest uniform deformation within the ECSZ (strike N23°W) (26 nstrain yr −1 extension normal to the zone and 39 nstrain yr −1 simple right-lateral shear across it) with no significant deformation in the two blocks (the Sierra Nevada mountains and southern Nevada) on either side. The deformation may be imposed by right-lateral slip at depth on the individual major fault systems within the zone if the slip rates are: Death Valley-Furnace Creek fault 3.2±0.9 mm yr −1 , Hunter Mountain-Panamint Valley fault 3.3±1.6 mm yr −1 , and Owens Valley fault 6.9±1.6 mm yr −1 . However, this estimate of the slip rate on the Owens Valley fault is 3 times greater than the geologic estimate.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

On the modified Mercalli intensities and magnitudes of the 1811-1812 New Madrid earthquakes

We reexamine original felt reports from the 1811–1812 New Madrid earthquakes and determine revised isoseismal maps for the three principal mainshocks. In many cases we interpret lower values than those assigned by earlier studies. In some cases the revisions result from an interpretation of original felt reports with an appreciation for site response issues. Additionally, earlier studies had assigned modified Mercalli intensity (MMI) values of V-VII to a substantial number of reports that we conclude do not describe damage commensurate with intensities this high. We investigate several approaches to contouring the MMI values using both analytical and subjective methods. For the first mainshock on 02∶15 LT December 16, 1811, our preferred contouring yields M ω7.2–7.3 using the area-moment regressions of Johnston [1996]. For the 08∶00 LT on January 23, 1812, and 03∶45 LT on February 7, 1812, mainshocks, we obtain M ω7.0 and M ω7.4–7.5, respectively. Our magnitude for the February mainshock is consistent with the established geometry of the Reelfoot fault, which all evidence suggests to have been the causative structure for this event. We note that the inference of lower magnitudes for the New Madrid events implies that site response plays a significant role in controlling seismic hazard at alluvial sites in the central and eastern United States. We also note that our results suggest that thrusting may have been the dominant mechanism of faulting associated with the 1811–1812 sequence.

Arkansas, Illinois, Kentucky, Missouri, Tennessee↗

USGS2001

Seismic mapping of shallow fault zones in the San Gabriel Mountains from the Los Angeles Region Seismic Experiment, southern California

During the Los Angeles Region Seismic Experiment (LARSE), a reflection/refraction survey was conducted along a profile (line 1) extending from Seal Beach, California, northeastward to the Mojave Desert and crossing the Los Angeles and San Gabriel Valley basins and San Gabriel Mountains. In most shot gathers from the southern and central San Gabriel Mountains, clear secondary arrivals are seen that merge, or appear to merge, with first arrivals at three locations, including the location of the Vincent thrust fault, an exposed late Mesozoic/early Cenozoic megathrust. These secondary arrivals are interpretable as reflections in the shallow crust (<5 km depth) from a concave-upward interface that projects to the surface in the north near the Vincent thrust fault, is offset in its central part at the San Gabriel fault (an old branch of the San Andreas fault), and terminates in the south at 1 to 2 km depth at the southern mountain front. The velocity structure above and below this interface strongly suggests it is the Vincent thrust fault: intermediate velocities (6.2 km/s), consistent with mylonites overlying the Vincent thrust fault, are observed above it; lower velocities (5.8 km/s), consistent with the Pelona Schist underlying the Vincent thrust fault, are observed below it. Problems arise, however, in attempting to match this reflector to the exposed Vincent thrust fault, which is seen in outcrops east of line 1. The Vincent thrust fault is shallower than the reflector in most places. An unmapped structure (steep fault, monocline, or thrust fault) is required between line 1 and the outcrops that either drops the Vincent thrust fault down to the depths of the reflector or repeats the Vincent thrust fault beneath line 1 in the footwall of another thrust fault. An alternative interpretation of the reflector is a deep greenstone horizon within the Pelona Schist, although this alternative is not favored by the velocity structure. Copyright 2001 by the American Geophysical Union.

California↗

USGS2000

Fracture process zone in granite

In uniaxial compression tests performed on Aue granite cores (diameter 50 mm, length 100 mm), a steel loading plate was used to induce the formation of a discrete shear fracture. A zone of distributed microcracks surrounds the tip of the propagating fracture. This process zone is imaged by locating acoustic emission events using 12 piezoceramic sensors attached to the samples. Propagation velocity of the process zone is varied by using the rate of acoustic emissions to control the applied axial force. The resulting velocities range from 2 mm/s in displacement-controlled tests to 2 μm/s in tests controlled by acoustic emission rate. Wave velocities and amplitudes are monitored during fault formation. P waves transmitted through the approaching process zone show a drop in amplitude of 26 dB, and ultrasonic velocities are reduced by 10%. The width of the process zone is ∼9 times the grain diameter inferred from acoustic data but is only 2 times the grain size from optical crack inspection. The process zone of fast propagating fractures is wider than for slow ones. The density of microcracks and acoustic emissions increases approaching the main fracture. Shear displacement scales linearly with fracture length. Fault plane solutions from acoustic events show similar orientation of nodal planes on both sides of the shear fracture. The ratio of the process zone width to the fault length in Aue granite ranges from 0.01 to 0.1 inferred from crack data and acoustic emissions, respectively. The fracture surface energy is estimated from microstructure analysis to be ∼2 J. A lower bound estimate for the energy dissipated by acoustic events is 0.1 J.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Earthquake stress drop and laboratory-inferred interseismic strength recovery

We determine the scaling relationships between earthquake stress drop and recurrence interval t r that are implied by laboratory-measured fault strength. We assume that repeating earthquakes can be simulated by stick-slip sliding using a spring and slider block model. Simulations with static/kinetic strength, time-dependent strength, and rate- and state-variable-dependent strength indicate that the relationship between loading velocity and recurrence interval can be adequately described by the power law V L ∝ t r n where n≈−1. Deviations from n=−1 arise from second order effects on strength, with n>−1 corresponding to apparent time-dependent strengthening and n<−1 corresponding to weakening. Simulations with rate and state-variable equations show that dynamic shear stress drop Δτ d scales with recurrence as d Δτ d / d ln t r ≤σ e ( b-a ), where σ e is the effective normal stress, μ=τ/σ e , and ( a-b )= d μ ss / d ln V is the steady-state slip rate dependence of strength. In addition, accounting for seismic energy radiation, we suggest that the static shear stress drop Δτ s scales as d Δτ s / d ln t r ≤σ e (1 +ζ)( b-a ), where ζ is the fractional overshoot. The variation of Δτ s with ln t r for earthquake stress drops is somewhat larger than implied by room temperature laboratory values of ζ and b-a . However, the uncertainty associated with the seismic data is large and the discrepancy between the seismic observations and the rate of strengthening predicted by room temperature experiments is less than an order of magnitude.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

Acoustic properties of a crack containing magmatic or hydrothermal fluids

We estimate the acoustic properties of a crack containing magmatic or hydrothermal fluids to quantify the source properties of long-period (LP) events observed in volcanic areas assuming that a crack-like structure is the source of LP events. The tails of synthetic waveforms obtained from a model of a fluid-driven crack are analyzed by the Sompi method to determine the complex frequencies of one of the modes of crack resonance over a wide range of the model parameters α/ a and ρ f /ρ s , where αis the P wave velocity of the rock matrix, a is the sound speed of the fluid, and ρ f and p s are the densities of the fluid and rock matrix, respectively. The quality factor due to radiation loss ( Q r ) for the selected mode almost monotonically increases with increasing α/ a , while the dimensionless frequency ( v ) of the mode decreases with increasing α/ a and ρ f /ρ s . These results are used to estimate Q and v for a crack containing various types of fluids (gas-gas mixtures, liquid-gas mixtures, and dusty and misty gases) for values of a , ρ f , and quality factor due to intrinsic losses ( Q i ) appropriate for these types of fluids, in which Q is given by Q −l = Q r −l + Q i −1 . For a crack containing such fluids, we obtain Q ranging from almost unity to several hundred, which consistently explains the wide variety of quality factors measured in LP events observed at various volcanoes. We underscore the importance of dusty and misty gases containing small-size particles with radii around 1 μm to explain long-lasting oscillations with Q significantly larger than 100. Our results may provide a basis for the interpretation of spatial and temporal variations in the observed complex frequencies of LP events in terms of fluid compositions beneath volcanoes.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Rapid fluid disruption: A source for self-potential anomalies on volcanoes

Self-potential (SP) anomalies observed above suspected magma reservoirs, dikes, etc., on various volcanoes (Kilauea, Hawaii; Mount Unzen, Japan; Piton de la Fournaise, Reunion Island, Miyake Jima, Japan) result from transient surface electric fields of tens of millivolts per kilometer and generally have a positive polarity. These SP anomalies are usually attributed to electrokinetic effects where properties controlling this process are poorly constrained. We propose an alternate explanation that contributions to electric fields of correct polarity should be expected from charge generation by fluid vaporization/disruption. As liquids are vaporized or removed as droplets by gas transport away from hot dike intrusions, both charge generation and local increase in electrical resistivity by removal of fluids should occur. We report laboratory observations of electric fields in hot rock samples generated by pulses of fluid (water) through the rock at atmospheric pressure. These indicate the relative amplitudes of rapid fluid disruption (RFD) potentials and electrokinetic potentials to be dramatically different and the signals are opposite in sign. Above vaporization temperatures, RFD effects of positive sign in the direction of gas flow dominate, whereas below these temperatures, effects of negative sign dominate. This suggests that the primary contribution to observed self-potential anomalies arises from gas-related charge transport processes at temperatures high enough to produce vigorous boiling and vapor transport. At lower temperatures, the primary contribution is from electrokinetic effects modulated perhaps by changing electrical resistivity and RFD effects from high-pressure but low-temperature CO2 and SO2 gas flow ripping water molecules from saturated crustal rocks. If charge generation is continuous, as could well occur above a newly emplaced dike, positive static potentials will be set up that could be sustained for many years, and the simplest method for identifying these hot, active regions would be to identify the SP anomalies they generate.

Hawai'i↗

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↗