Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Applied Geophysics”

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 109 records · Page 6Linked to original sources
USGS2009

Elastic wave speeds and moduli in polycrystalline ice Ih, si methane hydrate, and sll methane-ethane hydrate

We used ultrasonic pulse transmission to measure compressional, P, and shear, S, wave speeds in laboratory-formed polycrystalline ice Ih, si methane hydrate, and sll methane-ethane hydrate. From the wave speed's linear dependence on temperature and pressure and from the sample's calculated density, we derived expressions for bulk, shear, and compressional wave moduli and Poisson's ratio from -20 to 15??C and 22.4 to 32.8 MPa for ice Ih, -20 to 15??C and 30.5 to 97.7 MPa for si methane hydrate, and -20 to 10??C and 30.5 to 91.6 MPa for sll methane-ethane hydrate. All three materials had comparable P and S wave speeds and decreasing shear wave speeds with increasing applied pressure. Each material also showed evidence of rapid intergranular bonding, with a corresponding increase in wave speed, in response to pauses in sample deformation. There were also key differences. Resistance to uniaxial compaction, indicated by the pressure required to compact initially porous samples, was significantly lower for ice Ih than for either hydrate. The ice Ih shear modulus decreased with increasing pressure, in contrast to the increase measured in both hydrates ?? 2009.

Journal of Geophysical Research B: Solid Earth↗
USGS2022

Introduction to the special issue of the Consortium of Organizations for Strong Motion Observation Systems (COSMOS) international guidelines for applying noninvasive geophysical techniques to characterize seismic site conditions

Knowledge about local seismic site conditions provides critical information to account for site effects that are commonly observed in strong motion recordings. Certainly, other wave propagation effects can influence these observations, which are attributable to variations in material properties of the paths traveled by the waves, as well as the characteristics of the seismic source. However, local geologic conditions, particularly, when under shear-wave excitation, are known to have a strong influence on the behavior of ground shaking in the frequency range that is expected to directly affect the built environment. Thus, shear waves traveling in the shallow subsurface—defined here as tens to hundreds of meters beneath the ground surface—are the main foci for application and research in the earthquake engineering community.

Journal of Seismology↗
USGS1980

Driving- stress waveform and the determination of rock internal friction by the stress-strain curve method.

Harmonic distortion in the stress-time function applied to rock specimens affects the measurement of rock internal friction in the seismic wave periods by the stress-strain hysteresis loop method. If neglected, the harmonic distortion can cause measurements of rock internal friction to be in error by 3O% in the linear range. The stress-time function therefore must be recorded and Fourier analysed for correct interpretation of the experimental data. Such a procedure would also yield a value for internal friction at the higher harmonic frequencies.-Author

Geophysical Journal of the Royal Astronomical Soci↗
USGS2022

Quantifying large-scale surface change using SAR amplitude images: Crater morphology changes during the 2019-2020 Shishaldin Volcano eruption

Morphological processes often induce meter-scale elevation changes. When a volcano erupts, tracking such processes provides insights into the style and evolution of eruptive activity and related hazards. Compared to optical remote-sensing products, synthetic aperture radar (SAR) observes surface change during inclement weather and at night. Differential SAR interferometry estimates phase change between SAR acquisitions and is commonly applied to quantify deformation. However, large deformation or other coherence loss can limit its use. We develop a new approach applicable when repeated digital elevation models (DEMs) cannot be otherwise retrieved. Assuming an isotropic radar cross-section, we estimate meter-scale vertical morphological change directly from SAR amplitude images via an optimization method that utilizes a high-quality DEM. We verify our implementation through simulation of a collapse feature that we modulate onto topography. We simulate radar effects and recover the simulated collapse. To validate our method, we estimate elevation changes from TerraSAR-X stripmap images for the 2011–2012 eruption of Mount Cleveland. Our results reproduce those from two previous studies; one that used the same dataset, and another based on thermal satellite data. By applying this method to the 2019–2020 eruption of Shishaldin Volcano, Alaska, we generate elevation change time series from dozens of co-registered TerraSAR-X high-resolution spotlight images. Our results quantify previously unresolved cone growth in November 2019, collapses associated with explosions in December–January, and further changes in crater elevations into spring 2020. This method can be used to track meter-scale morphology changes for ongoing eruptions with low latency as SAR imagery becomes available.

Alaska↗
USGS2012

Use of weather types to disaggregate general circulation model predictions

General circulation models (GCMs) simulate climatic conditions with a grid cell resolution on the order of 100,000 km 2 . This resolution is inadequate to assess the effects of climatic change on water resources at a regional scale. A method has been developed that uses weather-type analysis as a tool to spatially disaggregate GCM predictions to make them useful for water resource studies. The method has been applied to the Delaware River basin to predict the effects of doubling atmospheric carbon dioxide on precipitation patterns in the region. An application of the technique to the Delaware River basin indicates that future climatic conditions will show minimal changes in weather-type frequency, implying that air circulation patterns will remain unchanged. Results of this study indicate that changes in regional precipitation patterns under a doubling of atmospheric carbon dioxide will be a result of within-type changes in weather characteristics.

Journal of Geophysical Research↗
USGS2013

Size distribution of Parkfield’s microearthquakes reflects changes in surface creep rate

The nucleation area of the series of M6 events in Parkfield has been shown to be characterized by low b-values throughout the seismic cycle. Since low b -values represent high differential stresses, the asperity structure seems to be always stably stressed and even unaffected by the latest main shock in 2004. However, because fault loading rates and applied shear stress vary with time, some degree of temporal variability of the b -value within stable blocks is to be expected. We discuss in this study adequate techniques and uncertainty treatment for a detailed analysis of the temporal evolution of b -values. We show that the derived signal for the Parkfield asperity correlates with changes in surface creep, suggesting a sensitive time resolution of the b-value stress meter, and confirming near-critical loading conditions within the Parkfield asperity.

Geophysical Journal International↗
USGS2005

Estimating hydrodynamic roughness in a wave-dominated environment with a high-resolution acoustic Doppler profiler

Hydrodynamic roughness is a critical parameter for characterizing bottom drag in boundary layers, and it varies both spatially and temporally due to variation in grain size, bedforms, and saltating sediment. In this paper we investigate temporal variability in hydrodynamic roughness using velocity profiles in the bottom boundary layer measured with a high-resolution acoustic Doppler profiler (PCADP). The data were collected on the ebb-tidal delta off Grays Harbor, Washington, in a mean water depth of 9 m. Significant wave height ranged from 0.5 to 3 m. Bottom roughness has rarely been determined from hydrodynamic measurements under conditions such as these, where energetic waves and medium-to-fine sand produce small bedforms. Friction velocity due to current u*c and apparent bottom roughness z0a were determined from the PCADP burst mean velocity profiles using the law of the wall. Bottom roughness kB was estimated by applying the Grant-Madsen model for wave-current interaction iteratively until the model u*c converged with values determined from the data. The resulting kB values ranged over 3 orders of magnitude (10-1 to 10-4 m) and varied inversely with wave orbital diameter. This range of kB influences predicted bottom shear stress considerably, suggesting that the use of time-varying bottom roughness could significantly improve the accuracy of sediment transport models. Bedform height was estimated from kB and is consistent with both ripple heights predicted by empirical models and bedforms in sonar images collected during the experiment. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
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↗
USGS2021

The finicky nature of earthquake shaking-triggered submarine sediment slope failures and sediment gravity flows

Since 2011, seafloor temperatures, pressures, and seismic ground motions have been measured by the seafloor cabled Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET) on the Nankai margin. These measurements, high-resolution bathymetry, and abundant contextual information make the DONET region seem ideally suited to provide constraints on seismic shaking-triggered sediment slope failures and gravity flows, particularly since numerous published studies have linked paleo- to modern earthquakes to failures and flows within the DONET. The occurrences of the local 2016 M6.0 Mie-ken and regional M7.0 Kumamoto earthquakes within and at regional distances, respectively, from the DONET data set provided an opportunity to explore this potential. We used DONET seismic recordings of the posited triggering shaking and to search for submarine slide signals and continuous temperature and pressure data to detect pulses of warm and densified water indicative of passing flows. We developed and applied a variety of analytical methods to eliminate signals generated by water column processes, while leaving slope failures and sediment gravity flow anomalies as residuals. Our explorations yielded no evidence that earthquake shaking initiated either phenomenon, which we suggest reflects the finicky nature both of the detection of and the physical processes that contribute to slope failures and flows (i.e., both require satisfying precise suites of conditions). Nonetheless, this negative result, our analyses, and the estimates of physical properties we derived for them, provide useful lessons and inputs for future studies.

Journal of Geophysical Research↗
USGS2024

Predicting characteristic length scales of barrier island segmentation in microtidal environments

Segmented barrier islands can be found in regions with small tidal ranges. In contrast to tidally dominated barriers, where inlet dynamics are thought to control island length scales, the controls on barrier island length scales in wave-dominated environments have not been quantified. These microtidal barriers typically have a curved shoreline, suggesting the influence of wave-driven alongshore sediment transport. Microtidal barriers are also typically hydrodynamically isolated from one another, as weak tidal flows limit interactions between adjoining barriers. To better understand the controls on and scales of barrier segmentation in the relative absence of tides, here we develop a theoretical framework to estimate the alongshore length scales at which a barrier will either breach or heal following a disturbance in the barrier morphology. The non-dimensional framework compares the timescales of overwash (advective) and alongshore sediment transport (diffusive) processes along barrier island chains. We then apply this framework to modern barrier islands in the microtidal Gulf of Mexico using wave hindcast data and the lengths, widths, heights, and lagoon depths measured from remotely sensed geospatial data and topobathymetric data. We find that most of these barriers are currently longer than their critical length scale, often as a result of coastal restoration efforts. Our critical length scale analysis suggests that most of the Gulf of Mexico barriers are vulnerable to segmentation despite coastal restoration efforts intended to protect fisheries and the mainland coasts.

Florida, Louisiana, Mississippi↗
USGS2011

Detecting aseismic strain transients from seismicity data

Aseismic deformation transients such as fluid flow, magma migration, and slow slip can trigger changes in seismicity rate. We present a method that can detect these seismicity rate variations and utilize these anomalies to constrain the underlying variations in stressing rate. Because ordinary aftershock sequences often obscure changes in the background seismicity caused by aseismic processes, we combine the stochastic Epidemic Type Aftershock Sequence model that describes aftershock sequences well and the physically based rate‐ and state‐dependent friction seismicity model into a single seismicity rate model that models both aftershock activity and changes in background seismicity rate. We implement this model into a data assimilation algorithm that inverts seismicity catalogs to estimate space‐time variations in stressing rate. We evaluate the method using a synthetic catalog, and then apply it to a catalog of M ≥ 1.5 events that occurred in the Salton Trough from 1990 to 2009. We validate our stressing rate estimates by comparing them to estimates from a geodetically derived slip model for a large creep event on the Obsidian Buttes fault. The results demonstrate that our approach can identify large aseismic deformation transients in a multidecade long earthquake catalog and roughly constrain the absolute magnitude of the stressing rate transients. Our method can therefore provide a way to detect aseismic transients in regions where geodetic resolution in space or time is poor.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

An effective medium inversion algorithm for gas hydrate quantification and its application to laboratory and borehole measurements of gas hydrate-bearing sediments

The presence of gas hydrate in marine sediments alters their physical properties. In some circumstances, gas hydrate may cement sediment grains together and dramatically increase the seismic P- and S-wave velocities of the composite medium. Hydrate may also form a load-bearing structure within the sediment microstructure, but with different seismic wave attenuation characteristics, changing the attenuation behaviour of the composite. Here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate-bearing sediments. The velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure. The attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures. We relate velocity and attenuation increases in hydrate-bearing sediments to their hydrate content, using an effective medium inversion algorithm based on the self-consistent approximation (SCA), differential effective medium (DEM) theory, and Biot and squirt flow mechanisms of fluid flow. The inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space. We applied our algorithm to a data set from the Mallik 2L–38 well, Mackenzie delta, Canada, and to data from laboratory measurements on gas-rich and water-saturated sand samples. Predictions using our algorithm match the borehole data and water-saturated laboratory data if the proportion of hydrate contributing to the load-bearing structure increases with hydrate saturation. The predictions match the gas-rich laboratory data if that proportion decreases with hydrate saturation. We attribute this difference to differences in hydrate formation mechanisms between the two environments.

Geophysical Journal International↗
USGS2010

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California↗
USGS2017

Accounting for orphaned aftershocks in the earthquake background rate

Aftershocks often occur within cascades of triggered seismicity in which each generation of aftershocks triggers an additional generation, and so on. The rate of earthquakes in any particular generation follows Omori's law, going approximately as 1/ t . This function decays rapidly, but is heavy-tailed, and aftershock sequences may persist for long times at a rate that is difficult to discriminate from background. It is likely that some apparently spontaneous earthquakes in the observational catalogue are orphaned aftershocks of long-past main shocks. To assess the relative proportion of orphaned aftershocks in the apparent background rate, I develop an extension of the ETAS model that explicitly includes the expected contribution of orphaned aftershocks to the apparent background rate. Applying this model to California, I find that the apparent background rate can be almost entirely attributed to orphaned aftershocks, depending on the assumed duration of an aftershock sequence. This implies an earthquake cascade with a branching ratio (the average number of directly triggered aftershocks per main shock) of nearly unity. In physical terms, this implies that very few earthquakes are completely isolated from the perturbing effects of other earthquakes within the fault system. Accounting for orphaned aftershocks in the ETAS model gives more accurate estimates of the true background rate, and more realistic expectations for long-term seismicity patterns.

Geophysical Journal International↗
USGS2016

A large refined catalog of earthquake relocations and focal mechanisms for the Island of Hawai'i and its seismotectonic implications

We present high-quality focal mechanisms based on a refined earthquake location catalog for the Island of Hawai'i, focusing on Mauna Loa and Kīlauea volcanoes. The relocation catalog is based on first-arrival times and waveform data of both compressional and shear waves for about 180,000 events on and near the Island of Hawai'i between 1986 and 2009 recorded by the seismic stations at the Hawaiian Volcano Observatory. We relocate all the earthquakes by applying ray tracing through an existing three-dimensional velocity model, similar event cluster analysis, and a differential-time relocation method. The resulting location catalog represents an expansion of previous relocation studies, covering a longer time period and consisting of more events with well-constrained absolute locations. The focal mechanisms are obtained based on the compressional-wave first-motion polarities and compressional-to-shear wave amplitude ratios by applying the HASH program to the waveform cross correlation relocated earthquakes. Overall, the good-quality (defined by the HASH parameters) focal solutions are dominated by normal faulting in our study area, especially in the active Ka'ōiki and Hīlea seismic zones. Kīlauea caldera is characterized by a mixture of approximately equal numbers of normal, strike-slip, and reverse faults, whereas its south flank has slightly fewer strike-slip events. Our relocation and focal mechanism results will be useful for mapping the seismic stress and strain fields and for understanding the seismic-volcanic-tectonic relationships within the magmatic systems.

Hawai'i↗
USGS2015

Resolving bathymetry from airborne gravity along Greenland fjords

Recent glacier mass loss in Greenland has been attributed to encroaching warming waters, but knowledge of fjord bathymetry is required to investigate this mechanism. The bathymetry in many Greenland fjords is unmapped and difficult to measure. From 2010 to 2012, National Aeronautics and Space Administration's Operation IceBridge collected a unique set of airborne gravity, magnetic, radar, and lidar data along the major outlet glaciers and fjords in Greenland. We applied a consistent technique using the IceBridge gravity data to create 90 bathymetric profiles along 54 Greenland fjords. We also used this technique to recover subice topography where warm or crevassed ice prevents the radar system from imaging the bed. Here we discuss our methodology, basic assumptions and error analysis. We present the new bathymetry data and discuss observations in six major regions of Greenland covered by IceBridge. The gravity models provide a total of 1950 line kilometers of bathymetry, 875 line kilometers of subice topography, and 12 new grounding line depths.

Journal of Geophysical Research B: Solid Earth↗
USGS2018

A model ensemble for projecting multi‐decadal coastal cliff retreat during the 21st century

Sea cliff retreat rates are expected to accelerate with rising sea levels during the 21 st century. Here we develop an approach for a multi‐model ensemble that efficiently projects time‐averaged sea cliff retreat over multi‐decadal time scales and large (>50 km) spatial scales. The ensemble consists of five simple 1‐D models adapted from the literature that relate sea cliff retreat to wave impacts, sea level rise (SLR), historical cliff behavior, and cross‐shore profile geometry. Ensemble predictions are based on Monte Carlo simulations of each individual model, which account for the uncertainty of model parameters. The consensus of the individual models also weights uncertainty, such that uncertainty is greater when predictions from different models do not agree. A calibrated, but unvalidated, ensemble was applied to the 475 km‐long coastline of Southern California (USA), with 4 SLR scenarios of 0.5, 0.93, 1.5, and 2 m by 2100. Results suggest that future retreat rates could increase relative to mean historical rates by more than two‐fold for the higher SLR scenarios, causing an average total land loss of 19 – 41 m by 2100. However, model uncertainty ranges from +/‐ 5 – 15 m, reflecting the inherent difficulties of projecting cliff retreat over multiple decades. To enhance ensemble performance, future work could include weighting each model by its skill in matching observations in different morphological settings

Journal of Geophysical Research F: Earth Surface↗
USGS2019

The MTPy software package for magnetotelluric data analysis and visualisation

The magnetotelluric (MT) method is increasingly being applied to a wide variety of geoscience problems. However, the software available for MT data analysis and interpretation is still very limited in comparison to many of the more mature geophysical methods such as the gravity, magnetic or seismic reflection methods. MTPy is an open source Python package to assist with MT data processing, analysis, modelling, visualization and interpretation. It was initiated at the University of Adelaide in 2013 as a means to store and share Python code amongst the MT community (Krieger & Peacock, 2014). Here we provide an overview of the software and describe recent developments to MTPy. These include new functionality and a clean up and standardization of the source code, as well as the addition of an integrated testing suite, documentation, and examples in order to facilitate the use of MT in the wider geophysics community.

Journal of Open Source Software↗