LETTER doi:10.1038/nature09417
Melting above the anhydrous solidus controls the location of volcanic arcs
Philip C. England1 & Richard F. Katz1
Segregation of magma from the mantle in subduction zones is one (,120 6 40 km; refs 1, 3 and 5). The sharpness of the volcanic fronts of the principal mechanisms for chemical differentiation of the implies that a key process in the generation or transport of magma is Earth. Fundamental aspects of this system, in particular the pro- similarly focused beneath the arcs, but there is no consensus as to what cesses by which melt forms and travels to the Earth’s surface, remain that process may be. A wide range of metamorphic and melting re- obscure. Systematics in the location of volcanic arcs, the surface actions, either in the slab or in the mantle wedge, have been proposed expression of this melting, are widely considered to be a clue to as candidate processes. Some authors suggest that the arcs lie above processes taking place at depth, but many mutually incompatible places where the degree of melting in the mantle wedge becomes high interpretations of this clue exist (for example, see refs 1–6). We enough for the melt to segregate from the solid2,4. Others suggest that discriminate between those interpretations by the use of a simple the locations of arcs are determined by the release of fluid near the top scaling argument derived from a realistic mathematical model of of the slab in reactions that are either strongly pressure-dependent1,3,5 heat transfer in subduction zones. The locations of the arcs cannot or strongly temperature-dependent6,9. None of these suggestions has, be explained by the release of fluids in reactions taking place near however, produced a successful quantitative prediction of the location the top of the slab. Instead, the sharpness of the volcanic fronts, of volcanic arcs. Here we take a different approach: starting with the together with the systematics of their locations, requires that arcs observed correlation between the descent speed of the slab and its must be located above the place where the boundary defined by the depth beneath the volcanic arcs8, we use a simple mathematical model anhydrous solidus makes its closest approach to the trench. We to fit the data and reveal the petrological processes responsible for the show that heat carried by magma rising from this region is sufficient locations of the arcs. to modify the thermal structure of the wedge and determine the Although calculations of the full temperature field in subduction pathway through which both wet and dry melts reach the surface. zones require numerical models10, their results can be encapsulated in Volcanic arcs are characterized by sharp fronts whose locations may simple scaling relations that show11 that temperatures within the mantle be described, with misfits of no more than a few kilometres, by small wedge and at the top of the slab depend upon a single parameter, Vrd2/ circles on the Earth’s surface (Fig. 1 and refs 7 and 8); furthermore, the k. Here V is the convergence rate across the plate boundary, d is the dip depth of the top of the slab beneath these fronts falls in a narrow range of the slab, r isthe radial distance from the wedgecorner (Fig. 1), and k is
a b Trench Arc front
δ δ z 100 Overriding plate w
D 500 D R 1,000
T
Region enlarged in b Wedge 500 Subducting slab 1,1,000 000 V Zero-stress boundary Zero-stress
Figure 1 | Idealized cross-sections of a subduction zone, drawn wedge is at a depth zw. b, Enlargement of temperature structure of the rectangle perpendicular to the trench and the island arc. a, Two plates converge at a in a: Isotherms are shown at intervals of 100 uC. A schematic isotherm, labelled speed V, with the slab of oceanic lithosphere being subducted at an angle d T, has its closest approach to the wedge corner (its ‘nose’) immediately beneath beneath the overriding plate. The arc front is a zone a few kilometres wide, the volcanic front, at a distance R from the corner (black circle). An open circle across which volcanic activity begins as one moves away from the trench; it lies marks the top of the slab directly below the corner of this isotherm, at a depth D at a distance D above the top of the slab. The relative motion between slab and below the surface. The distance R cannot be determined by observation, but is overriding plate generates a creeping flow in the wedge of mantle between similar to RD 5 (D 2 zw)/sind, the radial distance of the top of the slab (open them, which follows stream lines, shown as curved lines. The corner of the circle) from the wedge corner.
1Department of Earth Sciences, South Parks Road, Oxford, OX1 3AN, UK.
700|NATURE|VOL467|7OCTOBER2010 ©2010 Macmillan Publishers Limited. All rights reserved LETTER RESEARCH
a b c 1,400 800 0 1,300 700 50 1,200 600 T 1,100 500 Moho 100 40 mm yr–1
1,000 (°C) 400 (°C) r s
T –1 T 60 mm yr 150
900 300 Depth (km) 80 mm yr–1 800 200 100 mm yr–1 200 700 100 600 0 250 0 25 50 75 100 125 0 25 50 75 100 125 0 50 100 150 200 250 (Vrδ2)/κ (Vrδ2)/κ Horizontal distance (km) Figure 2 | Scaling relations for temperatures in the core of the mantle the top of the slab) for convergence rates V of 40–100 mm yr21 and a slab dip of wedge, and at the top of the slab. a, Maximum temperature in the wedge, Tr,as 40u. c, The temperature structure near the wedge corner for a calculation with a a function of dimensionless distance from the wedge corner, Vrd2/k, where V is convergence speed V of 80 mm yr21 and a slab dip of 40u. The location of the convergence speed, d is dip of the slab, and k is thermal diffusivity. Dots with 500 uC isotherm is shown by green lines and the 1,225 uC and 1,275 uC error bars indicate the averages and standard deviations of Tr, determined from isotherms are shown by red lines. Green arrowheads show the horizontal extent calculations in which V varies from 10 mm yr21 to 100 mm yr21, in steps of over which some part of the oceanic crust is at a temperature of 500 uC. Red 10 mm yr21, while d varies from 20u to 70u in steps of 10u. The red line arrowheads show the range over which the maximum temperature in the corresponds to the theoretical expression for Tr (equation (1)), with mantle wedge lies between 1,225 uC and 1,275 uC (a typical range of T0 5 1,420 uC, B 5 3.3 and b 520.8. b, As for a, but for temperatures within temperature represented by the error bars in Fig. 3b). Uncertainties in the the slab. The thick blue line corresponds to the theoretical expression for Ts,the temperatures arising from idealizations in the model are discussed in section B2 temperature on the top of the slab (equation (2)), with C 5 1.4 and c 520.06. of the Supplementary Information. Thin blue lines indicate the temperatures at the base of the crust (7 km below the thermal diffusivity of the mantle. In this scaling, the maximum variation rules out the hypothesis that the arcs are located above the temperature in the mantle wedge Tr is given by: place where the top of the slab reaches a critical pressure correspond- "# 1,3,5 b ing to a single dehydration reaction . We may also rule out the Vrd2 T