Thermal structure of a vanishing subduction system: An example of seismically-derived crustal temperature along the Italian peninsula

The active tectonic processes in convergent margins confer a high degree of complexity to the crust. Determining the thermal structure is, therefore, key to better elucidate the nature of those processes. In order to reconstruct the thermal structure of the crust beneath the Italian peninsula, we combine the most recent and accurate shear-wave velocity model that is currently available with thermodynamic modelling, assuming a global average crustal composition with no lateral variations. Our model, presented in terms of Moho temperature and crustal thermal gradients, shows a very good agreement with the known thermal anomalies associated with the backarc spreading related to the Apennine subduction. Importantly, we envisage a new anomalous region of high Moho temperatures in NW Italy (T > 800 °C at 30 km), at the transition between the Alps and Apennine orogens. The lowest temperatures of our model, corresponding to geothermal gradients <19 °C km-1, are obtained in the still active but slow-convergent portion of the northern Apennine. Moho temperatures increase moving southwards along the Apennine chain, an observation that is coherent with the evidence of ceasing subduction and consequent rebalancing of the depressed isotherms along the slab. Our results suggest that a thermal structure in different tectonic settings can be inferred with acceptable uncertainties based on absolute seismic velocity models. In this sense, our approach can be extended to any other region.