Satellite-based Interferometric Synthetic Aperture Radar Research Articles

Two-dimensional thermal and dynamical strain in landfast sea ice from InSAR: results from a new analytical inverse method and field observations

Abstract Observing continuous strain in sea ice at geophysical scales of tens of meters to kilometers requires displacement measurements made with millimeter-scale precision. Satellite-based interferometric synthetic aperture radar (InSAR) provides such precise measurements of relative surface displacement over broad spatial areas at regular intervals and, unlike point displacement measurements, it allows confident delineation of continuously deforming regions. However, InSAR only captures the 1-D component of surface displacement parallel to a radar's lines-of-sight. Additional analysis is required to translate between these 1-D observations and the horizontal or vertical displacements they arose from. Previous studies utilize an iterative inverse model to constrain estimates of horizontal surface displacement from InSAR. Here we build upon that work outlining a new analytical inverse modeling method for quantifying displacement and strain over continuous regions of sea ice and provide comparison between model results and independent displacement observations. We demonstrate the inverse method over both landfast and drifting ice along the Alaskan coastline. These intercomparisons highlight environments in which displacements inverted from interferograms may be used as an independent estimator of surface strain, as well as the potential for the outlined inverse methods to be used in conjunction with other observing methods.

Application of Sentinel-1 InSAR to monitor tailings dams and predict geotechnical instability: practical considerations based on case study insights

Tailings storage facilities (TSFs) impound mining waste behind dams to ensure public safety, but failure incidents have prompted calls for more robust monitoring programs. Satellite-based interferometric synthetic aperture radar (InSAR) has grown in popularity due to its ability to remotely detect millimeter-scale displacements in most urban and some natural terrains. However, there remains a limited understanding of whether InSAR can be as accurate or representative as on-the-ground instruments, whether failures can be predicted in advance using InSAR, and what variables govern the quality and reliability of InSAR results. To address these gaps, we analyze open-source, medium-resolution Sentinel-1 data to undertake a ground-truth assessment at a test site and a forensic analysis of five failure cases. We use a commercial software with an automated Persistent Scatterer (PS) workflow (SARScape Analytics) for all case study sites except one and a proprietary algorithm (SqueeSAR) with a dual PS and Distributed Scatterer (DS) algorithm for the ground-truth site and one forensic case. The main goal is to deliver practical insights regarding the influence of algorithm/satellite selection, environmental conditions, site activity, coherence thresholds, satellite-dam geometry, and failure modes. We conclude that Sentinel-1 InSAR can serve as a hazard-screening tool to help guide where to undertake targeted investigations; however, most potential failure modes may not exhibit InSAR-detectable accelerations that could assist with time-of-failure prediction in real time. As such, long-term monitoring programs should ideally be integrated with a combination of remote sensing and field instrumentation to best support engineering practice and judgment.

Interpretation of Bridge Health Monitoring Data from Satellite InSAR Technology

This paper presents a study on applying satellite Interferometric Synthetic Aperture Radar (InSAR) technology for the remote monitoring of road bridges and interpreting the results from a structural standpoint. The motivation behind this study arises from the widespread deterioration observed in many road bridges worldwide, leading to the need for large-scale, economic, and effective structural health monitoring (SHM) techniques. While traditional contact-type sensors have cost sustainability limitations, remote sensing techniques, including satellite-based InSAR, offer interesting alternative solutions. The objective of this study is three-fold: (i) to process InSAR data specifically for road bridges in operational conditions through the Multi-Temporal InSAR technique and extract displacement time series of reflective targets on their decks; (ii) to interpret the observed millimetric bridge displacements to verify the consistency with expected response to environmental loads and the possibility to detect unexpected behaviours; and (iii) to investigate the correlation between bridge displacements and environmental loads as temperature and river water flow variations. The study focuses on the multi-span prestressed concrete A22 Po River Bridge in Italy, utilising a dataset of X-Band HIMAGE mode Stripmap images acquired over eight years by the satellite constellation COSMO-SkyMed. The study demonstrates the effectiveness of InSAR-based SHM in detecting temperature-induced displacements and identifying different bridge spans simply by studying the sign of the correlation between displacements and temperature variation. It also reveals an unexpected behaviour in a portion of the bridge retrofitted to prevent scour issues a few years before the dataset start date. Furthermore, the correlation between pier displacements and river level variations underscores the importance of considering environmental factors and the geotechnical characteristics of the foundation soils in bridge monitoring. The results obtained from this study are significant with a view to using this satellite InSAR-based monitoring for early detection of anomalous bridge behaviour on a large scale.

Activity of Okgye Limestone Mine in South Korea Observed by InSAR Coherence and PSInSAR Techniques

The Okgye limestone mine, which is the largest open-pit limestone mine located in a mountainous area in Korea, suffered a collapse in 2012 that claimed four casualties. Restoration work on the rocky mined-out slopes, as well as mining and dumping activities, are still in progress. Monitoring slope stability is important to prevent the sudden collapse of slopes, which can be efficiently performed by satellite-based interferometric synthetic aperture radar (InSAR) techniques. Firstly, we obtained elevation changes using InSAR-generated Copernicus 30 m DEM in 2014 and an SRTM 1Sec DEM in 2000, through which the area was roughly classified into the mining area, tailings storage area, and the mined-out area. A time series of 12-day coherence images produced by Sentinel-1B SAR were averaged annually to produce an RGB-composite image to observe the change in mining activities during 2018, 2019, and 2020. We found many persistent scatterers (PS) when observing the ground displacement, both in the ascending and descending orbits, from which we decomposed this into the vertical and east components. The largest displacement of 63.6 mm/year was observed during 2019 and 2020 in the tailings storage area in the direction of the dumping slope. For the rocky outcrops and the transmission tower, we found a seasonal oscillation, which can be interpreted as the thermal expansion of limestone and iron. This paper demonstrated that the surface stability and deformation of open-pit mines could be effectively monitored by combining InSAR DEM, coherence, and PSInSAR techniques.

Surface Deformation Surrounding the 2021 Mw 7.2 Haiti Earthquake Illuminated by InSAR Observations

ABSTRACT Earthquakes pose a major threat to the people of Haiti, as tragically shown by the catastrophic 2010 Mw 7.0 earthquake and more recently by the 2021 Mw 7.2 earthquake. Both events occurred within the transpressional Enriquillo–Plantain Garden fault zone (EPGFZ), which runs through the southern peninsula of Haiti and is a major source of seismic hazard for the region. Satellite-based Interferometric Synthetic Aperture Radar (InSAR) data are used to illuminate the ground deformation patterns associated with the 2021 event. The analysis of Sentinel-1 and Advanced Land Observation Satellite (ALOS)-2 InSAR data shows (1) the broad coseismic deformation field; (2) detailed secondary fault structures as far as 12 km from the main Enriquillo–Plantain Garden fault (EPGF), which are active during and after the earthquake; and (3) postseismic shallow slip, which migrates along an ∼40 km unruptured section of the EPGF for approximately two weeks following the earthquake. The involvement of secondary faults in this rupture requires adjustments to the representation of hazard that assumes a simple segmented strike-slip EPGF. This work presents the first successful use of phase gradient techniques to map postseismic deformation in a vegetated region, which opens the door to future studies of a larger number of events in a wider variety of climates.

A country-wide assessment of Iran's land subsidence susceptibility using satellite-based InSAR and machine learning

Land subsidence (LS), which mainly results from poor watershed management, is a complex and nonlinear phenomenon. In the present study, LS at a country-wide assessment of Iran was mapped by using several geo-environmental conditioning factors (namely, altitude, slope degree and aspect, plan and profile curvature, distance from a river, road or fault, rainfall, geology and land use) into a machine learning algorithm-based artificial neural network (ANN), and a powerful group method of data handling (GMDH). The total dataset includes historical LS and non-LS locations, identified by the interferometric synthetic aperture radar (InSAR). The whole dataset was divided into two subsets at a ratio of 70:30 for training and validating the model, respectively. ANN- and GMDH-based LS maps were evaluated using receiver-operator characteristic (ROC) curves. The information gain ratio (IGR) was calculated to determine the relative importance of the conditioning factors. The results showed that all of the considered factors contributed significantly to the LS mapping in Iran, with geology having the strongest impact. According to the ROC curve analysis, both ANN and GMDH-based LS maps were accurate, but the map obtained by the GMDH approach had a higher accuracy than that of ANN. Southwestern, northeastern and some parts of the central region of Iran were shown to be susceptible to LS in the future. According to the GMDH susceptibility map, 10% of Iran exhibits high or very high susceptibility to LS in the future. The provinces of Hamedan and Khouzestan had the highest percentage of areas at risk of LS. According to the InSAR deformation map, 39%, 20%, 25%, 13% and 3% of the investigated areas are subject to a yearly LS of −1 to −2.5, −2.5 to −5, −5 to −7.5, −7.5 to −10 and −10 to −20 cm, respectively. The province of Razavi Khorasan in the northeast of Iran had the largest area (about 3500 km2) vulnerable to LS occurrence. Based on the LS susceptibility map, the provinces of Ardebil, Kurdistan, West and East Azerbaijan, Sistan and Baluchistan and Kermanshah, although not currently undergoing a high rate of LS, will be at high risk of severe LS in the future.

Spatio-Temporal Quality Indicators for Differential Interferometric Synthetic Aperture Radar Data

Satellite-based interferometric synthetic aperture radar (InSAR) is an invaluable technique in the detection and monitoring of changes on the surface of the earth. Its high spatial coverage, weather friendly and remote nature are among the advantages of the tool. The multi-temporal differential InSAR (DInSAR) methods in particular estimate the spatio-temporal evolution of deformation by incorporating information from multiple SAR images. Moreover, opportunities from the DInSAR techniques are accompanied by challenges that affect the final outputs. Resolving the inherent ambiguities of interferometric phases, especially in areas with a high spatio-temporal deformation gradient, represents the main challenge. This brings the necessity of quality indices as important DInSAR data processing tools in achieving ultimate processing outcomes. Often such indices are not provided with the deformation products. In this work, we propose four scores associated with (i) measurement points, (ii) dates of time series, (iii) interferograms and (iv) images involved in the processing. These scores are derived from a redundant set of interferograms and are calculated based on the consistency of the unwrapped interferometric phases in the frame of a least-squares adjustment. The scores reflect the occurrence of phase unwrapping errors and represent valuable input for the analysis and exploitation of the DInSAR results. The proposed tools were tested on 432,311 points, 1795 interferograms and 263 Sentinel-1 single look complex images by employing the small baseline technique in the PSI processing chain, PSIG of the geomatics division of the Centre Tecnològic de Telecomunicacions de Catalunya (CTTC). The results illustrate the importance of the scores—mainly in the interpretation of the DInSAR outputs.

Permafrost in monitored unstable rock slopes in Norway – new insights from temperature and surface velocity measurements, geophysical surveying, and ground temperature modelling

Abstract. The warming and subsequent degradation of mountain permafrost within alpine areas represent an important process influencing the stability of steep slopes and rock faces. The unstable and monitored slopes of Mannen (Møre and Romsdal county, southern Norway) and Gámanjunni-3 (Troms and Finnmark county, northern Norway) were classified as high-risk sites by the Norwegian Geological Survey (NGU). Failure initiation has been suggested to be linked to permafrost degradation, but the detailed permafrost distribution at the sites is unknown. Rock wall (RW) temperature loggers at both sites have measured the thermal regime since 2015, showing mean rock surface temperatures between 2.5 and −1.6 ∘C depending on site and topographic aspect. Between 2016 and 2019 we conducted 2D and 3D electrical resistivity tomography (ERT) surveys on the plateau and directly within the rock wall back scarp of the unstable slopes at both sites. In combination with geophysical laboratory analysis of rock wall samples from both sites, the ERT soundings indicate widespread permafrost areas, especially at Gámanjunni-3. Finally, we conducted 2D thermal modelling to evaluate the potential thermal regime, along with an analysis of available displacement rate measurements based on Global Navigation Satellite System (GNSS) and ground- and satellite-based interferometric synthetic aperture radar (InSAR) methods. Surface air and ground temperatures have increased significantly since ca. 1900 by 1 and 1.5 ∘C, and the highest temperatures have been measured and modelled since 2000 at both study sites. We observed a seasonality of displacement, with increasing velocities during late winter and early spring and the highest velocities in June, probably related to water pressure variations during snowmelt. The displacement rates of Gámanjunni-3 rockslide co-vary with subsurface resistivity and modelled ground temperature. Increased displacement rates seem to be associated with sub-zero ground temperatures and higher ground resistivity. This might be related to the presence of ground ice in fractures and pores close to the melting point, facilitating increased deformation. The study demonstrates and discusses the possible influence of permafrost, at least locally, on the dynamics of large rock slope instabilities.

An Integrated InSAR-Borehole Inclinometer-Numerical Modeling Approach to the Assessment of a Slow-Moving Landslide

ABSTRACT We use results of satellite-based interferometric synthetic aperture radar, Global Positioning System, and borehole inclinometer data to constrain numerical models that improve understanding of slope deformation at the Alexandria landslide, British Columbia, Canada. Surface monitoring data and borehole slope inclinometer measurements provide important insight into the slope failure mechanism. We initially analyzed the data in a geographic information system (GIS) to create thematic maps of the landslide area (hillshade, slope, and aspect), to identify key geological features, and to produce an engineering geomorphology map of the landslide. The monitoring data and the geological/engineering geomorphological models provide important constraints for two-dimensional landslide limit equilibrium and finite difference analyses. The initial limit equilibrium analysis improved understanding of the sliding surfaces. The finite difference models were then used to simulate and investigate the potential slope deformation mechanism. The combined slope monitoring/modeling results show that the Alexandria landslide is a slow-moving, multiple-block, retrogressive slope failure. The close agreement between the limit equilibrium and finite difference analyses, together with the satellite and ground-based slope monitoring and GIS data, highlight the importance of using a multidisciplinary/integrated approach in landslide studies.

Advanced analysis of satellite data reveals ground deformation precursors to the Brumadinho Tailings Dam collapse

Catastrophic failure of a tailings dam at an iron ore mine complex in Brumadinho, Brazil, on 25th January 2019 released 11.7 million m3 of tailings downstream. Although reportedly monitored using an array of geotechnical techniques, the collapse occurred without any apparent warning. It claimed more than 200 lives and caused considerable environmental damage. Here we present the Intermittent Small Baseline Subset (ISBAS) technique on satellite-based interferometric synthetic aperture radar (InSAR) data to assess the course of events. We find that parts of the dam wall and tailings were experiencing deformation not consistent with consolidation settlement preceding the collapse. Furthermore, we show that the timing of the dam collapse would have been foreseeable based on this observed precursory deformation. We conclude that satellite-based monitoring techniques may help mitigate similar catastrophes in the future.

On the use of adjoints in the inversion of observed quasi-static deformation

SUMMARY An adjoint-based conjugate gradient algorithm provides an efficient means for imaging sources of deformation within the Earth, such as volume stresses associated with fluid flow in aquifers and reservoirs. For time intervals over which the overburden deforms elastically, one can calculate the gradient elements for a single model update using just two numerical simulations. The first is a forward run that is used to compute the residuals associated with the given iteration. The second simulation is to evaluate the application of the adjoint operator to the residuals. In this adjoint calculation, the residual displacements are applied as sources at the measurement locations, driving the deformation in the simulation. The volume stress on the source grid blocks, in response to the residual displacements, provide the gradient components. We apply this technique to satellite-based interferometric synthetic aperture radar (InSAR) line-of-sight displacements that were observed over an oil reservoir in California’s Central Valley. We find that the adjoint-based gradient estimates, requiring 18 CPU seconds, agree with conventional numerical calculations that take over 3700 CPU seconds to compute. Conjugate gradient algorithms utilizing the conventional approach and adjoint-based gradient computations give roughly the same reductions in misfit and similar final estimates of reservoir volume change.

Very Local Subsidence Near the Hot Spring Region in Hakone Volcano, Japan, Inferred from InSAR Time Series Analysis of ALOS/PALSAR Data

Monitoring of surface displacement by satellite-based interferometric synthetic aperture radar (InSAR) analysis is an effective method for detecting land subsidence in areas where routes of leveling measurements are undeveloped, such as mountainous areas. In particular, InSAR-based monitoring around well-developed hot spring resorts, such as those in Japan, is useful for conserving hot spring resources. Hakone Volcano is one of the major hot spring resorts in Japan, and many hot spring wells have been developed in the Owakudani fumarole area, where a small phreatic eruption occurred in 2015. In this study, we performed an InSAR time series analysis using the small baseline subset (SBAS) method and ALOS/PALSAR scenes of the Hakone Volcano to monitor surface displacements around the volcano. The results of the SBAS-InSAR time series analysis show highly localized subsidence to the west of Owakudani from 2006–2011 when the ALOS/PALSAR satellite was operated. The area of subsidence was approximately 500 m in diameter, and the peak rate of subsidence was approximately 25 mm/year. Modeling using a point pressure source suggested that the subsidence was caused by a contraction at approximately 700 m above sea level (about 300 m below the ground surface). The rate of this contraction was estimated to be 1.04 × 104 m3/year. Hot spring water is collected from a nearby well at almost the same depth as the contraction source, and its main dissolved ion component is chloride ions, suggesting that the hydrothermal fluids are supplied from deep within the volcano. The land subsidence suggests that the fumarole activity is attenuating due to a decrease in the supply of hydrothermal fluids from deeper areas.

General Survey of Large-scale Land Subsidence by GACOS-Corrected InSAR Stacking: Case Study in North China Plain

Abstract. Satellite-based InSAR (Interferometric Synthetic Aperture Radar) provides an effective way to measure large-scale land surface motions. Currently, the atmospheric phase delay is one of the most critical issues in InSAR deformation monitoring. Generic Atmospheric Correction Online Service (GACOS) is a free, globally available and easy-to-implement tool to generate high-resolution zenith total delay maps, which could be used for InSAR atmospheric delay correction. The mean velocity could then be estimated by stacking multiple GACOS-corrected interferograms. We applied the proposed GACOS-corrected InSAR stacking method in the North China Plain and analysed its performance. Within the 549 interferograms, more than 85 % gained positive correction performances. The correlation between the phase-dZTD indicator and the performance reached 0.89, demonstrating a significant relationship. Deformation maps revealed by InSAR stacking with and without GACOS corrections showed that GACOS could mainly remove the topography-related and long wavelength signals.

Monitoring the fate of injected CO2 using geodetic techniques

Geodetic methods comprise one class of geophysical data that are sensitive to changes in effective pressure within operating reservoirs, albeit indirectly through induced deformation. Geodetic observations, which have observation intervals that vary from seconds to days, weeks, or months, generally provide more frequent sampling compared to existing geophysical methodologies (such as seismic time-lapse monitoring), which typically invoke repeat times of months to years. These differences in sampling intervals are primarily due to the extensive effort, and hence cost, of conducting geophysical field operations, which often precludes executing a large number of surveys. Satellite-based interferometric synthetic aperture radar (InSAR) is cost effective and used in many applications, including monitoring the injection of carbon dioxide (CO2) for both long-term storage and enhanced oil production. An application to the geologic sequestration of CO2 in Algeria revealed northwest migration along a fault/fracture zone intersected by the injection well. A study in a Texas field demonstrated that enhanced oil recovery utilizing CO2 leads to observable surface deformation that may be used to characterize the sequestered CO2 and to estimate the pressure changes within the reservoir induced by injection and production.

Spatiotemporal Ocean Tidal Loading in InSAR Measurements Determined by Kinematic PPP Solutions of a Regional GPS Network

The coastal crustal deformation caused by ocean tidal loading (OTL) varies spatially and temporally, and this spatiotemporal variation in satellite-based interferometric synthetic aperture radar (InSAR) measurements needs to be determined. In this article, we propose a spatiotemporal modeling method to estimate the OTL displacements in InSAR measurements using the kinematic precise point positioning (PPP) solutions of a regional GPS network. We tested the method through an experiment using 25 Sentinel-1B images and long-term observations of 172 GPS reference sites from Southern California. The experimental results suggest that there are significant OTL and solid Earth tide effects in the differential InSAR interferogram, which is greater than 40 mm. We find that the spatial characteristics of OTL variations can be expressed as a high-order polynomial in the two variables of latitude and longitude, and the spatiotemporally modeled PPP tidal estimates of the high-density GPS sites can provide high precision OTL correction for all the pixels in the interferogram. In the last part of the study, we show that the spatial large-scale signals in the differential interferograms of Sentinel-1B data are mainly atmospheric delay, solid Earth tidal, and OTL effect, and demonstrate the importance of the tidal correction in the InSAR measurements.

Satellite-Based InSAR Monitoring of Highway Bridges: Validation Case Study on the North Channel Bridge in Ontario, Canada

A space-borne bridge displacement monitoring technology was demonstrated on a major highway bridge in Cornwall (Ontario, Canada) over a 15-month period. Two major challenges had to be overcome and solutions were identified and implemented. The first challenge related to the proper selection and analysis of satellite imagery to optimize the number and signal quality of point targets on the bridge for adequate movement detection and analysis. The second challenge was the comparison and validation of independent sets of data having different measurement methods, reference baselines, and distributions in space and time. The key findings include: (i) the regression analysis of interferometric synthetic aperture radar (InSAR) data from the satellite produced the best coherence when fitted against the following three independent variables: height, ambient temperature, and time; (ii) the bridge railings appeared to be excellent natural reflectors with their sharp edges and regular spacing along the bridge, showing a return of over 3,000 clear point targets to the satellite; (iii) the InSAR displacement thermal sensitivity data was found to compare very well to numerical modeling thermal data. The results show great promise and value in applying satellite-based technology for the remote monitoring of highway and railway bridges to alert engineers of excessive movement and, in turn, will help optimize preventive maintenance management, extend structural lifespan, minimize traffic disruptions due to late repair, and ensure structural integrity.

Antarctic ice velocities from GPS locations logged by seismic stations

AbstractIn 2007–08, seismologists began deploying passive seismic stations over much of the Antarctic ice sheet. These stations routinely log their position by navigation-grade global positioning system (GPS) receivers. This location data can be used to track the stations situated on moving ice. For stations along the traverse from Zhongshan station to Dome A in East Antarctica and at the West Antarctic Ice Sheet divide the estimated velocities of the ice surface based on positions recorded by navigation-grade GPS are consistent with those obtained by high-accuracy geodetic GPS. Most of the estimated velocities have an angle difference of <28° with the steepest downhill vector of the ice surface slope at the stations. These results indicate that navigation-grade GPS measurements over several months provide reliable information on ice sheet movement of ≥1 m yr-1. With an uncertainty of ~0.3–1 m yr-1, this method is able to resolve both very slow ice velocities near Dome A and velocities of >100 m yr-1 on Thwaites Glacier. Information on ice velocity at three locations for which no data from satellite-based interferometric synthetic aperture radar are available have also been provided using this method.

The 2011 Hudson volcano eruption (Southern Andes, Chile): Pre-eruptive inflation and hotspots observed with InSAR and thermal imagery

Hudson is one of the most active volcanoes in the Southern Andes—it had one of the largest eruptions of the 20th century in 1991 (VEI = 5) and smaller eruptions in 1971 (VEI = 3), maybe 1973, and 2011 (VEI of 1-2). We use satellite-based interferometric synthetic aperture radar (InSAR) and thermal imagery to characterize the activity of Hudson between 2004 and 2011 and during the 2011 eruption. InSAR data show that the volcano inflated between 2004 and 2010 with a maximum change rate of between 2 and 3 cm/yr—about half of the deformation rate observed during a previous deformation episode from 1993–1999. Inversion for an inflating point source suggests magma accumulation beneath the SW part of the caldera at an average depth of 10 km. This inferred source is deeper than both the sources estimated for the magma chamber of the 1991 eruption (from petrology) and for the 1993–1999 deformation event. Also, the deformation from 2004–2010 is centered at a slightly different location and has a smaller volume change than that between 1993–1999—further indicating that there is either a large magma reservoir or several separate ones. While the deformation center is a few km from the eruption location near the caldera rim, the two are possibly linked since the predicted static Coloumb stress changes due to the inferred inflation source would encourage unclamping on potential faults in the caldera rim. We also analize nighttime satellite thermal images from MODIS and ASTER. While MODIS did not show any unambiguous evidence for hot spots, ASTER thermal imagery show that at least four months before the eruption there were locations with temperatures 7–8oK above background. Lahars observed by helicopter overflights on 4 March 2011 and October 2011 suggest that the hotspots may have been caused by lakes or subglacial melting. There is no InSAR data available for the months immediately preceding the eruption, but the ASTER thermal imagery results may indicate an increase in geothermal activity that could have been used to forecast the eruption.

Earthquake swarms in South America

We searched for earthquake swarms in South America between 1973 and 2009 using the global Preliminary Determination of Epicenters (PDE) catalogue. Seismicity rates vary greatly over the South American continent, so we employ a manual search approach that aims to be insensitive to spatial and temporal scales or to the number of earthquakes in a potential swarm. We identify 29 possible swarms involving 5–180 earthquakes each (with total swarm moment magnitudes between 4.7 and 6.9) within a range of tectonic and volcanic locations. Some of the earthquake swarms on the subduction megathrust occur as foreshocks and delineate the limits of main shock rupture propagation for large earthquakes, including the 2010 Mw 8.8 Maule, Chile and 2007 Mw 8.1 Pisco, Peru earthquakes. Also, subduction megathrust swarms commonly occur at the location of subduction of aseismic ridges, including areas of long-standing seismic gaps in Peru and Ecuador. The magnitude–frequency relationship of swarms we observe appears to agree with previously determined magnitude–frequency scaling for swarms in Japan. We examine geodetic data covering five of the swarms to search for an aseismic component. Only two of these swarms (at Copiapó, Chile, in 2006 and near Ticsani Volcano, Peru, in 2005) have suitable satellite-based Interferometric Synthetic Aperture Radar (InSAR) observations. We invert the InSAR geodetic signal and find that the ground deformation associated with these swarms does not require a significant component of aseismic fault slip or magmatic intrusion. Three swarms in the vicinity of the volcanic arc in southern Peru appear to be triggered by the Mw= 8.5 2001 Peru earthquake, but predicted static Coulomb stress changes due to the main shock were very small at the swarm locations, suggesting that dynamic triggering processes may have had a role in their occurrence. Although we identified few swarms in volcanic regions, we suggest that particularly large volcanic swarms (those that could be detected using the PDE catalogue) occur in areas of infrequent eruption and may be related to large regional fault zones.

Subsidence and uplift of Sidoarjo (East Java) due to the eruption of the Lusi mud volcano (2006–present)

Global positioning system (GPS) and satellite-based InSAR (Interferometric Synthetic Aperture Radar) measurements of the subsidence and uplift of a populated area of Sidoarjo, East Java are due to the eruption of the Lusi mud volcano (2006–present). These data are the first direct quantitative measurements of deformation due to the growth of a mud volcano edifice. The GPS data were recorded over periods of a few hours to several months and show that between June 2006 and September 2007, the earth’s surface has been subsiding at rates of 0.1–4 cm/day. Maximum rates of subsidence occurred in an area 300–400 m to the northwest of the main mud volcano vent. Horizontal displacements were 0.03–0.9 cm/day and were also towards this area. In general uplifts of up to 0.09 cm/day were recorded in areas outside of the edifice. Changes in elevation measured using satellite imagery (InSAR technique) provide regional datasets of subsidence and uplift. They confirm that during the first year a roughly circular area was undergoing sag-like subsidence centered to the northwest of the main vent and that uplift was occurring 3–4 months after the initiation of the eruption due to the movement Watukosek fault system. Subsidence occurred due to the weight of mud and man-made dams and the collapse of the overburden due to removal of mud from the subsurface. Assuming constant rates of subsidence of 4 cm/day, then in the centre of the edifice there would be up to 44 m of subsidence in 3 years, and up to 16 m in 10 years. The mud volcano is now in a self-organizing state with new fluid conduits forming as a result of the collapse.