Flat slab subduction

Flat slab subduction is characterized by a low subduction angle (<30 degrees to horizontal) beyond the seismogenic layer and a resumption of normal subduction far from the trench.^[1] A slab refers to the subducting lower plate. A broader definition of flat slab subduction includes any shallowly dipping lower plate, as in western Mexico. Flat slab subduction is associated with the pinching out of the asthenosphere, an inland migration of arc magmatism (magmatic sweep), and an eventual cessation of arc magmatism.^[2] The coupling of the flat slab to the upper plate is thought to change the style of deformation occurring on the upper plate's surface and form basement-cored uplifts like the Rocky Mountains.^[2][3] The flat slab also may hydrate the lower continental lithosphere^[2] and be involved in the formation of economically important ore deposits.^[4] During the subduction, a flat slab itself may deform or buckle, causing sedimentary hiatus in marine sediments on the slab.^[5] The failure of a flat slab is associated with ignimbritic volcanism and the reverse migration of arc volcanism.^[2] Multiple working hypotheses about the cause of flat slabs are subduction of thick, buoyant oceanic crust (15–20 km)^[6] and trench rollback accompanying a rapidly overriding upper plate and enhanced trench suction.^[7] The west coast of South America has two of the largest flat slab subduction zones.^[2] Flat slab subduction is occurring at 10% of subduction zones.^[3]

A diagram representing flat slab subduction

History of idea

The idea has its beginnings in the late 1970s.^[8] Seismic studies of the Andean margin seemed to show a zone of subhorizontal lower plate at a depth of 100 km. The Cornell-Carnegie debate between Cornell University geophysicists and workers at the Carnegie Institute of Washington centered on whether local deployments of seismometers would yield better results than looking at global (teleseismic) data. The Carnegie Institution seemed to have won the day with the local deployment imaging the flat slab where teleseismic data argued for a shallowing dipping slab with no near horizontal zone.^[9] The idea was taken up to explain the Laramide orogeny, as the flat slab subduction zones on the Andean margin are associated with more inboard surface deformation and magmatic gaps.^[2] Flat slab subduction is an active area of research; the causal mechanisms for its occurrence have not been sorted out.

Causal mechanisms and consequences of flat slab subduction

Causal mechanisms

There are several working hypotheses for the initiation of flat slab subduction. The buoyant ridge hypothesis seems to be favored at the moment.^[3]

Subduction of buoyant oceanic crust

The subduction of bathymetric highs such as aseismic ridges, oceanic plateaus, and seamounts has been posited as the primary driver of flat slab subduction.^[3] The Andean flat slab subduction zones, the Peruvian slab and the Pampean (Chilean) flat slab, are spatially correlated with the subduction of bathymetric highs, the Nazca Ridge and the Juan Fernandéz Ridge, respectively. The thick, buoyant oceanic crust lowers the density of the slab, and the slab fails to sink into the mantle after coming to a shallow depth (~100 km) due to the lessened density contrast.^[6] This is supported by the fact that all slabs are under ~50 Ma.^[10] However, there are cases where aseismic ridges on the same scale as the Nazca Ridge are subducting normally, and cases where flat slabs are not associated with bathymetric highs.^[11] There are few flat slabs in the Western Pacific in areas associated with the subduction of bathymetric highs.^[12] Geodynamic modeling has called into question whether buoyant oceanic crust alone can generate flat slab subduction.^[10]

Trenchward motion of overriding plate with cratonic keel

Another explanation for slab flattening is the lateral movement of the overriding plate in a direction opposite to that of the downgoing slab. The overriding plate is often equipped with a cratonic keel of thick continental lithosphere which, if close enough to the trench, can impinge upon the flow in the mantle wedge.^[7] Trench suction is included in this causal mechanism. Trench suction is induced by the flow of the asthenosphere in the mantle wedge area; trench suction increases with subduction velocity, a decrease of the mantle wedge thickness, or an increase in the mantle wedge viscosity.^[13] Trench retreat is motion of the trench in a direction opposite to that of plate convergence thought to be related to the position of the trench along the larger subduction zone with retreat occurring near the edges of subduction zones.^[14] Modeling experiments have shown that if the cratonic lithosphere is thick and the trench retreats, the shutdown of the mantle wedge increases trench suction to an extent that the slab flattens.^[7]

Consequences

Delay in the eclogitization

Eclogite is a dense (3.5 g/cu. cm), garnet-bearing rock that is formed as the oceanic crust subducts to zones of high pressure and temperature. The reaction that forms eclogite dehydrates the slab and hydrates the mantle wedge above. The now denser slab more effectively sinks.^[15] A delay in eclogitization could arise through the subduction of zone thicker oceanic lithosphere without deeply penetrating faults. Oceanic crust is normally faulted at the trench rise by the bending of the plate as it subducts. This may be an effect or a cause of flat slab subduction, but it seems as though it is more likely an effect. A resumption of normally dipping subduction beyond the flat slab portion is associated with the eclogite reaction, and the amount of time needed to accumulate enough eclogite for the slab to start sinking may be what limits temporal scale of flat slab subduction.^[6]

Magmatic gaps and adakitic volcanism

As the subducting plate flattens there is an inboard migration in the magmatic arc that can be tracked. In the Chilean flat slab region (~31–32 degrees S), around 7–5 Ma there was an eastward migration, broadening and gradual shutdown down of the volcanic arc associated with slab flattening.^[16] This occurs as the previous magmatic arc position on the upper plate (100–150 km above subducting plate) is no longer aligned with the zone of partial melting above the flattening slab.^[17] The magmatic arc migrates to a new location that coincides with the zone of partial melting above the flattening slab. Magmatism before the Laramide orogeny migrated all the way to western South Dakota.^[2] Eventually, the magmatic activity above the flat slab may completely cease as the subducting plate and upper plate pinch out the mantle wedge.^[2] Upon the failure of the flat slab, the mantle wedge can again start circulating hot asthenosphere (1300 degrees C) in an area that has been heavily hydrated, but that had not produced any melt; this leads to widespread ignimbritic volcanism, which is seen in both the Andean flat slab effected regions and the western United States.^[18]

Adakites are dacitic and andesitic magmas that are highly depleted in heavy rare-earth elements and high strontium/yttrium ratios and may be derived of melting of the oceanic crust.^[17] Adakites are thought to erupt or be emplaced during the transition from normally dipping subduction to flat subduction as the magmatic arc widens and migrates more inland.^[17] Adakitic rocks can be seen in modern Ecuador,^[19] a possible incipient flat slab zone, and in central Chile there are 10-5 Ma adakitic rocks.^[20] Thus, adakitic rocks could be used as marker of past episodes of flat slab subduction.

Surface deformation

See also: Andean foreland basins

Flat slabs are thought to result in zones of broad, diffuse deformation in the upper plate located far landward from the trench.^[3] Flat slab subduction is associated with basement-cored uplifts also known as "thick-skinned" deformation of the overriding plate like the Sierra Pampeanas in South America possibly associated with the subduction of the Juan Fernandéz Ridge.^[21] These areas of basement-cored uplifts are visually correlated with flat slab subduction zones.^[16] In contrast, "thin-skinned" deformation is the normal mode of upper plate deformation, and does not involve basement rock. Crustal shortening is observed to extend farther inland than in normally dipping subduction zones; the Sierra Pampeanas are over 650 km east of the trench axis.^[21] Flat slabs have been used as an explanation for the Laramide Orogeny^[18] and the central Altiplano-Puna region.^[22] Another interesting feature that may be associated with the flat slab subduction of the Nazca Ridge is the Fitzcarrald arch located in the Amazonian Basin. The Fitzcarrald arch is a long-wavelength, linear topographic feature extending from eastern Peru to western Brazil beyond the Subandean thrust front into an undeformed area and rising ~600 masl.^[23] The Fitzcarrald arch has the effect of splitting the Amazonian Basin into three subbasins: northern Amazonian foreland basin, southern Amazonian foreland basin, and the eastern Amazonian foreland basin.^[24][25]

Seismicity

The shape of the flat slab is constrained through earthquakes within the subducting slab and the interface between the upper plate and the subducting slab.^[16] Flat slab zones along the Andean margin release 3–5 times more energy through upper plate earthquakes than adjacent, more steeply dipping subduction zones.^[3] Upper plate earthquake focal mechanisms indicate that stress is aligned parallel with motion of the plate, and that stress is transmitted high into the upper plate from the lower.^[26] The reason for this enhanced seismicity is more effective coupling of the upper and lower plates. In normal subduction zones the coupling interface, the area in which the two plates are in close proximity, between the two plates is ~100–200 km long, but in flat slab subduction zones the coupling interface is much longer, 400–500 km.^[26] Although the lower lithosphere of the upper deforms plastically, numerical modeling has shown stress can be transmitted to crustal regions which behave in a brittle fashion.^[27] Along the subducting plate seismicity is more variable, especially intermediate-depth earthquakes. The variability may be controlled by the thickness of the crust and how efficiently it can release water. Thick crust that is not as deeply fractured by trench rise normal faulting may not dehydrate rapidly enough to induce intermediate-depth earthquakes.^[1] The Peruvian flat slab lacks significant intermediate-depth earthquakes and is associated with the subduction of the ~17 km thick Nazca Ridge.^[1]

Andean flat slabs

In the late 1970s early research recognized the unique nature of the two large flat slab subduction zones along the Andean margin of South America.^[28][29] Two large and one smaller current flat slab subduction segments exist along the Andean margin: the Peruvian, Pampean, and the Bucaramanga. Three Cenozoic flat slab segment are also known: Altiplano, Puna, and Payenia.

The Peruvian flat slab is located between the Gulf of Guayaquil (5 degrees S) and Arequipa (14 degrees S), extending ~1500 km along the strike of the subduction zone. The Peruvian flat slab is the largest in the world,^[3] and extends ~700 km inboard from the trench axis. The subducting plate starts at a dip of 30 degrees then flattens out at a depth of 100 km under the Eastern Cordillera and Subandean zone.^[30] The segment is visually correlated with the subduction of the Nazca Ridge, an aseismic ridge with thickened crust. The second highest zone in the Andes, Cordillera Blanca, is associated with the Peruvian flat slab segment and uplift of basement-cored blocks. Volcanism in the area ceased in the Late Miocene (11-5 Ma). Plate reconstructions time the collision of the Nazca Ridge with the subduction zone at 11.2 Ma at 11 degree S, which implies that the northern extent of the Peruvian flat slab may require some other subducted feature like an oceanic plateau. A putative subducted plateau, the Inca Plateau, has been argued for.^[31]

The Pampean or Chilean flat slab segment is located between 27 degrees S and 33 degrees S, extending ~550 km along the strike of the subduction zone. The Pampean flat slab similarly extends ~700 km inboard from the trench axis. The segment is visually correlated with the Juan Fernandez Ridge, and the highest peak in the Andes, the non-volcanic Aconcagua (6961 m). This area has undergone the same "thick-skinned" deformation, leading to the high mountain peaks.

The Bucaramanga segment was recognized in early eighties from limited seismological evidence.^[32] The segment is encompassed between 6 and 9 degrees N in Colombia, extending ~350 km along the strike of the subduction zone.

Other flat slabs

There are several other flat slab segments that warrant a mention:^[3]

• Alaskan: 145–150 degrees W along the Aleutian Trench associated with the Yakutat microplate
• Costa Rica: 82–84 degrees W associated with the Cocos Ridge
• Mexico: 96–100 degrees W associated with the Tenhuantepec Ridge
• Cascadian, United States: 46–49 degrees N associated with the subduction of young oceanic crust

Economic geology

Subduction of thick oceanic crust could be linked with the metallogenesis of copper and gold deposits.^[4] The 10 largest young (<18 Ma) gold deposits in South America are associated with flat slab segments.^[4] Enhanced metallogenesis may be caused by the cessation of magmatism in the arc allowing the conservation of sulfur-rich volatiles.^[4] The failure of the putative flat slab under western North America may have been vital in producing Carlin-type gold deposits.^[33]

Early Earth subduction

Early Earth's mantle was hotter and it has been proposed that flat slab subduction was the dominant style.^[34] Computer modeling has shown that an increase in oceanic plate buoyancy associated with enhanced oceanic crust production would have been counteracted by decreased mantle viscosity, so flat slab subduction would not have been dominant or non-existent.^[10]

References

1. Kumar, Abhash; Wagner, Lara S.; Beck, Susan L.; Long, Maureen D.; Zandt, George; Young, Bissett; Tavera, Hernando; Minaya, Estella (2016-05-01). "Seismicity and state of stress in the central and southern Peruvian flat slab". Earth and Planetary Science Letters. 441: 71–80. Bibcode:2016E&PSL.441...71K. doi:10.1016/j.epsl.2016.02.023.
2. Humphreys, Eugene; Hessler, Erin; Dueker, Kenneth; Farmer, G. Lang; Erslev, Eric; Atwater, Tanya (2003-07-01). "How Laramide-Age Hydration of North American Lithosphere by the Farallon Slab Controlled Subsequent Activity in the Western United States". International Geology Review. 45 (7): 575–595. Bibcode:2003IGRv...45..575H. doi:10.2747/0020-6814.45.7.575. ISSN 0020-6814. S2CID 15349233.
3. Gutscher, Marc-André; Spakman, Wim; Bijwaard, Harmen; Engdahl, E. Robert (2000-10-01). "Geodynamics of flat subduction: Seismicity and tomographic constraints from the Andean margin". Tectonics. 19 (5): 814–833. Bibcode:2000Tecto..19..814G. doi:10.1029/1999TC001152. ISSN 1944-9194.
4. Rosenbaum, Gideon; Giles, David; Saxon, Mark; Betts, Peter G.; Weinberg, Roberto F.; Duboz, Cecile (2005-10-30). "Subduction of the Nazca Ridge and the Inca Plateau: Insights into the formation of ore deposits in Peru". Earth and Planetary Science Letters. 239 (1–2): 18–32. Bibcode:2005E&PSL.239...18R. doi:10.1016/j.epsl.2005.08.003.
5. Li, Yong-Xiang; Zhao, Xixi; Jovane, Luigi; Petronotis, Katerina E.; Gong, Zheng; Xie, Siyi (2015-12-01). "Paleomagnetic constraints on the tectonic evolution of the Costa Rican subduction zone: New results from sedimentary successions of IODP drill sites from the Cocos Ridge". Geochemistry, Geophysics, Geosystems. 16 (12): 4479–4493. Bibcode:2015GGG....16.4479L. doi:10.1002/2015GC006058. ISSN 1525-2027.
6. Antonijevic, Sanja Knezevic; Wagner, Lara S.; Kumar, Abhash; Beck, Susan L.; Long, Maureen D.; Zandt, George; Tavera, Hernando; Condori, Cristobal (2015-08-13). "The role of ridges in the formation and longevity of flat slabs". Nature. 524 (7564): 212–215. Bibcode:2015Natur.524..212A. doi:10.1038/nature14648. ISSN 0028-0836. PMID 26268192. S2CID 205244754.
7. Manea, Vlad C.; Pérez-Gussinyé, Marta; Manea, Marina (2012-01-01). "Chilean flat slab subduction controlled by overriding plate thickness and trench rollback". Geology. 40 (1): 35–38. Bibcode:2012Geo....40...35M. doi:10.1130/G32543.1. ISSN 0091-7613.
8. Uyeda, S.; Sacks, I. Selwyn (1977-01-05). "Subduction zones, mid-ocean ridges, oceanic trenches and geodynamicsInterrelationships between volcanism, seismicity, and anelasticity in western South America". Tectonophysics. 37 (1): 131–139. doi:10.1016/0040-1951(77)90043-9.
9. Hasegawa, Akira; Sacks, I. Selwyn (1981-06-10). "Subduction of the Nazca Plate beneath Peru as determined from seismic observations". Journal of Geophysical Research: Solid Earth. 86 (B6): 4971–4980. Bibcode:1981JGR....86.4971H. doi:10.1029/JB086iB06p04971. ISSN 2156-2202. S2CID 53443213.
10. van Hunen, Jeroen; van den Berg, Arie P; Vlaar, Nico J (2004-08-16). "Various mechanisms to induce present-day shallow flat subduction and implications for the younger Earth: a numerical parameter study". Physics of the Earth and Planetary Interiors. Plumes and Superplumes. 146 (1–2): 179–194. Bibcode:2004PEPI..146..179V. doi:10.1016/j.pepi.2003.07.027.
11. Skinner, Steven M.; Clayton, Robert W. (2013-06-01). "The lack of correlation between flat slabs and bathymetric impactors in South America" (PDF). Earth and Planetary Science Letters. 371–372: 1–5. Bibcode:2013E&PSL.371....1S. doi:10.1016/j.epsl.2013.04.013.
12. Rosenbaum, Gideon; Mo, Won (2011-04-01). "Tectonic and magmatic responses to the subduction of high bathymetric relief". Gondwana Research. Island Arcs: Their role in growth of accretionary orogens and mineral endowment. 19 (3): 571–582. Bibcode:2011GondR..19..571R. doi:10.1016/j.gr.2010.10.007.
13. Stevenson, D. J.; Turner, J. S. (1977-11-24). "Angle of subduction". Nature. 270 (5635): 334–336. Bibcode:1977Natur.270..334S. doi:10.1038/270334a0. S2CID 4205429.
14. Schellart, W. P.; Freeman, J.; Stegman, D. R.; Moresi, L.; May, D. (2007-03-15). "Evolution and diversity of subduction zones controlled by slab width". Nature. 446 (7133): 308–311. Bibcode:2007Natur.446..308S. doi:10.1038/nature05615. ISSN 0028-0836. PMID 17361181. S2CID 4420049.
15. Pennington, Wayne D. (1984-02-20). "Geodynamics of Back-Arc Regions The effect of oceanic crustal structure on phase changes and subduction". Tectonophysics. 102 (1): 377–398. doi:10.1016/0040-1951(84)90023-4.
16. Alvarado, Patricia; Pardo, Mario; Gilbert, Hersh; Miranda, Silvia; Anderson, Megan; Saez, Mauro; Beck, Susan (2009-06-01). Flat-slab subduction and crustal models for the seismically active Sierras Pampeanas region of Argentina. Vol. 204. pp. 261–278. doi:10.1130/2009.1204(12). ISBN 9780813712048. ISSN 0072-1069. {{cite book}}: |journal= ignored (help)
17. Gutscher, Marc-André; Maury, René; Eissen, Jean-Philippe; Bourdon, Erwan (2000-06-01). "Can slab melting be caused by flat subduction?". Geology. 28 (6): 535–538. Bibcode:2000Geo....28..535G. doi:10.1130/0091-7613(2000)28<535:csmbcb>2.0.co;2. ISSN 0091-7613.
18. Humphreys, Eugene (2009-06-01). "Relation of flat subduction to magmatism and deformation in the western United States". Backbone of the Americas: Shallow Subduction, Plateau Uplift, and Ridge and Terrane Collision. Vol. 204. pp. 85–98. doi:10.1130/2009.1204(04). ISBN 9780813712048. ISSN 0072-1069. {{cite book}}: |journal= ignored (help)
19. Gutscher, M. -A; Malavieille, J; Lallemand, S; Collot, J. -Y (1999-05-15). "Tectonic segmentation of the North Andean margin: impact of the Carnegie Ridge collision". Earth and Planetary Science Letters. 168 (3–4): 255–270. Bibcode:1999E&PSL.168..255G. doi:10.1016/S0012-821X(99)00060-6.
20. Litvak, Vanesa D.; Poma, Stella; Kay, Suzanne Mahlburg (2007-09-01). "Paleogene and Neogene magmatism in the Valle del Cura region: New perspective on the evolution of the Pampean flat slab, San Juan province, Argentina". Journal of South American Earth Sciences. 24 (2–4): 117–137. Bibcode:2007JSAES..24..117L. doi:10.1016/j.jsames.2007.04.002. hdl:11336/77992.
21. Jordan, T. E.; Allmendinger, R. W. (1986). "Sign In". American Journal of Science. 286 (10): 737–764. doi:10.2475/ajs.286.10.737.
22. Cahill, Thomas; Isacks, Bryan L. (1986-04-01). "An apparent double-planed Benioff Zone beneath northern Chile resulting from misidentification of reflected phases". Geophysical Research Letters. 13 (4): 333–336. Bibcode:1986GeoRL..13..333C. doi:10.1029/GL013i004p00333. ISSN 1944-8007.
23. Espurt, N.; Baby, P.; Brusset, S.; Roddaz, M.; Hermoza, W.; Regard, V.; Antoine, P.-O.; Salas-Gismondi, R.; Bolaños, R. (2007-06-01). "How does the Nazca Ridge subduction influence the modern Amazonian foreland basin?". Geology. 35 (6): 515–518. Bibcode:2007Geo....35..515E. doi:10.1130/G23237A.1. ISSN 0091-7613.
24. Baby, P.; Guyot, J. L.; Deniaud, Y.; Zubieta, D.; Christophoul, F.; Rivadeneira, M.; Jara, F. (1999-01-01). "The High Amazonian Basin: tectonic control and mass balance" (PDF). International Symposium MANAUS 99, Hydrological and Geochemical Process in Large-scale Rivers Basins: Manaus (Brazil).
25. Kronberg, B.I.; Fralick, P.W.; Benchimol, R.E. (1998-09-01). "Late Quaternary sedimentation and palaeohydrology in the Acre foreland basin, SW Amazonia". Basin Research. 10 (3): 311. Bibcode:1998BasR...10..311K. doi:10.1046/j.1365-2117.1998.00067.x. ISSN 1365-2117. S2CID 140183791.
26. Gutscher, Marc-André (2002-04-01). "Andean subduction styles and their effect on thermal structure and interplate coupling". Journal of South American Earth Sciences. Flat-Slab Subduction in the Andes. 15 (1): 3–10. Bibcode:2002JSAES..15....3G. doi:10.1016/S0895-9811(02)00002-0.
27. Spencer, Jon E. (1994-01-01). "A numerical assessment of slab strength during high-and low-angle subduction and implications for Laramide orogenesis". Journal of Geophysical Research: Solid Earth. 99 (B5): 9227–9236. Bibcode:1994JGR....99.9227S. doi:10.1029/94jb00503.
28. Ramos, Victor A.; Folguera, Andrés (2009-01-01). "Andean flat-slab subduction through time". Geological Society, London, Special Publications. 327 (1): 31–54. Bibcode:2009GSLSP.327...31R. doi:10.1144/SP327.3. ISSN 0305-8719. S2CID 43604314.
29. Barazangi, Muawia; Isacks, Bryan L. (1976-11-01). "Spatial distribution of earthquakes and subduction of the Nazca plate beneath South America". Geology. 4 (11): 686–692. Bibcode:1976Geo.....4..686B. doi:10.1130/0091-7613(1976)4<686:sdoeas>2.0.co;2. ISSN 0091-7613.
30. Dorbath, L.; Dorbath, C.; Jimenez, E.; Rivera, L. (1991-01-01). "Seismicity and tectonic deformation in the Eastern Cordillera and the sub-Andean zone of central Peru" (PDF). Journal of South American Earth Sciences. 4 (1): 13–24. Bibcode:1991JSAES...4...13D. doi:10.1016/0895-9811(91)90015-D.
31. Gutscher, M. -A.; Olivet, J. -L.; Aslanian, D.; Eissen, J. -P.; Maury, R. (1999-09-15). "The "lost inca plateau": cause of flat subduction beneath peru?". Earth and Planetary Science Letters. 171 (3): 335–341. Bibcode:1999E&PSL.171..335G. doi:10.1016/S0012-821X(99)00153-3.
32. Pennington, Wayne D. (1981-11-10). "Subduction of the Eastern Panama Basin and seismotectonics of northwestern South America". Journal of Geophysical Research: Solid Earth. 86 (B11): 10753–10770. Bibcode:1981JGR....8610753P. doi:10.1029/JB086iB11p10753. ISSN 2156-2202.
33. Muntean, John L.; Cline, Jean S.; Simon, Adam C.; Longo, Anthony A. (2011-02-01). "Magmatic-hydrothermal origin of Nevada/'s Carlin-type gold deposits". Nature Geoscience. 4 (2): 122–127. Bibcode:2011NatGe...4..122M. doi:10.1038/ngeo1064. ISSN 1752-0894.
34. Abbott, Dallas; Drury, Rebecca; Smith, Walter H. F. (1994-10-01). "Flat to steep transition in subduction style". Geology. 22 (10): 937–940. Bibcode:1994Geo....22..937A. doi:10.1130/0091-7613(1994)022<0937:ftstis>2.3.co;2. ISSN 0091-7613.