A consistentmodel for fluid distribution, viscosity distribution, and flow-thermal structure in subduction zone

Water plays crucial roles in the subduction zone dynamics affecting the thermal-flow structure through the fluid processes. We aim to understand what controls the dynamics and construct a model to solve consistently fluid generation, fluid transport, its reaction with the solid and resultant viscosity, and thermal-flow structure. We highlight the effect of mechanical weakening of rocks associated with hydration. The viscosity of serpentinite (eta(serp)) in subduction zones critically controls the flow-thermal structure via extent of mechanical coupling between the subducting slab and overlying mantle wedge. When eta(serp) is greater than 10(21) Pa s, the thermal-flow structure reaches a steady state beneath the volcanic zone, and the melting region expands until C-in (initial water content in the subducting oceanic crust) reaches 3 wt %, and it does not expand from 3 wt %. On the other hand, when eta(serp) is less than 10(19) Pa s, the greater water dependence of viscosity (expressed by a larger f(v)) confines a hot material to a narrower channel intruding into the wedge corner from a deeper part of the back-arc region. Consequently, the overall heat flux becomes less for a larger f(v). When age(ba) (age of back-arc basin as a rifted lithosphere) = 7.5 Ma, the increase in fv weakens but shifts the melting region toward the trench side because of the narrow channel flow intruding into the wedge corner, where as it shuts down melting when age(ba) = 20 Ma. Several model cases (particularly those with eta(serp) = 10(20) to 10(21) Pa s and a relatively large f(v) for C-in = 2 to 3 wt %) broadly account for the observations in the Northeast Japan arc (i.e., location and width of volcanic chain, extent of serpentinite, surface heat flow, and seismic tomography), although the large variability of surface heat flow and seismic tomographic images does not allow us to constrain the parameter range tightly.