3.3. 2017–2025 LULC Change
An analysis of land use changes between 2017 and 2025 in the city centre of Elazığ reveals a 3.30% increase in built-up areas, expanding from 12,029.66 ha to 19,349.66 ha. This growth is primarily attributed to the development of new residential zones following the earthquake, driven by the need for safer housing options. Additionally, within the framework of the Murat River Basin Rehabilitation Project, which was completed in 2023, afforestation efforts within the central district led to an increase in forested (tree-covered) areas by 0.15%, rising from 241.10 ha to 571.58 ha. In contrast, bare lands experienced a notable decrease of 1.40%, shrinking from 3,905.83 ha to 790.05 ha (Table 4).
Table 4
LULC change between 2017–2025
| | Area (2017) | Area (2025) | ∆ | 2017% | 2024% | ∆ % |
| Bare | 3905.83 ha | 790.05 ha | -3115.78 ha | 1.76 | 0.36 | -1.40 |
| Built | 12029.66 ha | 19349.66 ha | 7320.00 ha | 5.42 | 8.71 | 3.30 |
| Crops | 85875.36 ha | 85348.91 ha | -526.45 ha | 38.68 | 38.44 | -0.24 |
| Range | 99617.03 ha | 93126.70 ha | -6490.33 ha | 44.87 | 41.94 | -2.92 |
| Trees | 241.10 ha | 571.58 ha | 330.48 ha | 0.11 | 0.26 | 0.15 |
| Water | 20361.53 ha | 22843.61 ha | 2482.08 ha | 9.17 | 10.29 | 1.12 |
With regard to other land use classes, agricultural land (crops) decreased by 0.24%, from 85,875.36 ha to 85,348.91 ha, while rangelands declined by 2.92%, from 99,617.03 ha to 93,126.70 ha. In contrast, water bodies increased by 1.12%, rising from 20,361.53 ha to 22,843.61 ha (Fig. 5). This increase is likely associated with the construction of irrigation and livestock watering ponds within the central district, implemented as part of the Murat River Basin Rehabilitation Project (Fig. 5).
How land use has changed is detailed in Fig. 6.
3.4. 2040 LULC Scenario
In determining the 2040 LULC scenario, several spatial variables were considered, including slope, aspect, distance to fault lines, distance to rivers, distance to railways, distance to roads, and distance to reservoirs and water bodies. To assess the relationships among these variables, Pearson's correlation method was employed. The correlation matrix summarizing the strength and direction of these relationships is presented in Table 5.
Table 5
Pearson correlation results
| | Slope | Distance to Faults | Distance to Rivers | Distance to Railways | Distance to Roads | Distance to Waterbodies |
| Aspect | 0.0815 | 0.0128 | 0.0646 | -0.0050 | -0.1067 | 0.1351 |
| Slope | | -0.0596 | 0.0470 | -0.0121 | 0.1894 | 0.1139 |
| Distance to Faults | | | 0.2997 | 0.8824 | 0.2491 | 0.1660 |
| Distance to Rivers | | | | 0.3174 | -0.0002 | 0.2979 |
| Distance to Railways | | | | | 0.2882 | 0.1501 |
| Distance to Roads | | | | | | -0.1396 |
Pearson correlation analysis was conducted to examine the relationships between the morphometric and spatial variables evaluated in the study. The resulting symmetrical correlation matrix offers valuable insights into the direction and strength of linear associations among the variables.
The analysis revealed a very strong positive correlation between the variables "Distance to Faults" and "Distance to Rivers" (r = 0.882), suggesting that fault lines and river systems within the study area are either spatially distant from one another or follow a parallel distribution pattern. Furthermore, a perfect correlation was observed between "Distance to Rivers" and "Distance to Railways" (r = 1.000), indicating a significant spatial overlap between these two features.
Among the moderate positive correlations, notable relationships were identified between "Slope" and "Distance to Roads" (r = 0.189), and between "Distance to Faults" and "Distance to Railways" (r = 0.249). These findings imply that topographic slope and proximity to fault lines may indirectly influence the spatial alignment of transportation networks such as roads and railways.
Negative correlations were generally weak. For instance, the slight negative correlation between "Aspect" and "Distance to Roads" (r = -0.106) may reflect subtle location preferences of transportation infrastructure in relation to terrain orientation. However, due to the low magnitude of this correlation, its statistical significance and interpretive value remain limited (Fig. 7).
An analysis of the learning curve of the artificial neural network (ANN) model during the training process reveals a marked decrease in the loss values for both the training and validation datasets in the initial stages. The rapid decline in loss, particularly within the first 20 epochs, indicates that the model effectively adapts to the learning process. Following approximately the 50th epoch, both loss curves reach minimal values and stabilize, suggesting that the model successfully converges by the end of the training phase. Moreover, the close alignment between training and validation loss values demonstrates the model’s strong generalization ability and indicates a low risk of overfitting. These findings collectively suggest that the ANN model exhibits stable and robust performance across both training and validation datasets (Fig. 8).
The accuracy and Kappa coefficient values obtained from the classification analysis indicate that the model yields highly reliable and consistent results. According to the analysis, the overall accuracy was calculated as 97.52%, demonstrating that a substantial proportion of the classified samples were correctly identified, and the classification process was largely successful.
The evaluation of Kappa coefficients further underscores the robustness of the classification, independent of chance agreement. The \(\:{K}_{o}\) coefficient of 0.96227 corresponds to an “almost perfect agreement,” suggesting that the model performs significantly better than random classification. Similarly, the \(\:{K}_{h}\) value (0.96315) confirms that the classification is consistently accurate across all land use classes. Notably, the \(\:{K}_{l}\) coefficient of 0.99909 reflects an exceptionally high level of spatial accuracy, indicating that the model performs extremely well in predicting the correct locations of land use types.
Collectively, these performance metrics highlight the high contextual and spatial reliability of the classification method and affirm the scientific robustness of the results (Basheer et al., 2022; Cohen, 1960; Sajan et al., 2022; Zhang, 2016). The projected land use changes under the 2040 urban sprawl scenario for Elazığ are presented in Table 6.
Table 6
2040 Elazığ province LULC datas
| Land Use | Area | Percent |
| Bare | 781.21 | 0.35 |
| Built | 19336.31 | 8.71 |
| Crops | 85362.58 | 38.45 |
| Range | 93138.45 | 41.95 |
| Trees | 567.57 | 0.26 |
| Water | 22844.40 | 10.29 |
| Total | 222030.51 | 100 |
In the development of the 2040 land use scenario, a correlation analysis was conducted to evaluate the relationships among key spatial variables, including slope, aspect, distance to fault lines, distance to rivers, and distance to railways. The detailed correlation matrix, which illustrates the associations between these variables and their potential influence on land use change, is presented in Fig. 9.
An examination of land use data for the period 2017–2025 in the Elazığ Central District reveals that spatial transformation in the region has primarily occurred along the axes of urban sprawl, loss of natural areas, and the conversion of agricultural lands. The most prominent change during this period is a 60.89% increase in built-up areas. This substantial growth indicates a pattern of urban expansion from the city center toward peripheral zones and reflects zoning policies that appear to favor construction activities. Moreover, this trend aligns with existing literature, which highlights that post-disaster housing policies in developing countries often promote settlement expansion into rural areas (Araya & Cabral, 2010; Seto et al., 2012). Critics contend that urban sprawl results in inefficient land use, escalates infrastructure development costs, and fosters automobile dependency, all of which collectively undermine the sustainability and vitality of urban centres (Palmieri & Serre, 2019). Similar patterns of urban expansion have been observed in many developing cities across Turkey (Güneş, 2021; Hayrullahoğlu et al., 2021).
The expansion of settlement areas has occurred in inverse proportion to the decline in natural land cover classes. Notably, bare lands experienced a substantial reduction of 79.97%, indicating that urban development has predominantly taken place on previously unused or vacant land. However, the decrease of 6.52% in rangelands and 0.61% in croplands highlights that urban growth is not confined to idle spaces but is also encroaching upon productive rural areas. This trend raises concerns regarding its potential long-term impacts on food security, rural development, and biodiversity conservation (Seto et al., 2012).
On the other hand, tree covered areas exhibited the highest proportional increase, rising by 137.17%. This significant change may be attributed to afforestation efforts undertaken within the city, the implementation of landscaping projects, or improvements in the accuracy of classification methods. However, the sustainability of such increases warrants further evaluation based on criteria such as the species of trees planted, the spatial extent of the afforested areas, and the continuity and ecological functionality of these landscapes (Barredo & Demicheli, 2003).
Similarly, the 12.21% increase observed in water bodies suggests that the hydrological infrastructure in Elazığ such as dams, ponds, and flood control systems has undergone notable changes in parallel with urban development. Elazığ, particularly recognized for its proximity to the Keban Dam, exhibits a spatial utilization of water resources that differs from the national average in Turkey.
According to the 2040 scenario maps generated through Artificial Neural Network (ANN) supported MOLUSCE modeling, it is projected that construction pressure on natural areas will intensify if current development trends persist. The model demonstrates a high level of predictive reliability, with an overall accuracy of 97.52% and an \(\:{K}_{o}\) coefficient of 0.96, underscoring the robustness of its estimations (Basheer et al., 2022; Sajan et al., 2022).
The 2040 scenario further suggests that if existing planning policies remain unchanged, the spatial development of Elazığ will continue in a manner inconsistent with the principles of sustainability. In this context, it becomes evident that to manage future urbanization trends, the establishment of ecological corridors, strict protection of agricultural lands, and the imposition of limits on construction pressure are essential measures (Lambin & Meyfroidt, 2011; Palmieri & Serre, 2019).
In light of these findings, it can be concluded that urban sprawl in Elazığ Central District has occurred in a low-density yet spatially impactful manner, producing observable effects on ecological systems, rural landscapes, and socio-economic structures. Comparable spatial dynamics have been reported in other rapidly developing cities across Turkey (Ewing & Hamidi, 2015). This pattern of urbanization highlights the urgent need for integrated green infrastructure strategies, a balanced rural–urban interface, and the adoption of sustainable land use policies in future planning initiatives.