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.     

2D Self-Potential tomographies for studying groundwater flows in the Varco d''Izzo landslide (Basilicata, southern Italy)

A geoelectrical monitoring activity has been carried out to improve the geological and hydrogeological knowledge about the Varco d'Izzo landslide (Potenza, Basilicata, Southern Apennine, Italy), an active rotational–translational slide evolved in earth-flow. In this work we have focused on the Self-Potential (SP) method by applying three different SP measuring techniques and combining modern technologies for data acquisition and new methods for tomographic inversion. A SP map and three static SP tomographies have been carried out to better analyse the groundwater circulation system and to better reconstruct the geometry of the landslide body. In the accumulation zone, which is the area more exposed to the geomorphological activity, a new SP measuring strategy has been applied. This strategy, based on time-continuous 2D SP tomographies, helps identify water flow changes in subsurface by studying the time series of SP tomographic images. The analysis of time-dependent changes of water infiltration in near surface is the key to better understand the hydrogeological processes underlying the ground instability phenomena. The time-lapse analysis of tomographic images has allowed us not only to investigate the correlation between the temporal changes of SP signals and rainfall events, but also to quantify the range of these changes. The modification of the distribution of the SP source accumulation zones is associated with the dynamics of the groundwater flows. These preliminary results allow us to consider the SP tomographic method as a tool for geophysical monitoring of landslide areas and encourage to develop new measuring systems for near-real time applications.     

Multi-sensor DInSAR applied to the spatiotemporal evolution analysis of ground surface deformation in Cerro Prieto basin,Baja California,Mexico, for the 1993–2014 period

The combined effects of active tectonics and anthropogenic activities, primarily geothermal resources exploitation for electricity production in Cerro Prieto geothermal field, influence the ground surface deformation in Cerro Prieto basin, Baja California, Mexico. In this study, a large set of multi-sensor C-band SAR images have been employed to reconstruct the spatiotemporal evolution of aseismic ground surface deformation that has affected Cerro Prieto basin from 1993 to 2014. Conventional DInSAR together with the interferograms stacking procedure was applied. The results showed that the study area presented considerable surface deformation (mainly subsidence) during the entire time of the investigation. The main changes in rate and pattern of surface deformation have a good correlation in time and space with the changes in production in the Cerro Prieto geothermal field. Comparison of LOS displacement maps from different viewing geometries, and decomposition (where possible) of LOS displacement into vertical and horizontal (east–west) components, revealed considerable horizontal displacement which mostly reflects the ground movement at and beyond the margin of the subsidence basin toward the areas of highest subsidence rates. In addition, the validation of the DInSAR results by comparing them against measurements from leveling surveys was performed, confirming the high reliably of satellite interferometry for the ground surface deformation rate mapping in the study area.     

Integrated geophysical and geomorphological approach to investigate the snowmelt-triggered landslide of Bosco Piccolo village (Basilicata, southern Italy)

Most of Basilicata region in the southern Italian Apennines is characterized by landslides often developing in clayey–marly formations. Many events have been triggered by extreme rainfall or snowmelt. The most important happened (on February–March 2005) at Bosco Piccolo 5 km far from Potenza. This landslide developed subsequently to rapid snowmelt occurred during alternating short periods of high temperatures and intense and continuous snowfalls. This complex landslide affected 4 ha of surface and reached a maximum depth of 20 m inducing damage and collapse of about 80% of the buildings in the village. An integrated multidisciplinary approach has been adopted to study the landslide. A multitemporal aerial photo interpretation and detailed geomorphological surveys have been carried out. Combined Electrical Resistivity Tomography (ERT) and Self-Potential (SP) measurements have been performed and calibrated with boreholes stratigraphy. Such an integrated approach allowed us to partially reconstruct the geometry of the investigated body and to evaluate the effectiveness of drainage system planned for the area.     

Application of electrical resistivity tomography (ERT) technique to detect underground cavities in a karst area of South Korea

Recently, the electrical resistivity tomography (ERT) technique has been increasingly applied to underground cavity detection filled with groundwater and/or clay. In this study, an ERT survey was conducted to examine the spatial distribution and shape of underground cavities in a karst area in Korea. Numerical ERT modeling and inversion were carried out to select appropriate field survey parameters. Two geological structures, similar to those that exist in the field survey area, were modeled: a sinkhole model and a vein-type structure. Electrode configurations using pole–pole, pole–dipole, dipole–pole, and dipole–dipole arrays were tested with the two geologic models and compared. The numerical modeling showed that the resistivity distribution of the ERT tomogram using the pole–dipole array produced the best representation of the two models. The pole–dipole and pole–pole array configurations were applied in an ERT field survey using six boreholes. The field ERT results show that cavity areas filled with clay or groundwater appeared as low-resistivity anomalies in the limestone formation. This was particularly evident in the inclined vein-type structures attributed to fracture zones, where the field data were similar to the modeling results. Therefore, the ERT technique could be effectively used for detecting underground limestone cavities under similar geological conditions as those at our study area.     

Combining field data with infrared thermography and DInSAR surveys to evaluate the activity of landslides: the case study of Randazzo Landslide (NE Sicily)

Infrared thermography and DInSAR surveying methodologies have been integrated and compared to study the state and distribution of activity of one of the largest landslides of northeastern Sicily, Italy, which is characterized by a slow retrogressive evolution. The Randazzo Landslide was triggered by heavy rainfalls in 1996 and its activity caused the disruption of a strategic transportation corridor, along with the formation of a landslide dam in the Alcantara Valley. After more than 20 years, the slope is still affected by instability and worrying signs of reactivations (e.g., deformation and cracking of road pavement) were surveyed in the field. The application of infrared thermography during the landslide survey allowed detecting peculiar features according to the different thermal pattern of the elements occurring along the slope. In particular, areas with different surface temperature were associated to vegetated spots, steep portions, bare sectors, and contact between different lithologies. Incipient portions and ancient landslide bodies were located, proving the utility of this experimental approach to the surveying of unstable slopes. Information on the entity, in terms of velocity, of the recent landslide displacements were provided by a DInSAR monitoring, which allowed ascertaining the presence of movements affecting several sectors of the landslide, mainly matching with the areas highlighted by thermal images. Achieved outcomes represent a scientifically relevant datum providing important information on the surveyed landslide from the risk management point of view and proving the utility of integrating infrared thermography and DInSAR methodologies for the study of complex movements.     

PS技术及其在地表形变监测中的应用现状与发展

合成孔径雷达干涉测量技术（InSAR）为地表形变监测提供了极具应用潜力的手段,具有大面积、高空间分辨率、全天候及成本低的优点。但由于大气条件变化、地表覆被等时间空间去相干的影响,其精度和普适性受到极大限制。 近几年发展的永久散射体（PS）技术在传统差分干涉测量（DInSAR）中引入时间维,分析长时间内保持稳定的像元集相位变化,获得毫米级的地表形变测量精度,同时有效地解决了时间空间去相关和大气非均质性影响的问题,目前在滑坡、地面沉降和地质灾害监测等领域得到了广泛的应用。PS技术具有高精度、高时间分辨率、能极大提高影像利用率的优点;但只适用较小区域、需要大量影像、且不适于分析快速突变的地表形变。为克服PS应用中的问题,近年来出现了三角反射体技术、多平台PS技术及相关性像元分析（CPT）技术,使PS技术应用具有更广泛的适应性。     

Spatio-temporal analysis of rockfall pre-failure deformation using Terrestrial LiDAR

We present a long-term spatio-temporal analysis of rock slope evolution using a Terrestrial LiDAR aiming to improve our understanding of the link between pre-failure deformation and the spatial prediction of rockfalls. We monitored the pilot study area located at the Puigcercós cliff (Catalonia, Spain) over a period of 1,705 days and detected the deformation of nine different cliff regions together with a high rockfall activity. An exact match was observed between the progressively deformed areas and the regions recently affected by three of the highest magnitude rockfall events, demonstrating a causal relationship between pre-failure deformation and rockfall occurrence. These findings allowed us to make a forward spatial prediction of future failures, hypothesizing a high probability of failure in the six remaining regions. We observed an exponential acceleration of the deformation close to failure, in accordance with tertiary creep theory. However, the temporal analysis of the deformed areas showed a complex and variable behavior, so no exact prediction of the date of failure can yet be made. Our findings have broadened our understanding of the pre-failure behavior of rockfalls and have clear implications for the future implementation of early warning systems.     

Complementing geotechnical slope stability and land movement analysis using satellite DInSAR

This paper explores the potential of using satellite radar inteferometry to monitor time-varying land movement prior to any visible tension crack signs. The idea was developed during dedicated geotechnical studies at a large open-pit lignite mine, where large slope movements (10–20 mm/day) were monitored and large fissures were observed in the immediate area outside the current pit limits. In this work, differential interferometry (DInSAR), using Synthetic Aperture Radar (SAR) ALOS images, was applied to monitor the progression of land movement that could potentially thwart mine operations. Early signs of land movements were captured by this technique well before their visual observation. Moreover, a qualitative comparison of DInSAR and ground geodetic measurements indicates that the technique can be used for the identification of high risk areas and, subsequently, for the optimization of the spatial distribution of the available ground monitoring equipment. Finally, quantitative land movement results from DInSAR are shown to be in accordance with simultaneous measurements obtained by ground means.     

Multiyear time-lapse ERT to study short- and long-term landslide hydrological dynamics

The geology of the “Vence” landslide (0.8 million m³, south-eastern France) explains the complex hydrology of the site which plays a key role in the destabilization of the slope (water circulation within the sliding mass, fluid exchanges between superficial layers and deep karstic aquifer through faults). To understand fluid circulations within the unstable slope, a 9.5-year multi-parametric survey was set up. The survey combines electrical resistivity tomography (daily acquisition), rainfall records since 2006 and boreholes monitoring groundwater level since 2009. The objective of this work is to present an automated clustering analysis applied to the ERT data enabled to locate geological units displaying distinct hydrogeological behaviours. Clustering analysis, based on a hierarchical ascendant classification (HAC), helped to simplify the ERT section isolating three groups of apparent resistivity values. Comparing the variations of these clusters’ behaviours in time to the variations of the groundwater levels on site, we identified hydrogeological units. The role of the faults cutting the substratum is thereby highlighted. It is the simultaneous analysis of such a large real dataset that allowed obtaining robust results characteristic of the long-term behaviour of the natural hydrogeological system. This type of qualitative information on the variability of the slope hydrogeological behaviour both spatially and temporally is crucial to help improving the conversion of resistivity data into hydrologic quantities. Indeed, the definition of petrophysical models to convert ERT measurements into hydrological measurements should be site-specific and take into account the spatial and temporal variability of the medium. In this work, we show a method that can also help to focus on the areas in depth that have different levels of permeability and observe how the saturation degree evolves in time. This can be used to optimize the location of additional instrumentation (such as temperature probes and chemical sampling) and, thus, help in the prevention of the risk in such problematic areas.     

Integrated analysis of SAR interferometric and geological data for investigating long-term reclamation settlement of Chek Lap Kok Airport,Hong Kong

We investigated the long-term reclamation settlement of Chek Lap Kok Airport, Hong Kong, one of the largest reclamation projects worldwide. An advanced InSAR technique was applied to two independent ASAR data sets along ascending (25 scenes) and descending (22 scenes) orbits acquired between 2003 and 2008, and allowed us to generate deformation velocity maps and displacement time series for the investigated area. Intercomparison of deformation velocity results from the two ASAR data sets and their validation against leveling survey measurements were carried out in order to assess accuracy and reliability of the InSAR observations. The intercomparison and validation results demonstrated the potential of elaborate InSAR processes on local-scale areas to monitor and detect a long-term residual settlement in such complex reclamation area. The InSAR-derived deformation velocity maps showed a fair stability in areas of the two original rocky islands of Chek Lap Kok and Lam Chau, and on most of the buildings at the site due to foundation on piles reaching the granite bedrock, whereas a generally moderate (3–7 mm/year) to strong (> 10 mm/year) and spatially variable settlement trend was detected in the majority of reclamation areas of the airport. The consistent relationship we found between the observed settlement pattern and the underlying geological conditions suggested that this recoded differential settlement was mainly related to the variations in thickness and soil properties of both the fill materials and the alluvial deposits below the reclamation. A geological and geotechnical interpretation of the lithostratigraphy of the alluvial deposits beneath the reclamation implied that the relatively large settlement rate measured up to 2008 was resulted from a still active primary consolidation process.     

Using differential SAR interferometry to map land subsidence: a case study in the Pingtung Plain of SW Taiwan

Synthetic aperture radar (SAR) interferometry (InSAR) is a geodetic tool widely applied in the studies of earth-surface deformation. This technique has the benefits of high spatial resolution and centimetre-scale accuracy. Differential SAR interferometry (DInSAR) is used to measure ground deformation with repeat-pass SAR images. This study applied DInSAR and persistent scatterers InSAR (PSInSAR) for detecting land subsidence in the Pingtung Plain, southern Taiwan, between 1995 and 2000. In recent years, serious land subsidence occurred along coastal regions of Taiwan as a consequence of over-pumping of underground water. Results of this study revealed that the critical subsidence region is located on the coast near the estuary of Linpien River. It is also found that subsidence was significantly higher during the dry season than the wet season. The maximum annual subsidence rate of the dry season is up to −11.51 cm/year in critical subsidence region and the vertical land movement rate is much slower during the wet season. The average subsidence rates in wet and dry seasons are −0.31 and −3.37 cm/year, respectively. As a result, the subsidence rate in dry seasons is about 3 cm larger than in wet seasons.     

高密度电法在输电线路塔基基础附近多年冻土探测中的应用北大核心CSCD

多年冻土区输电线路塔基基础附近活动层厚度和地下冰变化与基础稳定性密切相关,塔基施工的热扰动和混凝土基础的热效应使得基础周围冻土易发生退化,不利于基础的稳定。高密度电法是冻土工程环境研究中常用的地球物理方法,其探测结果的可靠性和分辨能力受数据采集方式、目标体地电结构影响。为减小对输电线路塔基附近冻土特征识别的不确定性,通过建立基础周围多年冻土地电模型的正反演模拟,发现活动层处于融化状态时各种装置方式数据采集均能较好地反映活动层厚度的起伏,但由于冻融锋面附近显著的电阻率差异,难以识别多年冻土层内的地下冰空间分布特征。而活动层处于冻结状态时进行探测能显著提高对多年冻土层内的地下冰空间分布特征识别精度,其中偶极-偶极装置可较好地识别高、低含冰量区域的发育位置和形态特征。在青藏直流输电线路塔基基础附近冻土探测中证实了方法的有效性,探测结果揭示了施工过程和基础热效应导致的塔基基础附近的地下冰退化。以上研究表明,通过正反演模拟,根据具体探测目标选择合适的探测时机和数据采集方式,能显著提高高密度电法探测结果的有效性和精度。     

Phase and amplitude analyses of SAR data for landslide detection and monitoring in non-urban areas located in the North-Eastern Italian pre-Alps

The main aim of this paper is to exploit information obtained from satellite SAR data to detect and monitor instability phenomena affecting hilly and scarcely urbanized areas, overtaking some of the restrictions due to the presence of thick vegetation. To this end, phase and amplitude analyses of COSMO-SkyMed SAR data were carried out on two landslides located in the North-Eastern Italian pre-Alps: Cischele roto-translational slide and Val Maso rotational slide—earth flow. In the first case, the careful choice of processing parameters allowed to evaluate landslide displacement fields considering the phase difference between SAR acquisitions. In the second case, the speed of movement and the deep changes in morphology and vegetation induced by the landslide did not allow to apply DInSAR techniques; in this case the variation in the amplitude between SAR acquisitions allowed to detect the area affected by the instability. Obtained results show that methods and techniques to analyse satellite SAR data could be further refined in order to provide useful tools for landslide mapping and monitoring.     

Investigating landslides with space-borne Synthetic Aperture Radar (SAR) interferometry

This paper is addressed to readers without advanced knowledge of remote sensing. It illustrates some current and potential uses of satellite Synthetic Aperture Radar interferometry (InSAR) for landslide assessment. Data acquired by SAR systems can provide 3D terrain models and be used to assist in regional scale investigations, e.g. aimed at evaluation of susceptibility of slopes to failure. Under favourable environmental conditions, the innovative Permanent Scatterers (PS) technique, which overcomes several limitations of conventional SAR differential interferometry (DInSAR) applications in landslide studies, is suitable for monitoring slope deformations with millimetric precision. The PS technique combines the wide-area coverage typical of satellite imagery with the capability of providing displacement data relative to individual image pixels. With the currently available radar satellites, however, only very slow ground surface displacements can be reliably detected and measured. The presented case study of a landslide from the Liechtenstein Alps indicates that the most attractive and reliable contribution provided by this remote sensing technique lies in the possibility of (i.) wide-area qualitative distinction between stable and unstable areas and (ii.) qualitative (relative) hazard zonation of large, slow landslides based on the identification of segments characterised by different movement rates. Since only the radar line of sight projection of the displacements can be detected, a quantitative exploitation of the PS data is possible only where sufficient ground truth is available. In site specific or single landslide investigations the PS data can represent a very useful complementary data source with respect to the information acquired through ground based observations and in situ surveying. However, the difficulties associated with the feasibility assessments of the applicability of SAR data to local scale problems, as well as with the interpretation of PS results, require a close collaboration between landslide experts and specialists in advanced processing of radar satellite data. The interpretation of the exact geotechnical significance of small, radar sensed ground surface deformations is challenging, especially where ground truth is lacking. Although any ground deformation is potentially of interest to an engineering geologist, detection of movements in both vertical and horizontal directions is needed in the case of landslides to evaluate slope failure mechanisms. With their high radar viewing angles, however, the current space-borne systems can detect only a fraction of the horizontal component of movement. It is expected that the upcoming SAR dedicated missions with new sensors and different acquisition geometries, combined with the rapid developments in the field of advanced radar data processing, will allow a full 3D reconstruction of deformation data and help to further reduce the current limitations of the PS and similar DInSAR approaches.     

DInSAR-based monitoring of land subsidence related to groundwater over-exploitation: example from developing urban center of Nairobi,Kenya

Over-exploitation of groundwater in many evolving urban settings causes ground subsidence and permanent loss of aquifer storage capacity. DInSAR (differential interferometric synthetic aperture radar) time series data from 2016 to 2019 were used to monitor and model the surface deformation around Nairobi, Kenya, where the water demand has exceeded the supply without capacity augmentation for over two decades. The aquifer system constitutes hard rock to semiconfined ash beds in volcanic terrain. The Small Baseline DInSAR technique identified the spatial pattern of subsidence and magnitude (line-of-sight (LOS) velocity), which exceeds 41 mm/year in the semiconfined aquifer towards the western-central part of Nairobi. The spatial distribution of subsidence is consistent with the groundwater level drop and probable compaction modeled using aquifer characteristics for 1950–2015. The Global Navigation Satellite System (GNSS) data at a station from 2007 to 2018 indicate a cumulative 4-cm subsidence which is comparable to ~2.5-cm LOS subsidence from the present study for 2016–2019. The correlation with other hydrological data suggests the aquifer is experiencing inelastic subsidence due to unsustainable groundwater extraction, putting a massive strain on Nairobi’s aquifer system. The present DInSAR based study establishes its effectiveness in the monitoring of groundwater over-exploitation-based subsidence and associated hazard to the aquifer in emerging urban centers.

     

Integrated near surface geophysics across the active Mount Marzano Fault System (southern Italy): seismogenic hints

Here, we describe an original geophysical multi-method approach applied to the Mount Marzano Fault System. This is one of the most hazardous seismogenic faults of the Apennines (Irpinia, southern Italy), and it was responsible for the 1980, Mw 6.9, earthquake, along with many others before. We carried out electrical resistivity tomography (ERT), ground penetrating radar (GPR) measurements, and horizontal-to-vertical spectral ratio (HVSR) microtremor analysis along several common transects designed across the potential and/or certain fault traces. The data obtained from these non-invasive, inexpensive, expeditious methods mutually integrate with and complement each other, providing a valuable subsurface image of the near surface fault architecture. ERT depicts the general shallow image of the fault zone and of the fault-controlled sedimentary basin, with the depth of the buried bedrock cross-correlated through ambient-noise HVSR results. GPR delineates the very shallow geometry of the fault and of the associated deformation. Coupled with previous paleoseismological studies, these data allow the evaluation of some fault parameters and the precise locating of the fault trace, to aid future paleoseismological investigations aimed at seismic risk reduction programs.     

Ground motion phenomena in Caltanissetta (Italy) investigated by InSAR and geological data integration

Urban areas are frequently affected by ground instabilities of various origins. The location of urban zones affected by ground instability phenomena is crucially important for hazard mitigation policies. Satellite-based Interferometric Synthetic Aperture Radar (InSAR) has demonstrated its remarkable capability to detect and quantify ground and building motion in urban areas, especially since the development of Advanced Differential Interferometric SAR techniques (A-DInSAR). In fact, the high density of reflectors like buildings and infrastructures in urban areas improves the quality of the InSAR signal, allowing sub-centimetric displacements to be reliably detected. The A-DInSAR techniques allow urban zones affected by ground deformation to be located and mapped, but clearly they are not able to point out the causes of the instability phenomena. These can only be highlighted by an integrated analysis of multidisciplinary data, like geological, geotechnical, SAR interferometric and historical data. The overlay of these data, which is possible within a Geographic Information System (GIS), is a useful tool to identify ground motion phenomena affecting urban zones. In this study we apply this kind of approach to Caltanissetta, a provincial capital in Sicily (Italy), where local damage has been detected. The reconstruction of the local near-surface geology shows the presence of zones affected by local natural hazard factors, essentially due to the local presence of soils with poor mechanical properties or swelling soils, high topographic gradients and steep slopes on loose soils. Processing 17 ASAR-ENVISAT SAR images covering the time interval October 2002–December 2005 by means of an A-DInSAR procedure, the Caltanissetta deformation map has been realized. It shows that most of the city is stable, with the exception of three zones, situated in the northwestern, northeastern and southern parts of the city, respectively. Two of them, characterized by high topographic gradients and steep slopes on sandy soils, are affected by subsidence ground motion. An uplift motion is recognized in the other zone, characterized by the local presence of expansible clays. Geotechnical swelling tests carried out on them have shown a swelling behavior. On site surveys have highlighted the presence of damage in the zones affected by ground motion.     

Investigation of waste disposal areas using electrical methods: a case study from Chania,Crete, Greece

A geophysical investigation of a landfill area was conducted as inextricable phase of a preliminary geotechnical, geological and hydrogeological study of the area under investigation. Electrical resistivity tomography (ERT), one of the most promising prospecting techniques mainly concerning its effective contribution to resolve several environmental problems, was applied for the geophysical modeling. ERT is a robust imaging method the theory and implementation of which are well documented in geophysical research literature. The survey site is located in Akrotirion Peninsula, about 7 km east of Chania city in Crete I. The geological setting comprises Trypalion limestones, marly limestones and a clayed weathered layer. Cavities are also identified at various scales which are recent tectono-karstic voids. Due to the above-mentioned geological singularity and in the light of the requirement for an environmentally safe construction of the landfill, an ERT survey was carried out. Specifically, seven geoelectrical tomographies were conducted reaching the prospecting depth of 36 m. The resulted images conduced, to overcome geotechnical problems since the spatial distribution of karstic voids was determined, to plan the future sites for waste disposal as the geological conditions were studied in detail and to reliably estimate the thickness of the already deposited wastes. The resulted images were confirmed using available borehole logs and ambient noise measurements.     

Magnetic and ground probing radar measurements for soil pollution mapping in the industrial area of Val Basento (Basilicata Region, Southern Italy): a case study

The joint application of electromagnetic techniques for near-surface exploration is a useful tool for soil pollution monitoring and can also contribute towards describing the spatial distribution of pollutants. The results of a geophysical field survey that was carried out for characterizing the heavy metal and waste disposal soil pollution phenomena in the industrial area of Val Basento (Basilicata region, Southern Italy) are presented here. First, topsoil magnetic susceptibility measurements have been carried out for defining the spatial distribution of superficial pollution phenomena in the investigated area. Second, detailed and integrated measurements based on a high-resolution magnetic mapping and ground probing radar (GPR) profiling have been applied to investigate the subsurface in two industrial areas located in more polluted sites that were identified during the first phase. Our monitoring strategy discloses the way to rapidly define the zone characterized by high pollution levels deriving from chemical industries and traffic emissions and to obtain the way information about the presence of local buried sources of contamination.