Next Article in Journal
Analyzing Land Use Changes in the Metropolitan Jilin City of Northeastern China Using Remote Sensing and GIS
Next Article in Special Issue
Tomographic Imaging of a Forested Area By Airborne Multi-Baseline P-Band SAR
Previous Article in Journal
Long-term Satellite NDVI Data Sets: Evaluating Their Ability to Detect Ecosystem Functional Changes in South America
Previous Article in Special Issue
Ship Detection in SAR Image Based on the Alpha-stable Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Atmospheric Effects on InSAR Measurements and Their Mitigation

1
Department of Land Surveying and Geo-Informatics, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong, China
2
School of Info-Physics and Geomatics Engineering, Central South University, Changsha 410083, Hunan, China
*
Author to whom correspondence should be addressed.
Sensors 2008, 8(9), 5426-5448; https://doi.org/10.3390/s8095426
Submission received: 28 July 2008 / Revised: 29 August 2008 / Accepted: 29 August 2008 / Published: 3 September 2008
(This article belongs to the Special Issue Synthetic Aperture Radar (SAR))

Abstract

:
Interferometric Synthetic Aperture Radar (InSAR) is a powerful technology for observing the Earth surface, especially for mapping the Earth's topography and deformations. InSAR measurements are however often significantly affected by the atmosphere as the radar signals propagate through the atmosphere whose state varies both in space and in time. Great efforts have been made in recent years to better understand the properties of the atmospheric effects and to develop methods for mitigating the effects. This paper provides a systematic review of the work carried out in this area. The basic principles of atmospheric effects on repeat-pass InSAR are first introduced. The studies on the properties of the atmospheric effects, including the magnitudes of the effects determined in the various parts of the world, the spectra of the atmospheric effects, the isotropic properties and the statistical distributions of the effects, are then discussed. The various methods developed for mitigating the atmospheric effects are then reviewed, including the methods that are based on PSInSAR processing, the methods that are based on interferogram modeling, and those that are based on external data such as GPS observations, ground meteorological data, and satellite data including those from the MODIS and MERIS. Two examples that use MODIS and MERIS data respectively to calibrate atmospheric effects on InSAR are also given.

1. Introduction

Synthetic Aperture Radar (SAR) Interferometry, commonly referred to as InSAR, IFSAR or SARI, is a synthesis of the SAR and the interferometry techniques [1]. InSAR is a powerful technology for topographic and ground surface deformation mapping due to its all-weather and day-and-night imaging capability, wide spatial coverage, fine resolution, and high measurement accuracy. Rogers and Ingalls [2] reported the first application of radar interferometry in Earth-based observations of Venus, while Graham [3] was regarded as the the first to apply an InSAR system to Earth topographic mapping. Airborne and spaceborne InSAR systems were then applied to Earth observation by Zebker and Goldstein [4] and Goldstein et al. [5], respectively. Gabriel et al. [6] first demonstrated the potential of differential InSAR (DInSAR) for centimeter or sub-centimeter level surface deformation mapping over a large area.
Significant progress has been made in further developing InSAR technology in the past two decades with the availability of a vast amount of globally covering SAR images from, e.g., ERS, Radarsat, JERS, Envisat, ALOS and TerraSAR sensors and with a wide range of applications of the technology (e.g., [7-19]). It is expected that InSAR will play a wider and more important role in both research and applications in the future with the advances of the technology and many ambitious SAR missions planed.
InSAR technology, however, has also limitations. One of the most intractable is the effect of the atmosphere (mainly the troposphere and the ionosphere) on repeat-pass InSAR. It is well known that electromagnetic waves are delayed (slowed down) when they travel through the troposphere. The effect often introduces significant errors to repeat-pass InSAR measurements. Massonnet et al. [8] first identified such effects. Since then, some intensive research has been carried out aiming to better understand and mitigate the effects. Zebker et al. [20] reported, for example, that spatial and temporal changes of 20% in the relative humidity of the troposphere could lead up to 10 to 14 cm errors in the measured ground deformations and 80 to 290 m errors in derived topographic maps for baselines ranging from 100 m to 400 m in the case of the SIR-C/X-SAR. A number of researchers have concluded that the tropospheric effects are a limiting factor for wide spread applications of repeat-pass InSAR (e.g., [11, 21-23]).
Contrary to the effects of the troposphere, the ionosphere tends to accelerate the phases of electromagnetic waves when they travel through the medium. The zenith ionospheric range error is proportional to the total electron content (TEC) in the ionosphere. For example, for C-band SAR, a TEC of 1 × 1016 m-2 causes a phase shift of about half a cycle [28]. The ionosphere is however a dispersive medium affecting the radar signals proportionately to the square of the wavelength [83]. For example, if the ionosphere causes 1.5 m range errors to the C-band (wavelength = 5.6 cm) signals, it will cause about 24 m range errors to the L-band (wavelength = 23 cm) signals if the same imaging geometry and atmospheric conditions are assumed. Since there are only very limited published works available on the ionospheric effects on InSAR, we will limit our discussions to the tropospheric effects hereafter. We will review systematically the work carried out in studying the atmospheric, especially the tropospheric effects on InSAR. The basic principles of the atmospheric effects on repeat-pass InSAR are first introduced. Research results on the properties of the atmospheric effects will then be examined. The various methods developed for mitigating the atmospheric effects will finally be studied.

2. Repeat-Pass SAR Interferometry

InSAR can be classified into across- and along-track interferometry according to the interferometric baseline formed, or single- and repeat-pass interferometry according to the number of platform passes involved. Two antennas are mounted on the same platform in along-track interferometry and a single platform pass suffices [24]. Across-track interferometry can be performed either with a one-antenna (e.g., ERS, Envisat) or a two-antenna (e.g., SRTM) SAR system. Revisit to the same scene is required for a one-antenna SAR system so that this is called repeat-pass SAR interferometry [25]. The atmospheric effects in the single-pass interferometry are basically removed completely in the interferometric computation as the effects are almost the same for the two SAR images. In repeat-pass interferometry, however, the atmospheric effects can become significant as the atmospheric conditions can vary considerably between the two SAR acquisitions. We will hereinafter limit our discussions to repeat-pass InSAR only.
The geometrical configuration of repeat-pass SAR interferometry is illustrated in Figure 1a. A1 and A2 are the positions of radar platforms corresponding to the two acquisitions. The phases, ψ1 and ψ2, measured at the two platform positions to a ground point are:
ψ 1 = 4 π λ L 1 , ψ 2 = 4 π λ L 2
where L1 and L2 are the slant ranges and λ is the wavelength of the radar signal. The interferometric phase ϕ is then
ϕ = ψ 1 ψ 2 = 4 π λ ( L 1 L 2 )
Under the far field approximation, one gets
ϕ = ψ 1 ψ 2 4 π λ B = 4 π λ B sin ( θ α )
where α is the orientation angle of the baseline and θ is the look angle.
When assuming a surface without topographic relief as illustrated in Figure 1b, the interferometric phase becomes [11]
ϕ 0 = 4 π λ B sin ( θ 0 α )
where θ0 is the look angle. If topographic relief is present, the look angle will differ from θ0 by δθ,
ϕ = 4 π λ B sin ( θ 0 + δ θ α )
Combining Equations (4) and (5), we get the “flattened” phase
ϕ flat = ϕ ϕ 0 4 π λ B cos ( θ 0 α ) δ θ = 4 π λ B δ θ
The relationship between the topographic height and δθ can be easily established (see Figure 1b)
h L δ θ 0 sin θ 0
Thus the topographic height can be expressed as
h = λ L 4 π B sin θ 0 cos ( θ 0 α ) ϕ flat
The aforementioned process of topography reconstruction is based on the assumption that the imaged surface is stationary during the acquisitions. The interferometric phase in repeat-pass interferometry in fact measures any ground displacement in addition to topography. DInSAR is the technique to extract displacement signature from a SAR interferogram over the acquisition period. In Figure 2, there is an exaggerated ground displacement Δd between the two acquisitions whose projection onto radar line-of-sight (LOS) direction is Δr.
The displacement will introduce a variation of interferometric phase which is proportional to Δr:
Δ ϕ flat = 4 π λ Δ r
Therefore, the interferometric phase includes topography information as well as deformation information,
ϕ flat 4 π λ B cos ( θ 0 α ) L sin θ 0 h + 4 π λ Δ r
To map the ground deformation between two SAR acquistions, the topographic contribution must be removed. According to the ways to remove the topographic contribution, three types of DInSAR configuration can be distinguished: (1) two-pass plus external DEM, (2) three-pass, and (3) four-pass. In two-pass plus external DEM formulation, a SAR interferogram (topographic interferogram thereinafter) is simulated based on the DEM and the imaging geometry of the “real” interferogram (deformation interferogram thereinafter) and is removed from the deformation interferogram. However, in three-pass and four-pass formulations, both the topographic and the deformation interferograms are generated from SAR images. The only difference between them is that in three-pass interferometry, one image is shared by both the topographic and the deformation interferograms. The two-pass plus external DEM and three-pass and four-pass configuration DInSAR can be expressed as:
Δ r two = λ 4 π ( ϕ d ϕ sim , t )
Δ r three , four = λ 4 π ( ϕ d B d B t ϕ t )
where ϕd and ϕt are phases of deformation and topography interferograms, respectively, and B d and B t are perpendicular baseline components of the deformation and topography interferograms, respectively.
The interferometric phase in Equation (10) may also include linear phase ramps caused by orbital errors that should be modeled and removed to derive the ground deformation [22, 28]. This can at times become a problem when the deformation or topography phases also have linear trends. We will however not discuss this problem further in this paper.

3. The Atmosphere and its Effects on Repeat-Pass InSAR

Atmospheric artifacts in SAR interferograms are mainly due to changes in the refractive index of the medium. These changes are mainly caused by the atmospheric pressure, temperature and water vapor. In most cases, the spatial variations of pressure and temperature are not large enough to cause strong, localized phase gradients in SAR interferograms. Their effects are generally smaller in magnitude and more evenly distributed throughout the interferogram when comparing with that of the water vapor, and sometimes difficult to be distinguished from errors caused by orbit uncertainties [22, 26]. The artifact caused by localized water vapor generally dominates the atmosphere induced artifacts in SAR interferograms. Water vapor is mainly contained in the near-ground surface troposphere (up to about 2 km above ground), where a strong turbulent mixing process occurs. Turbulent mixing can result in three-dimensional (3D) spatial heterogeneity in the refractivity and can cause localized phase gradient in both flat and mountainous regions [27, 28]. Besides turbulent mixing, another atmospheric process with clear physical origin is the stratification of the atmosphere. Stratification of the atmosphere into layers of different vertical refractivity causes additional atmospheric delays in mountainous regions [27, 28]. It should be noted that although water vapor is often considered the most important parameter causing the tropospheric delays, there are cases, e.g., in regions with strong topography, changes in pressure between two acquisitions can generate a bigger tropospheric delay signal than humidity variation.
Clouds are formed when the water vapor in the air condenses into visible mass of droplets or frozen crystals. Clouds are divided into two general categories, layered and convective. These are named stratus clouds and cumulus clouds respectively. The liquid water content in the stratiform clouds is usually low so that they do not cause significant range errors to SAR signals. The liquid water content in the cumulus clouds can however range from 0.5 to 2.0 g/m3 and cause zenith delays of 0.7 to 3.0 mm/km [26], significant to InSAR measurements.
Due to the propagation delay of radar signals, in repeat-pass SAR interferometry systems, the phase measurements corresponding to Equation (1) becomes:
ψ 1 = 4 π λ ( L 1 + Δ L 1 ) , ψ 2 = 4 π λ ( L 2 + Δ L 2 )
where ΔL1 and ΔL2 are atmospheric propagation delays of radar signals corresponding to the first and the second acquisitions. This gives the interferometric phase
ϕ = ψ 1 ψ 2 = 4 π λ ( L 1 L 2 ) + 4 π λ ( Δ L 1 Δ L 2 )
where 4 π λ ( L 1 L 2 ) are topography and surface deformation induced interferometric phase, and 4 π λ ( Δ L 1 Δ L 2 ) is the atmosphere induced interferometric phase. From Equation (14), we can see that the atmosphere induced phase errors are easily interpreted as topography or surface deformation.
It is obvious from Equation (14) that it is the relative tropospheric delay (ΔL1 − ΔL2) that causes errors in InSAR measurements. If the atmospheric profiles remain the same at the two acquisitions, the relative tropospheric delay will disappear. In addition, if ΔL1 − ΔL2 = constant for all the resolution cells in an area of interest, the atmospheric effects will also be cancelled out. The two conditions are, however, next to impossible to occur in practice. First, the troposphere, especially the tropospheric water vapor, varies significantly over periods of a few hours or shorter. It is, therefore, highly unlikely to have the same atmospheric profiles even over currently the shortest revisit interval of one day (for ERS-1/ERS-2). Second, it is also rather rare for the relative tropospheric delays to be constant for all the resolution cells due to local tropospheric turbulences, which affect flat terrain as well as mountainous terrain and to vertical stratification which only affects mountainous terrain [27-29].
The influences of the atmosphere induced phase errors on repeat-pass topographic and two-pass surface deformation measurements are straightforward [4, 9]
σ h = λ 4 π B L sin θ σ ϕ
σ Δ r , two = λ 4 π σ ϕ
where σϕ is the phase error in the interferogram; σh is the resultant height error; and σΔr,two is the deformation error for two-pass D-InSAR.
Assuming the same standard deviation σϕ on each interferogram, the covariance matrixes of ϕ = [ϕd ϕt] T for three-pass and four-pass interferometry are:
C o v ϕ , three = [ σ ϕ 2 1 2 σ ϕ 2 1 2 σ ϕ 2 σ ϕ 2 ]
C o v ϕ , four = [ σ ϕ 2 0 0 σ ϕ 2 ]
According to error propagation theorem, the effects of atmosphere induced phase errors on deformation mapping in three-pass and four-pass interferometry are:
σ Δ r , three = λ 4 π 1 B d B t + ( B d B t ) 2 σ ϕ
σ Δ r , four = λ 4 π 1 + ( B d B t ) σ ϕ

4. Properties of Atmospheric Signals in SAR Interferograms

4.1 Atmospheric Signals from SAR Interferograms

A SAR interferogram is a superposition of information on the topography, the surface deformation between the two SAR acquisitions, the differential atmospheric propagation delays between the two SAR acquisitions, and various noise (e.g., [22, 26]). The contribution from the topography can be removed by using a reference elevation model. That from the surface deformation can be neglected or removed if the surface deformation of the study area between the two SAR acquisitions is insignificant or the deformation is known. In addition, multi-looking operations and careful interferometric processing can help to suppress the noise. Therefore at the end an interferogram that contains only the atmospheric signature can be obtained [26]. The atmospheric signature thus obtained is very useful for studying the properties of atmospheric effects on InSAR. Besides, the atmospheric signals can be used to derive various atmospheric products. For example, Hanssen et al. [30] used atmospheric signals derived from SAR interferograms to map high-resolution water vapor.

4.2 Anisotropic Properties of Atmospheric Signals

Radon transform is the projection of image intensities along a radial line at a specified angle. A single Radon transform is a mapping of an image from two dimensions to one dimension where the image intensities collapse to a profile. Radon transform is therefore a tool to investigate anisotropy in images since systematic intensity variations in a particular direction will be visible as a profile [31].
Hanssen [26] first used Radon transform to examine the anisotropy of atmospheric signatures in SAR interferograms, while Jónsson [32] used it to characterize the anisotropy of the noise in SAR interferograms. Li et al. [33] used Radon transform to study the anisotropy of atmospheric signatures in four SAR interferograms over Shanghai. The results are shown in Figure 3.
The Radon transform of atmospheric signals showed varying degrees of anisotropy. For example, the first transform (Figure 3a) showed strong asymmetry especially for profiles of 0° to 90°. This implies that there are areas of very different atmospheric signals in the southwest and northeast corners of the interferogram. However, as the authors pointed out, none of the transforms showed complex variations in the signals, perhaps due to the fact that the studied region is very flat. The results are quite different from those obtained in mountainous regions where the atmospheric effects vary significantly (e.g., [34]) perhaps due to the vertical stratification or the “static” effect of the troposphere in the mountainous regions [29, 35, 43] and the effects of mountains on local weather conditions.

4.3 Gaussianity of Atmospheric Signals

It is important to examine the Gaussianity of atmospheric signals in SAR interferograms as different processing strategies must be applied for Gaussian and non-Gaussian signals. There are a number of hypothesis tests to study whether a signal is Gaussian or non-Gaussian. The Jarque-Bera test is based on classical measures of skewness and kurtosis and it examines whether the sample skewness and kurtosis are unusually different from their expected values [36]. The Hinich test is however a frequency domain test that examines the deviation of the bispectrum of the signal from zero as the bispectrum of a Gaussian signal is zero [37]. Li et al. [33] used both the Jarque-Bera and the Hinich methods to test the atmospheric signals in the four SAR interferograms over Shanghai. The results from both of the methods indicate that the atmospheric signals in all the interferograms are non-Gaussian.

4.4 Spectral Characteristics of Atmospheric Signals

The spectrum of atmospheric signals in a SAR interferogram reveals the energy distribution of the atmospheric effects at different spatial scales. Two-dimensional (2D) FFT is generally used to estimate the 2D power spectra of atmospheric signals. As the power spectra derived can be very noisy, the one-dimensional (1D) rotationally averaged power spectra are usually calculated from the 2D power spectra and used to study the energy distribution of atmospheric signals (e.g., [28, 38]).
Goldstein et al. [39] first calculated the power spectra of atmospheric signals in a SAR interferogram, and demonstrated that the spectra followed a power law distribution with a power exponent of −8/3. This feature is associated with Kolmogorov turbulences, indicating the nature of scale invariance or scaling [40]. Hanssen [26] analyzed the spectra of atmospheric signals in 26 SAR interferograms over Netherland. The results also showed the power law feature. Li et al. [33] calculated the power spectra of atmospheric signals for the four interferograms over Shanghai (Figure 4). It is very clear that the signals follow on the whole the power law distribution. The results are in good agreement with those presented for Mojave desert of California by Goldstein et al. [39] and for the Groningen and Flevoland area of Netherlands by Hanssen [28].
The power law spectral characteristics of the atmospheric signals are very useful. For example, Ferretti et al. [41] used the spectral characteristics to estimate the powers of thermal noise and atmospheric effects, and developed a method based on the results to combine SAR DEMs in wavelet domain. Ferretti et al. [42] also utilized the spectral characteristics to design filters to separate atmospheric effects from nonlinear subsidence. Williams et al. [43] considered that the low-frequency (long wavelength) components of atmospheric effects had larger energy so that sparse external data such as GPS and ground meteorological data can be used to calibrate the effects. Li et al. [44, 45] used the power law nature of the atmospheric effects in designing algorithms to model and correct the effects based on GPS and meteorological data.
The power law can be described by
E ( k ) k β
where E(k) is the power; k is the spatial frequency; and β is the power exponent. The power exponent is an important indicator of data stationarity. Theoretically, when 1 < β < 3, the data series are considered non-stationary but with stationary increments [28, 46]. The estimated spectral exponents range from 2.31 to 2.66 so that the signals have this property. Stationary increments lead to stationary structural function but do not imply that the variance and covariance of the atmospheric signals can be uniquely determined. Therefore, care should be taken when using InSAR data to constrain geophysical models, where the covariances of the noise are generally needed (e.g., [32, 38, 47]).
The 3D Kolmogorov tropospheric turbulence can occur within the region of up to several kilometers in elevation, usually referred to as the effective height of the wet troposphere. The LOS ranges of space-borne SAR systems are much larger than the effective height. The wet tropospheric delays can be therefore typically modeled as 2D turbulence [28].
The power exponent is an important parameter for estimating to what extent the atmospheric effects can be determined and removed with the help of external data. The accumulated energy of the atmospheric signals can be estimated based on the information for different scales by integrating the atmospheric power over the spatial frequencies. The spatial scales corresponding to 90% of the energy thus computed for the four interferograms mentioned above are 0.82, 1.01, 0.94, and 0.29 km, respectively [33]. The spatial scales can be considered as the lowest spatial resolution of external atmospheric data required to calibrate 90% of the atmospheric effects in the SAR interferograms. Therefore, to calibrate 90% of the atmospheric effects for the four interferograms, the spatial resolution of the external atmospheric data (assuming no measurement errors) must reach 0.82, 1.01, 0.94, and 0.29 km, respectively [33]. This can be a reference when one applies corrections for atmospheric effects on InSAR based on external data.

5. Mitigation of Atmospheric Effects on Repeat-Pass InSAR

5.1 Correction of Atmospheric Effects based on External data

5.1.1 Correction of atmospheric effects based on ground meteorological observations

Using ground meteorological data to calibrate tropospheric delays in radio ranging has been well documented [48-52]. Hanssen and Feijt [53] and Zebker et al. [20] proposed the use of the Saastamoinen model to assess the potential effects of the troposphere on InSAR measurements. Delacourt et al. [29] presented a case study of correcting atmospheric errors by using meteorological observations at a reference point together with tropospheric delay models, vertical gradient models of the meteorological parameters, and the DEM of the study area. The results showed that tropospheric corrections reached 2 fringes for some interferograms, and that on average the accuracy of the interferograms was about ±1 fringe after the corrections were applied. Bonforte et al. [54] demonstrated congruence between the tropospheric zenith delays estimated from GPS observations and from tropospheric models and meteorological data. The results confirmed that the meteorological data could be applied to calibrate InSAR measurements like what had been done to correct GPS observations, and suggested a possible integration of the two data sources for improving models of the atmospheric effects. Li et al. [45] studied InSAR atmospheric correction by using meteorological observations, GPS observations, and both types of observations. The results showed that the integration of the observations produced better results.
The difficulties in using meteorological data to correct atmospheric effects on InSAR include mainly the poor accuracy of the atmospheric delays estimated from empirical tropospheric models and the usually very sparse distribution of meteorological stations.

5.1.2 Correction of atmospheric effects based on GPS observations

The advances in GPS meteorology have enabled accurate estimation of tropospheric delays from GPS observations [55, 56], and have provided an opportunity to use GPS observations to evaluate and calibrate the atmospheric effects on InSAR measurements. However, the spatial resolution of GPS stations is in general much lower than that of InSAR data. This poses a potential limitation in applying GPS observations to correcting InSAR measurements.
Considering the power law nature of the atmospheric noise, Williams et al. [43] however dismissed the belief that the spatially sparse GPS observations (compared to the scales of the atmospheric irregularities and the resolutions of SAR data) were unsuitable for calibrating the atmospheric effects. Using simulated data, the authors demonstrated that in general it is possible to use sparsely distributed data to reduce the noise in a more densely distributed data set, and that in particular it is possible to use zenith delays estimated from GPS observations to reduce the atmospheric noise in InSAR measurements. Bock and Williams [57] reported through a cross validation analysis that using zenith delays estimated from GPS observations and the Kriging interpolator, more than 90% of the atmospheric delays at the unsampled points in a SAR image can be retrieved and therefore removed for the Los Angeles basin where fairly dense GPS stations had been in operation. On the other hand, only 40% of the atmospheric delays can be retrieved for regions outside the basin where the density of GPS stations is much lower. Also using cross validation analysis, Janssen et al. [58] tested the effectiveness of three interpolators, i.e., inverse distance weighting, Kriging and spline, in interpolating the GPS-derived atmospheric delays to the SAR resolution level and correcting the atmospheric effects on InSAR on a pixel-by-pixel basis. The results showed that the inverse distance weighting and Kriging interpolators are better than the spline interpolator. Webley et al. [59] proposed a procedure to use the water vapor delays derived from both the GPS observations and the non-hydrostatic three-dimensional (NH3D) meteorological model to calibrate the atmospheric effects on InSAR.
The research results in [43, 57-59] are mainly from cross validation analysis, but not from corrections to real SAR interferograms. Li et al. [44] recently proposed a new method and applied it to correct a SAR interferogram. In this method, an atmospheric delay map for each SAR acquisition is generated in two steps. First, a “mean” atmospheric delay map is calculated using the method adopted by Delacourt et al. [29]. Second, the “mean” atmospheric delay map is amended with the atmospheric zenith delays derived from a dense GPS network, mainly to calibrate the estimated “mean” atmospheric delays and to compensate their horizontal heterogeneity. Using 14 GPS stations over Mt. Etna, the authors corrected a SAR interferogram and achieved 27.2% overall improvement in the accuracy of the InSAR measurements. Based on the variance model of water vapor delays derived by Emardson et al. [47], a linear interpolator and the best linear unbiased estimator, Li et al. [60] developed a GPS topography-dependent turbulence model for InSAR atmospheric correction. Test results show that the model is much better than the inverse distance weighting interpolator.
There are many GPS networks around the world operated in continuous mode. If GPS observations prior to and after SAR acquisitions are available, they can also contribute to the correction of atmospheric effects on InSAR. Onn [62] and Onn and Zebker [61] applied this method to InSAR atmospheric correction based on Taylor's “frozen-flow” hypothesis. The results showed that additional improvement can be achieved when GPS observations prior to and after SAR acquisitions are added to the GPS-based InSAR atmospheric correction models.
The various methods proposed to date that use GPS observations to correct InSAR atmospheric effects differ primarily in the algorithms used to generate atmospheric delay maps from the spatially sparse GPS atmospheric delay measurements. Therefore, the accuracy of the corrections depends on how much atmospheric delays can be retrieved at the unsampled locations from the sparse GPS measurements. With the gradual increase in the density of GPS networks around the world, the method should become more and more useful.

5.1.3 Correction of atmospheric effects based on high-resolution meteorological models

Numerical meteorological modeling is an essential tool in atmospheric research. Numerical meteorological modeling can be carried out on global, regional or mesoscale. The global and regional numerical meteorological models are usually too coarse to model atmospheric effects on InSAR. The mesoscale numerical models can have a continuous time scale and a horizontal spatial scale of a few kilometers, and are therefore suitable for InSAR atmospheric correction. Integrated water vapor content along the radar paths can be retrieved from such models and used to calibrate atmospheric effects on InSAR.
Wadge et al. [59] used the local-scale non-hydrostatic three-dimensional models (NH3D) to simulate the path delays due to water vapor over Mount Etna, and found that the NH3D delays were in general agree well with the ERS-2 SAR interferogram and the GPS estimates. Webley [65] and Webley et al. [64] tested correcting atmospheric effects on a descending and two ascending SAR interferograms over Mount Etna by using the path delays derived from the NH3D models. The results showed that the correction can result in up to 28.6% improvements in terms of the phase standard deviations. The accuracy improvement is however highly dependent on the data used to initialize the NH3D models. Foster et al. [66] used the MM5 models (a non-hydrostatic mesoscale meteorological model produced by the National Center for Atmospheric Research (NCAR)/Pennsylvania State University) to predict the atmospheric delay maps and then to correct 44 SAR interferograms over Hawaii. The results showed that on average atmospheric effects with wavelengths of 30 km or greater can be significantly reduced, while those with wavelengths shorter than 30 km cannot be effectively reduced. More recently, Puysségur et al. [67] found that water vapor content estimated from Envisat Medium Resolution Imaging Spectrometer (MERIS) and from MM5 model were consistent and unbiased, and thus proposed to integrate MM5 model and MERIS data for InSAR atmospheric correction. Test results showed that about 43% of the atmospheric signals can be removed. High-resolution meteorological models have offered some promising opportunities for mitigating atmospheric effects on InSAR although further research needs to be carried out to enhance the accuracy and reliability of the method.

5.1.4 Correction of atmospheric effects based on MODIS data

The near-IR water vapor products provided by the Moderate Resolution Imaging Spectroradiometer (MODIS) have a spatial resolution of 1 km× 1 km (at nadir) and an accuracy of 5-10% [68]. The high-resolution water vapor products appear to be very useful for modeling and correcting atmospheric effects on InSAR although as an optical sensor MODIS measurements are sensitive to the presence of clouds. The resolutions of even the densest GPS networks in the world, e.g., the Southern California Integrated GPS Network (SCIGN), are more than ten times sparser than the resolution of MODIS.
Li et al. [69] first presented some results of using MODIS data to correct atmospheric effects on InSAR over Mount Etna and Los Angeles. Li et al. [70] proposed an integration of MODIS and GPS data for InSAR atmospheric correction, where the GPS data (more exactly the GPS precipitable water vapor (PWV) data) are mainly used to calibrate the MODIS PWV data. Experiments over Los Angeles area showed that the atmospheric singals in the SAR interferograms were significantly reduced with this method and the geophysical signals in the InSAR measurements became more prominent after the corrections were made.
Considering that all data interpolators unavoidably suffer from smoothing effects [71], Li [72] proposed a hybrid algorithm that jointly uses the Kriging interpolator and the conditional spectral simulation method to interpolate the MODIS PWV in correcting the atmospheric signals in SAR interferograms.
Figure 5 shows the original and the corrected ERS-2 interferogram over Los Angeles by using the MODIS data and the developed hybrid algorithm. Note that the original interferogram has been corrected for topographic phases with a known DEM and for deformation phases with GPS positioning results from SCIGN [72]. Thus, the signals left in the original interferogram can be considered solely from the atmospheric effects. After the corrections were applied, the negative atmospheric phases in the southwestern part of the interferogram and the positive atmospheric phases along the eastern margin of the interferogram were largely removed. The phases at the lower central part of the interferogram became however more significant, coupled with some under-modeled positive/negative residual phases in the corners of the interferogram. The phase standard deviation in the original deformation-free interferogram (Figure 5a) is 14.9 mm, while this becomes 10.6 mm in the corrected interferogram (Figure 5b), representing an improvement of 28.9% in the measurement accuracy.
The MODIS PWV measurements however are sensitive to the presence of clouds as noted earlier, which limits significantly the use of MODIS PWV in cloudy regions. In addition, systematic biases in space borne MODIS PWV measurements may exist and need to be calibrated with more accurate PWV measurements (e.g., GPS PWV).

5.1.5 Correction of atmospheric effects based on MERIS data

The MERIS onboard the Envisat satellite allows for global retrieval of PWV every three days, with two near infrared water vapor channels. It therefore can acquire water vapor data simultaneously with the Advanced SAR (ASAR). Its PWV measurements have a resolution as high as 300 m and accuracy higher than that of MODIS [69]. MERIS measurements therefore offer an opportunity for the atmospheric effects on ASAR measurements to be accurately modeled.
Li et al. [73] assessed the potential of using MERIS near-infrared water vapor products to correct ASAR interferometric measurements. The MERIS and the GPS/radiosonde water vapor products tested agreed to each other to within 1.1 mm (standard deviation) on average. It was also pointed out that the major limitation with the use of MERIS water vapor products is the low frequencies of cloud free conditions, i.e., about 25% globally although for certain areas like Easter Tibet and Southern California, the frequencies can be much higher.
Using the Los Angeles area as an example, Li et al. [74] showed that MERIS water vapor data could significantly reduce atmospheric effects in SAR interferograms. After corrections were made with the MERIS data, the RMS difference between GPS and InSAR range changes in the satellite LOS direction decreased from 0.89 cm to 0.54 cm in one interferogram, and from 0.83 cm to 0.59 cm in another. Puysségur et al. [67] proposed the integration of MM5 simulated water vapor data and MERIS data for InSAR atmospheric correction, as noted earlier. However, no significant improvements were found by adding the MERIS data to the MM5 model.
Figure 6 shows an example of correcting atmospheric effects on InSAR using MERIS data over Hong Kong region. The SAR images are acquired on 30 April 2006 and 11 March 2007, respectively. Topographic phases have been removed with a reference DEM. GPS positioning results have shown that there were no significant deformations during this period in region. The signals in Figure 6a can therefore be considered from atmospheric heterogeneity only. Figure 6b shows the corrected interferogram with reduced resolution (RR) MERIS water vapor data. It can be seen that the positive atmospheric phases in some of the areas have been significantly removed. The standard deviation of the phases in the original interferogram (Figure 6a) is 5.72 mm and that in the corrected interferogram (Figure 6b) is 4.51 mm, indicating an improvement of about 21% after the corrections were made.

5.2 Correction of Atmospheric Effects based on Correlation Analysis

There are mainly two types of correlation analysis adopted to reduce atmospheric effects on InSAR. The first type analyzes the correlation between interferograms, and the second the correlation between atmosphere-induced interferometric phases and elevation. Sarti et al. [75] proposed to characterize the atmospheric artifacts in SAR interferograms and to remove them through analyzing the correlation between interferograms. Fruneau and Sarti [76] proposed to separate the deformation signals from atmospheric artifacts by exploiting the correlation between interferograms. The method does not aim to remove atmospheric noise from a SAR interferogram, but manages to extract the common (correlated) deformation signals within two SAR interferograms through analyzing the correlation of the signals. Using this method, the authors successfully extracted the deformation signals from interferograms over Paris. Sarti et al. [77] compared this method with other methods for atmospheric effect mitigation like the stacking and the persistent scatterer (or permanent scatterer) InSAR (PSInSAR) method and pointed out the advantages of the correlation analysis method when the number of available SAR images is not large.
Beauducel et al. [35] proposed to separate deformation signals from atmospheric artifacts over Mount Etna by analyzing the correlation between the atmosphere-induced interferometric phases and the elevations. Using 238 interferograms over the area, the authors jointly estimated the deformations and the tropospheric delays. The results revealed that the estimated large-scale deformation and magma evolution from this study were much less than those from other studies, perhaps due to the fact that the atmospheric artifacts (ranging from -2.7 to +3.0 fringes) had been better accounted for in this study. Using ten SAR interferograms over Sakurajima volcano, Remy et al. [78] carefully investigated the relationship between the atmosphere-induced interferometric phases and the elevations, and found that the non-linear piecewise polynomial form of cubic splines was better in modeling the atmospheric delays than the linear models.
Chaabane et al. [81] suggested using the correlation between interferograms and that between atmosphere-induced interferometric phases and the elevation to correct for the atmospheric effects. In this approach, the global-scale atmospheric contribution is corrected by exploiting the correlation between the interferometric phases and the elevation, while the local atmospheric artifacts are corrected based on the correlation between interferograms containing a common acquisition. Test results with 81 differential interferograms covering the Gulf of Corinth (Greece) show that (1) the average uncertainty of the stacked deformation map has been decreased from ± 26 mm to ±12 mm, and (2) the RMS value of the differences between InSAR and GPS measurements at four stations has decreased from ±30 mm to ±19 mm after applying the correction.
The method of correlation analysis is advantageous in that no external data are needed. The method however strongly depends on the correlations between the deformations and between the atmospheric signals in different interferograms. Weak correlation may lead to insufficient atmospheric effect reduction.

5.3 Correction of Atmospheric Effects based on Pair-Wise Logic

The atmospheric signature in a SAR interferogram can be determined with the pair-wise logic method [21]. Atmospheric perturbations that are different from the pattern of local ground displacements can be identified by comparing interferograms spanning different time intervals. The method was used to find a 25×20 km kidney-shaped feature caused by ionospheric perturbations [8, 21]. Massonnet and Feigl [21] also found irregular patterns of up to three complete fringes resulted from tropospheric turbulences or increased water vapor over a 5×10 km area with this method. The qualitative nature of this method however makes it difficult to give exact values of the atmospheric effects. Hanssen [28] therefore suggested to sum or subtract two interferograms that use a common SAR image for removing atmospheric anomalies. The approach has also been referred to as the method of linear combination. It is effective when the atmospheric anomalies exist only in the common SAR image of the two interferograms.

5.4 Correction of Atmospheric Effects based on PSInSAR Technique

PSInSAR is a relatively new interferometic processing method [42, 79, 80]. It works on temporally stable coherent targets (permanent scatterers) only and can overcome the difficulties of coherence loss and atmospheric heterogeneities in conventional SAR interferometry. In PSInSAR, the atmospheric effects are modeled as linear phase ramps in the azimuth and the range directions for small ground areas or a more sophisticated model that includes, e.g., the linear ramps as well as the topography dependent term and the turbulence can be used for large rugged ground areas [82]. Parameters of an atmospheric model are estimated jointly with other unknowns such as the DEM errors and the LOS ground deformations at the permanent scatterers. The estimated atmospheric effects corresponding to each interferogram are then resampled onto the image grid with an interpolator and removed from the interferogram. Ferretti et al. [42, 79] reported that improved estimation of local topography and terrain motions was resulted over Ancona, Italy and Pomona, California with this method. Hooper et al. [80] modified the PSInSAR algorithms and applied the method to study the temporal and spatial deformation of volcanoes. A shortcoming of the method is that a significant number of SAR images over the same area, typically over 30, are needed to get reliable results.

5.5 Reduction of Atmospheric Effects with the Stacking Method

Stacking is a method that reduces the atmospheric effects on InSAR by averaging independent SAR interferograms. Assuming that atmospheric effects are uncorrelated between the interferograms, averaging N independent interferograms will reduce the atmospheric signals to 1 / N fold. The method was once regarded as the only viable solution to the problem of atmospheric effect mitigation [20]. Williams et al. [43] considered the method of interferogram stacking and that of atmospheric effect calibration with the assistance of external data (such as continuously operating GPS) to be complementary and suggested the two to be used simultaneously. Ferretti et al. [41] proposed a weighted averaging method by taking into account the spectral features of the thermal noise and the atmospheric component. Stacking in general degrades the temporal resolution of InSAR measurements, and the method works when there are only linear ground deformations as non-linear deformations can be lost in the process of stacking.
We have in this section looked through the various existing methods for mitigating the atmospheric effects on InSAR measurements. It should however be pointed out that in principle an optimal integration of some of the methods should yield the best results.

6. Conclusions

Atmospheric effects are one of the limiting error sources in repeat-pass InSAR measurements. They can introduce errors of over ten centimeters to ground deformations and of several hundred meters to DEMs measured with the conventional DInSAR method when considering the typical baseline geometries used. Studies have shown that atmospheric signals in SAR interferograms are anisotropic and non-Gaussian in distribution. The spectra of the atmospheric signals follow a power law distribution with the power exponent very close to -8/3. Various methods have been developed for mitigating the atmospheric effects on InSAR measurements based on external data such as ground meteorological observations, GPS data, satellite water vapor products such as those from MERIS and MODIS, and results from numerical meteorological modeling. These methods are typically able to reduce the atmospheric effects by about 20-40 percents. The other methods developed for mitigating the atmospheric effects are mainly based on simple data analysis or numerical solutions, including the pair-wise logic, the stacking, the correlation analysis, and the PSInSAR methods. Each of the methods developed has its pros and cons. The most suitable method should be chosen considering the number of SAR scenes acquired, the method used for InSAR processing, the atmospheric conditions (e.g., cloud conditions) and the external data available. Despite the progress already made in the research, further studies are still necessary in the area to develop more effective methods for the mitigation of the atmospheric effects.

Acknowledgments

The work presented was supported by the Research Grants Council of the Hong Kong Special Administrative Region (Project Nos.: PolyU 5157/05E and PolyU 5161/06E), the National Science Foundation of China (Project Nos.: 40774003 and 40404001), the Hong Kong Polytechnic University (Project No.: G-YX48), the National High-Tech. Program of China (Project No.: 2006AA12Z156), and West China 1:50000 Topographic Mapping Project. Some of the images used were provided by the European Space Agency under a Category 1 User Project (AO-4458, 4914).

References and Notes

  1. Rosen, P.A.; Hensley, S.; Joughin, I.R.; Li, F.K.; Madsen, S.N.; Rodriguez, E.; Goldstein, R.M. Synthetic Aperture Radar Interferometry. Proc. IEEE 2000, 88, 333–382. [Google Scholar]
  2. Rogers, A.E.; Ingalls, R.P. Venus: Mapping the surface reflectivity by radar interferometry. Science. 1969, 165, 797–799. [Google Scholar]
  3. Graham, L.C. Synthesis interferometric radar for topographic mapping. Proc. IEEE 1974, 62, 763–768. [Google Scholar]
  4. Zebker, H.A.; Goldstein, R.M. Topographic mapping from Interferometric Synthetic Aperture Radar Observations. J. Geophys. Res. 1986, 91, 4993–4999. [Google Scholar]
  5. Goldstein, R.M.; Zebker, H.A.; Werner, C.L. Satellite radar interferometry: Two-dimensional phase unwrapping. Radio Sci. 1988, 23, 713–720. [Google Scholar]
  6. Gabriel, A.K.; Goldstein, R.M.; Zebker, H.A. Mapping small elevation changes over large areas: differential radar interferometry. J. Geophys. Res. 1989, 94, 9183–9191. [Google Scholar]
  7. Massonnet, D.; Rossi, M.; Carmona, C.; Adragna, F.; Peltzer, G.; Feigl, K.L.; Rabaute, T. The displacement field of the Landers earthquake mapped by radar interferometry. Nature 1993, 364, 138–142. [Google Scholar]
  8. Massonnet, D.; Feigl, K.L.; Rossi, M.; Adragna, F. Radar interferometry mapping of deformation in the year after the Landers earthquake. Nature 1994, 369, 227–230. [Google Scholar]
  9. Zebker, H.A.; Werner, C.L.; Rosen, P.A.; Hensley, S. Accuracy of Topographic Maps Derived from ERS-1 Interferometric Radar. IEEE T. Geosci. Remote Sens. 1994, 33, 823–836. [Google Scholar]
  10. Zebker, H.A.; Rosen, P.A.; Goldstein, R.M.; Gabriel, A.; Werner, C.L. On the derivation of coseismic displacement fields using differential radar interferometry: The Landers Earthquake. J. Geophys. Res. 1994, 99, 19617–9634. [Google Scholar]
  11. Rosen, P.A.; Hensley, S.; Zebker, H.A.; Webb, F.H.; Fielding, E.J. Surface deformation and coherence measurements of Kilauea Volcano, Hawaii, from SIR-C radar interferometry. J. Geophys. Res. 1996, 101, 23109–23125. [Google Scholar]
  12. Lanari, R.; Fornaro, G.; Riccio, D.; Migliaccio, M.; Papathanassiou, K.P.; Moreira, J.R.; Schwabisch, M.; Dutra, L.; Puglisi, G.; Franceschetti, G.; Coltelli, M. Generation of digital elevation models by using SIR-C/X-SAR multifrequency two-pass interferometry: the Etna case study. IEEE T. Geosci. Remote Sens. 1996, 34, 1097–1114. [Google Scholar]
  13. Small, D. Generation of Digital Elevation Models through Spaceborne SAR Interferometry. Ph.D dissertation, University of Zurich, Zurich, 1998. [Google Scholar]
  14. Tobita, M.; Fujiwara, S.; Ozawa, S.; Rosen, P.A.; et al. Deformation of the 1995 North Sakhalin earthquake detected by JERS-1/SAR interferometry. Earth Planets Space 1998, 50, 313–335. [Google Scholar]
  15. Fujiwara, S.; Rosen, P.A.; Tobita, M.; Murakami, M. Crustal deformation measurements using repeat-pass JERS-1 Synthetic Radar Interferometry near the Izu Peninsula, Japan. J. Geophys. Res. 1998, 103(B2), 2411–2426. [Google Scholar]
  16. Rufino, G.; Moccia, A.; Esposito, S.D. EM generation by means of ERS Tandem Data. IEEE T. Geosci. Remote Sens. 1998, 36, 1905–1913. [Google Scholar]
  17. Liu, G.X.; Ding, X.L.; Chen, Y.Q.; Li, Z.L.; Li, Z.W. Monitoring land subsidence at the Chek Lap Kok Airport Using InSAR Technology. Chinese Sci. Bull 2001, 46, 1778–1782. [Google Scholar]
  18. Liu, G.X.; Ding, X.L.; Li, Z.L.; Li, Z.W.; Chen, Y.Q.; Yu, S.B. Pre- and Co-Seismic Ground Deformations of the 1999 Chi-Chi, Taiwan Earthquake, Measured with SAR Interferometry. Comput. Geosci. 2004, 30, 333–343. [Google Scholar]
  19. Ding, X.L.; Liu, G.X.; Li, Z.W.; Li, Z.L.; Chen, Y.Q. Ground subsidence monitoring in Hong Kong with satellite SAR interferometry. Photogramm. Eng. Rem. S. 2004, 70, 1151–1156. [Google Scholar]
  20. Zebker, H. A.; Rosen, P. A.; Hensley, S. Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps. J. Geophys. Res. 1997, 102, 7547–7563. [Google Scholar]
  21. Massonnet, D.; Feigl, K. Discrimination of geophysical phenomena in satellite radar interferograms. Geophys. Res. Lett. 1995, 22, 1537–1540. [Google Scholar]
  22. Tarayre, H.; Massonnet, D. Atmospheric propagation heterogeneities revealed by ERS-1. Geophys. Res. Lett. 1996, 23, 989–992. [Google Scholar]
  23. Li, Z.W.; Ding, X.L.; Liu, G.X.; Huang, C. Atmospheric Effects on InSAR Measurements - A Review. Geom. Res. Austr. 2003, 79, 43–58. [Google Scholar]
  24. Gens, R.; van Genderen, J.L. SAR interferometry- Issues, techniques, applications. Int. J. Remote Sens. 1996, 17, 1803–1835. [Google Scholar]
  25. Li, F.K.; Goldstein, R.M. Studies of multibaseline spaceborne interferometric synthetic aperture radars. IEEE T. Geosci. Remote Sens. 1990, 28, 88–97. [Google Scholar]
  26. Hanssen, R. Atmospheric heterogeneities in ERS tandem SAR interferometry; DEOS Report No.98.1; Delft University press: Delft, the Netherlands, 1998. [Google Scholar]
  27. Massonnet, D.; Feigl, K. Radar interferometry and its application to changes in the earth's surface. Rev. Geophys. 1998, 36, 441–500. [Google Scholar]
  28. Hanssen, R.F. Radar Interferometry: Data interpretation and Error Analysis; Kluwer Academic Publishers: Dordrecht, 2001. [Google Scholar]
  29. Delacourt, C.; Briole, P.; Achache, J. Tropospheric corrections of SAR interferograms with strong topography: application to Etna. Geophys. Res. Lett. 1998, 25, 2849–2852. [Google Scholar]
  30. Hanssen, R.F.; Wechwerth, T.M.; Zebker, H.A.; Klees, R. High-resolution water vapor mapping from interferometric radar measurements. Science 1999, 283, 1297–1299. [Google Scholar]
  31. Bracewell, R.N. Two-Dimensional Imaging; Prentice Hall: New Jersey, 1995. [Google Scholar]
  32. Jónsson, S. Modeling Volcano and Earthquake Deformation From Satellite Radar Interferometric Observations. Ph.D dissertation, Stanford University, 2002. [Google Scholar]
  33. Li, Z.W.; Ding, X.L.; Huang, C.; Zou, Z.R. Atmospheric effects on repeat-pass InSAR measurements over Shanghai region. J. Atmos. Terr. Phys. 2007, 69, 1344–1356. [Google Scholar]
  34. Ding, X.L.; Ge, L. L.; Li, Z.W.; Rizos, C. Atmospheric Effects on InSAR Measurements in Southern China and Australia: A Comparative Study. International Archives of the Photogrammetry, Remote Sensing and Spatial Information Science, Istanbul, Turkey, 12- 23 July. 2004; XXXV, pp. 70–75.
  35. Beauducel, B.; Briole, P.; Froger, J.L. Volcano-wide fringes in ERS synthetic aperture radar interferograms of ETNA (1992-1998): Deformation or tropospheric effect? J. Geophys. Res. 2000, 105, 16391–16402. [Google Scholar]
  36. Jarque, C.M.; Bera, A.K. Efficient tests for normality, homoscedasticity and serial independence of regression residuals. Econ. Lett. 1980, 63, 255–259. [Google Scholar]
  37. Hinich, M.J.; Wilson, G. R. Detection of Non-Gaussian Signals in Non-Gaussian Noise Using the Bispectrum. IEEE T. Acoust. Speech. 1990, 38, 1126–1131. [Google Scholar]
  38. Lohman, R.B.; Simons, M. Some thoughts on the use of InSAR data to constrain models of surface deformation: Noise structure and data downsampling. Geochem. Geophys. Geosys. 2005, 6, Q01007. [Google Scholar]
  39. Goldstein, R.M. Atmospheric limitations to repeat-track radar interferometry. Geophys. Res. Lett. 1995, 22, 2517–2520. [Google Scholar]
  40. Tatarski, V.I. Wave propagation in a turbulent medium; McGraw-Hill: New York, 1961. [Google Scholar]
  41. Ferretti, A.; Prati, C.; Rocca, F. Multibaseline InSAR DEM Reconstruction: The Wavelet Approach. IEEE T. Geosci. Remote Sens. 1999, 37, 705–715. [Google Scholar]
  42. Ferretti, A.; Prati, C.; Rocca, F. Nonlinear Subsidence Rate Estimation Using Permanent Scatters in Differential SAR Interferometry. IEEE T. Geosci. Remote Sens. 2000, 38, 2202–2212. [Google Scholar]
  43. Williams, S.; Bock, Y.; Fang, P. Integrated satellite interferometry: Tropospheric noise, GPS estimates and implications for interferometric synthetic aperture radar product. J. Geophys. Res. 1998, 103, 27051–27068. [Google Scholar]
  44. Li, Z. W.; Ding, X.L.; Huang, C.; Wadge, G.; Zheng, D.W. Modeling of atmospheric effects on InSAR measurements by incorporating terrain elevation information. J. Atmos. Terr. Phys. 2006, 66, 1189–1194. [Google Scholar]
  45. Li, Z. W.; Ding, X.L.; Liu, G.X. Modeling Atmospheric Effects on InSAR with Meteorological and Continuous GPS Observations: Algorithms and Some Test results. J. Atmos. Terr. Phys. 2004, 66, 907–917. [Google Scholar]
  46. Davis, A.; Marshak, A.; Wiscombe, W.; Cahalan, R. Scale invariance of liquid water distribution in marine stratocumulus. Part I: spectral properties and stationarity issues. J. Atmos. Sci. 1996, 53, 1538–1558. [Google Scholar]
  47. Emardson, T.R.; Simons, M.; Webb, F.H. Neutral atmospheric delay in interferometric synthetic aperture radar applications: Statistical description and mitigation. J. Geophys. Res. 2003, 108, 2231–2238. [Google Scholar]
  48. Hopfield, H. S. Tropospheric effect on electromagnetically measured range: Prediction from surface weather data. Radio Sci. 1971, 6, 357–367. [Google Scholar]
  49. Saastamoinen, J. Atmospheric correction for troposphere and stratosphere in radio ranging of satellites. The use of artificial satellites for Geodesy. Gephys. Monogr. Ser. 1972, 247–251. [Google Scholar]
  50. Askne, J.; Nordius, H. Estimation of tropospheric delay for microwaves from surface weather data. Radio Sci. 1987, 22, 379–386. [Google Scholar]
  51. Baby, H.B.; Gole, P.; Lavergnat, J. A model for the tropospheric excess path length of radio waves from surface meteorological measurements. Radio Sci. 1988, 23, 1023–1038. [Google Scholar]
  52. Li, Z.W.; Ding, X.L.; Chen, W.; Liu, G.X.; Shea, Y.K.; Emerson, N. Comparative study of tropospheric empirical model for Hong Kong region. Surv. Rev. 2008, in press. [Google Scholar]
  53. Hanssen, R.; Feijt, A. A first quantitative evaluation of atmospheric effects on SAR interferometry. Fringe 96′ Workshop on ERS SAR Interferometry, Zurich, Switzerland, 30 Sep.-2 Oct. 1996; pp. 277–282. [Google Scholar]
  54. Bonforte, A.; Ferretti, A.; Prati, C.; Puglisi, G.; Rocca, F. Calibration of atmospheric effects on SAR interferograms by GPS and local atmosphere models: first results. J. Atmos. Terr. Phys. 2001, 63, 1343–1357. [Google Scholar]
  55. Bevis, M.; Businger, S.; Herring, T.A.; Rocken, R.; Anthes, R.A.; Ware, R.H. GPS Meteorology: Remote sensing of atmospheric water vapor using the Global Positioning System. J. Geophys. Res. 1992, 97, 15787–15801. [Google Scholar]
  56. Liu, Y. Remote Sensing of Water Vapor Content Using GPS Data in Hong Kong Region. Ph.D dissertation, Hong Kong Polytechnic University, Hung Hom, 1999. [Google Scholar]
  57. Bock, Y.; Williams, S. Integrated Satellite Interferometry in Southern California. EOS 1997, 78(29). [Google Scholar]
  58. Janssen, V.; Ge, L.L.; Rizos, C. Tropospheric correction to SAR interferometry from GPS observations. GPS Solut. 2004, 8, 140–151. [Google Scholar]
  59. Webley, P.W.; Bingley, R.M.; Dodson, A.H.; Wadge, G.; Waugh, S.J.; James, I.N. Atmospheric water vapor correction to InSAR surface motion measurements on mountains: results from a dense GPS network on Mount Etna. Phys. Chem. Earth. 2002, 27, 363–370. [Google Scholar]
  60. Li, Z.H.; Fielding, E.J.; Cross, P.; Muller, J.-P. Interferometric synthetic aperture radar atmospheric correction: GPS topography-dependent turbulence model. J. Geophys. Res. 2006, 111, B02404. [Google Scholar] [CrossRef]
  61. Onn, F.; Zebker, H.A. Correction for interferometric synthetic aperture radar atmospheric phase artifacts using time series of zenith wet delay observations from a GOS network. J. Geophys. Res. 2006, 111, B09102. [Google Scholar] [CrossRef]
  62. Onn, F. Modeling water vapor using GPS with application to mitigating InSAR atmospheric distortions. Ph.D dissertation, Stanford University, 2006; p. 176. [Google Scholar]
  63. Wadge, G.; Webley, P.W.; James, I.N.; Bingley, R.; Dodson, A.; Waugh, S.; Veneboer, T.; Puglisi, G.; Mattia, M.; Baker, D.; Edwards, S.C.; Edwards, S.J.; Clarke, P.J. Atmospheric models, GPS and InSAR measurements of the tropospheric water vapor fields over Mount Etna. Geophys. Res. Lett. 2002, 29, 1905–1908. [Google Scholar]
  64. Webley, P.W.; Wadge, G.; James, I.N. Determining radio wave delay by non-hydrostatic atmospheric modeling of water vapour over mountains. Phys. Chem. Earth. 2004, 29, 139–148. [Google Scholar]
  65. Webley, P.W. Atmospheric water vapor correction to InSAR surface motion measurements on mountains: Case Study on Mount Etna. Ph.D dissertation, University of Reading, UK, Reading, 2003. [Google Scholar]
  66. Foster, J.; Brooks, B.; Cherubini, T.; Shacat, C.; Businger, S.; Werner, C. Mitigating atmospheric noise for InSAR using a high resolution weather model. Geophys. Res. Lett. 2006, 33, L16304. [Google Scholar]
  67. Puysségur, B.; Michel, R.; Avouac, J.-P. Tropospheric phase delay in interferometric synthetic aperture radar estimated from meteorological model and multispectral imagery. J. Geophys. Res. 2007, 112, B05419. [Google Scholar]
  68. Gao, B.-C.; Kaufman, Y.J. Water vapor retrievals using Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared channels. J. Geophys. Res. 2003, 108, 4389–4398. [Google Scholar]
  69. Li, Z.H.; Muller, J.P.; Cross, P. Tropospheric correction techniques in repeat-pass SAR interferometry. Proceedings of the FRINGE 2003 workshop, Frascati, Italy, 1- 5, December 2003.
  70. Li, Z.H.; Muller, J.-P.; Cross, P.; Fielding, E.J. Interferometric synthetic aperture radar (InSAR) atmospheric correction: GPS, Moderate Resolution Imaging Spectroradiometer (MODIS), and InSAR integration. J. Geophys. Res. 2005, 110, B03410. [Google Scholar]
  71. Journel, A.G.; Kyriakidis, P.C.; Mao, S. Correcting the Smoothing Effect of Estimators: A Spectral Postprocessor. Math. Geol. 2000, 32, 787–813. [Google Scholar]
  72. Li, Z. W. Modeling atmospheric effects on repeat-pass InSAR measurements. Ph.D Dissertation, The Hong Kong Polytechnic University, Hong Kong, 2005; p. 143. [Google Scholar]
  73. Li, Z. H.; Muller, J. P.; Cross, P.; Albert, P.; Fischer, J.; Bennartz, R. Assessment of the potential of MERIS near-infrared water vapour products to correct ASAR interferometric measurements. Int. J. Remote Sens. 2006, 27, 349–365. [Google Scholar]
  74. Li, Z.H.; Fielding, E.J.; Cross, P.; Muller, J.P. Interferometric synthetic aperture radar atmospheric correction: Medium Resolution Imaging Spectrometer and Advanced Synthetic Aperture Radar integration. Geophys. Res. Lett. 2006, 33, L06816. [Google Scholar]
  75. Sarti, F.; Vadon, H.; Massonnet, D. A method for the automatic characterization of atmospheric artifacts in SAR interferograms by correlation of multiple interferograms over the same site. Proceedings of IGARRS'99, Hamburg, Germany, 28 June- 2 July 1999.
  76. Fruneau, B.; Sarti, F. Detection of ground subsidence in the city of Paris using radar interferometry: isolation of deformation from atmospheric artifact using correlation. Geophys. Res. Lett. 2000, 27, 3981–3984. [Google Scholar]
  77. Sarti, F.; Fruneau, B.; Cunha, T. Isolation of atmospheric artifacts in differential interferometry for ground displacement detection: comparison of different methods. Proceedings of the European Space Agency ERS-Envisat Symposium, Gothenburg, Sweden, 16 - 20 October 2000.
  78. Remy, D.; Bonvalot, S.; Briole, P.; Murakami, M. Accurate measurements of tropospheric effects in volcanic areas from SAR interferometry data: application to Sakurajima volcano (Japan). Earth Planet. Sci. Lett. 2003, 213, 299–310. [Google Scholar]
  79. Ferretti, A.; Prati, C.; Rocca, F. Permanent Scatters in SAR interferometry. IEEE T. Geosci. Remote Sens. 2001, 39, 8–20. [Google Scholar]
  80. Hooper, A.; Zebker, H.; Segall, P.; Kampes, B. A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers. Geophys. Res. Lett. 2004, 31, L23611. [Google Scholar]
  81. Chaabane, F.; Avallone, A.; Tupin, F.; Briole, P.; Maître, H. A Multitemporal Method for Correction of Tropospheric Effects in Differential SAR Interferometry: Application to the Gulf of Corinth Earthquake. IEEE T. Geosci. Remote Sens. 2007, 45, 1605–1615. [Google Scholar]
  82. Ferreti, A.; Novali, F.; Passera, E.; Prati, C.; Rocca, F. Statistical analysis of atmospherical components in ERS SAR data. In Fringe'05; Esrin, Prascati; 28 November - 2 December 2005.
  83. Curlander, J.C.; McDonough, R.N. Synthetic Aperture Radar: Systems and Signal Processing.; Wiley: New York, 1991. [Google Scholar]
Figure 1. Interferometric geometry (from Li et al. [23]).
Figure 1. Interferometric geometry (from Li et al. [23]).
Sensors 08 05426f1
Figure 2. Geomtry of DInSAR.
Figure 2. Geomtry of DInSAR.
Sensors 08 05426f2
Figure 3. Radon transform of atmospheric signals in SAR pairs acquired on: (a) 19 and 20 February, 1996, (b) 25 and 26 March, 1996, (c) 3 and 4 June, 1996, and (d) 16 November and 21 December, 1999 (from Li et al. [33]).
Figure 3. Radon transform of atmospheric signals in SAR pairs acquired on: (a) 19 and 20 February, 1996, (b) 25 and 26 March, 1996, (c) 3 and 4 June, 1996, and (d) 16 November and 21 December, 1999 (from Li et al. [33]).
Sensors 08 05426f3
Figure 4. Power spectra of atmospheric signals for SAR pairs acquired on: (a) 19 and 20 February, 1996, (b) 25 and 26 March, 1996, (c) 3 and 4 June, 1996, and (d) 16 November and 21 December, 1999. (from Li et al. [33]).
Figure 4. Power spectra of atmospheric signals for SAR pairs acquired on: (a) 19 and 20 February, 1996, (b) 25 and 26 March, 1996, (c) 3 and 4 June, 1996, and (d) 16 November and 21 December, 1999. (from Li et al. [33]).
Sensors 08 05426f4
Figure 5. Atmospheric path delay corrections for interferometric pair of 29 July 2000 and 18 August 2001 over Los Angeles basin, South California. (a) original interferogram (deformation field has been modeled with GPS observations and removed from the interferogram); (b) interferogram corrected using MODIS data. (from Li [72]).
Figure 5. Atmospheric path delay corrections for interferometric pair of 29 July 2000 and 18 August 2001 over Los Angeles basin, South California. (a) original interferogram (deformation field has been modeled with GPS observations and removed from the interferogram); (b) interferogram corrected using MODIS data. (from Li [72]).
Sensors 08 05426f5
Figure 6. Atmospheric path delay correction for interferometric pair of 30 April 2006 and 11 March 2007 over Hong Kong area. (a) original interferogram; (b) interferogram corrected using MERIS data.
Figure 6. Atmospheric path delay correction for interferometric pair of 30 April 2006 and 11 March 2007 over Hong Kong area. (a) original interferogram; (b) interferogram corrected using MERIS data.
Sensors 08 05426f6

Share and Cite

MDPI and ACS Style

Ding, X.-l.; Li, Z.-w.; Zhu, J.-j.; Feng, G.-c.; Long, J.-p. Atmospheric Effects on InSAR Measurements and Their Mitigation. Sensors 2008, 8, 5426-5448. https://doi.org/10.3390/s8095426

AMA Style

Ding X-l, Li Z-w, Zhu J-j, Feng G-c, Long J-p. Atmospheric Effects on InSAR Measurements and Their Mitigation. Sensors. 2008; 8(9):5426-5448. https://doi.org/10.3390/s8095426

Chicago/Turabian Style

Ding, Xiao-li, Zhi-wei Li, Jian-jun Zhu, Guang-cai Feng, and Jiang-ping Long. 2008. "Atmospheric Effects on InSAR Measurements and Their Mitigation" Sensors 8, no. 9: 5426-5448. https://doi.org/10.3390/s8095426

APA Style

Ding, X. -l., Li, Z. -w., Zhu, J. -j., Feng, G. -c., & Long, J. -p. (2008). Atmospheric Effects on InSAR Measurements and Their Mitigation. Sensors, 8(9), 5426-5448. https://doi.org/10.3390/s8095426

Article Metrics

Back to TopTop