Search USGS⌕ Search

Geology topics

E.R. Engdahl

Publications and source records attributed to E.R. Engdahl.

11 recordsLinked to original sources
USGS2013

Changes in the Seismicity and Focal Mechanism of Small Earthquakes Prior to an MS 6.7 Earthquake in the Central Aleutian Island Arc

On November 4 1977, a magnitude M s 6.7 (m b 5.7) shallow-focus thrust earthquake occurred in the vicinity of the Adak seismographic network in the central Aleutian island arc. The earthquake and its aftershock sequence occurred in an area that had not experienced a similar sequence since at least 1964. About 13 1/2 months before the main shock, the rate of occurrence of very small magnitude earthquakes increased abruptly in the immediate vicinity of the impending main shock. To search for possible variations in the focal mechanism of small events preceding the main shock, a method was developed that objectively combines first-motion data to generate composite focal-mechanism information about events occurring within a small source region. The method could not be successfully applied to the whole study area, but the results show that starting about 10 1/2 months before the November 1977 earthquake, there was a change in the mechanism of small- to moderate-sized earthquakes in the immediate vicinity of the hypocenter and possibly in other parts of the eventual aftershock zone, but not in the surrounding regions.

Book chapter↗
USGS2012

Relocation of local earthquakes by seismic ray tracing

Seismic ray tracing is used to relocate a set of local earthquakes recorded by a dense seismic array in Bear Valley, California. The crustal velocity structure is two-dimensional and incorporates most of the known and inferred velocity differences in the region. Relocated hypocenters fall within the San Andreas fault zone, and P residuals, without elevation or station corrections, now appear to be more internally consistent in the rms sense. Ray takeoff angles and azimuths are distinctly different from those of a uniformly layered earth model. Focal mechanism solutions using only traced data are well determined and have few inconsistent data. Zones of convergence and divergence of rays suggest large variations in observed amplitudes across the local network.

California↗
USGS2012

Teleseismically recorded seismicity before and after the May 7, 1986, Andreanof Islands, Alaska, earthquake

The May 7, 1986, Andreanof Islands earthquake ( M w 8.0) is the largest event to have occurred in that section of the Aleutian arc since the March 9, 1957, Aleutian Islands earthquake ( M w 8.6). Teleseismically well-recorded earthquakes in the region of the 1986 earthquake are relocated with a plate model and with careful attention to the focal depths. The data set is nearly complete for m b ≥4.7 between longitudes 172°W and 179°W for the period 1964 through April 1987 and provides a detailed description of the space-time history of moderate-size earthquakes in the region for that period. Additional insight is provided by source parameters which have been systematically determined for M w ≥5 earthquakes that occurred in the region since 1977 and by a modeling study of the spatial distribution of moment release on the mainshock fault plane. A technically significant component of oblique convergence in the central Aleutian arc results in its breakup into clockwiserotating and westward translating blocks. The western part of the Andreanof block is distinct from and stronger than flanking regions. The greater strength of this block segment and strong coupling along the main thrust zone result in the accumulation of high levels of shear stress, which give rise to great earthquakes near its eastern boundary. The occurrence of the 1986 rupture only 29 years after the 1957 earthquake may indicate that in the central Aleutians M 7+ earthquakes ordinarily do represent the predominant mode of strain release. Segmentation of the main thrust zone into upper and lower planes is supported by spatial and temporal patterns of seismicity and by focal mechanism data. This and other lines of evidence indicate a downdip increase in fault strength and possibly in heterogeneity within the main thrust zone in subduction zone environments. Aftershocks of the 1986 earthquake were bounded to the west by the Andreanof block boundary (Adak Canyon) and to the east by an aseismic segment of the main thrust zone near the subducted extension of the Amlia fracture zone. The aftershock distribution was bounded to the south by a forearc shear zone and to the north by the base of the main thrust zone in the Hawley Ridge segment and by the downdip edge of the upper plane of the main thrust zone in the eastern segment. Aftershocks which occurred near the volcanic line at shallow crustal depths in the upper plate were triggered by the mainshock and manifest a partial decoupling of oblique slip in this region along a west-striking right-lateral fault with low shear strength. Aftershock clustering along the main thrust zone was very similar to the distribution of prior seismicity, suggesting a continuation of long-term processes and the existence of areas with distinct mechanical properties. Interconnecting regions of low seismicity during both the premainshock and aftershock periods coincided with areas of major moment release during the mainshock. Seismicity data prior to the 1986 mainshock gave few clues about the location of the mainshock nucleation point, the mainshock size, and its time of occurrence. However, a large part of the mainshock moment release did coincide with a zone of seismic quiescence monitored by the Adak local seismograph network.

Journal of Geophysical Research Solid Earth↗
USGS2011

Toward a comprehensive catalog of global historical seismicity

The United States Geological Survey (USGS) and the Colorado School of Mines (CSM) have initiated a project to locate more accurately all earthquakes recorded by instruments during the period 1900 to 1963. Seismicity for this period (hereafter referred to as historical seismicity, following Lee et al . [1988]) is still poorly understood, even for basic parameters such as earthquake locations (see Figure 1). In some cases this is the result of inherent limitations in the distribution, response characteristics, and timing of the instruments. However, locations for most of the pre-1964 earthquakes are poorly determined simply because modern data analysis techniques have yet to be applied to the available arrival-time observations, which are mainly preserved as printed bulletins and not in a computer-ready digital format. The arduous task of hand-entering these data has prevented the systematic analysis and relocation of historical seismicity.

Eos, Transactions, American Geophysical Union↗
USGS1999

Seismic subduction of the Nazca Ridge as shown by the 1996-97 Peru earthquakes

By rupturing more than half of the shallow subduction interface of the Nazca Ridge, the great November 12, 1996 Peruvian earthquake contradicts the hypothesis that oceanic ridges subduct aseismically. The mainshock's rupture has a length of about 200 km and has an average slip of about 1.4 m. Its moment is 1.5 x 1028 dyne-cm and the corresponding M(w) is 8.0. The mainshock registered three major episodes of moment release as shown by a finite fault inversion of teleseismically recorded broadband body waves. About 55% of the mainshock's total moment release occurred south of the Nazca Ridge, and the remaining moment release occurred at the southern half of the subduction interface of the Nazca Ridge. The rupture south of the Nazca Ridge was elongated parallel to the ridge axis and extended from a shallow depth to about 65 km depth. Because the axis of the Nazca Ridge is at a high angle to the plate convergence direction, the subducting Nazca Ridge has a large southwards component of motion, 5 cm/yr parallel to the coast. The 900-1200 m relief of the southwards sweeping Nazca Ridge is interpreted to act as a 'rigid indenter,' causing the greatest coupling south of the ridge's leading edge and leading to the large observed slip. The mainshock and aftershock hypocenters were relocated using a new procedure that simultaneously inverts local and teleseismic data. Most aftershocks were within the outline of the Nazca Ridge. A three-month delayed aftershock cluster' occurred at the northern part of the subducting Nazca Ridge. Aftershocks were notably lacking at the zone of greatest moment release, to the south of the Nazca Ridge. However, a lone foreshock at the southern end of this zone, some 140 km downstrike of the mainshock's epicenter, implies that conditions existed for rupture into that zone. The 1996 earthquake ruptured much of the inferred source zone of the M(w) 7.9-8.2 earthquake of 1942, although the latter was a slightly larger earthquake. The rupture zone of the 1996 earthquake is immediately north of the seismic gap left by the great earthquakes (M(w) ~8.8-9.1) of 1868 and 1877. The M(w) 8.0 Antofagasta earthquake of 1995 occurred at the southern end of this great seismic gap. The M(w) 8.2 deep-focus Bolivian earthquake of 1994 occurred directly downdip of the 1868 portion of that gap. The recent occurrence of three significant earthquakes on the periphery of the great seismic gap of the 1868 and 1877 events, among other factors, may signal an increased seismic potential for that zone.

Pure and Applied Geophysics↗
USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗
USGS1992

The international seismological observing period in Africa

The International Seismological Observing Period (ISOP) is a specific time interval designated for enhanced international cooperation in the collection and dissemination of observatory measurements from the global seismographic network. The primary purpose of the ISOP is to strengthen the international infrastructure that supports current seismological practice and increase the cooperation among nations that operate seismological observatories. Measurements, reported by the existing global network and compiled by agencies such as the International Seismological Centre (ISC), are providing new information about earthquakes and the structure of the Earth of fundamental importance to the Earth sciences. However, these data represent but a small fraction of the information contained in the seismograms. One of the goals of the ISOP is to collect improved sets of data. In particular, the measurement and reporting of later-arriving phases, during a fixed ISOP period, from earthquakes selected for detailed observation by the cooperating stations will be encouraged. The use of advanced, digital instrumentation provides an unprecedented opportunity for enhancing the methods of seismogram interpretation and seismic parameter extraction, by the implementation of digital processing methods at seismic observatories worldwide. It must be ensured that this new information will be available to the entire seismological community. It is believed that this purpose is best served with an ISOP that promotes increased on-site processing at digital stations in Africa and elsewhere. Improvements in seismology require truly international cooperation and the educational aspects of seismological practice form one of the goals of the ISOP. Thus, workshops will be needed in Africa to train analysts in ISOP procedures and to introduce them to modern techniques and applications of the data. Participants will, thus, benefit from theoretical results and practical experience that are of direct relevance to their own work.

Tectonophysics↗
USGS1987

Analysis of broadband seismograms from selected IASPEI events

Broadband seismograms of body waves that are flat to displacement and velocity in the frequency range from 0.01 to 5.0 Hz can now be routinely obtained for most earthquakes of magnitude greater than about 5.5. These records are obtained either directly or through multichannel deconvolution of waveforms from digitally recording seismograph stations. In contrast to data from conventional narrowband seismographs, broadband records have sufficient frequency content to define the source-time functions of body waves, even for shallow events for which the source functions of direct and surface-reflected phases may overlap. Broadband seismograms for selected IASPEI events are systematically analysed to identify depth phases and the presence of subevents. The procedure results in improved estimates of focal depth, identification of subevents in complex earthquakes, and better resolution of focal mechanisms. We propose that it is now possible for reporting agencies, such as the National Earthquake Information Center, to use broadband digital waveforms routinely in the processing of earthquake data. ?? 1987.

Physics of the Earth and Planetary Interiors↗
USGS1986

Interpretation of broad-band seismograms from central Aleutian earthquakes.

Broad-band Graefenberg (GRF) array data from 11 moderate-size shallow-depth earthquakes in the central Aleutians have been used to study the effects of focal depth and structure across the arc on observed waveforms. The theoretical results, primarily phase arrival times, suggest that arc structure is responsible for many of the complicated features seen on vertical-component summation seismograms simulated with different instrument responses from the broad-band array data. Except for one trench event, all the earthquakes studied occurred along the plate interface zone, had similar thrust focal mechanisms, and differed only in depth. As a result, the effects of depth phases on observed GRF waveforms across the arc were found to be systematically related to the increase in focal depth along the shallow-dipping seismic zone. -from Authors

Annales Geophysicae, Series B↗
USGS1982

Earthquake location in island arcs

A comprehensive data set of selected teleseismic P-wave arrivals and local-network P- and S-wave arrivals from large earthquakes occurring at all depths within a small section of the central Aleutians is used to examine the general problem of earthquake location in island arcs. Reference hypocenters for this special data set are determined for shallow earthquakes from local-network data and for deep earthquakes from combined local and teleseismic data by joint inversion for structure and location. The high-velocity lithospheric slab beneath the central Aleutians may displace hypocenters that are located using spherically symmetric Earth models; the amount of displacement depends on the position of the earthquakes with respect to the slab and on whether local or teleseismic data are used to locate the earthquakes. Hypocenters for trench and intermediate-depth events appear to be minimally biased by the effects of slab structure on rays to teleseismic stations. However, locations of intermediate-depth events based on only local data are systematically displaced southwards, the magnitude of the displacement being proportional to depth. Shallow-focus events along the main thrust zone, although well located using only local-network data, are severely shifted northwards and deeper, with displacements as large as 50 km, by slab effects on teleseismic travel times. Hypocenters determined by a method that utilizes seismic ray tracing through a three-dimensional velocity model of the subduction zone, derived by thermal modeling, are compared to results obtained by the method of joint hypocenter determination (JHD) that formally assumes a laterally homogeneous velocity model over the source region and treats all raypath anomalies as constant station corrections to the travel-time curve. The ray-tracing method has the theoretical advantage that it accounts for variations in travel-time anomalies within a group of events distributed over a sizable region of a dipping, high-velocity lithospheric slab. In application, JHD has the practical advantage that it does not require the specification of a theoretical velocity model for the slab. Considering earthquakes within a 260 km long by 60 km wide section of the Aleutian main thrust zone, our results suggest that the theoretical velocity structure of the slab is presently not sufficiently well known that accurate locations can be obtained independently of locally recorded data. Using a locally recorded earthquake as a calibration event, JHD gave excellent results over the entire section of the main thrust zone here studied, without showing a strong effect that might be attributed to spatially varying source-station anomalies. We also calibrated the ray-tracing method using locally recorded data and obtained results generally similar to those obtained by JHD. ?? 1982.

Physics of the Earth and Planetary Interiors↗