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.
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.
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.
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.
Abstract
The instability situation affecting the Middle East poses threats to preservation of cultural heritage. Mapping efforts based on satellite imagery currently concentrate more on recording ...human-induced damage than impacts of unforeseen natural events (e.g. floods). In 2018, the UNESCO World Heritage Site of Bosra and the landscape of Sergiopolis-Resafa in Syria were flooded after heavy rainfall. While the first incident was reported by heritage organisations (although information was limited to the main monument), the second event was completely unknown. Using optical and radar satellite images from the European Commission’s Copernicus Sentinels fleet and the Italian Space Agency’s COSMO-SkyMed constellation, we prove that these data are an enormous reservoir of information to assess comprehensively the duration, extent and severity of such natural events. In Bosra, several key assets were flooded besides the Roman Theatre, with waters taking from few days to several weeks to evacuate and dry out. In Sergiopolis, while the main ruins were sheltered by the fortification walls, the nearby floodplain was inundated. The floodwater-flow pattern resembled the simulations developed by archaeologists to prove the existence of an ancient system of embankments and dam. Our results suggest that many other events posing risk to heritage assets, otherwise unnoticed, may be unveiled if current satellite imagery archives, yet to analyse, are systematically screened.
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.
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.
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.
“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).
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.