1. Introduction
Optical satellite images are nowadays established Earth observation data for research teams and practitioners to monitor cultural and archaeological heritage that is exposed to hazards and anthropogenic threats to their conservation. If clouds do not hamper the visibility of the area of interest, the features on the ground are usually clearly seen and the surface changes are detected. Thanks to these properties, optical remote sensing is regarded by the archaeological and heritage community as an objective source of information allowing for a conservative estimate of the condition on the ground [
1]. However, the expertise of the individuals, their level of training in mapping and interpreting features, as well as a priori knowledge and expectations on what to see, are all factors influencing the level of accuracy that can be achieved in the satellite-based assessment.
In the last years, optical remote sensing has proved extremely valuable in reducing the uncertainty about the condition of heritage assets across the Middle East and North Africa (MENA) region and, more generally, in the Mediterranean countries. Several initiatives at both the international and national level are ongoing. For example, the UNESCO Emergency Safeguarding of the Syrian Cultural Heritage project [
2], the American School of Oriental Research Cultural Heritage Initiative (ASOR CHI) [
3], the Association for the Protection of Syrian Archaeology (APSA2011) [
4], and the Endangered Archaeology in the Middle East & North Africa (EAMENA) [
5]. These projects collect information from satellite imagery not only to map the extent of damage and destruction incurred by built and movable heritage during crises, but also to assess the impact of modern human activities in ordinary times, such as agricultural exploitation of land, urban expansion, and infrastructure construction.
As recently recalled by Casana and Laugier [
6], the above-mentioned projects demonstrate that the research and user communities are moving from site-based studies and single-incident substantiations to more systematic, regional-scale efforts over larger regions. However, such large spatial coverage and monitoring of a multitude of sites at the same time still pose operational challenges; first of all, the need to access large volumes of timely and updated satellite imagery at a suitable resolution. Satellite data uploaded onto open visualization platforms (for example, Google Earth and Bing maps) have partly contributed to solving this issue and are now Earth observation facilities incorporated into the image interpretation methodologies of some of these heritage initiatives (for example, Reference [
5]).
The optical images used for such studies mostly come from satellites operated by commercial companies (for example, WorldView and GeoEye) because of their sub-meter spatial resolution. For instance, the potential for looting detection has been recently discussed in Reference [
7]. To avoid confusion on terminology, in this work, we adopt the classification by which optical satellite imagery is referred to as ‘very high resolution’ (VHR) if characterised by a spatial resolution of 5 m to even less than 1 m, and ‘high resolution’ (HR) if the spatial resolution is between 5 and 30 m; see, for example References [
8,
9]).
From the point of view of research development, the emphasis on the exploitation of VHR images—mostly driven by the need for searching for very small features—has led to very few publications that specifically investigate the role that regularly acquired HR images from constellations such as the Landsat missions may play in this application domain. It is also to be noted that most of these publications rather concentrate on land use impacts in the surroundings of the sites and/or their buffer zones and conservation areas. These often exploit broad classification schemes for their assessment (for example, urban versus agricultural [
10,
11]), but seldom focus on the recognition of specific man-made features and tracking of the temporal evolution of their patterns.
On the other side, the Sentinel-2 multispectral constellation of the European Commission programme Copernicus [
12,
13,
14] is an Earth observation system that is increasingly gathering attention in archaeological remote sensing. The thirteen spectral bands (443–2190 nm) and HR imaging capabilities in visible and near-infrared bands at 10 m spatial resolution have been already tested for archaeological prospection (for example, References [
15,
16]). On the contrary, it is still to be assessed how Sentinel-2 can be used systematically to detect prominent features and changes in heritage sites during both ordinary times and crisis. To the best of our knowledge, no research has been published yet to test Sentinel-2 for this specific capability.
This is, therefore, the aim and overarching research question of this paper. In particular, we explore the value of Sentinel-2 time series for multi-temporal HR assessment for event verification and hotspot mapping. Our appraisal is structured around a selection of case studies in Syria and Libya that we investigated using Sentinel-2 and Google Earth time series (
Section 2) to examine two scenarios (
Section 3):
The local heritage has been impacted by event(s) known based on precise background information (that is, including the location, time of occurrence, and typology of the event). Therefore the features to search for in the satellite images can be anticipated, are very distinctive, and easy to interpret in the cloud- and fog-free visibility conditions (
Section 3.1);
The heritage site is known to be a hotspot of specific threats to conservation (for example, intense and repeated looting) and the history of incidents may suggest that further events can happen, and their manifestation can be captured even at HR (
Section 3.2).
The results of the Sentinel-2 time series analysis were validated based on integration with VHR imagery—that is, DigitalGlobe and CNES/Airbus images available in Google Earth—and background knowledge from published literature and incident reports. When no Sentinel-2 data were available before a certain date of interest, we used the pre-existing Google Earth time series to identify the trends and dynamics of the investigated studies for further verification in the more updated Sentinel-2 time series.
This was the case in the demonstration site of Cyrene where there was an opportunity to investigate an anthropogenic process of a potential threat to heritage conservation that manifests at the surface through distinctive features changing in time. We, therefore, discuss how systematic Sentinel-2 acquisitions taken with short revisiting time can be used for the dynamic detection of such features (
Section 3.3). We then test the regional scalability of the method by extending the analysis across the entire Libyan region of Cyrenaica with a focus on coastal heritage sites affected by urban sprawling (
Section 4). We finally outline future research directions in this field accounting for the emerging opportunities offered by machine learning (
Section 5).
2. Materials and Methods
To carry out our appraisal, we tested the Sentinel-2 HR imagery on the following heritage sites in Syria and Libya (
Figure 1):
the UNESCO World Heritage Site (WHS) of the Ancient City of Aleppo in Syria (centre coordinates: 36°13′59″N; 37°10′00″E);
the Hellenistic site of Apamea in Syria (35°25′11″N; 36°24′05″E), proposed for inscription on the UNESCO World Heritage List in 1999;
the cultural landscape of Cyrene WHS and the modern town of Shahat in Libya (32°48′12″N; 21°51′46″E);
the ancient Greek towns of Apollonia (modern Susah; 32°54′07″N; 21°58′11″E), Ptolemais (modern Ad Dirsiyah or Tolmeita; 32°42′24″N; 20°57′10″E), and Tocra (32°32′20″N; 20°34′10″E) in northern Cyrenaica, eastern Libya.
Table 1 summarizes the heritage category, status and type of damage, site disturbance or land cover–land use change processes affecting the above demonstration sites. The background knowledge about the processes under investigation is based on published and grey literature, outcomes of previous satellite-based assessments, and information extracted from the news, reports, photos, and videos found in the web and social media. Specific citations are reported throughout the text and summarized in
Table 1. Our selection of demonstration sites aims to cover processes of potential threat to heritage conservation that occur in ordinary times as well as during times of crisis. Sentinel-2 data were analysed against different sizes, spatial extents, severities, and temporal evolutions of human activities, as detailed in
Section 3.
Our motivation to test the Sentinel-2 constellation for the condition assessment of heritage sites relies on the following image properties:
- (i)
Systematic global acquisition of multispectral imagery with consistent imaging parameters;
- (ii)
Spatial resolution of 10 m in the visible bands [
12,
13];
- (iii)
Short revisiting time (up to 5 days) over Europe and the Mediterranean countries [
24].
We used the Copernicus Open Access Hub to access the Sentinel-2 products acquired by the Multispectral Instrument (MSI) onboard the twin satellites Sentinel-2A and Sentinel-2B. This is the dedicated platform where Top-Of-Atmosphere (TOA) reflectance in cartographic geometry (that is, Level 1C, or L1C) products are regularly published by the European Space Agency (ESA) very soon after acquisition. Since April 2017, the Bottom Of Atmosphere (BOA) reflectance in cartographic geometry (that is, Level 2A, or L2A) products are also published 48–60 h after the L1C products, based on automated processing of the latter with the Sen2Cor processor and PlanetDEM Digital Elevation Model [
25]. For each demonstration site, the temporal period of our observation spanned from the first image available in 2015 up to December 2017.
Of the 13 MSI spectral bands, we focused on the B2 (490 nm), B3 (560 nm), B4 (665 nm) bands and their Red-Green-Blue (RGB) colour composites at full resolution. We intentionally concentrated on these bands because the majority of archaeologists and practitioners assessing the condition of heritage sites based on optical satellite imagery mostly rely on the use of RGB information of VHR satellite images (for example, References [
1,
3,
5,
6,
19]). A proof of concept on the visible bands, without further sophisticated processing of these and other bands, is expected to ease the initial uptake of Sentinel-2 by a wider spectrum of practitioners. We acknowledge though, that future research and implementation could take advantage of recently developed methods. For example, the sharpening of Sentinel-2 bands with spatial resolutions of 20 m and 60 m to a spatial resolution of 10 m, in order to obtain a full set of visible and near-infrared (VNIR) and shortwave infrared (SWIR) spectral bands with a spatial resolution of 10 m [
26].
To carry out our analysis we discarded products where clouds hampered the visibility of the heritage sites.
Since Sentinel-2 is an HR satellite constellation, we did not expect to distinguish small to very small features (that is, <10 m). Instead, we aimed to demonstrate that Sentinel-2 data can be used for the systematic and regular screening at the site scale to identify larger single features, clusters of small features, and accumulations of progressive (small) changes that eventually manifest as new prominent surface features. For example, while individual looting features could not be resolved in the Sentinel-2 data, clusters of several adjacent looting holes are clearly detectable (see
Section 3.2).
We ran a multi-temporal analysis by using the full Sentinel-2 archive available for each of the demonstration sites that provides, at best, a 5-day revisiting time from one acquisition to the next. However, because Sentinel-2A was launched on 23 June 2015 and in consideration of the time needed for commissioning (that is, ~3 months from the end of launch and early orbit phase), in some cases, the earlier cloud-free images available for the demonstration sites were acquired in October 2015.
We spatially analysed the Sentinel-2 stacks against ancillary data including, but not limited to site plans, maps of heritage assets, and thematic maps. In particular, for each demonstration site, we defined the T
0 time reference, against which we detected the new features based on either their textural or surface reflectance properties. For this purpose, we applied the state-of-the-art change detection techniques currently used across the specialist community (for example, Reference [
27]). In practice, for each pair of Sentinel-2 images acquired at T
n−1 and T
n, the surface reflectance at bands B4, B3, B2 and the textural properties were compared and clusters of pixels indicating hotspots of change were highlighted and zoned through a conversion to polygons.
We interpreted the observed features following the same standard mapping methodologies that are commonly used by research teams and practitioners to assess the condition of a heritage site using optical satellite imagery [
3,
5,
6,
19]. To label the features and associate them to specific events of relevance for the site conservation, we browsed incident reports that are regularly issued by local heritage bodies, international organizations and initiatives (for example, ASOR, Syrian Directorate-General for Antiquities and Museums—DGAM, and Geospatial Technologies and Human Rights Project of the American Association for the Advancement of Science—AAAS), as well as the news in broadcasts and social media.
For validation and trend analysis, we accessed open-source VHR images through Google Earth (see
Appendix A). These data provided the sub-meter spatial resolution and allowed us to upgrade the spatial definition of the identified features. However, their temporal coverage of the demonstration sites was discontinuous and not regular. In some circumstances, this implied the existence of temporal shifts with respect to the date of the reported incidents, as well as to the more regular temporal coverage provided by the Sentinel-2 time series. On the other hand, since Google Earth VHR images were acquired before June 2015, they were extremely useful to define the T
0 time reference, that is, the site condition prior to the first available Sentinel-2 image. The above-mentioned limitations are well discussed in the literature not only compared with the regular provision of VHR optical data through direct purchase, funds, sponsorships and agreements [
6], but also with the regularly acquired time series of satellite imagery acquired by space-borne sensors operating with different wavelengths (for example, the Synthetic Aperture Radar—SAR [
20]). We accounted for these advantages and limitations brought by VHR imagery when we interpreted them in combination with the Sentinel-2 time series (
Section 3 and
Section 4).
4. Discussion
The demonstration on the 2016 wall collapse in the Citadel of Aleppo WHS is an example of systematic examination of the regularly acquired optical satellite imagery, even at HR, that can provide objective evidence to detect impacts of damaging events on the conservation of local heritage assets, and substantiate incident reports and in situ observations. It is, of course, understood that we were able to capture this damaging event through the systematic monitoring with Sentinel-2 against the baseline reference of the Citadel of Aleppo because the collapse feature was clearly visible to the satellite, and its size and extent were large enough to be resolved at a 10 m spatial resolution.
It can be argued that, under different conditions, no feature could have been detected with Sentinel-2 due to its spatial resolution. In this regard, we recall the previous collapse that occurred along the eastern section of the Citadel walls on 11 July 2015. This event was reported and documented in the broadcast media the following day [
32] and also by ASOR [
33]. Based on published photographs taken from outside the Citadel, it can be observed that the collapse affected a few metres of the walls and the rubble accumulated in the upper part of the slope, just at the toe of the former walls. Therefore, this collapse generated a smaller event footprint than the event in August 2016.
Google Earth does not provide an image taken on the day or soon after the 2015 collapse happened (see
Table A1 in
Appendix A). We can only compare two DigitalGlobe images acquired on 15 December 2014 and 26 October 2015 (
Figure 9a,b) to analyse the collapse at a VHR. Although the visibility and illumination conditions are not ideal, the 26 October 2015 image somehow allows for the identification of the area affected by the collapse that occurred in July (
Figure 9b). A concave feature is also visible at the toe of the collapsed wall section that was not visible previously.
The comparison between the Sentinel-2 images acquired on 8 July 2015 and 17 August 2015 suggests that the collapse may not have been detected with Sentinel-2 without a priori knowledge. However, if the lower resolution of Sentinel-2 was a limiting factor for the event detection, subsequent Sentinel-2 images proved the benefit of having systematic acquisitions to follow the cascading effects of the wall collapse on the condition of the tell slope. By screening the Sentinel-2 dataset, it can be observed that the condition of the eastern slope of the tell does not change until 16 February 2016 when a feature resembling the track and toe deposit of collapsed debris becomes distinctively and permanently visible (
Figure 9c). This suggests that, in the absence of remediation works and due to the exposure to weathering, the slope further eroded and some rubble moved further down the slope. DigitalGlobe images acquired on 20 March 2016 and 23 March 2016 validate our Sentinel-2-based assessment (
Figure 9d), thus, implicitly providing proof that Sentinel-2 is valuable to track relevant changes in the heritage site beyond the detection of a single event.
Multi-temporal tracking is, undoubtedly, the key functionality that Sentinel-2 provides for hotspot monitoring, thanks to the combined high revisiting time of the two satellites. Our test on Apamea not only contributes to updating the information on the condition of this archaeological site as the last records date back to June 2016 [
26], but more importantly, it demonstrates that the Sentinel-2 constellation allows for monitoring looting and similar anthropogenic activities at this scale. While single small holes can be resolved only by means of VHR satellite imagery, the experiment on Apamea suggests that we do not need exclusively VHR images to track new and expanding looting with temporal accuracy. Far from stating that Sentinel-2 or other HR optical space missions should and could be used as a replacement for VHR Earth observation solutions, we instead propose the joint use of these systems making the best of their respective advantages (that is, freely accessible regular acquisitions on one side; higher spatial resolution on the other).
In Cyrene, we also proved that dense cloud-free Sentinel-2 time series can allow changes affecting cultural landscapes to be monitored and mapped with similar spatial accuracy and temporal granularity that can be achieved with an ad hoc VHR time series. Cyrene was a fortunate case study owing to the wealth of VHR images available via Google Earth. However, in cases in which the temporal coverage is discontinuous or cases in which such a dense VHR is not available, regularly acquired Sentinel-2 images become an asset for WHS condition monitoring.
While supervised classification with broad classes (for example, ‘urban’ versus ‘undeveloped land’) is an agile mapping approach, our test of feature extraction proves that an insightful analysis can be carried out with Sentinel-2 data to geospatially assess the changing degree of exposure of vulnerable heritage. With this method, it was clear that not all the areas of the cultural landscape surrounding the WHS of Cyrene were vulnerable at the same time. Prior to 2014, the urban sprawl was predominant east of Shahat and could potentially impact only the east necropolis. Starting from 2014 it appears that the urban sprawl became an increasing threat also for the conservation of the south necropolis, west of the modern town. Our experiment suggests that, if a dynamic assessment was undertaken in near real-time as new Sentinel-2 images were ingested in the time series analysis, the satellite-based condition assessment may inform the decision-making process and help to monitor the compliance of local land use changes with planning and landscape policies and regulations.
As found in Cyrene, the distinctive features relating to land cover and land use changes may be local manifestations of more regional phenomena. These situations can be captured over wide regions by exploiting the large swath of single scenes or mosaics of multiple (but coeval) Sentinel-2 images.
To simulate how Sentinel-2 could be used for this purpose, we ran the same test to extract and map block features due to the urban sprawl across Cyrenaica in Libya. We analysed a multi-temporal mosaic that we made by combining Sentinel-2 images available over Cyrenaica, from 29 July 2015 to 4 October 2015 (
Figure 10a). Our analysis suggests that in 2015 a number of block features were already widespread and common features in all the major urban settlements of Cyrenaica. The comparison with the Google Earth VHR time series pre-dating the Sentinel-2 mosaic confirms that the number and density of these block patterns started increasing since 2013, according to a general trend that can be observed in many other locations.
By intersecting the block features extracted from the Sentinel-2 mosaic and the location of known archaeological sites, three hotspots are found in the ancient Greek towns of Apollonia, Ptolemais, and Tocra (
Figure 10a). These are coastal sites embedded in the urban environments of the modern towns that have historically developed nearby [
23]: Susah, Ad Dirsiyah, and modern Tocra, respectively (
Figure 10b–d).
The synoptic view achieved with only two Sentinel-2 images allowed us to highlight slightly different situations in these three sites. In 2015, in Tocra and Apollonia, recent urbanisation was not immediately affecting the archaeological sites, although this was already an anthropogenic process that was apparently re-shaping the wider landscape (
Figure 10b,d). In particular, Susah expanded through one large project of urban development in the south-east and through new roads and building blocks in the south-west, but no new block features were developed at that time in the proximity of the archaeological site of Apollonia (
Figure 10d). However, some properties were already built very close to the southern boundary of the site (see also
Figure 11d). Similarly, there was no evidence of new blocks built in the proximity of the ruins of Tocra in 2015, but outside the modern town, several road-blocks had already appeared (
Figure 10b). Conversely, in mid-2015 urbanisation was already affecting the archaeological site of Ptolemais (
Figure 10c). A new road-block enclosing land lots had been already created close to the western boundary of Ptolemais, in a portion of land west of the Wadi Khambis that was previously vacant (see also
Figure 11b). Other road-blocks had been laid down along the coast, north-east of the site, beyond the Wadi Ziwana (
Figure 10c).
By analysing Sentinel-2 data until the end of 2017, no further changes were found within the footprints or along the boundaries of the archaeological sites of Tocra and Ptolemais (
Figure 11a–c). In particular, the urban sprawl continued south and east of the modern town of Tocra (
Figure 11a), so only continuous monitoring could show the future direct impacts on heritage conservation. In Ptolemais, the block feature that was observed in 2015 (
Figure 11b) did not appear to have been populated yet with buildings or farming activities in 2017 (
Figure 11c), so it had still not developed into a mature stage according to the cycle discussed in
Section 3.3 (see also
Figure 5). Conversely, we detected clear changes in the surface reflectance along the southern boundary of Apollonia that matched with the development of previously vacant land lots from late 2015 to mid-2017 (
Figure 11d,e). Therefore, we could update our assessment of the degree of exposure of these three sites to the impacts due to urban expansion and land use changes.
5. Conclusions
In this paper, we prove that, at the right spatial scale and for features that are a few meters in size, Sentinel-2 can be successfully applied for the condition assessment of heritage sites. The results of our research suggest that the general idea that ‘high resolution’ Earth observation sensors are of limited use compared to sensors operating at ‘very high resolution’ should be reconsidered.
We tested Sentinel-2 on both single localised events (for example, wall collapses) and sequences of recurrent incidents (for example, illegal excavations of archaeological sites), also searching for features of different size, spanning from clusters of meter-scale features (for example, looting holes) to single medium-size urban features (for example, road and building blocks, with dimensions of a few tens of meters). Despite the obvious constraint in the spatial resolution that we acknowledge in the discussion of the results (see
Section 3 and
Section 4), Sentinel-2 allowed for the detection of events of potential relevance for heritage conservation in all the demonstration sites investigated in this paper.
In renowned hotspots, the availability of prior knowledge of the processes to expect and of the baseline condition of the local heritage increases the likelihood for Sentinel-2 to capture textural anomalies or changes of surface reflectance properties that could relate to site disturbances or incidents of potential damage. In this regard, the case studies of Aleppo and Apamea (Syria) provide a clear demonstration.
However, to investigate processes such as urban development that manifest through spatial and temporal patterns of distinctive features, approaches of feature extraction making the best use of the two properties that make Sentinel-2 quite advantageous compared to VHR satellite imagery—the very short revising time between two consecutive acquisitions and the large swath width—are promising. These conditions cannot always be achieved using VHR data, as acknowledged in the literature [
6]. Sentinel-2 could contribute to addressing this limitation in the framework of a multi-scale condition assessment of heritage sites. For instance, a multi-temporal HR satellite-based condition assessment could provide information to schedule a targeted survey with ad hoc VHR imagery.
The test that we ran on the block features associated with urban sprawl in Cyrenaica (Libya) suggests that Sentinel-2 can well become an extremely useful source of information to screen wider regions and undertake comparative analyses of the condition of different heritage sites.
The approach that we demonstrate in this paper with Sentinel-2 may be applied to other satellite missions, such as Landsat-8. However, it is worth mentioning that the spatial resolution of Landsat-8 is coarser than that of Sentinel-2 (15 m for the panchromatic band and 30 m for the visible to short-wave infrared channels), the revisiting time longer (16 days) and the swath smaller (185 km). Therefore, Sentinel-2 appears to be a good candidate among the current freely accessible space solutions to bridge the gap between HR and VHR optical imagery for the condition assessment of heritage sites.
We can easily anticipate that the implementation of automatic processing chains searching for and extracting features through long time series would be a technological development accelerator towards the systematic use of Sentinel-2 for the condition assessment and monitoring of local to wide areas of study. This automation should provide a complement to the analyst-driven methods that currently represent the state-of-the-art methods in this field [
6] and, ideally, help to mitigate the drawbacks of manual examination. The latter is still a common practice across the archaeological and heritage community. However, although this approach has been proved to lead to reliable and successful mapping results, the community has started to question it from a methodological point of view. The reason is that, at equal conditions of feature visibility and imaging definition, the mapping accuracy achievable strongly depends on the operators’ skills and expertise. Casana and Laugier [
6] have recently reviewed this aspect specifically based on published exercises.
However, even in situations when operators are fully trained, manual examination remains a time-consuming task and is therefore not sustainable over large areas or for analysis of long time series. It is beyond the scope of this paper to demonstrate how automated processing for the detection and extraction of features should be designed and deployed. It is instead worth spending a few words on the possible technical requirements and recommendations that future research could take into account to develop such automated processing tools. In particular, the preparation of training data that are required to develop an effective machine learning approach appears to be a critical step. The experience gained by the community on different geographic locations is now being translated into reference catalogues of disturbance/threat/damage types and associated databases of documented incidents (see the example presented in [
5]). These are repositories from which it is possible to derive categories of simple, distinctive and frequently ubiquitous features and to define, for example, the labelled training data necessary to instruct supervised learning for the automatic recognition of patterns due to a certain process relevant for heritage conservation. Training data that are built according to the above-mentioned criteria (and specifically the distinctiveness of the pattern from other unrelated objects) should help to reduce the occurrence of false positive errors and mismatching. Recent publications presenting the first attempts of automation on feature detection in VHR multi-spectral satellite images have addressed this technical requirement by proposing algorithms searching for repeated featural motifs and building upon existing methods of supervised classification and unsupervised localization [
34]. Although these methods have been applied at VHR level, they outline one of the possible directions along which to develop methods suited for HR time series such as those provided by Sentinel-2.