Monitoring of Changes in Land Use/Land Cover in Syria from 2010 to 2018 Using Multitemporal Landsat Imagery and GIS

Understanding the effects of socio-ecological shocks on land use/land cover (LULC) change is essential for developing land management strategies and for reducing adverse environmental pressures. Our study examines the impacts of the armed conflict in Syria, which began in mid-2011, and the related social and economic crisis on LULC between 2010 and 2018. We used remote sensing for change detection by applying a supervised maximum likelihood classification to Landsat images of the three target years 2010, 2014, and 2018. Based on the computed extent of our LULC classes and accuracy assessment, we calculated area-adjusted estimates and 95% confidence intervals. Our classification achieved an overall accuracy of 86.4%. Compared to 2010, we found an increase in spatial extent for bare areas (40,011 km2), forests (2576 km2), and urban and peri-urban areas (3560 km2), whereas rangelands (37,005 km2) and cultivated areas (9425 km2) decreased by 2018. It is not possible to determine whether the changes in LULC in Syria will be permanent or temporary. Natural conditions such as climate fluctuations had an impact on the uses of the natural environment and cultivated areas during the study period, especially in regions suffering from water stress. Although seasonal precipitation patterns and temperature affect LULC change, however, we could not identify a prevailing climate trend towards more drought-prone conditions. Our analysis focuses on (potential) direct and indirect implications of the Syrian conflict on LULC change, which most notably occurred between 2014 and 2018. Conflict-related main drivers were human activities and demographic changes, which are mainly attributable to large-scale population displacement, military operations, concomitant socio-economic status, and control of local resources. As the study provides quantitative and qualitative information on the dynamics of LULC changes in Syria, it may serve as a framework for further relevant conflict-related research and support planning, management practices, and sustainable development.


Introduction
In the face of a growing world population, changing demands of land-based products, environmental changes and feedbacks, and increasing global interconnectivity due to ongoing globalization processes, understanding the diverse and complex interactions of social-ecological land systems and linkages across different scales remains a major challenge. Changes in land use/land cover (LULC) can result in direct and indirect impacts on natural habitats, biodiversity, and terrestrial landscapes. Hence, the conversion processes of LULC are a major driver in global environmental change [1][2][3][4]. As land cover is defined as the combination of biotic and abiotic components on the Earth's surface, it includes water bodies, natural vegetation (grassland, forests), bare rock/soil, and other The current population estimate of Syria is roughly 13.3 million people, according to the international organizations and the Central Bureau of Statistics (CBS) in Syria [69,70]. Prior to the Syrian conflict in 2011, the population of Syria was more than 23.5 million people, of whom 47% lived in rural areas [71]. The conflict resulted in considerable devastation in Syria, which triggered a population loss of about 43%, which included displaced people, deaths, and missing persons [69,70]. The rate of natural increase of the population within the country was 2.5% for the period 2010-2018 [69,71]. It is now estimated that about half of the population lives in major cities [66,69]. LULC in Syria are classified into four major groups as follows [65,71]: (1) Cultivable land (6.2 million hectares, ha). The cultivated land area is estimated at 5.7 million ha or 94% of the cultivable land, which accounts for only 33% of the total area of Syria. About 62% of the cultivated area is located in the three northern governorates (Aleppo, Ar-Raqqa, and Al-Hasakeh); (2) uncultivable land (3.7 million ha, 20% of the total area of Syria) includes buildings, public utilities, bare land, rocky land, and rock and water bodies; (3) rangeland, steppe and desert land (8.2 million ha, 44% of the total area of Syria),

Development Stages of Influence and Control in Syria from the Outbreak of the Conflict in 2011 until 2018
As LULC changes detected in this study took place during the conflict years, they are therefore related in one way or another to conflict circumstances and areas of control of the conflicting parties. Thus, it is important to describe the development stages of influence and control across Syria from the outbreak of the conflict until 2018. Since the outbreak of the Syrian Revolution in 2011, leading to an open civil war and border disputes with external forces, the control of various regions has been subject to several changes as illustrated by the following chronology of events [72,73].  September 2018-A new agreement is reached between Russia and Turkey. The objective of this agreement is to freeze the combat operations and to maintain the ceasefire in these areas, the deescalation zones.  [72,73] at the end of 2010 (a), 2014 (b) and 2018 (c), and the share of land control between the warring parties adapted from [72] (d).

Figure 2.
Distribution of zones of influence in Syria adapted from [72,73] at the end of 2010 (a), 2014 (b) and 2018 (c), and the share of land control between the warring parties adapted from [72] (d).
• March 2011-The first protests against the regime quickly spread across the country and security forces respond with arrests and shootings.  September 2018-A new agreement is reached between Russia and Turkey. The objective of this agreement is to freeze the combat operations and to maintain the ceasefire in these areas, the de-escalation zones.

Data Source
To analyze changes in LULC between 2010 and 2018 in Syria, we used images of the Landsat-5 Thematic Mapper (TM) and Landsat-8 Operational Land Imager (OLI) from the Landsat archive [75] with a spatial resolution of 28.5 m and 30 m, respectively. These images contain radiometric information of three visible and three infrared bands with spectral characteristics suitable for monitoring changes in LULC [18,76,77]. The list of Landsat TM/OLI imagery we acquired is given in Table 1. In general, due to photosynthetic activity, the growing season is more convenient for monitoring and mapping vegetation cover than the dormancy period [78]. The main planting season in Syria is the rainy season from January to February [65]. In this season, most winter crops of cereals and legumes are grown all over Syria to avoid the risk of frost. We selected Landsat images of May and April to capture maximum vegetation cover, as winter crops reach their green peak during this reference period [65]. It is often difficult to obtain cloud-free scenes above the coastal mountain range due to high cloud cover in March and April. Due to the lack of suitable images, we had to resort to scenes from the years 2015 and 2017 in two cases.
Supplementary data can be integrated as external inputs for remotely sensed data in a classification process according to different concepts [79]. Despite a lack of field observations due to the conflict in Syria, there was sufficient auxiliary data to identify the characteristics of the natural environment and land uses in the study area. The valuable data included the Google-Earth archive, topographic maps of various scales ranging from 1/25,000 to 1/100,000 covering the whole of Syria, and an assessment study of the International Center for Agricultural Research in the Dry Areas (ICARDA) [80]. We used these data as ground control points (GCPs) and accordingly as reference points for the geometric correction of Landsat images. In addition, we utilized the data to train the classifier, to assess the accuracy of the classification and to better visually interpret the Landsat images by deepening the background knowledge on land use in the study area. Figure 3 illustrates the methodological processing steps we performed to analyze LULC dynamics changes during the study period. Although the Landsat imagery level 1G used is radiometrically and geometrically corrected by the USGS, difficulties may counter when applying the polygon-based mosaic classification due to different atmospheric conditions [65]. We accounted for this by implementing geometric corrections, atmospheric corrections, best bands composite selections, and final subsetting in the pre-processing. In addition, pre-processing procedures for coarse and medium spatial resolution data, such as Landsat images, have a great benefit in presenting the objects of these images and in increasing visual interpretability [81].

Pre-Processing of the Satellite Data
Land 2020, 9, 226 9 of 30 to ensure the correctness of the GCPs selected. We found the RMSE to fall below critical thresholds, assuring the quality level and excluding distortion. [81,85]. Primarily, we performed the atmospheric correction to eliminate radiation imbalances and atmospheric impacts caused by scattering, absorption, and reflection using the software ENVI 5.3 [86][87][88]. We further processed the mosaic datasets (images) in ArcGIS 10.4 to isolate spectral bands that are highly correlated to avoid data redundancy at the edges during image classification [89].

Classification Uncertainties and Controls
The medium spatial resolution of the Landsat imagery and the spatial overlaps among LULC classes with a small spatial extent can result in spectral overlap, i.e., mixed pixel problems, which in turn can lead to erroneous classifications.
One of the issues encountered is the classification of planted lands with trees, especially rainfed lands, in bare areas. Extensive cultivated lands are spread over light-colored soils. Images of these We rectified and coregistered the Landsat images using a second-order polynomial model. This procedure is the most commonly used for the geometric correction of Landsat images since it is simple and fast, does not change the digital number value of pixels nor does it require additional sensor geometry data [82][83][84].
Based on this model, we reconstructed raw Landsat images by using the nearest neighbor resampling algorithm for all bands. Furthermore, we calculated the root-mean-squared-error (RMSE) to ensure the correctness of the GCPs selected. We found the RMSE to fall below critical thresholds, assuring the quality level and excluding distortion. [81,85]. Primarily, we performed the atmospheric correction to eliminate radiation imbalances and atmospheric impacts caused by scattering, absorption, and reflection using the software ENVI 5.3 [86][87][88]. We further processed the mosaic datasets (images) in ArcGIS 10.4 to isolate spectral bands that are highly correlated to avoid data redundancy at the edges during image classification [89].

Classification Uncertainties and Controls
The medium spatial resolution of the Landsat imagery and the spatial overlaps among LULC classes with a small spatial extent can result in spectral overlap, i.e., mixed pixel problems, which in turn can lead to erroneous classifications.
One of the issues encountered is the classification of planted lands with trees, especially rainfed lands, in bare areas. Extensive cultivated lands are spread over light-colored soils. Images of these lands in Google Earth showed that each pixel of Landsat satellite data typically contains between 10 and 16 trees. Therefore, it may be difficult to classify such disperse vegetation features with satisfactory accuracy on Landsat images. The spatial resolution of Landsat imagery was sufficient to distinguish irrigated tree crops and bare lands. Despite the difference in growth between rainfed and irrigated winter crops in Syria, particularly wheat and barley, the areas of these crops are homogeneous and similar in relation to spatial and spectral discrimination [65]. Therefore, it was necessary to include a series of Landsat images to identify the maximum growth dates for these crops, which are from mid-April until the beginning of May. This makes it possible to classify areas planted with these crops as agricultural land, as well as to distinguish between cultivated areas and bare lands according to the seasonal growth periods. The bare areas in Syria show great diversity, reflecting the geological structure and the complex topographic situation. Errors may arise due to low spectral separation of LULC classes and consequent misclassifications, especially in heterogeneous areas with mixed pixels at a small scale. Similar uncertainties arise with spectral and spatial overlaps between urban areas (built-up areas) and agricultural areas. This mainly concerns lands that fell fallow during the conflict or regions where spectral and spatial mixing between military artificial surfaces, urban areas, and bare areas occurred. We applied mask function in ArcGIS to separate different LULC from each other. A spectral and spatial overlap may also occur between natural water bodies and agricultural lands in the river beds, especially along the river Euphrates. These agricultural lands may be flooded with water as a result of changes in water capacity in spring, both due to natural conditions and the amount of water being released from the operation of dams in Turkey [65]. In addition, the impacts of the conflict on the area of actually cultivated land and water resource management practices need to be considered. To minimize classification uncertainty, the size and amount of training data in those regions were extended to provide better constraints for the classifier. We accounted for the bandwidth difference between Landsat 5 TM and Landsat 8 OLI sensors by applying an atmospheric correction. However, a study by Roy et al. (2015) identified higher reflectance of OLI sensors and that although atmospheric correction reduced differences in the NIR and the SWIR bands, the magnitude of the mean sensor reflectance differences in the visible bands increased [90]. Thus, any remaining impact on classification and change detection cannot be ruled out completely.

Ground Truth Data, Image Classification, and Accuracy Assessment
Various reviews and comparisons of detecting LULC changes with remote sensing-based techniques have been conducted with little consensus on optimal method for this purpose [91][92][93][94][95], including unsupervised classification (ISODATA-method), supervised classification and multistage supervised approaches (maximum likelihood classifier, neural network classifier, support vector machines). The maximum likelihood classifier (MLC) is a supervised classification approach that has been successfully used in studies conducted in comparable semi-arid environments with similar remote sensing data [46,63,65,84]. Therefore, we selected the widely used MLC algorithm due to its robustness, reliability, and its application suitability regarding our study area. The MLC-classifier requires training samples based on individual or composite spectral bands and, when applied, assigns pixels based on the maximum likelihood to the corresponding class.
We trained the MLC classifier with representative data assigned to all LULC classes that we derived from the ground truth reference points, utilizing bands 3 (red), 4 (NIR) and 5 (SWIR-1) of Landsat 5 TM and bands 4 (red), 5 (NIR) and 6 (SWIR-1) of Landsat 8 OLI". These data included 420 field reference points (GPS-measurements) in five Syrian governorates: Aleppo, Lattakia, Dara, Ar-Raqqa, and Deir ez-Zor. The GCPs formed the core set of our training data. To ensure the accuracy and reliability of the classified LULC categories, additional training samples were also collected by visual interpretations of Landsat imagery for the three study years. Thereby, we brought in our background knowledge of past and current LULC in the study area to confirm each GCP land-use allocation, if possible, with corresponding or nearly corresponding high-resolution imagery in Google Earth Pro. Furthermore, we utilized the NDVI to derive ground-reference points, especially when detecting agricultural areas in the eastern regions at the edge of the AL-Badia, where most of the bare areas or fallow areas are spread. The study by Foody et al. (2006) indicated that the size and number of training samples and the method of sampling affect the results of the land-use classification [96]. According to studies that successfully dealt with the classification of LULC in semi-arid regions, we set the minimum extent of the training areas to 50 m * 50 m [65,89,97]. In some cases, we chose larger training sites, especially in homogeneous areas, such as cultivated areas, to reduce the impact of technical geometrical noise that could appear in GPS measurements [65]. We aimed to cover the largest possible spectral range to achieve the representability of each LULC class in our training samples. We extracted the samples in the center of the representative and homogeneous regions [65]. The number and class distribution of our training samples are presented in Table 2. Subsequently, we applied the MLC algorithm to each Landsat image containing the spectral bands in false colors composite (band 5-4-3 as RGB for Landsat 5 TM and band 6-5-4 as RGB for Landsat 8 OLI). For each time step, we classified the Landsat images into six classes of LULC (Table 3). Afterward, we calculated the area of each category based on the pixel count. Finally, we smoothed the results of the image classification with the "majority-filter" (3 × 3 pixels) in ArcGIS to reduce the "salt & pepper" effect [98]. Table 3. Definitions of land use/cover classes in Syria.

Urban & peri-urban areas
Built-up areas and settlements (cities, villages), non-built-up areas (open structures related to human activities, such as quarries and permanent and temporary military zones)

Cultivated areas
Areas cultivated with crops: horticultural crops (tree crops, field and vegetable crops grown on small plots receiving intensive inputs), tree crops (olive grove areas and Aleppo pistachios trees), non-irrigated and irrigated crops (mostly rainfed winter/spring field crops and irrigated spring field crops).

Forest and other wooded areas
Undifferentiated coniferous and broadleaf, deciduous and evergreen, forest areas, scrubland, and degraded forest, often interspersed with tree crops Rangelands Sparse grasslands and shrublands of rocky hills and arid lands comprising the predominant grasslands and shrub areas used for extensive grazing and natural vegetation of saline depressions, which includes higher parts of saline depressions that are less frequently flooded than salt lakes Bare areas with or without sparse cover Unvegetated land, exposed rocks and burnt out areas

Water bodies Freshwater lakes and reservoirs, rivers, salt lakes, periodically inundated salt flats
To evaluate errors that may appear in the classification approach, we created a confusion matrix and computed the accuracies of the change detection results (overall accuracy, user's accuracy, and producer's accuracy) [99,100]. We based the accuracy assessment on a total of 1012 validation points to evaluate the classification accuracy of years 2010, 2014, and 2018 separately. We drew a stratified random sample (slightly modified inclusion probabilities) outside the training areas to generate validation points for each LULC class. Thus, the number of selected points in each class approximately matches the relative size of the corresponding LULC category. According to a recommendation by FAO (2016), we set the minimum sample size per strata to 20 samples [101]. However, we limited the number of water body samples due to its unique and differentiable spectral signature. We evaluated the validation points with high-resolution Google Earth Pro imagery (2010, 2014, and 2018) that correspond to the date (±15 days) of the Landsat images. We identified the current LULC allocation, compared it to our classified maps of 2010, 2014, and 2018 and subsequently calculated the accuracy. However, the classification of LULC may contain uncertainties, which can result from misclassification, homogeneity, and autocorrelations within verification samples during the manual preparation induced by the interpreter [102][103][104]. Thus, based on the error matrix and results of accuracy assessment, we calculated error-adjusted area estimates for each LULC class as well as confidence intervals at the level of 95% according to the statistical approach recommended by Olofsson et al. (2013Olofsson et al. ( , 2014 [104,105]. We carried out the LULC change detection by comparing the classification results of 2010, 2014, and 2018 based on adjusted area estimates for each LULC category.

Classification Accuracy Assessment
The area-adjusted overall accuracy achieved in the final classification maps was 87.2% (2010), 85.2% (2014), and 86.7% (2018) ( Table 4). The results accuracy assessment indicates that the errors associated with the overall evaluation accuracy are likely to be attributed to the low spectral separation of bare areas. In addition, overlap among different land-use classes such as rangelands and agricultural areas can cause spectral mixing and thus mixed-pixel problems occur which leads to difficulties in separating these areas in the Landsat data and subsequently to errors in the classification results. In the following, we refer to area-adjusted accuracy statistics, which however show similar results to former accuracies, except differences in 2014.
Taking all target years and all LULC classes into account, we obtained on average higher user´s accuracies (UA) than producer´s accuracies (PA), which overall indicates higher commission errors than omission errors. In 2010, rangelands covered approximately 44% of the total country area and achieved high PA of 92% and UA of 86%. It is therefore an important category for monitoring land use during the years of study. It should be noted that many regions assigned to the bare areas class can also easily be classified as rangelands due to grazing during the transition seasons (spring and autumn). Throughout all target years, bare areas showed on average both a high PA (85.3%) and UA (93.3%). The water bodies class obtained the highest PA results of 97% on average (2010, 2014, 2018) due to its distinctive spectral signature. The lower UA of water bodies (88.9% on average) can be attributed to commission errors in 2014 and the small sample size that assigns a strong weight to single misclassifications. The cultivated areas class achieved relatively low UAs compared to other LULC classes, with the lowest UA (75%) in 2018. This low accuracy is due to classification errors, where bare lands were misclassified as cultivated areas. Conspicuously, these areas are mostly located in the west and north of the country and correspond quite well to areas where precipitation is >250 mm/year. It, therefore, can be assumed that they represent predominantly cultivated crops, mostly cultivated rainfed crops, although complementary irrigation cannot be completely excluded. However, on the classified Landsat images, it is not possible to distinguish fully between areas that are consistently cultivated from bare areas. Bare areas contain small areas that are intensively cultivated with crops and vegetables in short periods of the year. Urban areas overall exhibited decent UA (81.6%) but by far the lowest PA (62.9%) on average. The urban areas category also includes open structures for human activities, such as quarries and military sites. Due to the fairly limited resolution of Landsat images, only large non-built-up areas can be displayed on the classified image. During the period under study, urban areas were mostly misclassified as rangelands and bare areas. The roads in the study area did not appear on the classified images because of mixed pixel signatures. Agricultural areas (fruit-tree plantations) could also be classified as forest cover. This is reflected in the PA of the forest classification. UA for this category remained largely unchanged. Forest cover in Syria is considered a rare category compared to other land uses. Therefore, the area-adjusted estimate of this category based on a statistical approach will be more accurate and reliable [104]. The results of the Land 2020, 9, 226 13 of 30 accuracy assessment of the supervised classification indicate an overestimation of bare areas, urban areas, and forest, whereas cultivated areas and rangeland, and an underestimation of urban areas to some extent. Table 4. Confusion matrices containing the (adjusted and not-adjusted) accuracy assessment of the classified LULC categories in Syria using a (slightly modified) stratified sample for the years 2010, 2014, and 2018. Here, the producer´s accuracy is PA and the user´s accuracy is UA.

LULC Change Dynamics
The LULC classification results are presented in Figure 4 and the corresponding LULC data in Table 5. The results show that the trend of LULC change between 2010 and 2018 did not witness any change before and after calculating adjusted area estimates, except for the forest cover. Area-adjustment accounts for differences in class accuracy and thus results in adjusted area extents that we refer to Land 2020, 9, 226 14 of 30 in the following analysis, also taking into account the 95% confidence intervals. To ensure good readability, we hereafter refer to a thousand hectares (tha) as a unit to quantify the magnitude of change (1 km 2 = 100 ha = 0.1 tha). Change detection revealed remarkable differences in the size and the trends among the LULC classes observed in the study period. Cultivated areas showed a reduction from 2010 to 2014 and a growth from 2014 to 2018. The most prominent changes in LULC were the continuous decline in rangelands and the increase in bare areas during the study period.
Land 2020, 9, x FOR PEER REVIEW 15 of 32 The LULC classification results are presented in Figure 4 and the corresponding LULC data in Table 5. The results show that the trend of LULC change between 2010 and 2018 did not witness any change before and after calculating adjusted area estimates, except for the forest cover. Areaadjustment accounts for differences in class accuracy and thus results in adjusted area extents that we refer to in the following analysis, also taking into account the 95% confidence intervals. To ensure good readability, we hereafter refer to a thousand hectares (tha) as a unit to quantify the magnitude of change (1 km 2 = 100 ha = 0.1 tha). Change detection revealed remarkable differences in the size and the trends among the LULC classes observed in the study period. Cultivated areas showed a reduction from 2010 to 2014 and a growth from 2014 to 2018. The most prominent changes in LULC were the continuous decline in rangelands and the increase in bare areas during the study period.

LULC Change Drivers
Based on the consistency of the classification assessment and the findings of the change detection analysis, we found that the land-use system in Syria has undergone rapid spatial and temporal changes during the eight years of study. In the introduction of the research, we have mentioned "high-level" driving forces and the corresponding direct and indirect drivers that could be behind changes in LULC. Given the large size of the study area, the complexity of the political and military landscape, and the change of influence and control areas on the ground in the light of the ongoing war in Syria, it was not possible to analyze all these drivers and show their direct and indirect impact on land-use changes during the study period. In the following, based on our analysis, and by reporting and interpreting evidence from various other relevant studies, we discuss potential direct and indirect drivers that may have led to these changes in LULC during the eight years of study. Thereby, we primarily focus on possible implications of the conflict, but also the potential effects of prevailing climatic conditions and trends.

LULC Changes and Conflict Conditions
During the study period, cultivated land was converted to non-cultivated land-use types in many areas in Syria. Generally, cultivated areas declined by 943 tha per year between 2010 and 2018. The results indicate a decrease in the use of arable land by 21% between 2010 and 2014 ( Table 5). We hypothesize that population loss during the conflict, especially in governorates with large agricultural areas such as Ar-Raqqa, Al-Hasakeh, Deir-ez-Zor, and Aleppo, had tremendous impacts on agriculture between 2011 and 2014, leading to a remarkable decrease in agricultural land. Analysis of demographic data from the Central Bureau of Statistics as well as studies that dealt with the humanitarian and agricultural situation in Syria showed a loss of agricultural employment [70,71,106,107]. Correspondingly, they reported a strong increase in the number of unemployed people and indicated population involvement in either the regular army or groups of the armed opposition. Thereby, related key processes are population displacement, forced migrations, and deaths. The number of refugees and displaced persons inside Syria reached 6.1 million people, while 5.5 million people took refuge outside Syria [108,109]. In 2011, more than 80% of the population in rural areas worked in agriculture [110]. Due to mismanagement of state institutions in the agricultural sector and intensification of rural-urban migration, agricultural areas, and thus the livelihoods of the population suffered substantial losses [109]. The FAO study on the state of agricultural production in Syria in 2013 included the governorates of rural Damascus, Homs, and Dara, and highlighted additional reasons for the decrease in cultivated land. The study suggested the lack of fuel for irrigation, plowing service, harvesting, and transportation of crops, as well as huge conflict-related agricultural production costs to have major impacts [110]. Another FAO study revealed that the area planted with cereal crops in 2014 was very low compared to previous years. According to the study, the area planted with wheat in 2014 was estimated at 900 thousand hectares, compared to 1.5 million hectares before the crisis [111]. As a result, farmers throughout Syria were no longer able to make ends meet, forcing many of them to abandon agricultural activity in search of more viable livelihoods. Poor management of water resources in the light of the conflict may also have driven people out of the rural areas, and thus further contributed to the decrease in the extent of agricultural land [109]. Military operations of conflict parties have had a direct impact on the diminishing area of cultivated land in most parts of Syria [73,109]. Large areas of agricultural land were destroyed in the course of military operations -in some cases deliberately -as happened in northeastern Syria, especially in Al-Hasakeh Governorate [108,109]. Military operations also had a major impact on the decline of tree crops, since tree crops were destroyed in many Syrian governorates, especially in Hama, Idleb, Aleppo, and rural Damascus [109]. Official reports by the Central Bureau of Statistics in Syria indicate that more than 51.7 tha of tree crops were cleared between 2010 and 2017, representing 12% of the total area planted with tree crops in 2010 and 7% in 2017 [69, 71,106,112]. In addition, most of the front lines of military operations occurred in or around agricultural lands in many areas in Idleb, Ar-Raqqa, Hama, and Aleppo governorates. Furthermore, agricultural land in these regions was confiscated by parties to the conflict. Some opposition forces also imposed taxes on agricultural products to secure financing [109]. The interplay of all drivers resulted in the displacement of large numbers of rural families from, and consequently to the abandonment of cultivated land. As the frequency of conflicts increased in most areas of Syria, growing insecurity disrupted the trade in agricultural crops and inputs, which resulted in shortages of fertilizers, energy, and manpower in agricultural activity [109]. As a result, agricultural production revenues continued to decline remarkably, thereby contributing to the decline in the area of agricultural land. According to the World Bank report in 2018, the GDP of the agricultural sector in Syria decreased by 41% in 2014 compared to 2010 as a result of the conflict´s implications [113].
Our change detection showed a slight increase in the area of agricultural land from 2014 to 2018 (80 tha). Since the magnitude of change is lower than the 95% confidence intervals of area extent in 2014 and 2018, we assume a class stabilization rather than an opposite (increasing) trend. This may be explained by the fact that the intensity of the conflict has declined between 2014 and 2018 in some Syrian governorates, especially in de-escalation areas. The relative stability in one form or another in these areas helped to create favorable conditions for agricultural activity, as did the return of many rural dwellers. Our change detection supports this assumption, as we observed the strongest increase in agricultural lands in governorates with de-escalation zones. This finding is consistent with data of previous studies that indicated a decrease in the number of unemployed people and an increase in the number of people in de-escalation areas whose economy is mainly based on agricultural resources [70,71,[106][107][108]. Profitable trade of agricultural products may also be a factor that contributed to increasing/stabilizing cultivated areas. This mainly affects crops, such as wheat, cotton, and barley, that are traded between regions of the regime and those outside its control, e.g., ISIS-controlled areas in which most of the agricultural land is concentrated. This is consistent with our finding of increasing agricultural lands in areas controlled by the Syrian regime and those that are under the control of the various opposition forces. The increase in the area of agricultural land between 2014 and 2018 can also be attributed to changes in the control and management of agricultural land. Cultivated land in northeastern Syria has increased in the area controlled by ISIS. This region includes the governorates of Al-Hasakeh and Ar-Raqqa and large parts of Deir-ez-Zor Governorate and contains nearly one-third of the agricultural areas in Syria. It remained under ISIS control until the end of 2017 [114]. In some areas of these governorates, croplands expanded to previously uncultivated areas, e.g., rangelands, driven by the purpose of enhancing agricultural production, and to compensate for the loss of agricultural areas that ISIS had lost in favor of opposition forces backed by the international coalition. Moreover, ISIS demanded landowners to maintain land cultivation to control food production and secure sources of financing through taxes imposed on agricultural products [114,115].
The spatial extent of rangelands showed relative stability, ranging from 46% to 55% of the country's area between 2003 and 2010 [116,117]. Our classification of 2010 demonstrates that rangeland areas have not undergone notable changes in their spatial extent until then and accounted for approximately 44% of the country. Analysis of LULC changes indicates that rangelands decreased substantially (34%) between 2010 and 2014 ( Table 5). The LULC classification of 2018 reveals a continuing trend of rangeland cover loss, with large areas of rangeland being converted into bare areas. Rangelands decreased by 8% between 2014 and 2018 ( Figure 4, Table 5). We found these results to be consistent with the pathways of conflict-related events, particularly in northeastern Syria and in the Syrian Al-Badia. Distribution maps of control zones between 2013 and 2014 indicate extensive and intense military conflicts between Syrian regime forces and opposition forces in the AL-Badia (Figure 2). Opposition forces eventually took control of most of the AL-Badia in June 2013 [73,74]. ISIS then acquired control of large swaths of land in Syria and Iraq, including Al-Badia, by January 2014 [73,74,118,119]. A consequence of these conflicts was the displacement of the population in these areas which was accompanied by the severe deterioration of livestock breeding and rangelands development [111]. Camp construction, urban expansion, the extension of military structures, and overgrazing and the lack of capacity to manage rangelands during wartime can contribute to a decrease in the rangeland area [65,108,120]. With the intensification of the conflict after 2014 in specific regions and the expansion of military operations, conflicting parties purposely removed vast areas of rangelands to expose fronts, facilitate surveillance, protect transport routes or infrastructure and to ensure the arrival of military and oil supplies. These practices were mainly concentrated in the Syrian Badia and on the sides of the main and subsidiary transport routes between the cities of Aleppo, Hama, Ar-Raqqa, and Deir ez-Zor [119,121].
Our classification result of 2010 shows that forests have a very limited spatial extent and accounted for only 2% of the total study area. Forests in Syria are mainly located in the coastal provinces (Lattakia and Tartous) and the governorates of Hama, rural Damascus, and Idleb (81% of the total forest area). The results of spatial-temporal trend analysis show that the forest cover increased by 30% between 2010 and 2014 and by a further 14% between 2014 and 2018 (Table 5). An increase in forest cover over several years in some areas of Syria may be not surprising. Growth of forest cover is likely to occur in areas where forests remained somewhat isolated from the effects of conflict, such as coastal governorates and Hama Governorate, where 62% of forest cover is located. Our outcomes differ from the results of previous studies and reports that showed a decrease in the area of forests in some coastal areas and the northwest of Idleb governorate for reasons related to the ongoing conflict 2010 and 2018, such as heavy logging to produce fuelwood and charcoal, and forest fires a result of mutual bombardment between the conflict parties [63,122,123].
It seems that the best explanation for this increase in forest cover between 2010 and 2018 is misclassification, which is due to the spectral confusion between the degraded forest cover and tree crops planted on terraced hillsides. To verify this outcome, the results of the forest cover classification for the three years were compared with the corresponding Google Earth Pro imagery. We found through visual inspection that forest areas, especially at relatively low elevations, are confused with areas of tree crops (Olives and fruit-tree plantations). This is also illustrated by the additional classification of forest cover at a regional or local level, which includes Lattakia Governorate, on the Syrian coast, where 30% of the forests and about 40% of the fruit trees are concentrated [71]. We also found that when the area classified on the image increases from the local to the national level, the total area of the forest cover increases. This can be attributed to the size of the classified area and the fragmented nature of the forest cover [124]. This can also be related to the spatial accuracy of Landsat imagery and the classification of mixed pixels containing dense tree crops and small forest spots as forest cover. To validate this interpretation, additional classifications related to the size of the classified area using different pixel sizes are required.
LULC classifications for the study years revealed that bare areas continued to expand between 2010 and 2018. The 2014 LULC classification shows that bare land increased by 41% and that the most prominent change was the conversion from rangelands to bare areas with or without sparse cover (covering 43% of the total area of Syria) (Figure 4). Most of these areas are located in the three eastern governorates of the country (Ar-Raqqa, Deir-ez-Zor, Al-Hasakeh) and the governorates of Homs and rural Damascus. The increase in bare lands between 2010 and 2014 (41%) is the outcome of rangeland degradation and shrinking areas of agricultural land. As a result of massive displacement of population and changes in the practice of land management, large areas of rangelands and agricultural land changed to bare lands, especially in areas that were controlled by ISIS in eastern and central Syria [108]. The increase in bare areas can also be attributed to the widespread distribution of mines through military operations, where the World Health Organization has shown that more than 8 million Syrians are at risk of lethal mines and remnants of war that spread over large areas and constitute a major impediment to land cultivation and pasture development [125,126]. Unfortunately, there is no accurate data on this important LULC driver in the period of war.
Our analysis reveals that the spatial extent of bare areas continued to increase after 2014. The results indicate an increase of bare areas by 9% between 2014 and 2018, albeit lower rates than before (Table 5). In addition to the drivers mentioned above, other main drivers behind the continuation of this increase in rangelands may be changes in land use activities and the escalation of the conflict between the warring parties in large parts of Syria after 2014. According to a report by Reuters, many farmers in ISIS-controlled areas did not cultivate their land for the 2015 season [127]. This was due to low fertilizer, lack of fuel, and instability as a result of the conflict, as well as deterioration of irrigated agriculture [127]. Tracking the course of the conflict-related events and the distribution of the control map indicates that gradually, after 2017, stability was reestablished in those areas from which the ISIS organization withdrew in the north-east and central Syria, in areas where military escalation was reduced as well as in areas where combat fronts were calm under the de-escalation agreements. Consequently, illegal economic activities began to grow on a large scale, such as oil production, refining, extraction of raw materials, and construction materials, which negatively affected the development of rangelands. This potentially contributed to the conversion of rangelands to bare areas, particularly in the Ar-Raqqa and Deir ez-Zor governorates [128,129].
In general, LULC change detection indicates a major increase in the spatial extent of urban and peri-urban areas between 2010 and 2014. Compared to 2010, the urban areas doubled and expanded by 105% until 2014. (Table 5). However, when considering confidence intervals, urban area expansion is subject uncertainty and ranges from 66% (low estimate) and 144% (high estimate). One possible explanation is the increase in population. Aleppo governorate and the eastern governorates of Ar-Raqqa, Deir-ez-Zor, and Al-Hasakeh particularly were affected by large and uncontrolled urban growth between 2011 and 2012 [130,131]. Poor management, oversight of state institutions at the beginning of the conflict, availability of resources, and the need to build new homes could encourage urban growth within existing urban areas and at urban fringes. The same pattern in similar scenarios was also detected in Gbanie's study (2018) [28]. With the intensification of the conflict and due to the fierce battles, large parts of urban areas became uninhabitable. Therefore, the wake of destruction, fallen buildings, and the desire to leave the most affected urban areas reinforce the necessity for new living spaces for the Syrian population and thus increase urbanization processes in the rural outskirts of the cities or even at new sites. As the conflict worsened further during the first half of the year 2013, people fled to relatively more stable areas, especially near the border and urban areas increased rapidly. As a result of the prevailing situation, large numbers of refugees, as well as urban and peri-urban dwellers, practiced urban and peri-urban agriculture in new landscapes to earn a living. This could ultimately lead to a marked expansion in urban and peri-urban areas. This case is in line with the findings of other studies, which showed that during forced internal migration of populations urban and peri-urban agricultural activities occur on a large scale, ultimately leading to an observed increase in the spatial extent of urban areas [28,132,133]. Some studies also indicated that many camps for internally displaced persons (IDP) have been established on large areas of agricultural land, particularly in the northwest (Aleppo and Idelb governorates) and the northeastern part of Syria (Al-Hasakeh governorate) [108,120,134]. This is what we can conclude from analyzing Google Earth images of 2012, 2013, and 2014. The camps and urban areas have a similar spectral signature. The spatial extent of these urban areas might have changed over time according to the military control of the parties to the conflict and the prevailing military situation. During the conflict, closed military zones were established and the existing military infrastructure expanded in the Syrian Coast region (in the provinces of Lattakia and Tartous) and in eastern and southern Syria (in the provinces of Deir-ez-Zor and rural Damascus), especially by countries that supported the Syrian regime in its war against the opposition [135,136].
The spatial extent of urban and peri-urban areas decreased between 2014 and 2018 by 11%. However, the urban class is subject to uncertainty, as the area change related to confidence intervals ranges from −8% (low estimate) to 29% (high estimate). We assume that the potential decrease in urban areas might be due to the strong population decline and low income of the population. The flight of civilian population to new regions inside the country and massive displacement of the population outside the country (33.8% internally and externally 26.6%) [108,130], as well as the decline in economic activities due to the country's volatile social and political environment, led to the disappearance of urban and peri-urban agriculture practiced by the displaced and local low-income residents [70,71,106,108,130]. Furthermore, several reports and conflict-related news in Syria indicate that large parts of villages, towns, and cities were partially or completely destroyed [137,138]. Both the scenario for the decline in urban area growth and the scenario of destruction of parts of the existing urban areas were observed in studies that dealt with similar armed conflicts around the world [3,139]. In addition, massive destruction led to strong fragmentation within urban areas in Syria. This may have resulted in isolated urban spots being replaced by agricultural areas or the sprawling bare areas Land 2020, 9, 226 20 of 30 on the classified images. Such patterns are consistent with previous studies that examined land-use changes during conflicts [2,16,140].
Our analysis of the spatio-temporal trend of water bodies showed a reduction of 48 tha (36%) by 2014. One explanation may be the uncontrolled increase of irrigated land, especially in the first two years of the conflict, i.e., before the refugee crisis, due to weak control by state institutions in Syria. Irrigated agriculture depends mainly on surface water reservoirs (stored water) such as lakes and dams built on waterways [141]. Another reason for this decrease in water may be low water storage of dams due to changes in water resource management practices during the conflict, especially after the large internal displacement of the population. Moreover, the low water storage of dams due to changes in water management practices as a result of the conflict situation in the country may have supported the decline of water bodies, especially after the large internal displacement of the population. For example, studies that examined changes in water resource in Dara governorate in southern Syria during the conflict indicated that the decrease by 49% in the storage of reservoirs in the study area is largely due to changes in water management practices during the period of refugee displacement between 2013 and 2015 [62,142]. This uncommon decrease in water levels was detected on satellite images as well as in the statistical analysis of water changes in the Syrian governorates during the study period. This result is consistent with the findings of previous studies that dealt with land-use changes in other conflict zones around the world [ 3,16].
The 2018 classification showed an increase in water bodies by 76 tha (48%) compared to 2014. One possible explanation is a considerable decrease in the area of irrigated agricultural land and consequently, the associated decrease in the use of irrigation water. This interpretation is consistent with other studies [62,143,144]. The decrease in irrigated area coincided with the large internal and external displacement of the population [ 108,130]. There are other reasons for the increase in the area of water bodies which can be attributed to the major changes in water resource management practices that have grown with the escalation and widening of the conflict [ 62]. From 2013 onwards, opposition groups controlled most surface water reservoirs such as dams in Syria. In late 2013, ISIS controlled the Euphrates River dam near the city of Ar-Raqqa in northern Syria, where the country's largest water surface is located [ 72]. The opposition forces lacked proper information about these dams as well as technical expertise and the institutional work to manage them effectively in a coordinated manner. Moreover, there was a lack of energy to operate these facilities [145]. This led to dependence on alternative sources, the groundwater. At the end of January 2017, the United Nations warned that damage to the Euphrates dam could lead to widespread flooding of the Euphrates Basin caused by the high water level behind the dam [146][147][148]. In addition, the water level of the Euphrates River rose about 10 m at the end of January 2017, which led to the inundation of large areas of agricultural land on both sides of the river for a long time [147]. These measurements indicate an increased reliance on groundwater for agriculture and other uses compared to surface water whose spatial extent has increased after 2014. Concurrently, this change in water usage was accompanied by massive population displacements both internally and externally, which may have caused changes in water resources management practices.

LULC Changes and Climatic Conditions
Climate variability, i.e., the interplay of the yearly temporal and spatial distribution of precipitation and temperature fluctuations, is considered a major driver of observed LULC patterns. Thereby, droughts are the most influential phenomena in Syria. Overall, a gradually declining average rainfall pattern since 1940 and constantly increasing average temperatures with a drastic upward trend since the 1980s can be observed in Syria, concomitant with long-term warming in the eastern Mediterranean [149][150][151][152]. Challenges related to climate variability and changes related to water use during the conflict can also lead to alterations in land productivity and irrigation needs, the large-scale failure of agricultural crops [152]. Moreover, climatic conditions are linked to the civil war as well as ongoing violent conflicts and resulting in poor management practices and unsustainable Land 2020, 9, 226 21 of 30 agricultural and environmental policies [153][154][155][156][157]. This can amplify environmental implications. Commonly, climate change is referred to as a threat multiplier [158]. The interaction of increasing temperatures and decreasing precipitation amplifies water stress, leading to more frequent and severe dry periods. Kelley et al. (2015) indicate a strong impact of droughts in Syria on the agricultural systems, especially in water-stressed areas. Occurring droughts were documented between 2007 and 2010 and in 2014 [115,152].
Low precipitation is directly reflected in the availability and extent of water bodies and may explain the decline from 2010 to 2014. The study by Müller et al. (2016), which examined the impact of the refugee crisis on changes in land use and water resources in southern Syria and northern Jordan between 2011 and 2015 also indicated that the decrease in volume of surface water stored in reservoirs and dams corresponds to the regional drought that Syria experienced in 2014 [62]. However, the report by JRC (2019) indicates strong climate anomalies and dry conditions for 2017 and 2018 as well [159] which does not correspond to the increase of water bodies from 2014 to 2018 found in our analysis. Therefore, we assume that the change in water cover is not a consequence of natural climate variability. Nonetheless, dry periods may have contributed to the decrease in cultivated areas, forests, and rangelands and the contrary increase in bare lands observed from 2010 to 2014 (Table 5). As two-thirds of the cultivated land in Syria is rain-fed, there is a high vulnerability and dependency on climatic conditions and natural year-to-year variability. This decline in agricultural areas in relation to climatic conditions often occurs in north and northeastern Syria [149]. In addition, the spatial extent of the irrigated agricultural areas may be affected by droughts due to this type of cultivation being linked to the volume of surface water stored in reservoirs, dams, and groundwater-fed by rainfall. According to the study by Müller et al. (2016) the low rainfall period 2013-2015 was accompanied by a marked decrease in the area of irrigated agriculture [62].
From 2010 to 2018, bare areas increase mainly to the expense of rangelands. The study by Houmsi et al. (2017) showed that there was a climatic shift of areas from dry-subhumid to semi-arid in the north and northwest of Syria due to drought during the period between 1981 and 2010 [64]. This shift was likely enhanced between 2010 and 2018 due to dry conditions in 2014, 2017, and 2018 [115,159]. This conclusion also corresponds to the results of our study, which showed an increase in the area of bare lands during the study period, especially in the governorates of Al-Hasakeh, Ar-Raqqa, Aleppo, and Idleb, where most of the bare lands are located (Figure 4). Rangelands are located in the more arid and fragile regions that are not suitable for cultivation. These vulnerable areas rely on seasonal precipitation and small changes in rainfall patterns and temperature may entail degradation processes and trigger land conversion. Thus, the consequences of the occurring dry periods for rangelands are more striking due to less resilience regarding water retention. Overgrazing may be another reason for the increase in bare areas at the expense of rangelands. According to the study by Al-Fares (2012), overgrazing led to the degradation of rangelands in drier semi-arid plains and dry interior regions [65]. This type of rangeland degradation is often associated with droughts and periods of little rainfall. This is consistent with the results of our study, which showed a decrease in rangelands and an increase in bare areas in the governorates of Al-Hasakeh, Ar-Raqqa, Deir-Ez-Ezor, Homs, and rural Damascus where the Syrian Badia forms the largest part of these governorates. It is believed that the effect of overgrazing is limited to the first years of study period i.e., until 2014 and to a specific spatial scale because of the exacerbation of the conflict between 2014 and 2018, along with droughts, presumably had the greatest impact on diminishing area of rangelands [111].
Quantifying climate change impact is subject to uncertainty and a major challenge, as numerous direct and indirect interactions contribute to LULC changes that vary across the regions. Our examination has shown that there has been no significant trend towards less/more drought-prone conditions in the study period. Further analysis is needed, e.g., by incorporating spatially resolved temperature and precipitation data or by assessing useful aridity-indices. Generally, long-term trends of precipitation and temperature do not necessarily translate into short-term effects. However, it cannot be ruled out completely, whether a short-term effect occurred and contributed substantially to LULC Land 2020, 9, 226 22 of 30 change. Climatic conditions have to be taken into account, however, we conclude that climate impacts regarding the fluctuation of precipitation and temperature were not the main driver for LULC in the study period.

Conclusions
We used multi-temporal Landsat in this study to monitor and analyze land use/land cover (LULC) change in Syria from 2010 to 2018. This study not only examines the spatial and temporal changes in LULC in Syria but also attempts to identify the causes of these changes by linking it to relevant published literature. Our results of spatial, temporal, and statistical analyses of land-use changes indicate that the degree of LULC change from 2010 to 2014 was greater than the change from 2014 to 2018 for all LULC classes except water bodies. We find that this change entailed severe conversions in the natural environment and agricultural land during the study period. The results also indicate that the spatial extent of cultivated areas and rangelands declined during the study period while the other LULC classes showed an increase. At this stage, it is not possible to determine whether the specific changes will be permanent or temporary.
Given the large size of the study area, the complexity of the political and military landscape in Syria, population movements and economic consequences, and the impact of droughts linked to climate change, it is difficult to identify single drivers of change for specific LULC types and to analyze their direct and indirect impact on LULC changes during the study period. By discussing our findings on the causes of change in LULC dynamics, we concluded that population loss and social unrest due to population displacement and long-term or short-term population pressure are the main potential causes of these changes.
We also concluded by tracking land control pathways during the study years that the spatial extent of LULC classes, such as agricultural lands and rangeland, is substantially affected by the prevailing conflict situation in north-eastern and southern Syria (governorates of Al-Hassakeh, Ar-Raqqa, Deir-ez-zor, and Dara) and in the Al-Badia (Homs and rural Damascus governorates), whereas other regions located in western Syria (Lattakia, Tartous and Hama governorates) appear to be less impacted by the political situation and the state of conflict, especially after 2014. It has become apparent to us that this can be attributed to differences in local dominant political power and the changing balance of power between the warring parties as well as resource management practices at local and national levels.
We found evidence, substantiated by published literature, that the effects of droughts on LULC during the study period mainly occurred in water bodies, cultivated areas, and rangelands, especially in water-stressed areas. Consequently, those changes naturally led to an alteration of other LULC classes, such as bare areas. By linking our results to literature published we could not identify a significant trend towards less/more drought-prone conditions in the study period and we, therefore, we do not consider climate variability and fluctuating climatic conditions as a major driver of LULC in Syria.
Monitoring the current dynamics of LULC changes in Syria and identifying the reasons behind them shows the continued change that will probably extend in the future. This most probably will continuously cause both rapid and continued increase in erosion and degradation of agricultural soils which in turn has serious economic, social and environmental consequences and impacts on land systems and its management, especially at the local level. The future trajectory mainly depends on the political and military scene in Syria and long-term climatic conditions associated with climate change.
The proposed procedure for detecting LULC changes in Syria is straightforward to apply and allows us to monitor future developments also by using remote sensing data from other seasons. To generate a more holistic view, temporally corresponding proximate drivers of LULC change, demographic and socio-economic may be implemented in a spatial model in the future. This study shows that remote sensing technology is an effective tool for monitoring temporal and spatial changes in LULC in Syria at the country-scale. Our study provides quantitative and qualitative information on the dynamics of LULC changes and aims to establish a basis for understanding the drivers of these