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 19 recordsLinked to original sources

Designs and test results for three new rotational sensors

We discuss the designs and testing of three rotational seismometer prototypes developed at the Institute of Geophysics, Academy of Sciences (Prague, Czech Republic). Two of these designs consist of a liquid-filled toroidal tube with the liquid as the proof mass and providing damping; we tested the piezoelectric and pressure transduction versions of this torus. The third design is a wheel-shaped solid metal inertial sensor with capacitive sensing and magnetic damping. Our results from testing in Prague and at the Albuquerque Seismological Laboratory of the US Geological Survey of transfer function and cross-axis sensitivities are good enough to justify the refinement and subsequent testing of advanced prototypes. These refinements and new testing are well along.

Journal of Seismology↗

A review of near-surface QS estimation methods using active and passive sources

Seismic attenuation and the associated quality factor ( Q ) have long been studied in various sub-disciplines of seismology, ranging from observational and engineering seismology to near-surface geophysics and soil/rock dynamics with particular emphasis on geotechnical earthquake engineering and engineering seismology. Within the broader framework of seismic site characterization, various experimental techniques have been adopted over the years to measure the near-surface shear-wave quality factor ( Q S ). Common methods include active- and passive-source recording techniques performed at the free surface of soil deposits and within boreholes, as well as laboratory tests. This paper intends to provide an in-depth review of what Q is and, in particular, how Q S is estimated in the current practice. After motivating the importance of this parameter in seismology, we proceed by recalling various theoretical definitions of Q and its measurement through laboratory tests, considering various deformation modes, most notably Q P and Q S . We next provide a review of the literature on Q S estimation methods that use data from surface and borehole sensor recordings. We distinguish between active- and passive-source approaches, along with their pros and cons, as well as the state-of-the-practice and state-of-the-art. Finally, we summarize the phenomena associated with the high-frequency shear-wave attenuation factor (kappa) and its relation to Q , as well as other lesser-known attenuation parameters.

Journal of Seismology↗

Key elements of regional seismic velocity models for long period ground motion simulations

Regional 3-D seismic velocity models used for broadband strong motion simulations must include compressional-wave velocity (Vp), shear-wave velocity (Vs), intrinsic attenuation (Qp, Qs), and density. Vs and Qs are the most important of these parameters because the strongest ground motions are generated chiefly by shear- and surface-wave arrivals. Because Vp data are more common than Vs data, many researchers first develop a Vp model and convert it to a Vs model. I describe recent empirical relations between Vs, Vp, Qs, Qp, and density that allow velocity models to be rapidly and accurately calculated. ?? Springer Science+Business Media B.V. 2007.

Journal of Seismology↗

3D crustal structure and long-period ground motions from a M9.0 megathrust earthquake in the Pacific Northwest region

We have developed a community velocity model for the Pacific Northwest region from northern California to southern Canada and carried out the first 3D simulation of a M w 9.0 megathrust earthquake rupturing along the Cascadia subduction zone using a parallel supercomputer. A long-period (<0.5 Hz) source model was designed by mapping the inversion results for the December 26, 2004 Sumatra–Andaman earthquake (Han et al., Science 313(5787):658–662, 2006) onto the Cascadia subduction zone. Representative peak ground velocities for the metropolitan centers of the region include 42 cm/s in the Seattle area and 8–20 cm/s in the Tacoma, Olympia, Vancouver, and Portland areas. Combined with an extended duration of the shaking up to 5 min, these long-period ground motions may inflict significant damage on the built environment, in particular on the highrises in downtown Seattle.

Journal of Seismology↗

Uncertainty estimates in broadband seismometer sensitivities using microseisms

The midband sensitivity of a seismic instrument is one of the fundamental parameters used in published station metadata. Any errors in this value can compromise amplitude estimates in otherwise high-quality data. To estimate an upper bound in the uncertainty of the midband sensitivity for modern broadband instruments, we compare daily microseism (4- to 8-s period) amplitude ratios between the vertical components of colocated broadband sensors across the IRIS/USGS (network code IU) seismic network. We find that the mean of the 145,972 daily ratios used between 2002 and 2013 is 0.9895 with a standard deviation of 0.0231. This suggests that the ratio between instruments shows a small bias and considerable scatter. We also find that these ratios follow a standard normal distribution ( R 2 &thinsp;=&thinsp;0.95442), which suggests that the midband sensitivity of an instrument has an error of no greater than &plusmn;6 % with a 99 % confidence interval. This gives an upper bound on the precision to which we know the sensitivity of a fielded instrument.

Journal of Seismology↗

Was the Mw 7.5 1952 Kern County, California, earthquake induced (or triggered)?

Several recent studies have presented evidence that significant induced earthquakes occurred in a number of oil-producing regions during the early and mid-twentieth century related to either production or wastewater injection. We consider whether the 21 July 1952 M w 7.5 Kern County earthquake might have been induced by production in the Wheeler Ridge oil field. The mainshock, which was not preceded by any significant foreshocks, occurred 98 days after the initial production of oil in Eocene strata at depths reaching 3 km, within ~1 km of the White Wolf fault (WWF). Based on this spatial and temporal proximity, we explore a potential causal relationship between the earthquake and oil production. While production would have normally be expected to have reduced pore pressure, inhibiting failure on the WWF, we present an analytical model based on industry stratigraphic data and best estimates of parameters whereby an impermeable splay fault adjacent to the main WWF could plausibly have blocked direct pore pressure effects, allowing the poroelastic stress change associated with production to destabilize the WWF, promoting initial failure. This proof-of-concept model can also account for the 98-day delay between the onset of production and the earthquake. While the earthquake clearly released stored tectonic stress, any initial perturbation on or near a major fault system can trigger a larger rupture. Our proposed mechanism provides an explanation for why significant earthquakes are not commonly induced by production in proximity to major faults.

California↗

Temporal patterns of induced seismicity in Oklahoma revealed from multi-station template matching

Over the past decade, Oklahoma became the most seismically active region of the mid-Continental USA as a result of industry operations. However, seismic network limitations and completeness of earthquake catalogs have restricted the types of analyses that can be performed. By applying multi-station template matching on the 23,889 cataloged earthquakes in Oklahoma and Southern Kansas between late-2008 and 2016, we increased the number of detected earthquakes to 209,409 events. While the improved catalog produced an order of magnitude events than the original catalog, the frequency-magnitude distribution remains similar to the original catalog. We found that the coefficient of variation of interevent times in small spatial bins tends to spatially correlate with the location of M ≥ 4 earthquakes. The improved catalog reveals the pervasiveness of swarm-like patterns in seismicity across the entire study region. The rapid increase in seismicity rate of these swarms in 2013 coincided with a reduction in the calculated p values (power law decay rates) before and after larger events. We also used the catalog to revisit the temporal patterns in the four M ≥ 5 sequences, finding more active foreshock behavior than previously recognized and variations in aftershock behavior. When compared against poroelastic stress models for the Pawnee and Fairview sequences, the catalog shows an improved correlation with stress that accounts for variable-rate injection, supporting the conclusion that injection rate is an important contributor to seismic hazard.

Oklahoma↗

Seismic site characterization with shear wave (SH) reflection and refraction methods

Reflection and critically refracted seismic methods use traveltime measurements of body waves propagating between a source and a series of receivers on the ground surface to calculate subsurface velocities. Body wave energy is refracted or reflected at boundaries where there is a change in seismic impedance, defined as the product of material density and seismic velocity. This article provides practical guidance on the use of horizontally propagating shear wave (SH-wave) refraction and reflection methods to determine shear wave velocity as a function of depth for near-surface seismic site characterizations. Method principles and the current state of engineering practice are reviewed, along with discussions of limitations and uncertainty assessments. Typical data collection procedures are described using basic survey equipment, along with information on more advanced applications and emerging technologies. Eight case studies provide examples of the techniques in real-world seismic site characterizations performed in a variety of geological settings.

Journal of Seismology↗

Microtremor array method using spatial autocorrelation analysis of Rayleigh‑wave data

Microtremor array measurements (MAM) and passive surface wave methods in general, have been increasingly used to non-invasively estimate shear-wave velocity structures (Vs) for various purposes. The methods estimate dispersion curves and invert them for retrieving S-wave velocity profiles. This paper summarizes principles, limitations, data collection and processing methods. It intends to enable students and practitioners to understand the principles needed to plan a microtremor array investigation, record and process the data, and evaluate the quality of investigation result. The paper focuses on the spatial autocorrelation (SPAC) processing method among microtremor array processing methods because of its relatively simple calculation and stable applicability.

Journal of Seismology↗

A review of the microtremor horizontal-to-vertical spectral ratio (MHVSR) method

The single-station microtremor horizontal-to-vertical spectral ratio (MHVSR) method was initially proposed to retrieve the site amplification function and its resonance frequencies produced by unconsolidated sediments overlying high-velocity bedrock. Presently, MHVSR measurements are predominantly conducted to obtain an estimate of the fundamental site frequency at sites where a strong subsurface impedance contrast exists. Of the earthquake site characterization methods presented in this special issue, the MHVSR method is the furthest behind in terms of consensus towards standardized guidelines and commercial use. The greatest challenges to an international standardization of MHVSR acquisition and analysis are (1) the what — the underlying composition of the microtremor wavefield is site-dependent, and thus, the appropriate theoretical (forward) model for inversion is still debated; and (2) the how — many factors and options are involved in the data acquisition, processing, and interpretation stages. This paper reviews briefly a historical development of the MHVSR technique and the physical basis of an MHVSR (the what ). We then summarize recommendations for MHVSR acquisition and analysis (the how ). Specific sections address MHVSR interpretation and uncertainty assessment.

Journal of Seismology↗

Flexible multimethod approach for seismic site characterization

We describe the flexible multimethod seismic site characterization technique for obtaining shear-wave velocity ( V S ) profiles and derivative information, such as the time-averaged V S of the upper 30 m ( V S30 ). Simply stated, the multimethod approach relies on the application of multiple independent noninvasive site characterization acquisition and analysis techniques utilized in a flexible field-based approach through a decision process based on primary factors such as the investigator’s available equipment, field logistics (e.g., available array deployment space, urban versus rural), geologic setting, and funding among other primary factors. In a multimethod approach, sites ideally should be characterized using both active and passive noninvasive (i.e., no drilling and only minimal disturbance to the ground surface) methods. Almost without exception, we recommend the use of at least one active-source technique for analyzing surface waves, which in the current state of practice would include one or more of the following: spectral analysis of surface waves (SASW; commonly Rayleigh waves) and multichannel analysis of surface waves (MASW; Rayleigh and/or Love waves). In addition, passive-source microtremor array methods (MAMs) using linear (one-dimensional) and two-dimensional arrays may be appropriate or even required for characterizing a given site. Similarly, the microtremor horizontal-to-vertical spectral ratio (mHVSR) method can be valuable for identifying shallow rock sites. Finally, P- and SH-wave seismic refraction methods warrant consideration, especially at rock and shallow-rock sites.

Journal of Seismology↗

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↗

Site-specific amplifications in Northwestern Turkiye: A generic approach

In this study, we derive generic site amplification functions applicable to Northwestern Türkiye and follow guidance on site classifications established by the United States National Earthquake Hazards Reduction Program (NEHRP). We employ the one-dimensional (1-D) equivalent linear ground response analysis method and focus on recordings from a select suite of free-field strong motion (SM) stations. These particular SM stations are associated with a comprehensive geotechnical site conditions dataset comprising 76 shear-wave velocity profiles. Despite their varying resolutions and depths, this diversity allows for a more representative analysis across different site conditions, enhancing the robustness of our findings. For analyses, we utilize a robust set of global ground motion records, selected from both within and outside of Türkiye, to serve as strong and weak input bedrock motions. Site amplification is determined by propagating these bedrock records through 76 site-specific Afet ve Acil Durum Yönetimi Başkanlığı (AFAD) station profiles and calculating the spectral acceleration ratio between the ground surface and the input motion. Our analyses reveal substantial disparities between the amplification functions derived from 1-D site response analyses using strong and weak input motions. These disparities are central to our analysis and inform the subsequent comparison with site amplification functions reported in the literature. We compare the generic site amplification functions derived in this study with those reported in the literature for regions with similar geological and seismic conditions to Northwestern Türkiye. For example, site amplification functions for NEHRP Site Classes C and D have been documented for Greece and the United States utilizing frequency-dependent methods such as the quarter-wavelength approximation. Our findings reveal that the generic site amplification functions derived in this study exhibit differences in peak amplitudes and predominant frequencies compared to those in the literature, reflecting regional variations in geotechnical and seismic properties. This study represents an initial effort to formulate generic site amplification functions for Türkiye, with a particular emphasis on the Northwestern region. A thorough literature review indicated no prior studies have proposed such functions specifically for this area.

Northwestern Turkiye↗

The initial subevent of the 1994 Northridge, California, earthquake: Is earthquake size predictable?

We examine the initial subevent (ISE) of the M 6.7, 1994 Northridge, California, earthquake in order to discriminate between two end-member rupture initiation models: the ‘preslip’ and ‘cascade’ models. Final earthquake size may be predictable from an ISE's seismic signature in the preslip model but not in the cascade model. In the cascade model ISEs are simply small earthquakes that can be described as purely dynamic ruptures. In this model a large earthquake is triggered by smaller earthquakes; there is no size scaling between triggering and triggered events and a variety of stress transfer mechanisms are possible. Alternatively, in the preslip model, a large earthquake nucleates as an aseismically slipping patch in which the patch dimension grows and scales with the earthquake's ultimate size; the byproduct of this loading process is the ISE. In this model, the duration of the ISE signal scales with the ultimate size of the earthquake, suggesting that nucleation and earthquake size are determined by a more predictable, measurable, and organized process. To distinguish between these two end-member models we use short period seismograms recorded by the Southern California Seismic Network. We address questions regarding the similarity in hypocenter locations and focal mechanisms of the ISE and the mainshock. We also compare the ISE's waveform characteristics to those of small earthquakes and to the beginnings of earthquakes with a range of magnitudes. We find that the focal mechanisms of the ISE and mainshock are indistinguishable, and both events may have nucleated on and ruptured the same fault plane. These results satisfy the requirements for both models and thus do not discriminate between them. However, further tests show the ISE's waveform characteristics are similar to those of typical small earthquakes in the vicinity and more importantly, do not scale with the mainshock magnitude. These results are more consistent with the cascade model.

Journal of Seismology↗

Jack Boatwright’s last science

In the months before he passed away, Jack Boatwright asked us to review a manuscript on source properties, specifically stress drop, of earthquakes in northeastern North America (NENA). This manuscript originated in research funded by the U.S. Nuclear Regulatory Commission (NRC), described in his final report to NRC and published as USGS Open-File Report 2018-1073 (Boatwright, 2018). We wish to call attention to Boatwright (2018) for those who may not be aware of what Jack was unable to put into journal publications.

Seismological Research Letters↗

The Survey’s first venture into seismology

The occurrence of two small but widely felt earthquakes in the Eastern United States in 1884 led an editor of the research journal Science to suggest, in the October 3 issue of the magazine, that an "earthquake club" be formed. Its purpose was so that observers and students of "this branch of physical geography" could become acquainted with each other and thus insure the proper collection and collation of their observations. Because this was its intent, he though the club should be "unembarrassed by formal regulations"; consequently, there is no record of whether such an organziation ever existed. But the other suggestion, thata of collecting earthquake records was acted on promptly.

Earthquake Information Bulletin (USGS)↗

CRUST 5.1: A global crustal model at 5° x 5°

We present a new global model for the Earth's crust based on seismic refraction data published in the period 1948-1995 and a detailed compilation of ice and sediment thickness. An extensive compilation of seismic refraction measurements has been used to determine the crustal structure on continents and their margins. Oceanic crust is modeled with both a standard model for normal oceanic crust, and variants for nonstandard regions, such as oceanic plateaus. Our model (CRUST 5.1) consists of 2592 5° × 5° tiles in which the crust and uppermost mantle are described by eight layers: (1) ice, (2) water, (3) soft sediments, (4) hard sediments, (5) crystalline upper, (6) middle, (7) lower crust, and (8) uppermost mantle. Topography and bathymetry are adopted from a Mandard database (ETOPO-5). Compressional wave velocity in each layer is based on field measurements, and shear wave velocity and density are estimated using recently published empirical V p -V s and V p -density relationships. The crustal model differs from previous models in that (1) the thickness and seismic/density structure of sedimentary basins is accounted for more completely, (2) the velocity structure of unmeasured regions is estimated using statistical averages that are based on a significantly larger database of crustal structure, (3) the compressional wave, shear wave, and density structure have been explicitly specified using newly available constraints from field and laboratory studies. Thus - this global crustal model is based on substantially more data than previous models and differs from them in many important respects. A new map of the thickness of the Earth's crust is presented, and we illustrate the application of this model by using it to provide the crustal correction for surface wave phase velocity maps. Love waves at 40 s are dominantly sensitive to crustal structure, and there is a very close correspondence between observed phase velocities at this period and those predicted by CRUST 5.1. We find that the application of crustal corrections to long-period (167 s) Rayleigh waves significantly increases the variance in the phase velocity maps and strengthens the upper mantle , velocity anomalies beneath stable continental regions. A simple calculation of crustal isostacy indicates significant lateral variations in upper mantle density. The model CRUST 5.1 provides a complete description of the physical properties of the Earth's crust at a scale of 5° × 5° and can be used for a wide range of seismological and nonseismological problems.

Journal of Geophysical Research B: Solid Earth↗

Calving seismicity from iceberg-sea surface interactions

Iceberg calving is known to release substantial seismic energy, but little is known about the specific mechanisms that produce calving icequakes. At Yahtse Glacier, a tidewater glacier on the Gulf of Alaska, we draw upon a local network of seismometers and focus on 80 hours of concurrent, direct observation of the terminus to show that calving is the dominant source of seismicity. To elucidate seismogenic mechanisms, we synchronized video and seismograms to reveal that the majority of seismic energy is produced during iceberg interactions with the sea surface. Icequake peak amplitudes coincide with the emergence of high velocity jets of water and ice from the fjord after the complete submergence of falling icebergs below sea level. These icequakes have dominant frequencies between 1 and 3 Hz. Detachment of an iceberg from the terminus produces comparatively weak seismic waves at frequencies between 5 and 20 Hz. Our observations allow us to suggest that the most powerful sources of calving icequakes at Yahtse Glacier include iceberg-sea surface impact, deceleration under the influence of drag and buoyancy, and cavitation. Numerical simulations of seismogenesis during iceberg-sea surface interactions support our observational evidence. Our new understanding of iceberg-sea surface interactions allows us to reattribute the sources of calving seismicity identified in earlier studies and offer guidance for the future use of seismology in monitoring iceberg calving.

Journal of Geophysical Research F: Earth Surface↗