Gutscher Et Al., 2000
Gutscher Et Al., 2000
Marc-Andr Gutscher* Laboratoire GTS, UMR 5573, Universit Montpellier II, F-34095 Montpellier, France
Ren Maury IUEM/Universit Bretagne Occidentale, Place Nicolas Copernic, F-29280 Plouzan, France
Jean-Philippe Eissen
IRD Centre de Bretagne, B.P. 70, F-29280 Plouzan, France
Erwan Bourdon
ABSTRACT
Slab melting has been suggested as a likely source of adakitic arc magmas (i.e., andesitic
and dacitic magmas strongly depleted in Y and heavy rare earth elements). Existing numerical
and petrologic models, however, restrict partial melting to very young ( < 5 Ma) oceanic crust
(typically at 6080 km depth). Paradoxically, most of the known Pliocene-Quaternary adakite
occurrences are related to subduction of 1045 Ma lithosphere, which should not be able to melt
under normal subduction-zone thermal gradients. We propose an unusual mode of subduction
known as flat subduction, occurring in ~10% of the worlds convergent margins, that can produce 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
(<6 Ma) occurrences of adakitic magmas. Observations from Chile, Ecuador, and Costa Rica
suggest a three-stage evolution: (1) steep subduction produces a narrow calc-alkaline arc, typically ~300 km from the trench, above the asthenospheric wedge; (2) once flat subduction begins,
the lower plate travels several hundred kilometers at nearly the same depth, thus remaining in
a pressure-temperature window allowing slab melting over this broad distance; and (3) once flat
subduction continues for several million years, the asthenospheric wedge disappears, and a volcanic gap results, as in modern-day central Chile or Peru. The proposed hypothesis, which
reconciles thermal models with geochemical observations, has broad implications for the study
of arc magmatism and for the thermal evolution of convergent margins.
Keywords: flat subduction, thermal structure, slab melting, adakites.
INTRODUCTION
Petrologic models suggest that formation of
trondhjemite-tonalite-dacite (TTD) by partial
melting of the subducted slab was widespread
during Archean time, and these rocks represent a
major component of Precambrian gneiss terranes
(Martin, 1986; Drummond and Defant, 1990).
However, given the cooler mantle temperatures
during the Phanerozoic, these processes should
be uncommon today (Martin, 1986). Recent
models of magma genesis emphasize the role of
fluids released via dehydration reactions from the
subducting plate, thereby causing melting in the
overlying mantle wedge (Schmidt and Poli,
1998). In the past decade, numerous occurrences
of Cenozoic adakites (andesitic and dacitic magmas characterized by strong depletion in heavy
rare earth elements and high Sr/Y ratios) have
been documented and discussed as possible
examples of slab melting (Fig. 1) (Defant and
Drummond, 1990; Morris, 1995). These include
Mount St. Helens (Defant and Drummond, 1993)
and the recent major eruption of Mount Pinatubo
(Prouteau et al., 1999).
Numerical and petrological studies of pressuretemperature-time (P-T-t) paths in subduction
zones suggest that this process can only occur for
subduction of very young (5 Ma) lithosphere
(Peacock et al., 1994) (Fig. 2A, solid lines). Para*Current address: GEOMAR, Marine Geodynamics,
Wischhofstrasse 1-3, D-24148, Germany. E-mail:
[email protected].
535
us
v = 3 cm/a
= 0 MPa
4.0
steep slab
P-T-t paths
3.0
P (GPa)
0
2
10
Hb
50
flat slab
P-T-t paths
50
2.0
1.0
200
400
600
T (C)
800
150
1400
us
flat slab
(points from
Fig. 3B,C)
[km]
ou
Ecuador (0 Ma)
300
[km]
500
400
800
600
300
Moho
600
900
melting
1300
900
in "tongue"
1300
100
slab melting
asthenosphere
150
Ga
Am
400
600
T (C)
200
100
asthenosphere
1400
1400
1.0
200
asthenosphere
1400
50
Hb
Ec
partial melting in
asthenospheric
wedge
300
600
2.0
ut
1000
isotherms and P-T-t paths involved in flat subduction, slab melting can occur (Fig. 2). Prolonged flat subduction will sufficiently cool both
lithospheres such that the downgoing oceanic
crust dehydrates before the necessary temperatures for partial melting are attained and thus the
final stage, a volcanic gap, is achieved (Fig. 3D).
This complete evolution has been documented
for central Chile between 20 Ma and the present
(Kay and Abbruzzi, 1996). Ecuador is known to
have had a narrow calc-alkaline arc at 5 Ma
(Barberi et al., 1988) and has evolved to stage 2
today due to flattening of the slab in response to
subduction of Carnegie Ridge (Pennington,
1981; Gutscher et al., 1999b).
The Costa Rican margin also exhibits a similar magmatic evolution (Defant et al., 1992).
Before 6 Ma, magmatism was calc-alkaline in
nature and the subduction angle presumably
steep. Around 5 Ma uplift began in the Tala536
900
1300
1000
wet so
lid
steep slab
(points from
Fig. 3A)
3.0
P (GPa)
600
asthenosphere
t
ou
v = 6 cm/a
25 m.y. old lithosphere
Moho
900
4.0
1300
ou
300
Moh
100
Ga
Am
[km]
Ec
300
600
20
50
10
20
C
0
300
[km]
400
500
600
300
300
600
50
200
100
Moho
600
900
900
1300
1300
[km]
100
slab melting
asthenosphere
150
1400
200
100
300
[km]
400
no active volcanism
500
600
0
Moho
300
300
600
50
[km]
Moho
600
900
900
no asthenospheric wedge
1300
1300
100
asthenosphere
asthenosphere
150
1400
100
200
300
[km]
400
500
600
Figure 3.Three-stage evolution showing transition from steep to flat subduction style due to subduction of buoyant, overthickened oceanic crust (e.g., oceanic plateau), with impact on magmatic arc.Thermal structure is constructed based on adaptation of published numerical models
(Peacock et al., 1994; Peacock, 1996) and analytical solution (Davies, 1999) (see text). A: Steep or
normal subduction ( > 30 dip), with narrow calc-alkaline arc, ~300 km from trench produced by
slab dehydration and partial melting in overlying asthenospheric wedge. B: Early stage magmatic flat subduction; oceanic slab crosses 700 C geotherm at ~80 km depth (300400 km from
trench) allowing broad adakitic arc to develop. Partial melting continues in narrow intervening
tongue of asthenospheric material. C: Late stage magmatic flat subduction; intervening asthenospheric 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
eruption 4.7 Ma (Kay and Abbruzzi, 1996). D: Prolonged (amagmatic) flat subduction; oceanic
slab is cooler than 600 C, asthenospheric wedge has disappeared, and thus no partial melting
can occur, and volcanism ceases.
(Atherton and Petford, 1993; Petford and Atherton, 1996; Kay and Abbruzzi, 1996), and melting
of the base of an accreted arc or ophiolite complex,
as invoked for Ecuador (Arculus et al., 1999).
The slab window model in general applies to
trenchspreading center triple junctions, and
thus could succeed in explaining the Chile triple
junction and possibly the Panama occurrence
(assuming that a spreading center segment was
recently subducted). However, it is insufficient
to explain the numerous cases where moderately
old oceanic crust is subducting (e.g., Peru, central Chile, Japan, Alaska) and where no evidence
of lithospheric tears exists.
Melting at the base of an overthickened orogenic crust (garnet eclogite) has been suggested
for Chile (Kay and Abbruzzi, 1996), Peru (Petford
and Atherton, 1996), and Ecuador (Monzier
et al., 1997) based on isotopic imprints and trace
element features indicative of crustal components. However, in Ecuador the crustal imprint is
known to be weak (<15%) (Barragan et al.,
1998; Bourdon et al., 1999). Furthermore, it is
very difficult to distinguish whether these characteristics reflect melting of the base of the crust,
or slab melting followed by crustal contamination during ascent through the Andean crust
(e.g., the MASH model, Hildreth and Moorbath,
1988). Because these two processes often cannot
be distinguished geochemically, the reasons
commonly given for dismissing the slab-melting
model are thermal considerations, i.e., that crust
older than 20 Ma cannot melt (Petford and
Atherton, 1996; Kay and Abbruzzi, 1996).
The first argument against melting of the base
of a crust thickened by underplating of mafic
material is that, while this explanation may work
for parts of South America, it does not explain
Cascadia, southern Alaska, or southwestern
Japan, where continental crustal thicknesses do
not exceed 40 km (Oleskevich et al., 1999).
Moreover, the relatively high Mg, Cr, and Ni
contents of Andean adakites and related rocks
are best explained by interaction with mantle
peridotite (Sen and Dunn, 1995).
The strongest argument against crustal melting (either of continental or accreted oceanic
materials) is temporal. All well-studied adakitic
provinces are known from radiometric dating to
have been active for only a brief time span of
23 m.y. or less (e.g., Panama and Costa Rica,
Defant et al., 1992). Adakitic magmatism in central Chile ceased at 4.7 Ma (Kay and Abbruzzi,
1996), yet subduction beneath the thickened
crust continues today. The same is true for
Ecuador, where subduction has been occurring
and arc magmas have been ascending through
the accreted ophiolitic crust for at least 20 m.y.
Why is there no older (1020 Ma) adakitic magmatism and why does its onset correspond to the
change in subduction dip?
It is important to note that our flat subduction
slab melting model includes this temporal evolu537
tion. It is only during the early stages of flat subduction that the upper part and leading edge of the
slab can melt. Prolonged flat subduction will cool
both lithospheres and impede partial melting. It
has been repeatedly suggested that only the subduction of hot, very young lithosphere will allow
slab melting to occur (Defant and Drummond,
1993; Peacock et al., 1994). The emphasis here is
not on the initial temperature of the slab, but on
its buoyancy and thus on the time spent at a relatively shallow (~80 km) depth, allowing the slab
to heat. Therefore, oceanic crust as old as 50 Ma
can melt if it is given sufficient time to warm.
Young crust (510 Ma) subducts steeply at lat
42S in Chile, yet no adakites are observed between lat 42 and 44S.
ACKNOWLEDGMENTS
We thank the reviewers Richard Arculus and Marc
Hirschmann for constructive suggestions that helped
improve the original manuscript. M.-A. Gutschers research was supported by a Marie Curie Training and
Mobility of Researchers postdoctoral stipend from the
European Union.
REFERENCES CITED
Abbott, D., Drury, R., and Smith, W.H.F., 1994, Flat to
steep transition in subduction style: Geology,
v. 22, p. 937940.
Arculus, R.J., Lapierre, H., and Jaillard, E., 1999,
Geochemical window into subduction and accretion processes: Raspas metamorphic complex,
Ecuador: Geology, v. 27, p. 547551.
Atherton, M., and Petford, N., 1993, Generation of
sodium-rich magmas from newly underplated
basaltic crust: Nature, v. 362, p. 144146.
Barberi, F., Coltelli, M., Ferrara, G., Innocenti, F.,
Navarro, J.M., and Santacroce, R., 1988, PlioQuaternary volcanism in Ecuador: Geological
Magazine, v. 125, p. 114.
Barragan, R., Geist, D., Hall, M., Larson, P., and Kurz,
M., 1998, Subduction controls on the composition
of lavas from the Ecuadorian Andes: Earth and
Planetary Science Letters, v. 154, p. 153166.
Bourdon, E., Eissen, J.-P., Cotten, J., Monzier, M.,
Robin, C., and Hall, M.L., 1999, Les laves calcoalcalines et caractre adakitique du volcan Antisana (Equateur): Hypothse ptrogntique: Paris,
Academie des Sciences Comptes Rendus, v. 328,
p. 443449.
Cahill, T., and Isacks, B.L., 1992, Seismicity and shape
of the subducted Nazca plate: Journal of Geophysical Research, v. 97, p. 17,50317,529.
Davies, J.H., 1999, Simple analytic model for subduction zone thermal structure: Geophysical Journal
International, v. 139, p. 823828.
Davies, J.H., and Stevenson, D.J., 1992, Physical
model of source region of subduction zone volcanics: Journal of Geophysical Research, v. 97,
p. 20372070.
Defant, M.J., and Drummond, M.S., 1990, Derivation of
some modern arc magmas by melting of young
subducted lithosphere: Nature, v. 347, p. 662665.
Defant, M.J., and Drummond, M.S., 1993, Mount St.
Helens: Potential example of the partial melting
of the subducted lithosphere in a volcanic arc:
Geology, v. 21, p. 547550.
Defant, M.J., Jackson, T.E., Drummond, M.S., de Boer,
J.Z., Bellon, H., Feigenson, M.D., Maury, R.C.,
and Stewart, R.H., 1992, The geochemistry of
young volcanism throughout western Panama and
southeastern Costa Rica: An overview: Geological Society of London Journal, v. 149, p. 569579.
538
Printed in USA