Search

Optical and Synthetic Aperture Radar (SAR) remote sensing has a long history of use and reached a good level of maturity in archaeological and cultural heritage applications, yet further advances are viable through the exploitation of novel sensor data and imaging modes, big data and high-performance computing, advanced and automated analysis methods. This paper showcases the main research avenues in this field, with a focus on archaeological prospection and heritage site protection. Six demonstration use-cases with a wealth of heritage asset types (e.g. excavated and still buried archaeological features, standing monuments, natural reserves, burial mounds, paleo-channels) and respective scientific research objectives are presented: the Ostia-Portus area and the wider Province of Rome (Italy), the city of Wuhan and the Jiuzhaigou National Park (China), and the Siberian “Valley of the Kings” (Russia). Input data encompass both archive and newly tasked medium to very high-resolution imagery acquired over the last decade from satellite (e.g. Copernicus Sentinels and ESA Third Party Missions) and aerial (e.g. Unmanned Aerial Vehicles, UAV) platforms, as well as field-based evidence and ground truth, auxiliary topographic data, Digital Elevation Models (DEM), and monitoring data from geodetic campaigns and networks. The novel results achieved for the use-cases contribute to the discussion on the advantages and limitations of optical and SAR-based archaeological and heritage applications aimed to detect buried and sub-surface archaeological assets across rural and semi-vegetated landscapes, identify threats to cultural heritage assets due to ground instability and urban development in large metropolises, and monitor post-disaster impacts in natural reserves.

ABSTRACTSatellite Interferometric Synthetic Aperture Radar (InSAR) is widely used for topographic, geological and natural resource investigations. However, most of the existing InSAR studies of ground deformation are based on relatively short periods and single sensors. This paper introduces a new multi-sensor InSAR time series data fusion method for time-overlapping and time-interval datasets, to address cases when partial overlaps and/or temporal gaps exist. A new Power Exponential Knothe Model (PEKM) fits and fuses overlaps in the deformation curves, while a Long Short-Term Memory (LSTM) neural network predicts and fuses any temporal gaps in the series. Taking the city of Wuhan (China) as experiment area, COSMO-SkyMed (2011–2015), TerraSAR-X (2015–2019) and Sentinel-1 (2019–2021) SAR datasets were fused to map long-term surface deformation over the last decade. An independent 2011–2020 InSAR time series analysis based on 230 COSMO-SkyMed scenes was also used as reference for comparison. The correlation coefficient between the results of the fusion algorithm and the reference data is 0.87 in the time overlapping region and 0.97 in the time-interval dataset. The correlation coefficient of the overall results is 0.78, which fully demonstrates that the algorithm proposed in our paper achieves a similar trend as the reference deformation curve. The experimental results are consistent with existing studies of surface deformation at Wuhan, demonstrating the accuracy of the proposed new fusion method to provide robust time series for the analysis of long-term land subsidence mechanisms.

One of the most promising applications of satellite data is providing users in charge of land and emergency management with information and data to support decision making for geohazard mapping, monitoring and early warning. In this work, we consider ground displacement data obtained via interferometric processing of satellite radar imagery, and we provide a novel post-processing approach based on a Functional Data Analysis paradigm capable of detecting precursors in displacement time series. The proposed approach appropriately accounts for the spatial and temporal dependencies of the data and does not require prior assumptions on the deformation trend. As an illustrative case, we apply the developed method to the identification of precursors to a mud volcano eruption in the Santa Barbara village in Sicily, southern Italy, showing the advantages of using a Functional Data Analysis framework for anticipating the warning signal. Indeed, the proposed approach is able to detect precursors of the paroxysmal event in the time series of the locations close to the eruption vent and provides a warning signal months before a scalar approach would. The method presented can potentially be applied to a wide range of geological events, thus representing a valuable and far-reaching monitoring tool.

The quantification of the extent and dynamics of land-use changes is a key metric employed to assess the progress toward several Sustainable Development Goals (SDGs) that form part of the United Nations 2030 Sustainable Development Agenda. In terms of anthropogenic factors threatening the conservation of heritage properties, such a metric aids in the assessment of achievements toward heritage sustainability solving the problem of insufficient data availability. Therefore, in this study, 589 cultural World Heritage List (WHL) properties from 115 countries were analyzed, encompassing globally distributed and statistically significant samples of “monuments and groups of buildings” (73.2%), “sites” (19.3%), and “cultural landscapes” (7.5%). Land-cover changes in the WHL properties between 2015 and 2020 were automatically extracted from big data collections of high-resolution satellite imagery accessed via Google Earth Engine using intelligent remote sensing classification. Sustainability indexes (SIs) were estimated for the protection zones of each property, and the results were employed, for the first time, to assess the progress of each country toward SDG Target 11.4. Despite the apparent advances in SIs (10.4%), most countries either exhibited steady (20.0%) or declining (69.6%) SIs due to limited cultural investigations and enhanced negative anthropogenic disturbances. This study confirms that land-cover changes are among serious threats for heritage conservation, with heritage in some countries wherein the need to address this threat is most crucial, and the proposed spatiotemporal monitoring approach is recommended.

Early-season crop mapping provides decision-makers with timely information on crop types and conditions that are crucial for agricultural management. Current satellite-based mapping solutions mainly rely on optical imagery, albeit limited by weather conditions. Very few exploit long-time series of polarized synthetic aperture radar (SAR) imagery. To address this gap, we assessed the performance of COSMO-SkyMed X-band dual-polarized (HH, VV) data in a test area in Ponte a Elsa (central Italy) in January–September 2020 and 2021. A deep learning convolutional neural network (CNN) classifier arranged with two different architectures (1-D and 3-D) was trained and used to recognize ten classes. Validation was undertaken with in situ measurements from regular field campaigns carried out during satellite overpasses over more than 100 plots each year. The 3-D classifier structure and the combination of HH+VV backscatter provide the best classification accuracy, especially during the first months of each year, i.e., 80% already in April 2020 and in May 2021. Overall accuracy above 90% is always marked from June using the 3-D classifier with HH, VV, and HH+VV backscatter. These experiments showcase the value of the developed SAR-based early-season crop mapping approach. The influence of vegetation phenology, structure, density, biomass, and turgor on the CNN classifier using X-band data requires further investigations, along with the relatively low producer accuracy marked by vineyard and uncultivated fields.

Imaging radar has been dramatically developed over the past decades enabling a better understanding of cultural heritage from a microwave perspective. Nonetheless, a dedicated survey and analysis of the performance of such technology in cultural heritage monitoring and management is required. In order to fill this gap, we first review the technology advance of imaging radar, including ground penetration radar, ground-based and airborne/satellite radar, in the focused cultural applications to grasp the development trend of these technologies. We then analyse the performance and limitations of imaging radar technologies based on their respective characteristics to facilitate the technology service in practical applications. Finally, we propose a flexible solution of imaging radar in cultural heritage through technical integration with pilot synergy applications in archaeological prospection and cultural heritage diagnosis and conservation.

A general framework for the integration of multi-sensor data for dry and fresh biomass retrieval is proposed and tested in Alpine meadows and pastures. To this purpose, hyperspectral spectroradiometer (as simulation of hyperspectral imagery) and biomass samples were collected in field campaigns and Copernicus Sentinel-1 Interferometric Wide (IW) swath SAR backscattering coefficients were used. First, a genetic algorithm feature selection was performed on hyperspectral data, and afterwards the resulting most sensitive bands where combined with SAR data within a support vector regression (SVR) model. The most sensitive hyperspectral bands were mainly located in different regions of the SWIR range for both fresh and dry biomass, and in the red and near-infrared regions mainly for dry biomass, but with less influence for fresh biomass. The R2 correlation values between the sampled and the estimated biomass range from 0.24 to 0.71. The relatively low performances are mainly related to the saturation effect in the optical bands, as well as to the paucity of points for high values of biomass. The methodology allows a better understanding of the interaction between grassland systems and the electromagnetic spectrum by offering a model with a reduced number of narrow bands in the context of a multi-sensor integration.

In Near and Middle Eastern archaeology, satellite-derived digital elevation models (DEM) of medium spatial resolution (≥30 m) are mostly used to locate and map archaeological mounds (namely ‘tells’), whereas high resolution DEMs (≤10 m) are still poorly exploited. To fill this gap, the 5 m resolution CartoSat-1 Euro-Maps 3D Digital Surface Model (DSM) is assessed vs. the 30 m Shuttle Radar Topography Mission (SRTM) global DEM, the Advanced Land Observing Satellite (ALOS) World 3D–30 m (AW3D30) and a 10 m COSMO-SkyMed DEM, on a test area in Wasit, southern Iraq, where the high density of tells is yet to be exhaustively documented. A total of 344 sites was mapped, with one order of magnitude improvement compared to previous mapping exercises, existing databases and historical maps. The morphometric analysis not only highlights the reliability of CartoSat-1 DSM height and volume estimates, but also suggests that, in the test area, the volume of a tell can robustly be calculated based on the simple knowledge of its basal area, following a quadratic function. Morphology and elevation of at least 53% irregularly shaped tells were impacted by anthropogenic disturbances. Morphometric indices (e.g., Topographic Position Index, DEViation from mean elevation) are a viable automated method to ease tells detection. When integrated with other satellite datasets (e.g., CORONA, Google Earth, Sentinel-2 imagery), the CartoSat-1 DSM can unveil morphological changes and support condition assessment. In Wasit, agriculture and modern development are among the major threats for tells preservation, alongside looting.

Volcanic activity may influence climate parameters and impact people safety, and hence monitoring its characteristic indicators and their temporal evolution is crucial. Several databases, communications and literature providing data, information and updates on active volcanoes worldwide are available, and will likely increase in the future. Consequently, information extraction and text mining techniques aiming to efficiently analyze such databases and gather data and parameters of interest on a specific volcano can play an important role in this applied science field. This work presents a natural language processing (NLP) system that we developed to extract geochemical and geophysical data from free unstructured text included in monitoring reports and operational bulletins issued by volcanological observatories in HTML, PDF and MS Word formats. The NLP system enables the extraction of relevant gas parameters (e.g., SO2 and CO2 flux) from the text, and was tested on a series of 2839 daily and weekly bulletins published online between 2015 and 2021 for the Stromboli volcano (Italy). The experiment shows that the system proves capable in the extraction of the time series of a set of user-defined parameters that can be later analyzed and interpreted by specialists in relation with other monitoring and geospatial data. The text mining system can potentially be tuned to extract other target parameters from this and other databases.

Correct use of multi-temporal Interferometric Synthetic Aperture Radar (InSAR) datasets to complement geodetic surveying for geo-hazard applications requires rigorous assessment of their precision and accuracy. Published inter-comparisons are mostly limited to ground displacement estimates obtained from different algorithms belonging to the same family of InSAR approaches, either Persistent Scatterer Interferometry (PSI) or Small BAseline Subset (SBAS); and accuracy assessments are mainly focused on vertical displacements or based on few Global Navigation Satellite System (GNSS) or geodetic leveling points. To fill this demonstration gap, two years of Sentinel-1 SAR ascending and descending mode data are processed with both PSI and SBAS consolidated algorithms to extract vertical and horizontal displacement velocity datasets, whose accuracy is then assessed against a wealth of contextual geodetic data. These include permanent GNSS records, static GNSS benchmark repositioning, and geodetic leveling monitoring data that the National Institute of Statistics, Geography, and Informatics (INEGI) of Mexico collected in 2014−2016 in the Aguascalientes Valley, where structurally-controlled land subsidence exhibits fast vertical rates (up to −150 mm/year) and a non-negligible east-west component (up to ±30 mm/year). Despite the temporal constraint of the data selected, the PSI-SBAS inter-comparison reveals standard deviation of 6 mm/year and 4 mm/year for the vertical and east-west rate differences, respectively, thus reassuring about the similarity between the two types of InSAR outputs. Accuracy assessment shows that the standard deviations in vertical velocity differences are 9−10 mm/year against GNSS benchmarks, and 8 mm/year against leveling data. Relative errors are below 20% for any locations subsiding faster than −15 mm/year. Differences in east-west velocity estimates against GNSS are on average −0.1 mm/year for PSI and +0.2 mm/year for SBAS, with standard deviations of 8 mm/year. When discrepancies are found between InSAR and geodetic data, these mostly occur at benchmarks located in proximity to the main normal faults, thus falling within the same SBAS ground pixel or closer to the same PSI target, regardless of whether they are in the footwall or hanging wall of the fault. Establishing new benchmarks at higher distances from the fault traces or exploiting higher resolution SAR scenes and/or InSAR datasets may improve the detection of the benchmarks and thus consolidate the statistics of the InSAR accuracy assessments.

Sulfur dioxide (SO2) degassing at Strombolian volcanoes is directly associated with magmatic activity, thus its monitoring can inform about the style and intensity of eruptions. The Stromboli volcano in southern Italy is used as a test case to demonstrate that the TROPOspheric Monitoring Instrument (TROPOMI) onboard the Copernicus Sentinel-5 Precursor (Sentinel-5P) satellite has the suitable spatial resolution and sensitivity to carry out local-scale SO2 monitoring of relatively small-size, nearly point-wise volcanic sources, and distinguish periods of different activity intensity. The entire dataset consisting of TROPOMI Level 2 SO2 geophysical products from UV sensor data collected over Stromboli from 6 May 2018 to 31 May 2021 is processed with purposely adapted Python scripts. A methodological workflow is developed to encompass the extraction of total SO2 Vertical Column Density (VCD) at given coordinates (including conditional VCD for three different hypothetical peaks at 0–1, 7 and 15 km), as well as filtering by quality in compliance with the Sentinel-5P Validation Team’s recommendations. The comparison of total SO2 VCD time series for the main crater and across different averaging windows (3 × 3, 5 × 5 and 4 × 2) proves the correctness of the adopted spatial sampling criterion, and practical recommendations are proposed for further implementation in similar volcanic environments. An approach for detecting SO2 VCD peaks at the volcano is trialed, and the detections are compared with the level of SO2 flux measured at ground-based instrumentation. SO2 time series analysis is complemented with information provided by contextual Sentinel-2 multispectral (in the visible, near and short-wave infrared) and Suomi NPP VIIRS observations. The aim is to correctly interpret SO2 total VCD peaks when they either (i) coincide with medium to very high SO2 emissions as measured in situ and known from volcanological observatory bulletins, or (ii) occur outside periods of significant emissions despite signs of activity visible in Sentinel-2 data. Finally, SO2 VCD peaks in the time series are further investigated through daily time lapses during the paroxysms in July–August 2019, major explosions in August 2020 and a more recent period of activity in May 2021. Hourly wind records from ECMWF Reanalysis v5 (ERA5) data are used to identify local wind direction and SO2 plume drift during the time lapses. The proposed analysis approach is successful in showing the SO2 degassing associated with these events, and warning whenever the SO2 VCD at Stromboli may be overestimated due to clustering with the plume of the Mount Etna volcano.

The aim of this paper is to address the monitoring of the recovery phase in the aftermath of Hurricane Matthew (28 September–10 October 2016) in the town of Jérémie, southwestern Haiti. This is accomplished via a novel change detection method that has been formulated, in a data fusion perspective, in terms of multitemporal supervised classification. The availability of very high resolution images provided by last-generation satellite synthetic aperture radar (SAR) and optical sensors makes this analysis promising from an application perspective and simultaneously challenging from a processing viewpoint. Indeed, pursuing such a goal requires the development of novel methodologies able to exploit the large amount of detailed information provided by this type of data. To take advantage of the temporal and spatial information associated with such images, the proposed method integrates multisensor, multisource, and contextual information. Markov random field modeling is adopted here to integrate the spatial context and the temporal correlation associated with images acquired at different dates. Moreover, the adoption of a region-based approach allows for the characterization of the geometrical structures in the images through multiple segmentation maps at different scales and times. The performances of the proposed approach are evaluated on multisensor pairs of COSMO-SkyMed SAR and Pléiades optical images acquired over Jérémie, in the aftermath of and during the three years after Hurricane Matthew. The effectiveness of the change detection results is analyzed both quantitatively, through the computation of accuracy measures on a test set, and qualitatively, by visual inspection of the classification maps. The robustness of the proposed method with respect to different algorithmic choices is also assessed, and the detected changes are discussed in relation to the recovery endeavors in the area and ground-truth data collected in the field in April 2019.

“Tells” are archaeological mounds formed by deposition of large amounts of anthropogenic material and sediments over thousands of years and are the most important and prominent features in Near and Middle Eastern archaeological landscapes. In the last decade, archaeologists have exploited free-access global digital elevation model (DEM) datasets at medium resolution (i.e., up to 30 m) to map tells on a supra-regional scale and pinpoint tentative tell sites. Instead, the potential of satellite DEMs at higher resolution for this task was yet to be demonstrated. To this purpose, the 3 m resolution imaging capability allowed by the Italian Space Agency’s COSMO-SkyMed Synthetic Aperture Radar (SAR) constellation in StripMap HIMAGE mode was used in this study to generate DEM products of enhanced resolution to undertake, for the first time, a systematic mapping of tells and archaeological deposits. The demonstration is run at regional scale in the Governorate of Wasit in central Iraq, where the literature suggested a high density of sites, despite knowledge gaps about their location and spatial distribution. Accuracy assessment of the COSMO-SkyMed DEM is provided with respect to the most commonly used SRTM and ALOS World 3D DEMs. Owing to the 10 m posting and the consequent enhanced observation capability, the COSMO-SkyMed DEM proves capable to detect both well preserved and levelled or disturbed tells, standing out for more than 4 m from the surrounding landscape. Through the integration with CORONA KH-4B tiles, 1950s Soviet maps and recent Sentinel-2 multispectral images, the expert-led visual identification and manual mapping in the GIS environment led to localization of tens of sites that were not previously mapped, alongside the computation of a figure as up-to-date as February 2019 of the survived tells, with those affected by looting. Finally, this evidence is used to recognize hot-spot areas of potential concern for the conservation of tells. To this purpose, we upgraded the spatial resolution of the observations up to 1 m by using the Enhanced Spotlight mode to collect a bespoke time series. The change detection tests undertaken on selected clusters of disturbed tells prove how a dedicated monitoring activity may allow a regular observation of the impacts due to anthropogenic disturbance (e.g., road and canal constructions or ploughing).

Wuhan is an important city in central China, with a rapid development that has led to increasingly serious land subsidence over the last decades. Most of the existing Interferometric Synthetic Aperture Radar (InSAR) subsidence monitoring studies in Wuhan are either short-term investigations—and thus can only detect this process within limited time periods—or combinations of different Synthetic Aperture Radar (SAR) datasets with temporal gaps in between. To overcome these constraints, we exploited nearly 300 high-resolution COSMO-SkyMed StripMap HIMAGE scenes acquired between 2012 and 2019 to monitor the long-term subsidence process affecting Wuhan and to reveal its spatiotemporal variations. The results from the Persistent Scatterer Interferometric SAR (PSInSAR) processing highlight several clearly observable subsidence zones. Three of them (i.e., Houhu, Xinrong, and Guanggu) are affected by serious subsidence rates and non-linear temporal behavior, and are investigated in this paper in more detail. The subsidence in Houhu is caused by soft soil consolidation and compression. Soil mechanics are therefore used to estimate when the subsidence is expected to finish and to calculate the degree of consolidation for each year. The COSMO-SkyMed PSInSAR results indicate that the area has entered the late stage of consolidation and compression and is gradually stabilizing. The subsidence curve found for the area around Xinrong shows that the construction of an underground tract of the subway Line 21 caused large-scale settlement in this area. The temporal granularity of the PSInSAR time series also allows precise detection of a rebound phase following a major flooding event in 2016. In the southern industrial park of Guanggu, newly detected subsidence was found. The combination of the subsidence curve with an optical time-series image analysis indicates that urban construction is the main trigger of deformation in this area. While this study unveils previously unknown characters of land subsidence in Wuhan and clarifies the relationship with the urban causative factors, it also proves the benefits of non-linear PSInSAR in the analysis of the temporal evolution of such processes in dynamic and expanding cities.

The growing volume of synthetic aperture radar (SAR) imagery acquired by satellite constellations creates novel opportunities and opens new challenges for interferometric SAR (InSAR) applications to observe Earth’s surface processes and geohazards. In this paper, the Parallel Small BAseline Subset (P-SBAS) advanced InSAR processing chain running on the Geohazards Exploitation Platform (GEP) is trialed to process two unprecedentedly big stacks of Copernicus Sentinel-1 C-band SAR images acquired in 2014–2020 over a coastal study area in southern Italy, including 296 and 283 scenes in ascending and descending mode, respectively. Each stack was processed in the GEP in less than 3 days, from input SAR data retrieval via repositories, up to generation of the output P-SBAS datasets of coherent targets and their displacement histories. Use-cases of long-term monitoring of land subsidence at the Capo Colonna promontory (up −2.3 cm/year vertical and −1.0 cm/year east–west rate), slow-moving landslides and erosion landforms, and deformation at modern coastal protection infrastructure in the city of Crotone are used to: (i) showcase the type and precision of deformation products outputting from P-SBAS processing of big data, and the derivable key information to support value-adding and geological interpretation; and (ii) discuss potential and challenges of big data processing using cloud/grid infrastructure.

Geohazards associated with the dynamics of the liquid and solid water of the Earth’s hydrosphere, such as floods and glacial processes, may pose significant risks to populations, activities and properties [...]

The Methana volcano in Greece belongs to the western part of the Hellenic Volcanic Arc, where the African and Eurasian tectonic plates converge at a rate of approximately 3 cm/year. While volcanic hazard in Methana is considered low, the neotectonic basin constituting the Saronic Gulf area is seismically active and there is evidence of local geothermal activity. Monitoring is therefore crucial to characterize any activity at the volcano that could impact the local population. This study aims to detect surface deformation in the whole Methana peninsula based on a long stack of 99 Sentinel-1 C-band Synthetic Aperture Radar (SAR) images in interferometric wide swath mode acquired in March 2015–August 2019. A Multi-Temporal Interferometric SAR (MT-InSAR) processing approach is exploited using the Interferometric Point Target Analysis (IPTA) method, involving the extraction of a network of targets including both Persistent Scatterers (PS) and Distributed Scatterers (DS) to augment the monitoring capability across the varied land cover of the peninsula. Satellite geodetic data from 2006–2019 Global Positioning System (GPS) benchmark surveying are used to calibrate and validate the MT-InSAR results. Deformation monitoring records from permanent Global Navigation Satellite System (GNSS) stations, two of which were installed within the peninsula in 2004 (METH) and 2019 (MTNA), are also exploited for interpretation of the regional deformation scenario. Geological, topographic, and 2006–2019 seismological data enable better understanding of the ground deformation observed. Line-of-sight displacement velocities of the over 4700 PS and 6200 DS within the peninsula are from −18.1 to +7.5 mm/year. The MT-InSAR data suggest a complex displacement pattern across the volcano edifice, including local-scale land surface processes. In Methana town, ground stability is found on volcanoclasts and limestone for the majority of the urban area footprint while some deformation is observed in the suburban zones. At the Mavri Petra andesitic dome, time series of the exceptionally dense PS/DS network across blocks of agglomerate and cinder reveal seasonal fluctuation (5 mm amplitude) overlapping the long-term stable trend. Given the steepness of the slopes along the eastern flank of the volcano, displacement patterns may indicate mass movements. The GNSS, seismological and MT-InSAR analyses lead to a first account of deformation processes and their temporal evolution over the last years for Methana, thus providing initial information to feed into the volcano baseline hazard assessment and monitoring system.

The coastline environment and urban areas of Peru overlooking the Pacific Ocean are among the most affected by El Niño-Southern Oscillation (ENSO) events, and its cascading hazards such as floods, landslides and avalanches. In this work, the complete archives of the European Space Agency (ESA)’s European Remote-Sensing (ERS-1/2) missions and European Commission’s Copernicus Sentinel-1 constellation were screened to select synthetic aperture radar (SAR) images covering the most severe and recent ENSO-related flooding events that affected Lima, the capital and largest city of Peru, in 1997–1998 and 2017–2018. Based on SAR backscatter color composites and ratio maps retrieved from a series of pre-, cross- and post-event SAR pairs, flooded areas were delineated within the Rímac River watershed. These are mostly concentrated along the riverbanks and plain, where low-lying topography and gentle slopes (≤5°), together with the presence of alluvial deposits, also indicate greater susceptibility to flooding. A total of 409 areas (58.50 km2) revealing change were mapped, including 197 changes (32.10 km2) due to flooding-related backscatter variations (flooded areas, increased water flow in the riverbed, and riverbank collapses and damage), and 212 (26.40 km2) due to other processes (e.g., new urban developments, construction of river embankments, other engineering works, vegetation changes). Urban and landscape changes potentially contributing, either detrimentally or beneficially, to flooding susceptibility were identified and considered in the overall assessment of risk. The extent of built-up areas within the basin was mapped by combining information from the 2011 Global Urban Footprint (GUF) produced by the German Aerospace Center (DLR), the Open Street Map (OSM) accessed from the Quantum GIS (QGIS) service, and 2011–2019 very high-resolution optical imagery from Google Earth. The resulting flooding risk map highlights the sectors of potential concern along the Rímac River, should flooding events of equal severity as those captured by SAR images occur in the future.

Remote sensing data and methods are increasingly being embedded into assessments of volcanic processes and risk. This happens thanks to their capability to provide a spectrum of observation and measurement opportunities to accurately sense the dynamics, magnitude, frequency, and impacts of volcanic activity in the ultraviolet (UV), visible (VIS), infrared (IR), and microwave domains. Launched in mid-2018, the Special Issue “Remote Sensing of Volcanic Processes and Risk” of Remote Sensing gathers 19 research papers on the use of satellite, aerial, and ground-based remote sensing to detect thermal features and anomalies, investigate lava and pyroclastic flows, predict the flow path of lahars, measure gas emissions and plumes, and estimate ground deformation. The strong multi-disciplinary character of the approaches employed for volcano monitoring and the combination of a variety of sensor types, platforms, and methods that come out from the papers testify the current scientific and technology trends toward multi-data and multi-sensor monitoring solutions. The research advances presented in the published papers are achieved thanks to a wealth of data including but not limited to the following: thermal IR from satellite missions (e.g., MODIS, VIIRS, AVHRR, Landsat-8, Sentinel-2, ASTER, TET-1) and ground-based stations (e.g., FLIR cameras); digital elevation/surface models from airborne sensors (e.g., Light Detection And Ranging (LiDAR), or 3D laser scans) and satellite imagery (e.g., tri-stereo Pléiades, SPOT-6/7, PlanetScope); airborne hyperspectral surveys; geophysics (e.g., ground-penetrating radar, electromagnetic induction, magnetic survey); ground-based acoustic infrasound; ground-based scanning UV spectrometers; and ground-based and satellite Synthetic Aperture Radar (SAR) imaging (e.g., TerraSAR-X, Sentinel-1, Radarsat-2). Data processing approaches and methods include change detection, offset tracking, Interferometric SAR (InSAR), photogrammetry, hotspots and anomalies detection, neural networks, numerical modeling, inversion modeling, wavelet transforms, and image segmentation. Some authors also share codes for automated data analysis and demonstrate methods for post-processing standard products that are made available for end users, and which are expected to stimulate the research community to exploit them in other volcanological application contexts. The geographic breath is global, with case studies in Chile, Peru, Ecuador, Guatemala, Mexico, Hawai’i, Alaska, Kamchatka, Japan, Indonesia, Vanuatu, Réunion Island, Ethiopia, Canary Islands, Greece, Italy, and Iceland. The added value of the published research lies on the demonstration of the benefits that these remote sensing technologies have brought to knowledge of volcanoes that pose risk to local communities; back-analysis and critical revision of recent volcanic eruptions and unrest periods; and improvement of modeling and prediction methods. Therefore, this Special Issue provides not only a collection of forefront research in remote sensing applied to volcanology, but also a selection of case studies proving the societal impact that this scientific discipline can potentially generate on volcanic hazard and risk management.

Vehicle detection from satellite imagery can support different applications, such as security and situational awareness. In the civilian domain, it can provide quantitative evidence to investigate urban mobility and traffic patterns in cities. Satellite synthetic aperture radar (SAR) can help in detecting vehicles in (almost) all weather conditions and during the day and night. In this study, the capability of SAR StripMap imaging mode data to monitor traffic is analyzed using the case study of Wuhan, China. In ordinary times, the bridges crossing the Yangtze river are the key infrastructure allowing for urban mobility in Wuhan. More recently, the city has been the first in the world to be put in lockdown due to the outbreak of the Coronavirus Disease of 2019 (COVID-19). Using a very long time series of 294 COSMO-SkyMed StripMap HIMAGE mode scenes collected from 2011 to 2020, we detected vehicles on seven bridges, estimated their speed, and analyzed the traffic pattern over time. Vehicles are detected based on their azimuth shift caused by their across-track motion. Our goal is to monitor the variations in traffic instead of single-car detection. The results from 2011 to 2019 show a general increase in the number of vehicles crossing the bridges, as new infrastructure was built over the years. Variations in detected vehicle numbers were especially found during the two events of the 7th International Military Sports Council (CISM) Military World Games in October 2019, and the COVID-19 lockdown in early 2020. These events were therefore used for internal validation of our assessment of traffic patterns. On the other side, TomTom traffic index data were used for external validation. The results and their comparison with TomTom data prove the effectiveness of our method in detecting traffic patterns, but also demonstrate that mostly large vehicles (e.g., trucks or buses) are detected. Future work should be carried out to improve the detection rate of smaller vehicles.

This paper provides the first wide-area Interferometric Synthetic Aperture Radar (InSAR) survey of the whole eastern Trans-Mexican Volcanic Belt (42,200 km2). The aims are to identify ground deformation hotspots within major urbanized areas and rural valleys, establish baselines in geothermal exploration sites, and analyze deformation at geothermal exploitation sites and its relationship with energy production. The whole 2003−2010 ENVISAT C-band SAR archive available over the region was processed with the Small BAseline Subset (SBAS) InSAR method to retrieve over 840,000 coherent targets and estimate their ground displacement rates and time series. Land subsidence hotspots due to aquifer drawdown are found within the city of Puebla (up to −53 mm/year vertical rates, groundwater pumping for industrial use), Tlaxcala and Apizaco (−17 mm/year, industrial and public), the valley of Tecamachalco (−22 mm/year, agricultural), Tulancingo (−55 mm/year, public, industrial and agricultural), and in the eastern Mexico City metropolitan area (−44 mm/year, agricultural). The baseline for the Acoculco caldera complex shows widespread ground stability. Conversely, localized subsidence patterns of −5 to −10 mm/year exist around Las Derrumbadas and Cerro Pinto in the Serdán-Oriental basin, due to intense groundwater pumping for agriculture. A well-defined land subsidence area with −11 mm/year maximum rates is found at Los Humeros volcanic complex within Los Potreros collapse, correlating well with energy production infrastructure location and historical steam production rates. Field surveys carried out in Acoculco and Los Humeros in 2018 provide supporting evidence for the identification of hydrothermal manifestations, and understanding of the landscape and surface deformation patterns within the geothermal fields.

Illegal excavations in archaeological heritage sites (namely “looting”) are a global phenomenon. Satellite images are nowadays massively used by archaeologists to systematically document sites affected by looting. In parallel, remote sensing scientists are increasingly developing processing methods with a certain degree of automation to quantify looting using satellite imagery. To capture the state-of-the-art of this growing field of remote sensing, in this work 47 peer-reviewed research publications and grey literature are reviewed, accounting for: (i) the type of satellite data used, i.e., optical and synthetic aperture radar (SAR); (ii) properties of looting features utilized as proxies for damage assessment (e.g., shape, morphology, spectral signature); (iii) image processing workflows; and (iv) rationale for validation. Several scholars studied looting even prior to the conflicts recently affecting the Middle East and North Africa (MENA) region. Regardless of the method used for looting feature identification (either visual/manual, or with the aid of image processing), they preferred very high resolution (VHR) optical imagery, mainly black-and-white panchromatic, or pansharpened multispectral, whereas SAR is being used more recently by specialist image analysts only. Yet the full potential of VHR and high resolution (HR) multispectral information in optical imagery is to be exploited, with limited research studies testing spectral indices. To fill this gap, a range of looted sites across the MENA region are presented in this work, i.e., Lisht, Dashur, and Abusir el Malik (Egypt), and Tell Qarqur, Tell Jifar, Sergiopolis, Apamea, Dura Europos, and Tell Hizareen (Syria). The aim is to highlight: (i) the complementarity of HR multispectral data and VHR SAR with VHR optical imagery, (ii) usefulness of spectral profiles in the visible and near-infrared bands, and (iii) applicability of methods for multi-temporal change detection. Satellite data used for the demonstration include: HR multispectral imagery from the Copernicus Sentinel-2 constellation, VHR X-band SAR data from the COSMO-SkyMed mission, VHR panchromatic and multispectral WorldView-2 imagery, and further VHR optical data acquired by GeoEye-1, IKONOS-2, QuickBird-2, and WorldView-3, available through Google Earth. Commonalities between the different image processing methods are examined, alongside a critical discussion about automation in looting assessment, current lack of common practices in image processing, achievements in managing the uncertainty in looting feature interpretation, and current needs for more dissemination and user uptake. Directions toward sharing and harmonization of methodologies are outlined, and some proposals are made with regard to the aspects that the community working with satellite images should consider, in order to define best practices of satellite-based looting assessment.

We processed X-band COSMO-SkyMed 3-m resolution StripMap HIMAGE time series (March 2011–June 2013) with the Stanford Method for Persistent Scatterers (StaMPS), to retrieve an updated picture of the condition and structural health of the historic centre of Rome, Italy, and neighbouring quarters. Taking advantage of an average target density of over 2800 PS/km2, we analysed the spatial distribution of more than 310,000 radar targets against: (1) land cover; (2) the location of archaeological ruins and restoration activities; and (3) the size, orientation and morphology of historical buildings. Radar interpretation was addressed from the perspective of conservators, and the deformation estimates were correlated to local geohazards and triggering factors of structural collapse. In the context of overall stability, deformation was identified at the single-monument scale, e.g., for the Roman cistern and exedra in the Oppian Hill. Comparative assessment against InSAR processing of C-band imagery (1992–2010) published in the literature confirms the persistence of ground motions affecting monuments and subsidence in southern residential quarters adjacent to the Tiber River, due to the consolidation of compressible deposits. Vertical velocity estimated from COSMO-SkyMed PS exceeds −7.0 mm/y in areas of recent urbanization.

Synthetic aperture radar (SAR) imagery has long been used in archaeology since the earliest space radar missions in the 1980s. In the current scenario of SAR missions, the Italian Space Agency (ASI)’s COnstellation of small Satellites for Mediterranean basin Observation (COSMO-SkyMed) has peculiar properties that make this mission of potential use by archaeologists and heritage practitioners: high to very high spatial resolution, site revisit of up to one day, and conspicuous image archives over cultural heritage sites across the globe. While recent literature and the number of research projects using COSMO-SkyMed data for science and applied research suggest a growing interest in these data, it is felt that COSMO-SkyMed still needs to be further disseminated across the archaeological remote sensing community. This paper therefore offers a portfolio of use-cases that were developed in the last two years in the Scientific Research Unit of ASI, where COSMO-SkyMed data were analysed to study and monitor cultural landscapes and heritage sites. SAR-based applications in archaeological and cultural heritage sites in Peru, Syria, Italy, and Iraq, provide evidence on how subsurface and buried features can be detected by interpreting SAR backscatter, its spatial and temporal changes, and interferometric coherence, and how SAR-derived digital elevation models (DEM) can be used to survey surface archaeological features. The use-cases also showcase how high temporal revisit SAR time series can support environmental monitoring of land surface processes, and condition assessment of archaeological heritage and landscape disturbance due to anthropogenic impact (e.g., agriculture, mining, looting). For the first time, this paper provides an overview of the capabilities of COSMO-SkyMed imagery in StripMap Himage and Spotlight-2 mode to support archaeological studies, with the aim to encourage remote sensing scientists and archaeologists to search for and exploit these data for their investigations and research activities. Furthermore, some considerations are made with regard to the perspectives opened by the upcoming launch of ASI’s COSMO-SkyMed Second Generation constellation.

Geoinformation derived from Earth observation (EO) plays a key role for detecting, analyzing and monitoring landslides to assist hazard and risk analysis. Within the framework of the EC-GMES-FP7 project SAFER (Services and Applications For Emergency Response) a semi-automated object-based approach for landslide detection and classification has been developed. The method was applied to a case study in North-Western Italy using SPOT-5 imagery and a digital elevation model (DEM), including its derivatives slope, aspect, curvature and plan curvature. For the classification in the object-based environment spectral, spatial and morphological properties as well as context information were used. In a first step, landslides were classified on a coarse segmentation level to separate them from other features with similar spectral characteristics. Thereafter, the classification was refined on a finer segmentation level, where two categories of mass movements were differentiated: flow-like landslides and other landslide types. In total, an area of 3.77 km² was detected as landslide-affected area, 1.68 km² were classified as flow-like landslides and 2.09 km² as other landslide types. The outcomes were compared to and validated by pre-existing landslide inventory data (IFFI and PAI) and an interpretation of PSI (Persistent Scatterer Interferometry) measures derived from ERS1/2, ENVISAT ASAR and RADARSAT-1 data. The spatial overlap of the detected landslides and existing landslide inventories revealed 44.8% (IFFI) and 50.4% (PAI), respectively. About 32% of the polygons identified through OBIA are covered by persistent scatterers data.

Very high-resolution (VHR) optical satellite imagery (≤5 m) is nowadays an established source of information to monitor cultural and archaeological heritage that is exposed to hazards and anthropogenic threats to their conservation, whereas few publications specifically investigate the role that regularly acquired images from high-resolution (HR) satellite sensors (5–30 m) may play in this application domain. This paper aims to appraise the potential of the multispectral constellation Sentinel-2 of the European Commission Earth observation programme Copernicus to detect prominent features and changes in heritage sites, during both ordinary times and crisis. We test the 10 m spatial resolution of the 3 visible spectral bands of Sentinel-2 for substantiation of single local events—that is, wall collapses in the UNESCO World Heritage site of the Old City of Aleppo (Syria)—and for hotspot mapping of recurrent incidents—that is, the archaeological looting in the archaeological site of Apamea (Syria). By screening long Sentinel-2 time series consisting of 114 images for Aleppo and 57 images for Apamea, we demonstrate that changes of textural properties and surface reflectance can be logged accurately in time and space and can be associated to events relevant for conservation. VHR imagery from Google Earth was used for the validation and identification of trends occurring prior to the Sentinel-2 launch. We also demonstrate how to exploit the Sentinel-2 short revisiting time (5 days) and large swath (290 km) for multi-temporal tracking of spatial patterns of urban sprawl across the cultural landscape of the World Heritage Site of Cyrene (Libya), and the three coastal ancient Greek sites of Tocra, Ptolemais, and Apollonia in Cyrenaica. With the future development of tailored machine learning approaches of feature extraction and pattern detection, Sentinel-2 can become extremely useful to screen wider regions with short revisiting times and to undertake comparative condition assessment analyses of different heritage sites.

The “Lines and Geoglyphs of Nasca and Palpa” in Peru are among the most well-known UNESCO World Heritage Sites globally, and an exemplar of site where heritage assets cannot be separated from their natural and anthropogenic environment. The site is exposed to interactions with natural processes, as well as human presence. In this work, 3-m resolution synthetic aperture radar (SAR) StripMap HIMAGE HH-polarised scenes acquired by the X-band COSMO-SkyMed constellation are exploited to track two events of human-induced landscape disturbance that occurred in December 2014 and January 2018. Pre-, cross-, and post-event interferometric SAR (InSAR) pairs characterised by small temporal and normal baselines allow the detection of temporal decorrelation associated with the two events, the extent and time reference of which match with online photographic and video evidence, published literature, web news, and press releases by the Ministry of Culture in Peru. Further elements enhancing the understanding of the 2018 event come from 10-m resolution Sentinel-2B satellite data that reveal the occurrence of apparent changes of surface reflectance due to uncovering of the light grey-yellow clay underneath the darker pebble constituting the fragile surface of the Pampa de Jumana. This scientific study confirms that SAR imagery archives, such as those being built by COSMO-SkyMed for Nasca, prove valuable for the retrospective analysis and digital recording of human-induced landscape disturbance events from space. These archives therefore act as essential sources of geospatial information on the conservation history of heritage sites and assets.

With a wide spectrum of imaging capabilities—from optical to radar sensors, low to very high resolution, continental to local scale, single-image to multi-temporal approaches, yearly to sub-daily acquisition repeat cycles—Earth Observation (EO) offers several opportunities for the geoscience community to map and monitor natural and human-induced Earth hazards from space. The Special Issue “Observing Geohazards from Space” of Geosciences gathers 12 research articles on the development, validation, and implementation of satellite EO data, processing methods, and applications for mapping and monitoring of geohazards such as slow moving landslides, ground subsidence and uplift, and active and abandoned mining-induced ground movements. Papers published in this Special Issue provide novel case studies demonstrating how EO and remote sensing data can be used to detect and delineate land instability and geological hazards in different environmental contexts and using a range of spatial resolutions and image processing methods. Remote sensing datasets used in the Special Issue papers encompass satellite imagery from the ERS-1/2, ENVISAT, RADARSAT-1/2, and Sentinel-1 C-band, TerraSAR-X and COSMO-SkyMed X-band, and ALOS L-band SAR missions; Landsat 7, SPOT-5, WorldView-2/3, and Sentinel-2 multi-spectral data; UAV-derived RGB and near infrared aerial photographs; LiDAR surveying; and GNSS positioning data. Techniques that are showcased include, but are not limited to, differential Interferometric SAR (InSAR) and its advanced approaches such as Persistent Scatterers (PS) and Small Baseline Subset (SBAS) methods to estimate ground deformation, Object-Based Image Analysis (OBIA) to identify landslides in high resolution multi-spectral data, UAV and airborne photogrammetry, Structure-from-Motion (SfM) for digital elevation model generation, aerial photo-interpretation, feature extraction, and time series analysis. Case studies presented in the papers focus on landslides, natural and human-induced subsidence, and groundwater management and mining-related ground deformation in many local to regional-scale study areas in Austria, Belgium, Italy, Slovakia, Spain, and the UK.

In this paper, we investigate land motion and groundwater level change phenomena using differential interferometric synthetic aperture radar (DInSAR) over the Northumberland and Durham coalfield in the United Kingdom. The study re-visits earlier research that applied a persistent scatterers interferometry (PSI) technique to ERS (European Remote Sensing) and ENVISAT (Environmental Satellite) data. Here, the Intermittent Small Baseline Subset (ISBAS) DInSAR technique is applied to ERS, ENVISAT and Sentinel-1 SAR datasets covering the late 1990s, the 2000s and the mid-2010s, respectively, to increase spatial coverage, aid the geological interpretation and consider the latest Sentinel-1 data. The ERS data identify surface depressions in proximity to former collieries, while all three data sets ascertain broad areas are experiencing regional scale uplift, often occurring in previously mined areas. Uplift is attributed to increases in pore pressure in the overburden following the cessation of groundwater pumping after mine closure. Rising groundwater levels are found to correlate to ground motion measurements at selected monitoring sites, most notably in the surrounding area of Ashington. The area is divided by an impermeable EW fault; to the south, surface heave was identified as groundwater levels rose in the 1990s, whereas to the north, this phenomenon occurred two decades later in the 2010s. The data emphasize the complexity of the post-mining surface and subsurface environment and highlight the benefit that InSAR, utilizing the ISBAS technique, can provide in its characterization.

Satellite Interferometric Synthetic Aperture Radar (InSAR), geological data and Small Unmanned Aerial Vehicle (SUAV) surveying was used to enhance our understanding of ground movement at five areas of interest in Northern Ireland. In total 68 ERS-1/2 images 1992–2000 were processed with the Small Baseline Subset (SBAS) InSAR technique to derive the baseline ground instability scenario of key areas of interest for five stakeholders: TransportNI, Northern Ireland Railways, Department for the Economy, Arup, and Belfast City Council. These stakeholders require monitoring of ground deformation across either their geotechnical infrastructure (i.e., embankments, cuttings, engineered fills and earth retaining structures) or assessment of subsidence risk as a result of abandoned mine workings, using the most efficient, cost-effective methods, with a view to minimising and managing risk to their businesses. The InSAR results provided an overview of the extent and magnitude of ground deformation for a 3000 km2 region, including the key sites of the disused salt mines in Carrickfergus, the Belfast–Bangor railway line, Throne Bend and Ligoniel Park in Belfast, Straidkilly and Garron Point along the Antrim Coast Road, plus other urbanised areas in and around Belfast. Tailored SUAV campaigns with a X8 airframe and generation of very high resolution ortho-photographs and a 3D surface model via the Structure from Motion (SfM) approach at Maiden Mount salt mine collapse in Carrickfergus in 2016 and 2017 also demonstrate the benefits of very high resolution surveying technologies to detect localised deformation and indicators of ground instability.

Vajedian et al. [1] present an improved method for the derivation of deformation parameters using satellite Interferometric Synthetic Aperture Radar (InSAR) data. The method is a modification of the Small Baseline Subset (SBAS) method as implemented in the StaMPS (Stanford Method for Persistent Scatterers) software. The modification includes many steps including the filtering of the differential interferograms, integration with GPS data and advanced phase unwrapping “to overcome a lot of short- and long-wavelength artifacts that are clearly visible in StaMPS results” (cf. [1], p. 8331). The authors refer to this new approach as the Improved SBAS, or ISBAS, method. [...]