Search USGS⌕ Search

SEARCH · Search USGS

Results for “Bulletin of the Geological Society of 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 829 records · Page 46Linked to original sources

Effects of local geological conditions in the San Francisco Bay region on ground motions and the intensities of the 1906 earthquake

Measurements of ground motion generated by nuclear explosions in Nevada have been completed for 99 locations in the San Francisco Bay region, California. The recordings show marked amplitude variations in the frequency band 0.25 to 3.0 Hz that are consistently related to the local geological conditions of the recording site. The average spectral amplifications observed for vertical and horizontal ground motions are, respectively: (1, 1) for granite, (1.5, 1.6) for the Franciscan Formation, (3.0, 2.7) for the Santa Clara Formation, (3.3, 4.4) for alluvium, and (3.7, 11.3) for bay mud. Spectral amplification curves define predominant ground frequencies in the band 0.25 to 3.0 E for bay mud sites and for some alluvial sites. Amplitude spectra computed from recordings of seismic background noise at 50 sites do not generally define predominant ground frequencies. The intensities ascribed to various sites in the San Francisco Bay region for the California earthquake of April 18, 1906, are strongly dependent on distance from the zone of surface faulting and the geological character of the ground. Considering only those sites (approximately one square city block in size) for which there is good evidence for the degree of ascribed intensity, the intensities for 917 sites on Franciscan rocks generally decrease with the logarithm of distance as Intensity = 2.69 - 1.90 log (Distance in kilometers). (1) For sites on other geological units, intensity increments, derived from this empirical relation, correlate strongly with the Average Horizontal Spectral Amplifications (AHSA) according to the empirical relation Intensity Increment = 0.27 + 2.70 log (AHSA). (2) Average intensity increments predicted for the various geological units are −0.3 for granite, 0.2 for the Franciscan Formation, 0.6 for the Great Valley sequence, 0.8 for the Santa Clara Formation, 1.3 for alluvium, and 2.4 for bay mud. The maximum intensity map predicted on the basis of these data delineates areas in the San Francisco Bay region of potentially high intensity for large earthquakes on either the San Andreas fault or the Hayward fault. The map provides a crude form of seismic zonation for the region and may be useful for certain general types of land-use zonation.

California↗

Probability of rupture of multiple fault segments

Fault segments identified from geologic and historic evidence have sometimes been adopted as features limiting the likely extends of earthquake ruptures. There is no doubt that individual segments can sometimes join together to produce larger earthquakes. This work is a trial of an objective method to determine the probability of multisegment ruptures. The frequency of occurrence of events on all conjectured combinations of adjacent segments in northern California is found by fitting to both geologic slip rates and to an assumed distribution of event sizes for the region as a whole. Uncertainty in the shape of the distribution near the maximum magnitude has a large effect on the solution. Frequencies of individual events cannot be determined, but it is possible to find a set of frequencies to fit a model closely. A robust conclusion for the San Francisco Bay region is that large multisegment events occur on the San Andreas and San Gregorio faults, but single-segment events predominate on the extended Hayward and Calaveras strands of segments.

Bulletin of the Seismological Society of America↗

Ground motions recorded in Rome during the April 2009 L’Aquila seismic sequence: site response and comparison with ground‐motion predictions based on a global dataset

The mainshock and moderate‐magnitude aftershocks of the 6 April 2009 M 6.3 L’Aquila seismic sequence, about 90 km northeast of Rome, provided the first earthquake ground‐motion recordings in the urban area of Rome. Before those recordings were obtained, the assessments of the seismic hazard in Rome were based on intensity observations and theoretical considerations. The L’Aquila recordings offer an unprecedented opportunity to calibrate the city response to central Apennine earthquakes—earthquakes that have been responsible for the largest damage to Rome in historical times. Using the data recorded in Rome in April 2009, we show that (1) published theoretical predictions of a 1 s resonance in the Tiber valley are confirmed by observations showing a significant amplitude increase in response spectra at that period, (2) the empirical soil‐transfer functions inferred from spectral ratios are satisfactorily fit through 1D models using the available geological, geophysical, and laboratory data, but local variability can be large for individual events, (3) response spectra for the motions recorded in Rome from the L’Aquila earthquakes are significantly amplified in the radial component at periods near 1 s, even at a firm site on volcanic rocks, and (4) short‐period response spectra are smaller than expected when compared to ground‐motion predictions from equations based on a global dataset, whereas the observed response spectra are higher than expected for periods near 1 s.

Rome↗

Simulation of earthquake ground motions in the eastern United States using deterministic physics‐based and site‐based stochastic approaches

Earthquake ground‐motion recordings are scarce in the central and eastern United States (CEUS) for large‐magnitude events and at close distances. We use two different simulation approaches, a deterministic physics‐based method and a site‐based stochastic method, to simulate ground motions over a wide range of magnitudes. Drawing on previous results for the modeling of recordings from the 2011 M w 5.8 Mineral, Virginia, earthquake and using the 2001 M w 7.6 Bhuj, India, earthquake as a tectonic analog for a large magnitude CEUS event, we are able to calibrate the two simulation methods over this magnitude range. Both models show a good fit to the Mineral and Bhuj observations from 0.1 to 10 Hz. Model parameters are then adjusted to obtain simulations for M w 6.5, 7.0, and 7.6 events in the CEUS. Our simulations are compared with the 2014 U.S. Geological Survey weighted combination of existing ground‐motion prediction equations in the CEUS. The physics‐based simulations show comparable response spectral amplitudes and a fairly similar attenuation with distance. The site‐based stochastic simulations suggest a slightly faster attenuation of the response spectral amplitudes with distance for larger magnitude events and, as a result, slightly lower amplitudes at distances greater than 200 km. Both models are plausible alternatives and, given the few available data points in the CEUS, can be used to represent the epistemic uncertainty in modeling of postulated CEUS large‐magnitude events.

Bulletin of the Seismological Society of America↗

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States↗

Liquefaction and related ground failure from July 2019 Ridgecrest earthquake sequence

The 2019 Ridgecrest earthquake sequence produced a 4 July M 6.5 foreshock and a 5 July M 7.1 mainshock, along with 23 events with magnitudes greater than 4.5 in the 24 hr period following the mainshock. The epicenters of the two principal events were located in the Indian Wells Valley, northwest of Searles Valley near the towns of Ridgecrest, Trona, and Argus. We describe observed liquefaction manifestations including sand boils, fissures, and lateral spreading features, as well as proximate non‐ground failure zones that resulted from the sequence. Expanding upon results initially presented in a report of the Geotechnical Extreme Events Reconnaissance Association, we synthesize results of field mapping, aerial imagery, and inferences of ground deformations from Synthetic Aperture Radar‐based damage proxy maps (DPMs). We document incidents of liquefaction, settlement, and lateral spreading in the Naval Air Weapons Station China Lake US military base and compare locations of these observations to pre‐ and postevent mapping of liquefaction hazards. We describe liquefaction and ground‐failure features in Trona and Argus, which produced lateral deformations and impacts on several single‐story masonry and wood frame buildings. Detailed maps showing zones with and without ground failure are provided for these towns, along with mapped ground deformations along transects. Finally, we describe incidents of massive liquefaction with related ground failures and proximate areas of similar geologic origin without ground failure in the Searles Lakebed. Observations in this region are consistent with surface change predicted by the DPM. In the same region, geospatial liquefaction hazard maps are effective at identifying broad percentages of land with liquefaction‐related damage. We anticipate that data presented in this article will be useful for future liquefaction susceptibility, triggering, and consequence studies being undertaken as part of the Next Generation Liquefaction project.

California↗

Reflected seismic waves and their effect on strong shaking during the 1989 Loma Prieta, California, earthquake

Our data indicate that critical and postcritical reflections from crustal layers and the Moho produced increased shaking at discrete distances along the San Francisco Peninsula during the 1989 Loma Prieta earthquake. These reflections may have produced an increase in amplitude that is as much as 10 times greater than that of the direct arrival. Peak amplitude-distance patterns measured from explosive sources, synthetic seismograms, aftershocks, and the mainshock of the 1989 Loma Prieta earthquake indicate that (1) point sources, such as explosions, produce similar peak amplitude-distance relations as distributed, double-couple sources when the crustal structure is approximately a Poisson solid; (2) peak amplitudes from smaller point sources (explosions) may be scaled to those of larger sources (earthquakes); and (3) reflections caused a pattern of high amplitudes at specific distances along the San Francisco Peninsula that geographically correlates with areas of reported damage following from the Loma Prieta mainshock. Our study indicates that critical and postcritical reflections were stronger influences on the locations of strong shaking than local geology because a number of sites around the San Francisco Bay that are underlain by unconsolidated sediments experienced much less shaking than sites that were underlain by hard rock. Furthermore, some severely shaken hard-rock sites were farther from the epicenter than many of the less severely shaken “soft-sediment” sites. Models of the regional crustal structure and seismic-wave propagation paths may be useful in identifying sites around the San Francisco Bay region that are at risk due to reflected seismic energy, and this type of analysis may be useful in other seismically active regions.

California↗

Thin‐ or thick‐skinned faulting in the Yakima fold and thrust belt (WA)? Constraints from kinematic modeling of the saddle mountains anticline

The Yakima fold and thrust belt (YFTB) deforms the Columbia River Basalt Group flows of Washington State. The YFTB fault geometries and slip rates are crucial parameters for seismic‐hazard assessments of nearby dams and nuclear facilities, yet there are competing models for the subsurface fault geometry involving shallowly rooted versus deeply rooted fault systems. The YFTB is also thought to be analogous to the evenly spaced wrinkle ridges found on other terrestrial planets. Using seismic reflection data, borehole logs, and surface geologic data, we tested two proposed kinematic end‐member thick‐ and thin‐skinned fault models beneath the Saddle Mountains anticline of the YFTB. Observed subsurface geometry can be produced by 600–800 m of heave along a single listric‐reverse fault or ∼3.5 km of slip along two superposed low‐angle thrust faults. Both models require decollement slip between 7 and 9 km depth, resulting in greater fault areas than sometimes assumed in hazard assessments. Both models require initial slip much earlier than previously thought and may provide insight into the subsurface geometry of analogous comparisons to wrinkle ridges observed on other planets.

Washington↗

Status and performance of the ShakeAlert® earthquake early warning system: 2019-2023

The U.S. Geological Survey (USGS)‐operated ShakeAlert® system is the United States West Coast earthquake early warning system ( Given et al. , 2018 ). In this study we detail ShakeAlert’s performance during some of the largest events seen by the system thus far. Statewide public alerting using ShakeAlert messages was authorized in California in October 2019. Over the next few years, public alerts were expanded into Oregon and then into Washington ( U.S. Geological Survey, 2024 ). ShakeAlert source results are routinely compared to the USGS Comprehensive Catalog (ComCat; Guy et al. , 2015 ; U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which contains the earthquake location and magnitude determined using complete waveform data. M 4.5 and larger is the threshold used for public alerting and was deliberately set below the level where damage is likely to compensate for cases where the system underestimates the magnitude. Between 17 October 2019 and 1 September 2023, the ShakeAlert system created 95 events with maximum magnitude estimates of M ≥4.5, the public alerting threshold. 94 of the 95 events were due to real earthquakes. Seven were categorized “false” per ShakeAlert’s internal definition that there was no matching catalog event within 100 km and 30 s of origin time; however, all but one of these were real earthquakes that were poorly located, primarily because they were at the edges of the seismic network. Three detected events were labeled “missed” because they were very poorly located (>100 km location error). In addition, the system did not produce solutions for four ComCat events M ≥4.5 ( U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which were all at the edge of the alerting and network boundaries. The ShakeAlert system has accurately detected the majority of earthquakes that have occurred within the operational region since completing the public rollout, and alerts from the system have been delivered to millions of cell phone users throughout the West Coast.

California, Oregon, Washington↗

Upper-crustal structure of the inner Continental Borderland near Long Beach, California

A new P-wave velocity/structural model for the inner Continental Borderland (ICB) region was developed for the area near Long Beach, California. It combines controlled-source seismic reflection and refraction data collected during the 1994 Los Angeles Region Seismic Experiment (LARSE), multichannel seismic reflection data collected by the U.S. Geological Survey (1998-2000), and nearshore borehole stratigraphy. Based on lateral velocity contrasts and stratigraphic variation determined from borehole data, we are able to locate major faults such as the Cabrillo, Palos Verdes, THUMS-Huntington Beach, and Newport Inglewood fault zones, along with minor faults such as the slope fault, Avalon knoll, and several other yet unnamed faults. Catalog seismicity (1975-2002) plotted on our preferred velocity/structural model shows recent seismicity is located on 16 out of our 24 faults, providing evidence for continuing concern with respect to the existing seismic-hazard estimates. Forward modeling of P-wave arrival times on the LARSE line 1 resulted in a four-layer model that better resolves the stratigraphy and geologic structures of the ICB and also provides tighter constraints on the upper-crustal velocity structure than previous modeling of the LARSE data. There is a correlation between the structural horizons identified in the reflection data with the velocity interfaces determined from forward modeling of refraction data. The strongest correlation is between the base of velocity layer 1 of the refraction model and the base of the planar sediment beneath the shelf and slope determined by the reflection model. Layers 2 and 3 of the velocity model loosely correlate with the diffractive crust layer, locally interpreted as Catalina Schist.

Bulletin of the Seismological Society of America↗

A comparative ground response study near Los Angeles using recordings of Nevada nuclear tests and the 1971 San Fernando earthquake

A comparative ground response study at sites in the Los Angeles region is based on the extensive strong-motion data set recorded in the 1971 San Fernando earthquake and 159 three-component recordings of distant nuclear tests in Nevada. Amplitude spectral ratios computed for the nuclear test data over those frequency bands for which there is an adequate signal-to-noise ratio provide statistically stable estimates of the local ground response related to the type of local geologic conditions. Comparison of the strong-motion data recorded from the 1971 San Fernando earthquake with that recorded at 28 of the same sites for the nuclear tests show that the corresponding amplitude spectral ratios are statistically equivalent for most sites, provided reference stations are chosen to minimize effects on the earthquake data of azimuthal source variations and crustal propagation path. Statistical equivalence of the spectral ratios for the two types of data sources suggests that amplitude spectral ratios computed with respect to the appropriate reference station provide a first-order estimate of local ground response.

California↗

A terrain-based site characterization map of California with implications for the contiguous United States

We present an approach based on geomorphometry to predict material properties and characterize site conditions using the V S 30 parameter (time‐averaged shear‐wave velocity to a depth of 30 m). Our framework consists of an automated terrain classification scheme based on taxonomic criteria (slope gradient, local convexity, and surface texture) that systematically identifies 16 terrain types from 1‐km spatial resolution (30 arcsec) Shuttle Radar Topography Mission digital elevation models (SRTMDEMs). Using 853 V S 30 values from California, we apply a simulation‐based statistical method to determine the mean V S 30 for each terrain type in California. We then compare the V S 30 values with models based on individual proxies, such as mapped surface geology and topographic slope, and show that our systematic terrain‐based approach consistently performs better than semiempirical estimates based on individual proxies. To further evaluate our model, we apply our California‐based estimates to terrains of the contiguous United States. Comparisons of our estimates with 325 V S 30 measurements outside of California, as well as estimates based on the topographic slope model, indicate our method to be statistically robust and more accurate. Our approach thus provides an objective and robust method for extending estimates of V S 30 for regions where in situ measurements are sparse or not readily available.

California↗

Effects of thermal variability on broadband seismometers: Controlled experiments, observations, and implications

Isolating seismic instruments from temperature fluctuations is routine practice within the seismological community. However, the necessary degree of thermal stability required in broadband installations to avoid generating noise or compromising the fidelity in the seismic records is largely unknown and likely application dependent. To quantify the temperature sensitivity of seismometers over a broad range of frequencies, we artificially induced local temperature changes on three different models of seismometers to measure the effect of thermal variations on seismometer output. We found that diurnal temperature changes above 0.002°C root mean square (rms) showed significant changes in velocity and acceleration output in comparison to thermally stable reference measurements. We also found that sensor incoherent self‐noise increased with temperature variation; these increases in noise can be modeled as 1 / f "> 1 / f noise (pink noise), and are unlikely to be easily corrected for. These experimental results are compared with the data from Incorporated Research Institutions for Seismology (IRIS) U.S. Geological Survey (USGS) Global Seismographic Network (GSN) station TUC (Tucson, Arizona). This station is well instrumented with temperature sensors and has three different broadband seismometers, each of which uses a different method of thermal isolation. We show that the water bricks and borehole installations give ample temperature attenuation to thermally isolate seismometers from diurnal thermal variability that would compromise seismic data. We find that seismometer installations that provide thermal stability below 0.002°C rms could help to improve long‐period vertical seismic data across the GSN by decreasing temperature‐driven 1 / f "> 1 / f noise.

Bulletin of the Seismological Society of America↗

Effects of local geology on ground motion near San Francisco Bay

Measurements of ground motion generated by nuclear explosions in Nevada were made for 37 locations near San Francisco Bay, California. The results were compared with the San Francisco 1906 earthquake intensities and the strong-motion recordings of the San Francisco earthquake of March 22, 1957. The recordings show marked amplitude variations which are related consistently to the geologic setting of the recording site. For sites underlain by a layer of younger bay mud or artificial fill, maximum horizontal ground velocities generally increased with thickness of the layer and were as much as ten times greater than those recorded on nearby bedrock. The maximum vertical velocities for these sites were between 1 and 3.5 times greater. Spectral amplification curves clearly define a “dominant ground period” of about 1 second for sites underlain by younger bay mud. For sites underlain by older, more consolidated sediments, no clearly defined “dominant ground period” was found. Maximum ground velocities for the older bay sediment sites were about twice those recorded on bedrock.

California↗

Long Valley caldera and the UCERF depiction of Sierra Nevada range-front faults

Long Valley caldera lies within a left-stepping offset in the north-northwest-striking Sierra Nevada range-front normal faults with the Hilton Creek fault to the south and Hartley Springs fault to the north. Both Uniform California Earthquake Rupture Forecast (UCERF) 2 and its update, UCERF3, depict slip on these major range-front normal faults as extending well into the caldera, with significant normal slip on overlapping, subparallel segments separated by ∼10  km. This depiction is countered by (1) geologic evidence that normal faulting within the caldera consists of a series of graben structures associated with postcaldera magmatism (intrusion and tumescence) and not systematic down-to-the-east displacements consistent with distributed range-front faulting and (2) the lack of kinematic evidence for an evolving, postcaldera relay ramp structure between overlapping strands of the two range-front normal faults. The modifications to the UCERF depiction described here reduce the predicted shaking intensity within the caldera, and they are in accord with the tectonic influence that underlapped offset range-front faults have on seismicity patterns within the caldera associated with ongoing volcanic unrest.

Bulletin of the Seismological Society of America↗

Seismic velocity model of the central United States (Version 1): Description and simulation of the 18 April 2008 Mt. Carmel, Illinois, Earthquake

We have developed a new three‐dimensional seismic velocity model of the central United States (CUSVM) that includes the New Madrid Seismic Zone (NMSZ) and covers parts of Arkansas, Mississippi, Alabama, Illinois, Missouri, Kentucky, and Tennessee. The model represents a compilation of decades of crustal research consisting of seismic, aeromagnetic, and gravity profiles; geologic mapping; geophysical and geological borehole logs; and inversions of the regional seismic properties. The density, P ‐ and S ‐wave velocities are synthesized in a stand‐alone spatial database that can be queried to generate the required input for numerical seismic‐wave propagation simulations. We test and calibrate the CUSVM by simulating ground motions of the 18 April 2008 M w 5.4 Mt. Carmel, Illinois, earthquake and comparing the results with observed records within the model area. The selected stations in the comparisons reflect different geological site conditions and cover distances ranging from 10 to 430 km from the epicenter. The results, based on a qualitative and quantitative goodness‐of‐fit (GOF) characterization, indicate that both within and outside the Mississippi Embayment the CUSVM reasonably reproduces: (1) the body and surface‐wave arrival times and (2) the observed regional variations in ground‐motion amplitude, cumulative energy, duration, and frequency content up to a frequency of 1.0 Hz. In addition, we discuss the probable structural causes for the ground‐motion patterns in the central United States that we observed in the recorded motions of the 18 April Mt. Carmel earthquake.

Illinois↗

A V S30 map for California with geologic and topographic constraints

For many earthquake engineering applications, site response is estimated through empirical correlations with the time‐averaged shear‐wave velocity to 30 m depth ( V S 30 ). These applications therefore depend on the availability of either site‐specific V S 30 measurements or V S 30 maps at local, regional, and global scales. Because V S 30 measurements are sparse, a proxy frequently is needed to estimate V S 30 at unsampled locations. We present a new V S 30 map for California, which accounts for observational constraints from multiple sources and spatial scales, such as geology, topography, and site‐specific V S 30 measurements. We apply the geostatistical approach of regression kriging (RK) to combine these constraints for predicting V S 30 . For the V S 30 trend, we start with geology‐based V S 30 values and identify two distinct trends between topographic gradient and the residuals from the geology V S 30 model. One trend applies to deep and fine Quaternary alluvium, whereas the second trend is slightly stronger and applies to Pleistocene sedimentary units. The RK framework ensures that the resulting map of California is locally refined to reflect the rapidly expanding database of V S 30 measurements throughout California. We compare the accuracy of the new mapping method to a previously developed map of V S 30 for California. We also illustrate the sensitivity of ground motions to the new V S 30 map by comparing real and scenario ShakeMaps with V S 30 values from our new map to those for existing V S 30 maps.

California↗

Tectonics and seismicity of the southern Washington Cascade range

Geophysical, geological, and seismicity data are combined to develop a transpressional strain model for the southern Washington Cascades region. We use this model to explain oblique fold and fault systems, transverse faults, and a linear seismic zone just west of Mt. Rainier known as the western Rainier zone. We also attempt to explain a concentration of earthquakes that connects the northwest-trending Mount St. Helens seismic zone to the north-trending western Rainier zone. Our tectonic model illustrates the pervasive effects of accretionary processes, combined with subsequent transpressive forces generated by oblique subduction, on Eocene to present crustal processes, such as seismicity and volcanism.

Washington↗