Marco Polcari

Geodetic constraints to the source mechanism of the 2011–2013 unrest at Campi Flegrei (Italy) caldera

Campi Flegrei caldera (Italy) was affected by a new unrest phase during 2011–2013. We exploit two COSMO-SkyMed data sets to map the deformation field, obtaining displacement rates reaching 9 cm/yr in 2012 in the caldera center. The resulting data set is fitted in a geophysical inversion framework using finite element forward models to account for the 3-D heterogeneous medium. The best fit model is a north dipping mixed-mode dislocation source lying at ~5 km depth. The driving mechanism is ascribable to magma input into the source of the large 1982–1984 unrest (since similar source characteristics were inferred) that generates initial inflation followed by additional shear slip accompanying the extension of crack tips. The history and the current state of the system indicate that Campi Flegrei is able to erupt again, and the advanced techniques adopted provide useful information for short-term forecasting.

3D displacement field retrieved by integrating Sentinel-1 InSAR and GPS data: the 2014 South Napa earthquake

Abstract In this study the integration of Sentinel-1 InSAR (Interferometric Synthetic Aperture Radar) and GPS (Global Positioning System) data was performed to estimate the three components of the ground deformation field due to the Mw 6.0 earthquake occurred on August 24th, 2014, in the Napa Valley, California, USA. The SAR data were acquired by the Sentinel-1 satellite on August 7th and 31st respectively. In addition, the GPS observations acquired during the whole month of August were analyzed. These data were obtained from the Bay Area Regional Deformation Network, the UNAVCO and the Crustal Dynamics Data Information System online archives. The data integration was realized by using a Bayesian statistical approach searching for the optimal estimation of the three deformation components. The experimental results show large displacements caused by the earthquake characterized by a predominantly NW-SE strike-slip fault mechanism.

X- and C-Band SAR Surface Displacement for the 2013 Lunigiana Earthquake (Northern Italy): A Breached Relay Ramp?

In this paper, we applied the differential interferometric synthetic aperture radar (DInSAR) technique to investigate and measure surface displacements due to the Mw 5.3 (M1 5.2), June 21, 2013 earthquake, occurred north of the Apuan Alps (NW Italy), in the discontinuity zone between the Lunigiana and Garfagnana area. Two differential interferograms showing the coseismic displacement have been generated using X-band and C-band data, taken from COSMO-SkyMed and RADARSAT-2 satellites, respectively. Both interferograms highlighted a clear pattern of subsidence of few cm located between the Lunigiana and Garfagnana basins. We then modeled the observed SAR deformation fields using the Okada analytical formulation and found them to be consistent with an extensional fault plane dipping toward NW at about 50 . The integrated analysis of DInSAR, geological data, modeling, and historical seismicity suggest that the fault responsible for the June 2013 earthquake corresponds to a breached relay ramp connecting the Lunigiana and Garfagnana seismogenic sources.

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 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.

Exploitation of SAR data to detect burned areas in the Sila mountain area (southern Italy)

This study focuses on testing the SAR coherence changes from Sentinel-1 data to detect burned areas and to compare the results with optical Sentinel-2 derived burned area product to be used as validation. Visible Infrared Imaging Radiometer Suite (VIIRS) data at 350 m resolution was used to identify active fires locations. We focused on a sequence of wildfires that affected the Sila mountain area during the summer of the 2017. This area of the Calabria region (southern Italy) was interested by a range of fires for the second half of July and the whole month of August ([1], [2]) due also to an extremely dry and hot summer. We used a pair of optical images acquired from Sentinel- 2 satellites on 24 July 2017 (pre-events) and 23 August 2017 (post-events). Firstly, we computed the Normalized Difference Vegetation Index (NDVI) for both images and calculated the difference between these two (dNDVI) at 10m resolution; the results put in evidence several areas characterized by vegetation reduction, with dNDVI values up to 0.3-0.4. Concerning the SAR data, we evaluated the coherence changes by exploiting two pairs of Sentinel-1 SAR data over the same area. Both pairs were acquired along descending orbit, respectively before (on July, 19th and 31st) and after (on September, 5th and 17th) the fires occurred in the Sila mountain area. The coherence was computed separately for the first (γpre) and the second pair (γpost) and the difference γpost - γpre was calculated. In this way, we evaluated the difference in coherence between September, i.e. post-fires, and July, i.e. pre-fires expecting a higher coherence after burning, due to the vegetation reduction. In several areas, the coherence seems to be consistent with the fire events showing increments up to 0.20-0.25. However, the increasing of coherence difference could also be due to other reasons such as the soil moisture variations in the proximity of lakes/rivers or the seasonal cultivation changes. Further analysis integrating more information such as the SAR amplitude signal and the cross-polarized backscattering coefficient will be conducted in order to better evaluate and discriminate any contributions.

An improved data integration algorithm to constrain the 3D displacement field induced by fast deformation phenomena tested on the Napa Valley earthquake

Abstract In this work, we propose an improved algorithm to constrain the 3D ground displacement field induced by fast surface deformations due to earthquakes or landslides. Based on the integration of different data, we estimate the three displacement components by solving a function minimization problem from the Bayes theory. We exploit the outcomes from SAR Interferometry (InSAR), Global Positioning System (GNSS) and Multiple Aperture Interferometry (MAI) to retrieve the 3D surface displacement field. Any other source of information can be added to the processing chain in a simple way, being the algorithm computationally efficient. Furthermore, we use the intensity Pixel Offset Tracking (POT) to locate the discontinuity produced on the surface by a sudden deformation phenomenon and then improve the GNSS data interpolation. This approach allows to be independent from other information such as in-situ investigations, tectonic studies or knowledge of the data covariance matrix. We applied such a method to investigate the ground deformation field related to the 2014 Mw 6.0 Napa Valley earthquake, occurred few kilometers from the San Andreas fault system.

A possibile breached relay ramp causing the 2013 Lunigiana earthquake (Northern Italy)

We applied the Differential Interferometric Synthetic Aperture Radar (DInSAR) technique to investigate and model the Mw 5.3 (Ml 5.2), June 21st, 2013 earthquake, occurred in the discontinuity zone between the Lunigiana and Garfagnana areas, NW Italy. We analyzed the coseismic displacement from two differential interferogram generated using COSMO-SkyMed X-band and RADARSAT2 C-band data. Both interferograms show a pattern of subsidence of a few cm located between the Lunigiana and Garfagnana basins. We then modeled the observed SAR deformation fields using the Okada analytical formulation and found them to be consistent with an extensional fault plane dipping towards NW at about 50°. Our results suggest that the fault responsible for the June 2013 earthquake may correspond to a breached relay ramp connecting the Lunigiana and Garfagnana seismogenic sources.

Interpretation of the urban areas subsidence detected by multitemporal DInSAR: the Cassino test site

Lucchi, Renata G. ... et. al.-- 87° Congresso della Societa Geologica Italiana e 90° Congresso della Societa Italiana di Mineralogia e Petrologia, The Future of the Italian Geosciences - The Italian Geosciences of the Future, 10-12 September 2014, Milan, Italy.-- 1 pageThe Montellina Spring (370 m a.s.l.) represents an example of groundwater resource in mountain region. It is a significant source of drinking water located in the right side of the Dora Baltea Valley (Northwestern Italy), SW of Quincinetto town. This spring shows a morphological location along a ridge, 400 m from the Renanchio Torrent in the lower sector of the slope. The spring was investigated using various methodologies as geological survey, supported by photo interpretation, structural reconstruction, NaCl and fluorescent tracer tests, discharge measurements. This multidisciplinary approach, necessary due to the complex geological setting, is required for the importance of the Montellina Spring. It is interesting in the hydrogeological context of Western Alps for its high discharge, relatively constant over time (average 150 l/s), and for its location outside a fluvial incision and suspended about 40 m above the Dora Baltea valley floor (Lasagna et al. 2013). According to the geological setting, the hydrogeological reconstruction of the area suggests that the large amount of groundwater in the basin is essentially favoured by a highly fractured bedrock, covered by wide and thick bodies of glacial and gravitational sediments. The emergence of the water along the slope, in the Montellina Spring, is essentially due to a change of permeability between the deep bedrock and the shallow bedrock and/or surficial sediments. The deep bedrock, showing closed fractures and/or fractures filled by glacial deposits, is slightly permeable. The shallow bedrock, strongly loosened as result of gravitational phenomena, and the local gravitational sediments are, on the contrary, highly permeable. The concentration of water at the spring is due to several reasons. a) The spring is immediately downward a detachment niche, dipping towards the spring, that essentially drains the water connected to the change of permeability in the bedrock. b) It is along an important fracture, that carries a part of the losses of the Renanchio Torrent. c) Finally, it is favored by the visible and buried morphology. Although it is located along a ridge, the spring occurs in a small depression between a moraine and a landslide body. It also can be favored by the likely concave trend of buried base of the landslide. At last, tracer tests of the Renanchio Torrent water with fluorescent tracer are performed, with a continuous monitoring in the Montellina Spring. The surveys permit to verify and quantify the spring and torrent hydrogeological relationship, suggesting that only a small fraction of stream losses feeds the spring.