Open AccessThree‐dimensional dynamics of hydrous thermal‐chemical plumes in oceanic subduction zones
Author(s) -
Zhu Guizhi,
Gerya Taras V.,
Yuen David A.,
Honda Satoru,
Yoshida Takeyoshi,
Connolly James A. D.
Publication year - 2009
Publication title -
geochemistry, geophysics, geosystems
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.928
H-Index - 136
ISSN - 1525-2027
DOI - 10.1029/2009gc002625
Subject(s) - geology , slab , plume , subduction , geophysics , mantle convection , mantle wedge , trench , convection , rayleigh–taylor instability , buoyancy , mantle (geology) , viscosity , mechanics , seismology , instability , tectonics , materials science , thermodynamics , physics , layer (electronics) , composite material
Hydration and partial melting along subducting slabs can trigger Rayleigh‐Taylor‐like instabilities. We use 3‐D petrological‐thermomechanical numerical simulations to investigate small‐scale convection and hydrous, partially molten, cold plumes formed in the mantle wedge in response to slab dehydration. The simulations were carried out with the I3ELVIS code, which is based on a multigrid approach combined with marker‐in‐cell methods and conservative finite difference schemes. Our numerical simulations show that three types of plumes occur above the slab‐mantle interface: (1) finger‐like plumes that form sheet‐like structure parallel to the trench, (2) ridge‐like structures perpendicular to the trench, and (3) flattened wave‐like instabilities propagating upward along the upper surface of the slab and forming zigzag patterns parallel to the trench. The viscosity of the plume material is the main factor controlling the geometry of the plumes. Our results show that lower viscosity of the partially molten rocks facilitates the Rayleigh‐Taylor‐like instabilities with small wavelengths. In particular, in low‐viscosity models (10 18 –10 19 Pa s) the typical spacing of finger‐like plumes is about 30–45 km, while in high‐viscosity models (10 20 –10 21 Pa s) plumes become rather sheet‐like, and the spacing between them increases to 70–100 km. Water released from the slab forms a low‐viscosity wedge with complex 3‐D geometries. The computed spatial and temporal pattern of melt generation intensity above the slab is compared to the distribution and ages of volcanoes in the northeast Japan. Based on the similarity of the patterns we suggest that specific clustering of volcanic activity in this region could be potentially related to the activity of thermal‐chemical plumes.