Search USGSSearch

Geology topics

E.A. Okal

Publications and source records attributed to E.A. Okal.

5 recordsLinked to original sources

The 9 June 94 Bolivian Deep Earthquake: An exceptional event in an extraordinary subduction zone

We investigate the physical setting of the Bolivian shock based on the history of the subducting Nazca plate, intraslab seismicity, deep seismic moment release, and seismic tomography. South America has two broad regions of reverse arc curvature. Subduction constrained to this unique geometry produces slab kinking contortions that may cause unusual slab thickening as they sink to the bottom of the transition zone and encounter resistance to penetration into the lower mantle. Such contortions are observed at intermediate depths “upstream” from the slab source regions of both the great 1970 Colombian and 1994 Bolivian events. Thickening helps explain how the Nazca slab accommodates large seismic source dimensions at depths of 625–650 km.

Geophysical Research Letters

Near-field survey of the 1946 Aleutian tsunami on Unimak and Sanak Islands

The 1946 Aleutian earthquake stands out among tsunamigenic events because it generated both very high run-up near the earthquake source region and a destructive trans-Pacific tsunami . We obtained new data on the distribution of its tsunami in the near field along south-facing coasts between Unimak Pass on the west and Sanak Island on the east by measuring the height of driftwood and beach materials that were deposited by the tsunami above the extreme storm tide level. Our data indicate that (1) the highest measured run-up, which is at the Scotch Cap lighthouse, was 42 m above tide level or about 37 m above present storm tide elevation; (2) run-up along the rugged coast from Scotch Cap for 12 km northwest to Sennett Point is 12-18 m, and for 30 km east of Scotch Cap to Cape Lutke it is 24-42 m; (3) run-up along the broad lowlands bordering Unimak Bight is 10-20 m, and in-undation is locally more than 2 km; (5) run-up diminishes to 8 m or less at the southeast corner of Unimak Island ; (6) no evidence was found for run-up above present storm tides (about 4-5 m above MLLW) on the Ikatan Peninsula or areas along the coast to the west; and (7) run-up above storm tide level in the Sanak Island group is restricted to southwest-facing coasts of Sanak , Long, and Clifford Islands , where it is continuous and locally up to 24 m high. Generation of the tsunami by one or more major earthquake-triggered submarine landslides near the shelf edge south of Unimak Island seems to be the only viable mechanism to account for the data on wave arrival time, run-up heights, and distribution, as well as for unconfirmed anecdotal reports of local postquake increases in water depth and diminished bottom-fisheries productivity. A preliminary hydrodynamic simulation of the local tsunami propagation and run-up using a dipolar model of a possible landslide off Davidson Bank provides an acceptable fit to the characteristics of the distribution of local run-up, with a value at 34 m at the Scotch Cap lighthouse.

Bulletin of the Seismological Society of America

Deep earthquakes beneath the Fiji Basin, SW Pacific: Earth's most intense deep seismicity in stagnant slabs

Previous work has suggested that many of the deep earthquakes beneath the Fiji Basin occur in slab material that has been detached and foundered to the bottom of the transition zone or has been laid down by trench migration in a similar recumbent position. Since nowhere else in the Earth do so many earthquakes occur in slabs stagnated in the transition zone, these earthquakes merit closer study. Accordingly, we have assembled from historical and modern data a comprehensive catalogue of the relocated hypocenters and focal mechanisms of well-located deep events in the geographic area between the bottoms of the main Vanuatu and Tonga Wadati-Benioff zones. Two regions of deep seismogenesis are recognized there: (i) 163 deep shocks have occurred north of 15??S in the Vityaz Group from 1949 to 1996. These seismological observations and the absence of other features characteristic of active subduction suggest that the Vityaz group represents deep failure in a detached slab that has foundered to a horizontal orientation near the bottom of the transition zone. (ii) Another group of nearly 50 'outboard' deep shocks occur between about 450 and 660 km depth, west of the complexly buckled and offset western edge of the Tonga Wadati-Benioff zone. Their geometry is in the form of two or possibly three small-circle arcs that roughly parallel the inferred motion of Tonga trench migration. Earthquakes in the southernmost of these arcs occur in a recumbent high-seismic-wavespeed slab anomaly that connects both to the main inclined Tonga anomaly to the east and a lower mantle anomaly to the west [Van der Hilst, R., 1995. Complex morphology of subducted lithosphere in the mantle beneath the Tonga trench. Nature, Vol. 374, pp. 154-157.]. Both groups show complexity in their focal mechanisms. The major question raised by these observations is the cause of this apparent temporary arrest in the descent of the Tonga slab into the lower mantle. We approach these questions by considering the effects of buoyant metastable peridotite in cold slab material that was detached and rapidly foundered, or was buckled, segmented and laid out in the transition zone.

Physics of the Earth and Planetary Interiors

Metastable mantle phase transformations and deep earthquakes in subducting oceanic lithosphere

Earth's deepest earthquakes occur as a population in subducting or previously subducted lithosphere at depths ranging from about 325 to 690 km. This depth interval closely brackets the mantle transition zone, characterized by rapid seismic velocity increases resulting from the transformation of upper mantle minerals to higher-pressure phases. Deep earthquakes thus provide the primary direct evidence for subduction of the lithosphere to these depths and allow us to investigate the deep thermal, thermodynamic, and mechanical ferment inside slabs. Numerical simulations of reaction rates show that the olivine → spinel transformation should be kinetically hindered in old, cold slabs descending into the transition zone. Thus wedge-shaped zones of metastable peridotite probably persist to depths of more than 600 km. Laboratory deformation experiments on some metastable minerals display a shear instability called transformational faulting. This instability involves sudden failure by localized superplasticity in thin shear zones where the metastable host mineral transforms to a denser, finer-grained phase. Hence in cold slabs, such faulting is expected for the polymorphic reactions in which olivine transforms to the spinel structure and clinoenstatite transforms to ilmenite. It is thus natural to hypothesize that deep earthquakes result from transformational faulting in metastable peridotite wedges within cold slabs. This consideration of the mineralogical states of slabs augments the traditional largely thermal view of slab processes and explains some previously enigmatic slab features. It explains why deep seismicity occurs only in the approximate depth range of the mantle transition zone, where minerals in downgoing slabs should transform to spinel and ilmenite structures. The onset of deep shocks at about 325 km is consistent with the onset of metastability near the equilibrium phase boundary in the slab. Even if a slab penetrates into the lower mantle, earthquakes should cease at depths near 700 km, because the seismogenic phase transformations in the slab are completed or can no longer occur. Substantial metastability is expected only in old, cold slabs, consistent with the observed restriction of deep earthquakes to those settings. Earthquakes should be restricted to the cold cores of slabs, as in any model in which the seismicity is temperature controlled, via the distribution of metastability. However, the geometries of recent large deep earthquakes pose a challenge for any such models. Transformational faulting may give insight into why deep shocks lack appreciable aftershocks and why their source characteristics, including focal mechanisms indicating localized shear failure rather than implosive deformation, are so similar to those of shallow earthquakes. Finally, metastable phase changes in slabs would produce an internal source of stress in addition to those due to the weight of the sinking slab. Such internal stresses may explain the occurrence of earthquakes in portions of lithosphere which have foundered to the bottom of the transition zone and/or are detached from subducting slabs. Metastability in downgoing slabs could have considerable geodynamic significance. Metastable wedges would reduce the negative buoyancy of slabs, decrease the driving force for subduction, and influence the state of stress in slabs. Heat released by metastable phase changes would raise temperatures within slabs and facilitate the transformation of spinel to the lower mantle mineral assemblage, causing slabs to equilibrate more rapidly with the ambient mantle and thus contribute to the cessation of deep seismicity. Because wedge formation should occur only for fast subducting slabs, it may act as a “parachute” and contribute to regulating plate speeds. Wedge formation would also have consequences for mantle evolution because the density of a slab stagnated near the bottom of the transition zone would increase as it heats up and the wedge transforms to denser spinel, favoring the subsequent sinking of the slab into the lower mantle.

Reviews of Geophysics

Frequency-moment distribution of deep earthquakes; Implications for the seismogenic zone at the bottom of slabs

We present a systematic investigation of the variation with depth of the frequency of earthquake occurrence vs. seismic moment based on 16 years of Harvard Centroid Moment Tensor (CMT) solutions. We analyze depth variations of earthquake size distribution in terms of variations in the absolute value of the slope of the regression of the logarithm of the population vs. seismic moment, a quantity known as the β parameter. The shallowest earthquakes (0–50 km depth) exhibit a well-defined and robust size distribution regime characterized by a discontinuous increase in β with increasing moment. Others have shown that this increase probably represents the effects of a physical limit in the dimensions of the area of seismogenic slip of shallow earthquake sources. The population of deep earthquakes in the depth interval 500–600 km shows two markedly different distributions. The deep earthquakes in the Tonga region feature an initially high β value (0.92) at small moments and a lower β value (0.41) at high moments. In contrast, the size distribution of non-Tonga deep events shows the reverse of those changes (β = 0.41 at low moment and β = 1.17 at higher moment). To help explain these observations, we propose a model of deep seismogenesis that assumes three-dimensional earthquake source regions that vary principally in their transverse dimensions. The two-β segment behavior in the Tonga region and other subduction zones is thought to represent, in part, constraints owing to the threshold of completeness of the CMT catalog and to its short time interval of sampling. We interpret the differences between Tonga and other deep Wadati-Benioff zones as being a consequence of Tonga's markedly higher subduction rate and, hence, its colder thermal structure and presumably thicker region of seismogenesis. We interpret the critical moments at which β values change in terms of variations in the transverse thickness of deep seismogenic zones and estimate that it is about 11 km for the Tonga region and about 3 km for other zones at depths of 500–600 km. These results are generally consistent with deep earthquakes being restricted to wedge-shaped regions of peridotite persisting metastably to as deep as 700 km in old, rapidly descending and hence cold slabs. Failure is thought to occur in metastable peridotite by transformational faulting. Great deep earthquakes present special challenges to any theory of deep earthquakes based on slab thermal structure. For example, a continuing question is how such large events can fit in a thermally controlled seismogenic zone that is diminishing in its transverse dimensions with increasing depth. The very concept of a scale-invariant earthquake size distribution may be inappropriate for these rare events and the unusual settings in which they are found.

Physics of the Earth and Planetary Interiors