Polarimetric radar interferometry for improved mine detection and surface clutter rejection

A recently developed technique, polarimetric radar interferometry, is applied to tackle the problem of the detection of buried objects embedded in surface clutter. An experiment with a fully polarimetric radar in an anechoic chamber has been carried out using different frequency bands and baselines. The processed results show the ability of this technique to detect buried plastic mines and to measure their depth. This technique enables the detection of plastic mines even if their backscatter response is much lower than that of the surface clutter.

Ground deformation at Soufriere Hills Volcano, Montserrat during 1998-2000 measured by radar interferometry and GPS

We examine the motion of the ground surface on the Soufriere Hills Volcano, Montserrat between 1998 and 2000 using radar interferometry (InSAR). To minimise the effects of variable atmospheric water vapour on the InSAR measurements we use independently-derived measurements of the radar path delay from six continuous GPS receivers. The surfaces providing a measurable inter-ferometric signal are those on pyroclastic flow deposits, mainly emplaced in 1997. Three types of surface motion can be discriminated. Firstly, the surfaces of thick, valley-filling deposits subsided at rates of 150-120 mm/year in the year after emplacement to 50-30 mm/year two years later. This must be due to contraction and settling effects during cooling. The second type is the near-field motion localised within about one kilometre of the dome. Both subsidence and uplift events are seen and though the former could be due to surface gravitational effects, the latter may reflect shallow (< 1 km) pressurisation effects within the conduit/dome. Far-field motions of the surface away from the deeply buried valleys are interpreted as crustal strains. Because the flux of magma to the surface stopped from March 1998 to November 1999 and then resumed from November 1999 through 2000, we use InSAR data from these two periods to test the crustal strain behaviour of three models of magma supply: open, depleting and unbalanced. The InSAR observations of strain gradients of 75-80 mm/year/krn uplift during the period of quiescence on the western side of the volcano are consistent with an unbalanced model in which magma supply into a crustal magma chamber continues during quiescence, raising chamber pressure that is then released upon resumption of effusion. GPS motion vectors agree qualitatively with the InSAR displacements but are of smaller magnitude. The discrepancy may be due to inaccurate compensation for atmospheric delays in the InSAR data. (c) 2005 Elsevier B.V. All rights reserved.

Correction of atmospheric delay effects in radar interferometry using a nested mesoscale atmospheric model

The temporal variability of the atmosphere through which radio waves pass in the technique of differential radar interferometry can seriously limit the accuracy with which the method can measure surface motion. A forward, nested mesoscale model of the atmosphere can be used to simulate the variable water content along the radar path and the resultant phase delays. Using this approach we demonstrate how to correct an interferogram of Mount Etna in Sicily associated with an eruption in 2004-5. The regional mesoscale model (Unified Model) used to simulate the atmosphere at higher resolutions consists of four nested domains increasing in resolution (12, 4, 1, 0.3 km), sitting within the analysis version of a global numerical model that is used to initiate the simulation. Using the high resolution 3D model output we compute the surface pressure, temperature and the water vapour, liquid and solid water contents, enabling the dominant hydrostatic and wet delays to be calculated at specific times corresponding to the acquisition of the radar data. We can also simulate the second-order delay effects due to liquid water and ice.

Co-eruptive and inter-eruptive surface deformation measured by satellite radar interferometry at Nyamuragira volcano, D.R. Congo, 1996 to 2010.

Interferometric Synthetic Aperture Radar (InSAR) measurements of surface deformation at Nyamuragira Volcano between 1996 and 2010 reveal a variety of co-eruptive and inter-eruptive signals. During 7 of the 8 eruptions in this period deformation was measured that is consistent with the emplacement of shallow near-vertical dykes feeding the eruptive fissures and associated with a NNW-trending fissure zone that traverses the summit caldera. Between eruptions the caldera and the summit part of this fissure zone subsided gradually (b3–5 cm/year). We also find evidence of post-eruption subsidence around the sites of the main vents of some flank eruptions (2002, 2004, 2006, and 2010). In the 6 months prior to the 2010 eruption a10-km wide zone centred on the caldera inflated by 1–2 cm. The low magnitude of this signal suggests that the presumed magma reservoir at 3–8 km depth contains highly compressible magma with little stored elastic strain energy. To the north of the caldera the fissure zone splits into WNW and NE branches around a zone that has a distinct InSAR signal. We interpret this zone to represent an elevated, 'stable' block of basement rocks buried by lavas within the Rift Zone.

Registration And Visualisation Of Deformation Maps From Terrestrial Radar Interferometry Using Photogrammetry And Structure From Motion

This paper describes a general workflow for the registration of terrestrial radar interferometric data with 3D point clouds derived from terrestrial photogrammetry and structure from motion. After the determination of intrinsic and extrinsic orientation parameters, data obtained by terrestrial radar interferometry were projected on point clouds and then on the initial photographs. Visualisation of slope deformation measurements on photographs provides an easily understandable and distributable information product, especially of inaccessible target areas such as steep rock walls or in rockfall run-out zones. The suitability and error propagation of the referencing steps and final visualisation of four approaches are compared: (a) the classic approach using a metric camera and stereo-image photogrammetry; (b) images acquired with a metric camera, automatically processed using structure from motion; (c) images acquired with a digital compact camera, processed with structure from motion; and (d) a markerless approach, using images acquired with a digital compact camera using structure from motion without artificial ground control points. The usability of the completely markerless approach for the visualisation of high-resolution radar interferometry assists the production of visualisation products for interpretation.

Vehicular traffic surveillance and road lane detection using radar interferometry

Speed enforcement on public roadways is an important issue in order to guarantee road security and to reduce the number and seriousness of traffic accidents. Traditionally, this task has been partially solved using radar and/or laser technologies and, more recently, using video-camera based systems. All these systems have significant shortcomings that have yet to be overcome. The main drawback of classical Doppler radar technology is that the velocity measurement fails when several vehicles are in the radars beam. Modern radar systems are able to measure speed and range between vehicle and radar. However, this is not enough to discriminate the lane where the vehicle is driving on. The limitation of several vehicles in the beam is overcome using laser technology. However, laser systems have another important limitation: They cannot measure the speed of several vehicles simultaneously. Novel video-camera systems, based on license plate identification, solve the previous drawbacks, but they have the problem that they can only measure average speed but never top-speed. This paper studies the feasibility of using an interferometric linear frequency modulated continuous wave radar to improve top-speed enforcement on roadways. Two different systems based on down-the-road and across-the-road radar configurations are presented. The main advantage of the proposed solutions is they can simultaneously measure speed, range, and lane of several vehicles, allowing the univocal identification of the offenders. A detailed analysis about the operation and accuracy of these solutions is reported. In addition, the feasibility of the proposed techniques has been demonstrated with simulations and real experiments using a Ka-band interferometric radar developed by our research group.

Spatiotemporal characteristics of the Huangtupo landslide in the Three Gorges region (China) constrained by radar interferometry

The Huangtupo landslide is one of the largest in the Three Gorges region, China. The county-seat town of Badong, located on the south shore between the Xiling and Wu gorges of the Yangtze River, was moved to this unstable slope prior to the construction of the Three Gorges Project, since the new Three Gorges reservoir completely submerged the location of the old city. The instability of the slope is affecting the new town by causing residential safety problems. The Huangtupo landslide provides scientists an opportunity to understand landslide response to fluctuating river water level and heavy rainfall episodes, which is essential to decide upon appropriate remediation measures. Interferometric Synthetic Aperture Radar (InSAR) techniques provide a very useful tool for the study of superficial and spatially variable displacement phenomena. In this paper, three sets of radar data have been processed to investigate the Huangtupo landslide. Results show that maximum displacements are affecting the northwest zone of the slope corresponding to Riverside slumping mass I#. The other main landslide bodies (i.e. Riverside slumping mass II#, Substation landslide and Garden Spot landslide) exhibit a stable behaviour in agreement with in situ data, although some active areas have been recognized in the foot of the Substation landslide and Garden Spot landslide. InSAR has allowed us to study the kinematic behaviour of the landslide and to identify its active boundaries. Furthermore, the analysis of the InSAR displacement time-series has helped recognize the different displacement patterns on the slope and their relationships with various triggering factors. For those persistent scatterers, which exhibit long-term displacements, they can be decomposed into a creep model (controlled by geological conditions) and a superimposed recoverable term (dependent on external factors), which appears closely correlated with reservoir water level changes close to the river's edge. These results, combined with in situ data, provide a comprehensive analysis of the Huangtupo landslide, which is essential for its management.

Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain

Subsidence related to multiple natural and human-induced processes affects an increasing number of areas worldwide. Although this phenomenon may involve surface deformation with 3D displacement components, negative vertical movement, either progressive or episodic, tends to dominate. Over the last decades, differential SAR interferometry (DInSAR) has become a very useful remote sensing tool for accurately measuring the spatial and temporal evolution of surface displacements over broad areas. This work discusses the main advantages and limitations of addressing active subsidence phenomena by means of DInSAR techniques from an end-user point of view. Special attention is paid to the spatial and temporal resolution, the precision of the measurements, and the usefulness of the data. The presented analysis is focused on DInSAR results exploitation of various ground subsidence phenomena (groundwater withdrawal, soil compaction, mining subsidence, evaporite dissolution subsidence, and volcanic deformation) with different displacement patterns in a selection of subsidence areas in Spain. Finally, a cost comparative study is performed for the different techniques applied.

Imaging Land Subsidence Induced by Groundwater Extraction in Beijing (China) Using Satellite Radar Interferometry

Beijing is one of the most water-stressed cities in the world. Due to over-exploitation of groundwater, the Beijing region has been suffering from land subsidence since 1935. In this study, the Small Baseline InSAR technique has been employed to process Envisat ASAR images acquired between 2003 and 2010 and TerraSAR-X stripmap images collected from 2010 to 2011 to investigate land subsidence in the Beijing region. The maximum subsidence is seen in the eastern part of Beijing with a rate greater than 100 mm/year. Comparisons between InSAR and GPS derived subsidence rates show an RMS difference of 2.94 mm/year with a mean of 2.41 ± 1.84 mm/year. In addition, a high correlation was observed between InSAR subsidence rate maps derived from two different datasets (i.e., Envisat and TerraSAR-X). These demonstrate once again that InSAR is a powerful tool for monitoring land subsidence. InSAR derived subsidence rate maps have allowed for a comprehensive spatio-temporal analysis to identify the main triggering factors of land subsidence. Some interesting relationships in terms of land subsidence were found with groundwater level, active faults, accumulated soft soil thickness and different aquifer types. Furthermore, a relationship with the distances to pumping wells was also recognized in this work.

Forensic analysis of buildings affected by mining subsidence based on Differential Interferometry (Part III)

This work presents a forensic analysis of buildings affected by mining subsidence, which is based on deformation data obtained by Differential Interferometry (DInSAR). The proposed test site is La Union village (Murcia, SE Spain) where subsidence was triggered in an industrial area due to the collapse of abandoned underground mining labours occurred in 1998. In the first part of this work the study area was introduced, describing the spatial and temporal evolution of ground subsidence, through the elaboration of a cracks map on the buildings located within the affected area. In the second part, the evolution of the most significant cracks found in the most damaged buildings was monitored using biaxial extensometric units and inclinometers. This article describes the work performed in the third part, where DInSAR processing of satellite radar data, available between 1998 and 2008, has permitted to determine the spatial and temporal evolution of the deformation of all the buildings of the study area in a period when no continuous in situ instrumental data is available. Additionally, the comparison of these results with the forensic data gathered in the 2005–2008 period, reveal that there is a coincidence between damaged buildings, buildings where extensometers register significant movements of cracks, and buildings deformation estimated from radar data. As a result, it has been demonstrated that the integration of DInSAR data into forensic analysis methodologies contributes to improve significantly the assessment of the damages of buildings affected by mining subsidence.

Subsidence damage assessment of a Gothic church using differential interferometry and field data

The Santas Justa and Rufina Gothic church (fourteenth century) has suffered several physical, mechanical, chemical, and biochemical types of pathologies along its history: rock alveolization, efflorescence, biological activity, and capillary ascent of groundwater. However, during the last two decades, a new phenomenon has seriously affected the church: ground subsidence caused by aquifer overexploitation. Subsidence is a process that affects the whole Vega Baja of the Segura River basin and consists of gradual sinking in the ground surface caused by soil consolidation due to a pore pressure decrease. This phenomenon has been studied by differential synthetic aperture radar interferometry techniques, which illustrate settlements up to 100 mm for the 1993–2009 period for the whole Orihuela city. Although no differential synthetic aperture radar interferometry information is available for the church due to the loss of interferometric coherence, the spatial analysis of nearby deformation combined with fieldwork has advanced the current understanding on the mechanisms that affect the Santas Justa and Rufina church. These results show the potential interest and the limitations of using this remote sensing technique as a complementary tool for the forensic analysis of building structures.

Retrieval of grassland biophysical parameters using multitemporal optical and radar satellite data

The amount and quality of available biomass is a key factor for the sustainable livestock industry and agricultural management related decision making. Globally 31.5% of land cover is grassland while 80% of Ireland’s agricultural land is grassland. In Ireland, grasslands are intensively managed and provide the cheapest feed source for animals. This dissertation presents a detailed state of the art review of satellite remote sensing of grasslands, and the potential application of optical (Moderate–resolution Imaging Spectroradiometer (MODIS)) and radar (TerraSAR-X) time series imagery to estimate the grassland biomass at two study sites (Moorepark and Grange) in the Republic of Ireland using both statistical and state of the art machine learning algorithms. High quality weather data available from the on-site weather station was also used to calculate the Growing Degree Days (GDD) for Grange to determine the impact of ancillary data on biomass estimation. In situ and satellite data covering 12 years for the Moorepark and 6 years for the Grange study sites were used to predict grassland biomass using multiple linear regression, Neuro Fuzzy Inference Systems (ANFIS) models. The results demonstrate that a dense (8-day composite) MODIS image time series, along with high quality in situ data, can be used to retrieve grassland biomass with high performance (R2 = 0:86; p < 0:05, RMSE = 11.07 for Moorepark). The model for Grange was modified to evaluate the synergistic use of vegetation indices derived from remote sensing time series and accumulated GDD information. As GDD is strongly linked to the plant development, or phonological stage, an improvement in biomass estimation would be expected. It was observed that using the ANFIS model the biomass estimation accuracy increased from R2 = 0:76 (p < 0:05) to R2 = 0:81 (p < 0:05) and the root mean square error was reduced by 2.72%. The work on the application of optical remote sensing was further developed using a TerraSAR-X Staring Spotlight mode time series over the Moorepark study site to explore the extent to which very high resolution Synthetic Aperture Radar (SAR) data of interferometrically coherent paddocks can be exploited to retrieve grassland biophysical parameters. After filtering out the non-coherent plots it is demonstrated that interferometric coherence can be used to retrieve grassland biophysical parameters (i. e., height, biomass), and that it is possible to detect changes due to the grass growth, and grazing and mowing events, when the temporal baseline is short (11 days). However, it not possible to automatically uniquely identify the cause of these changes based only on the SAR backscatter and coherence, due to the ambiguity caused by tall grass laid down due to the wind. Overall, the work presented in this dissertation has demonstrated the potential of dense remote sensing and weather data time series to predict grassland biomass using machine-learning algorithms, where high quality ground data were used for training. At present a major limitation for national scale biomass retrieval is the lack of spatial and temporal ground samples, which can be partially resolved by minor modifications in the existing PastureBaseIreland database by adding the location and extent ofeach grassland paddock in the database. As far as remote sensing data requirements are concerned, MODIS is useful for large scale evaluation but due to its coarse resolution it is not possible to detect the variations within the fields and between the fields at the farm scale. However, this issue will be resolved in terms of spatial resolution by the Sentinel-2 mission, and when both satellites (Sentinel-2A and Sentinel-2B) are operational the revisit time will reduce to 5 days, which together with Landsat-8, should enable sufficient cloud-free data for operational biomass estimation at a national scale. The Synthetic Aperture Radar Interferometry (InSAR) approach is feasible if there are enough coherent interferometric pairs available, however this is difficult to achieve due to the temporal decorrelation of the signal. For repeat-pass InSAR over a vegetated area even an 11 days temporal baseline is too large. In order to achieve better coherence a very high resolution is required at the cost of spatial coverage, which limits its scope for use in an operational context at a national scale. Future InSAR missions with pair acquisition in Tandem mode will minimize the temporal decorrelation over vegetation areas for more focused studies. The proposed approach complements the current paradigm of Big Data in Earth Observation, and illustrates the feasibility of integrating data from multiple sources. In future, this framework can be used to build an operational decision support system for retrieval of grassland biophysical parameters based on data from long term planned optical missions (e. g., Landsat, Sentinel) that will ensure the continuity of data acquisition. Similarly, Spanish X-band PAZ and TerraSAR-X2 missions will ensure the continuity of TerraSAR-X and COSMO-SkyMed.

Multi-scale analysis of active landslides using two-pass differential interferometry

Landslides are common features of the landscape of the north-central Apennine mountain range and cause frequent damage to human facilities and infrastructure. Most of these landslides move periodically with moderate velocities and, only after particular rainfall events, some accelerate abruptly. Synthetic aperture radar interferometry (InSAR) provides a particularly convenient method for studying deforming slopes. We use standard two-pass interferometry, taking advantage of the short revisit time of the Sentinel-1 satellites. In this paper we present the results of the InSAR analysis developed on several study areas in central and Northern Italian Apennines. The aims of the work described within the articles contained in this paper, concern: i) the potential of the standard two-pass interferometric technique for the recognition of active landslides; ii) the exploration of the potential related to the displacement time series resulting from a two-pass multiple time-scale InSAR analysis; iii) the evaluation of the possibility of making comparisons with climate forcing for cognitive and risk assessment purposes. Our analysis successfully identified more than 400 InSAR deformation signals (IDS) in the different study areas corresponding to active slope movements. The comparison between IDSs and thematic maps allowed us to identify the main characteristics of the slopes most prone to landslides. The analysis of displacement time series derived from monthly interferometric stacks or single 6-day interferograms allowed the establishment of landslide activity thresholds. This information, combined with the displacement time series, allowed the relationship between ground deformation and climate forcing to be successfully investigated. The InSAR data also gave access to the possibility of validating geographical warning systems and comparing the activity state of landslides with triggering probability thresholds.