Two-Dimensional Geothermal Model of the Peruvian Andes above the Nazca Ridge Subduction

The aseismic Nazca Ridge produces localized flat-slab subduction beneath the South American margin at latitudes 10◦ to 15◦ S. The geological evolution and the spatio-temporal pattern of deformation of the upper plate have been strongly influenced by the presence of the flat slab. In this study, we investigated the lithospheric thermal structure of this region by elaborating a 2D geothermal model along a section across the top of the Nazca Ridge, the Peru–Chile trench, the Andean Cordillera, and the Amazonian Basin, for a total length of 1000 km. For the sake of modelling, the crust of the overriding plate was subdivided into two parts, i.e., a sedimentary cover (including the entire lithostratigraphic sequence) and a crystalline basement. Applying an analytical methodology, we calculated geotherms and isotherms by setting (i) thickness, (ii) density, (iii) heat production, and (iv) thermal conductivity for each geological unit and considering (v) heat flux at the Moho, (vi) frictional heating produced by faults, and (vii) plate convergence rate. The resulting model could make a significant advance in our understanding of how flat slab geometry associated with the Nazca Ridge subduction affects the thermal structure and hence the tectonic evolution of the region.