Explore the vast range of titles available.

We study the tectonic deformation from the February 2017 shallow earthquake sequence onshore Biga Peninsula (NW Turkey, NE Aegean region). We use InSAR interferograms (Sentinel-1 satellites) to identify the seismic fault (striking N110°E) and seismological data (parametric data and Moment Tensor solutions from NOA and KOERI catalogues) so as to refine its geometry and kinematics using inversion techniques. Despite the moderate magnitudes of the main events of the sequence (5.0 ≤ Mw ≤ 5.2), the total surface deformation is 2.2 fringes (or maximum 6.2 cm along LOS) and it is well visible with InSAR because of the shallow depth of the four main events (6–8 km) and the good coherence of the signal phase. Our geodetic inversion showed that the fault has normal-slip kinematics, dimensions of 6 by 6 km (length, width) and dips at 45°. The InSAR data are fitted by a uniform slip of 28 cm. In addition, 429 earthquakes were relocated with the HypoDD software and the use of a 1-D velocity model. The dip-direction of the fault is not retrievable from InSAR, but a south-dipping plane is clear from seismology and the aftershocks distribution. The spatial distribution of relocated events indicates the activation of one fault with a rupture zone length of about 10 km, a result of the occurrence of off-fault aftershocks along strike the main rupture. A stress inversion using 20 focal mechanisms (M ≥ 3.6; NOA solutions) indicates that faulting accommodates a N196°E extension. It is confirmed that moderate (5.0 ≤ M ≤ 5.2) shallow events can be traced in InSAR studies and can produce surface displacements that provide useful data in fault inversion.

We identify the source of the Mw = 6.4 earthquake that rocked north-central Albania on November 26, 2019 02:54 UTC. We use synthetic aperture radar interferograms tied to the time series of coordinates of two permanent Global Navigation Satellite System (GNSS) stations (DUR2 and TIR2). We model the source by inverting the displacement data. Assuming in our model a half-space elastic medium and uniform slip along a rectangular fault surface, we invert the 104 picked measurements on a couple of ascending and descending interferograms to calculate the parameters of the fault. All inversions made with different input parameters converge towards a stable and robust solution with root mean square (r.m.s.) residual of 5.4 mm, thus ~1/5 of a fringe. They reveal that the earthquake occurred deep in the crust on a low-angle fault (23°) dipping towards east with a centroid at 16.5 km depth. The best-fitting length and width of the fault are 22 and 13 km, and the reverse slip, 0.55 m. The seismic moment deduced from our model agrees with those of the published seismic moment tensors. This geometry is compatible with a blind thrust fault that may root on the main basal thrust, i.e., along the thrust front that separates Adria–Apulia from Eurasia. It is notable that there is a 123 ns yr−1 active shortening of the crust between the GNSS stations DUR2-TIR2 (equivalent to a shortening rate of 3.6 mm yr−1), and roughly in the east–west direction. Given this amount of strain the recurrence time of M6+ earthquakes along this fault should be of the order of 150 years.

The aerosol properties of Mount Etna’s passive degassing plume and its short-term processes and radiative impact were studied in detail during the EPL-RADIO campaigns (summer 2016–2017), using a synergistic combination of observations and radiative transfer modelling. Summit observations show extremely high particulate matter concentrations. Using portable photometers, the first mapping of small-scale (within ∼ 20 km from the degassing craters) spatial variability of the average size and coarse-to-fine burden proportion of volcanic aerosols is obtained. A substantial variability of the plume properties is found at these spatial scales, revealing that processes (e.g. new particle formation and/or coarse aerosols sedimentation) are at play, which are not represented with current regional scale modelling and satellite observations. Statistically significant progressively smaller particles and decreasing coarse-to-fine particles burden proportion are found along plume dispersion. Vertical structures of typical passive degassing plumes are also obtained using observations from a fixed LiDAR station constrained with quasi-simultaneous photometric observations. These observations are used as input to radiative transfer calculations, to obtain the shortwave top of the atmosphere (TOA) and surface radiative effect of the plume. For a plume with an ultraviolet aerosol optical depth of 0.12–0.14, daily average radiative forcings of - 4.5 and - 7.0 W/m 2 , at TOA and surface, are found at a fixed location ∼ 7 km downwind the degassing craters. This is the first available estimation in the literature of the local radiative impact of a passive degassing volcanic plume.

The 2015/11/17 Lefkada (Greece) earthquake ruptured a segment of the Cephalonia Transform Fault (CTF) where probably the penultimate major event was in 1948. Using near-source strong motion and high sampling rate GPS data and Sentinel-1A SAR images on two tracks, we performed the inversion for the geometry, slip distribution and rupture history of the causative fault with a three-step self-consistent procedure, in which every step provided input parameters for the next one. Our preferred model results in a ~70° ESE-dipping and ~13° N-striking fault plane, with a strike-slip mechanism (rake ~169°) in agreement with the CTF tectonic regime. This model shows a bilateral propagation spanning ~9 s with the activation of three main slip patches, characterized by rise time and peak slip velocity in the ranges 2.5–3.5 s and 1.4–2.4 m/s, respectively, corresponding to 1.2–1.8 m of slip which is mainly concentrated in the shallower (<10 km) southern half of the causative fault. The inferred slip distribution and the resulting seismic moment (M 0  = 1.05 × 10 19 N m) suggest a magnitude of M w 6.6. Our best solution suggests that the occurrence of large ( M w  > 6) earthquakes to the northern and to the southern boundaries of the 2015 causative fault cannot be excluded.

Synthetic Aperture Radar Interferometry (InSAR) is a space geodetic technique used for mapping deformations of the Earth’s surface. It has been developed and used increasingly during the last thirty years to measure displacements produced by earthquakes, volcanic activity and other crustal deformations. A limiting factor to this technique is the effect of the troposphere, as spatial and temporal variations in temperature, pressure, and relative humidity introduce significant phase delays in the microwave imagery, thus “masking” surface displacements due to tectonic or other geophysical processes. The use of Numerical Weather Prediction (NWP) models as a tropospheric correction method in InSAR can tackle several of the problems faced with other correction techniques (such as timing, spatial coverage and data availability issues). High-resolution tropospheric modelling is particularly useful in the case of single interferograms, where the removal of the atmospheric phase screen (and especially the highly variable turbulent component) can reveal large-amplitude deformation signals (as in the case of an earthquake). In the western Gulf of Corinth, prominent topography makes the removal of both the stratified and turbulent atmospheric phase screens a challenging task. Here, we investigate the extent to which a high-resolution WRF 1-km re-analysis can produce detailed tropospheric delay maps of the required accuracy by coupling its output (in terms of Zenith Total Delay or ZTD) with the vertical delay component in GNSS measurements. The model is operated with varying physical parameterization in order to identify the best configuration, and validated with GNSS zenithal tropospheric delays, providing a benchmark of real atmospheric conditions. We correct sixteen Sentinel-1A interferograms with differential delay maps at the line-of-sight (LOS) produced by WRF re-analysis. In most cases, corrections lead to a decrease in the phase gradient, with average root-mean-square (RMS) and standard deviation (SD) reductions in the wrapped phase of 6.0% and 19.3%, respectively. Results suggest a high potential of the model to reproduce both the long-wavelength stratified atmospheric signal and the short-wave turbulent atmospheric component which are evident in the interferograms.

In the last thirty years, Synthetic Aperture Radar interferometry (InSAR) and the Global Navigation Satellite System (GNSS) have become fundamental space geodetic techniques for mapping surface deformations due to tectonic movements. One major limiting factor to those techniques is the effect of the troposphere, as surface velocities are of the order of a few mm yr−1, and high accuracy (to mm level) is required. The troposphere introduces a path delay in the microwave signal, which, in the case of GNSS Precise Point Positioning (PPP), can nowadays be partly removed with the use of specialized mapping functions. Moreover, tropospheric stratification and short wavelength spatial turbulences produce an additive noise to the low amplitude ground deformations calculated by the (multitemporal) InSAR methodology. InSAR atmospheric phase delay corrections are much more challenging, as opposed to GNSS PPP, due to the single pass geometry and the gridded nature of the acquired data. Thus, the precise knowledge of the tropospheric parameters along the propagation medium is extremely useful for the estimation and correction of the atmospheric phase delay. In this context, the PaTrop experiment aims to maximize the potential of using a high-resolution Limited-Area Model for the calculation and removal of the tropospheric noise from InSAR data, by following a synergistic approach and integrating all the latest advances in the fields of remote sensing meteorology (GNSS and InSAR) and Numerical Weather Forecasting (WRF). In the first phase of the experiment, presented in the current paper, we investigate the extent to which a high-resolution 1 km WRF weather re-analysis can produce detailed tropospheric delay maps of the required accuracy, by coupling its output (in terms of Zenith Total Delay or ZTD) with the vertical delay component in GNSS measurements. The model is initially operated with varying parameterization, with GNSS measurements providing a benchmark of real atmospheric conditions. Subsequently, the final WRF daily re-analysis run covers an extended period of one year, based on the optimum model parameterization scheme demonstrated by the parametric analysis. The two datasets (predicted and observed) are compared and statistically evaluated, in order to investigate the extent to which meteorological parameters that affect ZTD can be simulated accurately by the model under different weather conditions. Results demonstrate a strong correlation between predicted and observed ZTDs at the 19 GNSS stations throughout the year (R ranges from 0.91 to 0.93), with an average mean bias (MB) of –19.2 mm, indicating that the model tends to slightly underestimate the tropospheric ZTD as compared to the GNSS derived values. With respect to the seasonal component, model performance is better during the autumn period (October–December), followed by spring (April–June). Setting the acceptable bias range at ±23 mm (equal to the amplitude of one Sentinel-1 C-band phase cycle when projected to the zenithal distance), it is demonstrated that the model produces satisfactory results, with a percentage of ZTD values within the bias margin ranging from 57% in summer to 63% in autumn.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.An amendment to this paper has been published and can be accessed via a link at the top of the paper.

We identify the source of the Mw = 5.6 earthquake that hit west-central Epirus on 21 March 2020 00:49:52 UTC. We use Sentinel-1 synthetic aperture radar interferograms tied to one permanent Global Navigation Satellite System (GNSS) station (GARD). We model the source by inverting the INSAR displacement data. The inversion model suggests a shallow source on a low-angle fault (39°) dipping towards east with a centroid depth of 8.5 km. The seismic moment deduced from our model agrees with those of the published seismic moment tensors. This geometry is compatible with reverse-slip motion along the west-verging Margariti thrust fault that accommodates part of the convergence within the collision zone between Apulia and Eurasia. We also processed new GNSS data and estimate a total convergence rate between Apulia and Eurasia of 8.9 mm yr−1, of which the shortening of the crust between the Epirus coastal GNSS stations and station PAXO in the Ionian Sea (across the Ionian Thrust) is equivalent to ~50% of it or 4.6 mm yr−1. By back-slip modelling we found that a 60-km wide deformation zone takes up nearly most of the convergence between Apulia-Eurasia, trending N318°E. Its central axis runs along the southwest coast of Corfu, along the northeast coast of Paxoi, heading toward the northern extremity of the Lefkada island. The island of Paxoi appears kinematically as part of the Apulian plate.

A strong, shallow earthquake occurred near Heraklion (Crete, Greece) on 27 September 2021. The earthquake produced significant ground deformation in the vicinity of Arkalochori village but without any evidence for surface ruptures of primary origin. We used geodetic (InSAR and GNSS) data to map motions of the Earth’s surface that occurred during and shortly after the earthquake. A 14 cm subsidence of the GNSS station ARKL and a maximum of 19 cm distance from the SAR satellite were recorded. The measured surface displacements were used to constrain the rupture geometry and slip distribution at depth. Our best-fitting inversion model suggests that the rupture occurred on a 13 km-long planar normal fault striking N195° E dipping 55° to the northwest, with major slip occurring to the east and updip of the hypocentre. The fault tip is located 1.2 km beneath the surface. The maximum coseismic slip occurred in the uppermost crust, in the depth interval of 4–6 km. A decrease in the fault offsets toward the Earth’s surface is likely caused by an increased frictional resistance of the shallow layers to rapid coseismic slip. Satellite observations made in the first month after the earthquake detected no post-seismic deformation (i.e., below one fringe or 2.8 cm). The seismic fault may be identified with the Avli (Lagouta) segment of the NNE-SSW striking, west-dipping, 23 km-long neotectonic Kastelli Fault Zone (KFZ). Part of the rupture occurred along the Kastelli segment, indicating a fault segment linkage and a history of overlapping ruptures along KFZ. Based on geological data and footwall topography we estimate an average slip rate between 0.17–0.26 mm/yr for the KFZ. The Arkalochori earthquake is a paradigm example for the on-going extension of Heraklion basin (central Crete) in the WNW-ESE direction, which is almost orthogonal to the E-W Messara graben and other active faults along the south coast of Crete.

GROUND-BASED measurements of volcano deformation can be used to assess eruptive hazard, but require the costly (and often hazardous) installation and maintenance of an instrument network. Here we show that spaceborne radar interferometry, which has already shown its utility in mapping earthquake-related deformation 1 , can be used to monitor long-term volcano deformation. Two families of synthetic aperture radar images, acquired from ascending and descending orbits by the satellite ERS-1, and looking at Mount Etna from opposite sides, cover the time period from 17 May 1992 to 24 October 1993, and include the second half of Etna's most recent eruption. Despite artefacts of the interferometric technique, we can observe a volcano-wide deflation, which is an expected consequence of the eruption, but which had not previously been appreciated. We quantify it using a simple model based on the change of pressure in a sphere located in an elastic half-space; the modelled deformation increases linearly with time until the end of the eruption. Our results show that it will be possible to use this technique to detect the inflation of volcanic edifices that usually precedes eruptions.