Surface displacement of the Mw 7 Machaze earthquake (Mozambique): Complementary use of multiband InSAR and radar amplitude image correlation with elastic modelling.

International audienceIn this paper we investigate the surface displacement related to the 2006 Machaze earthquake using Synthetic Aperture Radar Interferometry (InSAR) and sub-pixel correlation (SPC) of radar amplitude images. We focus on surface displacement measurement during three stages of the seismic cycle. First, we examined the co-seismic stage, using an Advanced SAR (ASAR) sensor onboard the Envisat satellite. Then we investigated the post-seismic stage using the Phase Array L-band SAR sensor (PALSAR) onboard the ALOS satellite. Lastly, we focussed on the inter-seismic stage, prior to the earthquake by analysing the L20 and JERS-1 SAR data. The high degree of signal decorrelation in the C-band co-seismic interferogram hinders a correct positioning of the surface rupture and correct phase unwrapping. The post-seismic L-band interferograms reveal a time-constant surface displacement, causing subsidence of the surface at a ~ 5 cm/yr rate. This phenomenon continued to affect the close rupture field for at least two years following the earthquake and intrinsically reveals a candidate seismogenic fault trace that we use as a proxy for an inversion against an elastic dislocation model. Prior to the earthquake, the JERS interferograms do not indicate any traces of pre-seismic slip on the sismogenic fault. Therefore, slip after the earthquake is post seismic, and it was triggered by the Machaze earthquake. This feature represents a prominent post-seismic slip event rarely observed in such a geodynamic context

Inflation of the Aira Caldera (Japan) detected over Kokubu urban area using SAR interferometry ERS data

International audienceNine ERS-1 and ERS-2 descending orbit data acquired over Aira Caldera between June 1995 and November 1998 were used to create 36 differential interferograms. Although the interferograms exhibit a relatively low level of coherence, even for couples sampling short time intervals (6 months), Synthetique Aperture Radar (SAR) observations reveal a distinct range change pattern over Kokubu urban area whose amplitude increases with the time separation between SAR images. The analysis of the ground deformation time series relative to the earliest ERS images showed a maximum uplift of about 20 mm between the north and the south of the urban area during the period covered by our satellite observations. Taking the reduced surface of the coherent area into account, we performed a simple modeling of the deformation field assuming an inflating spherical source within an elastic half-space medium located beneath the centre of the Aira Caldera. This simple model predicts a maximum volume increase of 20–30×106 m3 between 1995 and 1998, which would produce an inflation of about 70 mm at the centre of the Aira Caldera and 40 mm in the Kokubu south urban area. These results are in good agreement with other geophysical observations carried out on the Aira caldera during this period. Despite the limited spatial extent of the coherent areas around the Aira Caldera, this study shows that DinSAR method using data collected in C band can be successfully used to detect subtle ground displacement changes of the volcanic complex and thus provides complementary information to ground-based geodetic monitoring of dynamic processes of the Aira Caldera and Sakurajima volcano

Crustal strain in central Greece from repeated GPS measurements in the interval 1989-1997

A 66-station GPS network spanning central Greece, first observed in 1989, has been occupied fully on three occasions: June 1989, October 1991 and May 1993. Subsets of this network bounding the Gulf of Korinthos have also been occupied in June 1995, October 1995, May 1996 and September/October 1997. The first three occupations were processed using a fiducial GPS methodology, whereas later surveys were processed using CODE precise orbits. Combination of data from different surveys to yield smooth site velocities requires global network translations at each epoch to compensate for errors in the realization of the reference frame. This method provides a posteriori estimates of the relative coordinate errors and reference frame noise. Only one earthquake, the 1995 June 15 Egion event, has caused significant local coseismic displacement, and its effects on the interseismic velocity field are removed using an elastic dislocation model. We constrain the orientation of the 100 yr triangulation—GPS velocity estimates of Davies et al. (1997) using 14 sites common to the two networks. The goodness of fit of this transformation indicates that the short-term and 100 yr geodetic estimates of deformation are highly compatible. We infer that short-term geodetic studies are capable of determining longer-term deformation rates provided that transient, local effects can be modelled. From the combined velocity field, we estimate principal strains and rigid-body rotation rates at points on a regular grid using data from neighbouring sites. Strain rates are high within the Gulf of Korinthos and much lower elsewhere. The extension rate across the Gulf of Korinthos increases from east to west. Comparison of the extension rate with historical and recent rates of seismic release of strain reveals significant medium-term seismic hazard in the western Gulf of Korinthos, and may also indicate long-term aseismic strai

The April 3, 2010 earthquake along the Pernicana fault (Mt. Etna - Italy): analysis of satellite and in situ ground deformation data integrated by the SISTEM approach

Etna is worldwide known as one of the most studied and monitored active volcanoes. Flank instability along the eastern and southern portion of Mt. Etna has been observed and measured thanks to geodetic networks and InSAR data analysis. The spreading area is bordered to the north by the east-west Pernicana Fault System (PFS) which dynamic is often linked with the eruptive activity, as recently observed during the 2002-2003 eruption. A seismic sequence occurred since April 2-3, 2010, along the PFS with two very shallow (a few hundred meters) mainshocks of magnitude 3.6 and 3.5. Explosions and ash emissions at the summit craters followed this swarm and culminated some days later (April 7-8). Just after the earthquake, specific GPS surveys were carried out aimed at monitoring the eastern part of the Pernicana fault, and the leveling route on the northeastern flank of the volcano was also surveyed. Trying to investigate the deformation occurred along the PFS during the events of April 3rd 2010, we performed a DInSAR (Differential Interferometric Synthetic Aperture Radar) analysis of ascending and descending Envisat, and of ascending ALOS-PALSAR images encompassing the date of the earthquake. The Envisat interferograms show very intense but local deformation on the Envisat ascending data and a low signal for the descending geometry, close to the Pernicana fault trace. This is probably due to the oblique normal/leftlateral kinematics of the PFS (as deduced also by GPS and leveling data), indeed both vertical (lowering) and horizontal (eastwards) components of motion produce a strong stretching of the LOS (Line Of Sight) distance for ascending geometry, while the two components act in opposite ways for the descending geometry, resulting in lower LOS distance variations compared to the ascending data set. We analyzed also the ALOS pair referring to 21/02/2010 – 08/04/2010 time and acquired along the ascending track number 638. The ALOS interferogram clearly show three fringes corresponding to roughly 35 cm of LOS displacement. The preliminary modeling of the interferograms agree with the seismic information (very shallow faulting, seismic moment) and show that the medium behave elastically. In order to investigate the ground deformation pattern associated with this event, an application of the novel SISTEM (Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements) approach is presented here. To achieve higher accuracy and get better constraint of the 3D components of the displacements, we improved the standard formulation of SISTEM approach, based on the GPS and a single DInSAR sensor, in order to take into account all the available dataset (GPS, leveling, ascending and descending ENVISAT C-Band interferograms and the ALOS L-Band data). The 3D displacement maps obtained using the SISTEM approach well show the kinematics of the PFS, and are able to reconstruct also the ground deformation affecting the whole investigated area, defining the movements of the north-eastern flank of the volcano. These results, which provide an accurate spatial characterization of ground deformation, are hence promising for future studies aimed at improving the knowledge about the kinematics of the active faults of Mt. Etna

Synergistic use of Lagrangian dispersion and radiative transfer modelling with satellite and surface remote sensing measurements for the investigation of volcanic plumes: the Mount Etna eruption of 25–27 October 2013

Abstract. In this paper we combine SO2 and ash plume dispersion modelling with satellite and surface remote sensing observations to study the regional influence of a relatively weak volcanic eruption from Mount Etna on the optical and micro-physical properties of Mediterranean aerosols. We analyse the Mount Etna eruption episode of 25–27 October 2013. The evolution of the plume along the trajectory is investigated by means of the FLEXible PARTicle Lagrangian dispersion (FLEXPART) model. The satellite data set includes true colour images, retrieved values of volcanic SO2 and ash, estimates of SO2 and ash emission rates derived from MODIS (MODerate resolution Imaging Spectroradiometer) observations and estimates of cloud top pressure from SEVIRI (Spinning Enhanced Visible and InfraRed Imager). Surface remote sensing measurements of aerosol and SO2 made at the ENEA Station for Climate Observations (35.52° N, 12.63° E; 50 m a.s.l.) on the island of Lampedusa are used in the analysis. The combination of these different data sets suggests that SO2 and ash, despite the initial injection at about 7.0 km altitude, reached altitudes around 10–12 km and influenced the column average aerosol particle size distribution at a distance of more than 350 km downwind. This study indicates that even a relatively weak volcanic eruption may produce an observable effect on the aerosol properties at the regional scale. The impact of secondary sulfate particles on the aerosol size distribution at Lampedusa is discussed and estimates of the clear-sky direct aerosol radiative forcing are derived. Daily shortwave radiative forcing efficiencies, i.e. radiative forcing per unit AOD (aerosol optical depth), are calculated with the LibRadtran model. They are estimated between −39 and −48 W m−2 AOD−1 at the top of the atmosphere and between −66 and −49 W m−2 AOD−1 at the surface, with the variability in the estimates mainly depending on the aerosol single scattering albedo. These results suggest that sulfate particles played a large role in the transported plume composition and radiative forcing, while the contribution by ash particles was small in the volcanic plume arriving at Lampedusa during this event

GPS network monitor the Western Alps deformation over a five year period: 1993-1998

GPS surveys in the Western Alps, performed in the time span 1993-2003, estimated the current crustal deformation of this area.Published63-763.2. Tettonica attivaJCR Journalreserve

Domino-style earthquakes along blind normal faults in Northern Thessaly (Greece): kinematic evidence from field observations, seismology, SAR interferometry and GNSS

Here we present a joint analysis of the geodetic, seismological and geological data of the March 2021 Northern Thessaly seismic sequence, that were gathered and processed as of April 30, 2021. First, we relocated seismicity data from regional and local networks and inferred the dip-direction (NE) and dip-angle (38°) of the March 3, 2021 rupture plane. Furthermore, we used ascending and descending SAR images acquired by the Sentinel-1 satellites to map the co-seismic displacement field. Our results indicate that the March 3, 2021 Mw=6.3 rupture occurred on a NE-dipping, 39° normal fault located between the villages Zarko (Trikala) and Damasi (Larissa). The event of March 4, 2021 occurred northwest of Damasi, along a fault oriented WNW-ESE and produced less deformation than the event of the previous day. The third event occurred on March 12, 2021 along a south-dipping normal fault. We computed 22 focal mechanisms of aftershocks with M≥4.0 using P-wave first motion polarities. Nearly all focal mechanisms exhibit normal kinematics or have a dominant normal dip-slip component. The use of InSAR was crucial to differentiate the ground deformation between the ruptures. The majority of deformation occurs in the vertical component, with a maximum of 0.39 m of subsidence over the Mw=6.3 rupture plane, south and west of Damasi. A total amount of 0.3 m horizontal displacement (E-W) was measured. We also used GNSS data (at 30-s sampling interval) from twelve permanent stations near the epicentres to obtain 3D seismic offsets of station positions. Only the first event produces significant displacement at the GNSS stations (as predicted by the fault models, themselves very well constrained by InSAR). We calculated several post-seismic interferograms, yet we have observed that there is almost no post-seismic deformation, except in the footwall area (Zarkos mountain). This post-seismic deformation is below the 7 mm level (quarter of a fringe) in the near field and below the 1 mm level at the GNSS sites. The cascading activation of the three events in a SE to NW direction points to a pattern of domino-style earthquakes, along neighbouring fault segments. The kinematics of the ruptures point to a counter-clockwise change in the extension direction of the upper crust (from NE-SW near Damasi to N-S towards northwest, near Verdikoussa)

The Corinth Rift Laboratory, Greece (CRL): A Multidisciplinary Near Fault Observatory (NFO) on a Fast Rifting System

The western rift of Corinth (Greece) is one of the most active tectonic structures of the euro-mediterranean area. Its NS opening rate is 1.5 cm/yr ( strain rate of 10-6/yr) results into a high microseismicity level and a few destructive, M>6 earthquakes per century, activating a system of mostly north dipping normal faults. Since 2001, monitoring arrays of the European Corinth Rift Laboratory (CRL, www.crlab.eu) allowed to better track the mechanical processes at work, with short period and broad band seismometers, cGPS, borehole strainmeters, EM stations, …). The recent (300 kyr) tectonic history has been revealed by onland (uplifted fan deltas and terraces) and offshore geological studies (mapping, shallow seismic, coring), showing a fast evolution of the normal fault system. The microseismicity, dominated by swarms lasting from days to months, mostly clusters in a layer 1 to 3 km thick, between 6 and 9 km in depth, dipping towards north, on which most faults are rooting. The diffusion of the microseismicity suggests its triggering by pore pressure transients, with no or barely detected strain. Despite a large proportion of multiplets, true repeaters seem seldom, suggesting a minor contribution of creep in their triggering, although transient or steady creep is clearly detected on the shallow part of some majors faults. The microseismic layer may thus be an immature, downward growing detachment, and the dominant rifting mechanism might be a mode I, anelastic strain beneath the rift axis , for which a mechanical model is under development. Paleoseismological (trenching, paleoshorelines, turbidites), archeological and historical studies completed the catalogues of instrumental seismicity, motivating attempts of time dependent hazard assessment. The Near Fault Observatory of CRL is thus a multidisciplinary research infrastructure aiming at a better understanding and modeling of multiscale, coupled seismic/aseismic processes on fault systems.Grant for Researchers (CC) ID 188753

GPS constraints on deformation in northern Central America from 1999 to 2017, Part 1 – Time-dependent modelling of large regional earthquakes and their post-seismic effects

We use continuous and campaign measurements from 215 GPS sites in northern Central America and southern Mexico to estimate coseismic and afterslip solutions for the 2009 Mw = 7.3 Swan Islands fault strike-slip earthquake and the 2012 Mw = 7.3 El Salvador and Mw = 7.4 Guatemala thrust-faulting earthquakes on the Middle America trench. Our simultaneous, time-dependent inversion of more than 350 000 daily GPS site positions gives the first jointly consistent estimates of the coseismic slips for all three earthquakes, their combined time-dependent post-seismic effects and secular station velocities corrected for both the coseismic and post-seismic deformation. Our geodetic slip solutions for all three earthquakes agree with previous estimates that were derived via static coseismic-offset modelling. Our time-dependent model, which attributes all transient post-seismic deformation to earthquake afterslip, fits nearly all of the continuous GPS site position time-series within their severalmillimetre position noise. Afterslip moments for the three earthquakes range from 35 to 140 per cent of the geodetic coseismic moments, with the largest afterslip estimated for the 2012 El Salvador earthquake along the weakly coupled El Salvador trench segment. Forward modelling of viscoelastic deformation triggered by all three earthquakes for a range of assumed mantle and lower crustal viscosities suggests that it accounts for under 20 per cent of the observed post-seismic deformation and possibly under 10 per cent