skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Earthquake Focal Mechanisms and Stress Field for the Intermediate‐Depth Cauca Cluster, Colombia
Abstract Beneath southwestern Colombia, intermediate‐depth earthquakes in the Cauca cluster locate in the subducting Nazca plate and in two columns extending ~40‐km into the mantle wedge above the slab. To investigate the cluster, we determine focal mechanisms for 69 small‐to‐moderate‐sized (2.3 ≤ Ml ≤ 4.7) earthquakes in the cluster by fitting short‐period body wave waveforms using the cut‐and‐paste method. The focal mechanisms have various faulting types and variably oriented nodal planes. We invert the focal mechanisms of the intraslab earthquakes for the intraslab stress field but cannot fit the region with a homogeneous stress tensor. We find that the principal stress axes rotate with the slab geometry, which has a concave shape and increases in dip angle from north to south. The northern region has slab normal compression and similar‐magnitude maximum and intermediate principal stresses. The minimum stress axis is oriented ~41° counterclockwise from the downdip direction. In the steeper southern region, the intermediate stress axis orients in the downdip direction. Deviation from a typical downdip extensional stress field may result from a buoyant young slab, an eastward mantle flow push, and/or along‐strike compression from the concave shape of the slab. This stress field would allow slip along preexisting faults of various orientations, such as the trench‐perpendicular seafloor features presently observed offshore, and contribute to the apparent heterogeneity of the intraslab stress field. The mantle wedge earthquakes also have various focal mechanisms but tend to have a subvertical nodal plane that aligns with the earthquake locations.  more » « less
Award ID(s):
1760802 1141905
PAR ID:
10370746
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
124
Issue:
1
ISSN:
2169-9313
Page Range / eLocation ID:
p. 822-836
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The subducted old and cold Pacific Plate beneath the young Philippine Sea Plate at the Izu‐Bonin trench over the Cenozoic hosts regional deep earthquakes. We investigate slab morphology and stress regimes under different trench motion histories with mantle convection models. Viscosity, temperature, and deviatoric stress are inherently heterogeneous within the slab, which we link to the occurrence of isolated earthquakes. Models expand on previous suggestions that observed slab morphology variations along the Izu‐Bonin subduction zone, exhibited as shallow slab dip angles in the north and steeper dip angles in the south, are mainly due to variations in the rate of trench retreat from the north (where it is fast) to the south (where it is slow). Geodynamic models consistent with the regional plate tectonics, including oceanic plate age, plate convergence rate, and trench motion history, reproduce the seismologically observed principal stress direction and slab morphology. We suggest that the isolated ~680 km deep, 30 May 2015 Mw 7.9 Bonin Islands earthquake, which lies east of the well‐defined Benioff zone and has its principal compressional stress direction oriented toward the tip of the previously defined Benioff zone, can be explained by Pacific slab buckling in response to the slow trench retreat. 
    more » « less
  2. Abstract The extremely oblique Indo‐Burma subduction zone exhibits dextral strike‐slip faulting along the Sagaing, Kabaw, and Churachandpur‐Mao Faults as well as east‐west shortening between the Sagaing Fault and Bengal Basin. Through regional stress analysis, considering areas from central Tibet, around the eastern Himalaya Syntaxis, to Burma, it has been determined that the principal compressive stress directions align with the principal strain rates. The northeast‐southwest oriented compressive stress direction from the western Shan Plateau continues into Burma. Notably, P axes align with the topographic gradients, and T axes are sub‐parallel to the topographic contours in the Shan Plateau region south of 27°N. These stress patterns are consistent with a gravitational potential energy induced crustal and mantle flow. The alignment of the fast shear wave with the maximum strain rate and the colinear NW‐SE to E‐W fast direction of the SKS wave and T axis determined from focal mechanisms in the Shan Plateau suggest that the mantle lithosphere deforms in concert with the crust. We suggest crust and mantle flow south of the Red River Fault has resulted in widening of the lithosphere in the Shan Plateau in an east‐west direction. Therefore, the Sagaing Fault has bowed approximately 50–100 km westward if we assume that the Sagaing Fault was originally straight. Our results of regional stress inversion are consistent with late Miocene to present E‐W shortening in the Indo‐Burma subduction zone resulting from the release of gravitational potential energy from the central Tibetan Plateau. 
    more » « less
  3. Abstract The radial viscosity of the mantle is generally thought to increase by ∼10–100 times from the upper to lower mantle with a putative, abrupt increase at 660 km depth. Recently, a low viscosity channel (LVC) between 660 and 1,000 km has been suggested. We conduct a series of time‐dependent flow models with viscosity either increasing or decreasing at 660 km depth while tracking slab structure, state‐of‐stress, and geoid. We find that a LVC will lower the amplitude of long wavelength (>5,000 km) geoid highs over slabs, with amplitudes <10 m in height, while increasing the slab dip angle and downdip tension in the upper 300 km of slabs. A viscosity increase at 660 km gives rise to strong downdip compression throughout a slab and this pattern will largely go away with the introduction of the LVC. In addition, the endothermic phase change at 660 km depth can substantially affect the stress distribution within slabs but has a minor influence on the geoid. Models that fit the observed long wavelength geoid and observed focal mechanism in the western Pacific favor models without the presence of the LVC between 660 km and 1,000 km depths. 
    more » « less
  4. Hydration of the subduction zone forearc mantle wedge influences the downdip distribution of seismicity, the availability of fluids for arc magmatism, and Earth's long term water cycle. Reconstructions of present‐day subduction zone thermal structures using time‐invariant geodynamic models indicate relatively minor hydration, in contrast to many geophysical and geologic observations. We pair a dynamic, time‐evolving thermal model of subduction with phase equilibria modeling to investigate how variations in slab and forearc temperatures from subduction infancy through to maturity contribute to mantle wedge hydration. We find that thermal state during the intermediate period of subduction, as the slab freely descends through the upper mantle, promotes extensive forearc wedge hydration. In contrast, during early subduction the forearc is too hot to stabilize hydrous minerals in the mantle wedge, while during mature subduction, slab dehydration dominantly occurs beyond forearc depths. In our models, maximum wedge hydration during the intermediate phase is 60%–70% and falls to 20%–40% as quasi‐steady state conditions are approached during maturity. Comparison to global forearc H2O capacities reveals that consideration of thermal evolution leads to an order of magnitude increase in estimates for current extents of wedge hydration and provides better agreement with geophysical observations. This suggests that hydration of the forearc mantle wedge represents a potential vast reservoir of H2O, on the order of 3.4–5.9 × 1021 g globally. These results provide novel insights into the subduction zone water cycle, new constraints on the mantle wedge as a fluid reservoir and are useful to better understand geologic processes at plate margins. 
    more » « less
  5. Abstract Small‐scale intraslab heterogeneity is well documented seismically in multiple subduction zones, but its nature remains elusive. Previous efforts have been mostly focusing on the scattering strength at intermediate depth (<350 km), without constraining its evolution as a function of depth. Here, we illustrate that the inter‐source interferometry method, which turns deep earthquakes into virtual receivers, can resolve small‐scale intraslab heterogeneity in the mantle transition zone. The interferometric waveform observations in the Japan subduction zone require weak scattering (<1.0%) within the slab below 410 km. Combining with previous studies that suggest high heterogeneity level (∼2.5%) at intermediate depth, we conclude that the small‐scale intraslab heterogeneity weakens as slabs subduct. We suggest that the heterogeneities are caused by intraslab hydrous minerals, and the decrease in their scattering strength with depth reveals processes associated with dehydration of subducting slabs. 
    more » « less