TY - JOUR A1 - Toledo, T. A1 - Gaucher, E. A1 - Jousset, P. A1 - Jentsch, A. A1 - Haberland, C. A1 - Maurer, H. A1 - Krawczyk, C. A1 - Calò, M. A1 - Figueroa, A. T1 - Local Earthquake Tomography at Los Humeros Geothermal Field (Mexico) Y1 - 2020-12-07 VL - 125 IS - 12 JF - Journal of Geophysical Research: Solid Earth DO - 10.23689/fidgeo-4029 N2 - A passive seismic experiment using 25 broadband and 20 short‐period stations was conducted between September 2017 and September 2018 at Los Humeros geothermal field, an important natural laboratory for superhot geothermal systems in Mexico. From the recorded local seismicity, we derive a minimum 1‐D velocity model and obtain 3‐D Vp and Vp/Vs structures of Los Humeros. We improved the classical local earthquake tomography by using a postprocessing statistical approach. Several inversions were computed and averaged to reduce artifacts introduced by the model parametrization and to increase the resolution of the investigated region. Finally, the resulting Vp and Vp/Vs structures and associated seismicity were integrated with newly acquired geophysical and petrophysical data for comprehensive interpretation. The recorded seismicity is mainly grouped in three clusters, two of which seem directly related to exploitation activities. By combining new laboratory measurements and existing well data with our Vp model, we estimate possible geological unit boundaries. One large intrusion‐like body in the Vp model, together with neighboring high Vp/Vs anomalies, hints at a region of active resurgence or uplift due to the intrusion of new magma at the northern portion of the geothermal field. We interpret high Vp/Vs features as fluid bearing regions potentially favorable for further geothermal exploitation. Deep reaching permeable faults cutting the reservoir unit could explain fluid flow from a deeper local heat source in the area. N2 - Key Points: High‐quality earthquake data were collected to image the Vp and Vp/Vs models for the first time at Los Humeros geothermal field (Mexico). Inversions were performed by extending the classical earthquake tomography using a postprocessing statistical approach. Geological unit boundaries and fluid and gas bearing zones were interpreted considering new geological, geophysical, and petrophysical data. UR - http://resolver.sub.uni-goettingen.de/purl?gldocs-11858/8369 ER -