Examining forest cover change and deforestation drivers in Taunggyi District, Shan State, Myanmar

Myanmar has been experiencing a significant amount of deforestation and forest degradation in recent years. Being a developing country, people are heavily dependent on its forest for sustenance and livelihood. This study examines a methodology to identify potential drivers and their relative significance for deforestation. The study was tested in one of the districts but could be applied in other areas of the country. The forest and non-forest land cover maps from the Japan Aerospace Exploration Agency (JAXA) for the years 2008 and 2016 were used in the study. It was derived that 46.54% of study area is still covered with forest, but there has been a significant decrease in forest area by 7.29% between the years 2008 and 2016. We examined a number of spatially explicit potential drivers of deforestation such as infrastructure, elevation, slope, deforested land, and population. As informed prevention awareness of deforestation, we projected future forest conditions using a cellular automation modeling technique for the years 2020, 2025 and 2030. We found that major physical and socioeconomic driving factors of deforestation such as proximity to infrastructure (reservoirs and roads), certain elevation levels, slope, proximity to previously deforested area and population density are strongly associated with neighborhood deforestation. The future projection showed a decrease in forest area by 13.8% from 2016 to 2030. This work therefore provides crucial information on forest landscape for forest management in the district. The projective scenario of study area generated by the model highlights the need for forest conservation and planning while addressing the key drivers of deforestation, giving direction for future potential areas of REDD+ implementation in the region.


Introduction
Forests have extraordinary and far-reaching contribution for the wellbeing of humankind and biodiversity. It plays a fundamental role in providing food security, combating rural poverty and providing livelihood opportunities to a large population of forest dependent community. Forests provide a range of long-term ecosystem services like clean air and water, serve as carbon sinks and store, conserve biodiversity and mitigate climate change (FAO 2015). Globally, the world's natural forest area has undergone major change since 2010 by a decrease of 6.5 million ha per year (FAO 2015). Currently, deforestation and forest degradation is considered to be the second major contributor of greenhouse gas (GHG) emissions after the fossils fuels (IPCC 2008). Subsequently, Myanmar, a country with about 45% of land covered with forests supporting variety of distinct vegetation composition and wide range of faunal diversity, is undergoing unprecedented changes in forest quality and quantity in the recent years. Annually, 546 thousand ha of forest is estimated to be deforested in Myanmar and has been ranked to have the third highest forest loss in the world after Brazil and Indonesia (FAO 2015). Such substantial loss may bring devastating outcomes to human lives and the overall ecosystem in the country. With decreasing forest and population growing by one million annually, pressure for increase in food and firewood consumption is inevitable. Myanmar has increased 7.73 million ha of agricultural field from 1990 to 2011 which is the direct result of change in land use pattern. Shan State, a state covering 1/4th of Myanmar, had the largest net forest loss of 5647.7 km 2 , for the year 2001-2010 (Wang and Myint 2016) and was responsible for emitting 6.86 million tons of carbon annually (FAO/RECOFTC 2016).
To reduce the effect of GHG emissions from forests at the same time, contribute to reducing poverty and sustainable development, "Reducing emissions from deforestation and forest degradation" (REDD+) developed by parties to UNFCCC has been introduced in many developing countries including Myanmar. The Measurement, Reporting and Verification (MRV) activities of REDD+ include timely measurement of forest stocks, monitoring of forest loss and estimating forest-related emissions, among a few. For this, earth observation and geospatial technologies can provide spatially explicit datasets and products useful for monitoring forest cover, identification of drivers of deforestation and developing forest reference emission levels. Forest cover data products with global coverage like GlobeLand30 (Chen et al. 2015), the Global Tree Canopy Cover (Hansen et al. 2013), the JAXA ALOS-PALSAR forest/non-forest map (Shimada et al. 2014), and the Landsat Vegetation Continuous Fields (Sexton et al. 2013) could provide relevant secondary data for study and analysis (Estoque et al. 2018).Various geospatial modeling techniques used with these datasets can provide understanding of complex process of forest cover change, identification of the determinants/drivers of the change and projecting future scenarios (Thapa et al. 2013). Spatial modeling methods, logistic regression (Chowdhury 2006;Echeverria et al. 2008), genetic algorithm (Soares-Filho et al. 2013;Venema et al. 2005), weight of evidence (Maeda et al. 2011;Soares-Filho et al. 2010;Thapa et al. 2013), and artificial neural network (Khoi and Murayama 2011;Mas et al. 2004) have been used in the recent past to study complex process of forest change and associated biophysical and socioeconomic drivers of deforestation. One of the prominent methods includes Cellular automata (CA) with the weights of evidence (WofE) modeling framework that is based on Bayes' rule of probability with an assumption of conditional independence (Thapa 2012). It identifies the interaction between changes in forest cover and the various drivers such as physiography, socioeconomic factors, forest management and developmental policies. It has been widely used for modeling changes in land use and land cover (Bonham-Carter 1994).
Being high on deforestation rates (i.e., 1.8% for 2010-2015), Myanmar has significant prospects to develop conservation strategies and to study biophysical and socioeconomic drivers of forest change at local level (FAO 2010). Case studies of forest change in the country can be helpful for gaining an understanding of complex relationships between social and natural systems that drive landscape change, which is not previously done in the context of Myanmar. Therefore, in this research, we aim to study the forest cover changes and develop a method to examine the underlying drivers of deforestation in Taunggyi District of Shan State, Myanmar. Emphasis has been laid in spatial analysis of forest changes and examination of spatially explicit drivers of deforestation. In addition, future forest cover has been projected under Business As Usual (BAU) scenario for informed prevention awareness for forest management. Readily available global land cover maps were assessed; the one with highest comparative accuracy was used for the study. We employed CA model based on WofE framework to quantify the significant drivers of change and produce future trajectory projection which has the potential to be applied to other parts of the region and upscaled to country level to produce broader results, and will be useful in the study of deforestation trends, policy framing and forest management.

Study area
Taunggyi is a district located in the Southern Shan State of Myanmar (Fig. 1). The district covers an area of 24,239.10 km 2 extending from 96°9′46″E, 22°14′45″N to 97°30′38″E, 19°19′3″N. The district is home for 1.7 million people where 27.3% of population living in urban areas (Ministry of Immigration and Population 2015). Physiographically, the district is a part of Shan Plateau that ranges between 100 and 2500 m. The topography is divided into Shan Plateau, Myelat and Plains. Myelat and Plains cover 25% of the district, and the rest are high plateau with mountain ranges, hills and spurs. Plains cover Kalaw Township, Nyaungshwe and Hopone Townships. Ywangan, Pindaya and Pinlaung Townships are Myelat region, and the rest are Shan Plateau. Mountain ranges lie from north to south in Taunggyi District, and the rivers and streams also flow in the same direction. The forest type identified in the district was Moist Upper Mixed Deciduous Forest, Dry Upper Mixed Deciduous Forest, Evergreen Forest, Dry Hill forest, Indaing forest and Pine forest. As per land use status of the district in 1996 there are 18 Reserve Forest covering an area of 30.12% of the district followed by 11 protected public forests (2.98%), 11 protected area system (5.19%) and 12 unclass forest(19.60%) in the district (MONREC 2016). Forests in Taunggyi comprise some of tree species with high commercial importance such as Tectona grandis, Dalbergia cultrata, Shorea obtuse, also IUCN "near threatened species" of Dipterocarpus tuberculatus found in the open and dry forests in the district that is mainly being used for firewood purpose by the locals (Field survey, December 2017).

Assessment of land cover data
We accessed and compared three different land cover maps for the study area. These maps included land cover map provided by Japan Aerospace Exploration Agency (JAXA), forest-non-forest (FNF), SERVIR-Mekong regional land cover map (RLCM) and European Space Agency (ESA) Climate Change Initiative (CCI) for the year 2016. We reclassified these maps into forest and non-forest. As these maps are in different spatial resolutions, we have resampled them at 30 m for the assessment. For the preparation of reference data, 4000 random sample points at 99% confidence interval with 2% error margin were generated across the spatial extent of Taunggyi District. These points were then classified into the forest or non-forest category through inspection in Google Earth Pro. We used the 2016 image in Google Earth to categorize our reference points and used the image of the closest year in case of non-availability of image for the same year. JAXA FNF map depicted 78% overall accuracy which was highest as compared to 70% and 66% overall accuracy of ESA CCI Global land cover map and SERVIR-Mekong RLCM map.

Determining deforestation driver variables
For modeling deforestation trends, it is important to select best set of input variables in order to produce the most optimum fit between the empirical data and observed reality. Therefore, we identified a set of explanatory variables governing the changes from forest area to non-forest, for calibrating deforestation model and quantifying drivers leading to deforestation. Given this, it must be noted that there are no universal explanatory factors required for producing ideal simulation results and the set of variables used may not necessarily represent all the factors responsible for the event. The variables may differ from one landscape to another but similar factors found in many studies conducted by Thapa et al. (2013), Carodenuto et al. (2015), Chowdhury (2006), Kolb et al. (2013), Seto and Kaufmann (2003), Thapa et al. (2013), Elz et al. (2015), etc. We examined factors such as elevation, slope, proximity to previously deforested area, proximity to roads, proximity to waterways, proximity to newly constructed reservoirs, proximity to settlements and local-level population density (Table 1). Physical conditions such as elevation, slopes and rivers play an important role in land use changes. Elevation differences promote terrain differences that are associated with complete or selective clearance of forest in Myanmar. Flat land promotes clearance of forest for agricultural needs, while higher elevations associated with rugged topography in most cases make it difficult for clearance of forest stands, especially in the study area. Regions of 100-300 m altitude in this region are often associated with shifting cultivation and conversion of forest land to rubber and palm oil plantation (Myint 2018;Miettinen et al. 2014). Some slopes are subjected to waterlogged conditions, whereas others provide ideal drainage favoring cultivation in turn contributing to forest clearance as part of agricultural expansion. Construction of roads has opened up many intact forests for settlement and resource extraction mostly in the tropical areas. According to the study by Alves 2002, 90% of deforestation in Amazonian forest occurred within 100 km of road. Myanmar has been increasing its road network given, the various development projects, growth in foreign trade and population density. Most of the roads being constructed since 2005 in Myanmar tend to be linked with areas of new rubber plantations, mines and hydroelectric project sites built at the cost of forested areas, especially in states of Kachin and Shan. The further expansion in road network is projected to have adverse effect on forest, especially the ones being constructed around heavily forested border areas (UN REDD 2017). Hence, proximity to roads is an important determinant that encourages deforestation arising from higher human accessibility and chances of forest encroachment. Waterways of all forms like navigable rivers, streams and canals have positive correlation with deforestation. According to another study by Barber et al. (2014), nearly 95% of all deforestation in Brazilian Amazon occurred within 5.5 km of roads or 1 km of river. Waterways are largely used for the transportation of logged materials hence forested river banks in tropical areas are highly favorable for tree-felling activities. Construction of canals and related anthropogenic activities in the forested regions lead to deforestation and threaten the biodiversity of the region. As the population of Myanmar grows by approximately 1.5% annually, with fuel wood and charcoal consumption being 76.41% of total energy consumption in Myanmar, higher population density would directly put higher pressure on available land resources for food and livelihood opportunities (Sein et al. 2015). Hence, population density of the region is an important driver of deforestation. Model calibration and simulation were done in Dinamica EGO software that uses Markov chain approach for calculating landscape transition matrices between the years (https ://csr.ufmg.br/dinam ica/). The effects of the spatial factors on the landscape change and probability map of deforestation for those factors were calculated using the WofE method which is based on Bayesian probability theory. It takes into consideration the deferential effect of spatial factors for predicting deforestation (Soares-Filho et al. 2010). The WofE values represent the attraction or repulsion of deforestation (event) based on the effect of the driving factor to deforestation. Higher values of WofE (greater than 0) depict higher probability of deforestation, whereas negative values represent lower deforestation probability ranges. The variables were tested for spatial dependence using Cramer's V (Bonham-Carter 1994) which shows the values ranging from 0 to 1 depicting the degree of association between the driver factors. Cellular automata (CA) model integrated with WofE was then used to calibrate deforestation pattern across the district and build BAU scenarios for upcoming years.

Transition potential map calculation and simulation
The model was prepared through the standard calibration, simulation and validation method. For calibration, the WofE ranges and the coefficients for each of the driving factors concerning the event, i.e., change of land cover class from non-forest to forest, were calculated. The WofE coefficients calculated through transitional conditional probability are significant in understanding the influence of each these drivers. Due to lack of spatial information on infrastructure development plan, this research assumed that existing infrastructure will be constant and influence equally overtime. Therefore, all of these drivers except proximity to previously deforested area were set static during the model run. The simulation model was then produced by setting internal parameters of patch mean size, variance and isometry. Since the initial year of land cover map used was 2008, the model was iterated 8 times to produce a simulated land cover map of 2016. The model validation included finding change similarity between the observed 2016 and simulated 2016 map using exponential decay function in Dinamica.

Spatial analysis of forest change
JAXA forest-non-forest map has been analyzed for the year 2008 and 2016 to study the deforestation over the years. Analyzing the map, it was seen that forest cover in Taunggyi significantly decreased between the years 2008 and 2016 by 1769.28 km 2 which is 7.3% of the forest area. Non-forest category mainly consisting of agriculture, built-up land and tree cover < 10% has shown an increasing trend from 10,905.83 to 12,624.17 km 2 at the cost of forest land over the years. The water category also depicts an increase of 50.94 km 2 mainly because of built-up of new artificial reservoirs in this region in the last years ( Table 2).
Conversion of land cover from one class to another was also looked upon (Fig. 2). The highest area (3005.84 km 2 ) of inter-land cover class conversion occurred from forest to non-forest. Forest-to-non-forest conversion is observed, particularly in the southern areas of the district. Also, area of 54.27 km 2 of forest area is seen to be converted to water, significant portion of it because of the newly constructed reservoirs in the district.
We also analyzed the township-wise change in forest cover in Taunggyi District. The highest change from forest to non-forest category is seen in the township of Pinlaung where forest decreased by 893.345 km 2 at the same time non-forest type has increased to a significant amount (Table 3). The water category has also increased by a significant amount compared to 2008 year mainly accounting to construction of artificial reservoirs in the district. Reservoir in Thabyegon (Pinlaung Township) has resulted in the clearing of 63.50 km 2 of forest area. Townships of Pekhon, Sesai, Ywangan and Kalaw have observed a striking decrease in forest area over the years. Township of Hopone and Yaksauk surprisingly shows an increasing trend of forest cover and could be an indicator of good forest practices at least at the regional level. Hence, expansion of non-forest areas which mainly consist of agricultural land, urban settlement and barren land have been encroaching fast into the forest areas. Construction of new infrastructure such as reservoirs, dams and road network is also important reason behind such high decline of forest cover in the region.
At national scale, the deforestation trends are alarming. The net loss of forest in Myanmar was reported to be 546 K ha year −1 between 2010 and 2015 which was 25% higher than in the 1990s (FAO 2015). Myanmar has undergone agricultural expansion for growing edible crops like rice, fruits, pulses, etc., and plantation of palm oil and rubber which has been responsible for about 1 million ha of forest loss between 2002 and 2014 (UN REDD 2017). The deforestation rate of Myanmar has been ranked second among the southeastasian countries after Indonesia (684 K ha year −1 ) and third in the world. In Indonesia, much of the forested land has been converted to palm oil plantation. Evidence shows that 63% of new palm oil in the country has come at the expense of loss of biodiversity-rich tropical forest between 1990 and 2010. On the other hand, several countries of South and Southeast Asia have shown increasing forest trends for the year 2010-2015. The annual rate of forest gain has been 3.3%, 0.9%, 0.8% and 0.3% in Philippines , Vietnam, China and India that are included in the top 10 counties with highest net forest gain (FAO 2015).

Spatial analysis of deforestation drivers
The WofE coefficients for each of the 8 variables were derived, where the positive coefficients relate to probability of transition potential favoring deforestation, while the negative coefficients indicate that the corresponding variables have potential to repel deforestation at those value ranges. Table 4 shows the positive coefficients of the 8 variables and their ranges in descending order. The highest WofE coefficient of drivers variable being distance from the newly constructed reservoirs and could be labeled as  (Table 4). It is evident from the derived coefficients that these geophysical and socioeconomics drivers of change do have significance on observed deforestation over the concerned years. If we look at each of these variables separately, it is easier to notice the WofE at their particular ranges. It is clear from Fig. 3 that different factors have varying amounts of effect on deforestation in either promoting or repelling it at different range intervals. The peaks (positive values) depict deforestation potential, whereas the lows (negative values) show ranges that defy deforestation for each variable. Evaluating the drivers of deforestation, it was found that some drivers like agricultural expansion, infrastructural development and fuelwood consumption, forest fires, etc., have direct effect on deforestation. Other drivers that relate to social, cultural or technologies factors like population growth, poverty and policy barriers have indirect effect on forest loss. In the case of elevation (Fig. 3a) factor, there has been a significant evidence of deforestation in the altitude of 100-200 m ranges with high coefficient of 2.25. This range is associated with lowland depicting relatively flatter areas in Taunggyi ideal for cultivation, drained by waterways and has been observed to be used for development of hydropower projects. The weight then decreases midway between 500 and 1000 m elevation with again a steady increase of up to coefficient of 0.78 at 2500 m elevation. This implies that an increase in deforestation has been observed at elevation greater than 1000 m and could be a result of forest fragmentation due to shifting cultivation, soil erosion and forest degradation affecting mostly higher altitude of Taunggyi. In Shan State, one of the major causes of deforestation is shifting cultivation or slash and burn farming, traditionally known as "taungya" which has been practiced by ethnic minorities or a long time in hilly areas (Myint 2018).
The slope factor (Fig. 3b) pointed out deforestation to be active at slope range less than 5 degree depicting most of the flat areas of the regions. This major portion of townships of Taunggyi, Sesai and Yauksauk form ideal agricultural tracts that are constantly cleared to expand the existing cultivable land. On the other hand, in the regions hilly terrain small landholdings and population growth force farmers to overexploit the natural resources cutting more trees for fuelwood and clearing land on steep slopes for cultivation (Myint 2018). The probability of deforestation is calculated to be highest (0.6) at 60° slope. Traditional practice of shifting cultivation is mainly operative in the steeper slopes of the region. This practice is now being replaced by opium cultivation and rubber cultivation and commercial growth of such crops are related to low levels of development and higher risk to environmental challenges (Myint 2018;UNODC 2017). Natural or human-induced landslides and forest fires are also important factors for forest loss on steep slopes. Non-forest area mostly consists of agricultural land, urban centers and other land not covered by forest. Regions in proximity to non-forest areas (Fig. 3c) have shown high probability of deforestation, which decreases as the distance increases. The highest weight of 0.95 is nearest to distance of 250 m to non-forest land. Hence, expansion of agriculture over the years could be a cause of deforestation prevalent up to distance of 500 m after which it then gradually repels deforestation. It is generally a notion to clear the surrounding forest fringe for expansion of its existing agricultural land to feed the growing population. Most of the forested land in flatter regions of Taunggyi district is lost because of its conversion to cropland, mainly rice and maize. Many parts of the hilly regions in country have recently undertaken large-scale rubber plantation. There has been a significant increase in area covered by rubber plantation in Shan State from 4000 to 74,200 ha between 2004 and 2005. The southern region of state (Taunggyi District) is characterized by shifting cultivation making forest cover highly dynamic. It is estimated that 15,000 ha of forest in the country is affected each year due to shifting cultivation. These have been an important factor in the increase in non-forested area in the region (UN REDD 2017).
Infrastructural development that includes the construction of reservoirs and artificial waterways, roads and built-up areas has been on rise in the recent years. Myanmar has the largest potential for hydropower development in Southeast Asia. Owing to this, many hydropower projects have been constructed in the country at the expense of forest and more so within reserved forest and protected forest (UN REDD 2017). It was reported that 110,777 m 3 of timber was cleared for hydropower development between 2011 and 2012 (Woods 2015). The factor proximity to the reservoirs (Fig. 3e) shows high coefficient values until 900 m of distance from these newly built infrastructures. The positive WofE coefficient in the vicinity of the reservoirs is because of construction of a reservoir in area previously occupied by forest. This factor can be clearly observed in the case of Taunggyi District where the construction of reservoir at certain location has disturbed existing forest. Other commercial land use projects and infrastructural development such as mining and transportation development through forested areas have adverse effect on the forest health and in majority of the times result in forest loss. Proximity to waterways (Fig. 3d) would mean clearance of forest in river banks and easy transportation of logged materials through navigable rivers. It also relates to canals and other water channels which are built connecting reservoirs/rivers to agricultural land mostly through the forested areas whose construction is related to deforestation. Through the model calibration, deforestation has been observed up to a distance of 950 m from the waterway quantified by the WofE 0.015. Similarly, construction of roads in heavy forested regions of Myanmar is proven as one of the direct drivers of deforestation. It is the result of infrastructural extension, resource extraction, agricultural expansion and providing accessibility to settlement. Construction of highways and major roads (Fig. 3f) is observed to be promoting deforestation up to a distance of 550 m up until 1 km. Analyzing the effect of the proximity to settlement (Fig. 3h) has also proved that deforestation is observed in areas closer to settlements. This is one of the direct drivers of deforestation arising from expansion of urban built-up and increased settlement in rural areas. Myanmar has undergone rapid growth in urban and peri-urban settlement since the last decade, and its built-up areas have increased by 24% from 1992 to 2010, surely at the expense of its forest. Urban settlement increases risk of deforestation of nearby forest for extension of built-up areas and agricultural land. On the other hand, the increase in rural settlement areas would mean higher fuelwood dependency to nearby forest, leading to greater deforestation and forest degradation rates (Sein et al. 2015).
Population density (Fig. 3g) evaluated at an equal interval of 20 persons/km 2 showed high degree of deforestation at density greater than 100 persons/km 2 population density. It got highest at population density of 240 person/km 2 showing higher population density areas evidencing higher deforestation event. In order to meet the demand of growing population, large amount of forest land is under stress. Population growth in most cities like Nay Pi Taw, Yangon and Taunggyi has increased pressure on existing forested land for expansion of infrastructure and meeting agricultural needs. A large portion of rural, peri-urban and urban population is dependent on forest for their daily fuel consumption. Almost 60% of the population in the country still rely on forest for fuelwood needs, exerting pressure on forest resources. The average annual consumption of fuelwood per household is estimated to be 2.5 cubic tons. With the population projected to increase from 53.9 million in 2015 to 60.2 million by 2030, consumption of fuelwood and charcoal is projected to increase to 55 million m 3 by 2030 leading to higher deforestation and forest degradation rates in the country (UN REDD 2017).
The simulation model produces yearly landscape and probability maps for the given number of years. The landscape changes every year with the transition probability calculated for each year. The forest areas with high probability value have a higher chance of being deforested in the landscape map of simultaneous iterative year. The model can be run n number of times to produce result of landscape for the nth year. Therefore, to produce the simulated landscape 2016 map, the model was iterated 8 times based on the probability of transition for those many years. The spatial independency of the input factor maps was tested using Cramer's (Bonham-Carter 1994). This test shows the degree of association among the factors examined being full or independent ranging from values 0 to 1, respectively. No significant spatial dependency was detected, and the Cramer's was less than 0.25 in all comparison cases.

Model validation
Validation of a land change model is usually carried out to determine the prediction ability of the model by comparing the predicted result to the reference map. Spatial models are generally validated in the neighborhood context because maps may not always match pixel by pixel, but they may present similar spatial patterns and spatial agreement within a certain cell vicinity (Thapa et al. 2013). We used a two-way reciprocal fuzzy similarity method (Almeida et al. 2008) in which the similarity between two maps is considered based on the influence of the cell and also, to some extent, by its neighbor cell. We compared the spatial similarity between the reference map 2016 and the simulated map 2016 at different scales. The exponential decay function was used to incorporate fuzzy similarity method to validate the prediction power of the model. Since it is a two-way reciprocal fuzzy similarity method, the function compares the differences between the model generated changes in landscapes and the reference (observed) changes in landscapes.
The similarity of the landscape change is 55% at pixel level and increases to 75.90% at 9*9 window (Fig. 4). The window size is further increased to reach a similarity fitness of 87.79% at 21*21 window size. The input land cover resolution was 30 m and the window size 21*21 pixels, so spatial resolution was 630 m.
It is seen that with the increasing window size, the similarity between the maps also increases meaning that the model can predict spatial patterns more accurately at coarser spatial resolution. The model users may select the model in different spatial resolution depending on their scale requirement (Thapa et al. 2013). For example, regional level planners would require much more detailed results than the province level planners and could therefore opt for finer resolution results. Since the model depicts fairly accurate level of compatibility between real and simulated landscapes, it has been further used to develop a simulation model to depict future forest cover scenarios.

Forest cover projection
It is also possible to project the forest cover scenarios to future using the historical trend as part of the simulation model. Figure 5 shows the trend of forest cover for the observed year 2016 and simulated years 2020, 2025 and 2030.
Under BAU scenario, the forest area decreased to striking 33.72% in 2030 as compared to 46.54% in 2016. This decline of 13.72% (3349.84 km 2 ) in forest cover is the result of continuous operation of driving factors and rate of deforestation. From the total forest loss that occurred between 2016 and 2030, high deforestation tends to affect southern townships of Taunggyi. A closer look into the areas affected by deforestation in each of the 5-year intervals 2020, 2025 and 2030 gives an idea of potential of forest loss in that period in a greater detail (Fig. 6).
In the township-wise analysis of future forest loss scenarios (Fig. 7), it was projected that 35% of forest loss has taken place in Pinlaung followed by Hopone (12%), Sesai (11%), Nyaungshwe (10%) and Pekhon (10%). Deforestation tends to affect the southern townships of Taunggyi since these are the areas where the drivers are most actively operating. The townships ranked highest to lowest in the case area of forest loss in the projection This could be because of higher proximity to previously deforested land (non-forest in first iteration) owing to highly fragmented landscape in the township. The driver proximity to previously deforested area seems to be operating with high weightage in Hopone Township with high probability of deforestation in the township. The northern most township of Yauksauk accounts for least amount of forest loss for the simulated years 2016-2030 since this area consists of dense compact forest with little disturbance from human activities like construction of reservoirs and roads, fewer settlement and population density of only 35 persons/km 2 . High percentage of deforestation in the district is sure to have longterm effect on biodiversity and regional climate regime. Therefore, it is very important for immediate public and private intervention to conserve the forest and to discontinue the present trend of deforestation.
According to the Nationally Determined Contributions, 2016 (NDC) document of Myanmar, the country formulated several actions relevant to climate change mitigation taking 2030 as the target year. Forest preservation measures and the action area of "environment and natural resources" including REDD+ has been in the priority list of the government (MONREC 2018). Thus, keeping in mind the changes in the forest area and potential drivers causing the changes becomes extremely relevant in determining where to lay focus in the coming years while deciding priority zones and framing action plans.
This study mainly uses the available spatial variables as geophysical drivers of deforestation, but many researches, for example Kolb et al. (2013), used a large number of environmental, biophysical and social variables including landforms, hydrology, ethnicity, migration and socioeconomic status. This requires more number of data from field study, but as Pontius et al. (2008) pointed out that a complex model does not necessarily lead to higher predictive power. Generalized models that are flexible with provision for including increasing complexity are often preferred in planning process. Also, land change models can only provide projections based upon the quality of inputs provided to the model. There may be over-or underestimation of future land change due to fluctuations in the annual land change rates, afforestation programmes or changes in government policy that may impact deforestation rates (Elz et al. 2015;Angelsen and Wertz-Kanounnikoff 2008). In addition, there exist knowledge gaps and absence of data in several areas of work in Myanmar. The extent to which agricultural expansion is happening in non-forest areas or leads to the clearance of fairly intact natural forests is unknown (UN REDD 2017). Hence, complete understanding and inclusion of other socioeconomic drivers, geological factors, timber trading policies and forest management practices is important for the development of a more accurate model. Filling in these gaps would significantly help to more confidently develop measures and policies for a national REDD+ strategy.

Conclusion
The study used readily available JAXA FNF map to analyze the changes in forest cover and presented a simple and easily replicable model to study future changes for Taunggyi District. It was derived that the total forest cover declined from 53.84% in 2008 to 46.54% in 2016 and projected to 32.72% in 2030 in the study area. We quantified factors of deforestation and learnt that construction of new infrastructure like reservoirs, certain elevation levels, slope, population density and distance to previously deforested areas received the highest weights. Spatial probability map of deforestation was also simulated that helped demarcate areas with the highest probability to deforestation. A decline of 13.82% of forest area is simulated for the years 2016-2030.
Since this paper provides a closer look at the drivers and is one of the unique studies in Myanmar, the amount and quality of geographic data are very limited. This study can provide basis for other researches in the country. It could be of importance to policy makers while formation of baseline study, monitoring of forest and formulation of conservation strategies. As initiative like REDD+ gain momentum, this research can provide for better understanding of current forest loss, rate of deforestation and identification of critical areas to be brought under its effect. Since Myanmar has significant potential to reduce its forest carbon emissions and sustainably manage its forest carbon stocks, it is taking rapid steps in implementing REDD+ readiness phase. Any REDD+ Demonstration Activity in this district will need to enhance the protection of existing national parks, implementation of existing regulations and alternative land use opportunities for locals while addressing key drivers of deforestation.