1. Introduction
Although soil is acknowledged as a finite, non-reproducible resource [
1,
2], its utilization for new urbanization processes is still happening at fast rates [
3,
4,
5]. Furthermore, the soil is the most important pool of biodiversity [
6]. It provides a large variety of ecosystem services (supporting, regulative, provisioning, and aesthetic/cultural) [
7,
8,
9] while being a primary resource to be considered in adapting cities against multiple hazards [
10,
11].
Recent studies have claimed how the deep knowledge of soil characteristics (pedogenesis) could inform decision-making during urban planning [
12,
13]. It has been demonstrated that if well estimated, the soil characteristics play a vital role in diminishing the vulnerability of socio-ecological and technological systems [
14]. If new urban areas are built without consideration of soil characteristics (e.g., soil texture, composition, slope, depth, capability, and hydraulic conductibility), the probability is that urban systems will become more sensitive to multiple hazards (e.g., earthquakes, liquefaction, cloudburst, drought, landslides, urban heat island) [
15,
16,
17].
Moreover, some recent works on urban ecosystems show how soil typology significantly affects ecosystem performance in urban areas [
18]. To clarify this concept, it is essential to understand how the preliminary knowledge of soil can play a vital role both in limiting the ecological impacts in the processes of sealing (adding new urban areas) [
19] or, on the other hand, maximizing the biophysical performance while de-sealing [
20,
21]. For example, suppose an urban system is already vulnerable to pluvial flooding (cloudburst) [
22]; the impermeabilization of new areas on soil with good drainage can dramatically impact it. In contrast, sealing the same quantity of land in soil with limited infiltration can reduce the impact [
23,
24,
25]. The same is valid for the contrary situation: the de-sealing of already urbanized land is strongly influenced by hydraulic conductibility. Some modeling results show, for example, that when the de-sealing process occurs in soils with bad hydraulic conductibility, the results in terms of run-off containment can be dramatically limited or not perceptible at all [
26]. That is why soil knowledge can support so-called “performance-based solutions” [
27,
28,
29], which indeed can guide decision-makers while understanding exactly where the de-permeabilization can achieve its maximum efficiency in run-off containment [
30,
31]. Therefore, urbanization and greening processes can generate different degrees of vulnerability to urban systems.
The same concept can be applied to other hazards, such as earthquakes [
32,
33]. The construction of new anthropic areas near fault lines can exponentially increase the possibility of structural collapses and damages during an earthquake event with surface ruptures and seismic waves [
34]. Parameters such as fault lines, groundwater, lithology, and hydraulic conductivity are essential for cities in first-degree earthquake zones located in coastal areas. The damage and loss of life caused by the earthquake waves of the Samos earthquake, the epicenter of which was about 70 km away from the city center of Izmir in 2020, has been revealed in many studies. However, this earthquake caused a lot of damage even though it was far from the city center of Izmir. As a result of the rapid urban transformation campaign in the last two years, many buildings have been demolished, and new buildings are being built in the same places or newly opened settlements. In addition to the importance of parameters such as the structural features of the buildings and the number of floors, the main issue to be considered is the choice of location.
Another example can be applied to the reduction in biodiversity: the knowledge of soil chemical composition and physical structure can support green infrastructure design that maximizes the ecological connections in urban catchments [
35,
36].
With the introduction of advanced GIS techniques, land-use change (LUC) assessment becomes a fundamental tool to support analyzing the causes and consequences of land-use dynamics [
35]. LUC can be studied using geoprocessing treatment of digital spatial data and is applied in scientific research on land-use environmental effects at different scales [
36,
37,
38,
39]. Right after 1980, when the European Union launched a program aimed at monitoring land use and land-use change named “Corine Land Cover” (CLC), there was a flourishing utilization of LUC assessment for environment observation.
The 2004 special issue of the Journal of Environmental Management dedicated to “modelling land use change” (Veldkamp and Verburg 2004) demonstrated that after CLC implementation, LUC has been strongly employed in the academic and scientific literature to create site-specific land pattern analyses. LUC can be studied by geoprocessing analysis based on two methods:
The “differential” method, which allows for comparing statistical datasets of urban land uses in different time thresholds and extracting the absolute and relative variations [
3,
40,
41].
The “flows” method which is based on cross-tabulation analysis (matrix) of land uses between different time thresholds and extracts all possible land-use “substitutions” [
42,
43,
44].
In the study by Verburg et al. [
45], the authors show how the LUC models based on the flow method are useful to empirically detect the trend, variations and alterations of land which can generate different impacts on the urban landscape [
35,
46,
47]. Verburg concludes that among others, three factors of LUC analysis can lead the development of supporting decision-making tools: (i) the multi-scalar approach to analysis; (ii) the need for proximity/interaction of data within the context of the study, and (iii) the temporal dynamic as a fundamental approach to measuring fluxes rather than stock data. According to these factors, within this work we wanted to use a multidisciplinary, innovative approach, that allowed us to (i) deal the multi-scalarity of multiple datasets produced at different acquisition scales, (ii) address the interaction between local and metropolitan data while intersecting LUC features with different land characteristics and (iii) employ a cross-tabulation analysis approach that helps us to understand beyond potential future trends in land changes. Once extracted, the new urbanization features were superimposed with different soil layers, and then the soil characteristics were evaluated in the light of each anthropization process. Within this process, we wanted to check whether the urbanization process between 2012 and 2018 in the city of Izmir (Turkey) has taken place in favorable or unfavorable soils while employing spatial land-use change (LUC) analysis based on different geoprocessing phases [
48].
Even though a lot of studies have already explored Izmir’s characteristics of land [
49,
50], soil [
51,
52,
53], and plans or projects [
54,
55,
56], in this work these different domains are analyzed simultaneously and by interaction, thus furnishing a comprehensive picture of how geology information is used to steer land transformations in cities located in first-degree earthquake zones and coastal areas.
Our research hypothesis is that the quantity of information related to soil and geology characteristics nowadays are overall widespread and publicly accessible, but their systematic utilization in planning is still negligible and scarce, thus affecting the environmental sustainability and the vulnerability of urban areas to different typologies of hazards. Our assumption is that revealing the impact of past land-use changes through soil and geological features can shed some light on the importance of including this knowledge in the land transformation process, especially forwarding a more environmentally compatible urban system.
The paper is structured as follow: after the introduction (
Section 1) in
Section 2 (Material and Method) we introduced the catchment and defined the dataset and geoprocessing phases that shaped our geographic information system-based (GIS) methodology. Then, in
Section 3 (Results), we analytically reported the empirical results and some examples of our GIS analysis. In
Section 4 (Discussion), we discussed the results in light of the potential land-use policy implications. We outlined a list of possible actions to reduce the vulnerability of future urban systems, including a broad utilization of land-use change analysis to support landscape policies. Finally,
Section 5 (Conclusions) recalled the research question and briefly synthesized this study’s main innovations.
2. Materials and Methods
2.1. The Study Area
The İzmir Province is located in the Aegean region and spans 11,878 km
2 (
Figure 1). It is located on the western coast of Turkey while creating a peninsula in the Aegean Sea. The altitude ranges between 0 and 1200 m, and the climate is typically the Mediterranean, with a hot and dry summer and a rainy winter.
The average yearly precipitation recorded between 1991 and 2000 amounts to 675 mm, and is most concentrated during the wet, fall season. The average monthly temperature ranges from 3.8 degrees (°C) in January and 20.8 degrees (°C) in July and August [
57].
The landform is heterogenous, with high mountainous peaks directly surrounding Izmir’s gulf, where the dense urban settlement is distributed.
In recent years, urbanization occurred at the expense of plain, fertile agricultural areas (26 thousand hectares). Nonetheless, natural and seminatural Mediterranean soils (more than 10 thousand hectares) were also urbanized in the last 28 years (1990–2018) [
58].
The Aegean region, also located east of the Mediterranean earthquake belt, is one of Anatolia’s riskiest regions regarding earthquake potential [
59]. İzmir, the largest city in Western Anatolia in terms of population (4,425,789 in 2021) [
60], has been affected by earthquakes in its history and has witnessed a loss of life and property. When the historical and instrumental period records are examined [
61,
62,
63,
64,
65], in addition to the earthquakes that took place within the provincial borders of İzmir, many earthquakes in the surrounding provinces, basins, and even the Aegean Sea caused damage to the city. As can be seen especially in the 2020 Samos earthquake, damages occurred as a result of liquefaction in the alluvial soils in the structures in the city center, and in some districts, the earthquake intensity affected the buildings more strongly for a longer period of time with the amplification effect in the earthquake waves due to the basin structure and groundwater levels. For this reason, it is essential to consider seismicity and proximity of existing fault segments when examining soil properties in terms of urbanization. To properly analyze the region, we employed different datasets presented hereafter.
2.2. Data Structure: Urban Atlas
We employed the Urban Atlas database [
40] to analyze the most recent urbanization trend in Izmir (
Figure 2). The dataset results from a joint initiative of the Commission Directorate-General for Regional and Urban Policy and the Directorate-General for Enterprise and Industry in the EU Copernicus program (European Space Agency and the European Environment Agency) [
66].
Urban Atlas provides open-access pan-European comparable land-cover and land-use data that include all major functional urban areas (FUA) [
67,
68], including Turkey and the West Balkans. The Urban Atlas 2012 release contains 785 FUAs with more than 50,000 inhabitants, and the data can be downloaded in vector format (shapefile) together with the metadata for all the necessary information [
69]. Urban Atlas is mainly based on combining (statistical) image classification and visual interpretation of high-resolution satellite imagery. Multispectral SPOT 5 and 6 and Formosat-2 pan-sharpened imagery with a 2 to 2.5 m spatial resolution is used as input data. Built-up classification was conducted while employing the imperviousness layer and enriched by ancillary local sources [
19,
66]. The legend is articulated in 27 land-use and land-cover classes, and the minimum mapping units range from 0.25 ha for class one to 1 ha for classes two, three, and five (see
Table 1).
2.3. Data Structure: Soil Hydraulic Conductibility
The second dataset we employed is the map of saturated hydraulic conductibility (Ksat, mm/h, see
Table 2), defined as the soil’s ability to be vertically crossed by fluids according to the fundamental lithological, hydrogeological and geological data [
71].
Hydrological soil are generally examined in four main groups according to their approximate infiltration rate, considering their flow and run-off and flooding potentials with their vegetation and lithological characteristics [
72].
Group A: Soils (from over-drained sands or gravelly sands) have a high infiltration rate and permeability (low run-off potential) during precipitation.
Group B: Soils with moderate infiltration rate and permeability when saturated with groundwater. These mainly correspond to drainage units but contain more lithologically different materials.
Group C: These are generally impermeable, poorly conductive soils (clay content is present) with a slow infiltration rate when saturated with groundwater.
Group D: This group consists of almost impermeable, shallow soils on material with very high clay content. When saturated with groundwater, it has a very slow infiltration rate (high flow potential). The water transmission of these soils is prolonged.
Hydraulic conductibility is a function of soil porosity and, in turn, soil texture [
73,
74]: clay soils generally have a lower conductibility than sandy and gravel soils [
75] (see the relationship between geological units and hydrological characterization in
Table 3). Quaternary alluvium deposits (recent alluvium, talus, delta, and shoreline deposits) have the highest infiltration rate and lowest run-off potential because of included sand and gravel materials. With increasing clay content and small grain size, volcanic, volcano-sedimentary, and flysch units have relatively slow infiltration rates, so this situation causes unexpected surface run-off and flooding events. In addition, similar problems in the soils will be formed by the alteration and physical deformation of these units, which are predominantly in rock type (carbonates, andesite, dacite, etc.) in the field.
As it is known, there is an essential relationship between hydraulic conductivity and groundwater. Whether the soil is saturated with water and its relationship with the slope are parameters that should be considered in terms of engineering applications. Within the scope of this study, topographically inclined and high regions are observed as the areas where rock units are dominant. In addition, stream beds in these regions have been rehabilitated or taken into reinforced concrete channels. This situation directly affects the infiltration and conductivity processes.
Saturated hydraulic conductivity is significantly affected by surface infiltration on the slope. When the saturated soil hydraulic conductivity is high, the safety factor is significantly reduced as surface infiltration appears rapid [
75]. Permeability is a measure of the ability of water to pass through units. In situ saturated hydraulic conductivity means implicitly saturated permeability, considering the trapped air in the soil. This value is considered more appropriate than a valid unsaturated region’s fully saturated permeability coefficient.
For this reason, we see that the permeability values of Quaternary alluvial units may differ (detailed soil classification maps according to sand, clay, and silt material size have not been prepared yet), but the presence of water-bearing capacity is seen from the existing wells and groundwater levels. In volcanic (andesite, dacite) and carbonate units, the permeability parameter comes to the fore rather than hydraulic conductivity. Also, fracture-crack systems and faults can cause high permeability and conductivity values in these units. In
Table 3, these rock units were included in the classification by considering the characteristics of the altered and relatively low-thickness soil formations associated with the rock. In addition, when the relationship between seismicity and groundwater is analyzed based on rock units, there is no liquefaction in rock units. In this context, it seems that the Quaternary-aged alluvial deposits are the part that needs to be examined and paid attention to primarily for the province of İzmir.
2.4. Data Structure: Major Fault Segments
The active fault map database of Turkey on a 1:250,000 scale [
76] was revised as the last version of the tectonic structures within the country’s borders. In that study, active faults were used for faults that formed at least one surface rupture and offset in the Quaternary in geological times. Earthquake surface rupture, Holocene fault, and Quaternary fault definitions corresponding to the Quaternary time were mapped as neotectonic period faults indicating younger activity in geological times. Other faults, which are older and thought to be inactive, were evaluated as probable Quaternary faults or lineaments [
76].
Although there are no sharp definition boundaries in all these fault types, there are some limits in applications of real life. It is sometimes a common idea in academia and engineering practices that inactive faults are risk-free. However, mapping these faults with some morphological anomalies, lineaments, contours, and surface manifestations was possible. In addition, surface fractures/ruptures that occur after the earthquake can quickly disappear with surface erosion and atmospheric conditions in a short time. In other words, it cannot be determined today (unless paleoseismological studies are carried out) that these faults may have created surface ruptures in the past. For this reason, these undefined and undifferentiated faults should still be considered risk factors.
Furthermore, it is commonly thought that faults transfer their existing energy by triggering other fault segments near them during/after earthquakes. In this context, it is known that inactive fault segments already have deformation and fracture zones, and that these weak zones will become active again in the event of an energy–stress transfer, and this situation should not be ignored. For this reason, the active–inactive fault distinction was not made in this study, and it was considered that all kinds of faults and fractures could be re-activated in engineering studies.
The essential fault segments around İzmir and their types are given in
Table 4. While most faults have a normal fault character, others have a strike-slip character. The faults, which are thought to be located in the Izmir metropolitan area, the gulf, and especially the alluvial ground, have not been extensively investigated yet. These studies will be detailed in the coming years with the microzonation projects that started in 2022 under the execution of the İzmir Metropolitan Municipality. In this context, the critical fault segments (
Table 4) found in the 1/100,000 (personal data for earthquake prevention/avoidance zone example) and 1/250,000 scale maps are indicated in this study, and it is suggested they be considered in soil interactions for future studies with more detailed scales.
2.5. Data Structure: Land Capability Classification
A third dataset we employed is the map of land capability classification (LCC, see
Table 5), defined as the soil’s propensity to be used for agricultural purposes [
77,
78]. LCC is based on the 1/25,000 scale national soil database of Turkey based on large soil groups.
According to the physical and chemical characteristic, the LCC divides the soil into seven classes, ranging from one—a high propensity to use the soil for agricultural production to eight—no propensity for agricultural production [
69,
79]. According to the most common interpretation of the LCC, classes that range between one and four are considered generally “good” soils in terms of potential agricultural production (no specific limitation, from good to medium characterization). Classes five to seven represent soils with some or severe limitations for agricultural purposes, and while class eight presents several restrictions for productive purposes. From a perspective of sustainability, soils one, two or three are extremely precious [
80,
81]: they are the plainest, most fertile, deep, and biodiverse soils. At the same time, they are more favorable to being urbanized for the same reasons [
5,
82]: agricultural land represents the most suitable and cheap solution for new urban expansion. Therefore the trade-off between farm productivity and demand for urban land is always conflicting [
83,
84,
85].
2.6. Data Processing
The land-use change analysis was performed by using ESRI ArcGIS (ver. 10.8.1) [
86]. The geoprocessing phases included the intersection of vector polygons UA 2012 and 2018. Once intersected, the urbanization process was characterized by all the transitions from agricultural and natural or seminatural classes into urban areas (including the new urban green areas and green/sports facilities) [
43].
Once extracted, the quantitative evaluation of the urbanization was processed employing MS-Excel and a land-use change matrix was created using the pivot function. Then, the empirical analysis of the quantity and distribution of the most relevant changes was completed.
The impact of urbanization on soil hydraulic conductibility was processed by a spatial intersection between the features of the land-use changes with the hydrological soil group classification (
Table 3). Then the areas were recalculated, and the statistics were extracted using a MS Excel pivot table and analyzed.
Analogously, the impact of urbanization on LCC was processed by a spatial intersection between the features of the land-use changes with the large soil group map (
Table 5). The land capability classification layer only partially covered the land-use and land-cover catchment; thus, the spatial intersection between the land-use change and the LCC map accounted for only 73% of the features. Statistics were analyzed and extracted using a MS Excel pivot table.
The complete list of layers used in this research is reported in
Table 6, while in
Figure 3, there is a synthesis of the conceptual workflow.
3. Results
3.1. Land-Use Changes in the Metropolitan Area
Between 2012 and 2018, the process of urbanization occupied 6074.33 ha.
This process mainly happened at the expense of (see
Figure 4):
Arable land (annual crops) (2592.38 ha);
Herbaceous vegetation associations (natural grassland, moors...) (1498.71 ha);
Pastures (1171.06 ha);
Permanent crops (vineyards, fruit trees, olive groves) (491.35 ha);
Forests (239.85 ha).
Specifically, arable land (annual crops) has been mainly transformed into industrial, commercial, public, military, and private units. Herbaceous vegetation associations have been mainly transformed into construction sites. Pastures have been mostly converted into industrial or commercial areas, and herbaceous vegetation associations have been mainly transformed into mineral extraction and dump sites.
Industry or commercial areas, construction sites, or quarries/dump sites occupy the most abundant part of the urbanization process, which happens in agricultural or pasture land and herbaceous vegetation.
Figure 5 focuses on three different land take typologies. An abundant, spatially clustered amount of new transformation sites are distributed in the Menemen neighborhood (Gediz Delta urbanization
Figure 3A). Here, the construction of a new primary road network has been accompanied by the exponential growth of many residential blocks. In this northern cluster, the new construction sites occupy, on average, more than 4 ha. In comparison, the new residential areas occupy less than 2 ha, the average dimension of the new industrial and commercial sites.
Another spatially clustered land take concentration is distributed in the plain eastern corridor of the city, where, similarly to what happened in Menemen, the primary road construction created the precondition for higher accessibility of the industrial zone in the neighborhood (
Figure 3B). Here, the average dimension of construction sites is much higher (50 ha), and is dominated by the saturation and expansion of the vast industrial area of the city. The average size of the new production plants is 4 ha, thus emphasizing the massive impact of this extended kind of urbanization.
Finally, a third land take’s spatial cluster was recognized in the south-eastern corridor of the metropolitan area (
Figure 3C). Here, there is a mix of different land take patterns: a diffuse expansion of new high- and low-density residential built-up areas surrounding the city of Torbali with an average dimension of 1.5 ha, while there is a concentration of new industrial and commercial areas in the northern catchment, where the average dimension is over 5 ha.
Overall, the spatial characterization and distribution of the land take recorded during the five years demonstrate the heterogeneous character of the phenomena: the polarization of new residential and industrial/commercial areas often happened along the construction of new primary road networks, which, in turn, generated a higher real-estate attractiveness for land transformation. The size and dimension of the new urbanization varies according to the position, proximity with the existent settlement, and land take typology. Nevertheless, it is noteworthy that the major urban expansions are placed on the outskirts of the existing city, at the interface between urban and rural or natural land.
3.2. Impacts on Soil Hydraulic Conductibility
The land-use change impact matrix with hydrological soil group classification clarifies that more than 53% (4076 ha) of the total land take happened on A-type soil, thus causing the most problematic impact in terms of run-off control. It is followed by 13% for B-type, 14% for C-type, and 20% for D-type (
Table 7).
Impressively, the land-use classes that had the most significant impact on A-type soil is the industrial, commercial, public, military, and private units. More than 1673 ha of agricultural, natural, or seminatural soils with the highest hydrological capabilities in water absorption were urbanized and sealed mainly through new industries or commercial centers. Analogously, more than 800 ha of A-type soils were urbanized by new construction sites, and this means that the trajectories of transformation sites critically intercept the most hydraulically conductible. Finally, more than 356 ha of A-type soils were urbanized by new mineral extraction and dump sites.
It is worth mentioning that among the new urban development for new residential areas, the discontinuous dense urban fabric had the highest impact on A-type soils (289 ha). Again, this means that most of the new peripheral residential areas built during the last years are placed on highly hydraulically compatible soils. Furthermore, this type of soil degradation also affects the coastline while generating seawater intrusion and the consequent salinization of groundwater [
54,
56].
3.3. Impacts on Land Capability
The land-use change impact matrix with LCC shows that of the total amount of land taken by 2012–2018 which intersected the LCC dataset (4414 ha), 14.47% (638 ha) was in class one, 19.20% (847 ha) in class two, 9.24% (407 ha) in class three, 6.69% (736 ha) in class four, 9.13% (402 ha) in class six, 12.45% (549 ha) in class seven, 0.49% (21 ha) in class eight and 18.34% (809 ha) was located in urban land (no land capability classification, see
Figure 6).
In absolute terms, the class two soils were abundantly affected by urbanization, followed by urban soils and class four soils. It is worth mentioning that more than 638 ha of urbanized land was in class one, representing the most optimal and suitable location for agricultural purposes.
Looking at the land-use change impact matrix on LCC (
Table 8), it can be recognized that 406 ha of the land urbanized on class one soils were isolated structures, while 374 ha were on class two soils. This means that the new detached, peripheral structures located on the city’s border, without continuity and contiguity with the adjacent built-up land, are located on the better soils in terms of potential productivity, thus exacerbating the trade-off between competitive uses of the soil (urban vs. rural/productive) [
85,
87,
88].
Another relevant category is the mineral extraction and dumps sites, which occupied more than 178 ha of class two soil and 342 ha of class four soil. This result must also be interpreted in light of the landscape impact of extraction sites, which according to this data, are located in averagely good land in terms of potential productivity.
3.4. Impacts on Geological Units and Fault Lines
Figure 7 shows the distribution of urbanization changes (red polygons) on geological data. According to this data, most new areas are located on alluvium, volcanic units, and carbonates. In addition, newly built settlements are concentrated in the region between the İzmir Fault (IF) and Karşıyaka-Bornova faults (KBF). Faults, fractures, and crack systems also affect the hydrological characteristics of the units in this area. Especially in the areas of volcanic and carbonate rock, fissure and crack systems may accelerate the infiltration from the surface in the soil and the rock units below. For this reason, it is thought that there is infiltration at medium-slow levels in carbonate and volcanic units, mainly due to the fracture and fissure systems seen in field studies.
There are vital N–S, NE–SW, NW–SE, and E–W trending fault segments in the city of İzmir, located in a seismically very active region (
Figure 8).
As seen in
Figure 9, the regions with the lowest V
s30, the average shear wave velocity (V
s30: usually the average seismic shear wave velocity in the upper 30 m), coincide with the distributions of Quaternary alluvial soils around Bayraklı, Karşıyaka, and Konak [
52]. It is known that the basin/alluvium effects before, during, and after the earthquake in İzmir Bay and its surroundings cause more than expected deformations on the structures with groundwater level changes and soil amplification [
52]. As in the location distribution of the damaged buildings in the 2020 Samos earthquake, tremors affect rock masses and alluvium units differently on either side of the fault planes. For example, buildings in the hanging wall of normal faults are subject to more severe stress than the locations in the footwall. This situation is critical for buildings where the fault surface deformation has the maximum effect on the fault surface rupture. For this reason, some of the literature suggests earthquake prevention/avoidance zones [
90,
91].
4. Discussion
4.1. Soil Discussion
According to the analysis of land-use change, urbanization is happening regardless of the soil characteristics and the associated seismicity risk. Urbanization is constantly growing, especially in alluvial soils and where groundwater is near the surface.
As demonstrated by this comprehensive analysis, the shrinking of arable land is due to the uncontrolled diffusion of new infrastructure and commercial, productive, and residential areas, which creates spatial agglomerations of different sizes and forms while exposing the new settlements to vulnerable situations [
93,
94,
95].
The massive, uncontrolled utilization of land for the above mentioned uses is problematic for the following reasons [
96,
97,
98]:
The first is that without a proper metropolitan structural scheme for urban development, the massive diffusion of new anthropic areas is going to seal the most fertile soils in terms of potential biomass production while increasing, in the long term, the problem of food security at the local level;
The second is that, according to results, the reduction in biodiversity is accompanied by a significant impact on highly hydraulically conductible soils (A-type), which, in turn, generate a high risk of pluvial flooding in the newly developed urban areas [
99];
The third is that the reduction in the most fertile and hydraulically compatible soil is shrinking the cropland managed by extensive sustainable farming. The intensification of crop production in the remaining plain land is causing a change in the landscape management of rural land.
For the first two points, ecosystem vulnerability classifications based on soil characterization can be widely used in different study contexts to set policy targets and define site-specific regulations [
100,
101,
102]. By following comprehensive environmental recommendations, soil and geology characteristics can be used to determine new conservation areas through specific regulatory zoning and local planning activity.
If planned, land take should be limited, mitigated, or compensated [
103,
104]. The impact of land transformation on the ecosystem can be minimized while adopting green mitigation or compensative solutions. The new design and nature-based or technological solutions can drastically reduce the impermeabilization of the soil, increase the green areas, the shadow, and the vegetated biomass;
Compensation measures for the land take process can be considered at the landscape/metropolitan scale while designing specific public areas to host biodiversity [
105]. Land taxation can pay for new afforestation projects on public land in new development areas and re-balance the system’s ecosystem capability [
106]. Also, agroforestry projects can be optimal for connecting ecological corridors and existent biodiversity core sites.
Re-naturing solutions can also be considered in the existent built-up stock to reduce the anthropic pressure through ecological interventions, which include: de-sealing, building replacement, reduction in the urban footprint, and the adoption of nature-based solutions [
107].
The coupled integration of the above-mentioned recommendations (green mitigative solutions, compensation and afforestation, and re-naturing solutions) are all based on deep soil characterization, thus confirming the crucial role of soil knowledge in urban planning.
As for the impoverishment, shrinking, and consumption of rural areas described in
Section 3.1 (land-use change analysis), the situation is even more complicated, since the reduction in arable land (2592.38 ha), herbaceous vegetation associations (1498.71 ha) and mainly in pastures (1171.06 ha) is significantly affecting the entire primary production economy of the region, shifting from extensive local farm products to the massive plantation of monocultures such as corn production. Unfortunately, this trend has been evidenced by broad LUC analysis based on an extended temporal time-series (1990–2018) in Izmir’s metropolitan area [
58], where more than 34 thousand hectares of non-irrigated arable land became permanently irrigated land in less than thirty years.
The coupled action of the new city’s expansion at the expense of the most fertile land is generating tangible biophysical effects on ecosystems and social and ecological abandonment of more local, sustainable, extensive agriculture. These activities are based on mixed-use land (agroforestry), triggering people’s migration from inner farmland to the city center. In this case, alternative farming practices based on a stepwise introduction of more resilient, local, and less hydro-demanding farm production can be made possible only by long-term policies that envision the abandonment of intensive corn production and the shift to perennial cultivations. This requires an extensive vulnerability analysis of the water yield in each sub-catchment that supports the development of schemes and contracts with local farmers to favor sustainable, local productions.
Once more, the study of soil characteristics and the spatial evaluation of the ecosystem’s capacity can make a significant difference in planning. Indeed, land suitability analysis can immediately support the identification of different kinds of more sustainable rural products such as wheat, vegetables, and orchards while, in the long term, the introduction of perennial cultivations can represent optimal solutions also for the biodiversity of the region: olives, almonds, figs, grapes, pears, and chestnuts [
108].
The knowledge of the soil’s characteristics can play a significant role in planning for many reasons. For example, alluvial soils offer a perfect location for new buildings, but, at the same time, they are critical in terms of seismicity when considering geological and hydrogeological features. Important metropolitan cities such as Kathmandu (Nepal), Istanbul (Turkey), Mexico City (Mexico), and Los Angeles (USA) located on weak soils and fault lines in the world, have been exposed to loss of life and property due to repeated destructive earthquakes at specific periods. As viewed in the İzmir Earthquake Master Plan [
92], the average shear wave velocity V
s30 (usually the average seismic shear wave velocity in the upper 30 m) is the leading parameter for selecting the soil type (See
Figure 6). This parameter is also directly related to the expected amplification of the seismic shaking. Soil types A to D (vs. 30 more than 180 m/s), described in detail in [
109], range from rock or other rock-like formations to sandy clay soils. Type D and E soils are characterized by loose cohesionless soils and sharp contrast of hardness between a loose surface layer and a much more massive underlying foundation. In addition, potential liquefaction and wave amplification events during the earthquake threaten the existing structures as well as pose a risk factor in the expansion areas where new settlements are located.
4.2. Geologic Discussion
Earthquake surface rupture, another risk factor, is the appearance of the movement of the fault on the ground surface. Generally, these displacements occur as a result of destructive earthquakes along active faults. The probability of surface faulting and the amount of slip are related to earthquake magnitude, depth, fault geometry, and lithological features.
It is common for soil properties to have a destructive effect, such as earthquake wave amplification, and to rupture and damage directly under the structures [
52]. In this context, safety zones are defined at certain distances on both sides of the fault planes called earthquake prevention/avoidance zones. These zones at different distances have been proposed for different fault types worldwide. In addition, there are important studies on this subject in Turkey and the province of Izmir, but this protection method has not been applied much, except for a few certain regions [
110,
111].
While E–W trending fault segments are related to the major graben–horst systems of Western Anatolia [
59], approximately NE–SW trending fault segments are connected to the strike-slip-dominated İzmir-Balıkesir transfer zone [
51]. The critical fault segments are Gülbahçe, Seferihisar, Tuzla, and İzmir (
Figure 7). When we look around the city center, two critical tectonic zones affect the morphology of Izmir Bay. While the İzmir Fault bordered the Çatalkaya and Nif mountain elevations in the south, the Karşıyaka-Bornova faults bordered the volcanic Yamanlar Mountain in the north. Along these two fault segments, the ground forming the inner gulf of İzmir collapsed, and the waters of the Aegean Sea advanced to the internal parts [
51]. In these areas, the drainages are usually perpendicular to the faults. At the same time, these faults physically deformed the surrounding rocks. As a result of this deformation, soil formation in the character of permeable, unstable talus has occurred around the bay.
There were some studies conducted for the Izmir metropolitan area on earthquakes and soil liquefaction, and some risk mapping has been completed according to the results of this study [
52,
92].
Figure 10 shows a GIS application example of the region where one of the urbanization change areas determined in this study intersects with faults around Karşıyaka-Bornova. With this visualization, which is based on the principle of choosing a buffer zone on both sides of the fault in the range of 20 to 50 m (which may vary depending on the fault type: normal faults, strike-slip faults), we aim to limit destruction and casualties from potential surface ruptures while containing the urbanization [
112,
113]. Just as in this example, at many areas determined in the study area, fault segments, Quaternary alluvium, and areas with new urbanization spatially intersect, so there is a very risky situation in terms of surface rupturing.
4.3. Limits and Novelties
This study is mainly based on digital datasets (even though corrected/customized by direct field surveys), which were used to produce a large-scale assessment of the land take impacts on soil. This means that, in the case of hydraulic conductibility, the classification is based on a parametric association between soil types (materials) and conductivity. Even though this association is widely used for water run-off parametrization, it is well known that parametrization does not represent reality, especially in understanding the coupled integration between soil and vegetation, which can influence hydraulic characterization.
Also, the slope is another parameter that critically influences the infiltration capacity but also the urban expansion, which is limited in the mountain region. Indeed, the surface flow is dominant in high-sloping lands, while in the areas with flat topography, infiltration occurs at a higher rate compared to the sloping lands. Furthermore, not only the degree of the slope but also the kind of slope (concave, convex, etc.) affects infiltration, demonstrating that these parameters should be analyzed through a comprehensive characterization of the site. However, when we look at the findings of our field studies and the spread of geological formations in the region, permeability parameters are controlled by fracture-crack systems rather than soil hydraulic conductivity, since more rock units are located in the sloped regions.
We are aware that site-specific studies on environmental impact assessments of newly developed areas should be based on direct field observation and soil analysis developed at the local scale. Nevertheless, this study was not aimed at defining a new methodology for environmental impact assessment but rather demonstrating how much the past development trend has significantly impacted the most favorable soil in terms of water absorption and agricultural productivity, also resulting in deep contrast with the fault lines. In this context, GIS-based analysis revealed how the newly made site selections should be made by considering the soil characteristics. In this context, many important parameters not mentioned in this article based on engineering and architecture must be first investigated and then integrated into urban planning studies. Thus, it will be possible to create cities that are more resilient and less affected by future natural disasters. Furthermore, even though the economic benefit were not investigated here, the monetarization of the potential urbanization’s impact can play a vital role in easing the land transformations.
At the same time, the methodology adopted here can suggest essential changes in the urban planning decision-making phases. The dataset we employed is commonly diffused among different regions/countries; thus, we believe this approach can be easily replicated elsewhere. Specifically, during the preliminary phase of metropolitan planning, various development alternatives can support a suitability or conflict analysis while finding the less impacting and more resilient scenario for local communities. In light of climate change conditions, this assessment should be made before the final prefiguration of plans and projects to drastically reduce urban systemic vulnerability.
5. Conclusions
Are soil and geology characteristics considered in urban planning? This rhetorical question shaped the research article, where the competencies of geography, urban planning, topography, and soil science were simultaneously used to build a geodatabase for Izmir’s urban area.
First, we conducted a spatial analysis of the land-use change trends between 2012 and 2018 by employing a cross-tabulation matrix. We discovered the distribution and extension of the major land-use alterations in recent years. Then, to understand the impact of newly developed areas, new urbanizations were qualitatively analyzed by intersection with four ancillary datasets: the soil hydraulic conductibility, the land capability classification, the geological units, and fault lines.
The results confirmed a huge underestimation of the soil’s characteristics in urbanization processes: more than 53% of the total land take happened on A-type soil, thus causing the most problematic impact in terms of run-off control. Moreover, more than the 34% of new urbanizations were located on highly fertile soil, in terms of potential agricultural production (LCC). We discovered that new detached, peripheral structures were located on the city’s border, without continuity and contiguity with the adjacent built-up land.
Thirdly, we verified that new urban areas were built without considering the problem of soil liquefaction or on top of active fault lines.
All these limits were discussed while trying to address some policy recommendations, which included the effects of uncontrolled urban planning not only at the expense of the natural environment and the relative ecosystem services but also the long-term landscape degradation processes in the rural areas and the most problematic aspects of the exposure of citizens to earthquake effects.
Some specific recommendations were designed to address the problem of land take (green mitigative solutions, compensation or afforestation, and re-naturing areas), while other comprehensive solutions were discussed to reduce the intensification of corn cultivation in intensive farm productions.
However, in the current situation, it would be unrealistic to say that these measures can be immediately introduced to standard planning procedures in Turkey due to many procedural, bureaucratic and legal procedures that are not investigated by this study. Nevertheless, even though we know this practical limit, we demonstrated how geological characteristics can be employed to develop more rational and resilient urban planning.