Shear stresses on megathrusts: Implications for mountain building behind subduction zones

[ 1] Shear stresses tau on a subduction megathrust play an important role in determining the forces available for mountain building adjacent to a subduction zone. In this study, the temperatures and shear stresses on megathrusts in 11 subduction zones around the Pacific rim (Hikurangi, Tonga, Izu-Ogasawara, western Nankai, northeastern Japan, Aleutians, western Alaska, Cascadia, northern Chile, southern Chile) and SE Asia ( northern Sumatra) have been determined. The main constraint is that vertical normal stresses beneath the highlands behind the subduction zone are nearly equal to horizontal normal stresses, in the plane of a trench- or arc-normal section. For a typical brittle and ductile megathrust rheology, frictional shear stress tau = mu rho gz, for depth z, and ductile shear stress tau = A exp (B/RT) at temperature T, where mu, A, B are rheological parameters treated as constants. Rheological constants common to all the megathrusts ( m crust, m mantle, B) are determined by simultaneously solving for the force balance in the overlying wedge and megathrust thermal structure, using a simplex minimization algorithm, taking account of the induced mantle corner flow at depth (65 +/- 15 km (2 sigma)) and constant radiogenic heating (0.65 +/- 0.3 mu W m(-3) (2 sigma)) throughout the crust. The A constants are solved individually for each subduction zone, assuming that the maximum depth of interplate slip earthquakes marks the brittle-ductile transition. The best fit solution shows two groupings of megathrusts, with most subduction zones having a low mean shear stress in the range 7 - 15 MPa (mu(crust) = 0.032 +/- 0.006, mu(mantle) = 0.019 +/- 0.004) and unable to support elevations > 2.5 km. For a typical frictional sliding coefficient similar to 0.5, the low effective coefficients of friction suggest high pore fluid pressures at similar to 95% lithostatic pressure. Tonga and northern Chile require higher shear stresses with mu(crust) = 0.095 +/- 0.024, mu(mantle) = 0.026 +/- 0.007, suggesting slightly lower pore fluid pressures, at similar to 81% lithostatic. Ductile shear in the crust is poorly resolved but in the mantle appears to show a strong power law dependency, with B = 36 +/- 18 kJ mol(-1). A(mantle) values are sensitive to the precise value of B but are in the range 1 - 20 kPa. The power law exponent n for mantle flow is poorly constrained but is likely to be large ( n > 4). The brittle-ductile transition in the crust occurs at temperatures in the range 370 degrees C - 512 degrees C, usually close to the base of the crust and in the mantle at much lower temperatures ( 180 degrees C - 300 degrees C), possibly reflecting a marked change in pore fluid pressure or quasi ductile and subfrictional properties. In subduction zones where the subducted slab is older than 50 Ma, a significant proportion of the integrated shear force on the megathrust is taken up where it cuts the mantle and temperatures are <= 300 degrees C. In much younger subduction zones, the stress transmission is confined mainly to the crust. The shear stresses, particularly in the crust, may be kept low by some sort of lubricant such as abundant water-rich trench fill, which lowers the frictional sliding coefficient or effective viscosity and/or raises pore fluid pressure. The unusual high stress subduction zone in northern Chile lacks significant trench fill and may be poorly lubricated, with a mean shear stress similar to 37 MPa required to support elevations > 4 km in the high Andes. However, where the crust i thin in sediment-starved and poorly lubricated subduction zones, such as Tonga, the mean shear stress will still be low. Sediment may lubricate megathrusts accommodating underthrusting of continental crust, such as in the Himalayas or eastern central Andes, which have a low mean shear stress similar to 15 MPa.