War and Deforestation: Using Remote Sensing and Machine Learning to Identify the War-Induced Deforestation in Syria 2010–2019

: Armed con ﬂ icts and other types of violence are key drivers of human-induced landscape change. Since March 2011, Syria has been embroiled in a prolonged and devastating armed con ﬂ ict causing immense human su ﬀ ering and extensive destruction. As a result, over ﬁ ve million people have been forced to seek refuge outside the country’s borders, while more than six million have been internally displaced. This study focuses on examining the impact of this con ﬂ ict on forest cover by identifying the drivers of forest change. To assess this change, Landsat and PALSAR imagery were used to di ﬀ erentiate between forested and non-forested areas. Spectral information was syn-thetized using the Tasseled Cap transformation and the time series data was simpli ﬁ ed and despiked using the LandTrendr algorithm. Our results show that between 2010 and 2019 there was a substantial decrease of 19.3% in forest cover, predominantly concentrated in the northwestern region of Syria. This decline was induced by the armed con ﬂ ict, with several key drivers contributing to the decline, such as illegal logging activities conducted by both locals and refugees living in nearby forest areas. Drivers such as proximity to refugee camps, roads, and se tt lements played an important role in producing this change by facilitating access to forests. In addition, the occurrence of explosive events such as bombings and shelling near forests also contributed to this decline by causing forest ﬁ res. To mitigate further deforestation and reduce dependence on forests for fuel, it is crucial for local governments in the post-con ﬂ ict period to o ﬀ er sustainable alternatives for heating and cooking to both the local populations and refugees. Additionally, governments are recommended to enforce strict laws and regulations to protect forests and combat illegal logging activities. These measures are essential for preserving and restoring forests, promoting environmental sustainability, and ensuring the well-being of both displaced populations and local communities.


Introduction
Landscape dynamics are triggered by both ecological and anthropologic factors; as a result, they can be thought of as the combined effect of natural events and human activities on the environment. Changes in vegetation cover can generally be attributed to human-caused climate change, as well as land-use change drivers such as fires and overgrazing, both of which might have a substantial impact on the world's ecosystems [1][2][3][4][5].
Armed conflicts and warfare are among the most intensive and far-reaching human activities that can contribute to environmental damage and are key factors of human-caused landscape change [6,7]. Several studies have been conducted to highlight the impacts of armed conflicts and warfare on landscape change and to raise the concerns about the environmental repercussions of wars. For example, warfare activities such as the ruin of the landscape through carpet bombing, the slaughter of wildlife, and the use of biocides can all have disastrous environmental consequences and cause human suffering [8][9][10]. Warfare training activities, while essential for military readiness, can have severe environmental consequences, including vegetation removal and alteration, accelerated rates of soil erosion, increased soil compaction, and loss of wildlife habitats [8][9][10].
The role of armed conflicts in forest degradation is still not clear. Although wars almost always have devastating effects on people, they might also have negative and positive impacts on forests [11,12]. For example, armies burn forests to clear space and improve visibility to spot enemies and for soldiers to hunt wildlife for food [13,14]. Modern wars, particularly intranational conflicts, often play out in remote areas where armed factions seek the cover offered by deep forests, mountains, and other landscapes [12,15]. Protected-area boundaries lose effectiveness in this context, usually resulting in the evacuation of field staff and the suspension of conservation activities [16].
Some of the examples of positive consequences of wars on forests are that they may discourage the conversion of forests into pastures, as ranchers may fear being abducted or losing their livestock. Similarly, timber companies may limit their investment in such areas to avoid risking damage or loss to their expensive machinery [12,15].
Large concentrations of refugees and displaced people can put huge pressure on the local environment as they move to new areas to hunt, fish, and remove trees for building their houses and firewood, rapidly depleting such resources [13,17]. The impacts appear to be shaped by the use of natural resources near settlements, following large displacements of residents who had previously dispersed across the landscape during the war [18].
In Afghanistan and during the Afghan-Soviet war, numerous air raids caused massive and long-lasting wildfires in the woods of Paktia and Kunar [19,20]. Overall, the destruction of the forests was the worst environmental consequence of this conflict, with the forest area declining from 3.4% to 2.6% of the land area in less than ten years [19,21]. Refugees and displaced people in Sudan and South Sudan cut down vast areas of forest to obtain fuel supplies and construction materials in the subhumid zones in southern Darfur and Kordofan and over 1 million hectares of forest land were cleared between 2005 and 2010 [14,22,23].
The Mediterranean region has historically been a center of geopolitical warfare. Its distinctive geographical landscape at the crossroads of three continents has connected and divided the domains of influence of major powers in Europe, Africa, and Asia since ancient times. Securitization persisted beyond the East-West conflict and into the twentyfirst century, with its complex geopolitical landscapes and shifts connecting numerous scales, players, and security elements [24]. The Mediterranean region experienced about 18 armed conflicts between 1952 and 1995, with 12 of them taking place in North Africa and the Middle East [6]. Several long-term conflicts have made the Middle East and North Africa two of the most unstable areas in the world [25]. Some of these conflicts had an impact on landscape change. For example, Algeria lost about 2 million hectares of forests during the French occupation and 650 thousand hectares of forest were burned between 1956 and 1960 [26,27]. The repetitive armed conflicts in the North of Lebanon were associated with the degradation of 1020 hectares of vegetation. The rate of negative change in plant cover between 2006 and 2007 was at its highest [28].
In March 2011, the devastating armed conflict that erupted in Syria left more than 13 million people in need of assistance, causing the migration abroad of more than 5 million people and displacing more than 6 million people internally in the country [29]. About 1.4 million refugees fled to the Syrian coastline region [30], which comprises 90% of the country's forests, with forest ecosystems that are vital for carbon storage, timber, biodiversity, and recreational uses [31,32].
Remote sensing technology has become a vital tool in mapping, monitoring, and managing forest resources due to advanced data analysis techniques and improvements in remote sensing science and ecosystem modeling [33]. Prior to the launch of the first satellite in 1972, ecologists were restricted to field-based methods for addressing degradation in tropical forests [34]. Even in the early stages of satellite imagery programs, data accuracy was limited, but, by 1980, satellite imagery was able to differentiate between tree species and forest stands with similar spectral characteristics based on individual habitat preferences [35]. Machine learning approaches, such as the Random Forest (RF) and Support Vector Machine (SVM) algorithms, which can predict target variables with high accuracy, are powerful and commonly used ML techniques in this field. Combining remotely sensed data with field inventories provides the most consistent estimates [36][37][38][39][40].
Remote sensing can be used as a cost-efficient option for assessing degradation through proxies such as canopy cover percentage, with a decreasing trend indicating forest degradation [41]. In conflict areas such as Syria, where field-based methods are challenging due to war conditions, remote sensing by satellite images is an appropriate instrument for analyzing the impact of the conflict on the landscape. Land cover classification with satellite imaging time series has been shown to be effective for mapping forest, grassland, and agricultural dynamics [42][43][44][45][46].
In this study, our main objective is to analyze the change in forest cover in western Syria and quantify the recent changes in forest dynamics due to conflict-related drivers, answering the following questions: (a) Has the forest cover of Syria changed during the ongoing armed conflict (2010-2019)? (b) What are the patterns and spatial distribution of the forest cover change? (c) What impact do conflict-related drivers, such as refugee displacement, explosives, wildfires, and proximity to roads and settlements, have on the dynamics of forest cover?

Location
The study area includes the western part of Syria, which extends over an area of 37.512 km 2 and lies between 32°31′ to 36°84′ N and 35°71′ to 36°97′ E, with a coastline that extends for 183 km along the Mediterranean Sea ( Figure 1). It is important to note that the study area excludes the Golan Heights in the southwestern part. This region is recognized by the United Nations as Syrian territory occupied by Israel for over 39 years.

Climate
The climate in western Syria is considered a Mediterranean climate, characterized by dry hot summers, cold winters, and concentrated rain periods in spring and autumn [47]. Average temperatures in winter are moderate to cold, with precipitation between 100 and 1400 mm/year. In summer, the temperatures can rise up to 30 °C in most regions and, in some places, can go above 40 °C depending on the bioclimatic range [47].

Forests
According to data provided by the Syrian Ministry of Agriculture in 2005, the combined forest and woodland areas in the country cover approximately 0.6 million hectares, representing 3% of the total land use [48]. Among the various types of forest vegetation, oak forests dominate, with 59% of the total forested area, followed by coniferous forests, with approximately 27.5% [48]. The remaining species are considered associated species [48].
Ref. [49] classified the forest vegetation in five ecoregions according to the different geoclimatic zones in the country: the coastal mountains, the Baer and Basit Mountains, the Aleppo Mountain range, Hermon Mountain, and Jabal al-Arab Mountain; four of which fall within our study area: 1. Most of the forest cover in Syria is located in the western coastal mountains, which are divided into two parts due to the variation of precipitation that they receive annually: a.

Landsat Imagery
Models to identify forested and non-forested areas were trained from photointerpretation train samples and the spectral signal of the land cover. A set of Landsat TM, ETM+, and OLI surface reflectance images acquired during the spring season were retrieved and images were processed via the GEE environment [50]. We utilized the LandTrendr algorithm to aggregate data on a yearly basis using the annual medoid approach. By fitting a segmented model, we achieved a significant reduction in noise caused by variations in illumination conditions, phenological changes, atmospheric conditions, and geometric registration, offering an enhanced change detection capability [51,52]. We used the algorithm LandTrendr to aggregate data on a yearly basis, by means of annual medoid, and fit a segmented model. The LT consisted of a set of spectral temporal segmentation algorithms that isolated changes in time series. In general, the original time series trajectories were repeatedly simplified by identifying the most significant vertices of change. After potential vertices were identified and linked using flexible fitting rules, various segments of the updated trajectory were removed and remodeled to finally select the best model [51,53].
A smaller subset of variables was created by synthesizing the initial set of 6 spectral bands using the Tasseled Cap (TC) transformation [54]. TC is an orthogonal transformation of spectral data that is based on the three resultant axes, brightness (TCB), greenness (TCG), and wetness (TCW) and, concretely, TCW and TCB are highly correlated with forest structure [52,55]. The three components were transformed using [56] weighting coefficients. Also to mention is that the OLI data were homogenized with TM and ETM+ spectral data using the coefficients described in [57].

Palsar Data
For a better definition of forest areas, we used data from the L-band radar sensor PALSAR and PALSAR2, with 25 m spatial resolution. Concretely, the aggregated annual mosaics were available in GEE [58]. As only mosaics for the years 2010 and from 2015 to 2019 were available, then the base year was set as 2010 and we assessed the changes from 2015 to 2019. The product was provided in digital numbers (DN) and without speckle reduction. Image preprocessing was required to convert DN to gamma naught values in decibel units (dB) (Equation (1)) and a despeckling process using the refined Lee speckle filter [59].
Then, the available polarizations Horizontal-Horizontal (HH) and Horizontal-Vertical (HV) were selected and Gray Level Co-occurrence Matrix (GLCM) statistics were applied in order to compute textural metrics and identify patterns [60]. From the set of available metrics, we selected sum average, entropy, correlation, and variance applied in an array of 4 × 4 neighbors. Finally, all HH and HV original and textural bands were resampled to 30 m, using bilinear interpolation, and snapped with Landsat imagery.

Response Variable
In order to obtain "ground truth" information to discriminate forested and non-forested areas, we carried out a photointerpretation of the available imagery in Google Earth for the year 2019 and we restricted forest samples to where the trees were higher than 3 m, using a global forest height map [61]. The main goal was identifying those areas to serve as train and test samples for the evaluation of model performance. From the collection of imagery, we were able to identify a set of 3000 locations, including 3600 sites of forested areas and 2400 non-forested locations. Please see (Table A1, Appendix A) for more detailed information about the train and test subsamples.

Variable Selection and Modeling
To avoid redundancy and collinearity between spectral variables, we performed a variable selection analysis. To do so, we used the strategy based on the recursive variable elimination. This method refits Random Forest models [62] by eliminating variables one by one. Then, whichever number of variables yields the lowest error is evaluated and ranked by variable importance. The final selection results from a bootstrap process where the variable selection strategy is repeated and aggregated with 30 different subsamples to ensure the reliability of the variable selection.
To identify and detect forested and non-forested areas, we applied the Support Vector Machine (SVM) classification algorithm with previously selected spectral and derived variables from the year 2019. SVM is a non-parametric machine learning algorithm that has been specifically designed to analyze and recognize patterns based on building a series of hyperplanes in a high-dimensional space in an effort to find the greatest separability across classes [63,64]. SVM, concretely, the SVM-Radial Basis Function (RBF) kernel, is widely used in land cover classification models due to its good performance [65,66]. Cost and gamma parameters were adjusted by experimenting with all possible values between 1 and 700 and 0.001 and 1, respectively. A 10-fold cross validation was performed for model optimization using the "Caret" R-package [67]. Additionally, we produced spatial predictions and outputs using the "raster" package [68] and evaluated their performance using Cohen's Kappa index [69] and the accuracy percentage.

Deforestation Detection
Deforestation was detected by predicting with the model over the previous years and obtaining an annual mapping of forested and non-forested areas. We defined deforestation as a pixel where there was change from forest in the year 2010 to non-forested in a later year and, furthermore, continued to be detected as non-forested until the end of the time series in 2019. This filtering process avoided overestimating the phenomenon when detection was mainly produced by noise during data acquisition. This procedure provided a yearly mapping of deforestation.

Auxiliary Analyses
Deforestation is induced by human activity and is a consequence of war. For this reason, we described the spatial relationship between locations and infrastructures that can facilitate access to forests and lead to deforestation such as urban settlements, refugee camps, and roads by calculating the Euclidean distance between deforested areas and these areas. We summarized this information by means of the distribution of different ranges of distance and descriptive statistics. Additionally, forest fires from natural causes or induced by the explosions of war can contribute to deforestation. For that potential driver, we analyzed the relationship between deforested areas and forest fires. Furthermore, we analyzed whether deforestation occurred mostly in areas with higher or lower explosion activity. All of these auxiliary data were retrieved from different sources. Urban settlements and roads were extracted from OpenStreetMap [70] and refugee camp locations from data collected by U.S. Department of State (U.S. Department of State Humanitarian Information Unit, 2016), Humanitarian Information Unit [71]. For the forest fires database, we retrieved all images between 2015 and 2019 from MCD64A1 version 6 Modis product [72]. Explosion intensity data were retrieved from georeferencing, followed by an unsupervised 4-class classification that grouped the different colors in figure (Intensity of Explosive Incidents) published in the United Nations humanitarian needs report [73].

Land Cover Models Performance
SVM models usually perform well in deforestation monitoring; our most successful variable combination was obtained selecting all variables, but, as can be observed in Figure 2, the accuracy values were stabilized using five variables. To simplify and reduce the model complexity, we used the five most influential variables: TCW, SAR HV polarization GLCM correlation, SAR HH polarization GLCM correlation, SAR HH polarization GLCM variance, and HV polarization GLCM sum average (Table A3, Appendix A). SVM land cover model attained a predictive overall accuracy of 0.91 and a Kappa index of 0.83 (Table  A2, Appendix A). This high performance allowed us to make reliable predictions for different years and then to synthesize all the temporal information about deforestation.

Predictive Performance of the RF Models
The predictive performance of the model was satisfactory; the SVM model yielded a high AUC 0.84 (Figure 3). We observed a clear gap (Kruskal-Wallis test significant; p < 0.05) (Figure 4) that separated the predicted probability between deforested and non-deforested observations. The density plots corroborated the high contrast between categories ( Figure 4). In non-deforested locations, we found a frequency peak close to 0.1, after which the curve flattened. In deforested locations the probability peaked near 0.9, decreasing fast into low density values below the 0.8 threshold. Overall, all diagnoses provided proof of the satisfactory performance of the models, which supported our findings.

Overall Deforestation Patterns
The total loss of forest cover over the study period (2010-2019) was 19.3%, corresponding to a loss of 63.7 thousand hectares. The highest rate of forest decline was observed between 2010 and 2015, during which the total forest cover dropped by 11.5% (or 37.9 thousand hectares). During the period 2015-2018, a gradual but consistent decline was also observed, accounting for an accumulated reduction in forest area of 7.8% (See Table 1 and Figure 5). The analysis of satellite images from nine different years spanning from 2010 to 2019, with 2010 being the reference year, revealed a persistent decline in forest cover in the study area. Deforestation occurred mainly in the northwestern area, especially the coastal mountain range, and around Damascus in the South.

Continuous Variables: Distance to Roads, Urban Settlements, and Refugee Camps
We examined the relationship between road proximity to forests and deforestation. Results showed that the average distance between the deforested areas and the roads was 0.25 km and that the majority of deforested areas were within 0.04 to 0.39 km away from roads. A graph was created to show the marginal effect of road distance on deforestation probability ( Figure 6). According to the findings, the probability of deforestation was higher in areas where forests were closer to roads. For urban settlements, the average distance of the deforested areas was found to be 21.62 km and the majority of the deforested areas were between 12.6 and 30 km away from urban settlements. This indicated that the deforested areas were located much further away from urban settlements, with a wider variability in distance compared with roads.
For refugee camps, the average distance of deforested areas from the refugee camps was determined to be 28.55 km. However, higher rates of deforestation were observed to occur within a much smaller radius. The probability of deforestation was higher in areas that were located closer to refugee camps.

Discrete Variables: Intensity of Explosive Events
By overlaying the distribution of the deforested areas with the areas that experienced explosive events such as bombing and shelling, we found that over half of the total deforested land (53.7%) was located within the areas that experienced the highest level of explosive intensity, mainly in the north and northwestern parts of the study area. Additionally, 12.5% and 10.9% of the deforested area were found in regions with medium-intensity explosive events, while 20.3% of the deforested land was found in regions with low-intensity explosive events. Only 2.62% of the deforested land was found in areas that did not experience any explosive events during the study period (see Figure 7).
According to the analysis of MODIS fire data spanning from 2010 to 2019, a total of 15.24 thousand hectares of forest cover were burned, which constituted 23.4% of the overall decline in the forest cover in the study area (See Figure A2, Appendix A). The wildfire incidents were mainly concentrated in the northern and northwestern parts of the study area, corresponding to the regions with the highest deforestation rates, as shown in Figure  7.

Discussion
The models used in this study proved to be highly accurate in predicting deforestation, with an AUC value of 0.84 and density plots showing a clear distinction between deforested and non-deforested areas. The SVM land cover model achieved an overall predictive accuracy of 0.91 and a Kappa index of 0.83, which was consistent with previous studies that used machine learning algorithms such as Random Forests to predict deforestation. For instance, [74] used SVM and MD classifiers to detect the deforestation caused by oil palm plantations in two locations in Peninsular Malaysia. The results indicated that combined Landsat and PALSAR data provided better accuracies than either data source alone for both study areas. Additionally, [40] found that the SVM algorithm was more accurate than other machine learning models, as evidenced by its high AUC value of 0.90.
Our results indicated that, during the period of the armed conflict, the study area underwent a severe reduction in forest cover, resulting in a loss of 19.3% in 9 years from 2010 to 2019. In contrast, prior to the armed conflict, the forested areas in Syria exhibited a positive trend, with annual growth rates indicating a positive increase of 0.48% from 1997 to 2007, as reported by the Syrian Monitoring of Agricultural Resources [75]. The Food and Agriculture Organization also reported an increase of 119 thousand hectares in the forest cover in Syria from 1990 to 2010; this increase was mainly caused by the afforestation efforts made by the government to expand the forest cover in the country [76].
Before 2011, Syria used to produce around 30 million tree seedlings annually and the seedlings were dedicated to afforestation campaigns around the country. However, due to the conditions of the war, this number has decreased significantly to 1.5 million seedlings according to statistics from the Ministry of Agriculture [77].
The ongoing armed conflict has resulted in widespread poverty and increased financial strain on individuals, as well as a significant impact on the ability of the population to access basic necessities such as food and other essential items due to the sharp and unprecedented increase in prices [78]. This has created an economic crisis at both the local and national levels in Syria [78]. As a consequence of this economic crisis, sustained pressure has been applied on forest cover for several years, primarily due to the heavy logging needed to meet the high demand for firewood and charcoal from the local populations and surrounding urban areas. Oaks are the preferred tree species for logging according to loggers, as they prefer oak wood due to its high density and suitability for charcoal production [79]. This demand has increased due to the high cost of fuel and the severe cold during the winter months, which further has exacerbated the economic struggles of the population [80]. The relationship between poverty and deforestation has been extensively studied by researchers. Ref. [81] found that illegal logging in Central and Eastern Europe was explained by poverty, as well as weak law enforcement and land ownership reforms. In Cameroon, poverty was also identified as one of the factors underpinning illegal logging, along with conflicts and systemic corruption [82]. The informal sanctioning of logging activities of the rural poor, known as "humanizing the law", has been observed in the Philippines as a strategy to avoid aggravating rural poverty and fueling civil insurgency [83,84]. The lack of alternative livelihood projects has been identified as a factor contributing to illegal logging and providing such projects has been suggested as a means of controlling it [84].
Our results also suggested that the areas with higher deforestation rates were found to be closer to the roads. This result may have several possible explanations. One possible reason is that roads provide more access to forested areas, making it easier and more efficient for people to cut down trees for firewood, charcoaling, or agricultural expansion. Another reason could be that the presence of roads leads to increased development and population growth, which in turn leads to more demand for resources and land use. This could be particularly true in the western coastal part of the study area, which has received more than 1.4 million internal migrants since the beginning of the conflict [80].
Ref. [85] found that deforestation rates were significantly higher near roads and rivers compared with other areas in the Amazon, with nearly 95% of all deforestation occurring within 5.5 km of roads. This result was consistent with our study, but our areas with higher deforestation rates were located in the most accessible areas (0.25 km on average, 0.04 to 0.39 km away from roads). Similarly, [86] found a positive and significant correlation between road proximity and deforestation, indicating that the closer a forested area was to a road the more likely it was to be deforested. This relationship has often been driven by the reduced transportation costs associated with roads, which can incentivize increased economic activities and attract newcomers to the area, leading to further deforestation.
On the other hand, our study revealed that deforestation was more likely at greater distances from urban areas compared with roads. The highest rates of deforestation occurred between 12.6 and 30 km away from urban settlements. Interestingly, the majority of the deforestation was concentrated in the western part of the study area, where forests are located far from urban centers, especially on the western coastal mountain range. We suggest that this pattern may be attributed to the fact that loggers harvest wood in these remote forests and transport it to be sold in the city centers. It is also possible that the distance to urban settlements cannot be seen as an isolated factor contributing to the deforestation but is combined with other drivers such as roads that provide accessibility to transport from the forests to the wood markets.
The presence of refugees in camps can lead to deforestation, especially due to the demand for firewood. The collection of firewood for cooking and heating can result in the degradation of forest areas; firewood collection is often driven by income constraints and when the availability of substitutes for cooking fuel is limited [87]. For example, in the Rohingya camps in the southeastern border region of Bangladesh, the dense forest area decreased from 8.53 ha in 2016 to 4.49 ha in 2018 and the refugee settlement area increased from 271 ha in 2016 to 2679 ha in 2018 [88].
The use of firewood as an energy source has been a major issue in Syria, especially in besieged areas such as Idlib (northwest Syria). Nearly 90% of displaced people lack access to energy, which leads to numerous problems. The lack of access to energy makes it challenging for Internally Displaced People (IDPs) to meet their basic needs, including cooking and heating, which can lead to further environmental degradation as people turn to sources such as wood and use it in an unsustainable manner. The combination of these factors has resulted in widespread and illegal deforestation in the region, causing longterm environmental and humanitarian consequences [89]. The situation in Idlib and parts of Aleppo is complex and displacement appears to be a significant factor contributing to illegal logging. The massive influx of Internally Displaced People (IDPs) into the northwest provinces has resulted in the clearing of trees and the conversion of agricultural lands into temporary camps. Moreover, IDPs have resorted to cutting down trees to provide themselves with heating due to the soaring fuel prices and the bitterly cold winters that are prevalent in northern Syria [90].
The results of our study are consistent with previous research, indicating that deforestation rates are highest in areas closest to refugee camps (<9 km). Additionally, we found high deforestation values 28.5 km away from the refugee camps. Furthermore, we observed that refugee encroachment is unlikely to be the only factor causing deforestation, as the impacts extend beyond 28.5 km to up to 40 km from the camp. Hence, we suggest that other factors may contribute to deforestation, highlighting the complexity of this relationship between refugee camps and land-use change.
According to the United Nations Mine Action Service (UNMAS), Explosive Ordinance (EO) has been used in Syria through air and artillery strikes with rockets and mortars, as well as Improvised Explosive Devices (IEDs) and landmines, which have negatively impacted the civilian population [91]. Explosives not only have a negative impact on civilians, but they can also lead to forest fires. For example, the conflict in Ukraine has resulted in forest fires due to artillery shelling and bombing; in 2014, over 200 forest fires were caused by these actions in the Donbas region [92]. During the First World War, Germany, France, and Britain fired 1.45 billion shells, leading to massive deforestation in France, Belgium, and the surrounding regions [93]. Turkish forces have also set deliberate forest fires to diminish coverage for Kurdish fighters in the ongoing conflict with the Kurdistan Workers' Party (PKK) [94].
Ref. [95] also reported that forest fires in the northern Al-Kabeer River basin of the northern Province of Latakia in Syria in 2012-2014 were mainly caused by the war and the firing of shells and bombs in the area. These findings support our conclusion that shelling and bombing can be major causes of forest fires and deforestation.
In addition to these causes, charcoal production by locals in forested areas could be another contributor to forest fires. This process involves "charring" or partially burning tree branches and firewood in a shallow pit covered with insulating cloth, which can be dangerous as it may lead to the spread of fire due to sparks generated by the lit wood, especially during the dry season [79]. Over 1500 licensed and unlicensed charcoal fire pits are believed to exist in the countryside of the coastal region in Syria, with profits that can reach 300%. This illegal trade serves as a vital source of income for the local population [96].
Based on our findings, forest fires were responsible for 23.4% of the forest cover decline in the study area, mainly concentrated in the northwest. Over 50% of the deforestation occurred in high-intensity explosive event zones. Shelling, bombing, and charcoaling were possibly associated with forest fires and deforestation in the study area.
While this study identified the most important war-related drivers of forest cover change, it did not fully explore the social and political factors that contribute to this phenomenon, as it lacked the possibility to acquire field data due to the ongoing conflict.

Conclusions
We used remote sensing and machine learning methods to quantify the spatial and temporal change of the forest area over the course of nine years (2010-2019) of war in Syria. According to our results, a substantial decline in forest cover in the study area was found during the armed conflict, with a total loss of 19.3%.
The armed conflict has had both direct and indirect effects on deforestation. The direct effect of the conflict on deforestation was driven by an escalation of violence, resulting in a high number of internally displaced individuals who relied on the forest for heating and cooking, given the absence of other energy sources. The intensity of explosive events in forested areas contributed to deforestation by causing fires, resulting from artillery shelling and bombing.
Indirectly, the conflict drove deforestation by increasing poverty, which led to greater reliance on the forest as a more affordable alternative fuel source and an additional source of income for local communities, with the absence of forest guards and forest law enforcement during the war. Furthermore, the roads played an important role by facilitating the access to forest resources by linking people to the forest and local wood markets in adjacent urban settlements.
This study provides valuable insight into the impact of armed conflicts on forest cover. It highlights the need to understand the role of war in post-war restoration and conservation efforts. In the post-war period, it is crucial for local governments to provide sustainable substitutions for heating and cooking. This would reduce the reliance on forests for fuel and prevent further deforestation. Governments should also enforce laws and regulations to protect forests and combat illegal logging. Additionally, reforestation efforts should be prioritized to restore damaged areas and promote the growth of new forests. This could include the planting of native species and the establishment of protected areas. The involvement of local communities in these efforts is also important, as they have a vested interest in the preservation of their local environment.