The geodynamic origin of Los Humeros volcanic field in Mexico: insights from numerical simulations

Compared to normal arc-related volcanic eruptions, the formation of a volcanic caldera is a relatively atypical event. During caldera formation a series of large volumes of magma are erupted, reducing the structural support for the rock above the magma chamber and creating a large depression at the surface called caldera. Los Humeros volcanic field (LHVF) represents one of the largest volcanic calderas in Mexico. It is located some 400 km from the trench at the eastern edge of the Trans Mexican Volcanic Belt where the depth to the Cocos slab is more than 300 km. In this study we employ high-resolution two-dimensional thermomechanical numerical simulations of magma intrusions and a horizontal tectonic strain rate to better understand the influence of crustal deformation for the formation of Los Humeros caldera. A minimum number of three thermal anomaly pulses of hydrated mantle material (with diameter of 15 km or more) and a regional strain rate of 7.927 × 10^-16 s^-1 are required for magma to reach the surface. Modeling results show that regional extension coupled with deep thermal anomalies (with a temperature excess of ΔT ≥ 100 °C) that come in a specific chain-type sequence produce surface deformation patterns similar to LHVF. We propose an asthenospheric sub-slab deep source (> 300 km depth) for the thermal anomalies where previous studies showed the existence of a gap or tear in the Cocos slab.