Contribution of Earth Observation and Geospatial Information for Urban Planning of Historic Cities’ Centres: The Case Study of Nicosia, Cyprus

: The Sustainable Development Goals (SDGs) of the United Nations state that cities and human settlements need to be more inclusive, safe and resilient. In Europe cities have experienced dramatic physical, social and economic changes during the last decades while historic centres of European cities, among the most important assets of the European cultural heritage, are living paradoxes. They are deﬁned as “a collection of beauty, icon of well-being, model of sustainability, but abandoned”. This study investigates the changes in the urban landscape of Nicosia, a particular historical centre in the Mediterranean region (Cyprus). The city centre is characterised by exceptionally well-preserved Venetian fortiﬁcations. Due to political circumstances, the capital of Cyprus, Nicosia, is still divided and has been ruled by two different administrations for several decades. This study used optical multi-spectral satellite datasets processing, like the Landsat and the most recent Sentinel-2 products, to detect, identify and characterise signiﬁcant morphological transformations within the walled city and around it. This paper’s central thesis promotes a more systematic use of earth observation products and derivatives in decision-making processes that regard planning, use and management of urban resources in Europe, especially in support of urban planning strategies of historic cities.


Introduction
"Sustainable cities and communities", as Goal n. 11 of the United Nations' (UN) Sustainable Development Goals (SDGs), states that cities and human settlements need to be more inclusive, safe and resilient [1]. In particular, Goal 11.3 mentions that inclusive and sustainable urbanisation and the capacity for participatory, integrated and sustainable human settlement planning and management in all countries should be enhanced by the end of this decade (2030).
The current study focuses on the dimension of cities, as these are exposed to the various transformations and changes resulting from socio-economic and other dynamic activities. During the last few decades, European cities have experienced physical, social and economic changes due to relocations, conflicts and economic crises [2]. At the same time, historic centres of European cities are considered as important elements of the European cultural heritage, which are affected by tourist pressure [3]. The assumption is that cities' centres have been well-studied as they promote sustainability; however, the study of Pellegrini and Micelli [4] characterise them as living paradoxes, "a collection of beauty, icon of well-being, model of sustainability, but abandoned". surrounding areas are complex, a synoptic and spatial view offered from remote sensing is integral to measuring, modeling, and understanding these relationships." Therefore, including EO and other remotely sensed datasets within public SDIs can enable better analysis of urban management and promote more sustainabile interventions.
The current study aims to investigate and understand urban landscape changes and urban dynamics beyond traditional historical cities. Compared to many other European historical centers, the city centre under investigation, namely Nicosia (Cyprus), may be considered as an exception: the city has been divided and ruled by two different administrations for several decades (see Section 2). In addition, the interaction and relationship between Nicosia's historical city center and the rest of the modern sub-urban fabric is achieved though specific "gates", due to the radial alignment of the road network, since the city centre is characterised by an extremely well-preserved Venetian work of fortifications, an asset of high cultural value. This makes the case study unique, as we are aiming to understand how EO datasets and SDIs can support urban planning practices, given political constraints (the administrative division of the city) and the physical presence of an exceptional cultural asset (the Venetian fortifications).
In this study, we examine the changes in the urban fabric of Nicosia that occurred in the past. The study has been carried out using the Landsat and the most recent Sentinel-2 optical satellite datasets to detect, identify and characterise the significant morphological transformations within and around the walled city (Section 3). As discussed in Section 4, the objectives of the study, therefore, are the following: (i) Identify areas within and around the walls of Nicosia undergoing major changes occurring over the last 36 years; (ii) Identify the nature of these changes in specific areas; and (iii) Hence, observe possible city planning differences occurring within and around the Venetian walls (comparison of morphological changes). The previous objective would like to highlight the paper's central thesis for promoting a more systematic use of EO products and derivatives in decision-making processes regarding planning, use and management of urban resources. For this purpose, the National SDI of Cyprus was used as one of the primary references for evaluating the results of this paper (Section 5). The possibilities offered by the EO imagery analysis are hence discussed within a larger geo-spatial data framework as well as regarding their contribution to the urban planning strategies of historic cities.

Nicosia Historic City Centre
The main feature that characterises the Cypriot capital's shape and urban structure, namely Nicosia, is the Venetian walls that have maintained their geometric integrity for almost five centuries. Some historical details are provided below to describe the walls' setup and emphasise the significance of this structure.
The island of Cyprus was under Venetian rules in two distinct periods: during the first period from 1474 until 1489, the island was under the Venetian protectorate, while from 1489 until 1570, it was under Venetian direct sovereignty. Due to the rising power of the Ottoman Turks, Venice agreed to strengthen Nicosia's defence system; in 1567, the construction of the "royal fortress" began. For defence strategy purposes, the walled city's epicentre was shifted to the west while the circular perimeter, perhaps of Byzantine origin, was abandoned. Nicosia's Venetian walls were characterised by 11 equal bulwarks, "an impeccable topographic work", as stated by Marchesi in [19]. With that new defence system, the initial eight gates of the city were reduced to three: the gate of Saint Dominique (today known as Paphos gate), the gate of the Supervisor (today known as Kyrenia gate) and gate Gulia (today known as Famagosta gate). In the same period, the Pedieos river that was initially flowing through the centre of now fortified city, was diverted northwards by the Venetians for two reasons: to protect the residents from possible flooding and to fill the newly created moat encircling the new walls with water, therefore strengthening the defence system.
In the recent history of Cyprus, Nicosia was divided into the Greek and Turkish quarters in 1963. Since 1974, the whole island has been divided due to a Turkish military intervention. Today Nicosia remains a divided capital: the two parts of the city have been ruled by different administrations for more than four decades, and hence underwent dissimilar urban development changes. The central area of the city (as well as the whole island) is divided by a buffer zone (also known as the Green line), established in 1974 and under patrol of the "United Nations Peacekeeping Force in Cyprus", known as UNFI-CYP [20] (see Figure 1 in red). The Venetian walls surrounding the historic city, which have been exceptionally well-preserved for centuries, are under control of both parties on the north and on the south of the island [19]. example, in the case of road traffic, three gates are still in use as the main contemporary traffic hinges between the walled city and the surrounding areas. The Kyrenia gate is situated on the north side of the city. It is ruled by the Turkish Cypriot administration, while the Paphos and the Famagusta gates are situated respectively on west and east ends of the south part of the city that remained under the Government of the Republic of Cyprus. The three gates are illustrated in Figure 2.  The neighbourhood's division of the historic town also contemporarily divides the city's area beyond the fortifications. This fact influences some crucial city functions, such as the existing built environment, transportation, planning and other infrastructures. For example, in the case of road traffic, three gates are still in use as the main contemporary traffic hinges between the walled city and the surrounding areas. The Kyrenia gate is situated on the north side of the city. It is ruled by the Turkish Cypriot administration, while the Paphos and the Famagusta gates are situated respectively on west and east ends of the south part of the city that remained under the Government of the Republic of Cyprus. The three gates are illustrated in Figure 2.

National Spatial Data Infrastructures: Brief State of the Art in Cyprus
In Cyprus, the National SDI (NSDI) that fulfils the requirements of the European INSPIRE Directive [10] is coordinated by the Land and Surveyors Department (LSD). In 2011, the INSPIRE State of the art report evaluated the 'Cyprus Integrated Land Information System (CILIS)' as "an initiative considered to as an important step towards the development of a National SDI". The Cyprus LIS project was an umbrella program aiming to reduce duplication of land administration work among government agencies and to develop a Digital Cadastral Database (DCDB), a Survey Database (SDB) and a Topographical Database (TDB), suitable to support an integrated Land Information System. CILIS was defined as a tool for legal, administrative and economic decision-making and as an aid for planning and development [21]. The spatial data were documented and published in a catalogue, and these metadata were posted in January 2016 on the geoportal of the Republic of Cyprus [22]. As requested by the INSPIRE Directive, this geoportal is the main point for acquiring information regarding the implementation of INSPIRE in Cyprus [23].
It is also the main access point for e-services and geospatial data of the central government and semi-government organisations. The responsible Authority for the directive's implementation is the Ministry of Interior. At the same time, the INSPIRE Management Team in Cyprus is composed of the Department of Lands and Surveys, the Department of Environment and the Department of Information Technology Services. Town Planning and Housing Department are also included as key stakeholders and working groups.

National Spatial Data Infrastructures: Brief State of the Art in Cyprus
In Cyprus, the National SDI (NSDI) that fulfils the requirements of the European INSPIRE Directive [10] is coordinated by the Land and Surveyors Department (LSD). In 2011, the INSPIRE State of the art report evaluated the 'Cyprus Integrated Land Information System (CILIS)' as "an initiative considered to as an important step towards the development of a National SDI". The Cyprus LIS project was an umbrella program aiming to reduce duplication of land administration work among government agencies and to develop a Digital Cadastral Database (DCDB), a Survey Database (SDB) and a Topographical Database (TDB), suitable to support an integrated Land Information System. CILIS was defined as a tool for legal, administrative and economic decision-making and as an aid for planning and development [21]. The spatial data were documented and published in a catalogue, and these metadata were posted in January 2016 on the geoportal of the Republic of Cyprus [22]. As requested by the INSPIRE Directive, this geoportal is the main point for acquiring information regarding the implementation of IN-SPIRE in Cyprus [23]. It is also the main access point for e-services and geospatial data of the central government and semi-government organisations. The responsible Authority for the directive's implementation is the Ministry of Interior. At the same time, the IN-SPIRE Management Team in Cyprus is composed of the Department of Lands and Surveys, the Department of Environment and the Department of Information Technology Services. Town Planning and Housing Department are also included as key stakeholders and working groups.

Data Selection Over Time Period 1984-2020
The current study is focused on the changes observed in the urban fabric within and around the walled city of Nicosia over the last 36 years, using space-based observations. To minimise the impact due to seasonal vegetation changes, within and outside the historical centre, satellite images were selected over the year (the months May-June). The cloud coverage was set to be less than 5% for the scene of interest.
For these needs, the Landsat archive images of the NASA/USGS were explored. The Landsat science program offers the longest continuous global record of the Earth's surface, providing multi-spectral images at a medium resolution. In addition, the choice of Landsat imagery for urban changes referred to previous examples found in the literature, specifically on the consequences of urban sprawl [24][25][26][27].
The year selected referred to the less recent imagery available (1984), to the period of a few years before Cyprus has joined the European Union and to the year in which the study was conducted (2020). The imagery used was retrieved from the USGS data repository [28]. Table 1 shows some details regarding the images used in this paper.

Methodology
The workflow involved three processing phases. The first phase involves satellite data preparation. This includes the implementation of spectral indices relevant for the built environment and the data fusion using a literature-recommended band combination (Landsat 5 TM B7-B4-B2 combination). The second phase was focused on analysis and detection of areas subject to major changes using statistical methods such as the Standard Deviation (SD) and Principal Component Analysis (PCA). Finally, the third phase involved the assessment of the changes using high-resolution optical imagery from the Copernicus contributing missions and ground truth data from official publicly available INSPIRE data layers (3b). The overview of the workflow is illustrated in Figure 3. Regarding spectral indices in more detail, the normalised difference vegetation index (NDVI) takes into account the spectral difference between the red and the near-infrared (NIR) bands [29], while the NDBI exploits the short-wave infrared and NIR bands. Finally, the BUI is estimated as the difference between the NDVI and the NDB indices. The "built-up indices" were chosen as these can distinguish the urban areas  Regarding spectral indices in more detail, the normalised difference vegetation index (NDVI) takes into account the spectral difference between the red and the near-infrared (NIR) bands [29], while the NDBI exploits the short-wave infrared and NIR bands. Finally, the BUI is estimated as the difference between the NDVI and the NDBI indices. The "builtup indices" were chosen as these can distinguish the urban areas against the water bodies and vegetated area. This unsupervised classification brings the advantage of limiting the operator's subjective intervention in the mapping process [30]. The normalised difference built-up index (NDBI) [31] and the built-up index (BUI) [32] were used here are based on the spectral response of built-up areas due to their higher reflectance in the middle infrared wavelength range than in the near-infrared wavelength range. It should be stressed that all indices are normalised and provide results with a range of −1 to +1, thus minimising any seasonal differences between the different observation dates. Table 2 provides a summary and the equations of the indices used.

Spectral Signatures
Spectral signature diagrams were drawn based on the full spectrum of the images to identify potential spectral regions that could discriminate between the various targets (pixels) shown in the satellite images. Figure 4 presents indicative spectral signatures over six different targets as follow: vegetation (with green colour), urban areas (with red colour), soil areas (with brown colour), water bodies (with blue colour), asphalt targets (with black colour) and rivers (with cyan colour). The x-axis of Figure 4 shows the various wavelengths regions in terms of a micron (µm). As we see, water bodies tend to give relatively low pixel values, especially as we move on towards the near and the middle infrared part of the spectrum (after 0.80 µm). Vegetation shows a peak at the near-infrared part of the spectrum (at around 0.80 µm), and low values at the red part of the spectrum (0.60-0.70 µm), both indicative of spectral profile of healthy vegetated areas. A similar pattern to the one of vegetation was observed for the river target. This is because rivers are covered with dense vegetation. Urban areas show a decrease in the green part of the spectrum (at around 0.52-0.60 µm) and another decrease at the near-infrared part. An increase was also observed in the middle infrared part.
Findings from Figure 4 indicate that the various targets can be discriminated in the various parts of the spectrum, while the most appropriate to work with them are the green (0.52-0.60 µm), near-infrared (0.76-0.90 µm) and short-wave infrared (2.08-2.35 µm) wavelengths.
An overview of the three above mentioned indices in three distinct moments 1984, 2000 and 2020, is presented in Figure 5. More specifically, the NDVI index shows how non-built up areas within the walls (NDVI high values, i.e., white pixels) visible in 1984 (Figure 5a) increase over time and can be distinguished for the years 2000 and 2020 (Figure 5b,c). This is also evident in the "the buffer zone" area and around the Venetian walls of Nicosia: while in 2000, the southern part of the walls seemed to be more affected, in 2020, brighter pixels can be seen around the fortifications. As illustrated in Figure 5, in Nicosia, the land cover classes distinction considers mostly built areas and urban green areas.
infrared part of the spectrum (after 0.80 μm). Vegetation shows a peak at the near-infrared part of the spectrum (at around 0.80 μm), and low values at the red part of the spectrum (0.60-0.70 μm), both indicative of spectral profile of healthy vegetated areas. A similar pattern to the one of vegetation was observed for the river target. This is because rivers are covered with dense vegetation. Urban areas show a decrease in the green part of the spectrum (at around 0.52-0.60 μm) and another decrease at the near-infrared part. An increase was also observed in the middle infrared part. Findings from Figure 4 indicate that the various targets can be discriminated in the various parts of the spectrum, while the most appropriate to work with them are the green (0.52-0.60 μm), near-infrared (0.76-0.90 μm) and short-wave infrared (2.08-2.35 μm) wavelengths.

Comparison of Indices Used for Enhancement of Built-Up Areas
As already mentioned, for investigating land-use/land-cover changes and changes caused by urban sprawl, classification-based methods are commonly employed on remote sensing imagery. In [33], the authors compare environmental consideration and critical requirements between optical and non-optical sensors for urban change detection and monitor cities' expansion. Optical sensors used for change detection can provide information regarding detection and mapping of artificial changes, change rate, the spatial distribution of change types, analysis of urban change and growth, prediction of long-term trends of urban growth, etc. [33].
The second part of the analysis regarded the analysis using several indices, namely NDVI, NDBI and BUI. A comparison of existing built-up indices for urban area mapping, including NDBI and BUI, have been already examined on build-up areas of Mediterranean cities [18], while NDVI has been employed on multi-temporal Landsat TM imagery to derive overall trend changes in [29].

Comparison of Indices Used for Enhancement of Built-Up Areas
As already mentioned, for investigating land-use/land-cover changes and changes caused by urban sprawl, classification-based methods are commonly employed on remote sensing imagery. In [33], the authors compare environmental consideration and critical requirements between optical and non-optical sensors for urban change detection and monitor cities' expansion. Optical sensors used for change detection can provide information regarding detection and mapping of artificial changes, change rate, the spatial distribution of change types, analysis of urban change and growth, prediction of long-term trends of urban growth, etc. [33].
The second part of the analysis regarded the analysis using several indices, namely NDVI, NDBI and BUI. A comparison of existing built-up indices for urban area mapping, including NDBI and BUI, have been already examined on build-up areas of Mediterranean cities [18], while NDVI has been employed on multi-temporal Landsat TM imagery to derive overall trend changes in [29].
NDBI results show a similar interpretation to the NDVI. The dark pixels correspond to low values ranging from −0.1 to −0.4, indicating no water bodies in the area. The interpretation of the NDBI outcomes seems to confirm NDVI findings and aligned with the BUI index results: the buffer zone located in the middle of the walled city and surrounding Venetian walls have been subject to changes and an increase of (presumably urban) vegetation. However, in the southern part of the city, the vegetation tends to be visible around the walls; these are more clearly detected around the Kyrenia gate (see Section 5 for more details).
As far as the built-up areas, changes are visible over time in NDBI and BUI results. It could be argued that while in 1984 the urban district seemed to be more cohesive (Figure 5d,g), in the following years, the division due to the buffer zone is more evident. The built-up areas emerge as more highly concentrated in the centre of the historical city in the northern part and more distributed around main arteries and road connections in the city's southern part, even outside the city walls. In addition, some considerations can be made on Pedieos river as the riverbed used to cross the Nicosia city centre in a west-east direction (roughly corresponding to today's "buffer zone"): this area is today visible with high tones (white pixels) in NDVI figures (especially in the latest Sentinel

Standard Deviation Analysis over the Year-Pairs
A subsequent analysis involved a diachronic analysis, i.e., comparing the situ and changes over time in pairs of images (one per specific year) with a specific ban lection. In this case, the combination B7-B4-B2 (concerning the Landsat-5 sensor used. The literature suggests this pseudo-composite as the most promising one for u area enhancement and it was already tested on other urban areas [34]. The analysis h regarded artificially produced six-band multi-temporal images (see the workflo Figure 3 for details).  Figure 6a), changes could be observ It is important to note the difference in terms of spatial resolution highlighted in this step. Images acquired in 1984 and 2000 have a spatial resolution of 30 m (Landsat-5), while the 2020 Sentinel-2 provides a 10 m resolution (see Table 1). Such an increase in spatial resolution offers an added value for investigating both differences in land use/land cover and the morphological changes occurring within and around the city walls. Interesting features to observe in this sense include the road network and the distinction of single districts: in 1984 and 2000, a feature such as the Digeni Akrita street can be noticed in BUI (lighter pixels were surrounding the Venetian walls in the south-east, Figure 5g,h respectively), though it remains indistinct; in 2020, however, both the street and the existence of a dense urban fabric around it are visible in both NDBI and BUI (Figure 5f,i, respectively). Other important traffic arteries such as Larnakos Avenue, Pouliou-Kapota street (later Agiou Andreou) and Archiepiskopuo Makariou avenue can easily be identified, suggesting that a 10 m terrain pixels resolution is suitable for monitoring of the road network and hence possible respective changes in urban structure.

Standard Deviation Analysis over the Year-Pairs
A subsequent analysis involved a diachronic analysis, i.e., comparing the situation and changes over time in pairs of images (one per specific year) with a specific band selection. In this case, the combination B7-B4-B2 (concerning the Landsat-5 sensor) was used. The literature suggests this pseudo-composite as the most promising one for urban area enhancement and it was already tested on other urban areas [34]. The analysis hence regarded artificially produced six-band multi-temporal images (see the workflow in Figure 3 for details). Figure 6 shows the results of Standard Deviation (SD) analysis performed on these images for pairs 1984-2000 (Figure 6a), 2000-2020 ( Figure 6b) and 1984-2020 (Figure 6c). The results were expected to enhance the changes in the built-up areas over time. In the first period of 16 years considered (1984-2000, Figure 6a), changes could be observed in the city's historic centre and the south side. Also, some changes could also be noticed around the walls, such as the Famagusta gate (western part of the Walls) and in the south-west area around Digeni Akrita street. The second part of the time-lapse, from 2000 to 2020, seems to identify specific areas-one in the centre of the historic walled city and several cases around the Venetian walls. The central area could be referring to the change caused by the Ledra street passage between north and south Nicosia. In contrast, the changes around the Venetian walls seem to relate to specific interventions. the city's historic centre and the south side. Also, some changes could also be noticed around the walls, such as the Famagusta gate (western part of the Walls) and in the south-west area around Digeni Akrita street. The second part of the time-lapse, from 2000 to 2020, seems to identify specific areas-one in the centre of the historic walled city and several cases around the Venetian walls. The central area could be referring to the change caused by the Ledra street passage between north and south Nicosia. In contrast, the changes around the Venetian walls seem to relate to specific interventions.
Considering the period 1984 to 2020, the previously mentioned assumptions appeared to be confirmed: constant moderate changes occurred in the city within the walls with some "highlights" on both sides. The changes near the walls seem to have occurred on both sides on specific spots. In contrast, the urban fabric in a larger radius around the fortifications (circa 0.5 km) seemed to behave slightly differently. There are fewer changes that were observed by SD on the north, while on the south morphological changes seemed to follow some of the main city arteries.

Principal Component Analysis (PCA) over the Year-Pairs
In addition, the Principal Component Analysis (PCA) was performed on all three pairs: 1984-2000, 2000-2020 and 1984-2020. The paper examines the results of PCA in more detail on the pair 1984-2020 (Figure 7). Figure 7a shows the RGB natural colour composite of the Nicosia city centre, using the Sentinel-2 image. Similarly, Figure 7b shows the NIR-R-G pseudo colour composite of the same area. Pixels highlighted with red colour are used in the areas that are covered with healthy vegetation. From this view, the vegetated areas in the southern part of Nicosia are evident just outside of the Venetian walls. At the same time, the area of Pediaios river that crossed Nicosia is also clearly detectable. Figure 7c Figure 7d is the most important outcome of this analysis, as it highlights pixels with a red colour that have changed throughout the observation period. Considering the period 1984 to 2020, the previously mentioned assumptions appeared to be confirmed: constant moderate changes occurred in the city within the walls with some "highlights" on both sides. The changes near the walls seem to have occurred on both sides on specific spots. In contrast, the urban fabric in a larger radius around the fortifications (circa 0.5 km) seemed to behave slightly differently. There are fewer changes that were observed by SD on the north, while on the south morphological changes seemed to follow some of the main city arteries.

Principal Component Analysis (PCA) over the Year-Pairs
In addition, the Principal Component Analysis (PCA) was performed on all three pairs: 1984-2000, 2000-2020 and 1984-2020. The paper examines the results of PCA in more detail on the pair 1984-2020 (Figure 7). Figure 7a shows the RGB natural colour composite of the Nicosia city centre, using the Sentinel-2 image. Similarly, Figure 7b shows the NIR-R-G pseudo colour composite of the same area. Pixels highlighted with red colour are used in the areas that are covered with healthy vegetation. From this view, the vegetated areas in the southern part of Nicosia are evident just outside of the Venetian walls. At the same time, the area of Pediaios river that crossed Nicosia is also clearly detectable. Figure 7c shows the RGB pseudo colour composite using the first three principal components of the pair 1984-2020 as follows: a red colour covers areas with high values in the first principal component (PC1), while a green colour indicates pixels with high values in the second principal component (PC2); finally, a blue colour shows areas with high values at the third principal component (PC3). Individual results for the first three principal components (PC1-PC2 and PC3) are shown in Figures 7d, 7e and 7f, respectively. Figure 7d is the most important outcome of this analysis, as it highlights pixels with a red colour that have changed throughout the observation period. Having a closer look at this figure, we could extract areas of interest where this change was significant. As can be seen in Figure 7d, such hot-spot areas can be observed both within and outside the city walls of Nicosia.

Comparisson of the Results with Openily Available Geo-Spatial Information
To demonstrate how scientific findings in satellite remote sensing can be helpful, t study results include openly available repositories of other geo-spatial informatio namely: (i) digital globes repositories and (ii) Open Data found in geo-portals that pub administrations provide.

Evaluation of Results through HR Imagery and Digital Globes Repositories
Previous findings have revealed areas that changed dramatically throughout t

Comparisson of the Results with Openily Available Geo-Spatial Information
To demonstrate how scientific findings in satellite remote sensing can be helpful, the study results include openly available repositories of other geo-spatial information, namely: (i) digital globes repositories and (ii) Open Data found in geo-portals that public administrations provide.

Evaluation of Results through HR Imagery and Digital Globes Repositories
Previous findings have revealed areas that changed dramatically throughout the observation period. These changes can be due to land-use change as well as other seasonal changes or agricultural practices. To validate these findings, an archived high-resolution multi-spectral available from the Google Earth digital globe was used. Using the results obtained in PCA Component 1, six areas were selected for these purposes, as shown in Figure 8.   Table 3 shows various screenshots from these selected areas during 2011-2013 (middle column of Table 3) and recent satellite images of 2020 (right column of Table 3). For the area near point 1, we could not see any land-use changes; instead, there are some changes over the roof of the buildings. Small scale changes can be seen in the empty plot in the western part of this screenshot: today it is used as a parking place, while back in 2013 this was still an empty plot. Moving to the area around point 2, we can observe substantial changes, as the area was excavated (2013) and then was successively subject to a new construction intervention (2020). This is where the new Nicosia Municipality building is currently being constructed, preventing space observation of the excavated areas.        Table 4.  Table 4.  Table 4.

Evaluation of Results Using Open Geospatial Data and NSDI Services
Additional validation of results was done employing open geospatial data and services from the National Spatial Data Infrastructure (National INSPIRE geoportal). In particular, the data considered are shown in Table 4. The comparison of different layers is illustrated in Figure 9, specifically on six selected areas, the so-called "hot-spots" (Figure 8).
All images contain the layers provided by the NSDI with a WMS protocol. In particular, Natural Risk Zone layer, as proposed by the Water Development Department illustrates the Flood Plain boundary with 100 year return period (darker blue) and Flood Plain boundary with 500 year return period (lighter blue), as defined for Cyprus for the implementation of the EU Floods Directive [35]. The "buffer zone" area is represented with a grey-hatched layer.
Column A illustrates the hot-spots as seen on the Orthophoto WMS; column B shows the Land Use Land Cover layer (LULC); column C illustrates the Built-up index across all points for the year 2020; in column D, the close-up was made on the Standard Deviation (SD) image for the years 1984-2020. Images of the last two columns (C and D) were developed for this paper and have a resolution of 10 m (see Figures 5 and 6).
Looking into the six specific points, it can be noticed that points 1-4 are all near the areas identified as Natural Risk Zones, referring to flooding. In contrast, only points 5 and 6 are not in close proximity to these features. Points 2 and 4 are near the "buffer zone" area, while point 3 falls partially within these boundaries. In terms of land use, districts of points 1 and 5 refer to the LULC nomenclature "Industrial, commercial, military, public and private units"; points 2, 4 and 6 refer to a "Discontinuous dense urban fabric" (SL 50% to 80%); meanwhile, the district of point 3 refers to the category of "Continuous fabric" (with SL more than 80%). All points are well identifiable as white (light) pixels of column C (built-up areas) and red pixels of column D (significant changes in the urban fabric). Sustainability 2021, 13, x FOR PEER REVIEW 16 of 20 All images contain the layers provided by the NSDI with a WMS protocol. In particular, Natural Risk Zone layer, as proposed by the Water Development Department illustrates the Flood Plain boundary with 100 year return period (darker blue) and Flood Plain boundary with 500 year return period (lighter blue), as defined for Cyprus for the implementation of the EU Floods Directive [35]. The "buffer zone" area is represented with a grey-hatched layer.

Discussion
This paper supports the thesis that EO data and products can be used in urban planning and management practices, as well to contribute to socio-economic benefits in specific domains. In fact, public administrations have an increasing need for updated information on land use/land cover transformations in order to better tackle environmental policies.
For example, the results of this study could further investigate cluster analysis and spatial analysis to emphasise morphological changes, as have been conducted on other cities in Europe [36]. At the same time, a number of studies have been conducted regarding the use of space-based solutions for landscape pattern and vegetation fragmentation studies [37]. Additional applications on Nicosia could foresee the investigation of environmental "health" of the city and evaluation of quality of urban public green areas [38].
Observing the results of this paper, indices and statistics methods show specific changes over time, employing diachronic analysis of the same area across three different time thresholds. Diachronic analysis date back to the second half of the 18th century when the first cases of "diachronic" drawings of maps started to appear in urban cartography. Such drawings were illustrating the pre-existing territorial urban planning of a space and the projects intended for its reorganization, even if the latter was not yet implemented [39]. The experiments conducted here show how significant changes in the urban fabric ("hotspots") can be identified and compared to other sources of information to better understand their nature and purpose within future interventions.
The data evaluation confirms that the changes identified correspond to changes in the urban fabric (sometimes even in very densely urbanised districts) or a specific event with effects on the same fabric (e.g., excavation activities in the city center, point 4). Also, it could be noticed that many interventions involve very core districts of the Nicosia historic centre, on both the southern and northern side. In addition, two points (namely 1 and 6) were found in direct proximity to the Venetian walls, meaning that significant changes in the urban fabric have occurred respectively at a distance of 200 m and less than 100 m from this important cultural heritage asset and an urban landmark of the capital.
Specifically regarding the management of historic cities, these kinds of studies could be important when evaluating and planning the buffer-zone around Nicosia's Venetian walls. The term used here is not to be confused with the previously mentioned bufferzone of the overall Nicosia context but rather as intended by the UNESCO's definition for cultural heritage sites-an area to protect sites from negative influences [40]. In 2017, "The Venetian Works of Defence between 16th and 17th Centuries" was inscribed on the UNESCO World Heritage list. The existence of appropriate boundaries and buffer zones and protective measures were taken into account when evaluating integrity and authenticity of the selected monuments [41].

Conclusions
This study investigates the changes in the urban landscape of the historical centre of Nicosia, a Mediterranean capital characterised by an extremely well-preserved Venetian work of fortifications, an asset of high cultural value.
Some of the main conclusions of this paper are: - The study was conducted using free, full and open (FFO) earth observation data (namely Landsat and Sentinel-2 multi-spectral optical imagery). Free distribution satellite imagery can foster a better uptake of such information for targeted solutions and applications; - The study detects and characterises some major morphological transformations within the walled city and around it. Main changes in the urban fabric were identified, underlining that significant changes occur in already dense urban areas and in a close proximity to an extremely well preserved cultural heritage monument i.e., Nicosia Venetian walls; - The limitation of this research involves the evaluation of changes that was only conducted visually on previously collected datasets. Further investigations could employ the use of spectro-radiometers for ground-truth data collection, satellite image calibration and elaboration of supervised classification products; -In terms of drawbacks, in order to process and interpret satellite imagery, high technical capacity is needed both in terms of software licences and processor capability as well as skilled personnel; these could often be the main weaknesses when not employing space-based data more systematically within public administration bodies; - The results were further been compared to (i) to high-resolution imagery, and digital globes repositories and (ii) specific layers provided as WMS by the NSDI. These last steps highlight how satellite-based products can easily be consulted within GIS environments by professionals with no compulsory prior knowledge of remote sensing data elaboration; -In addition, the results of this paper can be provided as WMS within existing geospatial repositories in order to integrate already freely available geospatial information provided by public administrations.
As governments and local administrations are asked to promote and implement policies for more inclusive, safe and resilient cities in the future, integrated technological solutions should be put in practice to support decision-making processes that involve planning, use, and management of our urban resources. Inserting EO more frequently and more systematically within a larger geo-spatial data framework could provide a significant contribution to more sustainable urban planning strategies of the historic European cities, even in cities such as Nicosia, which are constrained by specific political situations but are rich in monuments with a high cultural importance and value.  Data Availability Statement: Publicly available datasets were analyzed in this study. This data can be found here: Satellite imagery: https://earthexplorer.usgs.gov/ (accessed on 30 April 2021). Open geospatial information on Cyprus: https://eservices.dls.moi.gov.cy/#/national/inspiregeoportalma pviewer (accessed on 30 April 2021).