Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Seismology”

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
USGS2015

A comparison of methods to estimate seismic phase delays--Numerical examples for coda wave interferometry

Time-shift estimation between arrivals in two seismic traces before and after a velocity perturbation is a crucial step in many seismic methods. The accuracy of the estimated velocity perturbation location and amplitude depend on this time shift. Windowed cross correlation and trace stretching are two techniques commonly used to estimate local time shifts in seismic signals. In the work presented here, we implement Dynamic Time Warping (DTW) to estimate the warping function – a vector of local time shifts that globally minimizes the misfit between two seismic traces. We illustrate the differences of all three methods compared to one another using acoustic numerical experiments. We show that DTW is comparable to or better than the other two methods when the velocity perturbation is homogeneous and the signal-to-noise ratio is high. When the signal-to-noise ratio is low, we find that DTW and windowed cross correlation are more accurate than the stretching method. Finally, we show that the DTW algorithm has better time resolution when identifying small differences in the seismic traces for a model with an isolated velocity perturbation. These results impact current methods that utilize not only time shifts between (multiply) scattered waves, but also amplitude and decoherence measurements. DTW is a new tool that may find new applications in seismology and other geophysical methods (e.g., as a waveform inversion misfit function).

Geophysical Journal International↗
USGS2023

Inconsistent citation of the Global Seismographic Network in scientific publications

The highly used Global Seismographic Network (GSN) is a pillar of the seismological research community and contributes to numerous groundbreaking publications. Despite its wide recognition, this survey found that the GSN is not consistently acknowledged in scientific literature and is underrepresented by roughly a factor of 3 in citation searches. Publication tracking is a key metric that factors into operational decisions and funding support for the network; thus, consistent and proper citation of the GSN is important. This study not only serves as a reminder for researchers using GSN observations to cite the network’s digital object identifiers (DOIs) but also promotes a community‐wide conversation among researchers, journal editors, network operators, and other stakeholders regarding more standardized policies and review processes to ensure seismic networks are properly and consistently recognized for their contributions to research.

Seismological Research Letters↗
USGS1991

Don’t fence us in

When I was a graduate student around 1950 I used to read the entire Bulletin of the Seismological Society of America. it was a pwoerful and inspiring educational experience, with an effect quite different from that of the more usual process of looking up a few articles in the chain of references in a subject of current interest. Reading the entire journal reveals how ideas, techniques, and seismologists appear and evolve. It is likely the best substitute for a firsthand personal experience with the early development of the field. And in spite of, or perhaps because of, the missteps, the wasted effort, and the lack of sophistication that those first volumes reveal, the reader can sense the opportunity and be inspired by the vibrancy of the young subject.

Earthquakes & Volcanoes (USGS)↗
USGS2009

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2011

Coulomb stress change sensitivity due to variability in mainshock source models and receiving fault parameters: A case study of the 2010-2011 Christchurch, New Zealand, earthquakes

Strong aftershocks following major earthquakes present significant challenges for infrastructure recovery as well as for emergency rescue efforts. A tragic instance of this is the 22 February 2011 M w 6.3 Christchurch aftershock in New Zealand, which caused more than 100 deaths while the 2010 M w 7.1 Canterbury mainshock did not cause a single fatality (Figure 1). Therefore, substantial efforts have been directed toward understanding the generation mechanisms of aftershocks as well as mitigating hazards due to aftershocks. Among these efforts are the prediction of strong aftershocks, earthquake early warning, and aftershock probability assessment. Zhang et al. (1999) reported a successful case of strong aftershock prediction with precursory data such as changes in seismicity pattern, variation of b-value, and geomagnetic anomalies. However, official reports of such successful predictions in geophysical journals are extremely rare, implying that deterministic prediction of potentially damaging aftershocks is not necessarily more scientifically feasible than prediction of mainshocks.

Christchurch↗
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↗
USGS2012

The new IASPEI standards for determining magnitudes from digital data and their relation to classical magnitudes

Why there is a need for measurement standards of magnitudes: In October 2005, the Commission on Seismic Observation and Interpretation of the International Association of Seismology and Physics of the Earth´s Interior (IASPEI) adopted the summary recommendations made by the IASPEI Working Group on Magnitudes on new measurement standards for widely used local, regional and teleseismic magnitude scales (IASPEI, 2005). These recommendations have recently been refined and detailed (IASPEI, 2013) and a final scientific report, to be published in a reputable international journal, is currently under preparation.

Book chapter↗
USGS2001

Deformation across the Pacific-North America plate boundary near San Francisco, California

We have detected a narrow zone of compression between the Coast Ranges and the Great Valley, and we have estimated slip rates for the San Andreas, Rodgers Creek, and Green Valley faults just north of San Francisco. These results are based on an analysis of campaign and continuous Global Positioning System (GPS) data collected between 1992 and 2000 in central California. The zone of compression between the Coast Ranges and the Great Valley is 25 km wide. The observations clearly show 3.8±1.5 mm yr −1 of shortening over this narrow zone. The strike slip components are best fit by a model with 20.8±1.9 mm yr −1 slip on the San Andreas fault, 10.3±2.6 mm yr −1 on the Rodgers Creek fault, and 8.1±2.1 mm yr −1 on the Green Valley fault. The Pacific-Sierra Nevada-Great Valley motion totals 39.2±3.8 mm yr −1 across a zone that is 120 km wide (at the latitude of San Francisco). Standard deviations are one σ. The geodetic results suggest a higher than geologic rate for the Green Valley fault. The geodetic results also suggest an inconsistency between geologic estimates of the San Andreas rate and seismologic estimates of the depth of locking on the San Andreas fault. The only convergence observed is in the narrow zone along the border between the Great Valley and the Coast Ranges.

California↗
USGS2024

Using open-science workflow tools to produce SCEC CyberShake physics-based probabilistic seismic hazard models

The Statewide (formerly Southern) California Earthquake Center (SCEC) conducts multidisciplinary earthquake system science research that aims to develop predictive models of earthquake processes, and to produce accurate seismic hazard information that can improve societal preparedness and resiliency to earthquake hazards. As part of this program, SCEC has developed the CyberShake platform, which calculates physics-based probabilistic seismic hazard analysis (PSHA) models for regions with high-quality seismic velocity and fault models. The CyberShake platform implements a sophisticated computational workflow that includes over 15 individual codes written by 6 developers. These codes are heterogeneous, ranging from short-running high-throughput serial CPU codes to large, long-running, parallel GPU codes. Additionally, CyberShake simulation campaigns are computationally extensive, typically producing tens of terabytes of meaningful scientific data and metadata over several months of around-the-clock execution on leadership-class supercomputers. To meet the needs of the CyberShake platform, we have developed an extreme-scale workflow stack, including the Pegasus Workflow Management System, HTCondor, Globus, and custom tools. We present this workflow software stack and identify how the CyberShake platform and supporting tools enable us to meet a variety of challenges that come with large-scale simulations, such as automated remote job submission, data management, and verification and validation. This platform enabled us to perform our most recent simulation campaign, CyberShake Study 22.12, from December 2022 to April 2023. During this time, our workflow tools executed approximately 32,000 jobs, and used up to 73% of the Summit system at Oak Ridge Leadership Computing Facility. Our workflow tools managed about 2.5 PB of total temporary and output data, and automatically staged 19 million output files totaling 74 TB back to archival storage on the University of Southern California's Center for Advanced Research Computing systems, including file-based relational data and large binary files to efficiently store millions of simulated seismograms. CyberShake extreme-scale workflows have generated simulation-based probabilistic seismic hazard models that are being used by seismological, engineering, and governmental communities.

Frontiers Earth Science Journal↗
USGS2015

Tidal and seasonal variations in calving flux observed with passive seismology

The seismic signatures of calving events, i.e., calving icequakes, offer an opportunity to examine calving variability with greater precision than is available with other methods. Here using observations from Yahtse Glacier, Alaska, we describe methods to detect, locate, and characterize calving icequakes. We combine these icequake records with a coincident, manually generated record of observed calving events to develop and validate a statistical model through which we can infer iceberg sizes from the properties of calving icequakes. We find that the icequake duration is the single most significant predictor of an iceberg's size. We then apply this model to 18 months of seismic recordings and find elevated iceberg calving flux during the summer and fall and a pronounced lull in calving during midwinter. Calving flux is sensitive to semidiurnal tidal stage. Large calving events are tens of percent more likely during falling and low tides than during rising and high tides, consistent with a view that deeper water has a stabilizing influence on glacier termini. Multiple factors affect the occurrence of mechanical fractures that ultimately lead to iceberg calving. At Yahtse Glacier, seismology allows us to demonstrate that variations in the rate of submarine melt are a dominant control on iceberg calving rates at seasonal timescales. On hourly to daily timescales, tidal modulation of the normal stress against the glacier terminus reveals the nonlinear glacier response to changes in the near-terminus stress field.

Alaska↗
USGS1990

Early spectral shaping boosts data quality

Reflection seismology is a method of conveying information about the interior of the earth to the surface of the earth. Most commonly the source and receivers are located at or near the surface. This paper is a conceptual discussion of how one may maximize the amount of information received from reflection paths between sources and receivers. Some of the basic tenets of information theory are discussed in simple terms for illustration purposes.

Oil & Gas Journal↗
USGS1994

The 1989-1990 eruptions of Redoubt Volcano: an introduction

Redoubt Volcano, located on the west side of Cook Inlet in south-central Alaska, erupted explosively on over 20 separate occasions between December 14, 1989 and April 21, 1990. Fourteen lava domes were emplaced in the summit area, thirteen of which were subsequently destroyed. The eruption caused economic losses estimated at over $160,000,000 making this the second most costly eruption in U.S. history. This economic impact provided the impetus for a integrated comprehensive account of an erupting volcano using both modern and classical research and modern techniques which in turn led to advances in eruption monitoring and interpretation. Research on such topics as dome formation and collapse and the resulting pyroclastic flows, elutriated ash, lightning, tephra, and flooding was blended with the rapid communication of associated hazards to a large user group. The seismology successes in predicting and monitoring eruption dynamics were due in part to (1) the recognition of long-period seismic events as indicators of the readiness of the volcano to erupt, and (2) to the development of new tools that allowed the seismicity to be assessed instantaneously. Integrated studies of the petrology of erupted products and volatile content over time gave clues as to the progress of the eruption towards completion.

Alaska↗
USGS2021

A subset of CyberShake ground-motion time series for response-history analysis

This manuscript describes a subset of CyberShake numerically simulated ground motions that were selected and vetted for use in engineering response-history analyses. Ground motions were selected that have seismological properties and response spectra representative of conditions in the Los Angeles area, based on disaggregation of seismic hazard. Ground motions were selected from millions of available time series and were reviewed to confirm their suitability for response-history analysis. The processes used to select the time series, the characteristics of the resulting data, and the provided documentation are described in this article. The resulting data and documentation are available electronically.

California↗
USGS2023

Strength recovery in quartzite is controlled by changes in friction in experiments at hydrothermal conditions up to 200°C

The rate of fault zone restrengthening between earthquakes can be influenced by both frictional and cohesive healing processes. Friction is dependent on effective normal stress while cohesion is independent of normal stress, potentially explaining—in part—the lack of depth dependence of earthquake stress drops. Although amenable to laboratory testing, few studies have systematically addressed the normal stress dependence of restrengthening rate. This is partially due to difficulty in separating relative contributions of friction and cohesion in recovery of fault strength. We present results from a series of slide-hold-slide tests on thin layers (≤10 𝜇m) of ultrafine quartz gouge that develop during shearing of initially bare-surface quartzite. Tests were conducted at 10 MPa constant pore pressure, 20–200 MPa constant effective normal stress, and temperatures of 22°–200°C. Restrengthening, defined as the difference between peak shear stress measured after resumption of sliding and steady-state sliding shear stress, increases with the log of hold duration. The 200°C healing rate, 0.014 per e-fold increase in time, is comparable to that determined from seismological observations along the Calaveras Fault, California. Construction of Mohr-Coulomb failure envelopes shows that changes in cohesion are small (<1 MPa) and independent of hold durations to 10 5 s, indicating that the increased strength is due to changes in the friction coefficient. These experimental results are inconsistent with the hypothesis that cohesive healing explains the depth independence of earthquake stress drop, but higher temperatures, longer time-scales, and more complex mineralogy could facilitate cohesive healing in natural fault systems.

Journal of Geophysical Research - Solid Earth↗
USGS2025

A crustal thermal model of the conterminous U.S. constrained by multiple data sets: A Monte-Carlo approach

The thermal structure of the continental crust plays a critical role in understanding its elastic and rheologic properties as well as its dynamic processes. Thermal parameter data sets on continental scales have been used to constrain the crustal thermal structure, including both the direct (e.g. temperature, heat flux and heat conductivity measured at the surface) and indirect (e.g. seismically derived Mohorovičić discontinuity (Moho) temperature, geomagnetically derived Curie depth) observations. In this study, we present a new continental scale crustal heat generation model with additional information from seismologically inferred crustal composition. Together with previous direct and indirect thermal parameter data sets in the conterminous United States, we use the new crustal heat generation model to construct a 3-D crustal temperature model under a newly developed Bayesian framework. Specifically, we first derive profiles of crustal heat generation based on an empirical geochemical relationship at 1683 locations where seismologically derived crustal composition information is available. Then for each of these locations, the average heat generation values in the upper, middle and lower crust are combined with other thermal parameters through a Markov Chain Monte-Carlo inversion for a conductive, vertically smooth temperature profile. The results, posterior distributions of temperature profiles, are used to generate a 3-D crustal thermal model with the uncertainties systematically assessed. The new temperature model overall exhibits similar patterns to that from the U.S. Geological Survey National Crustal Model, but also reduces possible biases and the model's dependence on a single thermal parameter.

conterminous United States↗
USGS2021

River terrace evidence of tectonic processes in the eastern North American plate interior, South Anna River, Virginia

We show that long-recognized seismicity in the central Virginia seismic zone of the eastern North American intraplate setting arises primarily from tectonic processes predicted by new, fully coupled plate tectonic geodynamic models. The study leverages much new geophysical and geologic data following the 2011 Mineral, Virginia, earthquake that ruptured a steeply dipping, northwest-verging reverse fault traversed by the South Anna River. The data are primarily assembled from a flight of six fluvial terrace geomorphic markers identified and correlated on texture, relative weathering, and numeric ages including one terrestrial cosmogenic nuclide (TCN) profile and 30 luminescence dates. Terrace thickness, stratigraphic age models, and incision rates downstream and upstream of the 2011 rupture are different. Long-term river incision rates of ∼25–30 m/My are superimposed on regional TCN-determined erosion rates of ∼8.5 m/My; however, there are at least 10 m of tectonically driven incision in the epicentral region at rates of ∼30–94 m/My. The inferred deformation resembles a hanging wall anticline above a blind reverse fault with a diffuse overlying carapace of minor brittle faults, an interpretation supported by seismology as well as bedrock and saprolite mapped across the epicentral region. These results are further supported by channel metrics that show nonuniform channel steepness ( k sn ) and a predicted steady-state channel elevation different from the actual channel elevation across the epicentral region. If all of the observed deformation is a consequence of the fault that ruptured in 2011, the recurrence interval of Mineral-sized events would be ∼5.5 ky.

Virginia↗
USGS2000

The crustal thickness of Australia

We investigate the crustal structure of the Australian continent using the temporary broadband stations of the Skippy and Kimba projects and permanent broadband stations. We isolate near-receiver information, in the form of crustal P-to-S conversions, using the receiver function technique. Stacked receiver functions are inverted for S velocity structure using a Genetic Algorithm approach to Receiver Function Inversion (GARFI). From the resulting velocity models we are able to determine the Moho depth and to classify the width of the crust-mantle transition for 65 broadband stations. Using these results and 51 independent estimates of crustal thickness from refraction and reflection profiles, we present a new, improved, map of Moho depth for the Australian continent. The thinnest crust (25 km) occurs in the Archean Yilgarn Craton in Western Australia; the thickest crust (61 km) occurs in Proterozoic central Australia. The average crustal thickness is 38.8 km (standard deviation 6.2 km). Interpolation error estimates are made using kriging and fall into the range 2.5–7.0 km. We find generally good agreement between the depth to the seismologically defined Moho and xenolith-derived estimates of crustal thickness beneath northeastern Australia. However, beneath the Lachlan Fold Belt the estimates are not in agreement, and it is possible that the two techniques are mapping differing parts of a broad Moho transition zone. The Archean cratons of Western Australia appear to have remained largely stable since cratonization, reflected in only slight variation of Moho depth. The largely Proterozoic center of Australia shows relatively thicker crust overall as well as major Moho offsets. We see evidence of the margin of the contact between the Precambrian craton and the Tasman Orogen, referred to as the Tasman Line.

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↗