Can Slab Melting Be Caused by Flat Subduction ?

Home , Adakite

Can slab melting be caused by flat subduction? GTS, 5573, F-34095 Marc-André Gutscher* Laboratoire UMR Université Montpellier II) Montpellier, France F-29280 René Maury , IUEM/Université Bretagne Occidentale, Place Nicolas CoDernic, Plouzané. France

e ae oretagne, Dar,IV, r-zyzuuriouzane, prance

ABSTRACT tion of adakitic magmas for these cases, related to Slab melting has been suggested as a likely source of adakitic arc magmas (i.e., andesitic an unusual subduction geometry known as flat and dacitic magmas strongly depleted in Y and heavy rare earth elements). Existing numerical subduction (Fig. 3) (Sacks, 1983; Pennington, and petrologic models, however, restrict partial melting to very young (15 Ma) oceanic crust 1984; Cahill and Isacks, 1992; Abbott et al., (typically at 60-SO km depth). Paradoxically, most of the known Pliocene-Quaternary adakite 1994; Gutscher et al., 1999a, 2000). Of the 10 occurrences are related to subduction of 10-45 Ma lithosphere,which should not be able to melt known flat slab regions worldwide (Fig. l), 8 are under normal subduction-zone thermal gradients. We propose an unusual mode of subduction linked to present or recent (<6 Ma) occurrences known as flat subduction, occurring in -10% of the world’s convergent margins, that can pro- of adakitic magmas (Table 1). duce the temperature and pressure conditions necessary for fusion of moderately old oceanic crust. Of the 10 known flat subduction regions worldwide, eight are linked to present or recent GEODYNAMIC MODEL (e6 Ma) occurrences of adakitic magmas. Observations from Chile, Ecuador, and Costa Rica Comparative observations from several well- suggest a three-stage evolution: (1) steep subduction produces a narrow calc-alkaline arc, typi- constrained flat slab regions along the Andean cally -300 km from the trench, above the asthenospheric wedge; (2) once flat subduction begins, and Central American margin suggest a three- the lower plate travels several hundred kilometers at nearly the same depth, thus remaining in stage evolution illustrating how a reduction of a pressure-temperature window allowing slab melting over this broad distance; and (3) once flat slab dip can lead to episodes of slab melting. In subduction continues for several million years, the asthenospheric wedge disappears, and a vol- the normal situation, a narrow calc-alkaline arc canic gap results, as in modern-day central Chile or Peru. The proposed hypothesis, which (volcanic line) develops above a steeply dipping reconciles thermal models with geochemical observations, has broad implications for the study slab, commonly between the 100 and 150 !un iso- of arc magmatism and for the thermal evolution of convergent margins, depth contours to the subducted plate (Fig. 3A). When the dip of the downgoing plate flattens, Keywords: flat subduction, thermal structure, slab melting, adakites. in response to a change in buoyancy (i.e., subduction of overthickened oceanic crust), the INTRODUCTION doxically, of the -20 known occurrences of arc typically widens, extending much farther Petrologic models suggest that formation of adakites worldwide, only 5 occur near Pliocene- (2400 km) from the trench (Fig. 3B and 3C). trondhjemite-tonalite-dacite(TTD) by partial Quatemary spreading center-trench triple junc- During this transitional stage, lasting several melting of the subducted slab was widespread tions, where very young oceanic lithosphere is million years, adakitic melts are generated during Archean time, and these rocks represent a being or was recently subducted. The remaining across a wide volcanic arc. This corresponds to major component of Precambrian gneiss terranes occurrences involve subduction of moderately modem Ecuador (arc 250-400 km from trench) (Martin, 1986; Drummond and Defant, 1990). old (10-45 Ma) lithosphere, which is not ex- (Monzier et al., 1997; Gutscher et al., 1999b; However, given the cooler mantle temperatures pected to melt under normal subduction-zone Bourdon et al., 1999) and to central Chile at during the Phanerozoic, these processes should pressure and temperature conditions. We present 10-5 Ma (arc 250-800 km from trench) (Kay be uncommon today (Martin, 1986). Recent a new geodynamic model to explain the forma- and Abbruzzi, 1996). Due to the unique set of models of magma genesis emphasize the role of fluids released via dehydration reactions from the subducting plate, thereby causing melting in the Figure 1. Global distribution of flat slab regions GOON overlying mantle wedge (Schmidt and Poli, (labels correspond to 1998). In the past decade, numerous occurrences Table 1) and adakitic mag- of Cenozoic adakites (andesitic and dacitic magmas (stars). Filled stars mas characterized by strong depletion in heavy are modern occurrences, unfilled stars are 1-6 Ma rare earth elements and high SrrY ratios) have old occurrences. Principal 30”N been documented and discussed as possible oceanic plateaus, hotspot examples of slab melting (Fig. 1) (Defant and tracks, and subducting Drummond, 1990; Moms, 1995). These include arcs are shaded graygray: no Hk-HikurancriHk-Hikurangi Plateau, 00 Mount St. Helens (Defant and Drummond, 1993) Lo-Louisville Ridge, and the recent major eruption of Mount Pinatubo Au-Austral Plateau, TÜ-Tu- (Prouteau et al., 1999). Tuamotu Plateau, Ma-Mq- Numerical and petrological studies of pressure- Marquesas Plateau, Mh- 30”s temperature-time (P-T-t) paths in subduction Manihiki Plateau, Li-Line Islands, OJ-Ontong Java zones suggest that this process can only occur for Plateau, ER-Euripik Rise, subduction of very young (15 Ma) lithosphere PK-Palau-Kyushu Ridge, (Peacock et al., 1994) (Fig. 2A, solid lines). Para- IB-lzu-Bonin Arc, Sh- 555 Shatsky Rise, MP-Mid- Pacific Mountains, He-Hess Rise, Ha-Hawaiian Chain, Em-Emperor Chain, YB-Yakutat “Current address: GEOMAR,Marine Geodynamics, Block, AS-Alaska Seamounts, Co-Cocos Ridge, GC-Galápagos-Carnegie Ridge, Nz-Nazca Wischhofstrassz- 1-3, D-24148, GeLmany. E-mail: Ridge, JF-Juan-Fernandez Ridge. mgssche&geomar.de. I Geo&g&JuLtgOO; v. 28; no. 6; p. 535-538; 3 hgures; 1 table. A varying age of subducting lithosphere (in my.)

4.0

3.0 slab aths P (GPa) 2.0

1.o

O O 1O0 200 300 [km] 400 500 600 O 200 400W) 600 800 1000 T B broad 'adakilic' arc R Central Chile (10-7 Ma) Ecuador (O Ma) O P-T-t paths based on isotherms in Figure 3 .._.____..._...... ______..__.300 __.__ I*I(I -~-LLLLL-2---DL-~~o-.~.=. v = 6 cm/a ...... ____ 7. ______y....:__.._.: _...... -...... ~ ...... ~ ...... 900 ....~

asthenosphere asthenosphere

150 -.-,,$QQ-""-- ..__.____~ --.-.._ - O 100 200 300 [km] 400 500 600

O 200 400W) 600 800 1000 T Figure 2. Pressure-temperature (f-T) meta- mbrphic reaction diagram showing expected f-%t (t = time) paths for varying conditions. Ec-eclogite, Am-amphibolite, Ga-garnet, Hb-hornblende. A P-Ft paths for normal 27" subduction dip (solid lines), with varying lithospheric ages, and subduction velocity of 3 cm/yr, with no shear heating, according to finite difference models (Peacock et al., 1994, O 100 200 300 tkml 400 500 600 Fig.8B). f-T-tpathsforflat subduction (dashed) transcunent faulting and/or block fault uplik (Sierras Pampeanas) are calculated using analytical solution for D Central Chile 10 Mal no aCliYB YO,caniSm -- - slab thermal structure (Davies, 1999) beneath 70 km depth.These paths are much flatter and intersect slab melting field (shaded) at 700 "C and 2.5 GPa for 10-50 Ma subducted lithosphere. B: P-T-t paths for steep (circles) and flat (diamonds and squares) subduction based on isotherms in Figure 3sampled every 50 km along slab surface. isotherms and P-T-t paths involved in flat sub- O 1O0 zoo 300 [kml 400 500 600 duction, slab melting can occur (Fig. 2). Pro- Figure 3.Three-stage evolution showing transition from steep to flat subduction style due to sub- longed flat subduction will sufficiently cool both duction of buoyant, overthickened oceanic crust (e.g., oceanic plateau), with impact on mag- lithospheres such that the downgoing oceanic matic arc.Therma1 structure is constructed based on adaptation of published numerical models crust dehydrates before the necessary tempera- (Peacocket al., 1994; Peacock, 1996) and analytical solution (Davies, 1999) (see text). A: Steep or tures for partial melting are attained and thus the normal subduction (230" dip), with narrow calc-alkaline arc, -300 km from trench produced by slab dehydration and partial melting in overlying asthenospheric wedge. B: Early stage mag- final stage, a volcanic gap, is achieved (Fig. 3D). matic flat subduction; oceanic slab crosses 700 "C geotherm at -80 km depth (300-400 km from This complete evolution has been documented trench) allowing broad adakitic arc to develop. Partial melting continues in narrow intervening for central Chile between 20 Ma and the present tongue of asthenospheric material. C: Late stage magmatic flat subduction; intervening astheno- (Kay and Abbruzzi, 1996). Ecuador is known to spheric tongue cools and retreats. Partial melting can only occur at great distance (-450- 600 km) from trench, and final pulse of magmatism occurs, as in Sierras Pampeanas, Pocho have had a narrow calc-alkaline arc at 5 Ma eruption 4.7 Ma (Kay and Abbruui, 1996). D: Prolonged (amagmatic) flat subduction; oceanic (Barben et al., 1988) and has evolved to stage 2 slab is cooler than 600 "C, asthenospheric wedge has disappeared, and thus no partial melting today due to flattening of the slab in response to can occur, and volcanism ceases. subduction of Camegie Ridge (Pennington, 1981; Gutscher et al., 1999b). The Costa %can margin also exhibits a simi- manca Cordillera, suggesting that the arrival of from slab melting (Defant et al., 1992). Today, lar magmatic evolution (Defant et al., 1992). Cocos Ridge began modifying the style of sub- the Talamanca segment of Costa Rica appears Before 6 Ma, magmatism was calc-alkaline in duction to the current flat slab configuration to have advanced to the final stage, a volcanic nature and the subduction angle presumably (Protti et al., 1994). Since 5 Ma, adakitic mag- gap, because no Quaternary stratovolcanoes are steep. Around 5 Ma uplift began in the Tala- matism developed, interpreted as originating present along this 100 km segment.

536 GEOLOGY, June 2000 TABLE I.FLAT SLAB REGIONS WORLDWIDE (Atherton and Petford, 1993; Petford and Ather- # Region' Length' Associated plateau@) Ages Volcanism ton, 1996; Kay and Abbruzzi, 1996), and melting (km) or blocks (Ma) Quaternary Age of adakites of the base of an accreted arc or ophiolite complex, arc (Ma) as invoked for Ecuador (Arculus et al., 1999). 1 Chile (28°S-330S) 550 Juan Femandez Ridge 43 No 4-7 Peru Nazca Ridge, Inca Plateau No The slab window model in general applies to 2 (2°S-150S)(lOS-ZON) 1500 3043 3-6 3 Ecuador 350 Carnegie Ridge 16-24 Yes 0-3 trench-spreading center triple junctions, and 4 Columbia (6"N-g0N) 350 Choco block 20 No (none known) thus could succeed in explaining the Chile triple 5 Costa Rica 250 Cocos Ridge 14-20 No' 1-3 6 Mexico 400 Tehuantepec Ridge 12-20 Yes (none known) junction and possibly the Panama occurrence 7 Cascadia (46ON-49'") 350 (none known) a Yes O 8 Southern Alaska 500 Yakutat terrane 45 No o" (assuming that a spreading center segment was 9 Southwestern Japan 600 Izu Bonin Arc, Pal.-Ky. Ridge 15-20 Yes O recently subducted). However, it is insufficient 10 New Guinea (135°E-1400E) 550 Euripik Ridge -20 No 2-41' ' 'Sources: (1)Kay and Abbruui, 1996, (2)Petford and Atherton, 1996. (3) Monzier et al., 1997, (4)Penninaton. 1981. to explain the numerous cases where moderately (5)Defant et al., 1992,(6) Suarez et al., 1990, (7) Defant and Drummond, 1993, (8) Preece, 1991, (9)Morri< 1995,(IO) old oceanic crust is subducting (e.g., Peru, cen- Pubellier et al., 1998. See also compilations (Pennington, 1984; McGeary et al., 1985; Abbott et al., 1994; Gutscher et al., 2000). tral Chile, Japan, Alaska) and where no evidence 'Cumulative length = 5400 km, which represents -10% of all subduction zones. of lithospheric tears exists. "ge of subducted lithosphere at trench. Melting at the base of an overthickened oro- "In Costa Rica there is a 100 km gap (Talamanca segment) with an absence of stratovolcanoes. "In eastern Alaska there is a 500-km-wide volcanic gap. The adjacent volcanoes Hayes and Wrangel1 show adakitic genic crust (gamet eclogite) has been suggested signatures. for Chile (Kay and Abbruzzi, 1996),Peru (Petford while SrN ratios for 4.4-2.6 Ma intrusions plot ìn the adakitic field, new Isotope data suggest strong contamination by Australian continental crust (Housh and McMahon, 2000). and Atherton, 1996), and Ecuador (Monzier et al., 1997) based on isotopic imprints and trace element features indicative of crustal compo- nents. However, in Ecuador the crustal imprint is THERMAL EVOLUTION AND for the thermal evolution resulting from a change known to be weak (45%) (Barragan et al., PETROLOGICAL CONSTRAINTS in subduction style from steep to flat, a thermal 1998; Bourdon et al., 1999). Furthermore, it is Although prolonged flat subduction generally model was constructed using an adaptation of very difficult to distinguish whether these char- produces a cooler thermal structure (i.e., both published numerical and analytical solutions. The acteristics reflect melting of the base of the crust, upper and lower lithospheres are cooler) (Fig. 3D), steep slab thermal strqcture is based on both finite or slab melting followed by crustal contamina- (Sacks, 1983; Vlaar, 1983; Henry and Pollack, difference models (Peacock et al., 1994; Peacock, tion during ascent through the Andean crust 1988; Dumitru et al., 1991), during the initiation 1996) and the analytical solution (Davies, 1999) (e.g., the MASH model, Hildreth and Moorbath, of flat subduction the leading edge of the under- and shows good agreement with earlier (Davies 1988). Because these two processes often cannot sliding slab can be anomalously heated. An and Stevenson, 1992; Stein and Stein, 1996) and be distinguished geochemically, the reasons analytic solution for the thermal structure of a recently published finite-element calculations commonly given for dismissing the slab-melting cold slab subducting through hot asthenosphere (Oleskevich et al., 1999; Peacock and Wang, model are thermal considerations, i.e., that crust demonstrates temperature to be a function of the 1999) for slabs with comparable subduction older than 20 Ma cannot melt (Petford and downdip length, the slab thickness, and the rate parameters (age, velocity, dip angle). The flat slab Atherton, 1996; Kay and Abbruzzi, 1996). of subduction, but independent of the subduction thermal structure below 70 km is calculated using The first argument against melting of the base dip, thus allowing application to a slab with vari- the analytical solution as described here (Davies, of a crust thickened by underplating of mafic able dip down its length (Davies, 1999). In other 1999). The subduction parameters and geometry material is that, while this explanation may work words, a slab moving a fixed distance through hot applied here are based loosely on typical Andean for parts bf South America, it does not explain asthenosphere will gain heat at the same rate values (25 Ma lithosphere with a 60 km thermal Cascadia, southern Alaska, or southwestern whether it descends vertically, at 60°, at 30°, or thickness, velocity = 6 cm/yr, and a 300 km sub- Japan, where continental crustal thicknesses do even if it moves horizontally (as in Fig. 3B and horizontal slab at 80 km depth). It assumes no not exceed 40 km (Oleskevich et al., 1999). 3C). This allows published steep-slab isotherms shear heating, no heat of fusion, and no advective Moreover, the relatively high Mg, Cr, and Ni (Peacock et al., 1994; Peacock, 1996) to be re- heat transport to the arc. This set of isotherms is contents of Andean adakites and related rocks calculated for a flat slab geometry for pressures sampled every 50 km along the downdip direction are best explained by interaction with mantle >2 GPa (-70 km depth) (Fig. 2A). Due to the un- in order to illustrate the P-T-r paths predicted by peridotite (Sen and Dunn, 1995). usual path traveled, first downward to 70 km this evolutionary model, shown in Figure 2B, by The strongest argument against crustal melt- depth, and then subhorizontally for several hun- the use of symbols (circles, diamonds, and ing (either of continental or accreted oceanic dred kilometers, the P-T-t paths differ greatly squares), each corresponding to a position in P-T materials) is temporal. All well-studied adakitic from steep subduction paths and cross over the space (Fig. 2B). Because the leading edge of the provinces are known from radiometric dating to wet solidus into the slab melting field. undersliding slab spends -5 m.y. at a nearly con- have been active for only a brief time span of To date, all published numerical models of sub- stant depth of 70-80 km (pressures of 2-2.5 GPa), 2-3 m.y. or less (e.g., Panama and Costa Rica, duction-zone thermal structure feature a constant it has sufficient time to heat to 700-800 OC, Defant et al., 1992). Adakitic magmatism in cen- (Davies and Stevenson, 1992; Furukawa, 1993; leading to the fluid-present partial melting of the tral Chile ceased at 4.7 Ma (Kay and Abbruzzi, Peacock et al., 1994; Peacock, 1996; Stein and oceanic crust (Drummond and Defant, 1990; 1996), yet subduction beneath the thickened Stein, 1996) or constantly increasing (Oleskevich Peacock et al., 1994; Prouteau et al., 1999). crust continues today. The same is true for et al., 1999; Peacock and Wang, 1999) slab dip, Ecuador, where subduction has been occurring and are typically run until thermal equilibrium is DISCUSSION and arc magmas have been ascending through reached. No numerical study has yet attempted to Alternate models have been proposed for pro- the accreted ophiolitic crust for at least 20 m.y. rigorously model a reduction in slab dip (flat sub- ducing adakitic-like magmas in a convergent Why is there no older (10-20 Ma) adakitic mag- duction), including its transient thermal structure. margin setting; melting of the edges of oceanic matism and why does its onset correspond to the Past attempts, either static (Vlaar, 1983) or one di- crust in the context of a slab window (Johnston change in subduction dip? mensional [Spencer, 1994), have proven unsatis- and Thorkelson, 1997). direct melting of the base It is important to note that our flat subduction factory. In order to illustrate a plausible sequence of an overthickened orogenic continental crust slab melting model includes this temporal evolu-

June 2000 GEOLOGY, 537

...... _ ...... C._ ...... - I ~.... tion. It is only during the early stages of flat sub- Drummond, M.S., andDefant,M.J., 1990,Amodel for Peacock, S.M., Rushmer, T., and Thompson, A.B., duction that the upper part and leading edge of the trondhjemite-tonalite-dacitegenesis and crustal 1994, Partial melting of subducting oceanic crust: growth via slab melting: Archean to modern com- Earth and Planetary Science Letters, v. 121, slab can melt. Prolonged flat subduction will cool parisons: Joumal of Geophysical Research, v. 95, p. 227-244. both lithospheres and impede partial melting. It p. 21,503-21,521. Pennington, W.D., 1981, Subduction of the Eastern has been repeatedly suggested that only the sub- Dumitru, T.A., Gans, P.B., Foster, D.A., and Miller, Panama basin and seismotectonics of northwest- duction of hot, very young lithosphere will allow E.L., 1991, Refrigeration of the western ern South America: Joumal of Geophysical Re- slab melting to occur (Defant and Drummond, Cordilleran lithosphereduring Laramide shallow- search, v. 86, p. 10,753-10,770. angle subduction: Geology, v. 19, p. 1145-1 148. Pennington, W.D., 1984, The effect of oceanic crustal 1993; Peacocket al., 1994). The emphasis here is Furukawa,Y., 1993, Depth of the decoupling plate inter- structure on phase changes and subduction: not on the initial temperature of the slab, but on face and thermal structure under arcs: Journal of Tectonophysics, v. 102, p. 377-398. its buoyancy and thus on the time spent at a rela- Geophysical Research, v. 98, p. 20.005-20,013. Petford, N., and Atherton, M., 1996, Na-rich partial tively shallow (-80 km) depth, allowing the slab Gutscher, M.-A., Olivet, J.-L., Aslanian, D., Maury, R., melts from newly underplated basaltic crust: The and Eissen, J.-P., 1999a, The "lost Inca Plateau": Cordillera Blanca batholith, Peru: Journal of to heat. Therefore, oceanic crust as old as 50 Ma Cause of flat subduction beneath Peru?: Earth and Petrology, v. 37, p. 1491-1521. can melt if it is given sufficient time to warm. Planetary Science Letters, v. 171, p. 335-341, Preece, S.J., 1991, Tephrostratigraphyof the lateCeno- Young crust (5-10 Ma) subducts steeply at lat Gutscher, M.-A., Malavieille, J., Lallemand, S., and zoic Gold Hill Loess, Fairbanks area, Alaska 42"s in Chile, yet no adakites are observed be- Collot, J.-Y, 1999b,Tectonic segmentation of the [Master's thesis]: Toronto, University of Toronto, tween lat 42" and North Andean margin: Impact of the Carnegie 164 p. 44"s. Ridge collision: Earth and Planetary Science Let- Protti, M., Guendel, F., and McNally, K., 1994, The ACKNOWLEDGMENTS ters, v. 168, p. 255-270. geometry of the Wadati-Benioff zone under We thank the reviewers Richard Arculus and Marc Gutscher, M.-A., Spakman, W., Bijwaard, H., and southern Central America and its tectonic signif- Hirschmann for constructive suggestions that helped Engdahl, E.R., 2000, Geodynamics of flat sub- icance: Results from a high-resolution local improve the original manuscript. M.-A. Gutscher's reduction: Seismicity and tomographic constraints seismographic network: Physics of the Earth and search was supported by a Marie Curie Training and from the Andean margin: Tectonics (in press). Planetary Interiors, v. 84, p. 271-287. Mobility of Researchers postdoctoral stipend from the Henry, S.G., and Pollack, H.N., 1988, Terrestrial heat Prouteau, G., Scaillet, B., Pichavant, M., and Maury, European Union. flow above the Andean subduction zone in Bolivia R.C., 1999, Fluid-present melting of oceanic and Peru: Journal of Geophysical_. Research, v. 93, crust in subduction zones: Geology, v. 27, REFERENCES CITED p. 15,153-15,162. p. 1111-1114. Abbott, D., Drury, R., and Smith, W.H.F., 1994, Flat to Hildreth.,. W.. andMoorbath. S.. 1988. Cmstal contribu- Pubellier, M., Prouteau, G., Maury, R.C., and Bellon, steep transition in subduction style: Geology, tion to arc magmatism in the Andes of central H., 1998, Le changement de régime plio-quater- Chile: Contributions to Mineralogy and Petrol- naire dans une zone d'&happement tectonique; V. 22, p. 937-940. Arculus, R.J., Lapierre, H., and Jaillard, E., 1999, ogy, v. 98, p. 455-489. contrôle par des intrusifs syntectoniques (Irian Geochemical window into subduction and accre- Housh, T., and McMahon, T.P., 2000, Ancient isotopic Jaya, Indonésie) [abs.]: Proceedings Reunion tion processes: Raspas metamorphic complex, characteristics of Neogene potassic magmatism des Sciences de la Terre Conference, Brest, Ecuador: Geology, v. 27, p. 547-551. in Western New Guine (Irian Jaya, Indonesia): France, p. 179-180. Atherton, M., and Petford, N., 1993, Generation of Lithos, v. 50, p. 217-239. Sacks, IS., 1983, The subduction of young litho- sodium-rich magmas from newly underplated Johnston, S.T., andThorkelson,D.J., 1997, Cocos-Nazca sphere: Journal of Geophysical Research, v. 88, basaltic crust: Nature, v. 362, p. 144-146. slab window beneath Central America: Earth and p. 3355-3366. Barberi, F., Coltelli, M., Ferrara, G., Innocenti, F., Planetary Science Letters, v. 146, p. 465-474. Schmidt, M.W., and Poli, S., 1998, Experimentally Navarro, J.M., and Santacroce, R., 1988, Plio- Kay, S.M., and Abbruzzi, J.M., 1996, Magmatic evi- based water budgets for dehydrating slabs and Quaternary volcanism in Ecuador: Geological dence for Neogene lithospheric evolution of the consequences for arc magma generation: Earth Magazine, v. 125, p. 1-14. central Andean "flat-slab" between 30"s and and Planetary Science Letters, v. 163, p. 361-379. Barragan, R., Geist, D.,Hall, M., Larson, P., and Kurz, 32"s: Tectonophysics, v. 259, p. 15-28. Sen, C., and Dum, T., 1995, Experimental modal M., 1998, Subduction controls on the composition Martin, H., 1986, Effect of steeper Archean thermal metasomatism of a spinel lherzolite and the pro- of lavas from the Ecuadorian Andes: Earth and gradient on geochemistry of subduction-zone duction of amphibole-bearing peridotite: Contri- Planetary Science Letters, v. 154, p. 153-166. magmas: Geology, v. 14, p. 753-756. butions to Mineralogy and Petrology, v. 119, McGeary, Nur, A., and Ben-Avraham, 1985, Bourdon, E., Eissen, J.-P., Cotten, J., Monzier, M., S., Z., p. 394-409. Robin, C., and Hall, M.L., 1999, Les laves calco- Spatid gaps in arc volcanism: The effect of colli- Spencer, J.E., 1994, A numerical assessment of slab sion or subduction of oceanic plateaus: Tectono- strength during high- and low-angle subduction alcalines et à caractère adakitique du volcan Anti- sana (Equateur): Hypothèse pétrogénétique: Paris, physics, v. 119, p. 195-221. and implications for Laramide orogenesis: Journal Academie des Sciences Comptes Rendus, v. 328, Monzier, M., Robin, C., HaU, M.L., Cotten, J., Mothes, P., of Geophysical Research, v. 99, p. 9227-9236. p. 443-449. Eissen, J.-P., and Samanigo, l?, 1997,Les adakites Stein, S., and Stein, C.A., 1996, Thermo-mechanical Cahill,T., andIsacks, B.L., 1992, Seismicity and shape d'Equateur: Modèle preliminaire: Paris, Academie evolution of oceanic lithosphere: Implicationsfor of the subducted Nazca plate: Journal of Geo- des SciencesComptes Rendus, v. 324, p. 545-552. the subduction process and deep earthquakes, in physical Research, v. 97, p. 17,503-17,529. Moms, P.A., 1995, Slab melting as an explanation of Bebout, G.E., et al., eds., Subduction: Top to bot- Davies, J.H., 1999, Simple analytic model for subduc- Quaternary volcanism and aseismicity in south- tom: American Geophysical Union Geophysical tion zone thermal structure: Geophysical Journal west Japan: Geology, v. 23, p. 395-398. Monograph 96, p. 1-17. International,v. 139, p. 823-828. Oleskevich, D.A., Hyndman, R.D., and Wang, K., Suarez, G., Montfret, T., Wittlinger, G., and David, C., Davies, J.H., and Stevenson, D.J., 1992, Physical 1999, The updip and downdip limits to great sub- 1990, Geometry of subduction and depth of the duction earthquakes: Thermal and structural seismogenic zone in the Guerrero gap, Mexico: model of source region of subduction zone vol- canics: Journal of Geophysical Research, v. 97, models of Cascadia, south Alaska, SW Japan, and Nature, v. 345, p. 336-338. p. 2037-2070. Chile: Journal of Geophysical Research, v. 104, Vlaar, N.J., 1983, Thermal anomalies and magmatism Defant, M.J., and Drummond, M.S., 1990, Derivation of p. 14,965-14,991. due to lithospheric doubling and shifting: Earth some modern arc magmas by melting of young Peacock, S.M., 1996, Thermal and petrologic struc- and Planetary Science Letters, v. 65, p. 322-330. subducted lithosphere: Nature, v. 347, p. 662-665. ture of subduction zones, in Bebout, G.E., et Defant, M.J., and Drummond, M.S., 1993, Mount St. al., eds., Subduction: Top to bottom: Ameri- Manuscript received October 5,1999 Helens: Potential example of the partial melting can Geophysical Union Geophysical Mono- Revised manuscript received March 14,2000 of the subducted lithosphere in a volcanic arc: graph 96, p. 119-133. Manuscript accepted March 28,2000 Geology, v. 21, p. 547-550. Peacock, S.M., and Wang, K., 1999, Seismic conse- Defant, M.J., Jackson, T.E., Drummond, M.S., de Boer, quences of warm versus cool subduction meta- J.Z., Bellon, H., Feigenson, M.D., Maury, R.C., morphism: Examples from southwest and north- and Stewart, R.H., 1992, The geochemistry of east Japan: Science, v. 286, p. 937-939. young volcanism throughout westem Panama and southeastern Costa Rica: An overview: Geologi- cal Society of London Joumal, v. 149, p. 569-579.

538 Printed in USA GEOLOGY, June 2000 .

j.. - - - - . , . ..

EDITORS David E. Fastovsky Ben A. van der Pluijm of University of Rhode Island University Michigan (401) 789-0238 (734) 615-3050 fax 401-783-61 55 fax 734-763-4690 E-mail: [email protected] E-mail:[email protected] Editor’s Assistant: Editor’s Assistant: Leslev A. Fastovskv Carol Traynor VOLUME 28 NO. 6 P. 481-576 June 2000

483 Life associated with a 2.76 Ga ephemeral 531 Channel migration and meander-bend pond?: Evidence from Mount Roe #2 paleosol curvature in the lower Mississippi River prior D, Rob Rye, Heinrich Holland to major human modification 487 Dating the Indian continental subduction Paul E Hudson, Richard H. Kesel and collisional thickening in the northwest 535 Can slab melting be caused by flat subduction? Himalaya: Multichronol.ogy of the Tso Morari Marc-AndréGutscher, RetJe‘ Maury, jean-PhilippeEissen, ecloaites ’ e o Erwan Bourdon lulia Ge Sigoyer, Valérie Chavagnac, Iai~neBlichert-Toft, 539 Weddell triple junction: The principal focus of /gor M. Villa, Béatrice Luais, Stéphane Guillot, Ferrar and Karoo magmatism during initial Michael Cosca, Georges Masde breakup ofTH. Gondwana 491 Reduced porphyry copper-gold deposits: A new D.H. Elliot, Fleming variation on an old theme 543 New high-pressure granulite event in the Stephen M. ROW~IJS Moine Supergroup, northern Scotland: 495 Earthquake focal depths, effective elastic Implications for Taconic (early Caledonian) thickness, and the strength of the continental crustal evolution K. lithosphere C, R.I. Friend, A. ]ones, I. M. BurIJs A. Maggi, I.A. lackson, D. McKenzie, K, Priestley 547 Overlapping volcanoes: The origin of Hilo Ridge, Hawaii 499 Dimethylsulfide production variations over T: the past 280 k.y. in the equatorial Atlantic: Robin Holcomb, Bruce K. Nelson, Peter W. Reiners, A first estimate Nuni-lynSawyer Anders S. Henriksson, Michael Sarr~thein, 551 Central American paleogeography controlled ’f Geoffrey Eglii~ton,1017 Poynter - i Pliocene Arctic Ocean moiiuscan migrations 1 503 Subaqueous lava flow lobes, observed on Louie Marincovich Ir. ROW KAIKO dives off Hawaii 555 Relationship between copper tonnage of Susuniu Umino, Peter W. Lipman, Sumie Obata Chilean base-metal porphyry deposits and 507 Lithium isotope evidence for light element Os isotope ratios decoupling in the Panama subarc mantle Ryan Mathur, loaquin Ruiz, Francisco Munizaga Paul B. Tomascak, )effrey G. Ryan, Marc). Defant 559 Alternative origin of aliphatic polymer In Coordinated strike-slip and normal faulting kerogen 51 1 D. in the southern Ozark dome of northern B.A. Stankiewicz, E.G. Briggs, R. Michels, Arkansas: Deformation in a late Paleozoic M.E. Collinson, M.B. Flannery, R.P. Evershed foreland 563 Re-Os isotope characteristics of postorogenic Mark R. Hudson lavas: Implications for the nature of young 51 5 Sclerosponges as a new potential recorder of lithospheric mantle and its contribution to enviro n ment a I ch a n g es: Lead in Ceratopor ella basaltic magmas P. nicholsoni Bruce E Schaefer, Simon Turner, Nick W, Rogers, Claire E. Lazareth, Philippe Willenz, lacques Navez, Chris ). Hawkesworfh, Helen M. Williams, Eddy Keppens, Frank Dehairs, Luc André D. Graham Pearson, Geoffrey M. Nowell 51 9 lntralithospheric differentiation and crustal 567 Names are better than numbers: growth: Evidence from the Borborema Nomenclature for sequences province, northeastern Brazil P. P. lake M. Hancock Sérgio Neves, Gorki Mariano, lgnez Guimarães, 569 Palynofacies prediction of distance from Adejardo E da Silva Filho, Silvana C. Melo sediment source: A case study from the 523 Dinosaur abundance was not declining in a Upper Cretaceous of the Pyrenees m “3 gap” at the top of the Hell Creek Richard V. Tyson, Ben Follows Formation, Montana and North Dakota Peter M. SheehaiJ, David E. Fastovsky, Claudia Barreto, Forum G. Raymond Hoffmann 572 Depth to pedogenic carbonate horizon 527 Crustal suture preserved in the Paleo- as a paleoprecipitation indicator? moterozoic Trans-Hudson oroaen, Canada Comment: Gregory). Retallack

(continued on next page) CONTENTS (continued) COVER 573 The Channeled Scabland: Back to Bretz? Optical photomicrograph, under cross-polarized light, of high- Comment: G. Komatsu, H. Miyamoto, K. /to, pressure metamorphosed pelite from the Tso Morari dome T: (eastern Ladakh, northwest Himalaya). This rock shows eclogitic H. Tosaka, Tokunaga foliation composed of stable Si phengite, kyanite, and quartz. Reply: lohn Shaw, Mandy Munro-Stasiuk, Brian Sawyer, Tso Because the Morari unit belongs to the Indian margin, the Claire Beaney, jerome-Etienne Lesemann, high-pressure, low-temperature assemblage implies the Alberto Musacchio, Bruce Rains, Robert R. Young E subduction of the Indian margin down to a depth of 70 km. Comment: Brian Atwater, Cary A. Smith, During the decompression, the eclogitic foliation is overprinted Richard 6. Waitt by recumbent folds and deformation leads to mineral re- Reply: lohn x25. Shaw, Mandy Munro-Stasiuk, crystallization in the hinge of the folds, Magnification See Brian Sawyer, Claire Beaney, lerome-EtienneR. Lesemann, “Dating the Indian continental subduction and collisional Alberto Musacchio, Bruce Rains, Robert Young thickening in the northwest Himalaya: Multichronology of the Tso Morari eclogites,” by Julia de Sigoyer et al., p. 487. Photo by Julia de Sigoyer. Cover design by Caynor Ann Bloom.

GSA‘s The full text of Geology(www.geosociety.org.pubs/journals.htm). articles is now available online at Manuscript submissions for Geology should be sent to Geology, Web page for Geology CSA, P.O. Box 9140, Boulder, CO 80301-91 40. Complete Articles are available as PDF files that may be viewed using information for contributors is posted on the GSA Web site at Adobe Acrobat Reader. www.geosociety.org.

GEOLOGY EDITORIAL BOARD Geological Society of America 1999-2000 Read, J. Fred Hacker, Bradley R. Bierman, Paul R. Andrews, lulian Ring, Uwe Harris, W. Burleigh Bird, Michael I. CHIEF EXECUTIVE OFFICER Boles, ]im Rodgers, David Hirschmann, Marc M. Birkeland, Peter W. Bottier, David Rollinson, Hugh Holliger, Klaus Cloetingh, Sierd Sara S. Foland Bowring, Sam Self, Steve lenkyns, Hugh Dean, Walter Brown, Mike ; Shinohara, Hiroshi Jones, Alan G. DeCeIles, Peter PRESIDENT VICE-PRESIDENT Canuel, Liz Sial, Alcides Kerrick, Derrill M. Dewey, John F. Mary LOUZoback Sharon Mosher Clift, Peter Tauxe, Lisa kharaka, Yousif Droser, Mary L. Dalziel, Ian Vallance, lames Koeberl, Christian Farley, Kenneth A. PAST PRESIDENT TREASURER DeMets, Chuck Wallace, Paul Koons, Peter C. Gardner, Thomas W. de Ronde, Cornel Wallis, Simon Maples, Christopher Hand, Martin Gail M. Ashley David E. Dunn D. Dipple, Greg Walter, Mike Marshall, lames Hovius, Niels Ferrari, Luca Watts, Anthony Meurer, William Kaufman, Darrell S. COUNCILORS Fisher, Andy Mills, Rachel A. Meyer, Grant A. 1999-2000 Fitzsimons, Ian 1999-2001 Mutti, Maria Mezger, Klaus K. Douglas W. Burbank Flessa, Karl W. Ahlberg, Anders Niemi, Tina Montgomery, David R. Furbish, Dave Alley, Richard B. Oppo, Delia Mountjoy, Eric Rhea Lydia Graham Allison, Peter D. Allison R. (Pete) Palmer Fyk, W. S. Ortiz, Joseph Normark, William R. Garrison, Robert Anderson, Robert S. Pomar, Luis Norris, R. D. Joaquin S. Ruiz Gerbe, Marie-Christine Baumgartner, Lukas P. Ross, Gerald M. Pain, Colin D. Hirth, J. Greg Blum, Joel Rowan, Elisabeth Lanier Pratson, Lincoln Frederick Bond, Gerard Roy, 1999-2001P. Hu, Feng Sheng I. Mousumi Ratschbacher, Lothar Noel James Keigwin, Lloyd Bralower, Timothy Rushmer, Tracy Rea, David K. Kelley, Patricia Brandon, Mark Rusmore, Margaret E. Ruiz, loaquin 5. Carol Simpson Buick, Ian Rob Van der Voo Mancktelow, Neil Santanach, Pere Ruppel, Carolyn D. G. Marzo, Mariano Bursik, Marcus Strecker, Manfred R. Sahagian, Dork L. Stephen Wells Mastin, Larry Chadwick, Oliver A. Stuewe, Kurt Shinn, Eugene A. Milliman, John Clemens, Steve Tait, Stephen Sloan, Lisa C. 2000-2002 Murray, Chris CoBabe, Emily White, Noel C. Spray, John G. Mary P. Anderson Nabelek, Peter Cuffey, Kurt Whitlock, Cathy Stern, Robert I. Rena M. Bonem Parmentier, Marc Defant, Marc J. Wilson, Paul A. Thomas, Ellen Ir. Parsons, Tom Delaney, Margaret .Willett, Sean Harry Y. (Hap) McSween Duncan, Robert 2000-2002 Barbara J. Tewksbury Phillips, Fred M. Zoback, Mark D. Powell, Chris Filippeili, Gabriel Arculus, Richard 1. Gray, David Baumiller, Tomasz Quade, lay R.

GEOLOGY (ISSN 0091-7613 USPS 994-580 CODEN GLGYB) is published monthly by the Geological Society of America, Inc., with offices at 3300 Penrose Place, Boulder, Colorado. Mailing address is P.O. Box 9140, Boulder, CO 80301-9140, U.S.A. Periodicals postage paid at Boulder, Colorado, and at additional mailing offices. Postmaster: Send address changes to Geology, Membership Services, P.O.Box 9140, Boulder, CO 80301-9140. Ccpyright O 2000, The Geological Society of America, Inc. (GSA). Ali rights reserved. Copyright not claimed on content prepared wholly by US. Government employees within the scope of their employment. Individual scientists are hereby granted permission, without fees or further requests to GSA, to: (1) use a single figure, a single table, and/or a brief paragraph of text in other subsequent works; (2) make unlimited photocopies of items in this journal for use in classrooms to further education and science; and (3) make up to five copies for distribution to associates in the furtherance of education and science. If $0.25 you need to make more than five photocopies for other noncommercial, http://www.copyrightcom;nonprofit purposes furthering science and education, remit a fee of per page-copy directly to the Copyright Clearance Center, 222 Rosewood Drive, Danvers, MA O1 923 USA, phone (978) 750-8400, when paying, reference Geology, ISSN 0091-761 3. Written permission is required from GSA for ail other forms of capture, reproduction, and/or distribution of any item in this journal by any means, including posting of part or all of the article on authors‘ or organizational Web sites, except that permission is granted to authors to post the abstracts only of their articles on their own or their organization’s Web site providing the posting includes this reference: ”The full paper was published ín the Geological Society of America’s journal Geology, [include year, month, and page number if known, where article appears or will appear].” 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. SUBSCRIPTIONS for 2000 calendar year: GSA Professional Members & Fellows 165 (2 years, $117); GSA Student Members and Associates, and K-12 Teachers $40 (2 years, $72); foreign members $77 (2 years, $139). Nonmembers and institutions: $355 for U.S./Canada/Mexico, add $12 overseas. Formats available: (1) monthly in paper-copy, as priced above; (2) full-text PDF file online at www. geosociety.org (no charge through 2000); (3) biannually on CD-ROMwith Bulletin, GSA Today, GSA Data Repository, and an Electronic Retrospective Index to journal articles from 1972: GSA Members $94 (2 years, 5169); GSA Student Members and Associates, and K-12 Teachers $64 (2years, $1 15); others, $355 if purchased alone, $255 with paper edition of Geology, $155 with paper editions of both Bulletin and Geology, add $1 2 for overseas, Call or e-mailfor information on 2-year CD-ROMsubscriptions. For member orders, call 1-888-443-4472 or (303) 447-2020, or email [email protected]. Nonmembers contact SubscriDtion Services at 1.800.472.1988. Claims: For nonreceiDt or damaaed copies,’ members contact Member Service Center, nonmembers contact Subscription Services. Claims are honored for one year; please alldw sufficient delivery time for overseas copies, up to six months. ADVERTISING: Contact Ann H. Crawford (acrawfordOgeosociety.org);phone (303) 447-2020 or fax 303-447-1 133. STAFF: Director of Publications, Peggy S. Lehr; Managing Editor, Jeanette Hammann; Editorial Assistant, Anika K. Burkard; Editorial Secretary, Marlene Mayer; Editorial Department: (303) 447-2020; Production, Gaynor Ann Bloom. GSA ONLINE: www.geosociety.org.