Inversion de polarité de subduction

L'inversion de polarité de subduction (ou SPR pour (en) Subduction Polarity Reversal) est un processus géologique dans lequel deux plaques convergentes changent de rôle : la plaque sus-jacente devient la plaque descendante, et vice versa. Deux éléments de base constituent une zone de subduction, une plaque de contrôle et d'une plaque de subduction[1]. Quand deux plaques se rapprochent l'une de l'autre sous l'effet des forces tectoniques, la plaque dite supérieure coulissera sur la plaque subductrice[1]. Ce type d'évènement tectonique se retrouve à de nombreuses limites de plaques[1].

Le concept d'inversion de la polarité de subduction

Certains géologues pensent que les rôles de plaque supérieure et subductrice ne restent pas indéfiniment les mêmes[2]. Une inversion de rôle peut intervenir, ce qui signifie que la plaque initialement subductée en dessous deviendra la plaque dominante[2]. Ce phénomène est appelé interrupteur de subduction[3], retournement de la polarité de subduction[4] ou inversion de polarité de subduction[2].

On trouve sur la carte tectonique des exemples de systèmes de subduction avec inversion de polarité de subduction :

Arrière-plan

modifier

Le phénomène d'inversion de polarité de subduction a été identifié lors de la collision d'un système de subduction intra-océanique[12], qui est la collision de deux plaques océaniques[1]. Lorsque deux plaques océaniques migrent l’une vers l’autre, l’une se subducte en dessous de l’autre. Généralement, la plaque océanique avec une densité plus élevée se subducte en dessous et l'autre l'emporte sur la dalle descendante[1]. Le processus se poursuit jusqu'à ce qu'une marge continentale flottante située au sommet de la plaque subductrice soit introduite dans la dalle descendante[2]. La subduction de la dalle devient plus lente et peut même cesser[2],[4]. Les géologues proposent différents modèles possibles pour prédire quelle sera la prochaine étape du système de subduction intra-océanique impliquant une croûte continentale flottante[2],[4]. L’un des résultats possibles est l’inversion de polarité de subduction[4],[11],[12],[13],[14],[15].

Modèles d'inversion de polarité de subduction

modifier

Même si de nombreux géologues s'accordent sur le fait qu'après l'implication d'une croûte continentale flottante, une inversion de polarité de subduction peut se produire, les opinions divergent sur les mécanismes conduisant au changement de direction de subduction. Il n’existe donc pas de modèle unique pour représenter l’inversion de polarité de subduction. La manière dont les géologues développent les modèles dépend des paramètres qu'ils prennent en compte[1]. Certains géologues tentent de construire des modèles d'inversion de subduction par le biais d'expériences en laboratoire[2],[12],[13] ou d'observations[4],[16]. Il existe trois modèles courants : la rupture de plaque[4], la double convergence[16] et la rupture lithosphérique[2].

Les modèles de fragmentation de plaque[4] et de double convergence sont basés sur des observations de géologues[16] et le modèle de fragmentation de lithosphère est basé sur une simulation expérimentale[2].

Les critères d'une inversion de polarité de subduction sont :

  1. Le système de subduction intra-océanique avec une plaque continentale flottante
  2. Le système de subduction qui cesse avec l'implication de la plaque continentale
  3. Une vieille dalle qui se brise[2],[4]

En plus des critères pour l'occurrence de l'inversion de polarité de subduction, certains géologues ont cherché à définir les facteurs déterminants de l'initiation de ce phénomène. Les géologues et sismologues Zhang et Leng affirment que « la résistance plastique et l’âge de la plaque océanique dominante dans le système de collision arc-continent contrôlent les modes d’initiation ». Lorsque le système comprend une plaque océanique dominante avec une faible résistance plastique et une plaque océanique plus jeune, la plaque supérieure "préfère" une « inversion spontanée de polarité de subduction »[17]. En effet, en l'absence de cette résistance à la déformation, la flottabilité négative prévaut et déclenche « spontanément » la subduction, tandis qu'une résistance plastique plus forte et une plaque océanique plus ancienne favorisent une « inversion de polarité induite par la subduction ». Une résistance plastique élevée de la plaque océanique entraîne une résistance à la subduction spontanée nécessitant ainsi une inversion de polarité de subduction induite par la compression.

Les différents modèles représentant l'inversion de polarité de subduction dépendent fortement des paramètres pris en compte par les géologues. En voici un tableau récapitulatif présentant les modèles de comparaison :

Différence Rupture de dalle Double convergence Rupture lithosphérique
Raisons de rupture de dalle Force de traction sur l'ancienne dalle Glissement latéral par la nouvelle dalle Une faille préexistante entraîne la pénétration d'une nouvelle dalle
Hébergement de la nouvelle dalle Fenêtre de manteau Un profond mouvement de décrochement La pénétration de la nouvelle dalle brise l'ancienne dalle

Rupture de dalle

modifier

Le modèle de la rupture de dalle a été développé en analysant la section transversale géologique le long de la collision entre la plaque eurasienne et la plaque marine philippine, où se produit une inversion continue de la polarité de subduction[4].

Lorsque deux plaques océaniques migrent l’une vers l’autre, une plaque l’emporte sur une autre formant un système de subduction. Plus tard, une marge continentale passive légère et flottante introduite dans ce système provoquera l'arrêt du système de subduction[4]. D'une part, la plaque flottante résiste à la subduction sous la plaque dominante[4]. D’un autre côté, la plaque océanique dense au niveau de la plaque subductrice préfère se déplacer vers le bas[4]. Ces forces opposées vont générer une force de traction ou une instabilité gravitationnelle sur la dalle descendante et conduire à sa rupture[18]. La cassure produit un espace qui formera une fenêtre à dans le manteau[4]. Par la suite, la marge continentale, la moins dense, forme la plaque dominante, tandis que la plaque océanique devient la plaque subductrice[4]. La direction du système de subduction change en raison de l'espace créé par la cassure, paramètre majeur de ce modèle[4].

Le diagramme d'évolution montrant comment l'inversion de la subduction initiée par une plaque cassante au niveau de la plaque subductrice : La couleur brune correspond aux croûtes continentales les moins denses ; La couleur blanche est la croûte océanique ; 1. Deux plaques se rapprochent l'une de l'autre ; 2. La croûte continentale, flottante, résiste à la subduction ; 3. La fenêtre du manteau est créée par l'instabilité gravitationnelle ; 4. Une nouvelle plaque de subduction se développe

Modèle de double convergence

modifier

Le modèle modèle de double convergence est développé sur la base de l'évolution géologique de la subduction alpine et apennine[16].

De même, lorsque deux plaques océaniques se rapprochent, le processus de subduction peut cesser avec l'arrivée d'un bloc continental flottant. Dans ces circonstances, une nouvelle plaque peut se former au niveau de la plaque dominante en raison de la compression régionale et de la différence de densité entre le bloc continental et la plaque océanique[16]. Un prisme orogénique se forme[16]. Cependant, un problème d'espace évident se pose quant à la manière d'accueillir deux plaques tectoniques. La solution réside dans le fait que la nouvelle plaque en formation se déplace non seulement verticalement, mais aussi latéralement, entraînant ainsi un mouvement de décrochement significatif[16]. Le développement de la coexistence de deux dalles opposées est décrit comme une subduction double face[19] ou prisme doublement convergent[16]. Finalement, le développement de la nouvelle dalle grandit et glisse sur l’ancienne dalle. L'ancienne dalle se brise et le prisme orogénique s'effondre. La nouvelle dalle arrête le mouvement latéral et s'enfonce en dessous[16]. La direction du système de subduction change[16].

Evolution du modèle à double convergence : La couleur marron représente la plaque Continentale ; La couleur blanche représente la plaque océanique ; C1. La plaque contenant à la fois la plaque continentale et la plaque océanique se subduit en dessous ; 2. Le bloc continental s'engage dans la subduction avec construction d'un prisme orogénique ; 3. La nouvelle dalle se développe et deux dalles présentent un mouvement de décrochement profond (double cercle signifie pointer hors de l'écran ; Croix à l'intérieur du cercle signifie pointer vers l'écran ; 4. La nouvelle dalle descend encore davantage ; 5. La vieille dalle se brise ; 6. La nouvelle dalle s'enfonce en dessous.

Rupture de la lithosphère

modifier

Le modèle de rupture de la lithosphère est simulé par des expériences sur les hydrocarbures en laboratoire[2]. Les chercheurs ont établi un modèle de zones de subduction qui est représenté de manière analogique par des hydrocarbures de différentes densités, symbolisant ainsi différentes couches présentes dans la zone de subduction[2].

Le modèle initial de la zone de subduction simulée utilise deux pistons. Le piston lié à la plaque supérieure est fixe, tandis que celui relié à la plaque subductrice est soumis à un taux de compression constant[2]. Plus significativement, il y a un arc magmatique relativement mince et une faille préexistante plongeant vers la plaque subductrice au niveau de la plaque dominante[2]. Le détachement de la faille préexistante se produit lorsque la marge continentale flottante entre en contact avec la plaque dominante[2]. En effet, la marge flottante oppose une résistance à la subduction, ce qui augmente considérablement la force de friction dans la région de contact. La subduction est ainsi stoppée. Ensuite, une nouvelle plaque subductrice se développe au niveau de la plaque dominante sous l'effet de la compression continue[2]. La nouvelle dalle en développement finit par pénétrer et briser l'ancienne dalle. Une nouvelle zone de subduction se forme avec une polarité opposée à la précédente[2].

En réalité, l'arc magmatique est une zone relativement faible au niveau de la plaque dominante car il possède une fine lithosphère et est encore affaibli par un flux de chaleur élevé[20],[21] et un fluide chaud[22],[23]. Les failles préexistantes dans cette simulation sont également courantes dans l'arc magmatique[24]. Cette expérience est une analogie réussie avec l'inversion de polarité de subduction qui se produit au Kamtchatka au début de l'Éocène[7],[25] et l'exemple actif dans la région de Taiwan[2],[11] ainsi qu'à Timor[26],[27].

Configuration de l'expérience d'A. Chemenda du modèle de rupture de la lithosphère : la couleur blanche indique la plaque océanique (densité plus élevée) ; La couleur marron indique la plaque continentale (densité inférieure) ; la couleur verte montre la faille préexistante ; Les plaques représentées par les hydrocarbures flottent au niveau de l'asthénosphère représentée par l'eau.
Le diagramme d'évolution montrant comment l'inversion de subduction initiée par une faille préexistante au niveau de la plaque dominante. 1 : Poussée par compression ; 2 : Une nouvelle dalle se développe avec la rupture de la faille ; 3 : La nouvelle dalle pénètre ; 4 : La nouvelle dalle casse l’ancienne dalle

L'exemple de Taïwan

modifier
La carte de Taiwan montre l'emplacement de la coupe géologique et les principales subductions

Le contraste saisissant des reliefs sur l'île de Taïwan a été l'objet de nombreuses recherches. La partie nord de Taiwan compte de nombreuses plaines plates telles que la plaine d'Ilan et la plaine de Pingtung[28], tandis que dans la partie sud de Taiwan se concentrent de nombreuses hautes montagnes comme Yushan atteignant environ 3 952 m. Cette énorme différence de topographie est la conséquence d'une inversion de la polarité de la subduction[4]. La plupart des modèles étudiant ce phénomène se concentreront sur une collision active à Taïwan qui semble révéler les premiers stades d’inversion de la subduction[4],[11],[12],[13],[14],[15].

La collision de l'arc de Luçon sur une ligne plus ou moins nord-sud dans la plaque marine des Philippines (PP) avec la plaque eurasienne (EP) à tendance est-ouest a commencé au milieu du Miocène[4], formant un système de subduction intra-océanique[12],[29]. L'île s'est formée par ce processus. La différence topographique entre sud et le nord de Taïwan forme un archive de l’évolution de la zone de subduction. La plaque de la mer des Philippines se subducte sous la plaque eurasienne dans la partie sud-ouest du WEP (bord ouest de la plaque de la mer des Philippines plongeant vers le nord)[4], et cette dernière remplace la première dans la partie nord-est de ce WEP[4]. La collision entre deux plaques a commencé au nord de Taïwan et s'est propagée vers le sud avec la région plus jeune dans la partie sud. Chaque étape naissante du processus d'inversion de subduction pourrait être étudiée en corrélant des coupes transversales dans diverses régions de l'île[30].

1) Coupe transversale A-A' (post-collision)[4] : La marge continentale passive de la plaque eurasienne, une croûte continentale flottante, dépasse la plaque maritime des Philippines. La plaque eurasienne subit un étirement lithosphérique, formant la fosse d'Okinawa .

2) Coupe transversale B-B'[4] : La plaque marine philippine s'enfonce sous la plaque eurasienne, et le recul de la tranchée Ryukyu conduit à l'effondrement extensionnel du coin orogénique de Taiwan[28]. La direction de la subduction change dans la section transversale C-C'.

3) Coupe transversale C-C'[4] : Une collision drastique entre deux plaques crée un coin d'accrétion et développe une ceinture orogénique. Les orogènes de Taiwan ont atteint leur hauteur maximale avec un taux d'érosion et de croissance égal[31]. L’angle de la dalle est de près de 80 degrés vers le bas[32].

4) Coupe transversale D-D'[4] : La plaque eurasienne est activement en subduction dans la plaque de la mer des Philippines à raison de 80 mm/an le long de la fosse de Manille[28]. La plaque pénètre dans le manteau et le volume de fusion dans le coin du manteau ne cesse d'augmenter. Pendant ce temps, l'angle de la dalle de subduction n'est pas aussi raide que dans la section transversale C-C'[32]. Le coin d’accrétion vient tout juste de se développer.
5) Coupe transversale E-E'[4] (pré-Collision) : La dalle pénètre sous la plaque de la mer des Philippines et amène des matériaux hydratés pour générer un prisme du manteau[4] et l'arc volcanique de Luçon.

Articles connexes

modifier

Références

modifier
  1. a b c d e et f (en) Dennis Brown, Arc-Continent Collision, Springer, coll. « Frontiers in Earth Sciences », (ISBN 9783540885573)
  2. a b c d e f g h i j k l m n o p q et r (en) A. I. Chemenda, R. -K. Yang, J. -F. Stephan et E. A. Konstantinovskaya, « New results from physical modelling of arc–continent collision in Taiwan: evolutionary model », Tectonophysics, vol. 333, nos 1–2,‎ , p. 159–178 (DOI 10.1016/S0040-1951(00)00273-0, Bibcode 2001Tectp.333..159C)
  3. (en) S. D. Willett et C. Beaumont, « Subduction of Asian lithospheric mantle beneath Tibet inferredfrom models of continental collision », Nature, vol. 369, no 6482,‎ , p. 642–645 (DOI 10.1038/369642a0, Bibcode 1994Natur.369..642W, S2CID 4239739)
  4. a b c d e f g h i j k l m n o p q r s t u v w x y z et aa (en) Louis S. Teng, C. T. Lee, Y. B. Tsai et Li-Yuan Hsiao, « Slab breakoff as a mechanism for flipping of subduction polarity in Taiwan », Geology, vol. 28, no 2,‎ , p. 155–158 (ISSN 0091-7613, DOI 10.1130/0091-7613(2000)28<155:sbaamf>2.0.co;2, Bibcode 2000Geo....28..155T, lire en ligne)
  5. (en) P. D. Ryan et J. F. Dewey, Arc-Continent Collision, Springer Berlin Heidelberg, coll. « Frontiers in Earth Sciences », , 373–401 p. (ISBN 9783540885573, DOI 10.1007/978-3-540-88558-0_13), « Arc–Continent Collision in the Ordovician of Western Ireland: Stratigraphic, Structural and Metamorphic Evolution »
  6. a et b (en) G. Molli et J. Malavieille, « Orogenic processes and the Corsica/Apennines geodynamic evolution: insights from Taiwan », International Journal of Earth Sciences, vol. 100, no 5,‎ , p. 1207–1224 (ISSN 1437-3254, DOI 10.1007/s00531-010-0598-y, Bibcode 2011IJEaS.100.1207M, S2CID 129517282)
  7. a et b (en) E. A Konstantinovskaia, « Arc–continent collision and subduction reversal in the Cenozoic evolution of the Northwest Pacific: an example from Kamchatka (NE Russia) », Tectonophysics, vol. 333, nos 1–2,‎ , p. 75–94 (DOI 10.1016/S0040-1951(00)00268-7, Bibcode 2001Tectp.333...75K)
  8. (en) Warren Bell Hamilton, Indonesia Departemen Pertambangan et United States Agency for International Development, Tectonics of the Indonesian region, U.S. Govt. Print. Off., (lire en ligne)
  9. (en) Robert McCaffrey, Peter Molnar, Steven W. Roecker et Yoko S. Joyodiwiryo, « Microearthquake seismicity and fault plane solutions related to arc-continent collision in the Eastern Sunda Arc, Indonesia », Journal of Geophysical Research: Solid Earth, vol. 90, no B6,‎ , p. 4511–4528 (ISSN 2156-2202, DOI 10.1029/JB090iB06p04511, Bibcode 1985JGR....90.4511M, S2CID 130977948)
  10. (en) « Structure and dynamics of subducted lithosphere in the Mediterranean region », Proceedings of the Koninklijke Nederlandse Akademie van Wetenschappen, vol. 95, no 3,‎ (ISSN 0924-8323, lire en ligne)
  11. a b c et d (en) A. I. Chemenda, R. K. Yang, C. -H. Hsieh et A. L. Groholsky, « Evolutionary model for the Taiwan collision based on physical modelling », Tectonophysics, an Introduction to Active Collision in Taiwan, vol. 274, no 1,‎ , p. 253–274 (DOI 10.1016/S0040-1951(97)00025-5, Bibcode 1997Tectp.274..253C)
  12. a b c d et e (en) Peter D. Clift, Hans Schouten et Amy E. Draut, « A general model of arc-continent collision and subduction polarity reversal from Taiwan and the Irish Caledonides », Geological Society, London, Special Publications, vol. 219, no 1,‎ , p. 81–98 (ISSN 0305-8719, DOI 10.1144/GSL.SP.2003.219.01.04, Bibcode 2003GSLSP.219...81C, S2CID 130378801, lire en ligne)
  13. a b et c (en) Serge Lallemand, Yvonne Font, Harmen Bijwaard et Honn Kao, « New insights on 3-D plates interaction near Taiwan from tomography and tectonic implications », Tectonophysics, vol. 335, nos 3–4,‎ , p. 229–253 (DOI 10.1016/S0040-1951(01)00071-3, Bibcode 2001Tectp.335..229L)
  14. a et b (en) Marzieh Baes, Rob Govers et Rinus Wortel, « Switching between alternative responses of the lithosphere to continental collision », Geophysical Journal International, vol. 187, no 3,‎ , p. 1151–1174 (ISSN 0956-540X, DOI 10.1111/j.1365-246X.2011.05236.x, Bibcode 2011GeoJI.187.1151B)
  15. a et b (en) Tetsuzo Seno, « The instantaneous rotation vector of the Philippine sea plate relative to the Eurasian plate », Tectonophysics, vol. 42, no 2,‎ , p. 209–226 (DOI 10.1016/0040-1951(77)90168-8, Bibcode 1977Tectp..42..209S)
  16. a b c d e f g h i et j (en) Gianluca Vignaroli, Claudio Faccenna, Laurent Jolivet et Claudia Piromallo, « Subduction polarity reversal at the junction between the Western Alps and the Northern Apennines, Italy », Tectonophysics, vol. 450, nos 1–4,‎ , p. 34–50 (DOI 10.1016/j.tecto.2007.12.012, Bibcode 2008Tectp.450...34V)
  17. (en) Shengxing Zhang et Wei Leng, « Subduction Polarity Reversal: Induced or Spontaneous? », Geophysical Research Letters, vol. 48, no 11,‎ (ISSN 0094-8276 et 1944-8007, DOI 10.1029/2021GL093201, lire en ligne, consulté le )
  18. (en) Alexander I. Shemenda, Subduction: Insights from Physical Modeling, Springer Science & Business Media, (ISBN 9780792330424, lire en ligne)
  19. (en) Winston C. Tao et Richard J. O'connell, « Ablative subduction: A two-sided alternative to the conventional subduction model », Journal of Geophysical Research: Solid Earth, vol. 97, no B6,‎ , p. 8877–8904 (ISSN 2156-2202, DOI 10.1029/91JB02422, Bibcode 1992JGR....97.8877T)
  20. (en) Claire A. Currie et Roy D. Hyndman, « The thermal structure of subduction zone back arcs », Journal of Geophysical Research: Solid Earth, vol. 111, no B8,‎ , B08404 (ISSN 2156-2202, DOI 10.1029/2005JB004024, Bibcode 2006JGRB..111.8404C)
  21. (en) C. A Currie, K Wang, Roy D Hyndman et Jiangheng He, « The thermal effects of steady-state slab-driven mantle flow above a subducting plate: the Cascadia subduction zone and backarc », Earth and Planetary Science Letters, vol. 223, nos 1–2,‎ , p. 35–48 (DOI 10.1016/j.epsl.2004.04.020, Bibcode 2004E&PSL.223...35C)
  22. (en) D. Arcay, M.-P. Doin, E. Tric et R. Bousquet, « Overriding plate thinning in subduction zones: Localized convection induced by slab dehydration », Geochemistry, Geophysics, Geosystems, vol. 7, no 2,‎ , Q02007 (ISSN 1525-2027, DOI 10.1029/2005GC001061, Bibcode 2006GGG.....7.2007A, S2CID 135199593, lire en ligne [PDF])
  23. (en) Satoru Honda et Takeyoshi Yoshida, « Application of the model of small-scale convection under the island arc to the NE Honshu subduction zone », Geochemistry, Geophysics, Geosystems, vol. 6, no 1,‎ , Q01002 (ISSN 1525-2027, DOI 10.1029/2004GC000785, Bibcode 2005GGG.....6.1002H, S2CID 134357121)
  24. (en) John Toth et Michael Gurnis, « Dynamics of subduction initiation at preexisting fault zones », Journal of Geophysical Research: Solid Earth, vol. 103, no B8,‎ , p. 18053–18067 (ISSN 2156-2202, DOI 10.1029/98JB01076, Bibcode 1998JGR...10318053T, lire en ligne [PDF])
  25. (en) Elena A Konstantinovskaia, « Geodynamics of an Early Eocene arc–continent collision reconstructed from the Kamchatka Orogenic Belt, NE Russia », Tectonophysics, vol. 325, nos 1–2,‎ , p. 87–105 (DOI 10.1016/S0040-1951(00)00132-3, Bibcode 2000Tectp.325...87K)
  26. (en) Eli A. Silver, Donald Reed, Robert McCaffrey et Yoko Joyodiwiryo, « Back arc thrusting in the Eastern Sunda Arc, Indonesia: A consequence of arc-continent collision », Journal of Geophysical Research: Solid Earth, vol. 88, no B9,‎ , p. 7429–7448 (DOI 10.1029/JB088iB09p07429, Bibcode 1983JGR....88.7429S)
  27. (en) D. B. Snyder, H. Prasetyo, D. J. Blundell et C. J. Pigram, « A dual doubly vergent orogen in the Banda Arc continent-arc collision zone as observed on deep seismic reflection profiles », Tectonics, vol. 15, no 1,‎ , p. 34–53 (ISSN 1944-9194, DOI 10.1029/95TC02352, Bibcode 1996Tecto..15...34S)
  28. a b et c (en) Jacques Angelier, Tsui-Yü Chang, Jyr-Ching Hu et Chung-Pai Chang, « Does extrusion occur at both tips of the Taiwan collision belt? Insights from active deformation studies in the Ilan Plain and Pingtung Plain regions », Tectonophysics, geodynamics and active tectonics in East Asia, vol. 466, nos 3–4,‎ , p. 356–376 (DOI 10.1016/j.tecto.2007.11.015, Bibcode 2009Tectp.466..356A, lire en ligne)
  29. (en) P. T. Leat et R. D. Larter, « Intra-oceanic subduction systems: introduction », Geological Society, London, Special Publications, vol. 219, no 1,‎ , p. 1–17 (ISSN 0305-8719, DOI 10.1144/GSL.SP.2003.219.01.01, Bibcode 2003GSLSP.219....1L, S2CID 131046715, lire en ligne)
  30. (en) Harm J. A. Van Avendonk, Kirk D. McIntosh, Hao Kuo-Chen et Luc L. Lavier, « A lithospheric profile across northern Taiwan: from arc-continent collision to extension », Geophysical Journal International, vol. 204, no 1,‎ , p. 331–346 (ISSN 0956-540X, DOI 10.1093/gji/ggv468, Bibcode 2016GeoJI.204..331V)
  31. (en) Suppe, J. (1984), « Kinematics of arc-continent collision, flipping of subduction, and backarc spreading near Taiwan », Mem. Geol. Soc. China, no 6,‎ , p. 21–33 (lire en ligne [PDF])
  32. a et b (en) Kamil Ustaszewski, Yih-Min Wu, John Suppe et Hsin-Hua Huang, « Crust–mantle boundaries in the Taiwan–Luzon arc-continent collision system determined from local earthquake tomography and 1D models: Implications for the mode of subduction polarity reversal », Tectonophysics, geodynamics and Environment in East Asia, vol. 578,‎ , p. 31–49 (DOI 10.1016/j.tecto.2011.12.029, Bibcode 2012Tectp.578...31U)