Multiscale dynamics of the Tonga—Kermadec subduction zone

Our understanding of mantle convection and the motion of plates depends intimately on our understanding of the viscosity structure of the mantle. While geoid and gravity observations have provided fundamental constraints on the radial viscosity structure of the mantle, the influence of short-wavelength variations in viscosity is still poorly understood. We present 2-D and 3-D finite-element models of mantle flow, including strong lateral viscosity variations and local sources of buoyancy, owing to both thermal and compositional effects. We first use generic 2-D models of a subduction zone to investigate how different observations depend on various aspects of the viscosity structure, in particular, the slab and lower-mantle viscosity and the presence of a low-viscosity region in the mantle wedge above the slab. We find that: (1) the strain rate provides a strong constraint on the absolute viscosity of the slab (10^(23) Pa s); (2) stress orientation within the slab is sensitive to the relative viscosity of the slab, lower mantle and the wedge; and (3) the stress state and topography of the overriding plate depend on the wedge viscosity and local sources of buoyancy. In particular, the state of stress in the overriding plate changes from compression to extension with the addition of a low-viscosity wedge. We then use observations of strain rate, stress orientation, dynamic topography and the geoid for the Tonga–Kermadec subduction zone as simultaneous constraints on the viscosity and buoyancy in a 3-D regional dynamic model. Together these observations are used to develop a self-consistent model of the viscosity and buoyancy by taking advantage of the sensitivity of each observation to different aspects of the dynamics, over a broad range of length-scales. The presence of a low-viscosity wedge makes it possible to match observations of shallow dynamic topography and horizontal extension within the backarc, and down-dip compression in the shallow portion of the slab. These results suggest that a low-viscosity wedge plays an important role in controlling the presence of backarc spreading. However, for a model with a low-viscosity and low-density region that provides a good fit to the observed topography, we find that a reduction of the slab density by a factor of 1.3 relative to the reference density model, is required to match the observed geoid. These results suggest that compensation of the slab by dynamic topography may be a much smaller effect at short to intermediate wavelengths than predicted by long-wavelength modelling of the geoid.

[1]  Robert W. Clayton,et al.  Constraints on the Structure of Mantle Convection Using Seismic Observations, Flow Models, and the Geoid , 1989 .

[2]  J. Royer,et al.  Revised plate motions relative to the hotspots from combined Atlantic and Indian Ocean hotspot tracks , 1993 .

[3]  Gregory M. Hulbert,et al.  Accurate determination of surface normal stress in viscous flow from a consistent boundary flux method , 1993 .

[4]  J. Jackson,et al.  Distributed deformation in the subducting lithosphere at Tonga , 1996 .

[5]  Louis Moresi,et al.  Role of faults, nonlinear rheology, and viscosity structure in generating plates from instantaneous mantle flow models , 1998 .

[6]  M. Bevis.  Seismic Slip and Down-Dip Strain Rates in Wadati-Benioff Zones , 1988, Science.

[7]  F. Richter.  On the driving mechanism of plate tectonics , 1977 .

[8]  Akira Hasegawa,et al.  Tomographic imaging of P and S wave velocity structure beneath northeastern Japan , 1992 .

[9]  H. Fukuyama,et al.  Generation of arc basalt magmas and thermal structure of the mantle wedge in subduction zones , 1983 .

[10]  D. McAdoo.  On the compensation of geoid anomalies due to subducting slabs , 1981 .

[11]  D. Stevenson,et al.  Physical model of source region of subduction zone volcanics , 1992 .

[12]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

[13]  S. Karato,et al.  Water, partial melting and the origin of the seismic low velocity and high attenuation zone in the upper mantle , 1998 .

[14]  D. Rubie,et al.  Experimental constraints on the depth of olivine metastability in subducting lithosphere , 2001 .

[15]  Shijie Zhong,et al.  Interaction of weak faults and non-newtonian rheology produces plate tectonics in a 3D model of mantle flow , 1996, Nature.

[16]  S. Karato.  Mapping water content in the upper mantle , 2013 .

[17]  B. Evans,et al.  Strength of the lithosphere: Constraints imposed by laboratory experiments , 1995 .

[18]  Magali I. Billen,et al.  A low viscosity wedge in subduction zones , 2001 .

[19]  A. Hasegawa,et al.  Deep structure of the northeastern Japan arc and its relationship to seismic and volcanic activity , 1991, Nature.

[20]  Variations in P wave speeds and outboard earthquakes: Evidence for a petrologic anomaly in the mantle transition zone , 2000 .

[21]  B. Hager,et al.  Effects of lateral viscosity variations on long-wavelength geoid anomalies and topography , 1989 .

[22]  H. Keppler,et al.  Solubility of water in the α, β and γ phases of (Mg,Fe)2SiO4 , 1996 .

[23]  D. Gubbins,et al.  DEFORMATION OF SUBDUCTED OCEANIC LITHOSPHERE , 1997 .

[24]  C. Scholz,et al.  On the mechanism of seismic decoupling and back arc spreading at subduction zones , 1995 .

[25]  A. E. Ringwood,et al.  Phase Transformations and Differentiation in Subducted Lithosphere: Implications for Mantle Dynamics, Basalt Petrogenesis, and Crustal Evolution , 1982, The Journal of Geology.

[26]  M. Brudzinski,et al.  Evidence for a large-scale remnant of subducted lithosphere beneath Fiji. , 2001, Science.

[27]  E. Boschi,et al.  Plate motion and dragging of the upper mantle: Lateral variations of lithospheric thickness and their implications for intraplate deformation , 1992 .

[28]  Bradford H. Hager,et al.  A simple global model of plate dynamics and mantle convection , 1981 .

[29]  C. Eloy,et al.  Free-surface formulation of mantle convection—II. Implication for subduction-zone observables , 1996 .

[30]  C. G. Chase.  Extension behind island arcs and motions relative to hot spots , 1978 .

[31]  T. Seno,et al.  Arc stresses determined by slabs: Implications for mechanisms of back‐arc spreading , 1998 .

[32]  H. Schmeling,et al.  On the effects of the lithosphere on mantle convection and evolution , 1982 .

[33]  Louis Moresi,et al.  Role of temperature‐dependent viscosity and surface plates in spherical shell models of mantle convection , 2000 .

[34]  Louis Moresi,et al.  The accuracy of finite element solutions of Stokes's flow with strongly varying viscosity , 1996 .

[35]  J. Morris,et al.  The subducted component in island arc lavas: constraints from Be isotopes and B–Be systematics , 1990, Nature.

[36]  D. Wiens,et al.  Seismic attenuation tomography of the Tonga‐Fiji region using phase pair methods , 1999 .

[37]  J. Z. Zhu,et al.  The finite element method , 1977 .

[38]  B. Hager,et al.  Effects of plate bending and fault strength at subduction zones on plate dynamics , 1999 .

[39]  Shijie Zhong,et al.  Analytic solutions for Stokes' flow with lateral variations in viscosity , 1996 .

[40]  S. Buiter,et al.  A modelling study of vertical surface displacements at convergent plate margins , 2001 .

[41]  B. Isacks,et al.  Lateral variations of seismic-wave attenuation in the upper mantle above the inclined earthquake zone of the Tonga Island Arc: Deep anomaly in the upper mantle , 1971 .

[42]  S. Stein,et al.  Thermo‐Mechanical Evolution of Oceanic Lithosphere: Implications for the Subduction Process and Deep Earthquakes , 1996 .

[43]  Sean C. Solomon,et al.  Localization of gravity and topography: constraints on the tectonics and mantle dynamics of Venus , 1997 .

[44]  S. Peacock.  Blueschist facies metamorphism, shear heating, and P-T-t paths in subduction shear zones , 1992 .

[45]  Y. Ogawa,et al.  Tension cracks on the oceanward slopes of the northern Japan and Mariana Trenches , 1997 .

[46]  J. G. Caldwell,et al.  An elastic-perfectly plastic analysis of the bending of the lithosphere at a trench , 1978 .

[47]  Ulrich R. Christensen,et al.  Some effects of lateral viscosity variations on geoid and surface velocities induced by density anomalies in the mantle , 1993 .

[48]  S. Newman,et al.  The role of water in the petrogenesis of Mariana trough magmas , 1994 .

[49]  J. D. Bremaecker.  Is the oceanic lithosphere elastic or viscous , 1977 .

[50]  B. Hager,et al.  Geoid Anomalies in a Dynamic Earth , 1984 .

[51]  D. Koch,et al.  The effect of lateral viscosity variations on surface observables , 1989 .

[52]  Kelin Wang,et al.  Seismic consequences of warm versus cool subduction metamorphism: examples from southwest and northeast japan , 1999, Science.

[53]  Shijie Zhong,et al.  Controls on trench topography from dynamic models of subducted slabs , 1994 .

[54]  A. Watts,et al.  On lithospheric flexure seaward of the Bonin and Mariana trenches , 1979 .

[55]  J. Stock,et al.  Morphology and origin of the Osbourn Trough , 2000 .

[56]  Bijaya B. Karki,et al.  Origin of lateral variation of seismic wave velocities and density in the deep mantle , 2001 .

[57]  Louis Moresi,et al.  Free-surface formulation of mantle convection—I. Basic theory and application to plumes , 1996 .

[58]  D. L. Anderson,et al.  Present-day plate motion constraint on mantle rheology , 1997 .

[59]  P. Kelemen,et al.  A review of melt migration processes in the adiabatically upwelling mantle beneath oceanic spreading ridges , 1997, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[60]  P. Molnar,et al.  Distribution of stresses in the descending lithosphere from a global survey of focal‐mechanism solutions of mantle earthquakes , 1971 .

[61]  B. Hager,et al.  Long-wavelength variations in Earth’s geoid: physical models and dynamical implications , 1989, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[62]  Bradford H. Hager,et al.  The distribution of earthquakes with depth and stress in subducting slabs , 1984 .

[63]  S. Karato,et al.  Effects of pressure on high-temperature dislocation creep in olivine , 2003 .

[64]  D. Rubie.  Solubility of water in the a , b and g phases of ( Mg , Fe ) 2 SiO 4 , 2022 .

[65]  C. Frohlich,et al.  The relationship between Wadati‐Benioff Zone geometry and P, T and B axes of intermediate and deep focus earthquakes , 1987 .

[66]  Greg Hirth,et al.  Water in the oceanic upper mantle: implications for rheology , 1996 .

[67]  M. Gurnis,et al.  Viscous flow model of a subduction zone with a faulted lithosphere: Long and short wavelength topography, gravity and geoid , 1992 .

[68]  M. Richards,et al.  The geoid constraint in global geodynamics: viscosity structure, mantle heterogeneity models and boundary conditions , 1997 .

[69]  Douglas A. Wiens,et al.  An empirical relationship between seismic attenuation and velocity anomalies in the upper mantle , 2000 .

[70]  M. Ravine,et al.  Geoid effects of lateral viscosity variation near the top of the mantle: A 2D model , 1993 .

[71]  R. Hyndman,et al.  An inverted continental Moho and serpentinization of the forearc mantle , 2002, Nature.

[72]  B. Hager.  Subducted slabs and the geoid: Constraints on mantle rheology and flow , 1983 .

[73]  Walter H. F. Smith,et al.  New version of the generic mapping tools , 1995 .

[74]  John A. Hildebrand,et al.  Depth Extent of the Lau Back-Arc Spreading Center and Its Relation to Subduction Processes , 1997 .

[75]  W. Holt.  Flow fields within the Tonga Slab determined from the moment tensors of deep earthquakes , 1995 .

[76]  Robert W. Clayton,et al.  Lower mantle heterogeneity, dynamic topography and the geoid , 1985, Nature.

[77]  R. V. D. Hilst.  Complex morphology of subducted lithosphere in the mantle beneath the Tonga trench , 1995, Nature.

[78]  H. Melosh.  Dynamic support of the outer rise , 1978 .

[79]  M. Gurnis,et al.  Constraints on the lateral strength of slabs from three-dimensional dynamic flow models , 1996 .

[80]  Kenneth E. Torrance,et al.  Thermal convection with large viscosity variations , 1971, Journal of Fluid Mechanics.

[81]  Shijie Zhong,et al.  Effects of plate and slab viscosities on the geoid , 1999 .

[82]  J. D. Zund,et al.  Geophysical Geodesy: The Slow Deformation of the Earth , 1991 .

[83]  R. Huene,et al.  Initial reports of the deep sea drilling project: National Science Foundation, Washington, D.C., 1969, 672 pp., U.S. $ 10.25 , 1971 .

[84]  Louis Moresi,et al.  Numerical investigation of 2D convection with extremely large viscosity variations , 1995 .

[85]  C. Froidevaux,et al.  Global plate motion and the geoid: a physical model , 1988 .