**Non-Destructive Techniques for Building Evaluation in Urban Areas: The Case Study of the Redesigning Project of Eleftheria Square (Nicosia, Cyprus)**

#### **Marilena Cozzolino 1,\*, Vincenzo Gentile 1, Paolo Mauriello <sup>1</sup> and Agni Peditrou <sup>2</sup>**


Received: 1 June 2020; Accepted: 16 June 2020; Published: 23 June 2020

**Abstract:** This paper deals with the application of non-destructive geophysical techniques of investigation in the urban environment of the city of Nicosia (Cyprus). The main aim of the research was, in the frame of the Eleftheria Square redesign project, to image subsurface properties in order to reduce the impact of hazards on the old buildings (therefore preserving the cultural heritage of the place), and on the new infrastructure under construction. Since 2008, electrical resistivity tomography (ERT), ground penetrating radar (GPR) and induced electromagnetic method (EMI) were employed during the different phases of the project to provide an understanding of geological stratigraphy, the detection of buried objects (archaeological structures and underground utilities) and the solution of unexpected events (such as water infiltration in the course of works). The geophysical results proved the efficiency of the adopted methods, adding scientific value to the knowledge of the studied area. The new gathered information helped the public administration technicians to plan direct and targeted interventions and to modify the original design of the project according to the discovery of archaeological findings.

**Keywords:** applied geophysics; urban geophysics; Eleftheria Square; ERT; probability-based ERT inversion (PERTI) method; GPR; EMI

#### **1. Introduction**

Applied geophysics, through the employment of different methodologies, represents a useful tool to map subsurface features in a non-invasive way, and up to date, successful documented researches in rural settings dealing with different research objectives are available in the literature. Geological investigations for characterizing faults [1–3], landslides [4,5], paleo-morphologies [6,7], litho-stratigraphies [8,9], acquifers [10], sinkoles [11,12], and seepage detection [13] are generally carried out through electrical resistivity tomography (ERT) [1–3,5,7,9], ground penetrating radar (GPR) [8,10], seismic [2,4–6], magnetic [14–17] and gravity [18] methods. Archae-geophysical prospections are particularly common in applications for the detection of buried structures, tombs, and channels through the application of electromagnetic methods [19,20], magnetometry [21], ERTs [22,23], GPR [24,25] or the combination of them. The diagnostics and the monitoring of buildings in architectural and engineering surveys are often reached by implementing linear variable displacement transducer [26], sonic tomography [27], infrared thermography [27–29] and GPR [30–32] because of the high resolution of results.

In contrast, the geophysical exploration in urban environments for different purposes is currently in constant development. The consistent noise, the low signal-to-noise ratio, the reduced spaces to operate in, and the presence of paved surfaces, asphalt or complex layers disturbed by anthropological buried shallow artifacts and metallic structures represent a problematic difficulty in geophysical prospections that can significantly influence data quality. The rising interest of those aspects are evidenced by the inclusion of the matter in ad-hoc sessions of international workshops and conferences [33–36] and by the increment of special issues of international journals aiming to provide an enrichment of showcases for urban geophysics [37–39]. Nevertheless, some applications can be found in current literature, providing information on cavity individuation [40], geological mapping [41,42], archaeological investigations [43–45] and underground pipe detection [46,47].

In this work, ground-sensing prospections were used to characterize the subsurface in the areas involved in the redesign project of Eleftheria Square (Nicosia, Cyprus) (Figure 1).

The activities lasted more than ten years, from 2008 until today; the geophysical studies were an integral part of the project and here we report the various steps. The action was performed in the framework of a collaboration between University of Molise and the Nicosia Master Plan and was partly supported by the Ministry of Foreign Affairs. In particular, the research was in the frame of a scientific

archaeological mission that had not only the function of scientific study but also represented a valuable tool for intercultural dialogue and development policies in the host country. The archeo-geophysical surveys performed at Famagusta Gate, Castellionissa Hall, Phafos Gate, the monumental complex of the Bedestan and Plateia Di Marche are within this project [49,50]. The choice of Eleftheria Square as a study case is justified by the majesty and prestige of the construction project. In addition, as it provided for a distortion of the original setting with an enormous impact on its soils, the landscape and the urban system, there was the real need to acquire information before direct interventions.

In this case, a careful evaluation of the best method to implement to solve a particular problem was carried out and the ERT, GPR and induced electromagnetic (EMI) techniques were used in this contest. Sometimes investigations were conducted in an emergency situation by providing quick answers directly on the investigation site for problem solving.

According to the architectural work being planned, the surveys were organized during the executive phases of the intervention with the following aims:


From a scientific point of view, we wanted to show, through the chosen case study, the efficiency of this type of investigation in very disturbed environments such as an urban context where stratigraphies are often disturbed by recent or ancient buried anthropogenic elements. In particular, for ERT data processing, this was the opportunity to test the probability-based ERT inversion (PERTI) method [51], applied here for the first time to model near-surface urban geoelectrical data. The algorithm was developed contextually with needs that the case study presented, having, among its main features, a good filter capacity and the possibility of performing real-time inversion directly in the field, thus allowing for fast modifications of the survey plan to better focus on the expected targets [51].

In the following sections, a brief description of the case study and its cultural value is given, then the methodological approach to the research is explained and the results are reported and discussed.

**Figure 2.** Google Earth™ satellite images (2008, 2011, 2020) of Eleftheria Square and surrounding areas (**a**–**c**) and location of geophysical surveys (**d**).

#### **2. Case Study**

Nicosia has been the capital of Cyprus since the 10th century. During the Middle Ages, walls larger than the actual fortification were built to protect the city. The engineer Giulio Savorgnano designed the Venetian walls that subsist today during the Venetian period (1489–1570) for defense from Ottoman attacks. They have a circular shape containing eleven pentagonal bastions and incorporate specific innovative military architecture (Figure 1). The fortification had three gates known today as Famagusta Gate, Paphos Gate and Kerynia Gate (Figure 1a). At the end of the 19th century, the city began to expand towards the southwest direction outside of the walls and new openings were realized by constructing bridges over the moat connecting the old town with the surrounding areas. The main square of Eleftheria Square today is located at the site of one of the first openings that was created on the walls towards the end of the 19th century in order to facilitate the direct connection of the old city with the newly developed areas around the walls (Figure 3). In 1882, a bridge in wood was built but later it was substituted with an in-filled solid bridge with stone walls. In 1930 during the English regime, the passage was widened while demolishing part of the bastion of the Venetian walls. The square has been a timeless gathering place for citizens linked to important historical moments of the country and the major celebrations of the capital city.

In 1974, the Turks invaded the island, occupying the 37% of it. The so-called "Green Line" (Figure 1a) separates the island in two parts; crossing the city of Nicosia actually represents today the unique divided capital in Europe. Due to this division, the historic city center suffered from depopulation, significant physical decay, loss of commercial activities and employment, and high concentration of social problems [52,53]. Since April 2003, it was possible to cross the Green Line and move between the two regions. In 2008, the first passage was unlocked inside the historical center of the capital through Ledra Street (Figure 1a).

In 1979, a significant agreement was reached between the Turkish and Cypriot community in order to prepare a bi-communal master plan with the aim to improve current and future living conditions for all residents of Nicosia. The Nicosia Master Plan was put under the auspices and the financial backing of the United Nations Development Programme (UNDP), incorporating the following objectives [53,54]:


**Figure 3.** Photographs of Eleftheria Square (Metaxa Square) in 1945 (Photo captured by J. Arthur Dixon) (**a**,**b**), 1957 (Photo captured by Felix Giaxis) (**c**), 1950 (**d**) and 1969 (**e**).

Thus, in the last decades the town has witnessed the growth of projects focusing on restoring historic buildings, preserving cultural heritage, and calling for new constructions in the periphery and along the historic fortification walls. An example is testified by the renovation of Elefteria Square, an important link between the old town and the modern city outside the walls. In 2005, an architectural competition assigned the project to Zaha Hadid Architects, an architectural office of international reputation, that proposed a modern solution which was blending the Venetian walls, the ditch and the limits between the two parts of the city [55] (Figure 4). The project, directed by Christos Passas, transformed the moat into a green park in combination with the restored Venetian walls, offering recreational facilities, spaces to be used for sport and cultural activities, and art exhibitions, while the upper part became the central core with a modern pedestrianized plaza. The eastern Solomon square was also restored with a new bus station (today completed) and a car park, which is currently under implementation under Omirou Avenue.

**Figure 4.** Eleftheria Square project: photographs of the work in progress captured in December 2019.

#### **3. Methodological Approach, Data Acquisition and Processing**

As a premise to the methodological approach used, it should be noted that the project of Eleftheria Square is a long-term plan that began in 2005 with its design and, today, its implementation is still ongoing. Furthermore, during the years the geophysical support was needed for the various purposes previously introduced, which cannot be separated from a case-by-case evaluation of the methodology to be adopted. Thus, the following workflow was considered:


(point-to-point) mode using frequencies in the range 5–16 kHz with vertically oriented dipoles. The investigation area was covered by profiles spaced 0.5 m inserted in a regular grid.

	- For ERT, the measured apparent resistivity datasets were processed in order to remove dragging effects that are typical of DD array and model the survey targets converting the values in real electrical resistivity values displayed as a function of depth below surface. To this end, a probability tomography approach has been applied for imaging the sources of anomalies into the analyzed grounds. The theory was first stated for the self-potential method [58], and then adapted to the resistivity method [59–61]. The primary approach, even if capable of distinguishing resistivity highs and lows in the field datasets considering a reference background resistivity, precluded the estimate of the intrinsic resistivities of the source bodies. The method was successfully used to delineate the geometry of the self-potential sources in the central volcanic area of Vesuvius (Naples, Italy) [62], to imaging through the resistivity method buried archaeological structures in the archeological site of Pompei (Naples, Italy) [63], at the Castle of Zena (Carpeneto Piacentino, Italy) [64], the prehistoric sites of Checua (Cundinamarca, Colombia) [65] and Grotta Reali (Rocchetta a Volturno, Italy) [66], the Archaeological Park of Aeclanum [67] and for fault detection in a site on the Matese Mountain (Italy) [3]. Successively, a data-adaptive probability-based ERT inversion (PERTI) method [51] was directly derived from the principles of the probability tomography in order to estimate the true resistivities. From a probabilistic point of view, the algorithm, being a non-linear approach, identifies, inside the set of possible solutions the most probable one, compatibly with the dataset acquisition scheme. The main features of the PERTI method, as reported in [51] are: (i) unnecessity of a priori information; (ii) full, unconstrained adaptability to any kind of dataset, including the case of non-flat topography; (iii) drastic reduction of computing time of even two orders of magnitude, with respect to the previous methods in complex 3D cases using the same computer; (iv) real-time inversion directly in the field, thus allowing for fast modifications of the survey plan to better focus the expected targets; (v) full independence from data acquisition techniques and spatial regularity, (vi) possibility to be used as an optimum starting model in standard iterative inversion processes in order to speed up convergence. As the PERTI method does not require a priori information and iterative processes, the computation of the route mean square (RMS) error between measured and modelled apparent resistivity values is useless. In fact, the resulting RMS error, whatever it is, can be lowered within the PERTI scheme. Many applications of the PERTI approach in near surface prospections are available in literature for solving archaeological research questions [49,68–71], for defining faults in Crete [72] and for imaging of the near-surface structure of the Solfatara crater, Campi Flegrei (Naples, Italy) [73]. In [49,51,68,72] the PERTI routine was tested using well-known commercial software of inversion of geoelectrical data and the comparison put in evidence for the coherence between the obtained results and the better filter capacity and great versatility of the PERTI algorithm. Here, for the first time, the method is applied to process field datasets acquired in urban environments, as its main features fit perfectly with the needs of the case study presented.
	- For GPR, as the interpretation of each section can lead to underestimation or overestimation of the reflected signals and makes it not easy to identify the effect of lateral bodies present in the subsoil, all sections were processed together using standard methodological approaches in order to obtain 2D horizontal maps at a different range of depth (GPR-SLICE 7.0 software) [74]. Data were converted by subtracting out the dc-drift (wobble) in the data, and at the same time adding a gain with time of 20. A time-zero correction was determined to designate the starting point of the wave and the center frequency of the antenna was matched. Then the bandpass filter and the background removal were respectively applied to reduce noise from oscillating

components that had a regular frequency cycle in the frequency domain and to remove striation noises that occurred at the same time. Processed radargrams were subsequently corrected with an automatic gain function applied to each trace based on the difference between the mean amplitude of the signal in the time window and the maximum amplitude of the trace. Thus, horizontal sections (time slices) were processed considering the whole dataset. Data were gridded using the inverse distance algorithm, which includes a search of all data within a fixed radius of 0.75 m of the desired point to be interpolated on the grid and a smoothing factor of two. Grid cell size was set to 0.01 m to produce high resolution images.


#### **4. Results**

#### *4.1. Survey 1*

In 2008, five ERT profiles, 36 m long, were conducted at Eleftheria Square with the aim to cross out possible remains of ancient walls and traces of portion of D'Avila Bastion (Figure 5). It clearly appears to have an asymmetrical shape, with the uncertainty that they were destroyed during the construction of the openings in the circuit at the beginning of the 19th century. The distance between the close profiles 2–4 was 1 m. The maximum reached depth was about 12 m at the center of the sections.

At 1 m in depth at the northeastern side of ERT 1 and ERT 5, two resistivity highs (about 300 Ohm m), vertically elongated, were detected (Figure 6a,d). They are precisely located were the portion of Venetian walls should be as pointed out with yellow dots in Figure 7a.

ERT profiles 2–5 highlight coupled medium resistive spots (about 100 Ohm) spaced 4 m; the location is reported with blue circles in Figure 6b–d. They were interpreted as the walls of the bastion. The conductive zones between them can be associated with the soil filling the fortress. In ERT 1 and ERT 5, the section of the ditch is represented by low resistivity values from the surface until a maximum depth of 8 m. It rests on a resistive layer, the probable deep bedrock.

**Figure 5.** Google Earth™ satellite images (2008) with the indication of gaps in Venetian wall (in blue) research subjects of Surveys 1 and 2.

**Figure 6.** Eleftheria Square: ERT profiles 1–5 (**a**–**d**) and the modelled horizontal resistivity map relative to 1.2 m in depth (**e**).

Figure 6e shows a horizontal resistivity map in which the resistivity highs are located on a technical plan of the area. In order to better highlight the medium resistive spots, the limits of color scale have been slightly modified.

An archaeological excavation performed in 2009 proved the existence of the Venetian walls in the point indicated at the end of ERT 5 and the bastion was entirely brought to light in the exact points indicated by geophysical prospections (Figures 6 and 7a). The ditch was also partially excavated. Subsequently, the original project of the new square needed to be modified, maintaining a distance of 3 m between the new construction and the northwestern side of the bastion as regulated by the International Charter for the Conservation of Monuments and Sites. The new construction is now completed, and it incorporates the ancient bastion (Figure 7b).

**Figure 7.** Eleftheria Square: Location of ERT profiles 1–5 and geophysical anomalies on a Google Earth™ satellite images (2011) after archaeological excavations at Eleftheria Square (**a**) and a picture of the excavated ancient bastion inserted in the new square (**b**).

#### *4.2. Survey 2*

Survey 2 was carried in 2008 and concerned Solomos Square, where the city's bus station is, at the intersection of Rigenis Street and Omirou Avenue (Figures 2 and 8a). Here, the presence of buried walls and archaeological remains were the target of the investigation. Five ERT profiles (6–10) with different lengths where acquired. ERT profiles 7–9 (Figure 8d–f), located at the center of the square, clearly put forth evidence of the filling of the ditch and the deep resistive substratum. With the exception of rare near-surface high-resistivity nuclei probably due to modern alteration, the unique dishomogeneities that can be attributable to ancient walls are the high-resistive features located at the right side of ERT profiles 6 and 10 (yellow dots in Figure 8a,c,f). In 2010, the square was remodeled, and direct archaeological excavations were not performed to verify geophysical results.

**Figure 8.** Solomon Square: location of ERT profiles 6–10 (**a**), data acquisition (**b**) and geophysical results (**c**–**f**).

#### *4.3. Surveys 3–4*

In 2009, twelve ERT profiles were carried out along Omirou Avenue with the aim of imaging the stratigraphy of the soil and to verify the presence of archaeological structures into the subsoil (ERT profiles 11–22 in Figure 9). As the planning of the survey was influenced by the necessity to ensure the traffic of cars, two main fragmented lines of investigation were conducted (blue and red lines in Figure 9). Profiles had variable lengths and the depth reached for all was 10 m. Figure 10 shows two tridimensional perspective views of the imaged ERT profiles. In general, a uniform stratigraphy is detected and characterized by materials with low resistivity values. In some cases, an increase of resistivity can be attested at about 8 m in depth. High resistivity features that can be associated to huge empty cavities were not detected, as well as it was not possible to distinguish small cavities filled up with loose conductive sediments.

**Figure 9.** Omirou Avenue: location of ERT profiles 11–22. Blue and red lines indicate the two main fragmented lines of investigation.

**Figure 10.** Omirou Avenue: 3D perspective views of ERT profiles 11–16 (**a**) and 17-22 (**b**).

Resistivity features with well-defined shapes similar to the ones that were detected on the ERT profiles of surveys 1 and 2 are not attested here, such as to hypothesize the presence of archaeological remains under the road plane. However, isolated shallow high resistivity values were highlighted in each section (marked with black arrows in Figure 10) that can suggest the existence of nuclei with a non-vanishing probability to find ancient remains. Noteworthy are the two high resistive nuclei that repeat on the ERT profiles 13, 14, 19 and 20, which seem to be on the same route.

In 2017, before the removal of the asphalt layer from the road, a GPR survey was performed. Figure 11a displays two examples of radargrams in which punctual targets, probable manholes, and horizontal discontinuities are marked with coloured arrows. The same marks are reported on the horizontal slice relative to the time window 10–15 ns, the one that better enhances the presence of targets into the subsoil (Figure 11b). The map shows the presence of different linear high amplitude values of the electromagnetic signal that probably refer to pipes buried into the shallow ground. Particularly evident is the presence of a high amplitude feature flanking the south side of the road attributable to the underground sewer pipes (with dotted lines in Figure 11b). Furthermore, the map clearly highlights the presence of manholes indicated with crosses in Figure 11b.

**Figure 11.** Omirou Avenue: two examples of radargrams (**a**) and GPR slice relative to the time window 10–15 ns. (**b**) Dotted lines and crosses indicate respectively underground pipes and manholes; the colored arrows mark the same anomalies on radargrams and slice.

#### *4.4. Survey 5*

In 2019, an unexpected problem occurred during the works: inside the underground substation where electrical cables arranged in six plastic pipes converge, water began to flow in a very abundant cascade directly from the pipes. To try to find out the cause, an excavation trench was carried out south of the cabin (Figure 12a), partially bringing to light the cables that were embedded inside a layer of reinforced concrete. The action was not conclusive as the leak could not be detected. The most reliable hypothesis of the origin of the water was the breakdown of the sewage pipes detected by the GPR in 2017 that pass to the center of the road junction, but no evidence was identified by direct verification. To this point, GPR and EMI investigations were performed in the very few spaces available (Figure 12b).

The resistivity maps obtained from the EMI survey allowed for locating the most probable wet areas. Figure 13 shows all the horizontal maps relating to the frequencies from 16 kHz to 5 kHz. At the edge of the sidewalk there is a conductive band that at lower frequencies (which give us information at greater depths) gradually flags and highlights the retaining wall of the underground rooms in the square. The slice that showed the highest conductivity is that related to the 11 kHz frequency (Figure 14b). The maps relative to lower frequencies put in evidence of a decreasing of the conductive values and highlight the resistive piles of reinforcement.

**Figure 12.** Corner between Evagorou Avenue and Omirou Avenue: the excavated trench (**a**), GPR and EMI data acquisition (**b**).

As expected, the GPR survey was not a solution in wet areas because the water completely absorbed the electromagnetic signal (Figure 14a). However, the buried electric cabin was precisely reported on the map and two nuclei with high amplitude of the electromagnetic signals, possible manholes, were detected.

A targeted excavation verification at that point confirmed the presence of conspicuous water into the soil. Although the origin of the water has not been identified, the investigation has averted the southern origin of the leak as hypothesized before the investigation. The most probable theory is that the problem should originate to the west of the trench where the pipeline of the city aqueduct is located. The issue is still being analyzed through inexpensive and non-invasive solutions.

**Figure 13.** EMI resistivity maps relative to the frequencies in the range 16–5 kHz.

**Figure 14.** GPR slice relative to the time window 14–18 ns (**a**), EMI resistivity map relative to the frequency of 11 kHz (**b**).

#### **5. Conclusions**

This paper concerned the contribution of non-destructive investigations within the redesigning project of Eleftheria Square and its cultural value protection according to its important role in the life of the city of Nicosia.

As a premise which forms an integral part of the discussion of the reported results, it is important to emphasize that the unique limitation of the methodological approach used in this work is in the general principle of geophysical diagnostics as it is an indirect methodology aimed at the spatial definition of anomalies present in the subsoil. Only a direct verification in the field can give certainty of the nature of the buried objects and in the case of ancient finds can define their age and the real archaeological value of the anomalies found. On the contrary, the strength of preventive investigations is the non-invasiveness of the survey, the speed of data acquisition and the possibility of being able guide direct investigations in a targeted and precise manner.

The Elephteria Square redesign project is an important case study for its majesty and prestige and for the historical and cultural context within which it is inserted. The activities lasted more than ten years, from 2008 until today, and the geophysical researches were an integral part of the project. Geophysical research combined with this type of initiative is important as it can be easily replicated in one of the many historical centers scattered around Cyprus or in other parts of the world interested in changes in new urban planning that can have an impact on modern and ancient structures.

The obtained results provided a precise answer to the research questions. In particular, the evaluation of the potential archaeological risk at Eleftheria Square and Solomon Square in 2008, before the beginning of the construction of the new structure, was successfully reached through the application of ERTs: the highs of resistivity interpreted as buried archaeological constructions were proved to be, after excavation, part of the Venetian wall and portions of D'Avila Bastion. Such a result allowed for preserving pieces of the local culture, avoiding its destruction. Furthermore, once discovered, they were inserted within the new complex of the square, adopting a modification to the project in order to highlight the archaeological findings.

The ERT surveys along Amirou Avenue in 2009 have permitted to delineate the geological stratigraphy of the soils, clarifying the absence of huge empty cavities and suggesting the existence of near surface resistive nuclei that probably can be attributable to ancient remains.

Regarding the methodological approach used to process electrical resistivity datasets, the data-adaptive probability-based ERT inversion (PERTI) method was proven to be appropriate for urban measurements because of its ability to remove dragging effects due to the dipole-dipole configuration used and to filter, sometimes, very noisy raw data. In general, within the first application of the probability tomography in urban settings, a good performance in drawing sections is proved.

The GPR surveys, using standard technique of data acquisition and processing, along Amirou Avenue in 2017 allowed the precise detection of underground utilities facilitating the technicians during works of excavation and avoiding unexpected damages to city supplies.

Finally, GPR and EMI investigations were carried out in order to understand the reason for an abundant flow of water that occurred inside the new underground electrical substation positioned under the new Eleftheria Square during excavations at the intersection of Evagorou Avenue and Omirou Avenue in 2009. As expected, the GPR survey was not able to overcome the limit of the presence of water but it was useful to report the precise location of the buried electric cabin on the map. The resistivity maps obtained from the EMI survey allowed the precise locating of the most probable wet areas. A targeted excavation verification at that point confirmed the presence of conspicuous water into the soil, averting the southern origin of the leak as hypothesized before the investigation.

In general, the geophysical results proved the efficiency of the adopted methods, adding scientific value to the status of knowledge of the subject area of study. In some cases, the verification of the supposed buried features was proved by excavation and a perfect correspondence between the geophysical previsions and the structures found in the subsoil was attested. In addition, the presented research supported the actions of the technicians involved in each phase of the project and often, in emergency, the results were provided in real time to resolve the problems encountered. A result of all the various urban activities is, in fact, the collaboration between geophysicists and urban decision makers, which should be encouraged in cases similar to that presented for an optimal city management.

**Author Contributions:** Conceptualization, M.C. and P.M.; methodology, M.C. and P.M.; software, M.C. and P.M.; validation, M.C. and P.M.; formal analysis, M.C. and P.M.; investigation, M.C., V.G. and P.M.; data curation, M.C., V.G. and P.M.; writing—original draft preparation, M.C., P.M. and A.P.; writing—review and editing, M.C. and P.M.; supervision, P.M. and A.P.; project administration, P.M. and A.P.; funding acquisition, P.M. and A.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** The Ministry of Foreign Affairs of Italy, in the frame of the Archaeological Missions, and Nicosia Master Plan funded part of this research. We thank the Embassy of Italy and the Ambassador Andrea Cavallari for the continuous support for the activities and His Excellency Gherardo La Francesca for being the promoters of this Italian-Cypriot collaboration. Special thanks to the technicians of Nicosia Municipality for administrative and technical support during the research activities.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Applied Sciences* Editorial Office E-mail: applsci@mdpi.com www.mdpi.com/journal/applsci

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel: +41 61 683 77 34

www.mdpi.com

ISBN 978-3-0365-5132-6