Rheological evolution of a Mediterranean subduction complex

We use field and microstructural observations, coupled to previously published P-T-time histories, to track the rheological evolution of an intracontinental subduction complex exposed in the Betic Cordillera in the western Mediterranean region. The body of rock we focus on, known as the Nevado-Filabride Complex (NFC), was originally part of the upper crust of the Iberian margin. It was subducted into hot asthenospheric mantle, then exhumed back toward the surface in two stages: an early stage of fast exhumation along the top of the subducting slab in a subduction channel, and a late stage of slower exhumation resulting from capture by a low-angle detachment fault rooted at the brittle-ductile transition. Each stage of deformation in the NFC was punctuated by changes in the dominant deformation mechanism. Deformation during initial subduction of the complex was accommodated by pressure-solution creep in the presence of a fluid phase – the grain sizes, stress magnitudes, and estimated strain rates for this stage are most consistent with a thin-film model for pressure solution in which the diffusion length scale is controlled by the grain size. During the early stages of exhumation within the subduction channel, deformation transitioned from pressure solution to dislocation creep due to increases in temperature, which resulted in increases in both water fugacity and grain size, each of which favor the dislocation creep mechanism. Differential stress magnitudes for this stage were ～10 MPa, and are consistent with simple models of buoyancy-driven channel flow. With continuing subduction-channel exhumation, deformation remained within the dislocation creep field because sequestration of free water into hydrous, retrogressive minerals suppressed the pressure-solution mechanism. Differential stresses progressively increased to ～100 MPa near the mouth of the channel during cooling as the rocks moved into mid-crustal levels. During the final, core-complex stage of exhumation, deformation was progressively concentrated into a narrow zone of highly localized strain beneath a mid-crustal detachment fault. Localization was promoted by a transition from dislocation creep to dislocation-creep-accommodated grain boundary sliding at temperatures of ～350–380 °C, grain sizes of ～4 μm and differential stress magnitudes of ～200 MPa. Peak differential stress magnitudes of ～200 MPa recorded just below the brittle-ductile transition are consistent with Byerlee's law in the upper crust assuming a vertical maximum principal stress and near-hydrostatic pore fluid pressures. Overall, the distribution of stress with temperature, coupled to independent constraints on strain rate from field observations and geochronology, indicate that the naturally calibrated Hirth et al. (2001) flow law for wet quartzite accurately predicts the rheological behavior of mid-crustal rocks deforming by dislocation creep.