Interferometric Synthetic Aperture Radar Velocity Research Articles

Tectonic deformation in El Salvador from combined InSAR and GNSS data

The country of El Salvador, located at the convergent boundary of the Cocos and Caribbean plates, is subject to frequent seismic events and complex surface deformation. This paper presents the first Interferometric Synthetic Aperture Radar (InSAR)-based velocity field of El Salvador produced by combining InSAR and new GNSS data, which allows us to gain a deeper understanding of the intricate crustal deformation processes in the region. We use L-Band ALOS PALSAR images from five ascending and two descending tracks, as well as 171 continuous and episodic GNSS stations across El Salvador and the Caribbean. The inversion of GNSS and InSAR velocities along fault-perpendicular profiles enable us to assess the accommodation of tectonic deformation along the El Salvador Fault Zone (ESFZ), providing a better constraint on the behaviour of its main structures. Additionally, we are able to better understand how the deformation is transmitted at the western (Jalpatagua-Ipala) and eastern (Fonseca Gulf) terminations of the ESFZ.

The 2022 Mw 6.6 Menyuan Earthquake in the Northwest Margin of Tibet: Geodetic and Seismic Evidence of the Fault Structure and Slip Behavior of the Qilian–Haiyuan Strike-Slip Fault

AbstractOn 8 January 2022, an Mw 6.6 earthquake occurred along the Lenglongling fault (LLLF) (Menyuan, Qinghai Province, China), in the northeast margin of the Tibetan plateau. Here, we use imaging geodesy observations from pixel-offset tracking (POT) and Interferometric Synthetic Aperture Radar (InSAR) to quantify the coseismic displacement features and then combine InSAR and teleseismic velocity waveforms to constrain a finite-fault model to examine the kinematic slip behavior of the LLLF. We observe maximum displacements of >50 cm in the line-of-sight InSAR observations and >1 m in the east–west POT observations. Our preferred uniform-slip model indicates that the causative fault of this 2022 event ruptured the middle segment of the Qilian–Haiyuan fault (i.e., the LLLF) with a strike of 108.0°, a dip angle of 83.6°, and slip concentrated within a 15×8 km2 area and reaching a maximum value of 4 m. The total scalar seismic moment was 1.24×1019 N·m, and ∼67% of the energy was released in the first 6 s. We speculate that the 2022 event was controlled by complex fault geometries and slip partitioning, thus raising potential seismic risks associated with the outward extrusion of the plateau.

Ground Displacement in East Azerbaijan Province, Iran, Revealed by L-band and C-band InSAR Analyses.

Iran, as a semi-arid and arid country, has a water challenge in the recent decades and underground water extraction has been increased because of improper developments in the agricultural sector. Thus, detection and measurement of ground subsidence in major plains is of great importance for hazard mitigation purposes. In this study, we carried out a time series small baseline subset (SBAS) interferometric synthetic aperture radar (InSAR) analysis of 15 L-band PALSAR-2 images acquired from ascending orbits of the ALOS-2 satellite between 2015 and 2020 to investigate long-term ground displacements in East Azerbaijan Province, Iran. We found that two major parts of the study area (Tabriz and Shabestar plains) are subsiding, where the mean and maximum vertical subsidence rates are −10 and −98 mm/year, respectively. The results revealed that the visible subsidence patterns in the study area are associated with either anthropogenic activities (e.g., underground water usage) or presence of compressible soils along the Tabriz–Shabestar and Tabriz–Azarshahr railways. This implies that infrastructure such as railways and roads is vulnerable if progressive ground subsidence takes over the whole area. The SBAS results deduced from L-band PALSAR-2 data were validated with field observations and compared with C-band Sentinel-1 results for the same period. The C-band Sentinel-1 results showed good agreement with the L-band PALSAR-2 dataset, in which the mean and maximum vertical subsidence rates are −13 and −120 mm/year, respectively. For better visualization of the results, the SBAS InSAR velocity map was down-sampled and principal component analysis (PCA) was performed on ~3600 randomly selected time series of the study area, and the results are presented by two principal components (PC1 and PC2).

Creep Along the Central San Andreas Fault From Surface Fractures, Topographic Differencing, and InSAR

AbstractImaging tectonic creep along active faults is critical for measuring strain accumulation and ultimately understanding the physical processes that guide creep and the potential for seismicity. We image tectonic deformation along the central creeping section of the San Andreas Fault at the Dry Lake Valley paleoseismic site (36.468°N, 121.055°W) using three data sets with varying spatial and temporal scales: (1) an Interferometric Synthetic Aperture Radar (InSAR) velocity field with an ~100‐km footprint produced from Sentinel‐1 satellite imagery, (2) light detection and ranging (lidar) and structure‐from‐motion 3‐D topographic differencing that resolves a decade of deformation over a 1‐km aperture, and (3) surface fractures that formed over the 3‐ to 4‐m wide fault zone during a drought from late 2012 to 2014. The InSAR velocity map shows that shallow deformation is localized to the San Andreas Fault. We demonstrate a novel approach for differencing airborne lidar and structure‐from‐motion topography that facilitates resolving deformation along and adjacent to the San Andreas Fault. The 40‐m resolution topographic differencing resolves a 2.5 ± 0.2 cm/yr slip rate localized to the fault. The opening‐mode fractures accommodate cm/yr of fault slip. A 90% ± 30% of the 1‐km aperture deformation is accommodated over the several meter‐wide surface trace of the San Andreas Fault. The extension direction inferred from the opening‐mode fractures and topographic differencing is 40°–48° from the local trend of the San Andreas Fault. The localization of deformation likely reflects the well‐oriented and mature fault.

Present-day tectonic activity along the central section of the Altyn Tagh fault derived from time series InSAR

The Altyn Tagh Fault (ATF) is a large-scale complex tectonic system. In this study, the present-day crustal deformation of the central section of the ATF (90.8°E−91.58°E) was obtained using 14 images on a descending track acquired between 2007 and 2010 from Advanced Land Observing Satellite (ALOS). To improve the accuracy of the interferograms, ALOS World 3D–30m (AW3D-30)Digital Surface Model (DSM) from the Japan Aerospace Exploration Agency was used in Small Baseline Subset (SBAS) Interferometric Synthetic Aperture Radar (InSAR) processing. The Line of Sight (LOS) deformation map show that there is an obvious zoning feature. With the main ATF as the boundary, the north block is concentrated around −35∼ −60 mm, and the south block is concentrated around −9 ∼ 11 mm. Based on the InSAR velocity map, we inverted for the strike-slip rate and locking depth of the fault using the two-dimensional strike-slip fault buried-dislocation model. The inversion results for the strike-slip rate at three selected cross-section locations perpendicular to the ATF were 6.1 mm/a, 5.3 mm/a and 7.9 mm/a from west to east; the corresponding locking depths were 9.5 km, 6.8 km and 12.3 km from west to east. The location and trend of the fault obtained by inversion are coincides with the Xorkol seismic belt.

Ground Deformation Analysis of Bolvadin (W. Turkey) by Means of Multi-Temporal InSAR Techniques and Sentinel-1 Data

Surface deformations were observed in Bolvadin, located in the province of Afyon (western Turkey) in the last decade which occurred without any destructive earthquakes. In this study, ground deformation of the Bolvadin region is analyzed by means of multi-temporal interferometric synthetic aperture radar (InSAR) technique with Sentinel-1 synthetic aperture radar (SAR) data. Sentinel-1 data acquired in ascending and descending orbits between October 2014 and October 2018 are processed with the Sentinel Application Platform (SNAP) and Stanford Method for Persistent Scatterers (StaMPS) open source software tools. Deformation velocity maps and line-of-sight (LOS) displacement time series are produced and compared with geology, groundwater level and the water surface area of Eber Lake nearby. Deformation velocity maps reveal significant subsidence in most of the town and surrounding regions, which is confirmed by field observations that show severe damage to the settlements and infrastructure of the town. The most severe subsidence is observed to be in the southern part of Bolvadin with rates up to 35 mm/year, which is characterized by the presence of soft alluvial deposits. Composed of slope debris/talus and conglomerate, the northeastern part of the deforming region experiences a relatively lower rate of subsidence. A strong correlation between LOS displacement time series and groundwater level exists both in the short and long term. Moreover, short term variations in LOS displacement time series also show high similarity with seasonal variations in the water surface area of Eber Lake located a few km southeast of the town. We conclude that the primary cause of subsidence is most probably the overexploitation of groundwater and hydrological changes because of (1) the strong correlation of subsidence with lithological units, (2) the similarity between deformation rate and groundwater level changes, (3) the correspondence of seasonal variations in water surface area and short-term deformation rate oscillations, and (4) the absence of InSAR velocity contrast across the active faults.

Strain Partitioning and Present‐Day Fault Kinematics in NW Tibet From Envisat SAR Interferometry

AbstractAn 8 year archive of Envisat synthetic aperture radar (SAR) data over a 300 × 500 km2 wide area in northwestern Tibet is analyzed to construct a line‐of‐sight map of the current surface velocity field. The resulting velocity map reveals (1) a velocity gradient across the Altyn Tagh fault, (2) a sharp velocity change along a structure following the base of the alluvial fans in southern Tarim, and (3) a broad velocity gradient, following the Jinsha suture. The interferometric synthetic aperture radar velocity field is combined with published GPS data to constrain the geometry and slip rates of a fault model consisting of a vertical fault plane under the Altyn Tagh fault and a shallow flat décollement ending in a steeper ramp on the Tarim side. The solutions converge toward 0.7 mm/yr of pure thrusting on the décollement‐ramp system and 10.5 mm/yr of left‐lateral strike‐slip movement on the Altyn Tagh fault, below a 17 km locking depth. A simple elastic dislocation model across the Jinsha suture shows that data are consistent with 4–8 mm/yr of left‐lateral shear across this structure. Interferometric synthetic aperture radar processing steps include implementing a stepwise unwrapping method starting with high‐quality interferograms to assist in unwrapping noisier interferograms, iteratively estimating long‐wavelength spatial ramps, and referencing all interferograms to bedrock pixels surrounding sedimentary basins. A specific focus on atmospheric delay estimation using the ERA‐Interim model decreases the uncertainty on the velocity across the Tibet border by a factor of 2.

Assessing the Feasibility of a National InSAR Ground Deformation Map of Great Britain with Sentinel-1

This work assesses the feasibility of national ground deformation monitoring of Great Britain using synthetic aperture radar (SAR) imagery acquired by Copernicus’ Sentinel-1 constellation and interferometric SAR (InSAR) analyses. As of December 2016, the assessment reveals that, since May 2015, more than 250 interferometric wide (IW) swath products have been acquired on average every month by the constellation at regular revisit cycles for the entirety of Great Britain. A simulation of radar distortions (layover, foreshortening, and shadow) confirms that topographic constraints have a limited effect on SAR visibility of the landmass and, despite the predominance of rural land cover types, there is potential for over 22,000,000 intermittent small baseline subset (ISBAS) monitoring targets for each acquisition geometry (ascending and descending) using a set of IW image frames covering the entire landmass. Finally, InSAR results derived through ISBAS processing of the Doncaster area with an increasing amount of Sentinel-1 IW scenes reveal a consistent decrease of standard deviation of InSAR velocities from 6 mm/year to ≤2 mm/year. Such results can be integrated with geological and geohazard susceptibility data and provide key information to inform the government, other institutions and the public on the stability of the landmass.

Interseismic strain accumulation across the central North Anatolian Fault from iteratively unwrapped InSAR measurements

AbstractThe North Anatolian Fault (NAF) is a major tectonic feature in the Middle East and is the most active fault in Turkey. The central portion of the NAF is a region of Global Navigation Satellite Systems (GNSS) scarcity. Previous studies of interseismic deformation have focused on the aseismic creep near the town of Ismetpasa using radar data acquired in a single line‐of‐sight direction, requiring several modeling assumptions. We have measured interseismic deformation across the NAF using both ascending and descending data from the Envisat satellite mission acquired between 2003 and 2010. Rather than rejecting incorrectly unwrapped areas in the interferograms, we develop a new iterative unwrapping procedure for small baseline interferometric synthetic aperture radar (InSAR) processing that expands the spatial coverage. Our method corrects unwrapping errors iteratively and increases the robustness of the unwrapping procedure. We remove long wavelength trends from the InSAR data using GNSS observations and deconvolve the InSAR velocities into fault‐parallel motion. Profiles of fault‐parallel velocity reveal a systematic eastward decrease in fault slip rate from 30 mm/yr (25–34, 95% confidence interval (CI)) to 21 mm/yr (14–27, 95% CI) over a distance of ∼200 km. Direct offset measurements across the fault reveal fault creep along a ∼130 km section of the central NAF, with an average creep rate of 8 ± 2 mm/yr and a maximum creep rate of 14 ± 2 mm/yr located ∼30 km east of Ismetpasa. As fault creep is releasing only 30–40% of the long‐term strain in the shallow crust, the fault is still capable of producing large, damaging earthquakes in this region.

Interseismic Deformation of the Altyn Tagh Fault Determined by Interferometric Synthetic Aperture Radar (InSAR) Measurements

The Altyn Tagh Fault (ATF) is one of the major left-lateral strike-slip faults in the northeastern area of the Tibetan Plateau. In this study, the interseismic deformation across the ATF at 85°E was measured using 216 interferograms from 33 ENVISAT advanced synthetic aperture radar images on a descending track acquired from 2003 to 2010, and 66 interferograms from 15 advanced synthetic aperture radar images on an ascending track acquired from 2005 to 2010. To retrieve the pattern of interseismic strain accumulation, a global atmospheric model (ERA-Interim) provided by the European Center for Medium Range Weather Forecast and a global network orbital correction approach were applied to remove atmospheric effects and the long-wavelength orbital errors in the interferograms. Then, the interferometric synthetic aperture radar (InSAR) time series with atmospheric estimation model was used to obtain a deformation rate map for the ATF. Based on the InSAR velocity map, the regional strain rates field was calculated for the first time using the multi-scale wavelet method. The strain accumulation is strongly focused on the ATF with the maximum strain rate of 12.4 × 10−8/year. We also show that high-resolution 2-D strain rates field can be calculated from InSAR alone, even without GPS data. Using a simple half-space elastic screw dislocation model, the slip-rate and locking depth were estimated with both ascending and descending surface velocity measurements. The joint inversion results are consistent with a left-lateral slip rate of 8.0 ± 0.7 mm/year on the ATF and a locking depth of 14.5 ± 3 km, which is in agreement with previous results from GPS surveys and ERS InSAR results. Our results support the dynamic models of Asian deformation requiring low fault slip rate.

InSAR bias and uncertainty due to the systematic and stochastic tropospheric delay

AbstractWe quantify the bias and uncertainty of interferometric synthetic aperture radar (InSAR) displacement time series and their derivatives, the displacement velocities, by analyzing the systematic and stochastic components of the temporal variation of the tropospheric delay. The biases due to the systematic seasonal delay depend on the SAR acquisition times, whereas the uncertainties depend on the standard deviation of the random delay, the number of acquisitions, the total time span covered, and the covariance of the time series of the stochastic delay between a pixel and the reference. We study the contribution of the wet delay to the InSAR observations along the western India plate boundary using (i) Moderate Resolution Imaging Spectroradiometer precipitable water vapor, (ii) stratified tropospheric delay estimated from the ERA‐I global atmospheric model, and (iii) seven Envisat InSAR swaths. Our analysis indicates that the amplitudes of the annual delay vary by up to ~10 cm in this region equivalent to a maximum displacement bias of ~24 cm in InSAR line of sight direction between two epochs (assuming Envisat IS6 beam mode). The stratified tropospheric delay correction mitigates this bias and reduces the scatter due to the stochastic delay. For ~7 years of Envisat acquisitions along the western India plate boundary, the uncertainty of the InSAR velocity field due to the residual stochastic wet delay after stratified tropospheric delay correction using the ERA‐I model is in the order of ~2 mm/yr over 100 km and ~4 mm/yr over 400 km. We discuss the implication of the derived uncertainties on the full variance‐covariance matrix of the InSAR data.

From quiescence to unrest: 20 years of satellite geodetic measurements at Santorini volcano, Greece

AbstractPeriods of unrest at caldera‐forming volcanic systems characterized by increased rates of seismicity and deformation are well documented. Some can be linked to eventual eruptive activity, while others are followed by a return to quiescence. Here we use a 20 year record of interferometric synthetic aperture radar (InSAR) and GPS measurements from Santorini volcano to further our understanding of geodetic signals at a caldera‐forming volcano during the periods of both quiescence and unrest, with measurements spanning a phase of quiescence and slow subsidence (1993–2010), followed by a phase of unrest (January 2011 to April 2012) with caldera‐wide inflation and seismicity. Mean InSAR velocity maps from 1993–2010 indicate an average subsidence rate of ~6 mm/yr over the southern half of the intracaldera island Nea Kameni. This subsidence can be accounted for by a combination of thermal contraction of the 1866–1870 lava flows and load‐induced relaxation of the substrate. For the period of unrest, we use a joint inversion technique to convert InSAR measurements from three separate satellite tracks and GPS observations from 10 continuous sites into a time series of subsurface volume change. The optimal location of the inflating source is consistent with previous studies, situated north of Nea Kameni at a depth of ~4 km. However, the time series reveals two distinct pressure pulses. The first pulse corresponds to a volume change (ΔV) within the shallow magma chamber of (11.56 ± 0.14) × 106 m3, and the second pulse has a ΔV of (9.73 ± 0.10) × 106 m3. The relationship between the timing of these pulses and microseismicity observations suggests that these pulses may be driven by two separate batches of magma supplied to a shallow reservoir. We find no evidence suggesting a change in source location between the two pulses. The decline in the rates of volume change at the end of both pulses and the observed lag of the deformation signal behind cumulative seismicity, suggest a viscoelastic response. We use a simple model to show that two separate pulses of magma intruding into a shallow magma chamber surrounded by a viscoelastic shell can account for the observed temporal variation in cumulative volume change and seismicity throughout the period of unrest. Given the similarities between the geodetic signals observed here and at other systems, this viscoelastic model has potential use for understanding behavior at other caldera systems.

Iceland rising: Solid Earth response to ice retreat inferred from satellite radar interferometry and visocelastic modeling

A broad uplift occurs in Iceland in response to the retreat of ice caps, which began circa 1890. Until now, this deformation signal has been measured primarily using GPS at points some distance away from the ice caps. Here, for the first time we use satellite radar interferometry (interferometric synthetic aperture radar) to constrain uplift of the ground all the way up to the edge of the largest ice cap, Vatnajökull. This allows for improved constraints on the Earth rheology, both the thickness of the uppermost Earth layer that responds only in an elastic manner and the viscosity below it. The interferometric synthetic aperture radar velocities indicate a maximum displacement rate of 24±4 and 31±4 mm/yr at the edge of Vatnajökull, during 1995–2002 and 2004–2009, respectively. The fastest rates occur at outlet glaciers of low elevation where ice retreat is high. We compare the observations with glacial isostatic adjustment models that include the deglaciation history of the Icelandic ice caps since 1890 and two Earth layers. Using a Bayesian approach, we derived probability density functions for the average Earth model parameters for three satellite tracks. Based on our assumptions, the three best fit models give elastic thicknesses in the range of 15–40 km, and viscosities ranging from 4–10× 1018 Pa s.

Interferometric synthetic aperture radar–GPS integration: Interseismic strain accumulation across the Hunter Mountain fault in the eastern California shear zone

The principal limitations of interferometric synthetic aperture radar (InSAR) to measure subtle, long‐wavelength deformation are uncertainties associated with the satellite orbits. We propose a method to remove orbital phase errors from the InSAR data by integrating InSAR and continuous GPS time series. We model the along‐track variation of the baseline errors as second‐order polynomials and estimate the coefficients using the continuous GPS measurements. We apply this method to a 600 km long region encompassing the Basin and Range and the eastern California shear zone. Comparison of the corrected InSAR velocities with independent GPS data shows that this method removes the long‐wavelength InSAR errors. The InSAR data reveal a region of sharp variation in the line‐of‐sight velocity across the Hunter Mountain fault. We model the deformation as interseismic elastic strain accumulation across a strike‐slip fault. The modeling suggests a fault slip rate of 4.9 ± 0.8 mm/yr and a locking depth of 2 ± 0.4 km. The shallow locking depth suggests that the Hunter Mountain fault is a transfer fault between low angle normal faults in the area.

Ice velocity measurements of Langjökull, Iceland, from interferometric synthetic aperture radar (InSAR)

AbstractWe have computed the velocity of Langjökull ice cap (LIC), Iceland, using InSAR data acquired by the European Remote-sensing Satellite in 1994 and an ice-surface elevation model derived from a GPS survey conducted in 1997. The velocity map covers 68% of the 925 km2 of LIC. Eight principal outlet glaciers were identified, along which ice velocities increased from zero at the ice divide to a maximum of 75 m a−1. When the InSAR velocity dataset is compared with surface GPS measurements of ice flow determined in summer 2001, there is generally good agreement, suggesting interannual and seasonal stability for much of the ice cap during this period. However, non-zero ice velocities at the termini of four of the eight glaciers surveyed suggest that these glaciers were advancing in 1994. Comparison of our GPS and InSAR velocities suggests that Sudurjökull experienced elevated velocity in late February 1994, which, when other studies are considered, appears to be consistent with a glacier surge.

Marine ice beneath the Filchner–Ronne Ice Shelf, Antarctica: a comparison of estimated thickness distributions

AbstractIn an earlier study, melt/freeze rates beneath most of the Filchner–Ronne Ice Shelf, Antarctica, were estimated using an assumption of a steady-state ice shelf applied to a velocity field derived from RADARSAT and ERS-1 interferometric synthetic aperture radar (InSAR) data and an ice-thickness map inferred from ERS-1 satellite radar altimeter data. Here, we use these basal accumulation rates and the InSAR velocity data to estimate the distribution and thickness of marine ice beneath the Filchner–Ronne Ice Shelf. These estimates are compared with a marine-ice thickness map derived from radio-echo sounding (RES) data and airborne radar altimetry. In general, we find close agreement between these estimates. In the few locations where there are significant differences, the discrepancies are largely attributable to artifacts in the radar-altimeter thickness map. With improvements such as a much smaller footprint, the next generation of altimeters should overcome many of these limitations, leading to improved marine-ice accumulation estimates. Overall, the good agreement with the RES data validates the InSAR-based estimates over much of the Ronne Ice Shelf.

Post-stagnation behavior in the upstream regions of Ice Stream C, West Antarctica

AbstractThe region where two active tributaries feed into the now stagnant Ice Stream C (ISC), West Antarctica, is thickening. In this region, we observe a correlation between faster ice flow (the tributaries) and elevated topography. We conclude that stagnation of ISC resulted in compression and thickening along the tributaries, eventually forming a “bulge” on the ice-sheet surface. Modern hydraulic potential gradients would divert basal meltwater from ISC to Ice Stream B (ISB). These gradients are primarily controlled by the bulge topography, and so likely formed subsequent to trunk stagnation. As such, we argue against “water piracy” as being the cause for ISC’s stagnation. Kinematic-wave theory suggests that thickness perturbations propagate downstream over time, but that kinematic-wave speed decreases near the stagnant trunk. This and modest diffusion rates combine to trap most of the tributary-fed ice in the bulge region. Using interfero-metric synthetic aperture radar velocity measurements, we observe that half of the ice within ISC’s southern tributary flows into ISB. That flow pattern and other observations of non-steady flow in the region likely result from stagnation-induced thickening along upper ISC combined with a longer period of thinning on upper ISB. If current trends in thickness change continue, more ice from upper ISC will be diverted to ISB.