A Multi Sensor Approach to Forest Type Mapping for Advancing Monitoring of Sustainable Development Goals (SDG) in Myanmar

Monitoring forests is important for measuring overall success of the 2030 Agenda because forests play an essential role in meeting many Sustainable Development Goals (SDG), especially SDG 15. Our study evaluates the contribution of three satellite data sources (Landsat-8, Sentinel-2 and Sentinel-1) for mapping diverse forest types in Myanmar. This assessment is especially important because Myanmar is currently revising its classification system for forests and it is critical that these new forest types can be accurately mapped and monitored over time using satellite imagery. Our results show that using a combination of Sentinel-1 and Sentinel-2 yields the highest accuracy (89.6% ± 0.16 percentage point(pp)), followed by Sentinel-2 alone (87.97% ± 0.11pp) and Landsat-8 (82.68% ± 0.13pp). The higher spatial resolution of Sentinel-2 Blue, Green, Red, Narrow Near Infrared and Short Wave Infrared bands enhances accuracy by 4.83% compared to Landsat-8. The addition of the Sentinel-2 Near Infrared and three Vegetation Red Edge bands further improve accuracy by 0.46% compared to using only Sentinel-2 Blue, Green, Red, Narrow Near Infrared and Short Wave Infrared bands. Adding the radar information from Sentinel-1 further increases the accuracy by 1.63%. We were able to map the two major forest types, Upper Moist and Upper Dry Mixed Deciduous Forest, which comprise 90% of our study area. Accuracies for these forest types ranged from 77 to 96% depending on the sensors used, demonstrating the feasibility of using satellite data to map forest categories from a newly revised classification system. Our results advance the ongoing development of the National Forest Monitoring System (NFMS) by the Myanmar Forest Department and United Nations-Food and Agriculture Organization (UN-FAO) and facilitates future monitoring of progress towards the SDGs.


Introduction
Forests play an essential role in addressing the Sustainable Development Goals (SDG) set by the United Nations in the 2030 Agenda, by providing goods and services necessary to meet the SDGs [1,2]. Specifically, SDG 15 (Life on Land) explicitly mentions forests and sustainable forest management in its targets 15.1 and 15.2, respectively. Afforestation and halting deforestation are recognized as key pathways for countries to reach their mitigation targets, because forests can serve as major carbon sinks (SDG 13). Given the significant contribution of forests to the SDGs, monitoring of forests is critical to measuring overall success of the 2030 Agenda.
In the context of SDG 15, accurate information on a country's forest area, as implied by indicator 15.1.1, is dependent on many factors including definitions of national classifications, data source and methodology [3]. Ensuring high quality data for monitoring and reporting, specifically accurate and up-to-date forest type baseline data, remains a logistic, financial and technical challenge in many countries, including Myanmar [4].
Myanmar lies in the Indo-Burma biodiversity hotspot [5] and has the highest forest cover in mainland Southeast Asia [6]. In 2011, the country started the process of political transition from a military government to a democracy accompanied by the opening of its economy. This led to a surge in foreign investments and resulted in increased deforestation. Myanmar had the third highest rate of deforestation globally between 2010 and 2015 [7]. In response, the Government is strengthening the existing National Forestry Master Plan which aims at demarcating 30% of the country's area for permanent forest cover and 10% of country's total area for protected areas [7] and has established ambitious targets for reducing emissions from land use change and forestry under the 2015 Paris Agreement and the Nationally Determined Contributions (NDCs), presently in the process of updating. However, different forest types differ in their conservation value, representation within protected areas and vulnerability to loss and degradation. A primary challenge for the management of Myanmar's unique forests is the absence of recent and accurate forest type maps. Furthermore, forest type definitions differ by national and international agencies and do not accurately represent the forest types at subnational levels. To address this issue, the Myanmar Forest Department is working with the United Nations-Food and Agriculture Organization (UN-FAO) to harmonize and modernize the forest type classifications under the United Nations Programme on Reducing Emissions from Deforestation and Forest Degradation (UN-REDD Programme) supported enhancement of the National Forest Monitoring System (NFMS). An important consideration while developing the revised forest type classification is the feasibility of using current satellite data sources to map the revised forest types accurately.
Studies evaluating Landsat-8 and Sentinel-2 imagery to map land-cover land-use change [26,32], forests attributes such as forest successional stages [33], canopy cover and leaf area index [27] and other forest variables such as stem volume, stem diameter, tree height and basal area [25], have found Sentinel-2′s performance comparable to or better than Landsat-8 depending on the band combination used in the evaluation process. Classification of land use and land cover in Burkina Faso [19] using all Sentinel-2 bands produced an overall accuracy that is 5% higher than using Landsat-8 bands only, while classification of land use and land cover using Sentinel-2 bands, similar to Landsat-8, improved accuracy by 4%. Sentinel-2 outperformed Landsat-8 in forest variable prediction on account of its Vegetation Red Edge band and better spatial resolution. In their study, [25] found that the average of the absolute value of the relative systematic error was 0.8% for Sentinel-2 and 1.2% for Landsat-8. However, research in boreal forests shows that the performance of Sentinel-2 was comparable or slightly better than Landsat-8 in discrimination of different successional forest stages (kappa =0.98 vs 0.9) [25] and estimation of canopy cover (root mean square error = 24% vs 24.5%) and leaf area index (root mean square error =19.6% vs 20.2%) [20]. Though studies have been conducted on the synergistic use of Landsat-8, Sentinel-2 and Sentinel-1 on the topics of land cover change, successional stages of tropical forests and forestry related variables, only one study has evaluated the potential of Sentinel-1 and -2 to classify detailed tropical forest types [34].
Random Forest is a popular machine learning algorithm for classifying large multiclass datasets as it effectively handles multi-collinearity and noise in multidimensional, non-parametric datasets, without overfitting. It works well even with datasets that have a large number of predictive variables [35,36]. Random Forest algorithm has been widely used in applications ranging from urban mapping [37,38] to forestry [35], exploiting information from spectral reflectance, texture and spatial metrics [39] derived from remotely-sensed datasets. Other methods of statistical analysis, such as discriminant analysis and logistic regression, have been used successfully to differentiate tillage classes in Nebraska from Landsat-5 TM [15]. A comparison of six base models of supervised classification, which included linear and quadratic discriminant analysis, logistic regression, classification trees, neural networks and support vector machines, found that Random Forest performed best in the prediction of Oak habitat [40]. Hence, we decided to use a Random Forest algorithm for our study. The availability of free moderate resolution satellite imagery and machine learning algorithms in open source cloud computing platforms provides an opportunity to analyze large satellite datasets quickly and with minimum computing infrastructure [41]. These advances now significantly aid broad-scale mapping for forest monitoring and management applications at subnational, national or continental scales.
Past satellite-based forest mapping studies in Myanmar have focused on assessing changes in forest cover at the national level [6,42] and in forest type at regional level [43] using Landsat images, but the utility of Sentinel-1 and -2 has not been investigated for forest type mapping in Myanmar. We used a Random Forest algorithm implemented in Google Earth Engine (GEE), a web-based cloud computing platform for geospatial analysis [41], to classify forest types from freely available, multitemporal Landsat-8 and Sentinel-1 and -2 images for the year 2018.
The overall goal of this study is to compare the contribution of multi-temporal images from Landsat-8 and Sentinel-1 and -2 for developing forest type maps focusing on a pilot study area, the Bago Region in Myanmar. This pilot area was selected because its landscape is dominated by Mixed Deciduous Forests which is assumed to represent the largest forest type by area (39%) in Myanmar [44]. Henceforth, Mixed Deciduous Forest is abbreviated to MDF in text, tables and figures.
Our objectives are: i) To provide a simple, open-source method to map forest types using satellite data sources that can be easily adopted by other countries. ii) Test the accuracy of mapping the revised forest classes using recent active and passive sensors-Sentinel-1, Sentinel-2 and Landsat-8. iii) Evaluate the contribution of each sensor and its metrics to successfully map the revised forest types. We hypothesize that: 1.
The accuracy of the Sentinel-2-based forest type classifications is higher compared to 2.
The accuracy of the combined Sentinel-1 and -2 based forest type classification is higher than the Sentinel-2 based forest type classification because the radar data adds more information on vegetation structure and moisture.

Study Area
The Bago Region is in the southern part of Central Myanmar between 46°45'N and 19°20'N and 94°35'E and 97°10'E ( Figure 1). The Bago Yoma mountain range (also known as Pegu Yoma or Pegu Range), runs in a north-south direction and influences the vegetation patterns in the landscape. The elevation ranges from 0 m to 1909 m with mean elevation of 146 m according to the analysis of elevation data from the Shuttle Radar Topographic Mission (SRTM) [45]. According to Geographic Information System (GIS) data shared by the Myanmar Forest Department, there are six major types of forests in Bago, including Bamboo, Indaing Forests, Lower MDF, Upper Dry MDF, Upper Moist MDF, and Swamps. Teak (Tectona grandis) is a commonly occurring species in the MDFs in this region and is commercially logged for its high-quality wood. Extensive logging of teak has resulted in degradation and soil erosion [46]. The Indaing Forests are Deciduous Dipterocarps forests, dominated mainly by Dipterocarpus tuberculatus.

Earth Observation Datasets and Processing
We used GEE [41] to select, process and classify dry season images (November 2017-April 2018) from Sentinel-1, Sentinel-2 and Landsat-8. Thus, our forest type map reflects the spatial distribution of different forests in 2018.

a) Landsat-8 and Sentinel-2 processing
To select images for analysis, we applied the image quality filter for all Landsat-8 images [47] and a cloud filter (< 50% cloud cover) to all Landsat-8 and Sentinel-2 'Top Of Atmosphere (TOA)' reflectance products. TOA reflectance products are corrected for differences in sun elevation and sunearth distance to account for spectral variability resulting from differences in acquisition date and time [48]. The image quality of the TOA reflectance products can be further improved by applying atmospheric corrections to yield Surface Reflectance products [49]. However, Surface Reflectance products for Sentinel-2 that cover our study area were not available for the year 2017 in GEE. Hence, to ensure that all our analyses can be transferrable via the GEE cloud computing facility and to maintain consistency between contrasting remote sensing models, we used TOA data products for both Landsat-8 and Sentinel-2. This is a common approach and has been widely implemented in land cover change products [50][51][52]. Cloud masks were applied to all Landsat-8 and Sentinel-2 TOA images based on information available in the Quality Assessment (QA) bands.
All the Landsat-8 bands used in our study have 30 m resolution. Sentinel-2 has four 10 m bands (R, B, G and NIR) and six 20 m bands (VRE1, VRE2, VRE3, NNIR, SWIR1, SWIR2). All the 10 m Sentinel-2 bands were resampled on-the-fly to 20 m using nearest neighbor in GEE to match the coarser resolution data (see https://developers.google.com/earth-engine/scale for on-the-fly resampling).

b) Sentinel-1 processing
We used Sentinel-1 Ground Range Detected (GRD) images available in GEE. These images were already calibrated and orthorectified using Sentinel-1 toolbox and have 10 m spatial resolution. For this analysis we chose images collected in Interferometric Wide (IW) swath mode. IW swath mode is the main acquisition mode over land for Sentinel-1. The IW mode collects data in dual polarization. It means that the microwave sensor transmits either horizontally polarized (H) or vertically polarized (V) waves. Depending on the transmission polarization, the received backscatter wave may be of one of the four types: HH (horizontal transmit and horizontal receive), HV (horizontal transmit and vertical receive), VV (vertical transmit and vertical receive) and VH (vertical transmit and horizontal receive) [53]. Details on data acquisition and processing algorithms for Sentinel-1 may be found in the Sentinel-1 User Guide Introduction (https://sentinel.esa.int/web/sentinel/user-guides/sentinel-1sar) and Sentinel-1 Technical Guide (https://sentinel.esa.int/web/sentinel/technical-guides/sentinel-1sar). Previous studies have found Sentinel-1 VV and VH polarizations to be useful in forest mapping applications [21,34]; hence, we used both the VV and VH backscatter in this study. As in the previous analysis, we resampled 10 m bands to 20 m.

c) SRTM
To include a topographic variable in our analysis we used 30 m resolution SRTM elevation data available in GEE.

Variables used in classification
To study the phenological characteristics of different forest types, we created six cloudless, monthly median composites (November(Nov), December(Dec), January(Jan), February(Feb), March(Mar) and April(Apr)) from the available images for Landsat-8 and Sentinel-2. To achieve this, we stacked images based on month and tile and then computed the monthly median value per pixel. We also calculated 18 texture metrics (Table 1, [54,55]) for each band and month for each sensor. All texture metrics were derived using a 3 × 3 window. Texture metrics [54,55] have been widely used to improve forest classifications [56,57] in complex landscapes [58] and tropical forests [59]. We also computed band ratios for all band combinations for each sensor. Finally, we include monthly Normalized Difference of all bands for Landsat-8 and Sentinel-2. Table 1. List of texture metrics and abbreviations used as variables in supervised forest classification from satellite imagery for the Bago Region of Myanmar based on [54,55].

No.
Name of Texture Metric Abbreviation 1.
Sum of Variance svar 8.
Sum of Entropy sent 9.
Cluster Prominence prom

Masks
To restrict our analysis to forest pixels, we applied three masks: a tree cover mask based on the percent tree cover threshold, a forest loss mask and Normalized Difference Vegetation Index (NDVI) mask.
a) The tree cover mask was created using the Treecover2000 layer from Global Forest Change (GFC) 2018 dataset [17]. The Treecover2000 layer from GFC 2018 dataset provides a continuous percent tree cover value for all land pixels for the year 2000. The pixel values of percent tree cover layer range from 0 to 100%. Depending on the biome, Hansen et al., 2010 recommended using a threshold of 25-30% tree cover to identify woody vegetation taller than 5 m as forest. We decided to adopt the threshold of 25% considering the broad range of variation of forest types in Myanmar ranging from dry forests to evergreen forests. Compared to the 10% threshold used by FAO, the 25% threshold is more conservative but is better suited for global scale applications of mid resolution satellite image [60]. Areas with tree cover less than 25% were excluded from analysis. c) Since the GFC 2018 product has a ± 1 year variation in the date of image used, and Bago region has high rate of deforestation [46], an additional NDVI mask was developed from the most recent Sentinel-2 image between November to April to ensure that the sampled pixels have vegetation. The NDVI value was computed as: where NIR is Sentinel-2, band 8A (Narrow Near Infrared band), and R is Sentinel-2, band 4 (Red band).
We selected pixels with NDVI values greater than 0.34 to further ensure that the samples have vegetation based on the recommendation available at https://custom-scripts.sentinelhub.com/custom-scripts/sentinel-2/ndvi/" and preliminary analysis of the value that had the best results classifying forests. A threshold of 0.34 was determined based on visual inspection of the landscape.

Training and Validation Datasets
The Myanmar Forest Department provided digital data from District Forest Inventories(DFI) for Bago Region. This data was collected by Department staff during field survey between 2010 and 2017 and it consists of 3180 georeferenced points with forest type information, spaced at 2 km intervals on a rectangular grid.
We overlaid the 3180 points on forest pixels identified earlier. Only points that overlay on forested pixels were selected for further analysis. For our training, we used randomly selected 984 data points that fell inside our forest mask and that represented six major forest types: Bamboo, Indaing The validation dataset was created by combining the remaining points of DFI and additional polygons groundtruthed using high resolution image on Google Earth. Additional validation polygons were created for Bamboo, Indaing, Lower MDF and Swamps that have a low number of points within DFI. Validation data for widespread forest types are randomly selected across the scene. However, because of the small coverage of a few forest types in the region, the validation data of those forest types are strategically searched and located around existing DFI data point. We searched high resolution images in Google Earth for rarer forest types around the sampled DFI points and matched them with characteristic spectral properties to identify validation data for rarer forest types. The process is a two-stage cluster sampling process, which is a valid and practical option when search cost and travel cost are prohibiting and when the extent of the study area is extensive [61]. The proportion of validation samples per classes in the validation data is: Swamp= 0.54% (40 pixels), Bamboo = 0.22% (16 pixels), Lower MDF = 0.50% (37 pixels), Upper Moist MDF = 28.9% (2125 pixels), Upper Dry MDF = 69.55% (5114 pixels) and Indaing Forest = 0.27% (20 pixels). The mixed approach of using existing data were possible and augmenting this data for validation, where needed, resulted in a split of 85% to 15% training to validation, respectively [62].

Analytical Approach and Flow
Our approach combines the time-series of moderate resolution passive (Sentinel-2, Landsat-8) and active (Sentinel-1) satellite datasets to extract spectral and structural information for each forest type and their response to seasonal variations ( Figure 2). To identify the important variables for our classification, we extracted the spectral and structural values for each training polygon across all bands in GEE and exported the resulting training dataset to R [63] and ran a Random Forest classifier. The 'importance' function in randomForest package [64] in R was used to identify the important variables. In general, increasing the number of variables increases model accuracy. However, since Random Forest models avoid overfitting, after a certain threshold, adding more variables does not necessarily improve model accuracy. This threshold can be determined by plotting a plot of Mean Decrease in Accuracy vs the model variables using the 'importance' function, run in R. For all our models we determined this threshold using the Mean Decrease in Accuracy vs the model variables plot to ensure that we choose correct number of variables for the models. A plot of the variables introduced vs the mean decrease in accuracy revealed a clear threshold, beyond which adding new variables did not improve classification accuracy much. The threshold values for the Landsat-8, Sentinel-2 all bands and Sentinel-1 and -2 models was 6, 5.6 and 12. Based on this threshold, we only employed important variables to run our final classification models in GEE. For example, for the Landsat-8 model, a plot of Mean Decrease in Accuracy vs the model variables revealed that beyond a threshold of 6, adding more model variables did not improve the accuracy much. To compare classification performance between sensors and to evaluate the contribution of different bands to the overall accuracy of the classification, we created four classification models ( Table 2)  During classification in Random Forest, we set the number of trees to 500 and the number of variables/split to the square root of the total number of variables in the model. We also used a 'seed' setting equal to 1 for each model rerun so that models could be reproduced reliably. To assess and compare the outputs and performance of our four different models (Landsat-8 Model, Sentinel-2 with Landsat-8 like bands Model, Sentinel-2 all bands Model and Sentinel-1 and -2 Model), we calculated the absolute and percent forest area for each forest type as well as area proportion weighed area estimates, user's accuracy and overall accuracy with 95% Confidence Intervals [24,65].

Areal Estimates
The four (Landsat-8 Model, Sentinel-2 with Landsat-8 like bands Model, Sentinel-2 all bands Model and Sentinel-1 and -2 Model) forest type classifications models produced consistent areal estimates with some minor variation among the forest types (  Table 3).
The swamp has the highest variability in the area, ranging from 2.3% (Landsat-8 model) to 9.05% (Sentinel-1 and -2 model) in the order Sentinel-2-all bands < Landsat-8 < Sentinel-2 with Landsat-8 like bands < Sentinel-1 and -2. The Sentinel-1 and -2 model predicts more Swamp area, almost by a factor of 4 ( Figure 3), when compared to the Landsat-8 model. The predicted Landsat-8 Swamp area estimate is 1.5 times smaller when compared to the predicted Swamp area by the Sentinel-2 with Landsat-8 like bands model, though they have same input variables derived from comparable bands. However, the area of Swamp predicted by the Landsat-8 model is 1.5 times more than the Sentinel-2 all band model.

Input predictors and Model Accuracy
Important predictors in the Landsat-8 model included bands B, G, R and SWIR1 ( Table 4). The Landsat-8 model had the lowest accuracy of all the models ( Table 5). The Sentinel-2 with Landsat-8 like bands model had the same input variables as the Landsat-8 model; however, it showed improved accuracy of 4.83%. More than 50% of the important variables in Sentinel-2 all band model were derived from the VRE bands and the Sentinel-2 NIR band ( Table 4). The addition of the VRE derived metrics improved the accuracy of the Sentinel-2 all band classification by 0.46% compared to the Sentinel-2 with Landsat-8 like bands model. Overall, the Sentinel-2 all bands model had an accuracy of 87.97% ± 0.11 pp, an improvement of 5.3% from the Landsat-8 model.
The savg texture variable appeared most often on the list of important variables (Table 4) in our classification models. Adding monthly composites and textures of VV, VH and the VV/VH bands (Table 5) in the Sentinel-1 and -2 model increased the accuracy by 1.63% and 6.92% compared to the Sentinel-2 all bands model and Landsat 8 model, respectively (Table 5).
The user's accuracy for most common and most widely distributed forest types Upper Moist MDF and Upper Dry MDF, together comprising 90% of our study area, varying from 77 to 96% depending on the data source (Table 6). Table 5. Accuracy of each model in separating forest types in the Bago Yoma of Myanmar.

Discussion
In this study we mapped the different forest types in the Bago Region, using satellite data and a revised forest classification provided by the Myanmar Forest Department. According to our maps, the Upper Dry MDF is the largest forest type by area, with its extent ranging between 7322.96 and 9040.27 km 2 , followed by Upper Moist MDF (5849.77-6551.12 km 2 ). A combination of both the forest types cover approximately 90% of the total forest area in Bago Region. Indaing is the least dominant forest type by area in the region.
Our four models (Landsat-8 Model, Sentinel-2 with Landsat-8 like bands Model, Sentinel-2 all bands Model and Sentinel-1 and -2 Model) showed overall agreement in total area and spatial extent (Table 3 and Figure 4). The disagreement between the four models on the extent of Swamp in the northeastern slopes of the Bago Yoma was investigated in detail due to availability of ground truth data from DFI. The DFI classifies this area as Upper Dry MDF, which is correctly predicted by the Landsat-8 model. Inspection of the area in Google Earth imagery shows it as a low-laying low forest cover region with presence of soil moisture evident from the rivers and streams running into the several lakes in this area. The variation in the model predictions may be explained by the sensitivity of the input variables in each model to soil moisture. The VV band of Sentinel-1 is highly sensitive to soil moisture and the sensitivity decreases with the increase in vegetation cover [66]. Hence, it is no surprise that the Sentinel-1 and -2 model with VV input variable classifies it as Swamp. The Sentinel-2 NIR band (Band 8) in Sentinel-2-all band model includes water absorption wavelengths that were intentionally avoided by the Landsat-8 NIR band design [67] hence, Landsat-8 classifies the same region correctly as Upper Dry MDF, while all the Sentinel-2 models show presence of varying degrees of Swamp pixels. Landsat-8 is known to underestimate the water areas, and often misclassified water with vegetation class in moist conditions [68]. This explains the low presence of Swamp pixels in Landsat Accuracy of the forest type classifications varied by sensor and forest type. Absence of overlap between the area weighted overall accuracy and associated 95% confidence intervals show that the differences between the four models are statistically significant (Table 5). Ranked in order of increasing overall accuracy, Landsat-8 (82.68% ± 0.13pp) < Sentinel-2 bands comparable to Landsat-8 bands (87.51% ± 0.12pp) < Sentinel-2 all bands (87.97% ± 0.11pp) < Sentinel-1 and -2 bands (89.6% ± 0.16pp), proving both our hypotheses to be true. Our first hypothesis holds true as the Landsat-8 based classification showed least accuracy (82.68% ± 0.13pp) compared to both Sentinel-2 based models (87.51% ± 0.12pp, 87.97% ± 0.11pp). Our second hypothesis also holds true as the Sentinel-1 and -2 based classification showed higher accuracy (89.6% ± 0.16pp) compared to the Sentinel-2 all band model (87.97% ± 0.11pp). These results can be further explained by considering the role played by important predictors. Important predictors in Landsat-8 model, included the B, G and SWIR bands. The Blue band is affected by atmosphere, and thus, widely known for its application in coastal and aerosol studies. However, our results are in agreement with previous studies [25,69,70] that found the Blue band to be useful in mapping broadleaf deciduous species. In leaf-off conditions the Blue and Green band of Landsat-8 is more sensitive to stand age in deciduous species [70]. The importance of Blue band can be explained based on the high absorption of blue light (~0.46 µm Landsat-8 Band 2) by chlorophyll b, but not by chlorophylls a, c or d. During the dry season, deciduous plants stop producing chlorophyll in preparation of dry/leaf-off season and, as a consequence, have lower absorption and higher reflectance in blue bands during the dry season when compared to non-deciduous plants. Landsat-8 B band was found to be among the top predictors of tree height and volume in boreal forests [25]. Combination of high-resolution drone imagery B band along with R and G has been shown to separate five deciduous species [69]. The SWIR1 band is sensitive to vegetation moisture [71]. Thus, it is expected that these bands will play an important role in identifying the moist and dry deciduous forests from other forest types.
The improvement of 4.83% accuracy of Sentinel-2 with Landsat-8 like bands model (87.51% ± 0.12pp) when compared to Landsat-8 model (82.68% ± 0.13pp) may be attributed to the difference in spatial resolution of Landsat-8 (30 m) and Sentinel-2 (20 m). Sentinel-2, with its higher spatial resolution, results in less mixed class pixels in the training and validation data. Efforts to discriminate grasses grown under different management practices exhibited higher accuracy for simulated Sentinel-2 MSI compared to Landsat-8 OLI (82% vs 75%) [72]. Additionally, small differences in the wavelengths covered by each band [73] may have affected accuracy.
Comparison of the classification by the Sentinel-2-all bands model to Sentinel-2 with Landsat-8 like bands model highlight the contribution of the VRE and NIR bands to forest type mapping. Adding the different VRE and NIR band metrics increases the accuracy by 0.46% to 87.97% ± 0.11 pp. Similar increases in accuracy have been reported by other studies when using VRE bands in a Random Forest classifier to map land cover [20,21,56].
Addition of the VV and VH bands in Sentinel-1 and -2 model improve the accuracy by 1.63% and show the contribution of the C band radar to forest type mapping over the optical imagery. The radar data complement the optical data and add information on structural components and phenology of the vegetation [74,75]. Both the VV and VH band are sensitive to stem volume [76,77] and were expected to help in separating hardwood forests from Bamboo forest. The Sentinel-1 and -2 model overclassifies the Swamp and Indaing Forests in the landscape. This might be due to Sentinel-1′s sensitivity to soil moisture [66] and limited ability of VH band to discriminate between forest types [34]. Sentinel-1 has been found to have limited improvements over Sentinel-2 in homogeneous forested landscapes [21]; hence, its limited contribution to forest type mapping in Bago, 90% of which is dominated by MDFs, is not unexpected.
Among the texture metrics, the most frequently occurring metric in the important variable list was savg. "savg" is related to tree volume and basal area [78] and its importance can be explained by the dominant presence of big hardwood deciduous trees in the MDFs of Bago Yoma region.
Overall, our analysis shows that a higher resolution, multi-sensor approach using Sentinel-1 and -2 achieves better overall accuracy in mapping forest types in Bago compared to using medium resolution Landsat-8. All our Sentinel-2 based models performed better than Landsat-8 models on account of higher resolution (20 m Sentinel-2 vs 30 m Landsat-8), added spectral bands (three VRE and one NIR in Sentinel-2) and radar (Sentinel-1) bands. Though the accuracy of the Sentinel based classifications is higher than Landsat-8, Sentinel-2 lacks the temporal legacy of Landsat. Landsat has provided satellite imagery for land cover analysis since 1972, while Sentinel-2 was launched in 2015. The Harmonized Landsat Sentinel-2 (HLS) dataset [79] combines Landsat and Sentinel-2 images to provide dense time series images, which provide improved information on vegetation phenology. The use of HLS datasets may be explored in future studies for forest type mapping. The HLS dataset is corrected for registration errors between Landsat-8 and Sentinel-2. In our study, we used Landsat-8 and Sentinel-2 images separately in the models and did not combine them into a harmonized time-series; thus, we did not co-register the Landsat-8 and Sentinel-2 images. Future studies may also consider using Sentinel-2 Surface Reflectance products, when it is available on GEE.
Our forest type map provides a transparent, unbiased and robust estimate of the distribution of different forest types in Bago Region. Identifying forest type distribution would enable the decision makers to have accurate inventory information on the unique forest resources. Forests play an important role in Myanmar's economy by providing direct and indirect benefits to many important sectors, such as forestry and agriculture. The combined value of the services that forests provide to the Myanmar economy has been estimated at US$ 7.3 billion each year [80]. Different forest types provide different types of values, goods and services and are exposed to different types of threats. For example, hardwood trees found in mature MDFs are valuable to Myanmar's economy. However, degradation of the MDFs due to overexploitation of the hardwood species it supports, is a serious environmental issue that needs to be addressed [46]. A distribution map of the MDF forest type is necessary to manage the remaining MDFs sustainably. Similarly, knowing distribution of the bamboo forest type would enable decision makers to formulate poverty alleviation policies, as bamboo and rattan are widely used as building material by people in rural areas. One-fourth of the country's total population falls under poverty and 85% of the total poor live in rural areas [81] close to forests, and protected areas depend on these resources for survival [82]. The Indaing forests in Central Dry Zone also face intensive anthropogenic pressure for conversion to agriculture [83]. Identification of the current distribution of the threatened Indaing forests from our maps can help the Government to develop policies for forest conservation. Sustainable forest management is, therefore, crucial for Myanmar to achieve the SDGs.
The recent socio-political reform provided the new Government an opportunity to draft new regulations to aid forest conservation. Forest management is recognized as a priority activity, implemented through three main policy and institutional arrangements-National Forestry Master Plan, National Biodiversity Strategies and Action Plan (NBSAP), the REDD+ readiness roadmap, the National REDD+ Strategy and Action Plan and the recently updated NDC for the land use change and forestry sector. In August 2018, the Government of Myanmar adopted the Myanmar Sustainable Development Plan (MSDP), which provides a blueprint of how Myanmar as a nation seeks to achieve sustainable development and contribute to the global SDGs. The MSDP is further divided into three pillars and five goals. Goal 5 addresses the role of forests in achieving sustainability. The absence of a recent forest type map of Myanmar makes it difficult for the Government to manage forest sustainability, and our study fills the gap by demonstrating a forest type mapping methodology for Myanmar that can be scaled to nationwide application. Our study contributes to the Government's efforts by addressing forest monitoring at a subnational level as most of the interventions contributing to these objectives operate at a local level. The forest type map developed provides a quantifiable baseline against which future changes in forest extent can be evaluated. Information on forest type distribution and quantification of changes in extent of different forest types are important for detailed assessment of carbon storage and international climate reporting requirements. Quantifying change in land cover and identifying the drivers of change also advance our understanding of land cover science.

Conclusions
This study demonstrates the utility of recent optical (Landsat-8, Sentinel-2) and radar (Sentinel-1) satellite data to map forest types in Myanmar using a revised forest type classification. Our analysis shows that using Sentinel-2 alone or a combination of Sentinel-1 and Sentinel-2 imagery greatly improves map accuracy (88-90%) compared to only using Landsat-8 (82.68% ± 0.13 pp) imagery. The higher spatial resolution and the characteristic wavelengths of Sentinel-2 improve the accuracy by 4.83% compared to Landsat-8. Addition of the VRE and NIR bands from Sentinel-2 and VV and VH bands from Sentinel-1 further increase the map accuracy by 0.46% and 1.63%, respectively. We were able to map the two major forest types, Upper Moist and Upper Dry MDF, which comprise 90% of our study area, with an accuracy ranging from 77 to 96% for the different sensor specific models, showing the feasibility of using recent satellite data to map the revised forest classifications.
Our approach used freely-available satellite data from the National Aeronautics and Space Administration (NASA) [84][85][86] and the European Space Agency (ESA) [87,88] and all the computation was done using free statistical software R and cloud computing environment of GEE [41], which makes our approach easily transferable across regions. This is especially useful for countries such as Myanmar, which lack recent accurate forest type map and have limited financial and technical resources to meet their international reporting obligations. Our efforts are expected to contribute to the ongoing National Forest Management System (NFMS) conducted by the Myanmar Forest Department and UN-FAO, and contribute to several SDGs, specifically SDG 15, target 15.1.