USGS Science⌕ Search

SEARCH · USGS Science

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

ETAS-positive: Extracting unbiased epidemic-type aftershock models from incomplete earthquake catalogs

Earthquake catalogs are inherently incomplete, and short-term fluctuations in catalog completeness can strongly bias estimates of seismicity parameters. Here, we extend the positive statistics framework recently introduced for estimating Gutenberg–Richter a and b values to the Epidemic-Type Aftershock Sequence (ETAS) model. The approach focuses on complete magnitude–time intervals defined by pairs of earthquakes in which the second event is larger than the first. Comprehensive synthetic tests demonstrate that the proposed ETAS + method retrieves nearly unbiased parameter estimates even from highly incomplete catalogs, whereas parameters obtained using classical estimation approaches are strongly biased. Application to Southern California reveals a similar trend, with the productivity scaling parameter α approaching the b value, consistent with self-similar earthquake triggering within estimation uncertainties. Owing to its robustness against catalog incompleteness and computational efficiency, the ETAS + method provides a powerful new tool for seismicity analysis and earthquake forecasting.

California↗

DevelNet: Earthquake detection on Develocorder films with deep learning: Application to the rangely earthquake control experiment

There exists over a century of instrumental seismic data; however, most seismograms recorded before the 1980s are only available in analog form. Although analog seismograms are of great value, they are underutilized due to the difficulties of making quantitative measurements on the original media and in converting them to digital time series. In this study, we present an alternative workflow, based on deep learning, to reconstruct an earthquake catalog from images of analog data without conversion to vector time series. We trained a convolutional neural network—DevelNet, using synthetic analog data to detect earthquakes on scanned multichannel Develocorder film images. We then developed an image‐based processing workflow to measure arrival times, locate, and determine the magnitudes of earthquakes in the data. We demonstrate the performance of this approach on two years of continuous Develocorder film recordings from the Rangely earthquake control experiment in the mid‐1970s. Our approach detects twice the number of events reported in the original catalog ( Raleigh et al. , 1976 ). This demonstrates that DevelNet efficiently detects earthquakes from Develocorder film scans, performs consistently over time, and is robust to changes in network geometry. Our locations generally agree with the original study, although the automatically measured arrival times are less precise than manual reading, leading to increased location scatter. Our automatic workflow of Develocorder films rivals the performance of skilled analysts in earthquake detection, but with minimal human intervention. This image‐based processing offers a new approach for effectively and efficiently extracting earthquake information from analog seismic data.

Bulletin of the Seismological Society of America↗

Earthquake of August 19, 1966, Varto Area, eastern Turkey

The earthquake of August 19, 1966, centered near Varto in eastern Turkey. It occurred at 12:22:09.6 G. M. T. and had a Richter magnitude of about 7. Intensities reached IX. Reportedly, 2,529 people were killed, 1,500 injured, and 19,013 buildings were demolished or heavily damaged. Surface fractures in three parallel fracture zones trending N65-70°W display both right-lateral strike slip and extension normal to the fractures. The trend and sense of displacement relate the fracturing to the North Anatolia fault system. Examples of geologic factors influencing damage include: (1) degree of saturation of the alluvium, (2) proximity to landslides, and (3) relation to surface faulting.

Near Varto in eastern Turkey↗

A method for estimating the uncertainty of seismic velocities measured by refraction techniques

Time residuals from 7 5-km segments of 18 crustal seismic-refraction profiles in the Basin and Range province are used to investigate the validity of the linear-regression model and to make large sample estimates of the variance in the travel time distributions. A formula for unbiased estimates of velocity uncertainty is derived, assuming a linear trend with distance for the variances of the travel-time distributions. If the recording units are symmetric about the center of the recording interval, this formula is equivalent to the one derived assuming the variances are equal. At the 95-per cent confidence level the chi-squared test implied 84 per cent of the time-residual samples were inconsistent with the hypothesis that their parent populations had Gaussian distributions. If the number of recording locations expceeds 8, confidence limits computed without the Gaussian assumption suggest the departures from normality are not significant for velocity uncertainty estimates. The large sample estimates of the time-residual populations may be applicable to other areas. This evidence motivated the development of a method, requiring very little numerical calculation, for estimating uncertainties in velocities. The method requires, in addition to the large sample estimates of the travel time variances, information on the quality of the data, the location of the recording interval, and the number of recording units. The method is useful for the design of new experiments and independent estimates of uncertainty reported in the literature.

Arizona, California, Idaho, Nevada, Oregon, Utah↗

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↗

Geological aspects of the May 31, 1970, Perú earthquake

Geological effects of the destructive May 31, 1970, Perú earthquake (Ms = 7.7) extended over roughly 65,000 km2 of west-central Perú. Earthquake-triggered slope failures of all types that occurred throughout the mountainous parts of the region extensively damaged transportation routes and irrigation canals and temporarily dammed some rivers and lakes. The geologically most important and spectacular of these, a cataclysmic debris avalanche from between 5,500 and 6,400 m altitude on the north peak of Huascarán Mountain, buried the city of Yungay and part of Ranrahirca (elevation about 2,500 m) with a loss of more than about 18,000 lives. The debris avalanche involved 50-100 million m3 of rock, ice, snow, and soil that traveled 14.5 km from the source to Yungay at an average velocity of between 280 and 335 km/hr. Exceptionally rapid movement of the avalanche is indicated by eyewitness accounts, by topographic irregularities as high as 140 m that were overridden, and locally by boulders weighing several tons that were hurled as much as 1,000 m beyond its margins. A pulse of muddy water from the debris avalanche that swept down the Río Santa 160 km to the sea inundated farms and small settlements, buried highway and railroad routes, and destroyed the diversion dam and access bridge to a major hydroelectric plant.

Bulletin of the Seismological Society of America↗

The earthquake sequence near Danville, California, 1970

Several thousand small earthquakes, with magnitudes ranging up to 4, occurred near Danville, California, during May, June, and July 1970. Seven temporary seismographs were installed near the epicentral region to augment an existing telemetered network within 1 day after the first felt earthquake. The dense concentration of 10 seismographs within 12 km of the seismic activity and the existence of a reversed seismic-refraction profile through the region permitted a very detailed study of this earthquake sequence. Over 400 events with magnitude greater than 0.5 were located with an average accuracy of ±1 km. The main earthquake region was approximately spherical: 1 km in radius, 6 km deep, and centered about 5 km southeast of Danville, 40 km east of San Francisco. Fault-plane solutions indicate that faulting associated with the earthquakes was of the right-lateral strike-slip type, with a strike of N 35°W.

California↗

Effects of the Bear Valley and San Juan Bautista earthquakes of 1972 on Geodimeter line lengths

Measurements of lines of the California Geodimeter network in the vicinity of the epicenters of four 1972 earthquakes (magnitudes 4.7, 4.7, 4.8, and 5.1) along the San Andreas fault system indicate that no significant anomalous changes in line length preceded or accompanied those earthquakes. Within the precision of measurement, the data are consistent with a linear change in line length with time. Measurements of two 20-km-long lines made 2 days before and 1 day after a magnitude 4.7 earthquake showed no significant change even though the epicenter of the earthquake was only about 15 km from each of the three stations at the end points of the lines.

California↗

Geodetic determination of strain at the Nevada Test Site following the Handley event

Repeated surveys of a trilateration network (aperture greater than 20 km) centered on ground zero for the HANDLEY event, a nuclear explosion at the Nevada Test Site with yield in excess of 1 megaton, suggest that the explosion induced an east-west extension of the network by more than 50 mm. In the year following the detonation, this deformation reversed such that the final configuration represented a small east-west contraction from the pre-HANDLEY state. In the subsequent 2-year period, only minor deformation was detected. Thus, the overall response of Pahute Mesa may be described as stable. The explosion-induced deformation is thought to be partly due to slip on faults driven by the large-amplitude seismic waves from the explosion. The mechanism of the postshot relaxation is not understood.

Nevada↗

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↗

Peak horizontal acceleration and velocity from strong-motion records including records from the 1979 imperial valley, California, earthquake

We have taken advantage of the recent increase in strong-motion data at close distances to derive new attenuation relations for peak horizontal acceleration and velocity. This new analysis uses a magnitude-independent shape, based on geometrical spreading and anelastic attenuation, for the attenuation curve. An innovation in technique is introduced that decouples the determination of the distance dependence of the data from the magnitude dependence.

California↗

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↗

Statistical relations among earthquake magnitude, surface rupture length, and surface fault displacement

In order to refine correlations of surface-wave magnitude, fault rupture length at the ground surface, and fault displacement at the surface by including the uncertainties in these variables, the existing data were critically reviewed and a new data base was compiled. Earthquake magnitudes were redetermined as necessary to make them as consistent as possible with the Gutenberg methods and results, which make up much of the data base. Measurement errors were estimated for the three variables for 58 moderate to large shallow-focus earthquakes. Regression analyses were then made utilizing the estimated measurement errors. The regression analysis demonstrates that the relations among the variables magnitude, length, and displacement are stochastic in nature. The stochastic variance, introduced in part by incomplete surface expression of seismogenic faulting, variation in shear modulus, and regional factors, dominates the estimated measurement errors. Thus, it is appropriate to use ordinary least squares for the regression models, rather than regression models based upon an underlying deterministic relation in which the variance results primarily from measurement errors. Significant differences exist in correlations of certain combinations of length, displacement, and magnitude when events are grouped by fault type or by region, including attenuation regions delineated by Evernden and others. Estimates of the magnitude and the standard deviation of the magnitude of a prehistoric or future earthquake associated with a fault can be made by correlating M s with the logarithms of rupture length, fault displacement, or the product of length and displacement. Fault rupture area could be reliably estimated for about 20 of the events in the data set. Regression of M s on rupture area did not result in a marked improvement over regressions that did not involve rupture area. Because no subduction-zone earthquakes are included in this study, the reported results do not apply to such zones.

Bulletin of the Seismological Society of America↗

Dynamic plane-strain shear rupture with a slip-weakening friction law calculated by a boundary integral method

A numerical boundary integral method, relating slip and traction on a plane in an elastic medium by convolution with a discretized Green function, can be linked to a slip-dependent friction law on the fault plane. Such a method is developed here in two-dimensional plane-strain geometry. The method is more efficient for a planar source than a finite difference method, and it does not suffer from dispersion of short wavelength components. The solution for a crack growing at constant velocity agrees closely with the analytic solution, and the energy absorbed at the smeared-out crack tip in the numerical calculation agrees with energy absorbed at the analytic singularity. Spontaneous plane-strain shear ruptures can make a transition from sub-Rayleigh to near-P propagation velocity. Results from the boundary integral method agree with earlier results from a finite difference method on the location of this transition in parameter space. The methods differ in their prediction of rupture velocity following the transition. The trailing edge of the cohesive zone propagates at the P-wave velocity after the transition in the boundary integral calculations.

Bulletin of the Seismological Society of America↗

Linear inversion of transmitted acoustic wave fields for three-dimensional modulus and density perturbations using a born-type approximation

A Born approximation is used to linearize the relationship, in the horizontal-wavenumber and frequency domains, between lateral perturbations of modulus and density in a layered half-space and the acoustic wave field observed at the surface when a plane wave is incident from below. The resulting equations can be used to perform a linear inversion of observed acoustic wave fields to obtain lateral perturbations in modulus and density. Since modulus and density effects are separated, gravity observations can be included in the inversion procedure without any assumptions about the relationship between density and acoustic velocity. Tests with synthetic data sets reveal that the inversion method gives useful results when the spatial scales of the inhomogeneities are smaller than several acoustic wavelengths. The inclusion of gravity observations in the inversion reduces the strong negative tradeoff between modulus and density perturbations.

Bulletin of the Seismological Society of America↗

A seismic-refraction profile across the San Andreas, Sargent, and Calaveras faults, west-central California

In 1981, the United States Geological Survey recorded a seismic-refraction profile across the southern Santa Cruz Mountains in west-central California to examine the shallow velocity structure of this seismogenic region. This 40-km-long profile, which consisted of three shotpoints, extended northeastward from near Watsonville, California, to Coyote Lake, crossing the San Andreas, Sargent, and Calaveras faults. This entire region is characterized by a highly heterogeneous upper crust. West of Watsonville, 1 km of alluvium with a velocity of 2.12 km/sec overlies a basement with a velocity 5.45 km/sec. The abrupt deepening of basement by 1.5 km just east of Watsonville at a subsurface fault suggests that the Zayante fault to the north and the Vergeles fault to the south are connected. The Tertiary sediment at the San Andreas fault is 2.5 km thick and has a velocity of 3.34 km/sec. The San Andreas fault is not marked by any distinctive seismic velocity features, whereas a 1-km-wide low velocity zone is found at both the Sargent and Calaveras faults. East of the Sargent fault, the basement surface forms a broad anticlinal structure, with velocities ranging from 5.4 to 6.0 km/sec. From the anticlinal crest, basement dips to the east beneath the Santa Clara Valley and reaches a maximum depth of 1 km on the east side; the overlying alluvium has a velocity of 2.7 km/sec. At the crest of the basement anticlinal structure, a vertical low-velocity zone coincident with exposed serpentine provides strong evidence of faulting.

California↗

Crustal structure of the southern Calaveras fault zone, central California, from seismic refraction investigations

A magnitude 5.7 earthquake on 6 August 1979, within the Calaveras fault zone, near Coyote Lake of west-central California, motivated a seismic-refraction investigation in this area. A northwest-southeast profile along the fault, as well as two fan profiles across the fault were recorded to examine the velocity structure of this region. The analysis of the data reveals a complicated upper crustal velocity structure with strong lateral variations in all directions. The near-surface layers consist of recent sediments with seismic P -wave velocities of 2.6 to 3.2 km/sec. These are underlain by rocks of the Great Valley Sequence which have an average velocity of 4.5 km/sec. The Great Valley Sequence is present along the whole profile; depths range from 4.3 to 4.8 km in the northwest near Anderson Lake and in the southeast of the profile line near Hollister. In the middle of the profile near Coyote Lake, however, this layer thins and we find a laterally limited higher velocity layer (5.1 km/sec) between depths of 2.8 to 4.8 km. The high-velocity zone, which coincides with a gravity high, also correlates spatially with the hypocentral area of the 6 August earthquake and its aftershocks and may therefore represent an asperity on the fault. Velocities within the fault zone were determined from the fan profiles. Near Anderson Lake, a pronounced delay of first arrivals on the fan records indicates a vertical 1- to 2-km-wide near-surface, low-velocity zone along the fault. Near Coyote Lake, the delays observed in the fan records correlate with two subsurface en-echelon fault planes which have been previously identified from lineations in the seismicity pattern. The structure of the lower crust is similar to the neighboring Diablo Range: a 8- to 9-km-thick upper crustal layer with a seismic velocity of 5.7 to 6.3 km/sec is underlain by a 3-km-thick layer with velocity 6.8 km/sec. In accordance with previous studies of the Diablo Range, there are indications of a pronounced lower crustal low-velocity zone between a depth of 17 and 23 km. The presence of this low-velocity zone suggests that higher velocity (igneous?) rocks of the middle crust have been thrust over sedimentary rocks.

California↗