The Pannonian basin is an extensional back-arc basin that has undergone neotectonic inversion and is currently shortening. The understanding and quantification of present-day deformation processes during this inversion are still incomplete. To this end, we investigate the active deformation of the Circum-Pannonian region via the interpolation of GNSS-derived velocity field and the derivation of the strain rate fields. For the interpolation of the velocity field, we use ordinary kriging, a strochastic interpolation method. Our results show that estimating a strain rate field that is virtually free of short-wavelength noise requires the scaling of the velocity uncertainties, i.e. assuming a minimum standard deviation of 1 mm/yr in our case. The deformation of the Circum-Pannonian region is defined by the 2–3 mm/yr, NNE-directed motion of the Dinarides, and by the 0.5–1.5 mm/yr, WSW to SSW directed motion of the eastern areas (European foreland, East Carpathians, South Carpathians, Transylvanian basin). These opposite-sense motions define a large-scale, on average NE-SW shortening and transpression-type deformation in the Dinarides as well as in the Pannonian basin, while the East and South Carpathians undergo regional N–S extension. Neotectonic structures generally show good agreement with the strain rate field, for example in the Dinarides, Eastern Alps, or in the western Pannonian basin. However, the presence of fault-parallel shortening or biaxial shortening along sinistral neotectonic structures in the central and eastern Pannonian basin show some discrepancy between current geodetic and observed neotectonic deformation. The vertical velocity field shows dominantly 100 and 1000 km wavelength signals, the former is probably related to the response of the Pannonian lithosphere-asthenosphere system to neotectonic basin inversion, while latter can possibly be explained by far-field subsidence patterns induced by the mantle response to melting of the Fennoscandian ice sheet during the current interglacial period.

We investigate ground movements induced by the 8 February 2016, Mw=4.2 earthquake at the Los Humeros Geothermal Field (Mexico) using Sentinel-1 radar interferometry. Previous estimated focal mechanism solution based on seismic data with a hypocentral depth of 1900 m could not resolve the measured coseismic surface deformation pattern. In this study, we applied inverse elastic dislocation models to estimate the source parameters of the seismic event. Our models suggest the reverse reactivation of the Los Humeros normal fault at a shallower depth (<1000 m), with a more significant left lateral component below ~400 m depth. The occurrence of such shallow events at Los Humeros pose increased risks for the neighboring communities and infrastructure. Therefore, continuous monitoring of seismicity and cautious planning of field operations are crucial.

A NNW-SSE striking fault swarm, including the Los Humeros fault, acts as a major boundary of the subsiding area observed by InSAR time-series between February 2016 and May 2019. A potential explanation of the reverse reactivation of the Los Humeros fault and following downward movement of the eastern fault block is the depressurization of the whole hydrothermal system. Such depressurization can occur due to the exploitation of the geothermal field and/or due to natural pressure/temperature changes related to magmatic activity.

On the 8th of February 2016, a Mw 4.2 earthquake was detected inside the Los Humeros caldera, located in the eastern sector of the Trans-Mexican Volcanic Belt. The event occurred after a sharp increase in the injection rate at the Los Humeros Geothermal Field and it was recorded by the seismic monitoring network of the power plant. The earthquake was felt by the local population and it caused damage in the power plant infrastructure. The focal mechanism solution of a previous study based on seismological data shows a reverse movement with a minor left-lateral component: Mw=4.2, depth=1500m, strike=169°, dip=61°, rake=42°. We have performed a geodetic and geomechanical analysis of the seismic source event based on ground deformation inferred from DInSAR. We used ascending and descending Sentinel-1 differential interferograms to retrieve the horizontal and vertical components of the co-seismic deformation. Subsequently, we inverted the estimated deformation to obtain the solution of an activated fault using the Okada model. These results shed light on the geomechanical aspects of the event and can help to understand the effects of field operations interacting with pre-existing structural features and active tectonic processes in the Los Humeros caldera.

Surface deformation due to fluid extraction can be detected by satellite-based geodetic sensors, providing important insights on subsurface geomechanical properties. In this study, we use Differential Interferometric Synthetic Aperture Radar (DInSAR) observations to measure ground deformation due to fluid extraction at the Los Humeros Geothermal Field (Puebla, Mexico). Our main goal is to reveal the pressure distribution in the reservoir and to identify reservoir compartmentalization, which can be important aspects for optimizing the production of the field. The result of the PS-InSAR (Persistent Scatterer by Synthetic Aperture Radar Interferometry) analysis shows that the subsidence at the LHGF was up to 8 mm/year between April 2003 and March 2007, which is small relative to the produced volume of 5×106 m3/year. The subsidence pattern indicates that the geothermal field is controlled by sealing faults separating the reservoir into several blocks. To assess if this is the case, we relate surface movements with volume changes in the reservoir through analytical solutions for different types of nuclei of strain. We constrain our models with the movements of the PS points as target observations. Our models imply small volume changes in the reservoir, and the different nuclei of strain solutions differ only slightly. These findings suggest that the pressure within the reservoir is well supported and that reservoir recharge is taking place.