Spatiotemporal Impacts of Urban Land Use/Land Cover Changes on Land Surface Temperature: A Comparative Study of Damascus and Aleppo (Syria)

: Monitoring the impact of changes in land use/land cover (LULC) and land surface temperature (LST) is of great importance in environmental and urban studies. In this context, this study aimed to analyze the dynamics of LULC and its impact on the spatiotemporal variation of the LST in the two largest urban cities in Syria, Damascus, and Aleppo. To achieve this, LULC changes, normalized difference vegetation index (NDVI), and LST were calculated from multi-temporal Landsat data for the period 2010 to 2018. The study revealed signiﬁcant changes in LULC, which were represented by a decrease in agricultural land and green areas and an increase in bare areas in both cities. In addition, built-up areas decreased in Aleppo and increased in Damascus during the study period. The temporal and spatial variation of the LST and its distribution pattern was closely related to the effect of changes in LULC as well as to land use conditions in each city. This effect was greater in Aleppo than in Damascus, where Aleppo recorded a higher increase in the mean LST, by about 2 ◦ C, than in Damascus, where it was associated with greater degradation and loss of vegetation cover. In general, there was an increasing trend in the minimum and maximum LST as well as an increasing trend in the mean LST in both cities. The negative linear relationship between LST and NDVI conﬁrms that vegetation cover can help reduce LST in both cities. This study can draw the attention of relevant departments to pay more attention to mitigating the negative impact of LULC changes in order to limit the increase in LST.


Introduction
Land use/land cover (LULC) changes in urban areas are a critical factor affecting the physical characteristics of the Earth's surface due to the unique characteristics that each LULC category possesses with respect to radiation and absorption energy [1][2][3]. Therefore, the occurrence of changes in LULC will lead to changes in the exchange of energy between the ground surface and the atmosphere, which, in turn, will lead to changes in the radiative and aerodynamic balance [1,[4][5][6]. These changes are enough to significantly alter the microclimate in urban areas [1,6,7]. This will, therefore, have consequences for local and regional weather variables [6]. This can eventually lead to sharp transformations in urban landscapes [6,8].
Studies have shown that one of the most well-known forms of modification of landscape involving LULC changes is the microclimate modification phenomenon in urban areas, which is represented by an increase in atmospheric and land surface temperature (LST) [5][6][7]9,10]. Since LST can be considered an effective measure for predicting radiation budgets for heat balance, it can be an important indicator in understanding the interactions that occur between the environment and humans in urban environments [11,12]; thus, this reflects the health status of the ecosystem [6]. In addition, LST can be an important factor controlling the physical and biological processes of land systems [11,13].
LST is obtained from direct measurements taken on the ground. Thanks to the ability of remote sensing techniques to capture the spatial patterns of the emitted thermal radiation, data of the thermal sensors of the satellites such as Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), operational land imager (OLI)/thermal infrared (TIR), ASTER, and MODIS are also used to calculate LST [14][15][16][17][18][19]. The data of these sensors, despite their low temporal accuracy compared to the data of field measurements of temperature, has wide spatial coverage in different time scales [6,16].
Studies on LST analyses have shown that LST that can be obtained from remote sensing data is the function of different land uses [6,11,20]; therefore, the LST can provide data with a high spatial resolution and continuous spatial distribution of an urban area as a large city or region [11]. This allows for adequate visualization of spatial relationships between temperature patterns and patterns of land use in urban areas, including infrastructure characteristics within these areas [5,16,[21][22][23][24].
Several studies have been conducted to correlate LST derived from remote sensing data and normalized difference vegetation index (NDVI) as a predictor of LST [1,6,10,19,[25][26][27][28][29] to identify surface urban heat island (UHI) distribution patterns. In view of the negative impacts of LULC changes as a result of the increase in urbanization and the loss of vegetation cover in urban areas, several studies have been conducted on the dynamics of LULC and its effects on LST based on remote sensing data [1,6,10,16,17,24,28,[30][31][32][33][34]. The results of previous studies have also shown that there is a strong relationship between LST and LULC patterns [1,6,16,28,30,33]. As a result, relationships between LULC changes and LST derived from remote sensing data can be studied. Thus, qualitative studies can be useful in land use planning in urban areas. In addition, the study of the relationship between LULC and LST may be of great value in urban climate studies.
Most of the studies based on remote sensing technology in studying the pattern, density, and distribution of UHI through LST and LULC changes confuse the LST, atmospheric UHI, and the surface UHI. The UHI is detected by measuring the temperature of the air for a daily cycle using data from meteorological stations [6,35,36], which cannot be obtained from remote sensing data. Surface UHI can also be detected using remote sensing thermal infrared data [6,37]. Therefore, in this research, in order to avoid uncertainty, the term surface temperature (LST) and its relationship to changes in LULC will be used.
Despite the fact that Damascus and Aleppo possessed strong centripetal forces at the regional level and were growing rapidly in the Middle East before the conflict began in mid-2011, they are located in different geographical environments. As a result, it is assumed that both cities witnessed a major change in land use and changes in the intensity of human activities during the devastating war. In addition, extreme hot weather has occurred in Syria with increasing frequency since that time. Given the critical importance of LST and LULC change studies in determining the local climate and health as well as in advocating urban planning policies in urban areas, the subject of this research becomes very important to understanding how LULC changes in urban areas can alter the surrounding thermal environment. Therefore, the main aim of this research was to understand the changes of LULC with trends of LST retrieved from the thermal infrared band of Landsat TM and OLI in Damascus and Aleppo during the 2010-2018 conflict period as well as analyze the effects of LULC changes on LST. This involves: (1) determining the extent, direction, and patterns of LULC changes; (2) calculating the LST values and analyzing its patterns; (3) exploring the spatial relationships between LST and different LULC classes and determining the role of LULC dynamics in LST temperature changes.

Study Area
Damascus and Aleppo are located in the interior of Syria (Figure 1). The city of Damascus is located in southwestern Syria at the geographic coordinates of 33 • 30 36 N and 36 • 17 28 E and covers an area of 116.75 km 2 with an average altitude of 680 m above mean sea level. Due to the existence of the Anti-Lebanon Mountains, which form the border between Syria and Lebanon, the Damascus region is sometimes subjected to drought. These mountains, which have high peaks of over 3048 m, block the passage of rain-producing weather systems from the Mediterranean Sea (rain shadow effect) [38][39][40]. Due to the concentration of intensive agriculture as irrigated agriculture in Al-Ghouta, which surrounds Damascus, these cultivated areas can affect the city's microclimate as a thermal modifier and a cooling factor within the city. To the east of Damascus lies the desert (steppe). Damascus city experiences the cold desert climate according to the Köppen-Geiger system [39]. The winters are cold and rainy, while the summers are dry and hot with a lower moisture content. The average annual rainfall was about 131 mm for the period 1985 to 2018, falling between October and May [38,40]. The average annual temperature was 17.5 • C for the period 1985 to 2018 [38,40]. Aleppo is located on a plateau with a height of 380 m in northwestern Syria at geographical coordinates of 36 • 12 04" N and 37 • 09 40" E, which is approximately 120 km from the Mediterranean Sea. Aleppo occupies an area of more than 144 km 2 . Aleppo is surrounded by agricultural land to the north and west, while it is close to the borders of the Syrian desert to the east. In the city of Aleppo, there is the cool steppe climate according to classification by Köppen [39]. The mountain range along the Mediterranean coast largely prevents the arrival of Mediterranean influences from the west on the climate of Aleppo. The average annual precipitation was 329.4 mm for the period 1985 to 2018, 80% of this precipitation being concentrated between October and March. The average annual temperature was 17.3 • C for the period 1985 to 2018 [41].
Damascus is the administrative capital of Syria and is, therefore, the center of a large number of administrative buildings of state institutions and residential areas. Aleppo was the largest urban area in Syria before the conflict, which included a large number of major industrial and commercial institutions and offices of major companies and was known as the industrial and economic capital of Syria. Both cities have witnessed and continue to experience a different population growth during the current conflict [42,43]. The population of Damascus in 2010 was 1,749,000, as compared to 2,375,000 people in Aleppo [1]. Aleppo was the largest city in the country and one of the fastest growing cities in the Middle East until the beginning of the conflict. However, as a result of the siege of the city of Aleppo and the long battle that took place in this city, there was tremendous destruction in many parts of the city, and hundreds of thousands of its inhabitants have been displaced. Aleppo has become the second largest city in Syria after the capital, Damascus, where the population of the city of Aleppo was estimated at about 1,850,000 people in 2017 [42,44,45]. Thus, Damascus became a densely populated metropolitan area with an estimated population of five million. The rate of growth in Damascus is currently the highest in the country, due to the war conditions that imposed new demographic and urban growth. However, the population of Damascus has declined in recent years as a result of the ongoing conflict, with a population of approximately 2,048,000 people in 2018 [1]. Population density in Damascus was about 17,541 people per km 2 in 2018, while the population density of Aleppo was lower, which amounted to 12,847 people per km 2 in the same year [1].

Data Acquisition
Landsat-5 Thematic Mapper (TM) and Landsat-8 operational land imager (OLI) images taken in late spring, from April, were used to monitor the changes in LULC, calculate the LSTs and NDVI, and evaluate the response of different LULC on LST in both cities for the years 2010, 2014, and 2018. These images, processed at level-1, have been obtained free of charge from Landsat archives available on the official United States Geological Survey (USGS) website [46]. These Landsat images have radiometric information of three visible and three infrared bands with spectral characteristics, which allows LULC monitoring and mapping as well as retrieval of LST [47][48][49][50].
Since the objective is to compare the changes in LULC and LSTs in both cities, images were chosen for each year of the study period from the same year and represent the same season. In addition, the scenes covered by these images were completely free of any coverage of clouds and captured exactly at one moment in time. Therefore, Landsat images are coherent enough to be interpreted as a single set when monitoring and analyzing changes in LULC and LSTs for the target years. Table 1 shows details of Landsat TM/OLI images used in the study.
Although there were no field observations to collect ground truth data due to the ongoing conflict in Syria, there were sufficient additional data to identify the characteristics of the natural environment and land uses in the two study areas. These data included high-resolution Google Earth Pro archive data and topographic maps of Aleppo and Damascus in scales ranging from 1:25,000 to 1:50,000 which were provided by the General Organization for Survey of Syria. Some aerial photos of small areas in both cities were also obtained from the military survey and are available online (photos of drones) [51]. Ground control points (GCPs) were derived from this additional data, which were used as reference points for the geometric correction of imagery. The additional data were also used to train the classifier, to validate the accuracy of the classification, and to better visually interpret the Landsat images by deepening the background knowledge on land use, which, in turn, provides a better understanding of LULC dynamics and the reasons behind these dynamics and changes in LSTs in both cities.

Preprocessing of the Landsat TM and OLI Imageries
The Landsat images of the level-1 are geometrically and radiometrically corrected by the USGS; nonetheless, there may be difficulties in applying the pixel-based classification of LULC, because the digital number of spectral images are often affected by various atmospheric conditions, such as atmospheric scattering, absorption, and reflection, as well as numerous factors such as sensor calibration, illumination geometry, sun angle, and visible near-infrared wavelengths [1,6,49,52,53]. Therefore, radiometric calibration is a preprocessing step to overcome radiation error caused by atmospheric scattering. This process aims to accurately determine the radiation brightness values at the sensor input and to further convert these values to the actual ground-reflected radiation [3,54,55] in order to derive good estimates of incoming radiation from the ground surface [16]. Furthermore, an atmospheric correction of all used bands is essential in techniques that use band ratios (such as NDVI and NIR/red ratio) [55]. Moreover, preprocessing procedures for the medium spatial resolution data of Landsat images have a great benefit in correcting different surface reflection differences, displaying the objects of these images, and increasing the possibility of visual interpretation during the classification process [1,56].
As a result, all the acquired images were subjected to preprocessing steps, which included geometric correction, atmospheric correction, and radiometric processing. The acquired Landsat images were preprocessed individually using ENVI 5.3 software of Harris Corporation in Melbourne, FL, USA, where these correction procedures have corresponding operation modules in the toolbox. To match the pixel sizes of all bands, the nearest-neighbor algorithm was used to resample the TM and OLI thermal bands to a 30 m pixel size. This step of preprocessing was performed using ENVI 5.3 [57]. The study area was then clipped from Landsat images for each study year using ArcGIS 10.5 [58]. The color combinations of Landsat imagery bands were also chosen, which take into account the contrast and distinction between different LULC, especially vegetation formations (band 5-4-3 as RGB for Landsat 5 TM and band 6-5-4 as RGB for Landsat 8 OLI) [50].

Image Classification and Accuracy Assessment
Classification of Landsat images of urban areas is a complex process due to the spectral heterogeneity shown by different LULC [1,59]. In previous studies, which addressed the detection of change in LULC patterns in semi-arid areas as comparable environments using similar data, the maximum likelihood classifier algorithm (MLC) achieved the best results [49,[60][61][62][63][64]. In addition, this classification method is known for its strong theoretical basis, its high classification accuracy, and its ability to accommodate changing spectral signatures of different LULC with appropriate selection of classifier training samples [1,49,64]. Therefore, the MLC algorithm was applied using ArcGIS 10.5 for Landsat image classification purpose for the two study areas. Classifier training samples were obtained as predefined spectral signatures for all classes of LULC from Google Earth archive data, topographic maps, and high-resolution aerial imagery (drone images). In order to achieve greater accuracy and reliability in the classification of LULC in the two study areas, additional training samples were also collected using visual interpretation from Landsat images and NDVI images, especially when detecting agricultural lands and green areas in the two study areas. Training samples were selected independently for each Landsat image in accordance with the size and location of LULC categories within the two study areas [1,49,65]. This ensures, as far as possible, to obtain different spectral signatures that represent all categories of LULC. Table A1 in Appendix A shows the number of training samples represented for the classified LULC categories. The MLC algorithm was then applied to each Landsat image containing the spectral bands in false colors composite. In order to detect LULC changes and better describe their patterns based on knowledge of the characteristics of land use, the density of development, and the social and economic levels of both cities, all images were classified into three categories of LULC ( Table 2). The area of each LULC category was calculated using the number of pixels. Finally, the LULC classification results were spatially smoothed using the "majority-filter" (3 × 3 pixels), to reduce the "salt&pepper" effect [66]. Table 2. Definitions of land use/land cover (LULC) categories in study area.

LULC Definitions
Urban or built-up area Built-up areas and settlements; non-built-up areas (pen structures related to human activities, such as quarries).

Agricultural land and green areas
Areas cultivated with crops: horticultural crops (tree crops, field and vegetable crops, which are grown on small plots receiving intensive inputs). Green and open areas in the urban areas (gardens parks, and green belts). Forest areas and other wooded areas.

Sparse grasslands and bare areas
Sparse grasslands and shrublands of rocky hills and arid lands. Unvegetated land, exposed rocks, bare areas and burnt out areas.
In order to assess the classification accuracy, a confusion matrix was applied to classified images [67,68]. For this purpose, a stratified random sampling approach was taken by allocating random sample points outside the training areas for each LULC category based on the minimum sample size for each stratum, which is 20 validation points [49,69,70]. The selected number of these random validation points in each category corresponds to the relative size of the corresponding classified LULC category [69], with a total of 218 and 230 points for the cities of Damascus and Aleppo, respectively (Table A2). The validation points for each LULC category in both cities were evaluated based on the high-resolution Google Earth Pro imagery that match the Landsat imagery dates (±15 days); then, the accuracy was calculated. The supervised classification of LULC can contain uncertainties, which can result from misclassification, homogeneity, and autocorrelations within the validation samples selected during manual preparation induced by the interpreter [71][72][73][74]. Therefore, based on the error matrix and accuracy assessment for the LULC classification, the error-adjusted area estimates for all LULC classes as well as the 95% confidence in-tervals were calculated according to the statistical approach recommended by Olofsson et al. [73,74]. Then, LULC change detection was performed based on adjusted area estimates in accordance with the three levels of development (urban planning) in both cities (Table 3, Figure 2). Table 3. Development Levels according to urban planning in Damascus and Aleppo in accordance with the traffic rings.

Development Levels Definitions
City center (inner-city core) within the inner traffic ring Which is characterized by very dense residential areas and commercial buildings.

Semi-urban/urban fringe between the inner and outer traffic rings
Which is characterized by dense residential areas, commercial and industrial structure and some small farming areas.

Rural outer fringe beyond the outer traffic ring
Which is characterized by low to moderate density residential neighborhoods, farming lands, and natural environment such as sparse grasslands and bare areas, except several major industrial areas such as the Aleppo industrial area in the northeast (industrial area Sheikh Najjar) and in the south of Aleppo (industrial area Ar-Ramouseh) and the industrial areas in the south of Damascus.

Derivation of NDVI
NDVI is calculated from the reflectance values of the visible (RED) and near-infrared bands (NIR) using Equation (1).
where NIR = band5 and NIR = band4 are spectral reflections in OLI, and NIR = band4 and NIR = band3 are spectral reflections in TM. The value of NDVI varies in the range −1 to 1. Negative values of NDVI correspond to non-vegetative cover or uncultivated surfaces. Positive values refer to cultivated soils or areas that contain some form of vegetation cover [25,62]. NDVI was applied to preprocessed Landsat images used in this study to determine proportion of vegetation coverage (P v ) in order to calculate the emission (ε) as input for the algorithm required for calculation of the LST and to analyze the relationship between vegetation cover and LST results. The images of NDVI were also useful for separating vegetation from other LULC categories during the classification process.

Retrieval of LST
Several studies dealing with changes in LSTs in arid and semi-arid environments have recommended the use of band 6 of Landsat 5 images because it does not show radiometric saturation in the hottest regions [6,75]. Band 10 from Landsat 8 is also recommended for the LST calculation because it has the best adjustment with respect to LST computation, compared to band 11, which shows an overestimation [6,76]. Therefore, the single band method was used to retrieve LST from thermal band 6 from TM images and the thermal band 10 from the OLI images. This method, which requires fewer parameters, is effective in obtaining the surface temperature from these thermal bands [16,36].
Given that the two study regions have limited spatial ranges and Landsat images of these regions were acquired in a cloudless sky, the degree of impact of atmospheric conditions will be consistent, so the relative temperature of the land surface and its spatial distribution across the study areas will also not be affected by these conditions. The following three steps show derivation of LST from Landsat TM and OLI sensors:

•
Conversion of pixel value from digital number (DN) value to spectral radiance value at the sensor entrance of the thermal band (L λ ) using Equation (2) [77].
where L λ is the spectral radiation value in W/(m 2 srµm); L min and L max are the radiance minimum and radiance maximum for the thermal band, which was acquired from a Landsat image metadata file; Q max is the maximum quantized calibrated pixel value, which was acquired from a Landsat image metadata file; and Q dn is the data value of image pixels, expressed as a digital number (for sensors TM = 255 and for Landsat 8 = 65,535).
where T b is brightness temperature at satellite (sensor) expressed in Kelvin; K 1 and K 2 are thermal conversion constants for thermal bands from a metadata file; and ln is a natural algorithm. K 1 and K 2 for band 10 OLI are 774.88 and 1321.08, respectively. K 1 and K 2 for band 6 TM are 607.76 and 1260.56, respectively.
where LST is the land surface temperature in Kelvin; λ is the wavelength of emission radiance (11.5 µm in band 6 for Landsat 5 TM and 10.8 µm in band 10 for Landsat 8 OLI); Due to the different ground surface materials, the corresponding pixels contain a variety of materials with different emissivity values, especially at resolutions of 100 m and 120 m [3]. This makes the calculation of emission very complicated due to its different values. To calculate the emission values for the study areas, the ground surface materials on Landsat images were divided into two types: the land of construction and the natural ground surface. The land surface emissivity of these two types was calculated using the NDVI thresholds method according to Equations (5)-(9), respectively [78,[80][81][82].
where ε surface is emissivity of the natural ground surface; ε construction is the built land emissivity; P v is the proportion of vegetation coverage; d ε is the effect of the geometrical distribution of the heterogeneous and rough surfaces as well as the internal reflections of these surfaces.
NDVI soil value can be approximately equal to the minimum value of NDVI, and NDVI veg value can be approximately equal to the maximum value of NDVI.
where F is a shape factor and is equal to 55 as average value, assuming that the Earth's surface has different geometrical distributions.

Correlation Analysis
To show a quantitative relationship between LST changes and LULC as urban parameters and to illustrate their spatial differences in both cities, a multiple linear regression analysis was used to calculate the correlation matrix between the different study variables, LULC, NDVI, and LST. To achieve this, the values of LST ranges calculated for each year and the size of the area for LULC extracted from the classification occupied by these ranges were used. To obtain the surface temperature ranges as an input parameter to calculate this correlation, the retrieved LSTs were divided into seven intervals according to the rules shown in Table 4, and then the corresponding LSTs were identified for each interval. Table 4. Rules for dividing the intervals and ranges of land surface temperature (LST) into seven categories for both study areas, adapted from [3].

Interval of LST Range of LST ( • C)
Low temperature Lst < Tavg + 2Std Sub-low temperature The study also employed quantitative analysis to determine independent changes in LST for LULCs in both cities based on LST derivation for different types of land use that was derived by visual interpretation on classified Landsat images: residential areas, non-residential areas, agricultural land, industrial areas, green areas, bare land, roads, airports, and military zones. This was achieved by converting pixels of LST into point data. Parametric values of these points were retrieved from LST maps derived from Landsat images. Sample points were then selected for comparison by overlaying the LST values of these points with the selected LULCs.

LULC Change Detection and Dynamics Analysis
Table A2 in Appendix A shows LULC classification accuracies in both cities based on the confusion matrix between 2010-2018. The area-adjusted accuracy of the final overall classification of LULC in Damascus varied between 87.1% and 89.4%, while this accuracy was between 86.1% and 90.0% in Aleppo. It is established that the standard overall accuracy for LULC classification ranges from 75% [83] to 90% [84]. Figure 3 presents the LULC pattern of Damascus and Aleppo. Table 5 also shows the results of LULC classification in both cities. The results of the LULC change revealed significant differences in the size and trajectory of land use during the study period, particularly in the city of Aleppo, which witnessed a siege and battles between 2012 and 2016 during the conflict. Analysis of LULC maps between 2010 and 2108 showed that cultivated areas and green areas in Damascus exhibited growth and decline during of conflict between 2010-2014 and 2014-2018, respectively, especially outside the outer traffic ring, while these areas in Aleppo recorded a continuous decline between 2010 and 2118 within and outside the outer traffic ring. The most prominent changes in LULC in Damascus were the continuous growth in urban or built-up areas and sparse grasslands with bare areas between 2010 and 2018 outside the outer traffic ring, while the city of Aleppo witnessed, for all the study years, a continuous decline in urban area as well as a continuous increase in sparse grasslands and bare areas inside and outside the outer traffic ring (Figure 3 and Table 5). The results indicate that the two study areas have undergone significant and different changes in LULC during the study period. This will lead to a corresponding change in the surface characteristics of these LULC, which implies subsequent changes in land surface albedo and surface energy balance.      From the above analysis, it is clear that there are qualitative differences in the intensity and distribution of LST in the cities of Damascus and Aleppo during the years of study. The results also indicate that the spatial pattern of LST in both cites has changed from a scatter pattern in 2010 to a more homogeneous pattern of LST in 2014. However, the scatter pattern reappeared to some extent in both cities in 2018, but with a larger proportion in Damascus.

Relationship between LST and LULC
The results shown in Tables 5 and 7 were used to validate if changes in LST are relevant to LULC in both cities. A 2010 linear regression analysis indicates a weak correlation coefficient between LST ranges and the corresponding LULC in both cities (34% for Damascus, 10% for Aleppo) ( Figure 6). The linear regression model between LST and LULC in 2014 shows a strong correlation coefficient in both cities, so that this coefficient explains 72% of the spatial relationship between these two variables in Damascus and 68% in Aleppo. For 2018, a weak correlation of 31% between LST and land cover is found for Damascus, compared with a strong correlation of 66% in Aleppo.  The results of the quantitative analysis to determine independent changes in LST for LULCs showed that LST values increased for almost all LULC types in both cities between 2010 and 2018, where this coincided with significant changes in LULC (Table 8).
For residential areas, under the influence of changes in LULC, an increase in LST was recorded in Damascus by 5 • C, while this increase was 4 • C in Aleppo. However, due to the nature of the various surface materials of these areas (different canopy structures: concrete, stone, etc.), their LST were not high compared to other urban areas. The calculated results also indicate that high urban areas in both the cities have shown higher LST from 2010 to 2018 than low urban areas. These areas in the city of Damascus recorded LST of 2 to 3 • C higher than in the city of Aleppo for all the study years. The increase ranged from 2 to 4 • C for non-residential urban areas in the city of Damascus and from 2 to 5 • C in the city of Aleppo. Non-residential urban areas associated with military-industrial activity also recorded significantly higher LST than other areas, particularly outside the outer traffic ring in both cities. The increase in LST in military sites such as airports and military barracks, which largely overlap with industrial zones, was about 5 • C in Aleppo and 4 • C in Damascus between 2010 and 2018. In contrast, the LST change was negative with about −8 • C in industrial areas in Aleppo for the same period. The industrial zones located behind the outer traffic ring in the city of Damascus recorded an increase in LST by 2 • C during the study period, which tended to increase in 2018. The results showed that bare areas in both cities recorded a slight increase in temperature between 2010 and 2018 (less than 1 • C). Although agricultural land is directly related to low LST, it has recorded high LST in the cities of Damascus and Aleppo. This increase was 2 • C in Damascus, while this was greater in Aleppo at 3 • C. The results of the LST analysis for the years 2014 and 2018 indicate that the agricultural lands that recorded a high LST are concentrated in the southwest of the cities of Damascus and Aleppo. In general, the green areas of both cities showed a relatively lower LST compared to other land uses, but this was in larger areas of Damascus. These areas in Damascus also showed higher temperatures, by approximately 2 • C, than Aleppo. The spatial distribution of LST within the green areas varied between 2010 and 2018, and this distribution was strongly influenced by the structure of land use and density of vegetation cover, especially inside and beyond the outer traffic ring in both cities.

Relationship between LST and NDVI
Linear regression analysis indicates a negative relationship between the values of NDVI and the LST values in the two cities over the entire study period (Figure 7). The regression analysis for 2010 showed a strong correlation coefficient between NDVI and LST. For the city of Damascus, 75% of the relationship is explained by NDVI, whereas in Aleppo, the relationship is stronger, at 83%. For the year 2014, the correlation coefficient in Damascus decreased significantly to 55%. By contrast, the coefficient of determination in Aleppo has increased slightly (2%), so that this coefficient explains 84% of the spatial relationship between these two variables. The low correlation between LST and NDVI in Damascus for the period 2010-2014 shows that there are other possible causes, which are responsible for LST rise. It is possible that these causes are related to uncontrolled industrial and human activities that have emerged in the city as well as local weather conditions. For 2018, the relationship between the two variables LST and NDVI increased in the case of Damascus compared with 2014, which recorded a strong correlation coefficient (70%). In the case of Aleppo, this relationship declined significantly, where the correlation became normal in 2018 (63%) after it was strong in 2014. Spatial analysis of the relationship between LST and NDVI led to the conclusion that the NDVI is largely responsible for the spatial variation patterns of LST in both the cities. This is also confirmed by the analysis of the relationship between spatial distribution of LST and LULC ( Figure 6).
The analysis of the relationship between LST and NDVI also showed that the surface temperature value decreases generally with the increasing value of NDVI in both cities during the study period. It has been concluded from the above that the relationship between these two variables is linear (Figure 7).

Discussion
In the years preceding the war in Syria, the cities of Damascus and Aleppo experienced a high rate of urbanization. This has been reflected in the continuous construction of buildings and structures as well as the significant expansion of infrastructure. The war broke out in mid-2011, where the consequences of this devastating war represented by population displacement, military operations, and socio-economic crisis had a significant impact on the change in LULC, which appeared in a controlled and uncontrolled increase in urban areas and a decrease in agricultural lands in Damascus, while in Aleppo, this was represented by a decrease in built-up areas, a significant loss of green areas and agricultural lands, and an increase in bare areas. These changes were especially large and rapid in the city of Aleppo between 2014 and 2018 compared to the city of Damascus, which can be attributed to the intense military operations, in addition to other possible causes such as industrial and human activities related to the conflict that appeared in the city of Aleppo, especially after 2016.
As mentioned in the introduction to the study, the degradation or loss of vegetation cover and the presence of unplanned urban activities in urban environments lead to a change in the characteristics of the land surface, which leads to changes in the thermal energy balance and modification of moisture distribution, which, in turn, leads to changes in the local climate, which is represented by a local increase in LST and changes in natural environment services. Therefore, this comparative study focused on monitoring the dynamics of LULC changes and their impact on LST in the cities of Damascus and Aleppo, the two largest cities in Syria.
Quantitative analysis to determine independent changes in LST for LULCs showed that LST values increased for almost all LULC types in both cities during the study period, where this coincided with LULC changes (Table 8). For the two study areas, statistical analysis of changes in LULC between 2010 and 2018 indicates that agricultural land and green areas decreased by 22 and 23% in Damascus and Aleppo, respectively. Additionally, urban or built-up areas increased by 4% between 2010 and 2018 in Damascus. In contrast, urban or built-up areas decreased by 10% over the same period in Aleppo. There was also a corresponding change in the surface properties of these areas, implying subsequent changes in the albedo of the land surface and a significant influence on local temperature modification in urban spaces of the two cities.
The results indicate that the city of Aleppo recorded a higher LST than the city of Damascus during the study period. The results of the study showed that the mean LST had an increasing trend in both cities throughout the study period. Quantitative analysis of LST also indicates that the difference in LST between urban areas within the outer traffic ring and the urban fringe with vegetation cover (agricultural land and green areas) was, on average, 2.8 • C for Damascus and 3.1 • C for Aleppo over the study period. The mean LST differences between urban and suburban areas in both cities increased to approximately 2 • C in Damascus and 3 • C in Aleppo during the study period. This is a characteristic of large cities that have experienced changes in LULC and changes in the functioning of urban ecosystem services [16,24,28,85,86]. The results also indicated that LST values differed according to land use in both cities and the increase in heterogeneity in land cover features, especially in Damascus (Figures 4 and 5).
An analysis of the relationship between LULC classes and LSTs indicates that changes in LST coincided significantly with changes in vegetation cover (agricultural land and green areas) and sparse grasslands with bare lands in both cities. Thus, the relationships between LULC and LST tend to change with changes in vegetation cover in Damascus and Aleppo. The values of the significant correlation coefficient between LST and LULC, and, in particular, in Aleppo, clearly reflect that both bare areas and non-residential areas, which increase at the expense of vegetation cover, cause a rise in LST. The results also showed that the changes in LSTs were closely related to the changes that occurred in non-residential areas in both cities as a result of the change in human activities and the accompanying changes in LULC in these areas during the study period. This can be seen by comparing LULC classes with LST imagery and values ( Figure 5 and Table 8). These changes in LULC were especially large and rapid in the city of Aleppo between 2014 and 2018 compared to the city of Damascus, which can be attributed to the intense military operations, in addition to other possible causes such as industrial and human activities related to the conflict that appeared in the city of Aleppo, especially after 2016. All this had a significant impact in modifying the conditions of the near-surface atmosphere as well as the surface temperature in these areas, which was reflected in the increasing trend in LST during the study period.
The results indicate that, after 2010 (as a base year), some non-residential urban areas, especially outside the outer traffic ring, recorded a significantly higher LST than other LULC due to the change in their use and their relation to industrial and military-industrial activity during the years of conflict in both cities. In contrast, other non-residential urban areas recorded lower LST for 2014 and 2018 compared to 2010. Industrial areas that have expanded at the expense of other LULC show a higher LST than other urban areas due to the presence of many artificial heat sources and the nature of building materials with high reflectivity, especially concrete and sheet metal roofs, unlike buildings in homogeneous urban areas (concrete, stone, etc.) that trap more heat during the day [27]. Industrial areas in Damascus city located behind the outer traffic ring recorded an increase in LST of 2 • C during the study period, which tended to increase in 2018 after the fighting subsided in the outskirts of the city, which led to an increase in industrial activity. On the other hand, for the same period, the industrial areas in Aleppo showed a negative change in the surface temperature of the earth, by approximately −8 • C. This can be attributed to the massive decline in industrial activity due to the conditions of the war, as Aleppo was considered the first industrial and economic center in Syria before the conflict (Table 8).
Military sites, such as airports and military barracks, which largely overlap with industrial areas, also recorded an increase in surface temperature during the study period, by approximately 5 • C in Aleppo and 4 • C in Damascus.
The results of the analysis of the spatial pattern of LST distribution in relation to the spatial extent of LULC patterns indicated that the high LST are distributed on the urban fringe and not in the center in both cities during the study period. This is associated with a significant loss of green area and scattered and dense agricultural land, an increase in bare areas, and fragmentation of peri-urban areas. This increase in surface temperature in the edges compared to the city center may appear in large cities whose urban fringes have sparse vegetation cover as well as an increase in unplanned urban activities. Here, local climate change processes may occur at different spatial scales as a result of environmental stress in these peripheral areas [6]. This urban periphery, which is considered exposed, is characterized by a higher reflection of thermal radiation, in contrast to the buildings inside and outside the inner traffic ring, which are largely homogeneous in building materials, and which can store a large part of solar radiation energy during the day [27]. Additionally, this observed variation in the LST distribution between the center and the urban periphery in both cities can be related to what is known as the phenomenon of the urban cold pool or urban cold island (UCI) [87]. This may explain the time of capture of the Landsat images used in this study.
The study also observed that the values of NDVI in the city of Damascus were higher than the values recorded in the city of Aleppo during the study period. However, the spatial extent of the vegetation cover (agricultural land and green areas) in the city of Damascus was much lower than in the city of Aleppo. There was a clear impact of the decrease in vegetation cover and the uncontrolled human activities associated with the conflict in both cities, especially in Aleppo, on the NDVI. The change in land use indicates that the spatial extent of green areas witnessed a significant decrease from 2010 to 2018 inside and outside the inner traffic ring in Aleppo (8 and 14%, respectively), compared to a slight decrease in Damascus for the same period, where this decrease amounted to 2.5 and 4%, respectively. While the area of agricultural land outside the outer traffic ring decreased by 18% in Damascus and 13% in Aleppo. The decline in vegetation cover in the city of Aleppo is primarily due to the acute water crisis, especially in the period between 2012 and 2016, which led to a significant lack of interest in green areas, in addition to the ongoing military operations and poor institutional management during the conflict and the displacement of the population to the countryside and other cities which experienced relatively greater stability, as well as outside the country [88][89][90]. It appears that the main driving factor for the decline in agricultural land in Damascus was the battles, which were concentrated, until 2017, on the outskirts of Damascus, especially in the areas of Eastern Ghouta and Western Ghouta adjacent to agricultural areas [91]. This decline in agricultural land in Damascus can also be attributed to the uncontrolled urban growth witnessed by the city, especially outside the outer traffic ring, as well as to the marginal expansion in urban areas, which was represented by the development of military and security infrastructure at the expense of agricultural land in the north and northwest of the city, especially between 2014 and 2018 [92].
The results of the NDVI obtained in this study for both study areas indicate that the values of this index tend to decrease in general during the study period. This decrease in NDVI is accompanied by a decrease in vegetation cover in both study areas, with an increase in the alteration in the natural environment to other LULC. The results of this study identified the areas of vegetation inside and outside the inner traffic ring in both study areas as areas with higher NDVI and lower LST. With regard to the low LST in the center of the two cities, inside and outside the inner traffic ring, they are distributed in areas with dense and natural vegetation cover and gardens, as well as in fields of dense cultivation, which show high NDVI values. Here, the vegetation, by process of evapo-transpiration, helps reduce LST [6,27]. This high relationship between the high NDVI values and the accompanying low LST within the inner and outer traffic ring of both cities is reinforced by the homogeneous distribution of vegetation cover. Regarding the vegetation cover outside the outer traffic ring of the urban area for both cities, there was no good consistency between the distribution of high NDVI values and low LST values. The urban fringes show different coverage and densities of vegetation cover and clear gradations of degradation due to changes in use over short periods of time. The increase in the values of the NDVI, which is accompanied by an increase in the LST in some areas in the urban area outside the outer traffic ring in both cities, especially in the city of Aleppo, can be attributed to local natural conditions such as high relative humidity in the air, where it is likely that vegetation cover did not affect the LST in these areas. A similar study indicated this scenario [93]. Part of the reason for this could also be related to the application of the NDVI, where it is difficult to extract information on specific types of land uses in urban environments using this index, especially bare lands and urban or built-up areas due to their relative similarities in terms of spectral signature [1,94]. Among the other possible reasons for the low relationship between NDVI and LST in the urban area outside the outer traffic ring in both cities is the concentration of industrial and commercial activities as well as the intense military activity after 2014 due to the war, most notably the city of Aleppo. Due to these intensive activities as well as the use of very poor fuels on a large scale as a result of the aggravation of the energy crisis, heat energy and pollutants, such as soot and black smoke, are released. These pollutants can absorb solar radiation as well as prevent the radiation emitted from the land surface from being released into the atmosphere. All this can contribute to an increase in the LST [27].
The results of the study showed a negative linear relationship between LST and NDVI in the cities of Aleppo and Damascus during the study period, and this means that the decrease in vegetation cover in both cities contributed greatly to the increase in the LST. These findings are consistent with the results reported in other similar studies that were conducted based on the analysis of Landsat data [14,15,[95][96][97][98]. In general, the detection of LULC changes and the values of the NDVI and its trend indicate that the vegetation cover within the two cities is not vigorous enough to maintain the energy balance of solar radiation and, therefore, does not contribute enough to maintain the ecosystem functions in the urban spaces of both cities.
There are no previous studies on LST changes in relation to LULC changes before and during the conflict in urban environments in Syria in order to compare the results. Additionally, there is no data on normalized LST to ensure that the LST values retrieved from the remote sensing data are comparable. However, it is possible that there are slight or significant differences in LST when compared with the results of similar studies on Syria. As mentioned in Section 3.1, the Landsat images used in this study were taken on a specific day and time. Since the surface characteristics and solar radiation balance change from time to time and from season to season, this may lead to differences in the spatial distribution of LST in both study areas. This could result in differences in the values and, thus, the power of the specific exponential relationships between LULC and LST as well as between the NDVI and LST, but these relationships are assumed to remain in the same trends. It can also be inferred that both cities are prone to increased LST for the reasons addressed in the analyses presented in this discussion.

Conclusions
This study is considered one of the first studies dealing with the spatial determination of land surface temperatures (LST) and their relationship with land use/land cover (LULC) changes using multi-temporal remote sensing data for Syrian cities. This study, as a comparative study on the two largest cities in Syria, Damascus, and Aleppo, which takes into account the conflict conditions and their potential effects on LST as well as LULC, can draw attention with its results. In this study, qualitative and quantitative analyses were conducted to study the relationship between land cover changes and land surface temperature between 2010 (as a base year before the war) and 2018.
In light of the advantages and disadvantages that were encountered in this research, several conclusions were reached: (1) the largest change in LULC was the decline in agricultural lands and green areas and the increase in bare areas in both the cities during the study period. Urban areas (build-up areas) also decreased in Aleppo, while these areas increased in Damascus. The main force for changing LULC in Damascus and Aleppo was the war and the consequent impact of politics, economy, and population distribution.
(2) The temporal analysis of the LST showed that the increase and decrease in the LST was largely related to the land use conditions of each city during the years of study. (3) Land cover patterns and changes have contributed to the modification of the local climate and also affected the value and distribution pattern of LST as well as its spatial variation between the center and the urban fringes, primarily through the degradation and loss of vegetation cover. The effect of this was greater in Aleppo than in Damascus. This was represented by the increasing trend in the minimum and maximum LST as well as the increasing trend of the mean LST. (4) In general, the spatial variation of LST was significant between urban areas in the center and suburban areas of both cities in 2010 (before the conflict). This variation was greater in Aleppo, but during the conflict, although the city of Aleppo has a larger vegetation cover, nearly similar LST were recorded in both cities. This means that similar conflict conditions in both cities had a close impact on LULC changes and thus on LST. (5) The relationships between LULC and LST tend to change with changes in vegetation (agricultural land and green areas) cover in Damascus and Aleppo. (6) The vegetation cover effectively contributes to mitigating the LST, especially in the center of the two cities, and this is confirmed by the negative linear relationship between the high normalized difference vegetation index (NDVI) values and the low accompanying LST. (7) The vegetation cover in its current extent cannot maintain the ecosystem services in both cities, as it is still subject to continuous degradation and loss as well as change in its density within short periods of time. This is also reinforced by changes in the characteristics of the surface as well as the local natural conditions for other LULC types.
The approach used in this study was important for detecting changes in LULC, NDVI, and LST. The results also showed that the Landsat images were ideal for monitoring LULC changes and LST and analyzing their spatial relationships. However, it may be problematic to choose images with similar conditions in terms of the environmental and climatic conditions of the urban area as well as the time of acquisition. This can be verified by further studies on the two study areas. An analysis of this relationship between LULC and LST can be carried out using high-resolution satellite images for different seasons of the year in future studies, taking into account local climatic conditions, such as wind speed, relative humidity, topography, as well as functional activities in both cities in greater detail. Therefore, more future investigations will also be required. In light of the current political, social and economic conditions related to the war in Syria, it is difficult to acquire data from field measurements. Therefore, the approach applied in this study, based on remote sensing technology, allows to monitor changes in the LST of cities quickly, accurately, and at low cost.
The persistence of these changes in the sensitive LULC of the ecosystem could make both cities environmentally more vulnerable, especially Damascus, which has more grasslands and bare lands that could, in the future, be converted to built-up areas. Ultimately, this study draws the attention of responsible departments to pay more attention to mitigate the negative impact of changes in LULC, at least to limit the reduction in vegetation cover and increase the number of green spaces in densely populated urban areas that would help reduce LST and contribute to the sustainability of both cities.  Table A2. Confusion matrices containing the (adjusted and not-adjusted) accuracy assessment of the LULC categories in Damascus and Aleppo 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.