Interferometric synthetic-aperture radar

From Wikipedia, the free encyclopedia

Interferometric synthetic aperture radar, abbreviated InSAR (or deprecated IfSAR), is a radar technique used in geodesy and remote sensing. This geodetic method uses two or more synthetic aperture radar (SAR) images to generate maps of surface deformation or digital elevation, using differences in the phase of the waves returning to the satellite[1][2][3] or aircraft. The technique can potentially measure millimetre-scale changes in deformation over spans of days to years. It has applications for geophysical monitoring of natural hazards, for example earthquakes, volcanoes and landslides, and in structural engineering, in particular monitoring of subsidence and structural stability.

Interferogram produced using ERS-2 data from 13 August and 17 September 1999, spanning the 17 August Izmit (Turkey) earthquake. (NASA/JPL-Caltech)

Technique[]

SAR amplitude image of Kīlauea (NASA/JPL-Caltech)

Synthetic aperture radar[]

Main article: Synthetic aperture radar

Synthetic aperture radar (SAR) is a form of radar in which sophisticated processing of radar data is used to produce a very narrow effective beam. It can be used to form images of relatively immobile targets; moving targets can be blurred or displaced in the formed images. SAR is a form of active remote sensing – the antenna transmits radiation that is reflected from the image area, as opposed to passive sensing, where the reflection is detected from ambient illumination. SAR image acquisition is therefore independent of natural illumination and images can be taken at night. Radar uses electromagnetic radiation at microwave frequencies; the atmospheric absorption at typical radar wavelengths is very low, meaning observations are not prevented by cloud cover.

Phase[]

Main article: Phase (waves)
Phase difference

SAR makes use of the amplitude and the absolute phase of the return signal data. In contrast, interferometry uses differential phase of the reflected radiation, either from multiple passes along the same trajectory and/or from multiple displaced phase centers (antennas) on a single pass. Since the outgoing wave is produced by the satellite, the phase is known, and can be compared to the phase of the return signal. The phase of the return wave depends on the distance to the ground, since the path length to the ground and back will consist of a number of whole wavelengths plus some fraction of a wavelength. This is observable as a phase difference or phase shift in the returning wave. The total distance to the satellite (i.e., the number of whole wavelengths) is known based on the time that it takes for the energy to make the round trip back to the satellite—but it is the extra fraction of a wavelength that is of particular interest and is measured to great accuracy.

In practice, the phase of the return signal is affected by several factors, which together can make the absolute phase return in any SAR data collection essentially arbitrary, with no correlation from pixel to pixel. To get any useful information from the phase, some of these effects must be isolated and removed. Interferometry uses two images of the same area taken from the same position (or, for topographic applications, slightly different positions) and finds the difference in phase between them, producing an image known as an interferogram. This is measured in radians of phase difference and, because of the cyclic nature of phase, is recorded as repeating fringes that each represent a full 2π cycle.

Factors affecting phase[]

Further information: Reflection (physics) and topography

The most important factor affecting the phase is the interaction with the ground surface. The phase of the wave may change on reflection, depending on the properties of the material. The reflected signal back from any one pixel is the summed contribution to the phase from many smaller 'targets' in that ground area, each with different dielectric properties and distances from the satellite, meaning the returned signal is arbitrary and completely uncorrelated with that from adjacent pixels. Importantly though, it is consistent – provided nothing on the ground changes the contributions from each target should sum identically each time, and hence be removed from the interferogram.

Once the ground effects have been removed, the major signal present in the interferogram is a contribution from orbital effects. For interferometry to work, the satellites must be as close as possible to the same spatial position when the images are acquired. This means that images from two satellite platforms with different orbits cannot be compared, and for a given satellite data from the same orbital track must be used. In practice the perpendicular distance between them, known as the baseline, is often known to within a few centimetres but can only be controlled on a scale of tens to hundreds of metres. This slight difference causes a regular difference in phase that changes smoothly across the interferogram and can be modelled and removed.

Corresponding interferogram of Kīlauea, showing topographic fringes (NASA/JPL-Caltech)

The slight difference in satellite position also alters the distortion caused by topography, meaning an extra phase difference is introduced by a stereoscopic effect. The longer the baseline, the smaller the topographic height needed to produce a fringe of phase change – known as the altitude of ambiguity. This effect can be exploited to calculate the topographic height, and used to produce a digital elevation model (DEM).

If the height of the topography is already known, the topographic phase contribution can be calculated and removed. This has traditionally been done in two ways. In the two-pass method, elevation data from an externally derived DEM is used in conjunction with the orbital information to calculate the phase contribution. In the three-pass method two images acquired a short time apart are used to create an interferogram, which is assumed to have no deformation signal and therefore represent the topographic contribution. This interferogram is then subtracted from a third image with a longer time separation to give the residual phase due to deformation.

Once the ground, orbital and topographic contributions have been removed the interferogram contains the deformation signal, along with any remaining noise (see Difficulties below). The signal measured in the interferogram represents the change in phase caused by an increase or decrease in distance from the ground pixel to the satellite, therefore only the component of the ground motion parallel to the satellite line of sight vector will cause a phase difference to be observed. For sensors like ERS with a small incidence angle this measures vertical motion well, but is insensitive to horizontal motion perpendicular to the line of sight (approximately north-south). It also means that vertical motion and components of horizontal motion parallel to the plane of the line of sight (approximately east-west) cannot be separately resolved.

One fringe of phase difference is generated by a ground motion of half the radar wavelength, since this corresponds to a whole wavelength increase in the two-way travel distance. Phase shifts are only resolvable relative to other points in the interferogram. Absolute deformation can be inferred by assuming one area in the interferogram (for example a point away from expected deformation sources) experienced no deformation, or by using a ground control (GPS or similar) to establish the absolute movement of a point.

Difficulties[]

A variety of factors govern the choice of images which can be used for interferometry. The simplest is data availability – radar instruments used for interferometry commonly don't operate continuously, acquiring data only when programmed to do so. For future requirements it may be possible to request acquisition of data, but for many areas of the world archived data may be sparse. Data availability is further constrained by baseline criteria. Availability of a suitable DEM may also be a factor for two-pass InSAR; commonly 90 m SRTM data may be available for many areas, but at high latitudes or in areas of poor coverage alternative datasets must be found.

A fundamental requirement of the removal of the ground signal is that the sum of phase contributions from the individual targets within the pixel remains constant between the two images and is completely removed. However, there are several factors that can cause this criterion to fail. Firstly the two images must be accurately co-registered to a sub-pixel level to ensure that the same ground targets are contributing to that pixel. There is also a geometric constraint on the maximum length of the baseline – the difference in viewing angles must not cause phase to change over the width of one pixel by more than a wavelength. The effects of topography also influence the condition, and baselines need to be shorter if terrain gradients are high. Where co-registration is poor or the maximum baseline is exceeded the pixel phase will become incoherent – the phase becomes essentially random from pixel to pixel rather than varying smoothly, and the area appears noisy. This is also true for anything else that changes the contributions to the phase within each pixel, for example changes to the ground targets in each pixel caused by vegetation growth, landslides, agriculture or snow cover.

Another source of error present in most interferograms is caused by the propagation of the waves through the atmosphere. If the wave travelled through a vacuum it should theoretically be possible (subject to sufficient accuracy of timing) to use the two-way travel-time of the wave in combination with the phase to calculate the exact distance to the ground. However, the velocity of the wave through the atmosphere is lower than the speed of light in a vacuum, and depends on air temperature, pressure and the partial pressure of water vapour.[4] It is this unknown phase delay that prevents the integer number of wavelengths being calculated. If the atmosphere was horizontally homogeneous over the length scale of an interferogram and vertically over that of the topography then the effect would simply be a constant phase difference between the two images which, since phase difference is measured relative to other points in the interferogram, would not contribute to the signal. However, the atmosphere is laterally heterogeneous on length scales both larger and smaller than typical deformation signals. This spurious signal can appear completely unrelated to the surface features of the image, however, in other cases the atmospheric phase delay is caused by vertical inhomogeneity at low altitudes and this may result in fringes appearing to correspond with the topography.

Persistent scatterer InSAR[]

Further information: Scattering

Persistent or permanent scatterer techniques are a relatively recent development from conventional InSAR, and rely on studying pixels which remain coherent over a sequence of interferograms. In 1999, researchers at Politecnico di Milano, Italy, developed a new multi-image approach in which one searches the stack of images for objects on the ground providing consistent and stable radar reflections back to the satellite. These objects could be the size of a pixel or, more commonly, sub-pixel sized, and are present in every image in the stack. That specific implementation is patented.

Some research centres and companies, were inspired to develop variations of their own algorithms which would also overcome InSAR's limitations. In scientific literature, these techniques are collectively referred to as persistent scatterer interferometry or PSI techniques. The term persistent scatterer interferometry (PSI) was proposed by European Space Agency (ESA) to define the second generation of radar interferometry techniques. This term is nowadays commonly accepted by scientific and the end user community.

Commonly such techniques are most useful in urban areas with many permanent structures, for example the PSI studies of European geohazard sites undertaken by the Terrafirma project.[5] The Terrafirma project provides a ground motion hazard information service, distributed throughout Europe via national geological surveys and institutions. The objective of this service is to help save lives, improve safety, and reduce economic loss through the use of state-of-the-art PSI information. Over the last 9 years this service has supplied information relating to urban subsidence and uplift, slope stability and landslides, seismic and volcanic deformation, coastlines and flood plains.

Producing interferograms[]

Further information: Interferogram

The processing chain used to produce interferograms varies according to the software used and the precise application but will usually include some combination of the following steps.

Two SAR images are required to produce an interferogram; these may be obtained pre-processed, or produced from raw data by the user prior to InSAR processing. The two images must first be co-registered, using a correlation procedure to find the offset and difference in geometry between the two amplitude images. One SAR image is then re-sampled to match the geometry of the other, meaning each pixel represents the same ground area in both images. The interferogram is then formed by cross-multiplication of each pixel in the two images, and the interferometric phase due to the curvature of the Earth is removed, a process referred to as flattening. For deformation applications a DEM can be used in conjunction with the baseline data to simulate the contribution of the topography to the interferometric phase, this can then be removed from the interferogram.

Once the basic interferogram has been produced, it is commonly filtered using an adaptive power-spectrum filter to amplify the phase signal. For most quantitative applications the consecutive fringes present in the interferogram will then have to be unwrapped, which involves interpolating over the 0 to 2π phase jumps to produce a continuous deformation field. At some point, before or after unwrapping, incoherent areas of the image may be masked out. The final processing stage involves geocoding the image, which resamples the interferogram from the acquisition geometry (related to direction of satellite path) into the desired geographic projection.

Hardware[]

Seasat (NASA/JPL-Caltech)

Spaceborne[]

Early exploitation of satellite-based InSAR included use of Seasat data in the 1980s, but the potential of the technique was expanded in the 1990s, with the launch of ERS-1 (1991), JERS-1 (1992), RADARSAT-1 and ERS-2 (1995). These platforms provided the stable, well-defined orbits and short baselines necessary for InSAR. More recently, the 11-day NASA STS-99 mission in February 2000 used a SAR antenna mounted on the space shuttle to gather data for the Shuttle Radar Topography Mission. In 2002 ESA launched the ASAR instrument, designed as a successor to ERS, aboard Envisat. While the majority of InSAR to date has utilised the C-band sensors, recent missions such as the ALOS PALSAR, TerraSAR-X and COSMO-SkyMed are expanding the available data in the L- and X-band.

Most recently, ESA launched Sentinel-1A and Sentinel-1B – two C-band sensors. Together, they provide InSAR coverage on a global scale and on a 6-day repeat cycle.

Airborne[]

Airborne InSAR data acquisition systems are built by companies such as the American Intermap, the German , and the Brazilian .[6]

Terrestrial or ground-based[]

A deformation plot showing slope instability using Terrestrial InSAR

Terrestrial or ground-based SAR interferometry (GBInSAR or TInSAR) is a remote sensing technique for the displacement monitoring of slopes,[7] rock scarps, volcanoes, landslides, buildings, infrastructures etc. This technique is based on the same operational principles of the satellite SAR interferometry, but the synthetic aperture of the radar (SAR) is obtained by an antenna moving on a rail instead of a satellite moving around an orbit. SAR technique allows 2D radar image of the investigated scenario to be achieved, with a high range resolution (along the instrumental line of sight) and cross-range resolution (along the scan direction). Two antennas respectively emit and receive microwave signals and, by calculating the phase difference between two measurements taken in two different times, it is possible to compute the displacement of all the pixels of the SAR image. The accuracy in the displacement measurement is of the same order of magnitude as the EM wavelength and depends also on the specific local and atmospheric conditions.

Applications[]

Rapid ground subsidence over the Lost Hills oil field in California. (NASA/JPL-Caltech)

Tectonic[]

InSAR can be used to measure tectonic deformation, for example ground movements due to earthquakes. It was first used for the 1992 Landers earthquake,[8] but has since been utilised extensively for a wide variety of earthquakes all over the world. In particular the 1999 Izmit and 2003 Bam earthquakes were extensively studied.[9][10] InSAR can also be used to monitor creep and strain accumulation on faults.

Volcanic[]

InSAR can be used in a variety of volcanic settings, including deformation associated with eruptions, inter-eruption strain caused by changes in magma distribution at depth, gravitational spreading of volcanic edifices, and volcano-tectonic deformation signals.[11] Early work on volcanic InSAR included studies on Mount Etna,[12] and Kilauea,[13] with many more volcanoes being studied as the field developed. The technique is now widely used for academic research into volcanic deformation, although its use as an operational monitoring technique for volcano observatories has been limited by issues such as orbital repeat times, lack of archived data, coherence and atmospheric errors.[14][15] Recently InSAR has been used to study rifting processes in Ethiopia.[16]

Subsidence[]

Ground subsidence from a variety of causes has been successfully measured using InSAR, in particular subsidence caused by oil or water extraction from underground reservoirs,[17] subsurface mining and collapse of old mines.[18] Thus, InSAR has become an indispensable tool to satisfactorily address many subsidence studies. Tomás et al.[19] performed a cost analysis that allowed to identify the strongest points of InSAR techniques compared with other conventional techniques: (1) higher data acquisition frequency and spatial coverage; and (2) lower annual cost per measurement point and per square kilometre.

Landslides[]

Although InSAR technique can present some limitations when applied to landslides,[20] it can also be used for monitoring landscape features such as landslides.[21][22][23]

Ice flow[]

Glacial motion and deformation have been successfully measured using satellite interferometry. The technique allows remote, high-resolution measurement of changes in glacial structure, ice flow, and shifts in ice dynamics, all of which agree closely with ground observations.[24]

Kamchatka Peninsula, Landsat data draped over SRTM digital elevation model (NASA/JPL-Caltech)

Infrastructure and building monitoring[]

InSAR can also be used to monitor the stability of built structures.[25] Very high resolution SAR data (such as derived from the TerraSAR-X StripMap mode or COSMO-Skymed HIMAGE mode) are especially suitable for this task. InSAR is used for monitoring highway and railway settlements,[26][27] dike stability,[28] forensic engineering [29] and many other uses.

DEM generation[]

Interferograms can be used to produce digital elevation maps (DEMs) using the stereoscopic effect caused by slight differences in observation position between the two images. When using two images produced by the same sensor with a separation in time, it must be assumed other phase contributions (for example from deformation or atmospheric effects) are minimal. In 1995 the two ERS satellites flew in tandem with a one-day separation for this purpose. A second approach is to use two antennas mounted some distance apart on the same platform, and acquire the images at the same time, which ensures no atmospheric or deformation signals are present. This approach was followed by NASA's SRTM mission aboard the space shuttle in 2000. InSAR-derived DEMs can be used for later two-pass deformation studies, or for use in other geophysical applications.

See also[]

Further reading[]

  • B. Kampes, Radar Interferometry – Persistent Scatterer Technique, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2006. ISBN 978-1-4020-4576-9

References[]

  1. ^ Massonnet, D.; Feigl, K. L. (1998), "Radar interferometry and its application to changes in the earth's surface", Rev. Geophys., vol. 36, no. 4, pp. 441–500, Bibcode:1998RvGeo..36..441M, doi:10.1029/97RG03139
  2. ^ Burgmann, R.; Rosen, P.A.; Fielding, E.J. (2000), "Synthetic aperture radar interferometry to measure Earth's surface topography and its deformation", Annual Review of Earth and Planetary Sciences, vol. 28, pp. 169–209, Bibcode:2000AREPS..28..169B, doi:10.1146/annurev.earth.28.1.169
  3. ^ Hanssen, Ramon F. (2001), Radar Interferometry: Data Interpretation and Error Analysis, Kluwer Academic, ISBN 9780792369455
  4. ^ Zebker, H.A.; Rosen, P.A.; Hensley, S. (1997), "Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps", Journal of Geophysical Research, vol. 102, no. B4, pp. 7547–7563, Bibcode:1997JGR...102.7547Z, doi:10.1029/96JB03804
  5. ^ "Terrafirma.eu.com: A pan-European ground hazard information service". Retrieved 22 January 2013.
  6. ^ "Revista Pesquisa Fapesp".
  7. ^ Longstaff, I.D. (2011). "Comparing real beam and synthetic aperture techniques for Slope Stability Radar" (PDF). Whitepaper, University of Queensland, Australia.[permanent dead link]
  8. ^ Massonnet, D.; Rossi, M.; Carmona, C.; Adragna, F.; Peltzer, G.; Feigl, K.; Rabaute, T. (1993), "The displacement field of the Landers earthquake mapped by radar interferometry", Nature, vol. 364, no. 6433, pp. 138–142, Bibcode:1993Natur.364..138M, doi:10.1038/364138a0, S2CID 4355142
  9. ^ "Envisat's rainbow vision detects ground moving at pace fingernails grow". European Space Agency. August 6, 2004. Retrieved 2007-03-22.
  10. ^ "The Izmit Earthquake of 17 August 1999 in Turkey". European Space Agency. Retrieved 2007-03-22.
  11. ^ Wadge, G. (2003), "A strategy for the observation of volcanism on Earth from space", Phil. Trans. Royal Soc.Lond., vol. 361, no. 1802, pp. 145–156, Bibcode:2003RSPTA.361..145W, doi:10.1098/rsta.2002.1117, S2CID 25985116
  12. ^ Massonnet, D.; Briole, P.; Arnaud, A. (1995), "Deflation of Mount Etna monitored by spaceborne radar interferometry", Nature, vol. 375, no. 6532, pp. 567–570, Bibcode:1995Natur.375..567M, doi:10.1038/375567a0, S2CID 4281294
  13. ^ Rosen, P. A.; Hensley, S.; Zebker, H. A.; Webb, F. H.; Fielding, E. J. (1996), "Surface deformation and coherence measurements of Kilauea Volcano, Hawaii, from SIR C radar interferometry", J. Geophys. Res., vol. 101, no. E10, pp. 23, 109–23, 126, Bibcode:1996JGR...10123109R, doi:10.1029/96JE01459
  14. ^ Pritchard, Matthew E.; Simons, Mark (2004). "Surveying Volcanic Arcs with Satellite Radar Interferometry: The Central Andes, Kamchatka, and Beyond". GSA Today. 14 (8): 4. doi:10.1130/1052-5173(2004)014<4:svawsr>2.0.co;2.
  15. ^ Stevens, N.F.; Wadge, G. (2004), "Towards operational repeat-pass SAR interferometry at active volcanoes", Natural Hazards, vol. 33, pp. 47–76, doi:10.1023/B:NHAZ.0000035005.45346.2b, S2CID 129026183
  16. ^ Wright, T.J.; Ebinger, C.; Biggs, J.; Ayele, A.; Yirgu, G.; Keir, D.; Stork, A. (2006), "Magma-maintained rift segmentation at continental rupture in the 2005 Afar dyking episode" (PDF), Nature, vol. 442, no. 7100, pp. 291–294, Bibcode:2006Natur.442..291W, doi:10.1038/nature04978, hdl:2158/1078052, PMID 16855588, S2CID 4319443
  17. ^ Tomás, R.; Márquez, Y.; Lopez-Sanchez, J.M.; Delgado, J.; Blanco, P.; Mallorquí, J.J.; Martínez, M.; Herrera, G.; Mulas, J. (2005). "Mapping ground subsidence induced by aquifer overexploitation using advanced Differential SAR Interferometry: Vega Media of the Segura River (SE Spain) case study". Remote Sensing of Environment. 98 (2–3): 269–283. Bibcode:2005RSEnv..98..269T. doi:10.1016/j.rse.2005.08.003. hdl:2117/13208.
  18. ^ Herrera, G.; Tomás, R.; Lopez-Sanchez, J.M.; Delgado, J.; Mallorqui, J.J.; Duque, S.; Mulas, J. (2007). "Advanced DInSAR analysis on mining areas: La Union case study (Murcia, SE Spain)". Engineering Geology. 90 (3–4): 148–159. doi:10.1016/j.enggeo.2007.01.001. hdl:2117/12906.
  19. ^ Tomás, R.; Romero, R.; Mulas, J.; Marturià, J.J.; Mallorquí, J.J.; Lopez-Sanchez, J.M.; Herrera, G.; Gutiérrez, F.; González, P.J.; Fernández, J.; Duque, S.; Concha-Dimas, A.; Cocksley, G.; Castañeda, C.; Carrasco, D.; Blanco, P. (2014). "Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain". Environmental Earth Sciences. 71: 163–181. doi:10.1007/s12665-013-2422-z. hdl:10045/36419. S2CID 128740704.
  20. ^ Colesanti, C.; Wasowski, J. (2006). "Investigating landslides with space-borne Synthetic Aperture Radar (SAR) interferometry". Engineering Geology. 88 (3–4): 173–199. doi:10.1016/j.enggeo.2006.09.013.
  21. ^ "Ground motion". European Space Agency. Archived from the original on 2008-05-21. Retrieved 2007-03-21.
  22. ^ Herrera, G.; Tomás, R.; Vicente, F.; Lopez-Sanchez, J.M.; Mallorquí, J.J.; Mulas, J. (October 2010). "Mapping ground movements in open pit mining areas using differential SAR interferometry". International Journal of Rock Mechanics and Mining Sciences. 47 (7): 1114–1125. doi:10.1016/j.ijrmms.2010.07.006.
  23. ^ Tomás, R.; Li, Z.; Liu, P.; Singleton, A.; Hoey, T.; Cheng, X. (2014). "Spatiotemporal characteristics of the Huangtupo landslide in the Three Gorges region (China) constrained by radar interferometry". Geophysical Journal International. 197 (1): 213–232. Bibcode:2014GeoJI.197..213T. doi:10.1093/gji/ggu017.
  24. ^ Goldstein, R.M.; Engelhardt, H.; Kamb, B.; Frolich, R. M. (1993), "Satellite radar interferometry for monitoring ice sheet motion: application to an Antarctic ice streamy", Science, vol. 262, no. 5139, pp. 1525–1530, Bibcode:1993Sci...262.1525G, doi:10.1126/science.262.5139.1525, PMID 17829380, S2CID 42622639
  25. ^ Tomás, R.; García-Barba, J.; Cano, M.; Sanabria, M.P.; Ivorra, S.; Duro, J.; Herrera, G. (November 2012). "Subsidence damage assessment of a gothic church using Differential Interferometry and field data". Structural Health Monitoring. 11 (6): 751–762. doi:10.1177/1475921712451953. hdl:10045/55037. S2CID 112142102.
  26. ^ Yu, B.; Liu, G.; Zhang, R.; Jia, H.; Li, T.; Wang, X.; Dai, K.; Ma, D. (2013). "Monitoring subsidence rates along road network by persistent scatterer SAR interferometry with high-resolution TerraSAR-X imagery". Journal of Modern Transportation. 21 (4): 236–246. doi:10.1007/s40534-013-0030-y.
  27. ^ Bianchini Ciampoli, L.; Gagliardi, V.; Clementini, C.; Latini, D.; Del Frate, F.; Benedetto, A. (2020). "Transport Infrastructure Monitoring by InSAR and GPR Data Fusion". Surveys in Geophysics. 41 (3): 371–394. doi:10.1007/s10712-019-09563-7.
  28. ^ Tomás, R.; Cano, M.; García-Barba, J.; Vicente, F.; Herrera, G.; Lopez-Sanchez, J.M.; Mallorquí, J.J. (2013). "Monitoring an earthfill dam using differential SAR interferometry: La Pedrera dam, Alicante, Spain". Engineering Geology. 157: 21–32. doi:10.1016/j.enggeo.2013.01.022. hdl:2117/19542.
  29. ^ Herrera, G.; Tomás, R.; Monells, D.; Centolanza, G.; Mallorquí, J.J.; Vicente, F.; Navarro, V.D.; Lopez-Sanchez, J.M.; Sanabria, M.; Cano, M.; Mulas, J. (2010). "Analysis of subsidence using TerraSAR-X data: Murcia case study". Engineering Geology. 116 (3–4): 284–295. doi:10.1016/j.enggeo.2010.09.010.

External links[]

Retrieved from ""