* Your assessment is very important for improving the work of artificial intelligence, which forms the content of this project
Download Geology
Oceanic trench wikipedia , lookup
Abyssal plain wikipedia , lookup
Cimmeria (continent) wikipedia , lookup
Supercontinent wikipedia , lookup
Great Lakes tectonic zone wikipedia , lookup
Mantle plume wikipedia , lookup
Cascade Volcanoes wikipedia , lookup
Algoman orogeny wikipedia , lookup
Baltic Shield wikipedia , lookup
Downloaded from geology.gsapubs.org on February 2, 2012 Geology New seismological constraints on growth of continental crust in the Izu-Bonin intra-oceanic arc Shuichi Kodaira, Takeshi Sato, Narumi Takahashi, Seiichi Miura, Yoshihiko Tamura, Yoshiyuki Tatsumi and Yoshiyuki Kaneda Geology 2007;35;1031-1034 doi: 10.1130/G23901A.1 Email alerting services click www.gsapubs.org/cgi/alerts to receive free e-mail alerts when new articles cite this article Subscribe click www.gsapubs.org/subscriptions/ to subscribe to Geology Permission request click http://www.geosociety.org/pubs/copyrt.htm#gsa to contact GSA Copyright not claimed on content prepared wholly by U.S. government employees within scope of their employment. Individual scientists are hereby granted permission, without fees or further requests to GSA, to use a single figure, a single table, and/or a brief paragraph of text in subsequent works and to make unlimited copies of items in GSA's journals for noncommercial use in classrooms to further education and science. This file may not be posted to any Web site, but authors may post the abstracts only of their articles on their own or their organization's Web site providing the posting includes a reference to the article's full citation. GSA provides this and other forums for the presentation of diverse opinions and positions by scientists worldwide, regardless of their race, citizenship, gender, religion, or political viewpoint. Opinions presented in this publication do not reflect official positions of the Society. Notes The Geological Society of America, Inc. Downloaded from geology.gsapubs.org on February 2, 2012 New seismological constraints on growth of continental crust in the Izu-Bonin intra-oceanic arc Shuichi Kodaira Takeshi Sato Narumi Takahashi Seiichi Miura Yoshihiko Tamura Yoshiyuki Tatsumi Yoshiyuki Kaneda Institute for Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technology, Showa-machi 3173-25, Kanazawa-ku, Yokohama, Japan Institute for Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technology, Natsushima 2-15, Yokosuka, Japan ABSTRACT The process by which continental crust has formed is not well understood, though such crust mostly forms at convergent plate margins today. It is thus imperative to study modern intra-oceanic arcs, such as those common in the western Pacific Ocean. New seismic studies along the representative Izu-Bonin intra-oceanic arc provide unique along-strike images of arc crust and uppermost mantle to complement earlier, cross-arc lithospheric profiles. These reveal two scales (1000–10 km scale) of variations, one at the scale of the Izu versus Bonin (thick versus thin) arc crust, the other at the intervolcano (~50 km) scale. These images show that: (1) the bulk composition of the Izu-Bonin arc crust is more mafic than typical continental crust, (2) the middle crust with seismic velocities similar to continental crust is predominantly beneath basaltic arc volcanoes, (3) the bulk composition beneath basaltic volcanoes changes little at thick and thin arc segments, and (4) a process to return lower crustal components to the mantle, such as delamination, is required for an arc crust to evolve into continental crust. Continued thickening of the Izu-Bonin crust, accompanied by delamination of lowermost crust, can yield velocity structure of typical continental crust. Keywords: intra-oceanic arc, continental crust, crustal evolution, seismic imaging, Izu-Bonin arc. INTRODUCTION The growth of continents is an important geoscientific issue. It appears that crust of this composition was accomplished by the episodic generation of intermediate or andesitic magmas, perhaps at intraoceanic arcs (Taylor, 1967; Rudnick, 1995; Rudnick and Fountain, 1995; Suyehiro et al., 1996), but that is questionable. The most problematic issue in understanding the evolution of continental crust is that presentday intra-oceanic arcs are dominated by basaltic magmatism, but the bulk composition of continental crust is andesitic (Arculus, 1981; Pearcy et al., 1990). Several petrological models have been proposed to address this paradox (Rudnick, 1995; Tatsumi, 2005), for example: (1) direct generation of andestic magma in the mantle by the reaction of slab-derived melts with mantle peridotites (e.g., Kelemen, 1995; Rapp et al., 1999) or by hydrous melting of mantle peridotite (e.g., Kushiro, 1974), and (2) fractional crystallization of basaltic magmas and/or remelting of basaltic (or mafic) initial arc crust and subsequent delamination of mafic lower crustal material (e.g., Kay and Kay, 1993; Jull and Kelemen, 2001; Tatsumi and Kogiso, 2003). The first model predicts the direct production of crust with both andesitic and basaltic components throughout the crustal evolution process (e.g., Tamura and Tatsumi, 2002), while the second model predicts a mafic bulk composition of crust before delamination, and that the composition is probably even more mafic in the initial stage of crustal evolution. Delamination is also required in the second model because even the mantle-derived magnesian andesites are more mafic than continental crust. Until arcs are deeply drilled, only geophysical investigations can constrain the composition and structure of arc crust. We can infer the bulk composition of crust if we have good information about its seismic velocity structure (e.g., Smithson et al., 1981; Kelemen and Holbrook, 1995). The temperature and composition of crust affect seismic velocities and thus calculations of average velocity of the crust. For example, seismic imaging, laboratory experiments, and thermodynamic modeling have shown that the seismic velocities at constant temperature depend primarily on the MgO and SiO2 content. The seismic velocity of bulk crust increases with increasing MgO content and decreasing SiO2 content (Kelemen and Holbrook, 1995; Behn and Kelemen, 2006), but addition of other oxides, such as FeO, CaO, and Al2O3, did not substantially affect seismic velocity, Vp, (Kelemen and Holbrook, 1995). This means that crust of higher bulk seismic velocity than typical continental crust has a more mafic or basaltic bulk composition (i.e., higher MgO content) than typical continental crust, and that crust of slower bulk seismic velocity is more silicic or andestic (i.e., higher SiO2 content). Temperature also influences seismic velocities in island arcs. However, we concluded that for the purpose of this study, the effect of temperature on velocity (i.e., ~0.2 km/s) is insignificant. Temperature (T ) variations of several hundreds of degrees Celsius are required to obtain a velocity variation of ~0.2 km/s if an average value of the temperature derivative for crustal rocks is used (dVp/dT = –4 × 10−4 km/s/°C) (Rudnick, 1995). However, even larger temperature variations (>1 × 103 °C) are required to explain the velocity variations discussed in this study if we apply a temperature derivative measured directly from rocks sampled at the Izu collision zone [dVp/dT = (–2.7 to –5.3) × 10−5 km/s/°C] (Kitamura et al., 2003). Temperature variations of this magnitude are unprecedented in crustal rocks. In this paper, we show new seismologically derived constraints on the evolution process of continental crust. For our study, we selected the Izu-Bonin arc, where the Pacific plate subducts beneath the Philippine Sea plate (Fig. 1). An advantage of the Izu-Bonin arc is that the effect of preexisting continental crust has been minimal since subduction of the Pacific plate initiated ca. 48 Ma (Stern et al., 2003). Moreover, different stages of a crustal evolution have been proposed for the Izu and Bonin parts of the arc system, on the basis of submarine topography, rock chemistry, distribution of hypocenters of deep earthquakes, and tectonic history (Stern et al., 2003; Yuasa, 1991). If we accept a tectonic scenario where the spreading center of the Parece Vela Basin once propagated to the south of the Izu arc ca. 25 Ma, an earlier stage of crustal evolution may be expected at the Bonin arc than at the Izu arc. This scenario proposes the generation of a wedge-shaped oceanic crust to the south of the Izu arc in the area of the present-day Bonin arc (Yuasa, 1991). In addition to the difference of the tectonic history between the Izu and the Bonin arc, a difference of magma production rate might be another factor controlling variation of crustal volume along the arc. DATA ACQUISITION AND PROCESSING We conducted two active-source wide-angle seismic experiments (Fig. 1) to image the structure immediately beneath the volcanic fronts of the Izu arc (in 2004) (Kodaira et al., 2007) and the Bonin arc (in 2005). © 2007 The Geological Society of America. For permission to copy, contact Copyright Permissions, GSA, or editing@geosociety.org. GEOLOGY, November 2007 Geology, November 2007; v. 35; no. 11; p. 1031–1034; doi: 10.1130/G23901A.1; 4 figures; Data Repository item 2007252. 1031 nT ren ch Downloaded from geology.gsapubs.org on February 2, 2012 35° Myk Mkr c Izu ar Krs Hcj Shc Ags Ssc Tsm Sfg G Nsi Bonin Kkt Ktk arc Kij Pacific Plate ar Parece Vela Basin na ia ar M Philippine Sea Plate c 135° 140° 145° Figure 1. Location map and topographic map of Izu-Bonin arc. Black and blue lines indicate wide-angle seismic profiles in 2004 and 2005, respectively. Os—Oh-shima; Nij—Nii-jima; Myk—Miyake-jima; Mkr— Mikura-jima; Krs—Kurose hole; Hcj—Hachijo-jima; Shc—South Hachijo caldera; Ags—Aoga-shima; Myn—Myojin knoll; Sms— South Sumisu; Ssc—South Sumisu caldera; Tsm—Torishima; Sfg— Sofugan; G—Getsuyo seamount; Ka—Kayo seamount; S—Suiyo seamount; Kn—Kinyo seamount; D—Doyo seamount; Nsi—Nishinoshima; Kkt—Kaikata seamount; Ktk—Kaitoku seamount; Kij—Kita Iou-jima; SFG-TL—Sofugan tectonic line. In the both experiments, a linear array of densely deployed (~5 km spacing) ocean-bottom seismographs (OBSs) and a large air-gun array (~197 L) were used. Our previous study (Kodaira et al., 2007) only described a seismic model obtained from the 2004 experiment, while this study processed the data from the 2005 experiment (Fig. 2; GSA Data Repository Appendix, Fig. DR11) and compiled structural models imaged from both data sets. Consequently, compiling the two data sets provides a deep seismic transect along the total of a 1050-km-long profile. We modeled the wide-angle seismic data in two ways: to image seismic velocity we used seismic refraction arrival tomography (Zhang and Toksöz, 1998), and to image seismic reflectivity we used a diffractionstack migration approach using picked reflection traveltimes (Fujie et al., 2006). Details of the modeling procedure are provided in the Appendix (Figs. DR1–DR3). In order to estimate a resolution of the model, we performed a checkerboard resolution test following a procedure described by Zelt and Barton (1998). The results of the test (Kodaira et al., 2007; Fig. DR4) show that the structure shallower than the dashed line in Figure 3A is well resolved by our data set. 1 GSA Data Repository item 2007252, Appendix, (modeling procedure, including reflectivity imaging, traveltime fitting, and resolution of the model), is available online at www.geosociety.org/pubs/ft2007.htm, or on request from editing@geosociety.org or Documents Secretary, GSA, P.O. Box 9140, Boulder, CO 80301, USA. 1032 S 8 4 Crus tal r efra ction 160 120 80 40 Offset (km) 00 ion act refr stal u r C 40 Mantle refraction refraction Mantle 80 120 Offset (km) 160 Ridge e 25° SF G-T L idg North American Plate Kn D Bonin uR ala S OBS50 Bonin arc S N Figure 2. Examples of observed record sections from OBS31 (oceanbottom seismometer) in the Izu profile and OBS50 in the Bonin profile. Because apparent velocities of the first arrivals were strongly affected by seafloor topography, we applied static corrections for the water layer in plotting data. Locations of OBSs are shown in Figure 3A. Reduction velocity used was 8 km/s. Horizontal axis indicates offset distance from OBS. A 5–15 Hz band-pass filter and 2 s automatic gain control were applied. See Fig. DR1 (see footnote 1) for more examples of data. nch -P Ka Eurasian Plate OBS31 Izu arc N 6 0 Myn Sms nin Tre hu us Ky 30° 10 2 Izu-Bo Shikoku Basin 2005 Jap a Os Nij 2004 T - D/8 (s) Japan SEISMIC IMAGE AND DISCUSSION The seismic data revealed marked structural differences between the Izu and Bonin arcs (Fig. 2). In the Izu section (left panel in Fig. 2), the observed range of the crustal refraction phase, which reflects the thickness of crust, is more than three times as long as that of the Bonin section (right panel in Fig. 2) (Fig. DR5). This difference indicates a considerable increase of crustal volume toward the Izu arc. The final seismic image using those data shows that the largest volume of crust (~32 km thick) is under the Izu arc beneath Hachijo-jima, and that the smallest volume of crust (~10 km thickness) is between the Kayo and Suiyo seamounts in the Bonin arc, where the major tectonic line (the Sofugan tectonic line) meets the volcanic front (Fig. 3). The 10-km-thick crust represents the thinnest arc crust among reported crustal structure of arcs on the Earth. The most important finding is the variation of the volume of the middle crust (Vp = 6.0–6.8 km/s). We interpreted the middle crust to be composed of plutonic rocks of felsic (Vp = 6.0–6.5 km/s) to intermediate (Vp = 6.5–6.8 km/s) composition on the basis of comparison to sonic-wave velocities of exposed rocks sampled at the collision zone between the Izu arc and central Japan (Kitamura et al., 2003). The felsic to intermediate composition is believed to be evidence of the continental composition of crust (Rudnick, 1995). Comparing the chemical composition of volcanic rocks (Bloomer et al., 1989; Yuasa and Nohara, 1992; Kodaira, et al., 2007) with the variation in volume of the middle crust (Figs. 3B, 3C) clearly demonstrates that the volume of the middle crust increases toward the centers of each of the basaltic volcanoes (red dots in Fig. 3B). Although the volume of middle crust in the Bonin arc is one-third of that in the Izu arc, the above structural characteristics (i.e., thickening the middle crust toward the centers of basalt volcanoes) are observed along the entire arc system (Fig. 3B). From these observations we conclude that the felsic to intermediate component of the crust was created beneath the basaltic volcanoes, and that this process took place even for the thinner arc crust. The model shows seismic velocities in the range 6.8–7.6 km/s underlying the middle crust. The seismic velocity of the upper half of this component (Vp = 6.8–7.2 km/s) corresponds to the velocity of gabbroic plutons exposed in the collision zone (Kitamura et al., 2003). The lower part of this component (Vp = 7.2– 7.6 km/s) is interpreted to represent mafic to ultramafic cumulates (Miller and Christensen, 1994; Kodaira et al., 2007), but exposures of these rocks have not been identified in the collision zone (Kitamura et al., 2003). To qualitatively evaluate the variation of the structure and composition of the crust, we calculated the average crustal velocity (Fig. 3B) using the seismic velocity between the top of the middle crust and the bottom of the mafic to ultramafic cumulate component. The crustal material with a velocity slower than 6.0 km/s (the upper crust) was excluded because upper GEOLOGY, November 2007 Downloaded from geology.gsapubs.org on February 2, 2012 Figure 3. A: Seismic velocity image along volcanic front from Izu to Bonin arc obtained by seismic refraction tomography. Seismic image of Izu section was reported on by Kodaira et al. (2007). Modeling procedure, including reflectivity imaging, traveltime fitting, and resolution of the model, is shown in the Appendix, Figures DR2–DR4 (see footnote 1). Checkerboard test (Fig. DR3) shows that the structure shallower than the dashed lines is well resolved. Layers A–E indicate geological interpretations of seismic image: A—Upper crust consisting of sediment, volcaniclastics, and volcanic rocks. B—Felsic composition plutons. C—Intermediate composition plutons. D—Mafic plutons. E—Mafic to ultramafic cumulates. F—Upper mantle. See also the Appendix. B: Average crustal seismic velocity (black line) and thickness of the middle crust (Vp = 6.0–6.8 km/s) (red line), which is interpreted to be plutonic rocks of felsic to intermediate composition. Black and red dots indicate average seismic velocities and thicknesses of the middle crust, respectively, beneath basaltic volcanoes. Blue dots show average crustal seismic velocities beneath basaltic volcanoes, but excluding the Vp = 7.2–7.6 km/s component. Orange shading shows the velocity range of typical continental crust (Christensen and Mooney, 1995). C: Average wt % SiO2 of volcanic rocks sampled and dredged from Quaternary volcanoes (Bloomer et al., 1989; Yuasa and Nohara, 1992; Kodaira, et al., 2007). Abbreviations as in Figure 1. crustal velocities are considered to be strongly affected by several parameters other than crustal composition, such as variable fracture distribution and porosity (e.g., Carlson and Gangi, 1985; Kelemen and Holbrook, 1995). We recognized the variation of average seismic velocity along the arc, which correlates well with the volume variation of the middle crust. However, the average seismic velocities beneath each basaltic volcano do not vary (Vp = ~6.8 km/s) from the thick Izu arc to the thin Bonin arc (black dots in Fig. 3B). This means that the volume ratios of each crustal component are equivalent in the both thick Izu arc and the thin Bonin arc beneath the basaltic volcanoes (Fig. 4A). It is important to note that those velocities are remarkably higher than the average velocity of the typical continental crust (Vp = 6.4 ± 0.21 km/s) (Christensen and Mooney, 1995). These observations provide two strong constraints on the growth process of continental crust: (1) the bulk chemical composition of the crust beneath the basaltic volcanoes is the same for thick and thin arc crust, and (2) even though felsic to intermediate crust has been formed beneath the Izu-Bonin arc, the bulk chemical composition of the crust beneath the basaltic volcanoes is still more mafic than that of typical continental crust. This latter observation suggests that to transform arc crust into continental crust, there must be a process to return the mafic to ultramafic cumulates to the mantle, such as delamination (Kay and Kay, 1993), foundering (Oliver et al., 2003), or transformation (Takahashi et al., 2007). The requirement for such a process is well demonstrated by the calcu- GEOLOGY, November 2007 lated average seismic velocities excluding the mafic to ultramafic cumulates (i.e., layer E in Fig. 3B; Vp = 7.2–7.6 km/s). The average velocities derived this way beneath the basaltic volcanoes (blue dots in Fig. 3B) are very close to the range of velocities for typical continental crust, with the exception of the two volcanoes (Nishino-shima and Kaikata seamount). This implies that continued thickening of the Izu-Bonin arc crust, accompanied by delamination of lowermost crust, can yield typical continental crust velocity structure. Vertical extension of the velocity-depth (V-D) profile (which represents crustal growth while maintaining constant volume ratios for each crustal component) (Fig. 4) also supports the proposed process. A 250% vertical extension of the V-D profile beneath the middle crust at the Suiyo seamount of the Bonin arc shows a similar pattern to that of the Izu arc at Aoga-shima. In addition, the 150% vertical extension of the V-D profile at the Izu arc agrees well with a typical continental crust (Christensen and Mooney, 1995; Rudnick and Fountain, 1995), except for ~10 km of the mafic to ultramafic cumulates layer. Another important outcome of this study is the finding of a unique structure beneath the rhyolite volcanoes that is predominantly observed between the large basalt volcanoes in the Izu arc among the Izu-Bonin arc. Peak to peak values of the variation curve of the average velocity in the Izu arc are larger (e.g., 0.25 km/s between Aoga-shima and Myojinknoll) (Fig. 3B) than those of the Bonin arc (e.g., 0.12 km/s near Kayo seamount). This difference is mainly attributed by the difference of vol- 1033 Downloaded from geology.gsapubs.org on February 2, 2012 Shields and Platforms B Erosion Izu arc Depth (km) 0 10 250 % 30 20 30 (7.1±0.2) 50 Delamination Continental crust C&M R&F 10 (6.3±0.2) 20 (6.6±0.1)150 % 40 0 Bonin arc Depth (km) A 40 1.8-5.8 Sediment, volcaniclastics, volcanic rocks 6.0-6.5 Felsic plutons 6.5-6.8 Intermediate composition plutons 6.8-7.2 Gabbroic plutons 7.2-7.6 Mafic to ultramafic cumulates layer Continental Moho 50 4 8 6 7 5 Velocity (km/s) Aoga-shiima Suiyo-smt. Vertical extension 150 % 250 % Figure 4. A: Schematic one-dimensional (1-D) structure and seismic velocities of Bonin arc, Izu arc, and predicted continental crust (left). Structure in red frame of left column represents typical continental structure for shields and platforms (Rudnick and Fountain, 1995). For thick arc crust to evolve to continental crust, a process is required to return the component interpreted as mafic to ultramafic cumulates (Vp = 7.2–7.6) to the mantle. The upper crustal component consisting of sediment, volcaniclastics, and volcanic rocks may be eroded during the evolution. B: 1-D velocity-depth profile at Aoga-shima in the Izu arc and at the Suiyo seamount in the Bonin arc. Dashed blue and red lines show 150% and 250% vertically extended profiles beneath the middle crust at Aoga-shima and the Suiyo seamount, respectively. Velocity-depth profiles of typical continental structures compiled by Christensen and Mooney (C&M, 1995) and Rudnick and Fountain (R&F, 1995) are superimposed. ume ratios of the middle crust, i.e., smaller volume ratios of the middle crust beneath the rhyolite volcanoes. Those observations imply that the bulk crustal composition beneath the rhyolite volcanoes in the Izu arc is even more mafic. Although seismic data do not constrain a cause of the rhyolite volcanism, our seismic image may suggest that there is a process for making more mafic crust beneath the rhyolite volcanoes due to reduction of volume ratio of the middle crust. ACKNOWLEDGMENTS This study was funded by the Institute for Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technology. REFERENCES CITED Arculus, R.J., 1981, Island arc magmatism in relation to the evolution of the crust and mantle: Tectonophysics, v. 75, p. 113–133, doi: 10.1016/ 0040-1951(81)90212-2. Behn, M.D., and Kelemen, P.B., 2006, Stability of arc lower crust: Insights from the Talkeetna Arc section, south-central Alaska, and the seismic structure of modern arcs: Journal of Geophysical Research, v. 111, doi: 10.1029/2006JB004327. Bloomer, S.H., Stern, R.J., Fisk, E., and Geschwind, C.H., 1989, Shoshonitic volcanism in the northern Mariana arc 1. Mineralogic and major and trace element characteristics: Journal of Geophysical Research, v. 94, p. 4469–4496. Carlson, R.L., and Gangi, A.F., 1985, Effect of cracks on the pressure dependence of P wave velocities in crystalline rocks: Journal of Geophysical Research, v. 90, p. 8675–8684. Christensen, N.I., and Mooney, W.D., 1995, Seismic velocity structure and composition of the continental crust: A global view: Journal of Geophysical Research, v. 100, p. 9761–9788, doi: 10.1029/95JB00259. Fujie, G., Ito, A., Kodaira, S., Takahashi, N., and Kaneda, Y., 2006, Confirming sharp bending of the Pacific plate in the northern Japan trench subduction zone by applying a traveltime mapping method: Physics of the Earth and Planetary Interiors, v. 157, p. 72–85, doi: 10.1016/j.pepi.2006.03.013. Jull, M., and Kelemen, P.B., 2001, On the conditions for lower crustal convective instability: Journal of Geophysical Research, v. 106, p. 6423–6446, doi: 10.1029/2000JB900357. Kay, R.W., and Kay, S.M., 1993, Delamination and delamination magmatism: Tectonophysics, v. 219, p. 177–189, doi: 10.1016/0040-1951(93)90295-U. Kelemen, P.B., 1995, Genesis of high Mg andesites and the continental crust: Contributions to Mineralogy and Petrology, v. 120, p. 1–19. 1034 Kelemen, P.B., and Holbrook, W.S., 1995, Origin of thick high-velocity igneous crust along the U.S. East Coast margin: Journal of Geophysical Research, v. 100, p. 10,077–10,094, doi: 10.1029/95JB00924. Kitamura, K., Ishihara, M., and Arima, M., 2003, Petrological model of the northern Izu-Bonin-Mariana arc crust: Constraints from high-pressure measurements of elastic wave velocities of the Tanzawa plutonic rocks, central Japan: Tectonophysics, v. 371, p. 213–221, doi: 10.1016/S0040-1951(03)00229-4. Kodaira, S., Sato, T., Takahashi, N., Ito, A., Tamura, Y., Tatsumi, Y., and Kanda, Y., 2007, Seismological evidence for variable growth of crust along the Izu intraoceanic arc: Journal of Geophysical Research, v. 112, doi: 10.1029/2006JB004593. Kushiro, I., 1974, Melting of hydrous upper mantle and possible generation of andesitic magma: Earth and Planetary Science Letters, v. 22, p. 294–299, doi: 10.1016/0012-821X(74)90138-1. Miller, D.J., and Christensen, N.I., 1994, Seismic signature and geochemistry of an island arc: A multidisciplinary study of the Kohistan accreted terrane, northern Pakistan: Journal of Geophysical Research, v. 99, p. 11,623– 11,642, doi: 10.1029/94JB00059. Oliver, S.B., Jones, C.H., and Sheehan, A.F., 2003, Foundering lithosphere imaged beneath the southern Sierra Nevada, California, USA: Science, v. 305, p. 660–662. Pearcy, L.G., DeBari, S.M., and Sleep, N.H., 1990, Mass balance calculations for two sections of island arc crust and implications for the formation of continents: Earth and Planetary Science Letters, v. 96, p. 427–442, doi: 10.1016/0012-821X(90)90018-S. Rapp, R.P., Shimizu, N., Norman, M.D., and Applegate, G.S., 1999, Reaction between slab-derived melts and peridotite in the mantle wedge: Experimental constraints at 3.8 GPa: Chemical Geology, v. 160, p. 3335–3356, doi: 10.1016/S0009-2541(99)00106-0. Rudnick, R.L., 1995, Making continental crust: Nature, v. 378, p. 571–578, doi: 10.1038/378571a0. Rudnick, R., and Fountain, D.M., 1995, Nature and composition of the continental crust: A lower crustal perspective: Reviews of Geophysics, v. 33, p. 267–309, doi: 10.1029/95RG01302. Smithson, S.B., Johnson, R.A., and Wong, Y.K., 1981, Mean crustal velocity: A critical parameter for interpreting crustal structure and crustal growth: Earth and Planetary Science Letters, v. 53, p. 323–332, doi: 10.1016/0012821X(81)90037-6. Stern, R.J., Fouch, M.J., and Klemperer, S.L., 2003, An overview of the Izu-BoninMariana subduction factory, in Eiler, J., ed., Inside the subduction factory: American Geophysical Union Geophysical Monograph 183, p. 175–222. Suyehiro, K., Takahashi, N., Ariie, Y., Yokoi, Y., Hino, R., Shinohara, M., Kanazawa, T., Hirata, N., Tokuyam, H., and Taira, A., 1996, Continental crust, crustal underplating and low-Q upper mantle beneath an oceanic island arc: Science, v. 272, p. 390–392, doi: 10.1126/science.272.5260.390. Takahashi, N., Kodaira, S., Klemperer, S.L., Tatsumi, Y., Kaneda, Y., and Suyehiro, K., 2007, Crustal structure and evolution of the Mariana intraoceanic island arc: Geology, v. 35, p. 203–206, doi: 10.1130/G23212A.1. Tamura, Y., and Tatsumi, Y., 2002, Remelting of an andesitic crust as a possible origin for rhyolitic magma in oceanic arcs: An example from the Izu-Bonin arc: Journal of Petrology, v. 43, p. 1029–1047, doi: 10.1093/ petrology/43.6.1029. Tatsumi, Y., 2005, The subduction factory; how it operates in the evolving Earth: GSA Today, v. 15, p. 4–10, doi: 10.1130/1052-5173(2005)015 [4:TSFHIO]2.0.CO;2. Tatsumi, Y., and Kogiso, T., 2003, The subduction factory: Its role in the evolution of the Earth’s crust and mantle, in Larter, R.D., and Leat, P.T., eds., Intraoceanic subduction systems: Tectonic and magmatic processes: Geological Society [London] Special Publication 219, p. 55–80. Taylor, S.R., 1967, The origin and growth of continents: Tectonophysics, v. 4, p. 17–34, doi: 10.1016/0040-1951(67)90056-X. Yuasa, M., 1991, The reason for the differences in geologic phenomena between northern and southern parts of the Izu-Ogasawara arc: Journal of Geography, v. 100, p. 458–463. Yuasa, M., and Nohara, M., 1992, Petrographic and geochemical along-arc variations of volcanic rocks on the volcanic front of the Izu-Ogawawara (Bonin) arc: Geological Survey of Japan Bulletin, v. 43, p. 421–456. Zelt, C.A., and Barton, P.J., 1998, 3D seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin: Journal of Geophysical Research, v. 103, p. 7187–7210, doi: 10.1029/97JB03536. Zhang, U., and Toksöz, M.N., 1998, Nonlinear refraction travel time tomography: Geophysics, v. 63, p. 1726–1737, doi: 10.1190/1.1444468. Manuscript received 19 March 2007 Revised manuscript received 11 May 2007 Manuscript accepted 29 June 2007 Printed in USA GEOLOGY, November 2007