Matteo Albano

Land subsidence, Ground Fissures and Buried Faults: InSAR Monitoring of Ciudad Guzmán (Jalisco, Mexico)

We study land subsidence processes and the associated ground fissuring, affecting an active graben filled by thick unconsolidated deposits by means of InSAR techniques and fieldwork. On 21 September 2012, Ciudad Guzman (Jalisco, Mexico) was struck by ground fissures of about 1.5 km of length, causing the deformation of the roads and the propagation of fissures in adjacent buildings. The field survey showed that fissures alignment is coincident with the escarpments produced on 19 September 1985, when a strong earthquake with magnitude 8.1 struck central Mexico. In order to detect and map the spatio-temporal features of the processes that led to the 2012 ground fissures, we applied InSAR multi-temporal techniques to process ENVISAT-ASAR and RADARSAT-2 satellite SAR images acquired between 2003 and 2012. We detect up to 20 mm/year of subsidence of the northwestern part of Ciudad Guzman. These incremental movements are consistent with the ground fissures observed in 2012. Based on interferometric results, field data and 2D numerical model, we suggest that ground deformations and fissuring are due to the presence of areal subsidence correlated with variable sediment thickness and differential compaction, partly driven by the exploitation of the aquifers and controlled by the distribution and position of buried faults.

Gravity-driven postseismic deformation following the Mw 6.3 2009 L’Aquila (Italy) earthquake

The present work focuses on the postseismic deformation observed in the region of L’Aquila (central Italy) following the Mw 6.3 earthquake that occurred on April 6, 2009. A new, 16-month-long dataset of COSMO-SkyMed SAR images was analysed using the Persistent Scatterer Pairs interferometric technique. The analysis revealed the existence of postseismic ground subsidence in the mountainous rocky area of Mt Ocre ridge, contiguous to the sedimentary plain that experienced coseismic subsidence. The postseismic subsidence was characterized by displacements of 10 to 35 mm along the SAR line of sight. In the Mt Ocre ridge, widespread morphological elements associated with gravitational spreading have been previously mapped. We tested the hypothesis that the postseismic subsidence of the Mt Ocre ridge compensates the loss of equilibrium induced by the nearby coseismic subsidence. Therefore, we simulated the coseismic and postseismic displacement fields via the finite element method. We included the gravitational load and fault slip and accounted for the geometrical and rheological characteristics of the area. We found that the elastoplastic behaviour of the material under gravitational loading best explains the observed postseismic displacement. These findings emphasize the role of gravity in the postseismic processes at the fault scale.

Subsidence Detected by Multi-Pass Differential SAR Interferometry in the Cassino Plain (Central Italy): Joint Effect of Geological and Anthropogenic Factors?

In the present work, the Differential SAR Interferometry (DInSAR) technique has been applied to study the surface movements affecting the sedimentary basin of Cassino municipality. Two datasets of SAR images, provided by ERS 1-2 and Envisat missions, have been acquired from 1992 to 2010. Such datasets have been processed independently each other and with different techniques nevertheless providing compatible results. DInSAR data show a subsidence rate mostly located in the northeast side of the city, with a subsidence rate decreasing from about 5–6 mm/yr in the period 1992–2000 to about 1–2 mm/yr between 2004 and 2010, highlighting a progressive reduction of the phenomenon. Based on interferometric results and geological/geotechnical observations, the explanation of the detected movements allows to confirm the anthropogenic (surface effect due to building construction) and geological causes (thickness and characteristics of the compressible stratum).

New geological data on the Cassino intermontane basin, central Apennines, Italy

This paper presents the outcome of photogeological analysis, field survey as well as borehole correlation, to add new data concerning the intermontane Cassino basin and define facies and thickness of its Quaternary deposits. The investigated area is located in the end sector of the Latina Valley, Southern Lazio and belongs to the Lazio–Abruzzi domain. With a maximum extension of 250 km2 and a NW–SE trend, this paleo-lake was created by lava flows from the Roccamonfina volcano, which barred the Latina Valley. The lacustrine facies were described during the 1960s but available bibliographic data about its spatial distribution and its relationship with the Mesozoic–Cenozoic bedrock are scarce. The thickness of Quaternary deposits has been deeply influenced by inherited paleo-morphologies. Indeed the geological reconstruction of the basin has outlined the presence of an underlying horst and graben structure. This latter obscured the contractional tectonics originated by the orogenetic deformation, even if new structures have been recognized at the foot of Montecassino Hill. We infer that a key role has been performed by the Roveto Valley–Atina–Caserta fault, between Atina and San Pietro Infine.

Discriminating between natural and anthropogenic earthquakes: insights from the Emilia Romagna (Italy) 2012 seismic sequence

The potential for oilfield activities to trigger earthquakes in seismogenic areas has been hotly debated. Our model compares the stress changes from remote water injection and a natural earthquake, both of which occurred in northern Italy in recent years, and their potential effects on a nearby Mw 5.9 earthquake that occurred in 2012. First, we calculate the Coulomb stress from 20 years of fluid injection in a nearby oilfield by using a poroelastic model. Then, we compute the stress changes for a 2011 Mw 4.5 earthquake that occurred close to the area of the 2012 mainshock. We found that anthropogenic activities produced an effect that was less than 10% of that generated by the Mw 4.5 earthquake. Therefore, the 2012 earthquake was likely associated with a natural stress increase. The probability of triggering depends on the magnitude of recent earthquakes, the amount of injected water, the distance from an event, and the proximity to the failure of the activated fault. Determining changes that are associated with seismic hazards requires poroelastic area-specific models that include both tectonic and anthropogenic activities. This comprehensive approach is particularly important when assessing the risk of triggered seismicity near densely populated areas.

Aftershocks, groundwater changes and postseismic ground displacements related to pore pressure gradients: Insights from the 2012 Emilia-Romagna earthquake

During the 2012 Emilia-Romagna (Italy) seismic sequence, several time-dependent phenomena occurred, such as changes in the groundwater regime and chemistry, liquefaction, and postseismic ground displacements. Because time-dependent phenomena require time-dependent physical mechanisms, we interpreted such events as the result of the poroelastic response of the crust after the mainshock. In our study, we performed a two-dimensional poroelastic numerical analysis calibrated with Cosmo-SkyMed interferometric data and measured piezometric levels in water wells. The simulation results are consistent with the observed postseismic ground displacement and water level changes. The simulations show that crustal volumetric changes induced by poroelastic relaxation and the afterslip along the mainshock fault are both required to reproduce the amplitude (approximately 4 cm) and temporal evolution of the observed postseismic uplift. Poroelastic relaxation also affects the aftershock distribution. In fact, the aftershocks are correlated with the postseismic Coulomb stress evolution. In particular, a considerably higher fraction of aftershocks occurs when the evolving poroelastic Coulomb stress is positive. These findings highlight the need to perform calculations that adequately consider the time-dependent poroelastic effect when modeling postseismic scenarios, especially for forecasting the temporal and spatial evolution of stresses after a large earthquake. Failing to do so results in an overestimation of the afterslip and an inaccurate definition of stress and strain in the postseismic phase.

New insights into earthquake precursors from InSAR

We measured ground displacements before and after the 2009 L’Aquila earthquake using multi-temporal InSAR techniques to identify seismic precursor signals. We estimated the ground deformation and its temporal evolution by exploiting a large dataset of SAR imagery that spans seventy-two months before and sixteen months after the mainshock. These satellite data show that up to 15 mm of subsidence occurred beginning three years before the mainshock. This deformation occurred within two Quaternary basins that are located close to the epicentral area and are filled with sediments hosting multi-layer aquifers. After the earthquake, the same basins experienced up to 12 mm of uplift over approximately nine months. Before the earthquake, the rocks at depth dilated, and fractures opened. Consequently, fluids migrated into the dilated volume, thereby lowering the groundwater table in the carbonate hydrostructures and in the hydrologically connected multi-layer aquifers within the basins. This process caused the elastic consolidation of the fine-grained sediments within the basins, resulting in the detected subsidence. After the earthquake, the fractures closed, and the deep fluids were squeezed out. The pre-seismic ground displacements were then recovered because the groundwater table rose and natural recharge of the shallow multi-layer aquifers occurred, which caused the observed uplift.

Did Anthropogenic Activities Trigger the 3 April 2017 Mw 6.5 Botswana Earthquake

On 3 April 2017, a Mw 6.5 earthquake occurred in Botswana, representing the second-strongest earthquake registered since 1949. Such an intraplate event occurred in a low seismic hazard area and was suspected to be an artificial earthquake induced by nearby anthropogenic activities (gas extraction). The possible relation between anthropogenic activities and the earthquake occurrence has been qualitatively investigated. We estimated the geometric and kinematic characteristics of the causative fault from the modeling of Sentinel-1 InSAR interferograms. Our best-fit solution for the main shock is represented by a normal fault located at a depth greater than 20 km, dipping 65° northeast, with a right-lateral component, and a mean slip of 2.7 m. The retrieved fault geometry and mechanism are incompatible with the hypothetical stress perturbation caused by the anthropogenic activities performed in the area. Therefore, the 3 April 2017 Botswana earthquake can be classified as a natural intraplate earthquake.

Application and analysis of geodetic protocols for monitoring subsidence phenomena along on-shore hydrocarbon reservoirs

Abstract In this study, we tested the “land-subsidence monitoring guidelines” proposed by the Italian Ministry of Economic Development (MISE), to study ground deformations along on-shore hydrocarbon reservoirs. We propose protocols that include the joint use of Global Positioning System (GPS) and multi-temporal Differential Interferometric Synthetic Aperture Radar (DInSAR) techniques, for a twofold purpose: a) monitoring land subsidence phenomena along selected areas after defining the background of ground deformations; b) analyzing possible relationships between hydrocarbon exploitation and anomalous deformation patterns. Experimental results, gathered along the Ravenna coastline (northern Italy) and in the southeastern Sicily (southern Italy), show wide areas of subsidence mainly related to natural and anthropogenic processes. Moreover, ground deformations retrieved through multi-temporal DInSAR time series exhibit low sensitivity as well as poor spatial and temporal correlation with hydrocarbon exploitation activities. Results allow evaluating the advantages and limitations of proposed protocols, to improve the techniques and security standards established by MISE guidelines for monitoring on-shore hydrocarbon reservoirs.

Cosesimic liquefaction phenomena from DInSAR after the May 20, 2012 Emilia earthquake

In this paper, we have investigated the capability of Differential Interferometry Synthetic Aperture Radar (DInSAR) technique to detect the ground effects induced by liquefaction phenomena occurred during the May 20, 2012 Emilia earthquake. To this aim, a set of COSMO-SkyMed (CSK) SAR images covering the coseismic phase has been used. The detected surface effects have been related to liquefaction of deep sandy layers. Thanks to the geological/geotechnical data in the area and a liquefaction susceptibility analysis of the subsoil, it has been identified a sandy layer between 9 and 13 m in deep, which probably liquefied during the earthquake. The estimated vertical displacements due to liquefaction are comparable with the values measured by DInSAR.