*Article* **Level of Contamination Assessment of Potentially Toxic Elements in the Urban Soils of Volos City (Central Greece)**

**Evangelia E. Golia 1,\*, Sotiria G. Papadimou 1,2, Christos Cavalaris <sup>3</sup> and Nikolaos G. Tsiropoulos <sup>2</sup>**


**Abstract:** A three-year study, designed to record the level of potentially toxic elements within the urban complex in the city of Volos, Greece, was carried out between 2018 and 2020. For the needs of the aforementioned study, 62 surface (0–15 cm) soil samples were collected each year (i.e., 186 samples in total) from an urban area of 3.65 km<sup>2</sup> , and the average value of pseudo-total metal concentration was measured. Soil pollution indices, such as the contamination factor (CF) and the geo-accumulation index (Igeo), were estimated regarding each of the metals of interest. The respective thematic maps were constructed, and the spatial variability of the contamination degree was displayed. Higher values of the CF and Igeo were obtained near the heavy traffic roads and beside the railway station, the bus stations, and the commercial port. The maps based on the pollution indices, along with the database that was constructed using the appropriate mathematical tools of geostatistical analysis, may be a useful tool for monitoring, prediction, and continuous verification of contamination in the urban soils of Volos city.

**Keywords:** contamination indices; geo-accumulation factor; heavy metals; GIS

#### **1. Introduction**

Important functions vital for agricultural, environmental, land preservation, perspective architecture, and urban activities usually take place in the soil [1,2]. Urban soils are inside or nearby cities, i.e., in urban and suburban areas that are significantly and persistently disrupted by human activity [3]. They include: (a) soils consisting of a mixture of materials of different origins, with inorganic or organic compositions, and from agricultural or forest areas, which are largely transformed by human intervention [4]; (b) soils in parks and green gardens that offer different compositions, uses, and management of agricultural land [5,6]; and (c) soils or mixtures resulting from various construction or metallurgical activities within cities [7].

Urban soils are created when human intervention takes place for many years, or after an immediate and abrupt concentration of human activities in a specific area and in a specific environmental background [8–10]. Compared to natural and agricultural soils, urban soil inherits certain characteristics of the parent rock, but it is also influenced by the local microenvironment as well as the process of its formation. Changes in the physical and chemical properties and in the biological cycles of nutrients are often observed.

Urban areas are often densely populated, so there is an over-consumption of natural resources and high production of pollutants. In recent decades, there has been an increase in urbanization and industrialization, as short-term and intensive human activities lead to large numbers of organic pollutants (such as polyhalogenated compounds), along with inorganic pollutants (such as potentially toxic metals) [11,12].

**Citation:** Golia, E.E.; Papadimou, S.G.; Cavalaris, C.; Tsiropoulos, N.G. Level of Contamination Assessment of Potentially Toxic Elements in the Urban Soils of Volos City (Central Greece). *Sustainability* **2021**, *13*, 2029. https://doi.org/10.3390/su13042029

Academic Editor: George D. Bathrellos Received: 18 January 2021 Accepted: 10 February 2021 Published: 13 February 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 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 (https:// creativecommons.org/licenses/by/ 4.0/).

It is well known that Potentially Toxic Elements (PTEs) accumulation in the soil leads to soil degradation, which can be a serious problem concerning human health, as metals can easily accumulate in the human body through food or ingestion of water as well as through respiration pathways [13]. Ecological risk assessments and their potential impact on human health are of great importance because they can point to unbearable dangers related to man and the environment [14]. In addition, they also require the identification of the specific driving factors that contribute to these risks [15]. The accumulation of metals in the upper soil can threaten health through consumption, digestion, and skin contact [3], and soil toxicity by heavy metals is highly correlated with metal levels and nature. In many cases, damage to human health is inevitable: headaches, insomnia, insanity, joint pain, and cancer. They can have a stronger effect on human health, plants, and soil organisms as well as in water filtration [16]. Problems due to urban soil systems are distinguished by other strongly affected soils, such as those found in quarries, mines, and airports far from cities [17]. There are different routes of exposure, such as direct ingestion through the food chain, inhalation through the mouth and nose, and direct skin contact. Therefore, it is important to assess the contamination and risk to human health of soil metals to improve the environment and protect human health [18,19].

Identifying potential sources of metals in soils is important for controlling priority pollutants. Recently, various mathematical models and methods of statistical analysis have been applied, such as the various indices usually called environmental factors (EFs) such as, the contamination factor (CF), the geo-accumulation index (Igeo), the availability factor (AF), principal component analysis (PCA), and factor analysis (FA) [20,21]. These tools are highly used to distinguish the natural and human origins of metals in soil systems [17] On the other hand, GIS in environmental pollution is a valuable tool for mapping out the contaminants in soil. Spatial display and pollutant distribution allow for a more efficient perspective to monitor soil contamination [22,23]. Spatial data management has a large range of applications along with great advantages [24,25] Methodologies incorporating GIS can be applied in soil contamination projects to facilitate the interpretation, estimation, and evaluation of information between stakeholders, which will improve stakeholder communication in the decision making process [26].The potential burden on the urban environment due to the presence of heavy metals in urban soils, along with the particles suspended in the air and deposited on the soil surface, has gained intense research interest in recent years [27,28]. However, limited studies are focused on risk assessment in urban and suburban areas, on densely populated areas, or on the industrial cities of Greece [29–31].

The objectives of the present study were (1) to evaluate heavy metal contamination levels in the urban soils of Volos, (2) to identify which PTEs are of most concern using soil contamination indices, and (3) to highlight the areas (sites) that are most at risk of pollution by constructing proper thematic maps to illustrate the spatial distribution in the study area. The novelty of the present research lies in the fact that, for the first time, a database has been created that enables the continuous monitoring of contamination assessment after recording the changes of the contamination index values in the area studied.

#### **2. Materials and Methods**

#### *2.1. Site Description*

Volos is a seaside port of Thessaly located in the middle of Greece. It was built inside the Pagasiticos Gulf and spreads along the plain beneath Pelion. It has a population of about 150,000 inhabitants. The area has a typical Mediterranean climate with wet and cold winters, dry and hot summers, a 14.4 ◦C average yearly temperature, and a mean yearly precipitation of 450 mm. At the west side, about 12–14 km from the city center, resides an industrial area that also includes a steel plant. Another large cement industry, producing seven types of cement, clinker, solid fuels, and aggregates, is located a distance of 4 km to the east part of the city (Figure 1). The present study is focused on the downtown city center, covering an area of 3.65 km<sup>2</sup> .

.

seven types of cement, clinker, solid fuels, and aggregates, is located a distance of 4 km to the east part of the city (Figure 1). The present study is focused on the downtown city

**Figure 1.** The city of Volos in Thessaly, Central Greece and the study area. **Figure 1.** The city of Volos in Thessaly, Central Greece and the study area.

#### *2.2. Soil Sampling 2.2. Soil Sampling*

center, covering an area of 3.65 km<sup>2</sup>

Sixty-two soil samples per year, consisting of three sub-soils each [32], were collected from diverse green spaces, including parks, playgrounds, squares close to the harbor, flower beds located close to main streets, and a neighborhood adjacent to a cement factory. The investigation began in June 2018. In the following two years (2019–2020), soil samples were collected from the same sampling points (i.e., 186 soil samples in total) of Volos. Sixty-two soil samples per year, consisting of three sub-soils each [32], were collected from diverse green spaces, including parks, playgrounds, squares close to the harbor, flower beds located close to main streets, and a neighborhood adjacent to a cement factory. The investigation began in June 2018. In the following two years (2019–2020), soil samples were collected from the same sampling points (i.e., 186 soil samples in total) of Volos.

#### *2.3. Chemical Analysis of Soil Samples 2.3. Chemical Analysis of Soil Samples*

as the maximum allowable concentrations.

The physicochemical analyses of the soil samples were conducted using the methods described by Page et al. [33]. The analyses were preceded by air-drying the soil samples for three days and sieving them using a 2 mm sieve. Clay, sand, and silt percentage, as well as cation exchange capacity (CEC), pH values, electrical conductivity, organic matter content (Table 1), and pseudo-total content (Table 2), were measured using the aqua regia method [34,35]. Metals were determined using an atomic absorption spectrophotometer (flame and/or graphite furnace) with the following detection limits: 0.1, 0.09, 0.09, 0.1, 0.08, 0.2, 0.1, 0.15 (mg L–<sup>1</sup> ) for Cu, Zn, Cd, Mn, Pb, Ni, Cr, and Co, respectively. A series of 11 standard solutions (Perkin Elmer® ) for each metal were used after conducting the proper solutions with 5% c. HNO3 [35]. For the verification of the accuracy of the analyses, a certified reference material (CRM) (No 141R, calcareous loam soil) from the Community Bureau of Reference (BCR) was analyzed with the soil samples. The results of the metal determination in the CRM are presented in the Supplementary Materials section. The European Council Directive 86/278/EEC [36] regulated the permissible limits for Cd, Cu, Co, Ni, Pb, and Zn. The Cr and Mn maximum values were reported in Kabata–Pendias [37] The physicochemical analyses of the soil samples were conducted using the methods described by Page et al. [33]. The analyses were preceded by air-drying the soil samples for three days and sieving them using a 2 mm sieve. Clay, sand, and silt percentage, as well as cation exchange capacity (CEC), pH values, electrical conductivity, organic matter content (Table 1), and pseudo-total content (Table 2), were measured using the aqua regia method [34,35]. Metals were determined using an atomic absorption spectrophotometer (flame and/or graphite furnace) with the following detection limits: 0.1, 0.09, 0.09, 0.1, 0.08, 0.2, 0.1, 0.15 (mg L−<sup>1</sup> ) for Cu, Zn, Cd, Mn, Pb, Ni, Cr, and Co, respectively. A series of 11 standard solutions (Perkin Elmer®) for each metal were used after conducting the proper solutions with 5% c. HNO3 [35]. For the verification of the accuracy of the analyses, a certified reference material (CRM) (No 141R, calcareous loam soil) from the Community Bureau of Reference (BCR) was analyzed with the soil samples. The results of the metal determination in the CRM are presented in the Supplementary Materials Section. The European Council Directive 86/278/EEC [36] regulated the permissible limits for Cd, Cu, Co, Ni, Pb, and Zn. The Cr and Mn maximum values were reported in Kabata–Pendias [37] as the maximum allowable concentrations.



Skewness Coefficient 0.710 0.747 –0.128 1.0096 –0.932 0.231


**Table 2.** Pseudo-total concentrations of potentially toxic elements (mean values of the three years of the study, n = 186).

<sup>a</sup> 10th percentile, <sup>b</sup> 50th percentile, <sup>c</sup> 90th percentile; <sup>d</sup> For Cr: maximum allowable concentrations ([37], p. 24). For Cd, Cu, Ni, Pb, and Zn: 86/278/EEC Directive [36]; <sup>e</sup> Average of all reported background reference values in Kabata–Pendias [37].

*2.4. Contamination Risk Assessment Indices of PTEs*

The following PTE contamination indices were estimated. [17,38]. The contamination factor (CF) is as follows:

$$\mathbf{C}F = \mathbf{C}\_{\mathbf{AR}} / \mathbf{C}\_{\mathbf{AR} \, nef} \tag{1}$$

where *C*AR is aqua regia extracted metal (mean values of the 3 years of the study) (mg kg−<sup>1</sup> soil). *C*AR *ref* is the background reference element concentration in uncontaminated areas (mg kg−<sup>1</sup> , PTE concentrations in pristine soils, "world soil average") [37].

The CF follows the classification: class I: CF < 1 (pristine soil), class II: CF = 1–3 (contamination is rated "moderate"), class III: CF = 3–6 ("considerable"), class IV: CF > 6 ("very high").

The geo-accumulation Index (Igeo) is as follows:

$$\text{I}\_{\text{geo}} = \log\_2(\text{C}\_{\text{AR}} \mid 1.5 \,\text{C}\_{\text{AR} \,\text{ref}}) \tag{2}$$

Igeo [39] is designated using Latin numbering [40]: class I: Igeo < 0, class II: Igeo = 0–1, class III: Igeo = 1–2, class IV: Igeo = 2–3, class V: Igeo = 3–4, class VI: Igeo = 4–5, class VII: Igeo > 5.

#### *2.5. Statistical and Geostatistical Procedures*

The 3-year data taken from the 62 sample sites were subjected to a one-way ANOVA at a *p* < 0.05 level using the SPSS-26 package, in order to identify potential timescale changes [23,41]. The results showed no significant differences for the mean three-year values, although an upward trend from year to year was detected (ANOVA results are listed in the Supplementary Materials Section). The mean values were used for further analysis and the construction of thematic maps. Initially, the mean values were checked for normality and homogeneity by means of the Kolmogorov–Smirnov and Shapiro–Wilks tests [42]. The results showed an abnormal distribution. Based on that distribution, a spatial interpolation process using the ordinary kriging along with data logarithmic transformation and a grid of 5 m spatial resolution, was performed in the SAGA 7.9.0 GIS (SAGA User Group Association, 2020) [43]. Prior to kriging, the datasets were also checked for anisotropy in two perpendicular directions. Showing a similar spatial autocorrelation pattern, an omnidirectional variogram was computed for every variable, and a suitable model was fitted to that. The geospatial statistics are presented in Tables 3 and 4. The ordinary kriging process resulted in raster images of pollution for each element. These rasters, along with the sampling point vector information and local points of interests (POIs), were integrated on the thematic maps using the QGIS 3.16.1 (QGIS Development Team, 2020) [44] and utilizing the web map service OpenStreetMap as a base layer (Figure 3).


**Table 3.** Descriptive statistics for the contamination factor (CF) of the potentially toxic elements and geostatistical parameters, for the creation and validation of the contour maps.

**Table 4.** Descriptive statistics for the Igeo of the potentially toxic elements and geostatistical parameters, for the creation and validation of the contour maps.


#### **3. Results**

#### *3.1. Physicochemical Properties of Soil Samples*

The soils ranged from sandy loam to clayey, with soil reaction values between 6.6 and 8.9 (Table 1). High clay content was found in 19% of the soil samples, while 57% had a sandy texture. According to the soil pH values, 82.3% of the soil samples were alkaline, while 17.7% had a slightly acid soil reaction (<7). Organic matter (OM) values ranged between 0.3% and 4.4%. Soil electrical conductivity (EC) varied between 1122 and 6958 µS cm−<sup>1</sup> , with the mean value at 3204 µS cm−<sup>1</sup> .

#### *3.2. Levels of Potentially Toxic Elements*

The mean values of metal concentrations during the three years' research along with their statistical characteristics are presented in Table 2. The mean values were found to be lower than maximum permitted values. The mean Cu concentration was 53.9, while the CEC directive value [36] was 140 mg kg−<sup>1</sup> . The other values were as follows: for Zn, the mean was 131.0, while the CEC limit was 300; for Pb, the mean was 35.9, while the CEC limit was 300; for Cd, the mean was 0.94, while the CEC limit was 3; for Cr, the mean was

93.2, while the Kabata–Pendias [37] value was 200; and for Ni, the average was 67.3, while the CEC limit was 75mg kg−<sup>1</sup> ). Average Ni concentration was lower than the EU limit, while the Co and Mn mean concentrations were higher than the background levels.

#### *3.3. Geostatistical analysis, Construction of Thematic Maps of Availability Indices of PTEs*

The thematic maps, presented in Figures 2 and 3, were created for the optimal depiction of spatial variability of the metal pollution assessment. The contamination factor (CF) and geo-accumulation index (Igeo) were used as useful tools to assess the pollution of the area. Descriptive statistics for the contamination factor (CF) and geo-accumulation index (Igeo) calculations along with the geostatistical parameters, for the creation and validation of the contour maps, are presented in Tables 3 and 4, respectively.

 **Figure 2.** Thematic maps of potentially toxic metals based on contamination factor (CF) values in Volos (mean values of the three years of the study).

1

**Figure 3.** Thematic maps of potentially toxic metals based on the geo-accumulation index (Igeo) values in Volos city (mean values of the three-year study).

2 The largest number of soil samples had CF values belonging to class II (1–3) and were characterized as "moderate contamination". In class II, the order of CFs was as follows: Cd = Zn > Ni > Co > Cu = Pb = Cr = Mn. In class III (CF: 3–6), characterized as "considerable contamination", the order changed as follows: Ni > Co > Zn = Cd. There were no soil samples with CF > 6, i.e., belonging to class IV with "very high contamination".

The Igeo values were also calculated. In class II (0–1), the order of the metal elements was as follows: Cd > Zn > Cr = Ni > Pb > Co > Mn > Cu. In class III (1–2), Ni was the most enriched element in the soil, followed by Co. Zn and Cd exhibited lower Igeo values in the decreasing sequence: Ni > Co > Zn = Cd. None of the metals studied had Igeo values > 2.

1

#### **4. Discussion**

#### *4.1. Physicochemical Properties of Soil Samples*

The soil samples were neutral to slightly alkaline. The presence of alkaline soils in the urban environment is rather common [7]. High clay content has usually been found in urban and peri-urban locations [28]. On the other hand, 19% of the samples had a sandy texture, which is probably due to the proximity of the samples to the coastal zone [45] or to the area near the river, which flows into the Pagasiticos Gulf. The organic matter content was rather low, but usually observed in Mediterranean soils with the characteristic climate, relatively high temperatures, and low humidity [46]. Anthropogenic activities create conditions of impaired salinity in soils, and especially in the capitals of the prefectures. A survey conducted in a Mediterranean site also revealed that soil salinity is controlled mainly by seawater intrusion and the usual application of salt-rich water [12]. The lower part of the study area was adjacent to the edges of the Pagasiticos Gulf, which naturally contributes to the increase in salinity [47].

#### *4.2. Levels of Potentially Toxic Elements and Contamination Indices*

The average CF values of the metals were >1, which indicates the strong anthropogenic influence in the area studied. The PTE contents were lower than those in soil samples in Athens, as described by Massas et al. [18] and Kelepertzis and Argyraki [5]. The soil samples seemed to be less contaminated than soils in areas with high anthropogenic activities [26,30,31,48]. The previous studies of Antoniadis et al. [48] and Kelepertzis et al. [30] made extensive reference to the impact of industrial activity on the peri-urban environment of Volos. The soil contamination indices (CF and Igeo) indicated low or moderate contamination. The elevated Ni Igeo may have been influenced by activities leading to enhanced Ni deposition [49]. Moreover, high Ni levels were previously reported in the nearby site, probably due to the geochemical composition of the parental material [48].

#### *4.3. Construction of Thematic Maps and Geostatistical Analysis*

The constructed maps based on the contamination indices [50] seemed to be a useful tool for the investigation of spatial distribution of metals in the city of Volos. To avoid misleading interpretations and alerts from visual observations, a uniform full-scale of the CF and Igeo values was used for the present maps' color classification. Despite reducing the spatial discrimination, this choice resulted in the maps depicting the potential toxicity of the studied metals [23]. In the present study, the cubic, exponential, and spherical models were used for the construction of the thematic maps in Figures 2 and 3, of the CF and Igeo, respectively. The cubic model seemed to be the proper one for the construction of maps based on the contamination factor of Cu along with the Pb, Ni, and Cr maps based on the Igeo indices [24]. The exponential and spherical models were preferred for the spatial distribution of Zn and Co in the CF and Igeo indices, respectively. Similar techniques and mathematical models have been used by other researchers [51,52] as they provide the best predictions. The low nugget/sill ratio (<0.25) obtained for most of the parameters indicated the spatial dependence of the variables [53].

#### 4.3.1. Maps of Contamination Factors

The Cu, Pb, Mn, and Cr contamination factor maps (Figure 2) revealed that the highest concentrations of these metals occurred near the heaviest traffic roads (Analipseos, K., Kartali, Demetriados), as well as in the areas adjacent to the railway station, the bus stations, and the commercial port, as mentioned in many studies [13,30,54]. The long-time and increased vehicle traffic using lead-based fuels, before the banning of lead-containing gasoline, may have led to increased PTE contents in urban soils [54,55]. Nevertheless, the CF values were high, and the ranking led to the rating of "moderate contamination" conditions. Cd and Zn revealed "moderate contamination" in all sampling points, with the exception of S40. The higher Cd concentrations were probably because of road construction works carried out during the sampling period (second year, 2019). The high CF value of Zn

in S1 near the commercial port led to "considerable contamination". Prolonged use and storage of metal objects and scrubs may have led to high Zn levels locally [29]. Ni and Co concentrations led to "considerable contamination" around the railway, the bus stations, and the port [56]. However, in the main area of Volos, the CF of Ni and Co indicated a "moderate contamination".

### 4.3.2. Maps of Geo-Accumulation Index

The maps (Figure 3), based on the Igeo distribution, revealed moderate contamination of Cu and Mn at the main roads and railway, and beside the bus stations and the city port (Igeo, 0–1), and the values suggested lower contamination level compared to other researches [28,29,37]. The Igeo values of Cd and Zn also belonged to class II (0–1), but in the S1 and S40 soil samples, an exaggeration was observed as higher values were recorded. High levels of metals beside ports are usually observed [56]. It is well known that metal concentration in aquatic ecosystems, such as ports, is measured in water and in sediment. Usually, metals exist in the lowest content in water and in higher concentrations in sea sediments [6,12]. The permanent presence of metals near ports is a problem of utmost concern as their constant contact with salts and the particularly high percentage of moisture creates oxidizing conditions. The metals therefore corrode and their concentrations in the soils accumulate [45]. In 16.1% of the soil samples, the Igeo values of nickel belonged to class III (1–2). Nickel has a large number of industrial uses and applications [57]. The occurrence of unexpectedly higher concentrations of nickel at the railway station and on the west side of the study area may be related to its geochemical origin, as recorded in a previous study [48]. The Igeo values of class II (0–1) were found for Pb, Co, and Cr near the main city roads, possibly affected by heavy traffic [57–60].

The factories located around the city of Volos may contribute to the pollution of the area, especially when the weather conditions are favorable for the spread of pollutants [30,48]. It is also worth mentioning that Volos is a rural area surrounded by the Pelion mountain to the north and the Pagasiticos Gulf to the south, creating a microenvironment favorable to the frequent entrapment of atmospheric pollution [29,30]. The burning of materials in the fireplaces of the houses to the right of the study area may also contribute to the increase in PTE pollution [61]. There is also evidence that the coal or crude oil used for household heating, gasoline, and diesel vehicle exhaust is of the major problems in urban soil pollution [29]. Concluding, the contamination seems to be significantly increasing near the commercial port of the city of Volos, close to the railway station, and in the intercity bus station along with the big and heaviest traffic roads of Volos [62]. The present study could give rise to further and more thorough research in the future. The identification of contamination sources at the local and national level could control or even reduce the levels of PTEs in the area of interest.

#### **5. Conclusions**

The level of PTEs in the soil samples of Volos were found to be lower than the maximum permitted values set by the European Union. Soil contamination indices (CF and Igeo) also confirmed the low to moderate level of contamination with PTEs. The formation of a database with three-year records of the PTE level is an innovation, as it has not been described before. Map formation, based on the values of the contamination indices, is a valuable, useful, and easy-to-handle tool for recording, monitoring, managing, and predicting contamination levels in the area of interest. Nevertheless, continuous long-term monitoring for the assessment of PTE levels is of paramount importance, as high PTE levels could ultimately have negative effects on the environment, food, and consequently, on the health of the citizens of Volos and the adjacent communities. Spatial environmental monitoring could be very useful for the assessment of environmental conditions and trends in the cities of Greece, and even worldwide, as it leads to an evidence-based approach for reporting to national policy makers. After all, maintaining good urban soil environmental quality is crucial for several socio-economic reasons along with human health and well-being.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/2071-105 0/13/4/2029/s1, (Table S1(SM) Pseudo-total concentrations of heavy metals of the reference material BCR CRM 141 R, Table S2(SM) Results of the analysis of variance (ANOVA) for the three-year data.

**Author Contributions:** Conceptualization, E.E.G. and N.G.T.; methodology, E.E.G. and C.C.; software, C.C.; validation, E.E.G. and N.G.T.; formal analfiguysis, E.E.G., S.G.P. and C.C.; investigation, E.E.G. and S.G.P.; resources, E.E.G., C.C. and N.G.T.; data curation, E.E.G., C.C. and S.G.P.; writing original draft preparation, E.E.G., S.G.P. and C.C.; writing—review and editing, E.E.G., C.C. and N.G.T.; visualization, C.C. and S.G.P.; supervision, E.E.G. and N.G.T.; All authors have read and agreed to the published version of the manuscript.

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

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data that support the findings of this study are available from the corresponding author upon reasonable request.

**Acknowledgments:** The authors would like to express their gratitude to the Postgraduate Program of the University of Thessaly, entitled: Environmental Changes Management and Circular Economy. Part of the results come from the post-graduate dissertation of student Sotiria G. Papadimou.

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

#### **References**

