The Izu-Bonin-Mariana Subduction Factory Robert J. Stern, Geosciences Dept., U. Texas at Dallas, Box 830688, Richardson TX 75083-0688, [email protected]

The IBM arc system lies in the Western Pacific and defines the eastern margin of the Plate (Fig. 1). Because IBM is an endmember arc system in a number of ways, it is an outstanding natural laboratory for studying earth’s largest geodynamic system, the Subduction Factory. IBM is the largest intra-oceanic convergent margin, being constructed entirely within and upon oceanic lithosphere. Lithosphere produced during the ~48 million years that this Subduction Factory has operated makes up a region about the size of India, but the presently active part extends for about 300 km west of the IBM trench. IBM manifests subduction of Pacific lithosphere, with the Pacific moving NW relative to IBM, at rates that vary from about 20 mm/a south of Guam to almost 60mm/a near Japan [1]. The Philippine Sea Plate itself is moving rapidly northwestward. Subducted lithosphere varies in age from mid-Jurassic (~170 Ma) outboard of the Marianas to early Cretaceous (~130Ma) adjacent to Japan (Fig. 2). The combination of a retreating upper plate and an extremely old (hence dense) subducting plate results in a strongly extensional convergent plate boundary, particularly in the southern IBM where these effects are maximized. Off-ridge volcanism was common during mid- and Late Cretaceous time on the Pacific Plate now outboard of southern IBM. As a result, there are significant differences in the bulk composition of subducted sedimentary columns in the north and south, although the thickness of sediments is relatively constant at about 500m [2]. This modest sediment thickness is well below the ~1km thickness required for development of an accretionary prism [3] and is completely subducted. Plate convergence is oblique over most of the IBM arc, approaching pure sinistral strike-slip motion in the northern Mariana arc (Fig. 3). The Wadati-Benioff zone is well-defined beneath IBM, dipping about 45° in the north and nearly vertically in the south (Fig. 4); some of the deepest seismic activity in earth – down to 700 km - is found beneath the central Mariana arc. These and other physical constraints need to be identified and incorporated in the construction of realistic, 4-D models of mantle and fluid flow, thermal evolution, and fluid/melt generation, migration, and storage. Earth’s only T-T-T triple junction defines the northern end of IBM, where collision continues with southern Honshu at a rate of about 40mm/a. This provides an opportunity to study terrane accretion in action [4] as well as exhuming IBM arc middle crust [5]. These exposures of arc crust are correlatable with crustal structure inferred for in situ IBM crust (Fig. 5). Not only does northern IBM provide an unparalleled opportunity to examine the formation of juvenile continental crust, societal considerations compel it: continuing convergence between buoyant IBM crust and Honshu presents an imposing earthquake hazard to the greater Tokyo metropolitan area, with far-reaching implications for the global economy. An analogous opportunity to study deeper arc crust and upper mantle lies on the wall of the 11km-deep Challenger Deep at the south end of the IBM arc system. Three combined forces – retreat of the Philippine Sea Plate, subduction of unusually dense lithosphere, and oblique convergence – combine to make IBM the most strongly extended convergent margin on the planet. This evisceration provides unique opportunities to monitor the Subduction Factory. Extension in the southern IBM is oriented trench-parallel in the forearc and trench-normal in the back-arc [6]. This strongly extensional regime provides three opportunities - forearc, active arc, and back- arc basin - to sample Subduction Factory fluids and melts, more than any other convergent margin on the planet. The Mariana forearc contains the only known occurrences of subduction-related serpentinite mud volcanoes (SMV) [7] on this planet. Flows from Mariana SMV contain fragments of blueschist from the subducted slab [8] as

- 1 - well as abundant mantle fragments [9]; several SMV are actively venting slab-derived fluids, some of which support active chemosynthetic communities. Magmatic activity is expressed differently within and along the IBM arc. A true back-arc basin (BAB), with seafloor spreading, is developed only in the Mariana Trough [10], although well-developed inter-arc basins are developed in the Bonin Arc Rifted Zone farther north (Fig. 6) [11, 12]. The active magmatic arc is largely submarine but with abundant subaerial volcanoes [13, 14]; strong variations in magmatic compositions affect both submarine and subaerial edifices. BAB lavas are often aphyric or have pillow rim glass that can readily be separated, and a substantial proportion are primitive basalts; these characteristics suit BAB lavas for analysis using approaches perfected for MORB. Arc lavas are predominantly porphyritic so that bulk compositions generally do not correspond to magmatic liquids, and should not be studied using the techniques appropriate for studying aphyric or glassy samples. Accumulation of plagioclase phenocrysts in particular has led to a misperception that mafic members are dominantly high-Al basalts when in fact aphyric samples or glass inclusions in phenocrysts are tholeiites. Primitive compositions (Mg#>65) are uncommon among IBM arc lavas, so fractionation conditions and history need to be resolved. Old techniques for studying petrography should be revived and new techniques will have to be perfected if we are to understand the magmatic evolution of porphyritic arc lavas. In future studies, petrographic descriptions should be reported with major and trace element data. The need to find new ways to study evolved, porphyritic samples promises to revive traditional petrography as well as stimulate developments in quantitative petrography/image analysis and microbeam analytical techniques. These techniques will also aid investigations of abundant cumulate xenoliths, found in the lavas of several IBM volcanoes. One issue that awaits resolution is the abundance of felsic material in the IBM arc. IBM has traditionally been thought to have a basaltic bulk composition, but recent evidence from geophysics [15], exposures in the collision zone [5], glass inclusions in phenocrysts [16], and the abundance of felsic tephra in DSDP cores [17] indicates that felsic rocks comprise an important part of the IBM arc. In spite of the different eruptive styles and extent of fractionation for arc and BAB, there are strong compositional affinities between arc and BAB suites, which provide different perspectives on important controversies and enigmas. The trace element signatures of these lavas strongly manifest the ‘subduction component’: enrichments in large-ion lithophile elements and depletion in high-field strength elements, both compatible and incompatible (Fig. 7). All arc and most BAB lavas have elevated water contents [12, 18] such that it is controversial the extent to which melts are generated by decompression [19] or fluxing by hydrous fluids [20]. One abiding mystery concerns how water gets into the source of Mariana BAB when the subducted slab does not lie beneath the spreading ridge? In contrast to widespread recognition that water in IBM melts is recycled from subducted materials, the source of other elements in IBM melts is less clear. In particular, controversy continues regarding the extent to which the IBM ‘subduction component’ manifests fractionations imposed when elements are transferred from the subducted plate to the overlying mantle wedge as opposed to being developed during re-equilibration of hydrous fluids and melts with convecting mantle. U-Th disequilibria studies indicate that strong fractionation of radionuclides occurred within the last 30 kyr [21], providing timescales for fluid-mediated fractionation but where and how this occurs are unresolved. There is abundant evidence that some subducted

components7 are recycled and11 can be found34 in young207 lavas,204 especially elevated water contents, ∂ Li [22], B/Be, ∂ B [23], ∂ S [24], and Pb/ Pb [25]. These studies differently emphasize roles of subducted sediments and altered oceanic crust. Rare gas

data is lacking3 for4 arc lavas but BAB lavas contain recycled atmospheric Ar in spite of mantle-like He/ He10 [26]; similar datasets for arc lavas promise to provide important constraints. The Be signal is muted probably because subducted sediment are much older than the half-life of this isotope [27]. Other isotopic and trace element data sets do

- 2 - not readily allow recycled components to be identified: O, Sr, Nd, and Pb isotopic compositions as well as K/Rb and K/Ba are remarkably constant in arc lavas despite being extremely heterogeneous in subducted components. This homogeneity requires efficient mixing of subducted components or effective re-equilibration of ascending fluids with convecting mantle. It is critical to resolve the nature of the mantle source beneath IBM; there are strong indications from Nd, Hf, Pb, and Os data that this mantle has affinities to that beneath the Indian Ocean [28-30]. Disparate geochemical and isotopic data sets and controversial conclusions provide important constraints for understanding how the IBM Subduction Factory operates. It will be an exciting challenge for scientists with different talents and perspectives to synthesize these observations and use these to develop and test hypotheses to better understand the operation and budget of this outstanding example of the Subduction Factory. A critical aspet will be to collect and distribute representative samples to analysts and to communicate these results to modelers, and for modellers to tell analysts what kinds of information they need. Another critical effort will be the seismic imaging of the IBM Subduction Factory, with as high a resolution as possible.

1. Seno, T., S. Stein, and A.E. Gripp, A model for the motion of the Philippine Sea Plate consistent with NUVEL-1 and geological data. Journal of Geophysical Research, 1993. 98(B10): p. 17,941-17,948. 2. Plank, T. and C. Langmuir, The chemical composition of subducting sediment and its consequence for the crust and mantle. Chemical Geology, 1998. 145: p. 325-394. 3. Le Pichon, X., P. Henry, and S. Lallemant, Accretion and erosion in subduction zones: The role of fluids. Annual Reviews of Earth and Planetary Sciences, 1993. 21: p. 307-331. 4. Taira, A., S. Saito, and others, Nature and growth rate of the Northern Izu-Bonin (Ogasawara) arc crust and their implications for continental crust formation. The Island Arc, 1998. 7: p. 395-407. 5. Kawate, S. and M. Arima, Petrogenesis of the Tazawa plutonic complex, central Japan: Exposed felsic middle crust of the Izu-Bonin-Mariana arc. The Island Arc, 1998. 7: p. 342-358. 6. Wessel, J.K., et al., Extension in the northern Mariana inner forearc. Journal of Geophysical Research, 1994. 99(B8): p. 15,181-15,203. 7. Fryer, P., et al., Serpentine bodies in the Forearcs of Western Pacific Convergent Margins: Origin and Associated Fluids, in Active Margins and Marginal Basins of the Western Pacific Convergent Margins, B. Taylor and J. Natland, Editors. 1995, American Geophysical Union: Washington DC. p. 259-270. 8. Fryer, P., C.G. Wheat, and M. Mottle, Mariana blueschist mud volcanism; implications for conditions within the subduction zone. Geology, 1999. 27: p. 103- 106. 9. Parkinson, I.J. and J.A. Pearce, Peridotites from the Izu-Bonin-Mariana forearc (ODP Leg 125); evidence for mantle melting and melt-mantle interaction in a supra- subduction zone setting. Journal of Petrology, 1998. 39: p. 1577-1618. 10. Fryer, P., Geology of the Mariana Trough, in Back-arc Basins, B. Taylor, Editor. 1996. 11. Taylor, B., et al., Structural development of Sumisu , Izu-Bonin Arc. Journal of Geophysical Research, 1991. 96: p. 113-129. 12. Stolper, E. and S. Newman, The role of water in the petrogenesis of Mariana Trough magmas. Earth and Planetary Science Letters, 1994. 121: p. 293-325. 13. Bloomer, S.H., R.J. Stern, and N.C. Smoot, Physical volcanology of the submarine Mariana and Volcano Arcs. Bulletin of Volcanology, 1989. 51: p. 210-224.

- 3 - 14. Yuasa, M., et al., Submarine topography of seamounts on the volcanic front of the Izu-Ogasawara (Bonin) Arc. Bulletin of the Geological Survey of Japan, 1991. 42: p. 703-743. 15. Suyehiro, K., et al., Continental crust, crustal underplating, and low-q upper mantle beneath an oceanic island arc. Science, 1996. 272: p. 390-392. 16. Lee, J. and R.J. Stern, Glass inclusions in Mariana Arc phenocrysts; A new perspective on magmatic evolution in a typical intra-oceanic arc. Journal of Geology, 1998. 106: p. 19-33. 17. Lee, J., R.J. Stern, and S.H. Bloomer, Forty million years of magmatic evolution in the Mariana Arc: The tephra glass record. Journal of Geophysical Research, 1995. 100: p. 17,671-17,687. 18. Newman, S., E. Stolper, and R.J. Stern, H2O and CO2 in magmas from the Mariana arc and back arc system. Geochemistry, Geophysics,Geosystems, 2000. 1. 19. Plank, T. and C.H. Langmuir, An evaluation of the global variations in the major element chemistry of arc basalts. Earth and Planetary Science Letters, 1988. 90: p. 349-370. 20. Eiler, J.M., et al., Oxygen Isotope Geochemistry of Oceanic-Arc lavas. Journal of Petrology, 2000. 41: p. 229-256. 21. Elliott, T., et al., Element transport from slab to volcanic front at the Mariana arc. Journal of Geophysical Research, 1997. B102: p. 14,991-15,019. 22. Moriguti, T. and E. Nakamura, Across-arc variation of Li isotopes in lavas and implications for crust/mantle recycling at subduction zones. Earth and Planetary Science Letters, 1998. 164: p. 167-174. 23. Ishikawa, T. and F. Tera, Two isotopically distinct fluid components involved in the Mariana Arc; evidence from Nb/B ratios and B, Sr, Nd, and Pb isotope systematics. Geology27, 1999: p. 83-86. 24. Alt, J.C., W.C. Shanks III, and M.C. Jackson, Cycling of sulfur in subduction zones; the geochemistry of sulfur in the Mariana island arc and back-arc trough. Earth and Planetary Science Letters, 1993. 119: p. 477-494. 25. Woodhead, J.D. and D.G. Fraser, Pb, Sr, and 10Be isotopic studies of volcanic rocks from the Northern Mariana Islands. Implications for genesis and crustal recycling in the Western Pacific. Geochimica et Cosmochimica Acta, 1985. 49: p. 1925-1930. 26. Ikeda, Y., et al., Noble gases in pillow basalt glasses from the northern Mariana Trough back-arc basin. The Island Arc, 1998. 7: p. 471-478. 27. Tera, F., et al., Sediment incorporation in island-arc magmas: Inferences from 10Be. Geochimica et Cosmochimica Acta, 1986. 50: p. 535-550. 28. Hickey-Vargas, R., J.M. Hergt, and P. Spadea, Ocean-type isotopic signature in Western Pacific marginal basins: Origin and significance, in Active Margins and Marginal Basins of the Western Pacific, B. Taylor and J. Natland, Editors. 1995, American Geophysical Union: Washinton DC. p. 175-197. 29. Parkinson, I.J., C.J. Hawkesworth, and A.S. Cohen, Ancient mantle in a modern arc: Osmium isotopes in Izu-Bonin-Mariana forearc peridotites. Science, 1998. 281: p. 2011-2013. 30. Pearce, J.A., et al., Hf-Nd Element and Isotope Perspective on the Nature and Provenance of Mantle and Subduction Components in Western Pacific Arc-Basin systems. Journal of Petrology, 1999. 40: p. 1579-1611. 31. Nakanishi, M., K. Tamaki, and K. Kobayashi, Magnetic anomaly lineations from Late Jurassic to Early Cretaceous in the west-central Pacific Ocean. Geophysical Journal International, 1992. 109: p. 701-719.

- 4 - 32. McCaffrey, R., Estimates of modern arc-parallel strain rates in fore arcs. Geology, 1996. 24: p. 27-30. 33. Chiu, J.-M., B.L. Isaaks, and R.K. Cardwell, 3-D configuration of subducted lithosphere in the western Pacific. Geophysical Journal International, 1991. 106: p. 99-111. 34. Peate, D.W. and J.A. Pearce, Causes of spatial compositional variations in Mariana arc lavas: Trace element evidence. The Island Arc, 1998. 7: p. 479-495.

- 5 - Figure Captions: Figure 1: IBM arc system, showing extent of crust generated over the 48 Ma life of the arc (fossil plus active) as well as those components which comprise the presently active IBM Subduction Factory. Also shown is the location of the deepest place on the face of the earth, the Challenger Deep (~11km deep) and the IBM collision zone. Fig. 2: Seafloor feeding into the IBM Subduction Factory, modified after [31]. Arrows are relative velocities of the Pacific Plate with respect to the Philippine Sea Plate, in mm/a, after [1]. DSDP and ODP sites sampling units being subducted beneath IBM are shown as well. Note that the sedimentary section being subducted beneath northern IBM has fewer volcanics and volcaniclastics than that being subducted beneath the southern IBM. Fig. 3: Obliquity of convergence between the Pacific and Philippine Sea plates, as inferred from earthquake slip vectors and modified after [32]. Note that convergence is highly oblique over much of the IBM arc system. Fig. 4: Generalized topology of IBM Wadati-Benioff Zone, modified after [33]. Two perspectives are shown, with contours colored at every 100 km depth. Figure 5: Structure of IBM arc crust at 32°15’N, modified after [15]. Vertical exaggeration is about 10x. Note the crustal thickness of 22km is nearly four times that of oceanic crust but only about half that of normal continental crust. Figure 6: Along-strike profiles of the IBM arc system, from Japan (left) to Guam (right). The thick solid line shows the bathymetry and topography along the volcanic axis of the active arc, with the thin dashed horizontal line marking sea level. The approximate locations of the principal island groups (Izu, Bonin-Volcano, and Mariana) are shown. Submarine volcanoes (and the Sofugan Tectonic Line, STL) are given as italicized abbreviations: Ku, Kurose;Ms, Myojin-sho; Do, Doyo; Kk, Kaikata; Kt, Kaitoku;F, Fukutoku-oka-no-ba; HC, Hiyoshi Volcanic Complex, Nk, Nikko; Fj, Fukujin, Ch, Chamorro, D, Diamante; R, Ruby, E, Esmeralda; T; Tracy. Subaerial volcanoes are given as normal abbreviations: O, Oshima; My, Miyakejima; Mi, Mikurajima; H, Hachijojima; A, Aogashima; Su, Sumisujima, T, Torishima; Sg, Sofugan; Nishinoshima; KIJ, Kita Iwo Jima; IJ, Iwo Jima; MIJ, Minami Iwo Jima; U, Uracas; M, Maug; As, Asuncion; Ag, Agrigan; P, Pagan; Al, Alamagan; G; Guguan; S, Sarigan; An, Anatahan. Dominant compositions of arc segments are also indicated. Locations of important zones of intra- arc and back-arc extension in the north (Bonin Arc Rifted Zone) and south (Mariana Trough Back-Arc Basin) are marked. The thick dashed line shows the maximum depth in the trench along its strike. Frontal arc elements are not shown, but consist of the Bonin or Ogasawara Islands between 26° and 28°N and the Mariana frontal arc islands between 13° and 16°N. ICZ = IBM collision zone. Fig. 7: ‘Spider’ diagram for Mariana arc lavas. Elements are listed in order of increasing compatibility in mantle minerals; data for typical Maraian arc lavas is from [34]. Notice strong enrichments in LIL and depletions in HFSC, including Nb and Ta.

- 6 - IBM 40¡ ARC Pacific SYSTEM Tokyo Plate Japan* . ... IBM . Collision . . . Zone Shikoku Izu Arc Basin . 30¡

. . . Bonin Arc .

Ridge . . West Philippine Basin . 20¡

. Mariana Arc . . Kyushu .

- . . . Parece Vela . Basin . Palau Guam

10¡

.. Yap Challenger Deep

Belau 140¡

Active Spreading Ridge Extinct Spreading Ridge (Mariana Trough) Active IBM system Fossil IBM system Age of Western Pacific Seafloor

Hawaii-Emperor Smts.

Shatsky 57 Rise 195 *194 1149 * * Cretaceous 197** 196 307 *Jurassic 198* 40

I-B-M 800 * 801 Mid-PacificMid-Pacific 452v v v v v vv vv vv Mountains v v * v v v v * v v v v v v v v v 21 v v v v v v v v v v v v v 802 v v v v v vvv v * v v Caroline Ridge v v

v v462 v v v * v v v Ontong Java Jurassic v Plateau Cretaceous 140°E 160 180°

v v Oligocene and younger crust v Mid-Cretaceous Flood Basalts Oceanic Plateaus (mostly Cretaceous) DSDP & ODP sites Trenches * 60 North South Izu Bonin Mariana 30

orthogonal convergence 0 ) or rate (mm/y) -30 °

-60

sinistral 34.1° 25.2° 16.9°N 142°E Obliquity ( strike-slip -90 0 1000 2000 3000 km plate convergence obliquity (from earthquake slip vectors) arc-parallel fore-arc slip rate 128°E 138° 148°

Izu 28°N Trench

N

18° N

Mariana Trench 8°

500 km Volcanic Fore-arc Arc 0 8.5 Pacific 8.0 Plate 7.5 5 2 Trench 7.0 6.5

10 6.0 Crust 5.5 5.0 15 4.5 Depth (km) 4.0 oho 3.5 Mantle M Crust (km/sec) Velocity Seismic P-wave 20 Mantle 3.0 Moho 2.5 100 km 2.0 25 1.5 138 139 140 141 142 143 Longitude ( ¡E) Low-K Tholeiite Shoshonite Medium-K Basalt Izu Islands Bonin &Volcano Islands Mariana Islands

ICZ My Ku Ms STL IJ F HC M Ch Al An R +1 O Mi H A Su T Sg Do N KkKt KIJ MIJ Nk Fj U As Ag P G S D E T SL

-2 Active Arc -4 Bonin Arc Rifted Zone Mariana Trough Back-Arc Basin (Km) -6 500 km Trench -8

-10 34 32 30 28 26 24 22 20 18 16 14 Degrees N latitude 100

Mariana Shoshonites

10

Mariana Trough

Sample/MORB 1

Mariana Islands

Cs Ba U TaLa Pb Sr Nd Hf Eu Ti HREE RbTh K Nb Ce Pr P Zr Sm HREE