The Global Range of Subduction Zone Thermal Structures from Exhumed Blueschists and Eclogites: Rocks Are Hotter than Models

The maximum-pressure P–T conditions ( P _max –T ) and prograde P–T paths of exhumed subduction-related metamorphic rocks are compared to predictions of P–T conditions from computational thermal models of subduction systems. While the range of proposed models encompasses most estimated P _max –T conditions, models predict temperatures that are on average colder than those recorded by exhumed rocks. In general, discrepancies are greatest for P _max < 2 GPa, where only a few of the highest- T model paths overlap petrologic observations and model averages are 100–300 ◦C colder than average conditions recorded by rocks. Prograde P–T paths similarly indicate warmer subduction than typical models. Both petrologic estimates and models have inherent biases. Petrologic analysis may overestimate temperatures at P _max where overprinting occurs during exhumation, although P–T paths suggest that relatively warm conditions are experienced by rocks on the prograde subduction path. Models may underestimate temperatures at depth by neglecting shear heating, hydration reactions and fluid and rock advection. Our compilation and comparison suggest that exhumed high- P rocks provide a more accurate constraint on P–T conditions within subduction zones, and that those conditions may closely represent the subduction geotherm. While exhumation processes in subduction zones require closer petrologic scrutiny, the next generation of models should more comprehensively incorporate all sources of heat. Subduction-zone thermal structures from currently available models appear to be inaccurate, and this mismatch has wide-reaching implications for our understanding of global geochemical cycles, the petrologic structure of subduction zones, and fluid–rock interactions and seismicity within subduction zones.