Kenneth C. Creager, Ling-Yun Chiao, John P. Winchester, Jr., E. Robert Engdahl
Geophysical Research Letters, Vol.22, No. 16, Pages 2321-2324, August 15, 1995
Geometric constraints associated with the concave oceanward bend in the trench and Andes Mountains, known as the Bolivian Orocline, force along-strike compression in the subducting slab. Our calculations of membrane strain rates, in an assumed continuous slab, demonstrate that this can be accommodated by either 10% along-strike compressive strain, or by geometric buckling. In nearly all other subduction zones exhibiting a similar bend in the trench, the compression is accommodated by slab-arching landward of the trench bend. The slab geometry is markedly different in South America. Two regions of nearly flat-lying slab north and south of the Bolivian Orocline result in the 150-km slab depth contour being much straighter than the trench profile. Our calculations show that this geometry helps to relieve along-strike compression. Along the deeper part of the subduction zone, the along-strike compression seems to be accommodated by pronounced along-strike buckling. We infer the approximate geometry of this buckle, in an otherwise aseismic region, from the great Mw 8.3 Bolivian deep earthquake on June 9, 1994, and its aftershocks. The deep slab appears to have a nearly east-west strike in the vicinity of this earthquake, whereas the strike is more nearly north-south to the north and south. The focal mechanisms of deep earthquakes exhibit consistent down-dip compression and along-strike null axes, tracking the contortions in slab geometry. However, our calculations indicate that the particle paths themselves are not down dip everywhere. Instead, where the deep-slab strike is oriented nearly east-west, the particle flow is more nearly along strike than down dip. In general this is not a major problem, but if the slab is anchored in the lower mantle, and moving very slowly, this produces large deformation rates near the base of the upper mantle where the slab is assumed to move at near plate rates. This increased, complex deformation in the deep slab may contribute to the conditions for very large earthquakes in two ways. First, this may dramatically increase deformation rates, and secondly, the increased deformation may lead to anomalous advective thickening of a proposed overdriven olivine wedge. This provides a mechanism for locally thickening the olivine wedge and allowing transformational faulting to occur over a much larger volume than allowed for normal subduction.
Figure 1. Map of South America showing depth contours of the WBZ at 100 km intervals (black lines) and 25 km intervals (red lines). Earthquake hypocentral depths are given in the key. Lower-focal hemisphere projections of focal mechanisms show compressive axes (white) and tension (red). Selected particle paths are projected from the slab surface straight up to the surface (solid lines). Grey lines show paths that would be taken if the oceanic plate were to continue rotating across the Earth's surface. Ticks mark time at 2 Ma intervals.
Figure 3. Map views showing theoretical slab contours (a and d, at 100 km intervals starting with 0 km), membrane strain rate tensor orientations (b and e), and effective strain rates (c and f). Lines (b and e) show directions of principal axes of strain-rate tensor. Thick bars indicate in-plane compression, thin lines are extension. Size of the strain rate is proportional to line length, but is truncated at a value of 10.