Projected Impacts of Climate Change on the Protected Areas of Myanmar

: Protected areas are the backbone of biodiversity conservation but are ﬁxed in space and vulnerable to anthropogenic climate change. Myanmar is exceptionally rich in biodiversity but has a small protected area system. This study aimed to assess the potential vulnerability of this system to climate change. In the absence of good biodiversity data, we used a spatial modeling approach based on a statistically derived bioclimatic stratiﬁcation (the Global Environmental Stratiﬁcation, GEnS) to understand the spatial implications of projected climate change for Myanmar’s protected area system by 2050 and 2070. Nine bioclimatic zones and 41 strata were recognized in Myanmar, but their representation in the protected area system varied greatly, with the driest zones especially underrepresented. Under climate change, most zones will shift upslope, with some protected areas projected to change entirely to a new bioclimate. Potential impacts on biodiversity include mountaintop extinctions of species endemic to isolated peaks, loss of climate specialists from small protected areas and those with little elevational range, and woody encroachment into savannas and open forests as a result of both climate change and rising atmospheric CO 2 . Myanmar needs larger, better connected, and more representative protected areas, but political, social, and economic problems make this di ﬃ cult.


Climate Change and Protected Areas
Tropical East Asia supports 15-25% of global terrestrial biodiversity in around 4% of Earth's total land area [1], but is also one of the most threatened regions of the planet [2]. Major threats include the world's highest deforestation rates [3], widespread habitat degradation [4], hunting [5], wildlife consumption and trade [6], and unsustainable land use [7]. In the near future, climate change is likely to become an additional major driver of habitat degradation and species loss [8]. Warming in the tropics has been large relative to natural climate variability, so most regions are already experiencing thermal climates that were unknown in the 19th century [9]. In contrast, changes in precipitation in tropical East Asia have generally been small relative to natural variability. In other parts of the world, many species are already shifting their distributions towards higher altitudes and latitudes [10], but there have been few observations of this in tropical Asia [1,8].
concessions in southern Myanmar and rubber plantations in south-eastern Myanmar (Tanintharyi Division, Mon State, and Kayin State) have had significant impacts on biodiversity [23,25].
Existing threats to biodiversity in Myanmar can be made more serious by climate change, through both direct impacts, like the loss of suitable habitat for species and reduced resilience in ecosystems, and through indirect impacts on humans and their dependence on the products and services produced by natural ecosystems [26]. Climate change impacts on biodiversity in Myanmar are still understudied, however. Between 1981 and 2010, the mean annual temperature increased by 0.25 • C per decade, within inland regions warming faster than coastal regions, and the mean annual rainfall also increased in most regions [18]. According to the Global Climate Risk Index 2020, Myanmar was the country second most affected by extreme weather events between 1999 and 2018, although 95% of the damage and casualties were caused by a single event, Cyclone Nargis, in 2008 [27]. The number of drought and flood events has increased, as has the intensity and frequency of cyclones. Cyclones Mala (2006), Nargis (2008), and Giri (2010) were the most severe and damaging cyclones Myanmar has experienced. Moreover, during the summer of 2010, 1482 cases of heat-related disorders and 260 heat-related deaths due to extreme high temperature were reported, and in July and August 2015, 1.6 million people were displaced as the result of flooding and landslides [28].

Myanmar's Protected Area System
Myanmar's first protected areas were established as wildlife sanctuaries around Buddhist monasteries since the 11th Century and some of these were transformed into formal protected areas to protect endangered wildlife species during the colonial period [29]. However, the major expansion in the number of PAs and the total area protected has occurred since 1980. There are currently 45 protected areas, covering 5.85% of the country's area, under the management of the Nature and Wildlife Conservation Division (NWCD) within the Ministry of Natural Resources and Environmental Conservation [30]. Among the 45 current PAs, eight (Hkakaborazi National Park, Htamanthi Wildlife Sanctuary, Indawgyi Lake Wildlife Sanctuary, Alaungdaw Kathapa National Park, Inlay Lake Wildlife Sanctuary, Meinmahla Kyun Wildlife Sanctuary, Lampi Marine National Park, and Natmataung National Park) have been recognized as ASEAN Heritage Parks (AHPs) for their particular biodiversity value or uniqueness within ASEAN countries [30].
The PA system is intended to represent the full complement of the country's biogeographic regions, but in practice, their effectiveness in conserving biodiversity is reduced by additional factors related to size, geographic representation, inadequate funding and management capacity, weak policy, and the regulatory framework [26]. Some PAs are close to national borders, like Khakaborazi National Park, Lenya National Park (with Namtok Huay Yang in Thailand), Tanintharyi Nature Reserve (with Kaengkrachan Forest Complex in Thailand), but still lack transboundary protected area management, which can play a crucial role in preserving biodiversity. Three of the 45 protected areas (Htaung Wi Taung, Thamihla Kyun Wildlife Sanctuary, and Eaisarthaya Cave-Geographic-Features Significant Area Myanmar) were omitted from the analyses in this study because they are very small.
Optimization of the spatial configuration of protected area networks would ideally be based on species distribution data [17], but such data are often missing, very patchy, or inaccessible in Myanmar. In this paper, therefore, we used a spatial modeling approach based on a statistically derived bioclimatic stratification [31,32] to understand the spatial implications of projected climate change for Myanmar's protected area network by 2050 and 2070. Changes in the distribution of bioclimatic zones and strata are used as indicators of potential impacts on species and ecosystems. This approach avoids the "cascade of model uncertainties" associated with the modeling of biological impacts, even in areas where the species distribution data are available [32]. The main objectives were: 1.
To summarize the temporal and spatial patterns of climate change projected for Myanmar by 2050 and 2070.

2.
To classify and map the climate types currently present in Myanmar and project the changes in the areas and spatial distributions of these climate types by 2050 and 2070.

3.
To assess the representativeness of the protected area system in terms of the current climate types present in Myanmar.

4.
To assess the impact of projected climate change on the representation of climate types within individual protected areas and in the protected area system as a whole.

5.
To identify the likely impacts of climate change on the protection of biodiversity in Myanmar's protected areas.

Bioclimatic Stratification
A spatial modeling approach based on a statistically derived bioclimatic stratification was used to predict and understand the spatial implications of projected climate conditions for Myanmar by the year 2050 and 2070. The spatial shifts of bioclimatic zones and strata were evaluated to estimate the potential impacts of climate change across Myanmar. The construction of the Global Environmental Stratification (GEnS), used here, is described in detail by Metzger et al. [33] and its use in the interpretation, understanding, and communication of global climate change projections is discussed by [32]. The GEnS is intended to be a globally consistent classification of land into relatively homogenous units, based on a statistical clustering of climate variables [33]. The GEnS classifies the world's land surface into 125 relatively homogeneous bioclimatic strata, aggregated into 18 zones, based upon a statistical analysis of current climate data . The zones have descriptive names while the strata have unique codes. To produce the original GEnS, Metzger et al. identified a subset of 36 biophysically relevant bioclimatic variables based on a statistical screening of 42 variables available from various climate datasets [33]. Principal Components Analysis revealed that the first three principal components, explaining greater than 99% of the total variation, were determined by only four variables.
Degree days >0 • C [34] • Daily sum of annual degrees of temperature above 0 • C, reflecting latitudinal and altitudinal temperature gradients, and plant growth periods [34].
Monthly mean temperature seasonality [34] • Standard deviation of the monthly mean temperature distribution Potential evapotranspiration (PET) seasonality [35] • Standard deviation of the monthly mean PET distribution.

Modeling of Projected Future Climate Conditions
Projected future conditions for the years 2050 (average for 2041-2060) and 2070 (average for 2061-2080), also derived from WorldClim 1.4, were selected to represent a short to medium-term time frame relevant to the needs of ecosystem managers, planners, and other policy and decision makers. The environmental stratification of Myanmar, which is based on climate data from 1960 to 2000, was reconstructed for 2050 and 2070, based upon projections for future climate conditions from three CMIP5 Earth System Models (ESMs), HadGEM2-ES, CNRM-CM5 and GFDL-CM3, which have been previously recommended for the region [37,38]. We used two Representative Concentration Pathways (RCPs), RCP2.6 and RCP8.5, representing the low and high greenhouse gas concentration scenarios, respectively [39]. The socio-economic assumptions on which these scenarios were originally based are no longer realistic in 2020, but they still serve to bracket the potential range of radiative forcing by the end of the century. RCP2.6 is consistent with meeting the Paris Agreement's 2 • C global warming target. The ESM runs were downscaled using the Delta method to 1 km 2 resolution, as with the current climate conditions. Note that the stratification is based on the current climate and that future climates are assigned to the most similar current climate zone and stratum.

Projected Climate Change by 2050 and 2070
The mean annual temperature for Myanmar in 1960-2000 was 23.2 • C, mean maximum temperature of the warmest month was 32.  Table S1). The projected increase in mean annual temperature by 2050 is greatest in the eastern parts of Myanmar, approaching and exceeding 3.0 • C under the RCP8.5 scenario with GFDL-CM3 and HadGEM2-ES (Supplementary Material Figure S1). By 2070, mean annual temperature is predicted to increase by 1. In contrast to the generally consistent pattern of warming, there is a wider spread among models for projections of mean annual precipitation (Supplementary Material Table S1). By 2050, annual precipitation is predicted to increase by 1.  Figure S6). Note that these changes are small relative to current interannual variation, and to rainfall variability over the last 200 years reconstructed from a tree-ring chronology [40].

Bioclimatic Stratification of Myanmar under Current Conditions
Nine bioclimatic zones (Table 1) and 41 bioclimatic strata (Table 2) were identified in Myanmar (Figure 1), ranging from the Extremely Cold and Wet zone at the highest elevations in the north, with a single stratum, to the Extremely Hot and Xeric zone at low elevations in the central dry zone, with three strata, of which Q4 is the hottest and driest ( Table 1). The mean annual temperatures for these zones are correlated with their average elevation and range from −4.1 • C, for the coldest zone at an average elevation of 5409 m to 26.7 • C, for the hottest zone at an average elevation of 155 m. The maximum temperatures are also correlated with their average elevation and range from 7.0 • C for the coolest zone to 38.1 • C for the hottest zone. However, annual precipitation shows no consistent relationship with elevation. Of the nine zones, the three at the highest elevations (Extremely Cold and Wet, Extremely Cold and Mesic, and Cold and Mesic) each cover <1% of the total land area of Myanmar. The zones with the largest areas are Extremely Hot and Moist, covering 224,377 km 2 (34% of the total area), mostly in the south, and Hot and Mesic, covering 219,324 km 2 (33%), mostly further north.

Projected Changes in the Spatial Distribution of Bioclimatic Zones and Strata by 2050 and 2070
By the year 2050, substantial spatial displacement of the bioclimatic zones is seen under both RCPs (Supplementary Material Figure S7). Under RCP2.6, there is a large expansion in the extents of the Extremely Hot and Xeric zone (from 77,215 km 2 to 93,411-137,028 km 2 ) and Extremely Hot and Moist zone (from 224,377 km 2 to 245,484-267,000 km 2 ), and also in the Hot and Mesic zone (from 219,324 km 2 to 223,454-233,226 km 2 ), except with HadGEM2-ES, which predicts a small decrease (from 219,324 km 2 to 218,834 km 2 ), while all other zones decrease with all models (Supplementary Material Table S2).  Table S2). Changes in strata are larger than those for zones, reflecting their narrower definitions, but are model-dependent (Supplementary Material Table S4). Four new strata appear by 2050, but most only with GFDL-CM3, and several disappear. Overall, 9-13% of the total area of Myanmar will change zone under RCP2.6 and 11-23% under RCP8.5, while 21-33% under RCP2.6 and 27-39% under RCP8.5 will change stratum (Supplementary Material Table S3).
In general, projected changes by 2070 are in the same direction as 2050, but larger for both zones and strata (Supplementary Material Figure S8). Under RCP2.6, there is a large expansion in the extent of the Extremely Hot and Xeric zone (from 77,215 km 2 Table S4). Changes in strata are similar to but larger than those by 2050 (Supplementary Material Table S5). Overall, 8-14% of the total area of Myanmar will change zone under RCP2.6 and 17-32% under RCP8.5, while 20-35% under RCP2.6 and 34-52% under RCP8.5 will change stratum ( Table 3). Table 3. Land area and percentage of total land area in Myanmar that is projected to change its bioclimatic zone or stratum by 2050 and 2070, with three earth system models and two RCPs.

Climates and Climate Change in the Protected Areas
The bioclimatic zones and strata were overlain on a map of the protected areas. All nine zones and 37 out of 41 bioclimatic strata are found in PAs. However, representation of the zones (Table 4) and strata (Supplementary Material Table S6) varied greatly. Of the whole area protected, 47% (19,513 km 2 ) is in the Hot and Mesic zone, 30% (12,564 km 2 ) in the Warm Temperate and Moist zone, and 11% (4555 km 2 ) in the Extremely Hot and Moist zone ( Table 4). All other zones contributed <5% each. The five coolest zones had >30% of their total area protected, while the two driest zones had <1%. Projected changes in climate within the protected areas are similar to those for Myanmar as a whole, with most zones shifting upslope (Supplementary Material Table S7). By 2050, 64-70% of all protected areas are projected to shift at least partly to different bioclimatic zones under RCP2.6 and 67-75% under RCP8.5, while 81-82% are projected to shift at least partly to new strata under RCP2.6 and 80-82% under RCP8.5 (Supplementary Material Table S9). In most cases, the majority of each protected area remains in the same zone, but 0-5% of all protected areas are predicted to shift completely to different zones under RCP2.6 and 2-12% under RCP8.5, and 10-19% are predicted to shift completely to different strata under RCP2.6 and 17-26% under RCP8.5 (Supplementary Material  Table S10). The projected changes are similar but larger by 2070 (Supplementary Material Tables S8, S11 and S12). Several new strata are also projected to appear within the Cold and Mesic, Hot and Dry, Extremely Hot and Moist, and Extremely Hot and Xeric zones under both RCPs by 2050 and 2070 (Supplementary Material Tables S13 and S14). In the protected area system as a whole, under both RCPs, there are large declines in the areas of the four coolest zones and large increases in the areas of the two hottest zones by 2050 (Supplementary Material Table S7), but with the amount of change depending on the model. These changes are generally projected to increase further by 2070 (Supplementary  Material Table S8). The spread among individual reserves is very broad: 0-100% of each PA shifting to a new zone and/or stratum (Supplementary Material Tables S15-S18). Both the smallest and largest percentage changes are in PAs with a small elevational range and/or total area, where the entire reserve either switches or does not.

Discussion
The bioclimatic stratification of Myanmar highlights the climatic diversity in the country, from the Extremely Cold and Mesic zone, above the alpine treeline on the tallest mountains in the north of the country, to the Extremely Hot and Moist zone, until recently largely occupied by tropical lowland rainforest, in the south, and the Extremely Hot and Xeric zone, with remnant patches of savanna and deciduous forests, in the center of the country. Representation of the bioclimatic zones and strata in protected areas is currently very uneven, with >30% of the areas of four coolest zones being protected and <1% of the two driest. To a large extent this reflects human pressures within these zones, with the cool high mountains sparsely populated, while the drier areas have had dense populations for centuries. The hyperdiverse Extremely Hot and Moist zone in the south is also underrepresented in terms of percentage area (2%) and is now under threat from expanding plantations [25]. The percentage of the total land area protected in Myanmar (<6%) is also small, by both regional and global standards, and well below both the 17% area target (Aichi target 11) agreed by the Convention on Biological Diversity for 2020 and the proposed 30% area target for 2030 [41]. The three earth system models used in this study all project an acceleration of recent warming trends across the whole of Myanmar, but with a fairly large spread (<1.5 • C) among models in the amount of warming. This spread is even wider for rainfall, in terms of both the amount and the spatial pattern of increase and decreases. In most projections and over most of Myanmar, rainfall is projected to increase, but the projected increases are small relative to both current interannual variation and variability over the last 200 years [40]. Changes in the bioclimatic stratification are therefore dominated by the increases in temperature, resulting in an upwards shift in average elevation for all zones and strata. The hottest zones increase in area while the cooler zones decline, with the coldest disappearing with two models. Changes in the strata are greater, reflecting their narrower bounds, but mostly model-dependent. Up to a third of Myanmar's land area will change bioclimatic zone by 2070, depending on the model and RCP, while up to half will change stratum. Projected changes within the protected area system are similar to those in the country as a whole, but individual reserves are highly variable, with a complete switch of bioclimatic zone or stratum in some small reserves, as well as some larger ones with a low elevational range.
The consequences of these climatic changes for biodiversity depend on how effective the bioclimatic stratification is as a proxy for species and ecosystems, both now and in the future. Too little biodiversity data is available in Myanmar to validate this assumption for the present day, but there is support from studies in similar ecosystems in southwest China [31] and the transboundary Kailash Sacred Landscape of China, India, and Nepal [42], as well as studies in other parts of the world. Validation of future predictions is not possible, but theory, paleoecological evidence, and some observations of responses to recent climate change suggests that the populations of many well-dispersed species will track changes in climate across the landscape [14,15]. However, poorly dispersed species and those with long life-cycles will not be able to keep up. In particular, most of the individual trees that will dominate Myanmar's forests in 2050 and 2070 are already growing and cannot move, although a majority of 20 tree species studied in Natma Taung National Park had a higher proportion of juveniles at the upper end of their ranges, suggesting that their populations will eventually shift upslope [43]. Failure to track rapid climate change creates ecosystems that are not in equilibrium with the climate of the time, with consequences that are currently unclear, but are likely to be include slower growth and increased vulnerability to pests and diseases [14].
An additional complication comes from the increase in carbon dioxide concentrations, which is not only the largest single driver of climate change, but also has a direct impact on plant physiology and thus on plant growth, competition, and vegetation [44]. This means that bioclimate alone cannot predict future vegetation structure and species composition, which will also depend on the CO 2 concentration. In other words, future analogues of modern climates are not necessarily ecologically equivalent. Rising atmospheric CO 2 does not impact animals directly, but they will be impacted indirectly through changes in vegetation structure and composition. A recent modeling study which simulated the impacts of climate change on vegetation in South Asia (including Myanmar), with and without increasing CO 2 , found that simulations with increasing CO 2 resulted in transitions from savanna into forest and deciduous forest into evergreen forest which did not occur in the absence of elevated CO 2 [44]. The vegetation model used (aDGVM2) does not include nutrient limitation, so the impacts of elevated CO 2 may be overestimated, but woody invasion of savannas in other parts of the world has been attributed, in part, to this mechanism [45].
The disappearance from Myanmar of the coldest bioclimatic zone, Extremely Cold and Wet, will have little direct impact on biodiversity, since this represents the summit zone of Mt Khakaborazi, which is barren rock and ice. In contrast, the large declines in the areas of the next three coolest zones, in both the country as a whole and the protected area system, will substantially reduce the area available for species adapted to high-mountain forest and alpine habitats in Myanmar. Upward shifts of several hundreds meters in steep topography, where they represent horizontal movements of a kilometer or two, may be within the dispersal capacities of most plant and animal species, but the area available declines with altitude on most mountains, and reaches zero at the summit. On isolated high mountains, such as Mt Victoria (Natma Taung) (3074 m) in southwest Myanmar, endemic species found only near the summit face potential mountain-top extinction. At the other extreme, species occurring in protected areas with little or no elevational range, because of flat topography (such as Chatthin and Shwesettaw Wildlife Sanctuaries) or small size (such as Chungponkan Wildlife Sanctuary, Lawkanada Sanctuary, and Wetthikan Bird Sanctuary), are threatened by the total loss of the bioclimatic zones or strata to which they are adapted, as the entire protected area undergoes a shift. Species adapted to open forests and savanna, such as the endangered Eld's deer (Rucervus eldi thamin) may be particularly vulnerable to woody encroachment and canopy closure, as a result of climate change and/or rising atmospheric CO 2 (see above).
The use of climatic data as a surrogate for biodiversity is not ideal, since bioclimatic zones and strata are not, in themselves, targets for conservation. This approach was necessitated by the patchy availability of biodiversity data in Myanmar. As more such data becomes available, it should be possible to calibrate the bioclimatic stratification in a way that makes it more useful for conservation planning [46]. Where biodiversity data is lacking, the addition of geological information would be an improvement on using just climate as a surrogate. Myanmar's extensive karsts, for example, support numerous narrow-range endemic species whose presence could not be predicted from climate alone.
It would be possible to make recommendations for additional protected areas and the expansion of existing ones based on this study, although recommendations based on climate variables alone should only be a first step. The vulnerability of the existing protected areas depends not only on their exposure to climate change, as assessed here, but also on their resilience (indicated by size, isolation, topographic variability, etc.) and capacity for adaptation [47]. Clearly, both the total area and the representativeness of the protected area system need to be increased, and connectivity across climatic gradients should be enhanced to permit species movements [11]. Extensive restoration of degraded vegetation, both passively (by removing the causes of degradation) and actively (by planting), may be needed, particularly in some lowland and drier areas [48]. Reintroduction of locally extirpated animal species may be practical where hunting can be controlled.
However, the protected area system in Myanmar is not currently limited by technical knowledge, but rather reflects, to a large extent, the legacies of decades of armed internal conflicts, some of which continue at a lower level, despite cease fires and peace agreements. These conflicts have limited the collection of biodiversity data and continue to make it very difficult to create new protected areas agreed by both the central and regional governments. As in many other countries, biodiversity protection in Myanmar is intimately linked with a variety of political, social, and economic issues, and progress in conservation will depend on progress in solving all these. Experience in some of these countries suggests that the best way forward is to take the technical knowledge-in this case, from climate change science-as a starting point and then to focus on policy, planning, and management issues in dialogue with major stakeholders [49].
Finally, we focus in this paper on protected areas, but the same climatic changes will also impact agricultural and urban areas, both directly and through their impacts on the supply of water and other services from natural ecosystems. Biodiversity conservation is easier in remote, unpopulated, areas, but arguably most important near to where most people live. Natural and restored ecosystems can not only provide a reliable source of water, but also reduce the risk from floods, cyclones, and other extreme weather events, regulate local climates, and provide accessible recreational and tourism opportunities, and associated economic benefits for local people [50].

Conclusions
Myanmar's current protected area system is small for such a biodiversity-rich country and not representative of the country's great climatic diversity. Projected climate change will result in a general upward shift in climate zones to higher altitudes, as well as model-dependent changes in rainfall, which are mostly small compared with current variability. Threats to biodiversity are expected for mountain-top endemics on isolated peaks, for the biotas of small protected areas and those with a low elevational range, and, in combination with rising CO 2 levels, to species dependent on open habitats vulnerable to woody encroachment. Biodiversity data are needed to refine predictions based on climate alone, but the major factors currently limiting the needed extensions to the protected area system are social and political, not scientific.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2225-1154/8/9/99/s1, Table S1: Projected changes in mean annual temperature and precipitation for Myanmar between 1960-2000 and 2050 and 2070, Table S2: Projected change in areal extent and mean elevation of bioclimatic zones and their upward shift by 2050, Table S3: Projected change in areal extent and mean elevation of bioclimatic zones and their upward shift by 2070, Table S4: Projected change in areal extent and mean elevation of bioclimatic strata and their upward shift by 2050, Table S5: Projected change in areal extent and mean elevation of bioclimatic strata and their upward shift by 2070, Table S6: Representation of bioclimatic strata within protected areas, Table S7: Projected change in areal extent and mean elevation of bioclimatic zones in protected areas and their upward shifts by 2050, Table S8: Projected change in areal extent and mean elevation of bioclimatic zones in protected areas and their upward shifts by 2070, Table S9: Percentage of all protected areas shifting to a different zone by 2050, Table S10: Percentage of all protected areas shifting to different strata by 2050, Table S11: Percentage of all protected areas shifting to a different zone by 2070, Table S12: Percentage of all protected areas shifting to different strata by 2070, Table S13; Projected change in areal extent and mean elevation of bioclimatic strata in protected areas and their upward shifts by 2050, Table S14. Projected change in areal extent and mean elevation of bioclimatic strata in protected areas and their upward shifts by 2070, Table S15: Percentage shift of each protected area to different bioclimatic zone for 2050 and 2070, Table S16: Percentage shift of each protected area to different bioclimatic strata  for 2050 and 2070, Table S17: Projected change in areal extent of bioclimatic zones in each protected area by 2050, Table S18: Projected change in areal extent of bioclimatic zones in each protected area by 2070, Figure S1: Change in mean annual temperature as projected for the year 2050, Figure S2: Change in mean annual temperature as projected for the year 2070, Figure S3: Change in maximum temperature of the warmest month as projected for the year 2050, Figure S4: Change in maximum temperature of the warmest month as projected for the year 2070, Figure S5: Change in mean annual precipitation as projected for the year 2050, Figure S6: Change in mean annual precipitation as projected for the year 2070, Figure S7: Bioclimatic stratification of Myanmar based on spatially interpolated weather station data and as projected for 2050, Figure S8: Bioclimatic stratification of Myanmar based on spatially interpolated weather station data and as projected for 2070.
Author Contributions: Conceptualization, T.N. and R.T.C.; methodology, T.N., R.J.Z., and R.T.C.; formal analysis, T.N. R.J.Z., and R.T.C.; writing-original draft preparation, T.N. and R.T.C.; writing-review, editing, and revision, all authors; supervision, R.T.C. All authors have read and agreed to the published version of the manuscript.