A finite-difference numerical program has been constructed to model the thermal structure of a subduction zone, including the combined effects of frictional heating and the heat sink created by dehydration of the subducted oceanic crust. Frictional heating is assumed to begin at the surface and to decrease linearly to zero at the depth at which dehydration begins. Dewatering of the oceanic crust occurs in the 80-125 km depth range and removes approximately 50 cal/gm of oceanic crust dehydrated. In young (<60 m.y. old) island arcs, isotherms are sufficiently elevated by frictional heating to melt the slab, perhaps accounting for the early tholeiitic volcanism recognized in some arcs. In old subduction zones, however, the frictional heating term must be increased to the equivalent of > 10 Kb shear stress (at the surface) before isotherms below the dehydration zone are elevated above the level predicted by conduction models. For lower stresses, the slab at $$\leq 150$$ km depth never becomes hotter than 600°C and thus does not melt. Alternative heat sources for arc volcanism are convection of asthenosphere above the slab and lowering of the melting point of overlying peridotite by release of water from dehydration of the subducted oceanic crust. In either case, the released water is probably a major transporter of Si, Na, K, Rb, Sr, Ba, and REE from the oceanic crust to the overlying wedge.
Read full abstract