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 37 records · Page 2Linked to original sources
USGS1999

Analysis of stress-induced oval fractures in a borehole at Deep Sea Drilling Project Site 504, eastern equatorial Pacific

Deep Sea Drilling Project (DSDP) Hole 504B is located in the eastern equatorial Pacific Ocean and extends to a total depth of 2111 m beneath the seafloor (mbsf). Several acoustic televiewer logs have been obtained in this well during successive stages of drilling, and the resulting digital images have revealed numerous oval-shaped fractures seemingly etched into the borehole wall. A theoretical examination of these stress-induced features identifies a unique and ephemeral set of stress distributions and magnitudes that are necessary for their production. Consequently, the ovals provide a basis for quantifying the magnitudes and orientations of the maximum and minimum horizontal principal stresses, S H and S h at this site. Vertical, truncated breakouts and horizontal tensile fractures define the spatial boundaries of the ovals. Explicit criteria for their occurrence are combined with estimates for various physical properties of the rock to yield a range of possible values for the horizontal principal stresses. The conspicuous oval geometry is completed by a curved fracture that joins the vertical and horizontal components. Its degree of curvature is delineated by the modified Griffith failure criterion and is directly related to the principal stress difference ( S H − S h ). Matching a series of type curves corresponding to specific values for ( S H − S h ) with the actual undistorted well bore images allows the magnitude of the stress difference to be further constrained. With a value for ( S H − S h ) of 45 ± 5 MPa the individual magnitudes of S H and S h are determined more precisely. Final estimates for the horizontal principal stresses in DSDP Hole 504B at a depth of 1200 mbsf are 141 MPa ≤ S H ≤ 149 MPa and 91 MPa ≤ S h ≤ 109 MPa. Stress magnitudes derived from this approach rely heavily upon the values of a variety of physical properties, and complementary laboratory measurements performed on relevant rock samples provide critical information. Uncertainties in estimating these properties translate into less precise predictions of principal stresses.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem

Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second‐order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well‐constrained, probabilistic model of storm and flood occurrence, has first‐order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment‐limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high‐frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.

California↗
USGS2019

Delayed dynamic triggering of disposal-induced earthquakes observed by a dense array in Northern Oklahoma

Recent increases in earthquake occurrence rates in Oklahoma have been linked to the injection of large volumes of saltwater, a byproduct of oil and gas extraction. Here we present a detailed study of remote earthquake triggering in an area of active injection‐induced seismicity in northern Oklahoma using data from the LArge‐n Seismic Survey in Oklahoma (LASSO) temporary array and nearby permanent broadband seismic stations. We estimate changes in earthquake rates and calculate the Coulomb failure stress changes on potential receiver faults due to passing teleseismic surface waves. A statistically significant increase in seismicity is observed ∼8 hr after the 16 April 2016 M w 7.8 Ecuador earthquake. The Coulomb stress changes associated with the Ecuador earthquake are on the order of ∼1 kPa. Physical mechanisms consistent with the observed dynamic stress threshold include failure driven by activation of aseismic slip or hydrological response of the fault system.

Oklahoma↗
USGS2000

Three-dimensional seismic velocity structure of the San Francisco Bay area

Seismic travel times from the northern California earthquake catalogue and from the 1991 Bay Area Seismic Imaging Experiment (BASIX) refraction survey were used to obtain a three-dimensional model of the seismic velocity structure of the San Francisco Bay area. Nonlinear tomography was used to simultaneously invert for both velocity and hypocenters. The new hypocenter inversion algorithm uses finite difference travel times and is an extension of an existing velocity tomography algorithm. Numerous inversions were performed with different parameters to test the reliability of the resulting velocity model. Most hypocenters were relocated <2 km from their catalogue locations. Large lateral velocity variations at shallow (<4 km) depth correlate with known surface geology, including low-velocity Cenozoic sedimentary basins, high-velocity Cenozoic volcanic rocks, and outcrop patterns of the major Mesozoic geologic terranes. Salinian arc rocks have higher velocities than the Franciscan melange, which in turn are faster than Great Valley Sequence forearc rocks. The thickess of low-velocity sediment is defined, including >12 km under the Sacramento River Delta, 6 km beneath Livermore Valley, 5 km beneath the Santa Clara Valley, and 4 km beneath eastern San Pablo Bay. The Great Valley Sequence east of San Francisco Bay is 4-6 km thick. A relatively high velocity body exists in the upper 10 km beneath the Sonoma volcanic field, but no evidence for a large intrusion or magma chamber exists in the crust under The Geysers or the Clear Lake volcanic center. Lateral velocity contrasts indicate that the major strike-slip faults extend subvertically beneath their surface locations through most of the crust. Strong lateral velocity contrasts of 0.3-0.6 km/s are observed across the San Andreas Fault in the middle crust and across the Hayward, Rogers Creek, Calaveras, and Greenville Faults at shallow depth. Weaker velocity contrasts (0.1-0.3 km/s) exist across the San Andreas, Hayward, and Rogers Creek Faults at all other depths. Low spatial resolution evidence in the lower crust suggests that the top of high-velocity mafic rocks gets deeper from west to east and may be offset under the major faults. The data suggest that the major strike-slip faults extend subvertically through the middle and perhaps the lower crust and juxtapose differing lithology due to accumulated strike-slip motion. The extent and physical properties of the major geologic units as constrained by the model should be used to improve studies of seismicity, strong ground motion, and regional stress.

Journal of Geophysical Research B: Solid Earth↗
USGS2023

Rock-to-metal ratios of the rare earth elements

The relative quantities of ore mined and waste rock (i.e., overburden) removed to produce the rare earth elements—their rock-to-metal ratios—were calculated for 21 individual operations or regions covering nearly all mine production in 2018. The results indicate that the rock-to-metal ratios for the total rare earth elements ranged from a low of 1.6 × 10 1 to a high of 3.6 × 10 3 , with operations in Brazil and Russia having the lowest ratios and ion-adsorption clays operations in China and Myanmar having the highest. For comparison, the global average rock-to-metal ratio for the total rare earth elements (9.8 × 10 2 ) fell between that of cobalt (8.6 × 10 2 ) and tungsten (1.1 × 10 3 ). Driven by their relative abundance in the ore and unit prices that were used in the economic allocation of the environmental burdens, the global rock-to-metal ratio for individual rare earth elements was lowest for cerium (2.3 × 10 1 ) and lanthanum (7.7 × 10 1 ) and highest for dysprosium (1.7 × 10 4 ), terbium (3.7 × 10 4 ), and lutetium (6.4 × 10 4 ). Like the rock-to-metal ratios for the total rare earth elements, rock-to-metal ratios for individual rare earth elements varied by roughly two orders of magnitude among the various operations examined. An alternative perspective of only accounting for the overburden that is physically removed in ion-adsorption clays in-situ operations yielded global rock-to-metal ratios that were an order of magnitude lower or less for many of the rare earth elements.

Journal of Cleaner Production↗
USGS1999

Formation of natural gas hydrates in marine sediments 1. Conceptual model of gas hydrate growth conditioned by host sediment properties

The stability of submarine gas hydrates is largely dictated by pressure and temperature, gas composition, and pore water salinity. However, the physical properties and surface chemistry of deep marine sediments may also affect the thermodynamic state, growth kinetics, spatial distributions, and growth forms of clathrates. Our conceptual model presumes that gas hydrate behaves in a way analogous to ice in a freezing soil. Hydrate growth is inhibited within fine-grained sediments by a combination of reduced pore water activity in the vicinity of hydrophilic mineral surfaces, and the excess internal energy of small crystals confined in pores. The excess energy can be thought of as a "capillary pressure" in the hydrate crystal, related to the pore size distribution and the state of stress in the sediment framework. The base of gas hydrate stability in a sequence of fine sediments is predicted by our model to occur at a lower temperature (nearer to the seabed) than would be calculated from bulk thermodynamic equilibrium. Capillary effects or a build up of salt in the system can expand the phase boundary between hydrate and free gas into a divariant field extending over a finite depth range dictated by total methane content and pore-size distribution. Hysteresis between the temperatures of crystallization and dissociation of the clathrate is also predicted. Growth forms commonly observed in hydrate samples recovered from marine sediments (nodules, and lenses in muds; cements in sands) can largely be explained by capillary effects, but kinetics of nucleation and growth are also important. The formation of concentrated gas hydrates in a partially closed system with respect to material transport, or where gas can flush through the system, may lead to water depletion in the host sediment. This "freeze-drying" may be detectable through physical changes to the sediment (low water content and overconsolidation) and/or chemical anomalies in the pore waters and metastable presence of free gas within the normal zone of hydrate stability.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Timescales and processes of methane hydrate formation and breakdown, with application to geologic systems

Gas hydrate is an ice-like form of water and low molecular weight gas stable at temperatures of roughly -10ºC to 25ºC and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one-sixth of Earth’s methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean-atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean-atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth’s near-surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some fields studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

Journal of Geophysical Research-Solid Earth↗
USGS2024

Subsurface redox interactions regulate ebullitive methane flux in heterogeneous Mississippi River deltaic wetland

As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive-transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient-driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

Journal of Advances in Modeling Earth Systems↗
USGS2021

Temporal seismic velocity variations: Recovery following from the 2019 Mw 7.1 Ridgecrest, California earthquake

We investigated seismic velocity changes ( dv / v ) associated with the 2019 Ridgecrest earthquake sequence with high‐frequency autocorrelations of ambient seismic noise data. Daily autocorrelation functions were computed for the entirety of 2019 and the first quarter of 2020 for broadband stations within the region, including the temporary broadband stations installed during the aftershock deployment. Travel time shifts in the daily autocorrelation functions, relative to the mean autocorrelation waveform, were computed to produce dv / v time series, which are sensitive to the evolving material properties of the shallow crust surrounding the Ridgecrest fault zone (RFZ). A short‐term velocity drop follows the M w 7.1 earthquake at stations in the vicinity of the rupture surface, while those greater than 50 km away showed no such drop. The maximum, absolute changes in seismic velocity are proportional to the logarithm of distance from the fault rupture and to the peak dynamic strain experienced during the earthquake. Near the areas of the highest coseismic slip within the RFZ, seismic velocities recovered over 3 months. However, in the vicinity of the nearby Garlock fault, where triggered slip manifested, and north of the RFZ, seismic velocities recovered within a month. We interpret the seismic velocity changes and their recovery to be largely due to changes in the physical properties of the shallow crust, such as fault zone damage recovery caused by the earthquake rupture process and in response to the large dynamic stresses of passing seismic waves from the mainshock.

California↗
USGS2021

High-resolution soil-moisture maps over landslide regions in northern California grassland derived From SAR backscattering coefficients

Slow-moving landslides are destabilized by accumulated precipitation and consequent soil moisture. Yet, the continuous high-resolution soil-moisture measurements needed to aid the understanding of landslide processes are generally absent in steep terrain. Here, we produce soil-moisture time-series maps for a seasonally active grassland landslide in the northern California coast ranges, USA, using backscattering coefficients from NASA's uninhabited aerial vehicle synthetic aperture radar at 6-m resolution. A physically based radar scattering model is used to retrieve the near-surface (5-cm depth) soil moisture for the landslide. Both forward modeling (backscattering estimation) and the retrieval (soil-moisture validation) show good agreement. The root-mean-square errors (RMSE) for vertical transmit vertical receive (VV) and horizontal transmit horizontal receive (HH) polarizations in forward model comparison are 1.93 dB and 1.88 dB, respectively. The soil-moisture retrieval shows unbiased RMSE of 0.054 m 3 /m 3 . Our successful retrieval benefits from the surface and double-bounce scattering, which is common in grasslands. The retrieved maps show saturated wetness conditions within the active landslide boundaries. We also performed sensitivity tests for incidence angle and found that the retrieval is weakly dependent on the angle, especially while using copolarized HH and VV together. Using the two copolarized inputs, the retrieval is also not sensitive to the change of orientation angles of grass cylinders. The physical model inversion presented here can be generally applied for soil-moisture retrieval in areas with the same vegetation cover types in California.

California↗
USGS1997

The physical basis of glacier volume-area scaling

Ice volumes are known for only a few of the roughly 160,000 glaciers worldwide but are important components of many climate and sea level studies which require water flux estimates. A scaling analysis of the mass and momentum conservation equations shows that glacier volumes can be related by a power law to more easily observed glacier surface areas. The relationship requires four closure choices for the scaling behavior of glacier widths, slopes, side drag and mass balance. Reasonable closures predict a volume-area scaling exponent which is consistent with observations, giving a physical and practical basis for estimating ice volumes. Glacier volume is insensitive to perturbations in the mass balance scaling, but changes in average accumulation area ratios reflect significant changes in the scaling of both mass balance and ice volume.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Introduction and digest to the Special Issue on Chemical Effects of Water on the Deformation and Strengths of Rocks

The important role of pore pressure in promoting such brittle processes as cataclasis, hydraulic fracturing, large‐scale faulting, and earthquakes within the crust is widely accepted in geology and geophysics [, 1957; , 1959; , 1958; ., 1963; , 1968; ., 1968; ., 1976; , 1973, 1980; , 1981]. Provided that fluid pressure is fully communicated with rock pore space, the effective normal stresses that control crack growth, macroscopic fracture, and friction are reduced by the magnitude of the fluid pressui'e. Beyond this physical effect of pore fluids, there are chemical effects of water on the strength of rocks that are also important in governing differential stresseg and flow in the continental crust. Some of these chemical effects of water on rock deformation have long been recognized.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Landsat yesterday and today: An American vision and an old challenge

Since the late 1960s, the United States government has invested more than $1 billion in designing, launching, and operating the Landsat (land satellite) series of Earth-observing satellites. Global change researchers, geologists, and environmental scientists have used images gathered by the satellites for purposes ranging from human health research, energy exploration, and pollution detection to agricultural assessments, urban growth monitoring, and earthquake lineament studies. The earliest data were captured on a digital medium called wide-band video tape (WBVT). However, two decades of unsound media storage conditions and a poorly maintained processing system have left the physically deteriorating WBVTs with no mechanism for interpretation. A national treasure was in jeopardy. With seed money from the National Aeronautics and Space Administration (NASA), the U.S. Geological Survey (USGS) began a project to rescue the data. More than 21,000 tapes from the 1970s have been transcribed to stable, archival media, preserving the data for future studies in Earth System Science.

Journal of Map & Geography Libraries↗
USGS2012

Advances in the simulation and automated measurement of well-sorted granular material: 1. Simulation

1. In this, the first of a pair of papers which address the simulation and automated measurement of well-sorted natural granular material, a method is presented for simulation of two-phase (solid, void) assemblages of discrete non-cohesive particles. The purpose is to have a flexible, yet computationally and theoretically simple, suite of tools with well constrained and well known statistical properties, in order to simulate realistic granular material as a discrete element model with realistic size and shape distributions, for a variety of purposes. The stochastic modeling framework is based on three-dimensional tessellations with variable degrees of order in particle-packing arrangement. Examples of sediments with a variety of particle size distributions and spatial variability in grain size are presented. The relationship between particle shape and porosity conforms to published data. The immediate application is testing new algorithms for automated measurements of particle properties (mean and standard deviation of particle sizes, and apparent porosity) from images of natural sediment, as detailed in the second of this pair of papers. The model could also prove useful for simulating specific depositional structures found in natural sediments, the result of physical alterations to packing and grain fabric, using discrete particle flow models. While the principal focus here is on naturally occurring sediment and sedimentary rock, the methods presented might also be useful for simulations of similar granular or cellular material encountered in engineering, industrial and life sciences.

Journal of Geophysical Research F: Earth Surface↗
USGS2023

Tree-ring derived avalanche frequency and climate associations in a high-latitude, maritime climate

Snow avalanches are a natural hazard in mountainous areas worldwide with severe impacts that include fatalities, damage to infrastructure, disruption to commerce, and landscape disturbance. Understanding long-term avalanche frequency patterns, and associated climate and weather influences, improves our understanding of how climate change may affect avalanche activity. We used dendrochronological techniques to evaluate the historical frequency of large magnitude avalanches (LMAs) in the high-latitude climate of southeast Alaska, United States. We collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between 1720 and 2018, which allowed us to reconstruct 82 years with LMA activity across three sub-regions. By combining this tree-ring-derived avalanche data set with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found that February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of LMA activity in the study area. Specifically, LMA activity occurred during winters with substantial February and March precipitation and neutral or negative (cold) ONI values, while years not characterized by LMAs occur more frequently during warm winters (positive ONI values). Our examination of the climate-avalanche relationship in southeast Alaska sheds light on important climate variables and physical processes associated with LMA years. These results can be used to inform long-term infrastructure planning and avalanche mitigation operations in an urban area, such as Juneau, where critical infrastructure is subject to substantial avalanche hazard.

Alaska↗
USGS1973

The chemistry of five accessory rock-forming apatites

Chemical and physical data are given for five samples of rock-forming apatite from diverse geologic environments in Nevada and Colorado. Four of these apatites contain rare-earth assemblages in which the cerium group is well represented but the yttrium group predominates. The fifth apatite contains a highly fractionated assemblage of the lighter (cerium group) rare earths similar to the assemblage typical of alkulic rocks.

Nevada, Colorado↗
USGS2010

The relationship between noise correlation and the Green's function in the presence of degeneracy and the absence of equipartition

Recent derivations have shown that when noise in a physical system has its energy equipartitioned into the modes of the system, there is a convenient relationship between the cross correlation of time-series recorded at two points and the Green's function of the system. Here, we show that even when energy is not fully equipartitioned and modes are allowed to be degenerate, a similar (though less general) property holds for equations with wave equation structure. This property can be used to understand why certain seismic noise correlation measurements are successful despite known degeneracy and lack of equipartition on the Earth. No claim to original US government works Journal compilation ?? 2010 RAS.

Geophysical Journal International↗
USGS2012

Interseismic uplift at the Nankai subduction zone, southwest Japan, 1951–1990

Uplift as a function of time from 1951 through 1990 has been deduced from annual mean sea level measured at 15 tide gages along the Nankai subduction zone, southwest Japan. The recurrence interval for rupture of the Nankai subduction zone is about 100 years, and the most recent rupture was in late 1946. Thus the 1951–1990 uplift record covers most of the first half of the earthquake cycle. The precision of the uplift record is better than could be obtained currently by annual Global Positioning System (GPS) measurements. The pre‐1960 uplift record shows rapid deformation that appears to terminate in 1959. After 1959 the uplift record is remarkably linear in time: Significant curvature in the uplift‐versus‐time plots is detected at only three of the 15 tide gages. The inferred uplift rates are not quantitatively consistent with the predictions of either the viscoelastic coupling or elastic half‐space models of subduction, but the agreement is sufficient to suggest that the causative physical mechanisms have been identified. The immediate postseismic response is controlled by the propagation of slip downward along the downdip extension of the coseismic rupture, and the interseismic response is caused by the accumulation of a slip deficit on the main thrust zone. Asthenosphere relaxation is not required to explain the observations.

Journal of Geophysical Research B: Solid Earth↗