USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Physics of the Earth”

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 73 records · Page 4Linked to original sources
USGS1977

Spectral characteristics of geomagnetic field variations at low and equatorial latitudes

Geomagnetic field spectra from eight standard observations at geomagnetic latitudes below 30?? were studied to determine the field characteristics unique to the equatorial region. Emphasis was placed upon those variations having periods between 5 min and 4 hr for a selection of magnetically quiet, average, and active days in 1965. The power spectral density at the equator was about ten times that the near 30?? latitude. The initial manifestation of the equatorial electrojet as evidenced by the east-west alignment of the horizontal field or the change in vertical amplitudes occurred below about 20?? latitude. Induced current effects upon the vertical component from which the Earth conductivity might be inferred could best be obtained at times and latitudes unaffected by the electrojet current. Values of about 1.6 ?? 103 mhos/m for an effective skin depth of 500-600 km were determined. The spectral amplitudes increased linearly with geomagnetic activity index, Ap. The spectral slope had a similar behavior at all latitudes. The slope changed systematically with Ap-index and showed a diurnal variation, centered on local noon, that changed form with geomagnetic activity.

Journal of Atmospheric and Terrestrial Physics↗
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↗
USGS2009

Numerical models of caldera deformation: Effects of multiphase and multicomponent hydrothermal fluid flow

Ground surface displacement (GSD) in large calderas is often interpreted as resulting from magma intrusion at depth. Recent advances in geodetic measurements of GSD, notably interferometric synthetic aperture radar, reveal complex and multifaceted deformation patterns that often require complex source models to explain the observed GSD. Although hydrothermal fluids have been discussed as a possible deformation agent, very few quantitative studies addressing the effects of multiphase flow on crustal mechanics have been attempted. Recent increases in the power and availability of computing resources allow robust quantitative assessment of the complex time-variant thermal interplay between aqueous fluid flow and crustal deformation. We carry out numerical simulations of multiphase (liquid-gas), multicomponent (H 2O-CO2) hydrothermal fluid flow and poroelastic deformation using a range of realistic physical parameters and processes. Hydrothermal fluid injection, circulation, and gas formation can generate complex, temporally and spatially varying patterns of GSD, with deformation rates, magnitudes, and geometries (including subsidence) similar to those observed in several large calderas. The potential for both rapid and gradual deformation resulting from magma-derived fluids suggests that hydrothermal fluid circulation may help explain deformation episodes at calderas that have not culminated in magmatic eruption.

Journal of Geophysical Research B: Solid Earth↗
USGS1974

Thermal inertia mapping from satellite – Discrimination of geologic units in Oman

The Nimbus III and IV satellites provide reflectance and emittance data from the earth's surface at 8-km resolution. These data have been used to derive a physical property of geologic materials termed '"thermal inertia" which appears to have great promise for discriminating surficial units. A thermal inertia map of part of Oman was produced from the Nimbus satellite measurements. Correlation of this map with a reconnaissance geologic map showed gross agreement with the major units but also suggested the need for some modifications of the reconnaissance map. Some of the anomalies were verified by comparison with a later, more detailed map; others remain unexplained and may indicate previously undiscriminated geologic units.

Journal of Research of the U.S. Geological Survey↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 2. Experimental tests

Scaling considerations indicate that miniature experiments can be used to test models of granular avalanches in which the effects of intergranular fluid and cohesion are negligible. To test predictions of a granular avalanche model described in a companion paper, we performed bench top experiments involving avalanches of dry sand across irregular basal topography that mimicked the complexity of natural terrain. The experiments employed a novel method of laser-assisted cartography to map the three-dimensional morphology of rapidly moving avalanches, thereby providing high-resolution data for comparison with model output. Model input consisted of two material properties (angles of internal and basal Coulomb friction of the sand), which were measured in independent tests, and of initial and boundary conditions that characterized the geometry of the experimental apparatus. Experimental results demonstrate that the model accurately predicts not only the gross behavior but also many details of avalanche motion, from initiation to deposition. We attribute this accuracy to a mathematical and computational formulation that conserves mass and momentum in three-dimensional physical space and satisfies the Coulomb equation in three-dimensional stress space. Our results support the hypothesis that a Coulomb proportionality between shear and normal stresses applies in moderately rapid granular flows and that complicated constitutive postulates are unnecessary if momentum conservation is strictly enforced in continuum avalanche models. Furthermore, predictions of our Coulomb continuum model contrast with those of a Coulomb point mass model, illustrating the importance of multidimensional modeling and model testing.

Journal of Geophysical Research Earth Surface↗
USGS2005

Magmatic effects of the Cobb hot spot on the Juan de Fuca Ridge

The interaction of the Juan de Fuca Ridge with the Cobb hot spot has had a considerable influence on the magmatism of the Axial Segment of the ridge, the second-order segment that overlies the hot spot. In addition to the construction of the large volcanic edifice of Axial Seamount, the Axial Segment has shallow bathymetry and a prevalence of constructional volcanic features along its 100-km length, suggesting that hot spot-derived magmas supplement and oversupply the ridge. Lavas are generally more primitive at Axial Seamount and more evolved in the Axial Segment rift zones, suggesting that fractional crystallization is enhanced with increasing distance from the hot spot because of a reduced magma supply and more rapid cooling. Although the Cobb hot spot is not an isotopically enriched plume, it produces lavas with some distinct geochemical characteristics relative to normal mid-ocean ridge basalt, such as enrichments in alkalis and highly incompatible trace elements, that can be used as tracers to identify the presence and prevalence of the hot spot influence along the ridge. These characteristics are most prominent at Axial Seamount and decline in gradients along the Axial Segment. The physical model that can best explain the geochemical observations is a scenario in which hot spot and mid-ocean ridge basalt (MORB) magmas mix to varying degrees, with the proportions controlled by the depth to the MORB source. Modeling of two-component mixing suggests that MORB is the dominant component in most Axial Segment basalts. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Spectral element modelling of fault-plane reflections arising from fluid pressure distributions

The presence of fault-plane reflections in seismic images, besides indicating the locations of faults, offers a possible source of information on the properties of these poorly understood zones. To better understand the physical mechanism giving rise to fault-plane reflections in compacting sedimentary basins, we numerically model the full elastic wavefield via the spectral element method (SEM) for several different fault models. Using well log data from the South Eugene Island field, offshore Louisiana, we derive empirical relationships between the elastic parameters (e.g. P-wave velocity and density) and the effective-stress along both normal compaction and unloading paths. These empirical relationships guide the numerical modelling and allow the investigation of how differences in fluid pressure modify the elastic wavefield. We choose to simulate the elastic wave equation via SEM since irregular model geometries can be accommodated and slip boundary conditions at an interface, such as a fault or fracture, are implemented naturally. The method we employ for including a slip interface retains the desirable qualities of SEM in that it is explicit in time and, therefore, does not require the inversion of a large matrix. We performa complete numerical study by forward modelling seismic shot gathers over a faulted earth model using SEM followed by seismic processing of the simulated data. With this procedure, we construct post-stack time-migrated images of the kind that are routinely interpreted in the seismic exploration industry. We dip filter the seismic images to highlight the fault-plane reflections prior to making amplitude maps along the fault plane. With these amplitude maps, we compare the reflectivity from the different fault models to diagnose which physical mechanism contributes most to observed fault reflectivity. To lend physical meaning to the properties of a locally weak fault zone characterized as a slip interface, we propose an equivalent-layer model under the assumption of weak scattering. This allows us to use the empirical relationships between density, velocity and effective stress from the South Eugene Island field to relate a slip interface to an amount of excess pore-pressure in a fault zone. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS2017

Decaying lava extrusion rate at El Reventador Volcano, Ecuador measured using high-resolution satellite radar

Lava extrusion at erupting volcanoes causes rapid changes in topography and morphology on the order of tens or even hundreds of meters. Satellite radar provides a method for measuring changes in topographic height over a given time period to an accuracy of meters, either by measuring the width of radar shadow cast by steep sided features, or by measuring the difference in radar phase between two sensors separated in space. We measure height changes, and hence estimate extruded lava volume flux, at El Reventador, Ecuador, between 2011 and 2016, using data from the RADARSAT-2 and TanDEM-X satellite missions. We find that 39 new lava flows were extruded between 9 February 2012 and 24 August 2016, with a cumulative volume of 44.8M m 3 dense rock equivalent and a gradually decreasing eruption rate. The average dense rock rate of lava extrusion during this time is 0.31 ± 0.02 m 3 s −1 , which is similar to the long-term average from 1972 to 2016. Apart from a volumetrically small dyke opening event between 9 March and 10 June 2012, lava extrusion at El Reventador is not accompanied by any significant magmatic ground deformation. We use a simple physics-based model to estimate that the volume of the magma reservoir under El Reventador is greater than 3 km 3 . Our lava extrusion data can be equally well fit by models representing a closed reservoir depressurising during the eruption with no magma recharge, or an open reservoir with a time-constant magma recharge rate of up to 0.35 ± 0.01 m 3 s −1 .

El Reventador Volcano↗
USGS2012

Casual instrument corrections for short-period and broadband seismometers

Of all the filters applied to recordings of seismic waves, which include source, path, and site effects, the one we know most precisely is the instrument filter. Therefore, it behooves seismologists to accurately remove the effect of the instrument from raw seismograms. Applying instrument corrections allows analysis of the seismogram in terms of physical units (e.g., displacement or particle velocity of the Earth’s surface) instead of the output of the instrument (e.g., digital counts). The instrument correction can be considered the most fundamental processing step in seismology since it relates the raw data to an observable quantity of interest to seismologists. Complicating matters is the fact that, in practice, the term “instrument correction” refers to more than simply the seismometer. The instrument correction compensates for the complete recording system including the seismometer, telemetry, digitizer, and any anti‐alias filters. Knowledge of all these components is necessary to perform an accurate instrument correction. The subject of instrument corrections has been covered extensively in the literature (Seidl, 1980; Scherbaum, 1996). However, the prospect of applying instrument corrections still evokes angst among many seismologists—the authors of this paper included. There may be several reasons for this. For instance, the seminal paper by Seidl (1980) exists in a journal that is not currently available in electronic format and cannot be accessed online. Also, a standard method for applying instrument corrections involves the programs TRANSFER and EVALRESP in the Seismic Analysis Code (SAC) package (Goldstein et al. , 2003). The exact mathematical methods implemented in these codes are not thoroughly described in the documentation accompanying SAC.

Seismological Research Letters↗
USGS2003

Why earthquakes correlate weakly with the solid Earth tides: Effects of periodic stress on the rate and probability of earthquake occurrence

We provide an explanation why earthquake occurrence does not correlate well with the daily solid Earth tides. The explanation is derived from analysis of laboratory experiments in which faults are loaded to quasiperiodic failure by the combined action of a constant stressing rate, intended to simulate tectonic loading, and a small sinusoidal stress, analogous to the Earth tides. Event populations whose failure times correlate with the oscillating stress show two modes of response; the response mode depends on the stressing frequency. Correlation that is consistent with stress threshold failure models, e.g., Coulomb failure, results when the period of stress oscillation exceeds a characteristic time tn; the degree of correlation between failure time and the phase of the driving stress depends on the amplitude and frequency of the stress oscillation and on the stressing rate. When the period of the oscillating stress is less than tn, the correlation is not consistent with threshold failure models, and much higher stress amplitudes are required to induce detectable correlation with the oscillating stress. The physical interpretation of tn is the duration of failure nucleation. Behavior at the higher frequencies is consistent with a second-order dependence of the fault strength on sliding rate which determines the duration of nucleation and damps the response to stress change at frequencies greater than 1/tn. Simple extrapolation of these results to the Earth suggests a very weak correlation of earthquakes with the daily Earth tides, one that would require >13,000 earthquakes to detect. On the basis of our experiments and analysis, the absence of definitive daily triggering of earthquakes by the Earth tides requires that for earthquakes, tn exceeds the daily tidal period. The experiments suggest that the minimum typical duration of earthquake nucleation on the San Andreas fault system is ???1 year.

Journal of Geophysical Research B: Solid Earth↗
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↗
USGS2007

The instantaneous rate dependence in low temperature laboratory rock friction and rock deformation experiments

[1] Earthquake occurrence probabilities that account for stress transfer and time-dependent failure depend on the product of the effective normal stress and a lab-derived dimensionless coefficient a . This coefficient describes the instantaneous dependence of fault strength on deformation rate, and determines the duration of precursory slip. Although an instantaneous rate dependence is observed for fracture, friction, crack growth, and low temperature plasticity in laboratory experiments, the physical origin of this effect during earthquake faulting is obscure. We examine this rate dependence in laboratory experiments on different rock types using a normalization scheme modified from one proposed by Tullis and Weeks [1987]. We compare the instantaneous rate dependence in rock friction with rate dependence measurements from higher temperature dislocation glide experiments. The same normalization scheme is used to compare rate dependence in friction to rock fracture and to low-temperature crack growth tests. For particular weak phyllosilicate minerals, the instantaneous friction rate dependence is consistent with dislocation glide. In intact rock failure tests, for each rock type considered, the instantaneous rate dependence is the same size as for friction, suggesting a common physical origin. During subcritical crack growth in strong quartzofeldspathic and carbonate rock where glide is not possible, the instantaneous rate dependence measured during failure or creep tests at high stress has long been thought to be due to crack growth; however, direct comparison between crack growth and friction tests shows poor agreement. The crack growth rate dependence appears to be higher than the rate dependence of friction and fracture by a factor of two to three for all rock types considered.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Physical property changes in hydrate-bearing sediment due to depressurization and subsequent repressurization

Physical property measurements of sediment cores containing natural gas hydrate are typically performed on material exposed, at least briefly, to non-in situ conditions during recovery. To examine the effects of a brief excursion from the gas-hydrate stability field, as can occur when pressure cores are transferred to pressurized storage vessels, we measured physical properties on laboratory-formed sand packs containing methane hydrate and methane pore gas. After depressurizing samples to atmospheric pressure, we repressurized them into the methane-hydrate stability field and remeasured their physical properties. Thermal conductivity, shear strength, acoustic compressional and shear wave amplitudes, and speeds of the original and depressurized/repressurized samples are compared. X– ray computed tomography images track how the gas-hydrate distribution changes in the hydrate-cemented sands owing to the depressurizaton/repressurization process. Because depressurization-induced property changes can be substantial and are not easily predicted, particularly in water-saturated, hydrate-bearing sediment, maintaining pressure and temperature conditions throughout the core recovery and measurement process is critical for using laboratory measurements to estimate in situ properties.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Process-based, morphodynamic hindcast of decadal deposition patterns in San Pablo Bay, California, 1856-1887

This study investigates the possibility of hindcasting-observed decadal-scale morphologic change in San Pablo Bay, a subembayment of the San Francisco Estuary, California, USA, by means of a 3-D numerical model (Delft3D). The hindcast period, 1856-1887, is characterized by upstream hydraulic mining that resulted in a high sediment input to the estuary. The model includes wind waves, salt water and fresh water interactions, and graded sediment transport, among others. Simplified initial conditions and hydrodynamic forcing were necessary because detailed historic descriptions were lacking. Model results show significant skill. The river discharge and sediment concentration have a strong positive influence on deposition volumes. Waves decrease deposition rates and have, together with tidal movement, the greatest effect on sediment distribution within San Pablo Bay. The applied process-based (or reductionist) modeling approach is valuable once reasonable values for model parameters and hydrodynamic forcing are obtained. Sensitivity analysis reveals the dominant forcing of the system and suggests that the model planform plays a dominant role in the morphodynamic development. A detailed physical explanation of the model outcomes is difficult because of the high nonlinearity of the processes. Process formulation refinement, a more detailed description of the forcing, or further model parameter variations may lead to an enhanced model performance, albeit to a limited extent. The approach potentially provides a sound basis for prediction of future developments. Parallel use of highly schematized box models and a process-based approach as described in the present work is probably the most valuable method to assess decadal morphodynamic development. Copyright ?? 2011 by the American Geophysical Union.

California↗
USGS2018

Temporal stress changes caused by earthquakes: A review

Earthquakes can change the stress field in the Earth’s lithosphere as they relieve and redistribute stress. Earthquake-induced stress changes have been observed as temporal rotations of the principal stress axes following major earthquakes in a variety of tectonic settings. The stress changes due to the 2011 Mw9.0 Tohoku-Oki, Japan, earthquake were particularly well documented. Earthquake stress rotations can inform our understanding of earthquake physics, most notably addressing the long-standing problem of whether the Earth’s crust at plate boundaries is “strong” or “weak.” Many of the observed stress rotations, including that due to the Tohoku-Oki earthquake, indicate near-complete stress drop in the mainshock. This implies low background differential stress, on the order of earthquake stress drop, supporting the weak crust model. Earthquake stress rotations can also be used to address other important geophysical questions, such as the level of crustal stress heterogeneity and the mechanisms of postseismic stress reloading. The quantitative interpretation of stress rotations is evolving from those based on simple analytical methods to those based on more sophisticated numerical modeling that can capture the spatial-temporal complexity of the earthquake stress changes.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Postseismic relaxation and aftershocks

[1] Perfettini et al. (2005) suggested that the temporal dependence of surface displacements u ( t ) measured in the epicentral area following an earthquake is related to N ( t ), the cumulative number of aftershocks, by the equation u ( t ) = a + bt + cN ( t ) + d (1 − e − αt ), where a , b , c , d , and α are constants chosen to fit the data and t is the postearthquake time. N ( t ) appears in the expression for u ( t ) because both the aftershocks and a portion of u ( t ) are thought to be driven by the same source, postseismic fault creep at subseismogenic depths on the downdip extension of the coseismic rupture. We show that this equation with the actually observed N ( t ) fits the postseismic displacements recorded on several baselines following each of five earthquakes: 1999 M7.6 Chi-Chi (Taiwan), 1999 M7.1 Hector Mine (southern California), 2002 M7.9 Denali (central Alaska), 2003 M6.5 San Simeon (central California), and 2004 M6.0 Parkfield (central California) earthquakes. Although there are plausible physical interpretations for each of the terms in the expression for u ( t ), the large number of adjustable constants ( a , b , c , d , and α ) involved in fitting the rather simple postseismic displacements diminishes the significance of the fit. Because the observed N ( t ) is well fit by the modified Omori's law, fault creep at depth presumably exhibits the same temporal dependence. That dependence could be explained if the rheology of the fault downdip from the coseismic rupture is consistent with ordinary transient creep. Montesi (2004) demonstrated that power law creep across a shear zone at depth would also produce that temporal signal.

Journal of Geophysical Research B: Solid Earth↗