Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Seismological 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 37 records · Page 2Linked to original sources
USGS1990

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↗
USGS2011

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↗
USGS1985

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↗
USGS2004

Geophysical data reveal the crustal structure of the Alaska range orogen within the aftershock zone of the Mw 7.9 Denali fault earthquake

Geophysical information, including deep-crustal seismic reflection, magnetotelluric ( mt ), gravity, and magnetic data, cross the aftershock zone of the 3 November 2002 M w 7.9 Denali fault earthquake. These data and aftershock seismicity, jointly interpreted, reveal the crustal structure of the right-lateral-slip Denali fault and the eastern Alaska Range orogen, as well as the relationship between this structure and seismicity. North of the Denali fault, strong seismic reflections from within the Alaska Range orogen show features that dip as steeply as 25° north and extend downward to depths between 20 and 25 km. These reflections reveal crustal structures, probably ductile shear zones, that most likely formed during the Late Cretaceous, but these structures appear to be inactive, having produced little seismicity during the past 20 years. Furthermore, seismic reflections mainly dip north, whereas alignments in aftershock hypocenters dip south. The Denali fault is nonreflective, but modeling of mt , gravity, and magnetic data suggests that the Denali fault dips steeply to vertically. However, in an alternative structural model, the Denali fault is defined by one of the reflection bands that dips to the north and flattens into the middle crust of the Alaska Range orogen. Modeling of mt data indicates a rock body, having low electrical resistivity (>10 Ω·m), that lies mainly at depths greater than 10 km, directly beneath aftershocks of the Denali fault earthquake. The maximum depth of aftershocks along the Denali fault is 10 km. This shallow depth may arise from a higher-than-normal geothermal gradient. Alternatively, the low electrical resistivity of deep rocks along the Denali fault may be associated with fluids that have weakened the lower crust and helped determine the depth extent of the aftershock zone.

Alaska↗
USGS2007

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↗
USGS2011

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–1.9-kyr estimated two-sigma [2δ] 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↗
USGS2018

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

Empirical improvements for estimating earthquake response spectra with random‐vibration theory

The stochastic method of ground‐motion simulation is often used in combination with the random‐vibration theory to directly compute ground‐motion intensity measures, thereby bypassing the more computationally intensive time‐domain simulations. Key to the application of random‐vibration theory to simulate response spectra is determining the duration ( D rms ) used in computing the root‐mean‐square oscillator response. Boore and Joyner (1984) originally proposed an equation for D rms , which was improved upon by Liu and Pezeshk (1999) . Though these equations are both substantial improvements over using the duration of the ground‐motion excitation for D rms , we document systematic differences between the ground‐motion intensity measures derived from the random‐vibration and time‐domain methods for both of these D rms equations. These differences are generally less than 10% for most magnitudes, distances, and periods of engineering interest. Given the systematic nature of the differences, however, we feel that improved equations are warranted. We empirically derive new equations from time‐domain simulations for eastern and western North America seismological models. The new equations improve the random‐vibration simulations over a wide range of magnitudes, distances, and oscillator periods.

Bulletin of the Seismological Society of America↗
USGS2002

Magnitude estimates of two large aftershocks of the 16 December 1811 New Madrid earthquake

The three principal New Madrid mainshocks of 1811-1812 were followed by extensive aftershock sequences that included numerous felt events. Although no instrumental data are available for either the mainshocks or the aftershocks, available historical accounts do provide information that can be used to estimate magnitudes and locations for the large events. In this article we investigate two of the largest aftershocks: one near dawn following the first mainshock on 16 December 1811, and one near midday on 17 December 1811. We reinterpret original felt reports to obtain a set of 48 and 20 modified Mercalli intensity values of the two aftershocks, respectively. For the dawn aftershock, we infer a Mw of approximately 7.0 based on a comparison of its intensities with those of the smallest New Madrid mainshock. Based on a detailed account that appears to describe near-field ground motions, we further propose a new fault rupture scenario for the dawn aftershock. We suggest that the aftershock had a thrust mechanism and occurred on a southeastern limb of the Reelfoot fault. For the 17 December 1811 aftershock, we infer a M W of approximately 6.1 ± 0.2 . This value is determined using the method of Bakun et al. (2002), which is based on a new calibration of intensity versus distance for earthquakes in central and eastern North America. The location of this event is not well constrained, but the available accounts suggest an epicenter beyond the southern end of the New Madrid Seismic Zone.

Arkansas, Illinois, Missouri, Tennessee↗
USGS2001

Triggered earthquakes and the 1811-1812 New Madrid, central United States, earthquake sequence

The 1811-1812 New Madrid, central United States, earthquake sequence included at least three events with magnitudes estimated at well above M 7.0. I discuss evidence that the sequence also produced at least three substantial triggered events well outside the New Madrid Seismic Zone, most likely in the vicinity of Cincinnati, Ohio. The largest of these events is estimated to have a magnitude in the low to mid M 5 range. Events of this size are large enough to cause damage, especially in regions with low levels of preparedness. Remotely triggered earthquakes have been observed in tectonically active regions in recent years, but not previously in stable continental regions. The results of this study suggest, however, that potentially damaging triggered earthquakes may be common following large mainshocks in stable continental regions. Thus, in areas of low seismic activity such as central/ eastern North America, the hazard associated with localized source zones might be more far reaching than previously recognized. The results also provide additional evidence that intraplate crust is critically stressed, such that small stress changes are especially effective at triggering earthquakes.

Bulletin of the Seismological Society of America↗
USGS2010

Intensity, magnitude, location and attenuation in India for felt earthquakes since 1762

A comprehensive, consistently interpreted new catalog of felt intensities for India (Martin and Szeliga, 2010, this issue) includes intensities for 570 earthquakes; instrumental magnitudes and locations are available for 100 of these events. We use the intensity values for 29 of the instrumentally recorded events to develop new intensity versus attenuation relations for the Indian subcontinent and the Himalayan region. We then use these relations to determine the locations and magnitudes of 234 historical events, using the method of Bakun and Wentworth (1997). For the remaining 336 events, intensity distributions are too sparse to determine magnitude or location. We evaluate magnitude and location accuracy of newly located events by comparing the instrumental- with the intensity-derived location for 29 calibration events, for which more than 15 intensity observations are available. With few exceptions, most intensity-derived locations lie within a fault length of the instrumentally determined location. For events in which the azimuthal distribution of intensities is limited, we conclude that the formal error bounds from the regression of Bakun and Wentworth (1997) do not reflect the true uncertainties. We also find that the regression underestimates the uncertainties of the location and magnitude of the 1819 Allah Bund earthquake, for which a location has been inferred from mapped surface deformation. Comparing our inferred attenuation relations to those developed for other regions, we find that attenuation for Himalayan events is comparable to intensity attenuation in California (Bakun and Wentworth, 1997), while intensity attenuation for cratonic events is higher than intensity attenuation reported for central/eastern North America (Bakun et al. , 2003). Further, we present evidence that intensities of intraplate earthquakes have a nonlinear dependence on magnitude such that attenuation relations based largely on small-to-moderate earthquakes may significantly overestimate the magnitudes of historical earthquakes.

Bulletin of the Seismological Society of America↗
USGS2020

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↗
USGS1989

The dependence of peak horizontal acceleration on magnitude, distance, and site effects for small-magnitude earthquakes in California and eastern North America

One-hundred and ninety free-field accelerograms recorded on deep soil (> 10 m deep) were used to study the near-source scaling characteristics of peak horizontal acceleration for 91 earthquakes (2.5 ≦ ML ≦ 5.0) located primarily in California. Of the six attenuation relationships developed in this study, the one considered most reliable is given by the expression.

Bulletin of the Seismological Society of America↗
USGS2013

Geometry and earthquake potential of the shoreline fault, central California

The Shoreline fault is a vertical strike‐slip fault running along the coastline near San Luis Obispo, California. Much is unknown about the Shoreline fault, including its slip rate and the details of its geometry. Here, I study the geometry of the Shoreline fault at seismogenic depth, as well as the adjacent section of the offshore Hosgri fault, using seismicity relocations and earthquake focal mechanisms. The Optimal Anisotropic Dynamic Clustering (OADC) algorithm (Ouillon et al., 2008) is used to objectively identify the simplest planar fault geometry that fits all of the earthquakes to within their location uncertainty. The OADC results show that the Shoreline fault is a single continuous structure that connects to the Hosgri fault. Discontinuities smaller than about 1 km may be undetected, but would be too small to be barriers to earthquake rupture. The Hosgri fault dips steeply to the east, while the Shoreline fault is essentially vertical, so the Hosgri fault dips towards and under the Shoreline fault as the two faults approach their intersection. The focal mechanisms generally agree with pure right‐lateral strike‐slip on the OADC planes, but suggest a non‐planar Hosgri fault or another structure underlying the northern Shoreline fault. The Shoreline fault most likely transfers strike‐slip motion between the Hosgri fault and other faults of the Pacific–North America plate boundary system to the east. A hypothetical earthquake rupturing the entire known length of the Shoreline fault would have a moment magnitude of 6.4–6.8. A hypothetical earthquake rupturing the Shoreline fault and the section of the Hosgri fault north of the Hosgri–Shoreline junction would have a moment magnitude of 7.2–7.5.

California↗
USGS2015

A long-term earthquake rate model for the central and eastern United States from smoothed seismicity

I present a long-term earthquake rate model for the central and eastern United States from adaptive smoothed seismicity. By employing pseudoprospective likelihood testing (L-test), I examined the effects of fixed and adaptive smoothing methods and the effects of catalog duration and composition on the ability of the models to forecast the spatial distribution of recent earthquakes. To stabilize the adaptive smoothing method for regions of low seismicity, I introduced minor modifications to the way that the adaptive smoothing distances are calculated. Across all smoothed seismicity models, the use of adaptive smoothing and the use of earthquakes from the recent part of the catalog optimizes the likelihood for tests with M ≥2.7 and M ≥4.0 earthquake catalogs. The smoothed seismicity models optimized by likelihood testing with M ≥2.7 catalogs also produce the highest likelihood values for M ≥4.0 likelihood testing, thus substantiating the hypothesis that the locations of moderate-size earthquakes can be forecast by the locations of smaller earthquakes. The likelihood test does not, however, maximize the fraction of earthquakes that are better forecast than a seismicity rate model with uniform rates in all cells. In this regard, fixed smoothing models perform better than adaptive smoothing models. The preferred model of this study is the adaptive smoothed seismicity model, based on its ability to maximize the joint likelihood of predicting the locations of recent small-to-moderate-size earthquakes across eastern North America. The preferred rate model delineates 12 regions where the annual rate of M ≥5 earthquakes exceeds 2×10 −3 . Although these seismic regions have been previously recognized, the preferred forecasts are more spatially concentrated than the rates from fixed smoothed seismicity models, with rate increases of up to a factor of 10 near clusters of high seismic activity.

Bulletin of the Seismological Society of America↗
USGS2003

Empirical ground-motion relations for subduction-zone earthquakes and their application to Cascadia and other regions

Ground-motion relations for earthquakes that occur in subduction zones are an important input to seismic-hazard analyses in many parts of the world. In the Cascadia region (Washington, Oregon, northern California, and British Columbia), for example, there is a significant hazard from megathrust earthquakes along the subduction interface and from large events within the subducting slab. These hazards are in addition to the hazard from shallow earthquakes in the overlying crust. We have compiled a response spectra database from thousands of strong-motion recordings from events of moment magnitude ( M ) 5-8.3 occurring in subduction zones around the world, including both interface and in-slab events. The 2001 M 6.8 Nisqually and 1999 M 5.9 Satsop earthquakes are included in the database, as are many records from subduction zones in Japan (Kyoshin-Net data), Mexico (Guerrero data), and Central America. The size of the database is four times larger than that available for previous empirical regressions to determine ground-motion relations for subduction-zone earthquakes. The large dataset enables improved determination of attenuation parameters and magnitude scaling, for both interface and in-slab events. Soil response parameters are also better determined by the data.

British Columbia, California, Oregon, Washington↗
USGS2011

Monitoring the Earthquake source process in North America

With the implementation of the USGS National Earthquake Information Center Prompt Assessment of Global Earthquakes for Response system (PAGER), rapid determination of earthquake moment magnitude is essential, especially for earthquakes that are felt within the contiguous United States. We report an implementation of moment tensor processing for application to broad, seismically active areas of North America. This effort focuses on the selection of regional crustal velocity models, codification of data quality tests, and the development of procedures for rapid computation of the seismic moment tensor. We systematically apply these techniques to earthquakes with reported magnitude greater than 3.5 in continental North America that are not associated with a tectonic plate boundary. Using the 0.02-0.10 Hz passband, we can usually determine, with few exceptions, moment tensor solutions for earthquakes with M w as small as 3.7. The threshold is significantly influenced by the density of stations, the location of the earthquake relative to the seismic stations and, of course, the signal-to-noise ratio. With the existing permanent broadband stations in North America operated for rapid earthquake response, the seismic moment tensor of most earthquakes that are M w 4 or larger can be routinely computed. As expected the nonuniform spatial pattern of these solutions reflects the seismicity pattern. However, the orientation of the direction of maximum compressive stress and the predominant style of faulting is spatially coherent across large regions of the continent.

Bulletin of the Seismological Society of America↗
USGS2002

Seismic structure of the crust and uppermost mantle of North America and adjacent oceanic basins: A synthesis

We present a new set of contour maps of the seismic structure of North America and the surrounding ocean basins. These maps include the crustal thickness, whole-crustal average P -wave and S -wave velocity, and seismic velocity of the uppermost mantle, that is, Pn and Sn . We found the following: (1) The average thickness of the crust under North America is 36.7 km (standard deviation [s.d.] ±8.4 km), which is 2.5 km thinner than the world average of 39.2 km (s.d. ± 8.5) for continental crust; (2) Histograms of whole-crustal P - and S -wave velocities for the North American crust are bimodal, with the lower peak occurring for crust without a high-velocity (6.9–7.3 km/sec) lower crustal layer; (3) Regions with anomalously high average crustal P -wave velocities correlate with Precambrian and Paleozoic orogens; low average crustal velocities are correlated with modern extensional regimes; (4) The average Pn velocity beneath North America is 8.03 km/sec (s.d. ± 0.19 km/sec); (5) the well-known thin crust beneath the western United States extends into north-west Canada; (6) the average P -wave velocity of layer 3 of oceanic crust is 6.61 km/sec (s.d. ± 0.47 km/sec). However, the average crustal P -wave velocity under the eastern Pacific seafloor is higher than the western Atlantic seafloor due to the thicker sediment layer on the older Atlantic seafloor.

Bulletin of the Seismological Society of America↗