Synthetic aperture radar interferometry is an imaging technique for measuring the topography of a surface, its changes over time, and other changes in the detailed characteristic of the surface. By exploiting the phase of the coherent radar signal, interferometry has transformed radar remote sensing from a largely interpretive science to a quantitative tool, with applications in cartography, geodesy, land cover characterization, and natural hazards. This paper reviews the techniques of interferometry, systems and limitations, and applications in a rapidly growing area of science and engineering.

/pdf/synthetic-aperture-radar-interferometry-28lm7jf0if.pdf

Synthetic aperture radar interferometry

Geophysical applications of radar interferometry to measure changes in the Earth's surface have exploded in the early 1990s. This new geodetic technique calculates the interference pattern caused by the difference in phase between two images acquired by a spaceborne synthetic aperture radar at two distinct times. The resulting interferogram is a contour map of the change in distance between the ground and the radar instrument. These maps provide an unsurpassed spatial sampling density (∼100 pixels km−2), a competitive precision (∼1 cm), and a useful observation cadence (1 pass month−1). They record movements in the crust, perturbations in the atmosphere, dielectric modifications in the soil, and relief in the topography. They are also sensitive to technical effects, such as relative variations in the radar's trajectory or variations in its frequency standard. We describe how all these phenomena contribute to an interferogram. Then a practical summary explains the techniques for calculating and manipulating interferograms from various radar instruments, including the four satellites currently in orbit: ERS-1, ERS-2, JERS-1, and RADARSAT. The next chapter suggests some guidelines for interpreting an interferogram as a geophysical measurement: respecting the limits of the technique, assessing its uncertainty, recognizing artifacts, and discriminating different types of signal. We then review the geophysical applications published to date, most of which study deformation related to earthquakes, volcanoes, and glaciers using ERS-1 data. We also show examples of monitoring natural hazards and environmental alterations related to landslides, subsidence, and agriculture. In addition, we consider subtler geophysical signals such as postseismic relaxation, tidal loading of coastal areas, and interseismic strain accumulation. We conclude with our perspectives on the future of radar interferometry. The objective of the review is for the reader to develop the physical understanding necessary to calculate an interferogram and the geophysical intuition necessary to interpret it.

/pdf/radar-interferometry-and-its-application-to-changes-in-the-2xi2itqotn.pdf

Radar interferometry and its application to changes in the Earth's surface

We investigate angular velocity vectors of the Philippine Sea (PH) plate relative to the adjacent major plates, Eurasia (EU) and Pacific (PA), and the smaller Caroline (CR) plate. Earthquake slip vector data along the Philippine Sea plate are inverted, subject to the constraint that EU-PA motion equals that predicted by the global relative plate model NUVEL-1. The resulting solution fails to satisfy geological constraints along the Caroline-Pacific boundary: convergence along the Mussau Trench and divergence along the Sorol Trough. We then seek solutions satisfying both the CR-PA boundary conditions and the Philippine Sea slip vector data, by adjusting the PA-PH and EU-PH best fitting poles within their error ellipses. We also consider northern Honshu to be part of the North American plate and impose the constraint that the Philippine Sea plate subducts beneath northern Honshu along the Sagmi Trough in a NNW-NW direction. Of the solutions satisfying these conditions, we select the best EU-PH as 48.2 deg N, 157.0 deg E, 1.09 deg/my, corresponding to a pole far from Japan and south of Kamchatka, and PA-PH, 1.2 deg N, 134.2 deg E, 1.00 deg/my. Predicted NA-PH and EU-PH convergence rates in central Honshu are consistent with estimated seismic slip rates. Previous estimates of the EU-PH pole close to central Honshu are inconsistent with extension within the Bonin backarc implied by earthquake slip vectors and NNW-NW convergence of the Bonin forearc at the Sagami Trough.

A model for the motion of the Philippine Sea Plate consistent with NUVEL‐1 and geological data

Synthetic aperture radar interferometry (InSAR) from Earth-orbiting spacecraft provides a new tool to map global topography and deformation of the Earth's surface. Radar images taken from slightly different viewing directions allow the construction of digital elevation models of meter-scale accuracy. These data sets aid in the analysis and interpretation of tectonic and volcanic landscapes. If the Earth's surface deformed between two radar image acquisitions, a map of the surface dis- placement with tens-of-meters resolution and subcentimeter accuracy can be con- structed. This review gives a basic overview of InSAR for Earth scientists and presents a selection of geologic applications that demonstrate the unique capabilities of InSAR for mapping the topography and deformation of the Earth.

/pdf/synthetic-aperture-radar-interferometry-to-measure-earth-s-1nugx1vz24.pdf

Synthetic aperture radar interferometry to measure earth's surface topography and its deformation

Many large earthquakes are preceded by one or more foreshocks, but it is unclear how these foreshocks relate to the nucleation process of the mainshock. On the basis of an earthquake catalog created using a waveform correlation technique, we identified two distinct sequences of foreshocks migrating at rates of 2 to 10 kilometers per day along the trench axis toward the epicenter of the 2011 moment magnitude (M(w)) 9.0 Tohoku-Oki earthquake in Japan. The time history of quasi-static slip along the plate interface, based on small repeating earthquakes that were part of the migrating seismicity, suggests that two sequences involved slow-slip transients propagating toward the initial rupture point. The second sequence, which involved large slip rates, may have caused substantial stress loading, prompting the unstable dynamic rupture of the mainshock.

https://science.sciencemag.org/content/sci/335/6069/705.full.pdf

Propagation of slow slip leading up to the 2011 M(w) 9.0 Tohoku-Oki earthquake.

Detailed analysis of Global Positioning System data from Japan's Geospatial Information Authority network provides a record of coseismic and postseismic slip distribution on the megathrust fault where the magnitude-9.0 Tohoku-Oki earthquake occurred on 11 March 2011. The coseismic slip area stretches some 400 kilometres along the Japan trench, matching the area of the preseismic locked zone. Afterslip is now overlapping the coseismic slip area and expanding into the surrounding regions. The authors conclude that such geodetic data could help to improve the forecasting of earthquake potential along other subduction zones. In the accompanying News & Views, Jean-Philippe Avouac discusses current models for assessing seismic hazard. Most large earthquakes occur along an oceanic trench, where an oceanic plate subducts beneath a continental plate. Massive earthquakes with a moment magnitude, Mw, of nine have been known to occur in only a few areas, including Chile, Alaska, Kamchatka and Sumatra. No historical records exist of a Mw = 9 earthquake along the Japan trench, where the Pacific plate subducts beneath the Okhotsk plate, with the possible exception of the ad 869 Jogan earthquake1, the magnitude of which has not been well constrained. However, the strain accumulation rate estimated there from recent geodetic observations is much higher than the average strain rate released in previous interplate earthquakes2,3,4,5,6. This finding raises the question of how such areas release the accumulated strain. A megathrust earthquake with Mw = 9.0 (hereafter referred to as the Tohoku-Oki earthquake) occurred on 11 March 2011, rupturing the plate boundary off the Pacific coast of northeastern Japan. Here we report the distributions of the coseismic slip and postseismic slip as determined from ground displacement detected using a network based on the Global Positioning System. The coseismic slip area extends approximately 400 km along the Japan trench, matching the area of the pre-seismic locked zone4. The afterslip has begun to overlap the coseismic slip area and extends into the surrounding region. In particular, the afterslip area reached a depth of approximately 100 km, with Mw = 8.3, on 25 March 2011. Because the Tohoku-Oki earthquake released the strain accumulated for several hundred years, the paradox of the strain budget imbalance may be partly resolved. This earthquake reminds us of the potential for Mw ≈ 9 earthquakes to occur along other trench systems, even if no past evidence of such events exists. Therefore, it is imperative that strain accumulation be monitored using a space geodetic technique to assess earthquake potential.

Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake

[1] We estimated the spatial and temporal evolution of the preceding aseismic slip from January 2003 to January 2011, the coseismic slip of the Tohoku earthquake, and the postseismic slip after the earthquake based on global positioning system (GPS) data. Time-dependent analysis indicates aseismic slip offshore of Miyagi and Fukushima prefectures from 2004 associated with a series of subduction earthquakes that overlap the aseismic slip area. These preceding aseismic and coseismic slip areas are centered between the centers of the coseismic and afterslip areas of the Tohoku earthquake offshore of Miyagi prefecture, while they overlap the coseismic and afterslip areas of the Tohoku earthquake off Fukushima prefecture. The timing of moment magnitude nine (Mw9) -class earthquakes appears to be controlled by multiple preceding slip events, smaller earthquakes and their afterslip. The preceding aseismic and coseismic slip decreased the coupling rate off the Tohoku coast, suggesting the possibility that the preceding slip represented a precursive stage of the Tohoku earthquake. The afterslip of the Tohoku earthquake occurred in an area where the coseismic slip was not large, complementing the large coseismic slip zone. The afterslip along Iwate-Miyagi extends up to 80 km in depth and is currently the sole mechanism of strain release in this depth range. The source region of the anticipated Miyagi-Oki earthquake shows small postseismic slip after the Tohoku earthquake, reflecting the energy release at the time of the earthquake. Aftershock activity is roughly governed by an afterslip process.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2011JB009120

Preceding, coseismic, and postseismic slips of the 2011 Tohoku earthquake, Japan

The 2005 northern Pakistan earthquake (magnitude 7.6) of 8 October 2005 occurred in the northwestern part of the Himalayas. We interpreted landslides triggered by the earthquake using black-and-white 2.5-m-resolution System Pour l’Observation de la Terre 5 (SPOT 5) stereo images. As a result, the counts of 2,424 landslides were identified in the study area of 55 by 51 km. About 79% or 1,925 of the landslides were small (less than 0.5 ha in area), whereas 207 of the landslides (about 9%) were large (1 ha and more in area). Judging from our field survey, most of the small landslides are shallow rock falls and slides. However, the resolution and whitish image in the photos prevented interpreting the movement type and geomorphologic features of the landslide sites in detail. It is known that this earthquake took place along preexisting active reverse faults. The landslide distribution was mapped and superimposed on the crustal deformation detected by the environmental satellite/synthetic aperture radar (SAR) data, active faults map, geological map, and shuttle radar topography mission data. The landslide distribution showed the following characteristics: (1) Most of the landslides occurred on the hanging-wall side of the Balakot–Garhi fault; (2) greater than one third of the landslides occurred within 1 km from the active fault; (3) the greatest number of landslides (1,147 counts), landslide density (3.2 counts/km2), and landslide area ratio (2.3 ha/km2) was found within Miocene sandstone and siltstone, Precambrian schist and quartzite, and Eocene and Paleocene limestone and shale, respectively; (4) there was a slight trend that large landslides occurred on vertically convex slopes rather than on concave slopes; furthermore, large landslides occurred on steeper (30° and more) slopes than on gentler slopes; (5) many large landslides occurred on slopes facing S and SW directions, which is consistent with SAR-detected horizontal dominant direction of crustal deformation on the hanging wall.

Interpretation of landslide distribution triggered by the 2005 Northern Pakistan earthquake using SPOT 5 imagery

An analysis of Global Positioning System (GPS) data revealed south-southeastward transient deformation caused by the expected slow-slip near the Boso peninsula, central Japan, approximately 5 years after a similar event in 2002. An area of aseismic slip with a moment magnitude of 6.6 was estimated off the Boso peninsula, adjacent to the area of associated seismic activity. The 2007 aseismic slip started from around August 10, expanded slightly to the north, and gradually ceased activity over a period of about 10 days. This rupture process is different from those of the last two events in which the slow-slip area moved from north to south. However, the three slow-slip events detected by the GPS array verify the hypothesis that the Boso slow-slips occur quasi-periodically in a certain area, accompanied by seismic swarms. The factors causing the slight differences in the rupture process, magnitude, recurrence time, and slip area among the observed slow-slip events remain unclear.

https://link.springer.com/content/pdf/10.1186%2FBF03352072.pdf

Occurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan

The two-dimensional displacement vector field on the Earth's surface associated with a large earthquake near Iwate volcano in Japan, was mapped by combining satellite radar interferometry images from two different track directions. The displacement reveals not only large main ruptures, but details of many small faults movement of aftershocks and volcanic earthquakes. The displacements of the aftershock faults coincide with a stress field reorganized by the main earthquake that accelerated the aftershock fault movements. The displacements mimic the topography, suggesting that repeated movements of the same faults in the Quaternary have led to the presently existing topography.

https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/1999GL011291

Mikio Tobita

Papers

Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake

Preceding, coseismic, and postseismic slips of the 2011 Tohoku earthquake, Japan

Interpretation of landslide distribution triggered by the 2005 Northern Pakistan earthquake using SPOT 5 imagery

Occurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan

2.5-D surface deformation of M6.1 earthquake near Mt Iwate detected by SAR interferometry