USGS ScienceSearch

SEARCH · USGS Science

Results for “Bulletin Seismological Society America”

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

16th century European earthquakes described in some contemporary woodcuts

In the middle of the 15 th century, shortly after Gutenberg’s invention of printing using individual lead type, the first illustrated broadsides (or broadsheets) began appearing in southern Germany. Usually printed on one side of a sheet of paper, they consisted of a woodcut illustration, sometimes colored, either by hand or by stencil, and several columns of text, often in verse. The appeal of the publications to the mass reader was often stimulated by sensationalism in word and picture, somewhat like that purveyed today’s tabloid press. What follows are woodcuts and portions of the text from 16 th to 17 th century broadsheets that describe earthquakes and volcanic eruptions. This article is an excerpt from Deresiewicz, H., 1982, Some 16th century European earthquakes as depicted in contemporary sources: Bulletin Seismological Society of America, vol. 72, p 507-523.

Earthquake Information Bulletin (USGS)

Shear-wave velocity site characterization in Oklahoma from joint inversion of multi-method surface seismic measurements: Implications for central U.S. Ground Motion Prediction

We analyze multimethod shear (SH)‐wave velocity ( ⁠ V S "> V S ) site characterization data acquired at three permanent and 25 temporary seismograph stations in Oklahoma that recorded M 4+ earthquakes within a 50 km hypocentral distance of at least one of the 2016 M 5.1 Fairview, M 5.8 Pawnee, or M 5.0 Cushing earthquakes to better constrain earthquake ground‐motion modeling in the region. We acquired active‐source seismic data for time‐averaged V S "> V S to 30 m depth ( V S "> V S V S 30 "> 30 ) at 28 seismograph stations near the Fairview, Pawnee, and Cushing epicentral areas. The SH‐wave refraction travel times coupled with Rayleigh‐ and Love‐wave phase velocity dispersion were extracted and modeled in a nonlinear least‐squares (L2) joint inversion to obtain a best‐fit 1D V S "> V S versus depth profile for each site. At a subset of sites where the preferred L2 inverse model did not optimally fit each of the Love, Rayleigh, and SH travel‐time datasets, we explore application of simulated annealing in a joint inversion to find a more global solution. V S "> V S V S 30 "> 30 values range from 262 to 807    m / s "> 807 m / s for the preferred measured (in situ) V S "> V S profiles, or National Earthquake Hazards Reduction Program (NEHRP) site class D to B, and are broadly comparable with estimates from previous data reports in the region. Site amplification estimates were calculated next from 1D SH transfer functions of the preferred V S "> V S VS profiles and then compared against observed horizontal‐to‐vertical spectral ratios (HVSRs) from nearby seismograph stations. We generally see good agreement between the predicted in situ model and the observed HVSR resonant frequencies, with nominal amplifications between 2 and 10 within the 2–15 Hz frequency band. Next, using 40 known in situ V S "> V S V S 30 "> 30 measurements in the region, we demonstrate that the in situ V S "> V S V S 30 "> 30 values improve the fit for selected suites of ground‐motion models (GMMs) for M 4+ earthquakes within a 50 km hypocentral distance when compared with proxy methods, arguing for future development of GMMs implementing in situ V S "> V S profiles.

Oklahoma

Epistemic uncertainty in California-wide synthetic seismicity simulations

The generation of seismicity catalogs on synthetic fault networks holds the promise of providing key inputs into probabilistic seismic-hazard analysis, for example, the coefficient of variation, mean recurrence time as a function of magnitude, the probability of fault-to-fault ruptures, and conditional probabilities for foreshock–mainshock triggering. I employ a seismicity simulator that includes the following ingredients: static stress transfer, viscoelastic relaxation of the lower crust and mantle, and vertical stratification of elastic and viscoelastic material properties. A cascade mechanism combined with a simple Coulomb failure criterion is used to determine the initiation, propagation, and termination of synthetic ruptures. It is employed on a 3D fault network provided by Steve Ward (unpublished data, 2009) for the Southern California Earthquake Center (SCEC) Earthquake Simulators Group. This all-California fault network, initially consisting of 8000 patches, each of ∼12 square kilometers in size, has been rediscretized into Graphic patches, each of ∼1 square kilometer in size, in order to simulate the evolution of California seismicity and crustal stress at magnitude M∼5–8. Resulting synthetic seismicity catalogs spanning 30,000 yr and about one-half million events are evaluated with magnitude-frequency and magnitude-area statistics. For a priori choices of fault-slip rates and mean stress drops, I explore the sensitivity of various constructs on input parameters, particularly mantle viscosity. Slip maps obtained for the southern San Andreas fault show that the ability of segment boundaries to inhibit slip across the boundaries (e.g., to prevent multisegment ruptures) is systematically affected by mantle viscosity.

California

Performance of Irikura recipe rupture model generator in earthquake ground motion simulations with Graves and Pitarka hybrid approach

We analyzed the performance of the Irikura and Miyake (Pure and Applied Geophysics 168(2011):85–104, 2011) (IM2011) asperity-based kinematic rupture model generator, as implemented in the hybrid broadband ground motion simulation methodology of Graves and Pitarka (Bulletin of the Seismological Society of America 100(5A):2095–2123, 2010), for simulating ground motion from crustal earthquakes of intermediate size. The primary objective of our study is to investigate the transportability of IM2011 into the framework used by the Southern California Earthquake Center broadband simulation platform. In our analysis, we performed broadband (0–20 Hz) ground motion simulations for a suite of M6.7 crustal scenario earthquakes in a hard rock seismic velocity structure using rupture models produced with both IM2011 and the rupture generation method of Graves and Pitarka (Bulletin of the Seismological Society of America, 2016) (GP2016). The level of simulated ground motions for the two approaches compare favorably with median estimates obtained from the 2014 Next Generation Attenuation-West2 Project (NGA-West2) ground motion prediction equations (GMPEs) over the frequency band 0.1–10 Hz and for distances out to 22 km from the fault. We also found that, compared to GP2016, IM2011 generates ground motion with larger variability, particularly at near-fault distances (<12 km) and at long periods (>1 s). For this specific scenario, the largest systematic difference in ground motion level for the two approaches occurs in the period band 1–3 s where the IM2011 motions are about 20–30% lower than those for GP2016. We found that increasing the rupture speed by 20% on the asperities in IM2011 produced ground motions in the 1–3 s bandwidth that are in much closer agreement with the GMPE medians and similar to those obtained with GP2016. The potential implications of this modification for other rupture mechanisms and magnitudes are not yet fully understood, and this topic is the subject of ongoing study. We concluded that IM2011 rupture generator performs well in ground motion simulations using Graves and Pitarka hybrid method. Therefore, we recommend it to be considered for inclusion into the framework used by the Southern California Earthquake Center broadband simulation platform.

Pure and Applied Geophysics

Software for inference of dynamic ground strains and rotations and their errors from short baseline array observations of ground motions

In two previous articles we presented a formulation for inferring the strains and rotations of the ground beneath a seismic array having a finite footprint. In this article we derive expressions for the error covariance matrices of the inferred strains and rotations, and we present software for the calculation of ground strains, rotations, and their variances from short baseline array ground-motion data.

Bulletin of the Seismological Society of America

Source parameters for small events associated with the 1986 North Palm Springs, California, earthquake determined using empirical Green functions

Using small events as empirical Green functions, which include path, site, and instrument effects on the P waveforms, source parameters were estimated for 25 M L 3.4 to 4.4 events associated with the 1986 North Palm Springs earthquake. The static stress drops ranged from 3 to 80 bars, for moments of 0.7 to 11 × 10 21 dyne-cm. There was a spatial pattern to the stress drops of the aftershocks which showed increasing values along the fault plane toward the northwest compared to relatively low values near the hypocenter of the mainshock. The highest values were outside the main area of slip, and are believed to reflect a loaded area of the fault that still has an higher level of stress which was not released during the main shock. The relatively high stress drop of a preshock in 1983, located close to the hypocenter of the main shock, is also interpreted as evidence that stress drops of magnitude 3 and 4 events can be used as indicators of the stress level in a particular area.

California

Regional correlations of V S 30 averaged over depths less than and greater than 30 meters

Using velocity profiles from sites in Japan, California, Turkey, and Europe, we find that the time-averaged shear-wave velocity to 30 m ( V S 30 ), used as a proxy for site amplification in recent ground-motion prediction equations (GMPEs) and building codes, is strongly correlated with average velocities to depths less than 30 m ( V Sz , with z being the averaging depth). The correlations for sites in Japan (corresponding to the KiK-net network) show that V Sz is systematically larger for a given V Sz than for profiles from the other regions. The difference largely results from the placement of the KiK-net station locations on rock and rocklike sites, whereas stations in the other regions are generally placed in urban areas underlain by sediments. Using the KiK-net velocity profiles, we provide equations relating V S 30 to V Sz for z ranging from 5 to 29 m in 1-m increments. These equations (and those for California velocity profiles given in Boore, 2004b) can be used to estimate V S 30 from V Sz for sites in which velocity profiles do not extend to 30 m. The scatter of the residuals decreases with depth, but, even for an averaging depth of 5 m, a variation in log V S 30 of ±1 standard deviation maps into less than a 20% uncertainty in ground motions given by recent GMPEs at short periods. The sensitivity of the ground motions to V S 30 uncertainty is considerably larger at long periods (but is less than a factor of 1.2 for averaging depths greater than about 20 m). We also find that V S 30 is correlated with V Sz for z as great as 400 m for sites of the KiK-net network, providing some justification for using V S 30 as a site-response variable for predicting ground motions at periods for which the wavelengths far exceed 30 m.

Bulletin of the Seismological Society of America

Determination of elastic wave velocity and relative hypocenter locations using refracted waves. I. Methodology

An arrival time difference method utilizing refracted arrivals from earthquakes in a homogeneous, layered earth model has been developed for the simultaneous determination of near-source (in situ) velocity and relative locations of earthquakes. The method is particularly applicable when analyzing data from arrays in which most of the recording stations are far (i.e., several focal depths) from a group of events. This iterative scheme locates earthquakes relative to a master event and performs an inversion for in situ velocity using a generalized inverseleast squares estimation procedure. Direct arrivals, when available, may be included to stabilize the inversion and increase the accuracy of the event locations. We tested this scheme on artificial data contaminated by random and systematic arrival time errors, gaps in azimuthal coverage, and inaccuracies in the assumed velocity model. As usual, depth is the least well-resolved hypocenter coordinate, but this scheme yielded accurate locations of most events while converging to the correct velocity model.

Bulletin of the Seismological Society of America

Holocene slip rate and revised characteristic earthquake parameters for the western segment of the Castle Mountain fault, Alaska

The western segment of the Castle Mountain fault poses a significant seismic hazard to the most populated region of south-central Alaska. We identify a previously unrecognized margin of a postglacial outwash channel that is offset right laterally 36 ± 4 m across the western segment of the Castle Mountain fault. This offset occurred after glaciers withdrew from the lowland 11,300–15,380 cal yr b.p. and after outwash channel margins were cut and stabilized 11,210–13,470 cal yr b.p. Using these ages and the measured separation, we obtain a maximum slip rate of 3.0 ± 0.6 mm yr −1 and a minimum slip rate of 2.8 ± 0.7 mm yr −1 . These are the first lateral slip rates for the Castle Mountain fault established by a field measurement. Based on timing of the most recent earthquake, 670 ± 60 yr b.p. , the Castle Mountain fault could have accumulated an average single-event slip of about 1.9 m (extremes range from 1.3 to 2.6 m). The fault consists of two segments; a surface-rupturing earthquake likely will be limited to the 62-km-long western segment. Area-magnitude regression calculations suggest that such an earthquake on the western Castle Mountain fault would have a moment magnitude of 6.9 to 7.3.

Alaska

BSSA: Worth thinking about

The Bulletin of the Seismological Society of America ( BSSA ) is a powerful community project that has helped us share the information necessary to keep our field moving forward since 1911. In some ways, BSSA is much like it has always been, and each issue provides us with a collection of research that has been improved by the peer review process and copyedited, typeset, and printed to make it easily readable. But BSSA has also constantly evolved, and the online edition with full-text searching of all content back to volume 1 and linked references that speed our navigation through the related literature is unlike anything J. C. Branner could have imagined when he edited the first issue.

Seismological Research Letters

Recent advances in rotational seismology

Rotational seismology is an emerging field of study concerned with all aspects of rotational motions induced by earthquakes, explosions, and ambient vibrations. Two recent monographs (Teisseyre et al. 2006; Teisseyre et al. 2008) and a Bulletin of the Seismological Society of America special issue on Rotational Seismology and Engineering Applications (Lee, Celebi et al. 2009) are useful starting points. Rotational seismology is of interest to a wide range of disciplines, including various branches of seismology, earthquake engineering, and geodesy, as well as to physicists using Earth-based observatories for detecting gravitational waves generated by astronomical sources.

Seismological Research Letters

Topographic Slope as a Proxy for Seismic Site-Conditions (VS30) and Amplification Around the Globe

Executive Summary It is well-known that large global earthquakes can have a dramatic effect on local communities and the built environment. Moreover, ground motions amplified by surficial materials can exacerbate the situation, often making the difference between minor and major damage. For a real-time earthquake impact alert system, such as Prompt Assessment of Global Earthquakes for Response (PAGER) (Wald and others, 2006), we seek to rapidly evaluate potential ground shaking in the source region and subsequently provide an estimate of the population exposure to potentially fatal levels of ground shaking in any region of the world. The contribution of surficial geology (particularly soft sediments) to the amplification of ground shaking is an important component in predicting the levels of ground motion observed at any site. Unfortunately, the availability of information regarding seismic siteconditions is only available at a few sites around the globe. Herein, we describe a methodology for deriving maps of seismic site-conditions anywhere in the world using topographic slope as a proxy. Average shear-velocity down to 30 m (or VS30) measurements are correlated against topographic slope to develop two sets of coefficients for predicting VS30: one for active tectonic regions that possess dynamic topographic relief, and one for stable continental regions where changes in topography are more subdued. These coefficients have been applied to the continental United States, in addition to other regions around the world. They are subsequently compared to existing site-condition maps based on geology and observed VS30 measurements, where available. The application of the topographic slope method in regions with abundant VS30 measurements (for example California, Memphis, and Taiwan) indicates that this method provides site condition-maps of similar quality, or in some cases, maps superior to those developed from more traditional techniques. Having a first-order assessment of seismic site-conditions anywhere in the world provides a valuable tool to rapidly estimate ground motions following any global earthquake, the primary motivation for this research. These VS30 maps will enable us to better quantify possible ground shaking and rapidly deliver these predictions to emergency managers and responders. In addition, the VS30 maps for the globe will also have practical applications for numerous related probabilistic- and scenario-based studies. To date, several researchers have requested maps or have used the approach outlined herein for their own applications (for example Cagnan and Kariptas, written commun., 2007; Harmandar and others, 2007). Given that we anticipate a significant demand for these products, we have developed an internet delivery service so that users can download maps and grids of seismic site-conditions for specified regions. To some extent, these grids can also be customized by the user if they disagree with the predefined correlations derived using the methodologies described within this report. Finally, this report represents a more comprehensive account of this technique and provides a more fully illustrated global description of results than that given in Wald and Allen (2007), which has been accepted for publication in the Bulletin of the Seismological Society of America.

Open-File Report

Overview of the Southern San Andreas Fault Model

This appendix summarizes the data and methodology used to generate the source model for the southern San Andreas fault. It is organized into three sections, 1) a section by section review of the geological data in the format of past Working Groups, 2) an overview of the rupture model, and 3) a manuscript by Biasi and Weldon (in review Bulletin of the Seismological Society of America) that describes the correlation methodology that was used to help develop the ?geologic insight? model. The goal of the Biasi and Weldon methodology is to quantify the insight that went into developing all A faults; as such it is in concept consistent with all other A faults but applied in a more quantitative way. The most rapidly slipping fault and the only known source of M~8 earthquakes in southern California is the San Andreas fault. As such it plays a special role in the seismic hazard of California, and has received special attention in the current Working Group. The underlying philosophy of the current Working Group is to model the recurrence behavior of large, rapidly slipping faults like the San Andreas from observed data on the size, distribution and timing of past earthquakes with as few assumptions about underlying recurrence behavior as possible. In addition, we wish to carry the uncertainties in the data and the range of reasonable extrapolations from the data to the final model. To accomplish this for the Southern San Andreas fault we have developed an objective method to combine all of the observations of size, timing, and distribution of past earthquakes into a comprehensive set of earthquake scenarios that each represent a possible history of earthquakes for the past ~1400 years. The scenarios are then ranked according to their overall consistency with the data and then the frequencies of all of the ruptures permitted by the current Working Group?s segmentation model are calculated. We also present 30-yr conditional probabilities by segment and compare to previous results. A distinctive aspect of the current model is that the probability is higher at both ends of the fault and that the ends have a much greater fraction of smaller events. There is a significant difference in the likelihood of large (M 7.7-8.0) earthquakes along the fault from north to south, with large 1857-like events common on the northern half of the southern San Andreas fault but relatively few M 7.7-8.0 expected on the southern half.

Open-File Report

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)

Integration of paleoseismic data from multiple sites to develop an objective earthquake chronology: Application to the Weber segment of the Wasatch fault zone, Utah

We present a method to evaluate and integrate paleoseismic data from multiple sites into a single, objective measure of earthquake timing and recurrence on discrete segments of active faults. We apply this method to the Weber segment (WS) of the Wasatch fault zone using data from four fault-trench studies completed between 1981 and 2009. After systematically reevaluating the stratigraphic and chronologic data from each trench site, we constructed time-stratigraphic OxCal models that yield site probability density functions (PDFs) of the times of individual earthquakes. We next qualitatively correlated the site PDFs into a segment-wide earthquake chronology, which is supported by overlapping site PDFs, large per-event displacements, and prominent segment boundaries. For each segment-wide earthquake, we computed the product of the site PDF probabilities in common time bins, which emphasizes the overlap in the site earthquake times, and gives more weight to the narrowest, best-defined PDFs. The product method yields smaller earthquake-timing uncertainties compared to taking the mean of the site PDFs, but is best suited to earthquakes constrained by broad, overlapping site PDFs. We calculated segment-wide earthquake recurrence intervals and uncertainties using a Monte Carlo model. Five surface-faulting earthquakes occurred on the WS at about 5.9, 4.5, 3.1, 1.1, and 0.6 ka. With the exception of the 1.1-ka event, we used the product method to define the earthquake times. The revised WS chronology yields a mean recurrence interval of 1.3 kyr (0.7&ndash;1.9-kyr estimated two-sigma [2&delta;] range based on interevent recurrence). These data help clarify the paleoearthquake history of the WS, including the important question of the timing and rupture extent of the most recent earthquake, and are essential to the improvement of earthquake-probability assessments for the Wasatch Front region.

Utah

Ground motions from the 7 and 19 September, 2017 Tehuantepec and Puebla-Morelos, Mexico, earthquakes

The 2017 M "> M 8.2 Tehuantepec and M "> M 7.1 Puebla‐Morelos earthquakes were deep inslab normal‐faulting events that caused significant damage to several central‐to‐southern regions of Mexico. Inslab earthquakes are an important component of seismicity and seismic hazard in Mexico. Ground‐motion prediction equations (GMPEs) are an integral part of seismic hazard assessment as well as risk and rapid‐response products. This work examines the observed ground motions from these two events in comparison to the predicted median ground motions from four GMPEs. The residuals between the observed and modeled ground motions allow us to study regional differences in shaking, the effects of each earthquake, and basin effects in Mexico City, Puebla, and Oaxaca. We find that the ground motions from these two earthquakes are generally well modeled by the GMPEs. However, the Tehuantepec event shows larger than expected ground motions at greater distances and longer periods, which suggests a waveguide effect from the subduction zone geometry. Finally, Mexico City and the cities of Puebla and Oaxaca exhibit very large ground motions, indicative of well‐known site and basin effects that are much stronger than the basin terms included in some of the GMPEs. Simple and rapid ground‐motion parameter estimates that include site effects are key for hazard and real‐time risk assessments in regions such as Mexico, where the vast majority of the population lives in areas where the aforementioned effects are relevant. However, GMPEs based on site correction terms dependent on topographic slope proxies underestimate, at least in the three cities tackled in this work, the observed amplification. Therefore, there is a need to improve models of seismic amplification in basins that could be included in GMPEs.

Bulletin of the Seismological Society of America

Spectral inversion for seismic site response in central Oklahoma: Low-frequency resonances from the Great Unconformity

We investigate seismic site response by inverting seismic ground‐motion spectra for site and source spectral properties, in a region of central Oklahoma, where previous ground‐motion studies have indicated discrepancies between observations and ground‐motion models (GMMs). The inversion is constrained by a source spectral model, which we computed from regional seismic records, using aftershocks as empirical Green’s functions to deconvolve site and path effects. Site spectra across the region exhibit multiple, strong, low‐frequency ( ⁠ f <2 Hz ⁠ ) resonances. Modeling of vertically propagating SH waves reproduces the mean amplitudes and frequencies of the site spectra and requires a deep ( ⁠ ∼1–2 km ⁠ ) impedance contrast. Comparison of regional seismic velocity models and geologic profiles indicates that the seismic impedance contrast is, or is in proximity to, the Great Unconformity, which marks the interface between Precambrian basement rocks and overlying Paleozoic sedimentary rocks. Depth to Precambrian basement increases to the southwest across the study region ( ⁠ ∼1500–4500 m ⁠ ), and the fundamental frequencies of the site spectra are anticorrelated with basement depth. The first higher‐mode resonance also exhibits dependence on basement depth; although modeling suggests that the second higher mode should depend on basement depth, site spectra do not support this. The low‐frequency resonances in central Oklahoma are not represented in the GMMs used in current seismic hazard analyses for tectonic earthquakes, though approaches to account for such features are under consideration in other regions of the central and eastern United States. Given the broad spatial extent of the Great Unconformity underlying eastern North America, it is likely that similar effects on seismic site response also occur in other areas. This study highlights the impact of regional geologic structure on earthquake ground motions and reiterates the need for modeling regional effects to improve ground‐motion predictions and seismic hazard assessments.

Oklahoma