Paul A. Rosen

Synthetic aperture radar interferometry

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.

Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps

Interferogram images derived from repeat-pass spaceborne synthetic aperture radar systems exhibit artifacts due to the time and space variations of atmospheric water vapor. Other tropospheric variations, such as pressure and temperature, also induce distortions, but the effects are smaller in magnitude and more evenly distributed throughout the interferogram than the wet troposphere term. Spatial and temporal changes of 20% in relative humidity lead to 10 cm errors in deformation products, and perhaps 100 m of error in derived topographic maps for those pass pairs with unfavorable baseline geometries. In wet regions such as Hawaii, these are by far the dominant errors in the Spaceborne Imaging Radar-C and X Band Synthetic Aperature Radar (SIR-C/X-SAR) interferometric products. The unknown time delay from tropospheric distortion is independent of frequency, and thus multiwavelength measurements, such as those commonly used to correct radar altimeter and Global Positioning System (GPS) ionospheric biases, cannot be used to rectify the error. In the topographic case, the errors may be mitigated by choosing interferometric pairs with relatively long baselines, as the error amplitude is inversely proportional to the perpendicular component of the interferometer baseline. For the SIR-C/X-SAR Hawaii data we found that the best (longest) baseline pair produced a map supporting 100 m contouring, whereas the poorest baseline choice yielded an extremely noisy topographic map even at this coarse contour interval. In the case of deformation map errors the result is either independent of baseline parameters or else very nearly so. Here the only solution is averaging of independent interferograms, so in order to create accurate deformation products in wet regions many multiple passes may be required. Rules for designing optimal data acquisition and processing sequences for interferometric analyses in nondesert parts of the world are (1) to use the longest radar wavelengths possible, within ionospheric scintillation and Faraday rotation limits, (2) for topography, maximize interferometer baseline within decorrelation limits* and (3) for surface deformation, use multiple observations and average the derived products. Following the above recipe yields accuracies of 10 m for digital elevation models and 1 cm for deformation maps even in very wet regions, such as Hawaii.

Updated repeat orbit interferometry package released

RO1_PAC V2.3, a Repeat Orbit Interferometry package that allows topographic and surface change researchers to apply Interferometric Synthetic Aperture Radar (InSAR) methods, is now freely available to the community InSAR is the synthesis of conventional SAR and interferometry techniques that have been developed over several decades in radio astronomy and radar remote sensing. In recent years, it has opened entirely new application areas for radar in the Earth system sciences, including topographic mapping and geodesy. RO1_PAC, developed primarily to work with European Remote Sensing (ERS) satellite radar data, currently supports ERS-1, ERS-2, and Japanese Earth Resources Satellite (JERS) radar data, and is configurable to work with “strip-mode” data from all existing satellite radar instruments. The first release of RO1_ PAC (V1.0) was made quietly in 2000, and roughly 30 groups in the academic and research community currently use it.

Accuracy of topographic maps derived from ERS-1 interferometric radar

An interferometric radar technique for topographic mapping of surfaces promises a high-resolution approach to the generation of digital elevation models. The authors present analyses of data collected by the synthetic aperture radar instrument on-board the ERS-1 satellite on successive orbits. Use of a single satellite in a nearly repeating orbit is attractive for reducing cost and spaceborne hardware complexity; also it permits inference of changes in the surface from the correlation properties of the radar echoes. The data have been reduced to correlation maps and digital elevation models. The correlation maps show that temporal correlation decreases significantly with time, but not necessarily at a constant well-defined rate, likely depending on environmental factors. When correlation among passes remains high, however, it is possible to form digital elevation models. Analyses of noise expected in ERS-1 interferometric data collected over Alaska and the southwestern United States indicate that maps with relative errors less than 5 m rms are possible in some regions. However, orbit uncertainties imply that tie points are required in order to reduce absolute height errors to a similar magnitude. The authors find that about 6 tie points per 40/spl times/40 km scene with 5 m rms or better height accuracy are needed to keep systematic map height errors below 5 m rms. The performance of the ERS-1 radar system for topographic applications, though useful for a variety of regional and local discipline studies, may be improved with respect to temporal decorrelation errors and absolute height acuity by modifying the orbit repeat period and incorporating precise orbit determination techniques. The resulting implementation will meet many, but not all, objectives of a global mapping mission. >

DETECTION OF AQUIFER SYSTEM COMPACTION AND LAND SUBSIDENCE USING INTERFEROMETRIC SYNTHETIC APERTURE RADAR, ANTELOPE VALLEY, MOJAVE DESERT, CALIFORNIA

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with high spatial detail and resolution of range displacement (±10 mm in change of land surface elevation) were developed for a groundwater basin (∼103 km2) in Antelope Valley, California, using radar data collected from the ERS-1 satellite. These data allow comprehensive comparison between recent (1993–1995) subsidence patterns and those detected historically (1926–1992) by more traditional methods. The changed subsidence patterns are generally compatible with recent shifts in land and water use. The InSAR-detected patterns are generally consistent with predictions based on a coupled model of groundwater flow and aquifer system compaction. The minor inconsistencies may reflect our imperfect knowledge of the distribution and properties of compressible sediments. When used in conjunction with coincident measurements of groundwater levels and other geologic information, InSAR data may be useful for constraining parameter estimates in simulations of aquifer system compaction.

Three-dimensional deformation caused by the Bam, Iran, earthquake and the origin of shallow slip deficit

Our understanding of the earthquake process requires detailed insights into how the tectonic stresses are accumulated and released on seismogenic faults. We derive the full vector displacement field due to the Bam, Iran, earthquake of moment magnitude 6.5 using radar data from the Envisat satellite of the European Space Agency. Analysis of surface deformation indicates that most of the seismic moment release along the 20-km-long strike-slip rupture occurred at a shallow depth of 4–5 km, yet the rupture did not break the surface. The Bam event may therefore represent an end-member case of the ‘shallow slip deficit’ model, which postulates that coseismic slip in the uppermost crust is systematically less than that at seismogenic depths (4–10 km). The InSAR-derived surface displacement data from the Bam and other large shallow earthquakes suggest that the uppermost section of the seismogenic crust around young and developing faults may undergo a distributed failure in the interseismic period, thereby accumulating little elastic strain.

Transient strain accumulation and fault interaction in the Eastern California shear zone

Satellite synthetic aperture radar interferometry reveals transient strain accumulation along the Blackwater‐Little Lake fault system within the Eastern California shear zone. The surface strain map obtained by averaging eight years (1992‐2000) of Earth Resource Satellite (ERS) radar data shows a 120-km-long, 20-kmwide zone of concentrated shear between the southern end of the 1872 Owens Valley earthquake surface break and the northern end of the 1992 Landers earthquake surface break. The observed shear zone is continuous through the Garlock fault, which does not show any evidence of left-lateral slip during the same time period. A dislocation model of the observed shear indicates right-lateral slip at 7 6 3 mm/yr on a vertical fault below ;5 km depth, a rate that is two to three times greater than the geologic rates estimated on northwest-trending faults in the eastern Mojave area. This transient slip rate and the absence of resolvable slip on the Garlock fault may be the manifestation of an oscillatory strain pattern between interacting, conjugate fault systems.

Postseismic Rebound in Fault Step-Overs Caused by Pore Fluid Flow

Near-field strain induced by large crustal earthquakes results in changes in pore fluid pressure that dissipate with time and produce surface deformation. Synthetic aperture radar (SAR) interferometry revealed several centimeters of postseismic uplift in pull-apart structures and subsidence in a compressive jog along the Landers, California, 1992 earthquake surface rupture, with a relaxation time of 270 ± 45 days. Such a postseismic rebound may be explained by the transition of the Poissons ratio of the deformed volumes of rock from undrained to drained conditions as pore fluid flow allows pore pressure to return to hydrostatic equilibrium.

Surface Displacement of the 17 May 1993 Eureka Valley, California, Earthquake Observed by SAR Interferometry

Satellite synthetic aperture radar (SAR) interferometry shows that the magnitude 6.1 Eureka Valley earthquake of 17 May 1993 produced an elongated subsidence basin oriented north-northwest, parallel to the trend defined by the aftershock distribution, whereas the source mechanism of the earthquake implies a north-northeast-striking normal fault. The �3-millimeter accuracy of the radar-observed displacement map over short spatial scales allowed identification of the main surface rupture associated with the event. These observations suggest that the rupture began at depth and propagated diagonally upward and southward on a west-dipping, north-northeast fault plane, reactivating the largest escarpment in the Saline Range.

The 2003 Bam (Iran) earthquake: Rupture of a blind strike-slip fault

A magnitude 6.5 earthquake devastated the town of Bam in southeast Iran on 26 December 2003. Surface displacements and decorrelation effects, mapped using Envisat radar data, reveal that over 2 m of slip occurred at depth on a fault that had not previously been identified. It is common for earthquakes to occur on blind faults which, despite their name, usually produce long-term surface effects by which their existence may be recognised. However, in this case there is a complete absence of morphological features associated with the seismogenic fault that destroyed Bam.