Next Article in Journal
The Relationship between the Frontal QRS-T Angle on ECG and Physical Activity Level in Young Adults
Next Article in Special Issue
Developing an Evaluation System Suitable for Coastal Rural Houses’ Characteristic Style and Its Inspiration for Rural Revitalization: Case Study of Rongcheng in Shandong Province
Previous Article in Journal
Level of Physical Activity, Sleep Quality and Use of Personal Protective Equipment of Students at Wroclaw Medical University during the COVID-19 Pandemic
Previous Article in Special Issue
Spatiotemporal Changes of Cultivated Land System Health Based on PSR-VOR Model—A Case Study of the Two Lake Plains, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Myanmar’s Land Cover Change and Its Driving Factors during 2000–2020

1
State Key Laboratory of Resources and Environmental Information System, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
School of Geosciences, Yangtze University, Wuhan 430100, China
3
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(3), 2409; https://doi.org/10.3390/ijerph20032409
Submission received: 5 December 2022 / Revised: 23 January 2023 / Accepted: 25 January 2023 / Published: 29 January 2023
(This article belongs to the Special Issue Land Use Change and Its Effects on Regional Sustainable Development)

Abstract

:
Land use/cover change (LUCC) research occupies an important place in the study of global change. It is important for the ecological protection and long-term development of a place. Current research is lacking in the study of dynamic changes at the national level in Myanmar over long time periods and sequences. Quantitative research on the driving factors of LUCC is also lacking. This paper uses the GLC_FCS30 (Global Land-Cover product with Fine Classification System) dataset and socio-economic statistical data in Myanmar to conduct the study. The dynamic change process of LUC (land use/cover) was investigated using the land use dynamic degree, land use transfer matrix, and Sankey diagram. Principal component analysis was used to derive the main drivers of LUCC. The drivers were quantified using multiple linear stepwise regression analysis and specific factors were analyzed. The spatial scope of the study is Myanmar, and the temporal scope is 2000–2020. Results: (1) In 2020, the spatial distribution of LUC in Myanmar shows predominantly forests and croplands. Forests account for 56.64% of the country’s total area. Agricultural land accounts for 25.59% of the country’s total area. (2) Over the time scale of the study, the trend of LUCC in Myanmar showed significant shrinkage of evergreen broad-leaved forest and deciduous broad-leaved forest (a total shrinkage of −3.34 × 104 km2) and expansion of the other land types. (3) Over the time scale of the study, the dynamic changes in LUCC in Myanmar most occurred as an interconversion between two land types, such as between cropland and deciduous broad-leaved forest, evergreen broad-leaved forest and shrubland, deciduous broad-leaved forest and shrubland, evergreen broad-leaved forest and evergreen needle-leaved forest, and evergreen broad-leaved forest and deciduous broad-leaved forest. (4) The dynamics of LUC in Myanmar is mainly influenced by the socio-economic level of the country. Among them, the impact of agricultural level is the most obvious. Specifically, Myanmar’s LUCC is mainly driven by urban population, urbanization rate, industrial value added, food production, and total population. Our research will enable the Myanmar government to make more scientific and rational land management and planning and to make more informed decisions. After understanding the basic situation of LUCC in Myanmar, the hydrological effects, biodiversity changes, and ecological service function changes due to land change in the region can be explored. This is the direction of future research.

1. Introduction

LUCC studies are indispensable in global change research. LUCC is capable of exhibiting climate change characteristics. The natural ecosystems on the Earth’s land surface can be affected by human activities. This impact can also be manifested through LUCC. Therefore, LUCC plays the role of a bridge and a link in the process of interaction between humans and nature [1]. As early as 1995, in order to increase understanding of the dynamic processes of LUCC and study their relationship with total environmental change of the Earth, IGBP and IHDP established the LUCC program. In 2006, IGBP and IHDP initiated the Global Land Project, which took over and replaced the LUCC program. Its main purpose is to synthesize and organize the ideas, knowledge, and methods from the research conducted by scientists in the land systems science community [2]. In 2014, ICSU and ISSC launched the 10-year Future Earth Initiative in collaboration with several international organizations such as UNESCO and UNEP. The program can provide society with the critical knowledge it needs to respond to the problems posed by environmental change on all of Earth’s surfaces. At the same time, it can provide the world with the opportunity to make the transition to sustainable development. It was the result of the gradual succession of global environmental change research strategies such as LUCC and GLP and met the core needs of the international community to emphasize linking land science research with addressing social challenges [3].
The basis of LUCC research is the LULC (Land Use & Land Cover) dataset. In the 1980s, Wilson and Henderson-Sellers of the University of Liverpool [4] produced a global land use/cover dataset for use in GCM climate models. After entering the 21st century, remote sensing observation technology and land mapping technology have developed rapidly. As a result, more and more LULC dataset products with global scale, higher spatial resolution, and shorter update cycles have appeared in the world, such as the IGBP DISCover product of the US Geological Survey [5] (https://daac.ornl.gov/, accessed on 30 September 2022), the GLC product of the EU Joint Research Center [6] (https://forobs.jrc.ec.europa.eu/, accessed on 29 September 2022), and the European Space Agency’s GLOBCOVER product [7,8] (https://www.esa.int/, accessed on 30 September 2022). These LULC products are mostly land mapping products with a spatial resolution of 300 m, 500 m, or even 1 km, and an update cycle of 5 years or even 10 years. In the past ten years, with the emergence of higher temporal and spatial resolution earth observation platforms and the accumulation of historical observation data, especially the further development of remote sensing cloud technology, higher spatial resolution (30 m, even 10 m), and higher time, the resolution (annual scale) LULC dataset was produced. Other datasets include the Globeland30 dataset produced by the Ministry of Natural Resources of China [9] (http://www.globallandcover.com/, accessed on 30 September 2022), the GLC_FCS30 dataset produced by the Academy of Aerospace Technology, Chinese Academy of Sciences [10,11,12] (https://data.casearth.cn/, accessed on 30 September 2022), the LSV10 data product produced by the European Space Agency [13] (https://esa-worldcover.org/, accessed on 30 September 2022), and ESRI 2020 Land Cover 10 m (ESRI10) data product produced by ESRI Corporation of the United States [14] (https://www.arcgis.com/, accessed on 30 September 2022).
Currently, the global LULC is changing rapidly under the influence of both the natural environment and human activities. This may break the balance of ecosystems and reduce biodiversity [15]. Therefore, when conducting research on LUCC, the analysis of how its drivers affect the sustainable development of the region is the focal point. Many studies have shown that the key element driving change in LULC is the economic activity of humans in society. Liu et al. [16] analyzed LUCC in China during the 1990s. Cropland area increased the most (44.88 million mu). The grassland area decreased the most (51.55 million mu). Policy regulation and economic drivers are the main causes of LUCC and its spatial and temporal variation. Moisa et al. [17] studied the LUCC in Addis Ababa municipality in Ethiopia. With the increase in population, vegetation cover and other land cover types have been transformed into building land. From 1991 to 2021, the agricultural land cover decreased by 66.8 km2 and the vegetated land cover decreased by 25.7 km2. The built-up area of the city was 96.6 km2 (18.3%) in 1991 and increased to 277.2 km2 (52.6%) in 2021. The development of urban society was the main reason for LUCC. Hu et al. [18] analyzed the LUCC in Ha Long City, Vietnam, since the founding of the country. The share of forest land shrank by 26.3%. The share of both urban land and industrial and mining storage has expanded by 4.3 times. The main reason for this change was the improvement in socioeconomic level. These studies provide examples that can be used for the extraction, research framework, and indicator methodology of LULC change impact factors for a particular country or region.
In studying land use change, we often use land use dynamic degree to analyze the extent of land type change and the land use transfer matrix to analyze the land type transformation process. These methods were used by Talebi et al. [19], Zhao et al. [20], and Xing et al. [21] in their study of land use dynamic degree in the east Meskine region of Ardabil province, Manwan Dam of the Lancang River, and low elevation areas of Ethiopia. When performing the driving mechanism of LUCC, there are generally multiple driving factors. The advantage of performing analysis with multiple drivers is that it makes the analysis results more comprehensive. The disadvantages are that it will generate data redundancy and a computationally intensive and complex analysis process. Moreover, there may be a correlation between multiple drivers. This may lead to less reliable final analysis results. Jiang et al. [22] and Liu et al. [23] used both principal component analysis and linear regression models in their study of the driving factors of LUCC in the Cuno River Basin and Chinese Jiangxi province.
Myanmar has long suffered from instability, sectarian violence, corruption, weak infrastructure, and a long history of colonial exploitation. The country’s economic and cultural level is generally relatively backward. According to a study by the United Nations Development Programme, Myanmar is extremely backward in human development in 2020, ranking 147th out of 189 countries in the world [24]. Previously, LULC studies related to Myanmar were usually conducted at smaller regional scales (such as watersheds, protected areas, counties, cities, and other administrative regions). For example, the study by McGinn et al. [25] studied LUCC in the Chindwin River Basin, Myanmar. During 1999–2019, the forest area in this area decreased by about 2%. The mining area increased by 0.38%. The area of agricultural land decreased by 2.1% overall. The increase in deforestation and mining activities has led to major problems such as water pollution, sedimentation, and river channel changes. Zaehringer et al. [26] studied LUCC in southeastern Myanmar from 1990 to 2017 and showed that the area shifted from smallholder rotation to natural rubber, betel nut, cashew, and oil palm trees. From this, it could be concluded that commercial interests were the core force for the huge loss of forest resources. The building land area increased from 1988 (730.64 ha) to 2018 (3607.20 ha). The urban forest area decreased from 2008 (4901.63 ha) to 2018 (3558.36 ha).
Although there are many studies on LUCC in Myanmar, there are still many problems. Among the many LULC data products that exist, the classification of various land types has mixed results. For example, ESRI has released the 2020 Land Cover product. It has good classification results for agricultural land, forests, water bodies, and building sites, but poor classification results for shrubs, wetlands, or fallow lands [27]. Additionally, Myanmar is relatively backward. Its socio-economic data and climate data are inadequate. This has prevented scholars from conducting fine-grained LUCC process portrayal and pattern description studies on it. For the analysis of the driving factors, most of the studies in Myanmar also stay at the stage of qualitative analysis. There are fewer quantitative analyses. Helen et al. [28] from the University of Tokyo, Japan, analyzed the LUCC process in the Myanmar city of Pyin Oo Lwin between 1988 and 2018. It was largely driven by demographic and socio-economic factors. This study provides evidence for urban green space protection, green infrastructure development planning, etc., but the study does not quantify the influencing factors. Moreover, the spatial scope of LUCC research is not large, and the time scale is also short. In addition, there is also a lack of in-depth discussions on the analysis of driving mechanisms based on mathematical statistics for national sustainable development and ecological protection.
In this paper, we use the long time series, high-resolution GLC_FCS30 dataset as the LUC study data. We use climate and national socio-economic statistics as the driving factor research data. The dynamic processes of spatial and land type change of LUC in Myanmar were analyzed using land use dynamic degree, Sankey diagram, and land use transfer matrix. Principal component analysis was applied to derive the main drivers of LUCC in Myanmar. Multiple linear stepwise regression in statistics was used to quantify the drivers of LUCC in Myanmar and to analyze the specific drivers. This study can enable the Myanmar government to conduct land planning and designing in a more scientific and practical way. Additionally, it can enable Myanmar to achieve sustainable development faster and protect the ecological environment better.

2. Data and Methods

2.1. Study Area Overview

Myanmar is in the northwestern part of the Indochina Peninsula (9°58′–28°31′ N, 92°20′–101°11′ E) (Figure 1). Myanmar is bordered by China in the north. It is connected to Laos and Thailand in the east. It is bordered to the southwest by the Andaman Sea and the Bay of Bengal in the Indian Ocean. It is bordered by Bangladesh and India to the northwest. The total land area is 67.66 × 104 km2, and the coastline is 0.28 × 104 km. Myanmar has high terrain in the north and low terrain in the south. It is surrounded by mountain ranges in the north, west, and east. It has high mountainous areas in the north (the Savage Mountains, also known as the Kachin Mountains, the Kumaon Ridge, and the Hukawng Valley Mountains). The Rakhine Mountains and the Naga Hills are in its western part. The Shan State plateau is to the east of it. Khyakabo Peak, near the Chinese border, is 5881 meters above sea level and is the highest peak in Myanmar. The Irrawaddy alluvial plain lies between its western mountains and eastern plateau. The terrain is low and flat, with an average elevation of 200 meters to 600 meters.
Myanmar is a tropical monsoon climate region. Its temperature does not vary much throughout the year. The whole year of Myanmar can be divided into three seasons: hot season, rainy season, and cool season. The hot season lasts from March to mid-May. The rainy season is from mid-May to October. The cool season is from November to February. The coldest month is January. It has an average temperature of 20–25 °C. The hottest months are April and May. It has an average temperature of 25–30 °C. Myanmar receives abundant rainfall. The most rainfall occurs in June, July, and August. The southwest monsoon is prevalent during this period. Precipitation in Myanmar is followed by May, September, and October. The rainfall in various places from May to October accounts for about 90–95% of the annual rainfall. Most of Myanmar has an annual rainfall of 4000 mm or more. The central part of Myanmar (Magway) is a rain shadow area with relatively low annual rainfall, but it also reaches 1000 mm.
Myanmar is divided into 7 provinces, 7 states, and 2 centrally administered municipalities. The seven provinces are Yangon, Mandalay, Bago, Magway, Sagaing, Ayeyarwady, and Tanintharyi. The seven states are Chin State, Shan State, Kayah State, Mon State, Kachin State, Karen State, and Rakhine; the two municipalities are Naypyidaw and Yangon. The province is the main inhabited area of the Burmese, and the state is the inhabited area of various ethnic minorities. Before 2005, the capital of Myanmar was Yangon, which is located in the southwest and borders the Andaman Sea. In 2005, the Myanmar government moved the capital to Naypyidaw in the central part of the Sittang Valley between the Bago Mountains and the Bonnong Mountains.
The social and economic level of Myanmar is relatively backward. The industry and service industries are underdeveloped. It is an agricultural country. The rural population accounts for about 70% of the country’s total population. Rural people make their living mainly from agricultural farming. The main crops grown include conventional rice, wheat, corn, and peas. It also grows a large number of cash crops, such as rubber, sugarcane, cotton, and palm. In recent years, peas have surpassed rice as the most important agricultural product exported by Myanmar. Myanmar is rich in forestry resources. Forests cover more than 56% of the country’s total land area. Myanmar produces and exports a large amount of teak, and the output of teak ranks first in the world. Myanmar produces and exports large quantities of teak. The production of teak ranks first in the world. In Myanmar’s industry, the labor-intensive processing and manufacturing industry represented by the textile and garment industry and some light industry sectors are developing fast. The development of all kinds of heavy industries is seriously lagging behind.

2.2. Data Sources

We downloaded the GLC_FCS30 dataset for the study from the China Earth Big Data Science Project Data Sharing Service (https://data.casearth.cn/, accessed on 30 October 2022). The spatial resolution of this dataset is 30 meters. It covers the period 1985–2020. Among them, every 5 years is 1 time period. It has a total of 8 time periods [12]. There are two levels of land cover types. Among them, there are 9 first-level types and 30 second-level types (Table 1), with classification accuracies of 82.5% and 68.7%, respectively [10]. There are 22 second-class land types in Myanmar (Table 2). In Myanmar, forests cover more than half of the total area of the country. We would like to study it in more detail and conclude the distribution and changes of each specific type of forest in Myanmar. So, only forests used Tier 2 land types in the study. All other land types used Tier 1 land types. The areas of bare areas and permanent snow and ice were extremely small. We merged bare areas into grasslands. We merged permanent snow and ice into wetlands & water bodies.
We calculated meteorological data (X1, X2) for Myanmar for 2000–2020 from day-by-day data from GLDAS-2.1 [29] and the PERSIANN-CDR dataset [30]. Their spatial resolution is 0.25 radians. We downloaded data on socio-economic development (X3–X9) from the World Bank (https://www.shihang.org/zh/home, accessed on 1 November 2022). We downloaded the rest of the data (X10–X16) from the World Food and Agriculture Organization (http://www.fao.org/, accessed on 3 November 2022).
Wood production includes wood fuel, sawn wood, veneer logs, veneer sheets, plywood, and other products. Vegetable production includes the production of green garlic, peppers, onions and shallots, other vegetables, etc. Oil crop production includes the production of sesame seeds, groundnuts, sunflower seeds, castor oil seeds, mustard seeds, cotton seeds, cashew nuts, etc., of which groundnuts production accounts for about 50%. Food production includes the production of wheat, soya beans, sorghum, rice, potatoes, millet, maize, lentils, peas, beans, cassava, sweet potatoes, cereals, etc., among which rice accounts for about 80%. Fruit production includes the production of areca nuts, coconuts, sugar cane, bananas, mangoes, guavas, mangosteens, other fruits, etc., of which sugar cane accounts for about 75%.
By analyzing various statistics of Myanmar, we established a system of driving factors of LUCC in Myanmar (Table 3).

2.3. LUCC Analysis Methods

Generally, there are two types of land use dynamic degree. One is the integrated land use dynamic degree. The other is the single land use dynamic degree.
As the name implies, the integrated land use dynamic degree (D) provides an integrated representation of the overall degree of change in all land use types within a study area [31]. It is able to express the quantitative characteristics of the variation. The mathematical expression is:
D = { i = 1 n ( Δ U i U i ) } × 1 T × 100 %
In the expression, D is the integrated land use dynamic degree in the research period; Δ U i is the area of change of class i LUC throughout the study; U i is the area of class i LUC at the beginning of the study; T is the research period; and n is the number of LUC types.
The single land use dynamic degree (K) can be used to quantify the changes in the number of a particular land type within a certain time frame in different subdivisions of the study area. Therefore, it can compare the changes in the number of land types in different subdivisions, and then predict its development trend [32]. The formula is:
K = U b U a U a × 1 T × 100 %
In the formula, K is the dynamic degree of a particular LUC type during the time period of this study; Ua denotes the area of a particular LUC at the start of the study; Ub denotes the area of the same LUC at the end of the study; and T denotes the time period studied.
The Markov model is a statistical model. It has been widely applied to various scientific fields. The land use transfer matrix is built on its basis. The land use transfer matrix allows the transfer between different land types to be quantified [33]. Moreover, the land use transfer matrix has to be used with a specific time interval when it is carried out. The land use transfer matrix can reveal the process by which various land types transform into each other during this time interval [34]. The land use transfer matrix is expressed mathematically as follows:
S i j = [ S 11   S 12 S 1 n S 21 S 22 S 2 n S n 1 S n 2 S n n ]
where S denotes the area of various types of land; i denotes the LUC type at the beginning of the study; j denotes the LUC type at the end of the study; and n denotes the number of LUC types.
We obtained areas of all land types within Myanmar by using the GLC_FCS30 dataset from 2000 to 2020. We calculated two land use dynamic degrees in Excel by Equations (1) and (2). We imported the land use data of Myanmar from 2000 to 2020 in the GLC_FCS30 dataset into ArcGIS 10.2. The land use transfer matrix was obtained by raster calculation.

2.4. Driving Factors Analysis Methods

Principal component analysis (PCA) uses an orthogonal transformation. It transforms the original random variable with correlated components into a new variable with uncorrelated components. Algebraically speaking, it converts the covariance array of the original variables into a diagonal array. Geometrically speaking, it transforms the original system of variables. After this operation, it makes it a new orthogonal system. In the case of a scattered distribution of sample points, it refers to the most open orthogonal direction of the scattering. In turn, the dimensionality reduction of the multidimensional variable system is achieved [35]. It converts multiple indicators into several comprehensive indicators. At the same time, it retains as much data and information in it as possible. Additionally, it is also able to eliminate redundancy arising from multiple data [36]. Through the principal component analysis method, the extracted principal component factors can be given new names according to the unique meanings reflected by professional knowledge and indicators. Thus, reasonable explanatory variables are obtained [37]. The selection method for each main factor is as follows:
{ F 1 = c 11 Z 1 + c 12 Z 2 + + c 1 n Z n F 2 = c 21 Z 1 + c 22 Z 2 + + c 2 n Z n F n = c n 1 Z 1 + c n 2 Z 2 + + c n n Z n
where Fi denotes principal component i, i = 1, 2, …, n; c denotes the eigenvector corresponding to the eigenvalues of the covariance matrix; and Z denotes the value obtained after normalization of the original variables. In this formula, for each i, there is c i 1 2 + c i 2 2 + + c i n 2 = 1 .
There are multiple drivers that affect LUCC. We use a multiple linear stepwise regression approach to conduct the analysis. In this method, variables are introduced in and out. It introduces variables individually. When each of the variables is entered, the variables that have been selected for inclusion are checked individually. When the previously entered variables become no longer significant due to the later introduced variables, they are removed. This ensures that only variables with significant effects are included in the regression model before each new variable is added. The introduction and elimination of variables was repeated in this manner. This operation is performed until neither influential significant independent variables are entered into the regression model nor influential insignificant independent variables are removed from the regression model [38]. Multiple linear stepwise regression analysis methods can build the most optimum or appropriate regression model. This allows a more in-depth study of the dependencies between variables [39]. This method can be used to make predictions for the dependent variable. It can improve the accuracy of the results. We also used this method to conduct the analysis of each driving factor presented in the paper. The expressions are as follows:
Y = β + α 1   X 1 + α 2   X 2 + + α n X n
where α1, α2, …, αn denote relevant coefficients; β denotes the constancy term; Y denotes the dependent variable; and X 1 X n denote the independent variables.
We imported Myanmar land use data and driving indicator data from 2000 to 2020 into SPSS software. After setting the relevant parameters, we let SPSS perform the corresponding operations. The final results of principal component analysis and multiple linear stepwise regression analysis were obtained. From the results, we can see the main driving forces and specific drivers of LUCC in Myanmar from 2000 to 2020.

3. Results Analysis

3.1. Spatial Distribution

In 2020, it can be seen from the statistical map of various types of LUC areas in Myanmar (Figure 2) that Myanmar is dominated by forests (accounting for 56.64% of the area) and cropland (accounting for 25.59% of the area). The combined area of the two land types accounts for more than 82% of the total national land area.
Specifically, forests have the widest distribution of area (43.18 × 104 km2, 56.64%). Among them, most are evergreen broad-leaved forests (19.06 × 104 km2, 25.01%) and deciduous broad-leaved forests (18.00 × 104 km2, 23.62%), and a few are evergreen needle-leaved forests (6.11 × 104 km2, 8.02%). The second largest LUC type is cropland (19.51 × 104 km2, 25.59%). The third largest LUC type is shrubland (10.96 × 104 km2, 14.37%). The area of grassland, impervious surfaces, wetlands and water bodies, etc., is small (total 2.59 × 104 km2, 3.39%). They contain grassland (0.48 × 104 km2, 0.63%), impervious surfaces (0.48 × 104 km2, 0.64%), and wetlands and water body (1.62 × 104 km2, 2.13%).
Cropland is mostly located in the central Irrawaddy alluvial plain, such as Sagaing, Magway, Mandalay, etc. There is also a small amount of distribution in the southwestern coastal areas, such as Ayeyarwady, Yangon, Mon State, etc.
Forests are mostly located in the Shan State plateau in the east, the northern region, and the mountainous region in the west. Among them, the evergreen broad-leaved forest is mostly located in Kachin State in the northern mountain area and Tanintharyi, Bago, and other places in the south. There is also a small amount of distribution in Chin State, Magway, Rakhine, and other places in the mountainous and hilly area in the west and Shan State area in the eastern Shan State plateau. The deciduous broad-leaved forest is mostly located in the alluvial plains of the central Ayeyarwady River, such as Bago, Sagaing, Magway, etc., and there is also a small amount of distribution in the Shan State area of the eastern Shan State plateau. The evergreen needle-leaved forest is mostly located in Kachin State in the northern mountain region, and the rest are scattered in the western mountain region and the eastern Shan State plateau.
In addition, the shrubland is mostly located in the Rakhine and Chin State areas in the western mountains, and the rest are scattered in the eastern Shan State plateau, northern mountainous areas, southern Mon State, Karen State, Tanintharyi, and other places. The grasslands are scattered in Kachin State, Sagaing, and other places in the northern mountains, and in the Shan State area of the eastern Shan State plateau. Impervious surfaces are concentrated in large cities, such as Sagaing and Mandalay in the central part, Yangon and Moulmein in the southern part, and others are distributed in the capitals or larger cities of various provinces or states. Wetlands and water bodies are generally distributed in points in Kachin State and Sagaing in the northern mountains, Shan State in the eastern Shan State plateau, Rakhine in the western mountains, and Bago, Yangon, Ayeyarwady, Tanintharyi, and other places in the south.

3.2. Analysis of Land Dynamic Change

During 2000–2020, the area of forest in Myanmar has shrunk (−3.34 × 104 km2, −0.36%), and the area of other land types has expanded. Among them, deciduous broad-leaved forest (−2.43 × 104 km2, −0.59%) and evergreen broad-leaved forest (−1.86 × 104 km2, −0.44%) showed a decreasing trend, and evergreen needle-leaved forest (+0.94 × 104 km2, +0.91%) showed an increasing trend. Cropland (+0.98 × 104 km2, +0.26%), wetlands and water bodies (+0.24 × 104 km2, +0.85%), and impervious surfaces (+0.94 × 104 km2, +0.91%) all showed a trend of area expansion. The most expanded area is shrubland (+1.82 × 104 km2, +1.00%), and the least expanded area is grassland (+0.09 × 104 km2, +1.19%).
In terms of the spatial distribution of land use dynamic degree, the main LUC changes in the first ten years (2000–2010) occurred in the central region, and the LUC changes in the eastern and western regions were weaker. In the next ten years, the LUC changes were most concentrated in the western and southern regions, and then in the eastern region, the intensity of LUC changes also climbed rapidly (Figure 3).
During 2000–2005, the spatial distribution of integrated land use dynamic degree in Myanmar was high in the central plains, low in the eastern plateaus and western mountains, and basically followed the law of decreasing with the elevation of the terrain. The integrated land use dynamic degree in most provinces is below 33%. Among them, the lowest is Rakhine (10.87%), and the highest is Yangon (56.06%).
During 2005–2010, the integrated land use dynamic degree increased in the mountains of western Myanmar and the plains of central Myanmar and decreased in southern Tanintharyi. The integrated land use dynamic degree in half of the provinces is above 33%. Among them, the lowest is Shan State (8.79%). The highest is Mon State (2156.12%), followed by Karen State (755.98%).
During 2010–2015, the integrated land use dynamic degree in Myanmar declined as a whole, and most provinces were below 33%. Among them, the lowest is Sagaing (9.98%). The highest is Yangon (1141.52%), followed by Kayah State (260.41%).
During 2015–2020, the integrated land use dynamic degree in Myanmar increased as a whole, showing a trend of increasing with the elevation of the terrain. Only Mandalay (19.09%) is below 20%. There are four provinces with a dynamic degree exceeding 100%, the highest being Chin State (1269.01%), followed by Kayah State (730.66%), Shan State (208.89%), and Rakhine (105.38%).

3.3. Source and Destination

During 2000–2020, the land area where transfer between different LUC types occurred in Myanmar was 18.51 × 104 km2. From 2000 to 2010, the land area where LUC types were transferred in Myanmar was 14.39 × 104 km2, and from 2010 to 2020, the area of transferred land was 16.84 × 104 km2. The conversion process mainly occurs between two land use/cover types (Figure 4), including cropland and deciduous broad-leaved forest, evergreen broad-leaved forest and shrubland, deciduous broad-leaved forest and shrubland, evergreen broad-leaved forest and evergreen needle-leaved forest, and evergreen broad-leaved forest and deciduous broad-leaved forest.
The cropland continues to show a trend of expansion. During 2000–2020, converted cropland to other land was a total of 2.60 × 104 km2, mainly converted into the deciduous broad-leaved forest (−55.82%, −1.45 × 104 km2), wetlands and water bodies (−15.29%, −0.40 × 104 km2), and shrubland (−12.95%, −0.34 × 104 km2). The total area of other land converted to cropland is +3.58 × 104 km2, mainly from the deciduous broad-leaved forest (+67.13%, +2.40 × 104 km2) and shrubland (+13.99%, +0.50 × 104 km2). Cropland shrinkage mainly occurs in the central plains such as Mandalay, Magway, Ayeyarwady, and Yangon; cropland expansion mainly occurs in the northern mountainous Sagaing, Kachin State, etc., the southern coastal Karen State, Mon State, Tanintharyi, etc., and the eastern plateau Shan State (Figure 5a).
The overall forest showed a shrinking trend. During 2000–2020, the total area of forest converted to other land was 7.16 × 104 km2. The total area of other land converted to forest was 3.40 × 104 km2.
Among them, the total area of evergreen broad-leaved forest transformed into other land types was 4.72 × 104 km2, mainly transformed into shrubland (−48.10%, −2.27 × 104 km2), evergreen needle-leaved forest (−26.36%, −1.24 × 104 km2), and deciduous broad-leaved forest (−19.36%, −0.91 × 104 km2). The total area of other land converted to the evergreen broad-leaved forest was 2.84 × 104 km2, mainly from shrubland (+37.51%, +1.06 × 104 km2), deciduous broad-leaved forest (+32.16%, +0.91 × 104 km2), and evergreen needle-leaved forest (+27.25%, +0.77 × 104 km2). The shrinking evergreen broad-leaved forests are mainly distributed in Kachin State, Sagaing, Rakhine, etc., in the northern and western mountains and Tanintharyi, Karen State, etc., in the south. The expansion is mainly distributed in Chin State, etc., in the western mountainous areas, and Shan State, etc., in the eastern plateau.
The total area of deciduous broad-leaved forest transformed into other land types was 6.35 × 104 km2, mainly transformed into cropland (−37.83%, −2.40 × 104 km2), shrubland (−28.37%, −1.80 × 104 km2), and evergreen needle-leaved forest (−16.93%, −1.08 × 104 km2). The total area of other land converted to the deciduous broad-leaved forest was 3.91×104 km2, mainly from cropland (+37.02%, +1.45 × 104 km2), shrubland (+24.58%, +0.96 × 104 km2), evergreen broad-leaved forest (+23.34%, +0.91 × 104 km2), and evergreen needle-leaved forest (+13.90%, +0.54 × 104 km2). The shrinkage of deciduous broad-leaved forests mainly occurred in Shan State and Kayah State in the eastern plateau, Sagaing and Chin State in the northern and western mountains, and Bago, Mon State, and Karen State in the south-central region. The expansion mainly occurred in Magway in the central plain, in the northern Sagaing and Kachin State in the mountains, Shan State in the eastern plateau, etc.
The total area of evergreen needle-leaved forest transformed into other land types was 1.55 × 104 km2, mainly transformed into the evergreen broad-leaved forest (−49.93%, −0.77 × 104 km2) and deciduous broad-leaved forest (−35.14%, −0.54 × 104 km2). The total area of other land converted to the evergreen needle-leaved forest was 2.51 × 104 km2, mainly from the evergreen broad-leaved forest (+49.62%, +1.24 × 104 km2) and deciduous broad-leaved forest (+42.91%, +1.08 × 104 km2). The shrinkage of evergreen coniferous forests mainly occurred in the Kachin State in the northern mountainous area and the Shan State in the eastern plateau. The expansion mainly occurred in the Chin State in the western mountainous area and the Shan State in the eastern plateau (Figure 5b).
The shrubland showed a continuous expansion trend. During 2000–2020, the total area of shrubland converted to other land was 2.68 × 104 km2, mainly transformed into the evergreen broad-leaved forest (−39.74%, −1.06 × 104 km2), deciduous broad-leaved forest (−35.93%, −0.96 × 104 km2), and cropland (−18.72%, −0.50 × 104 km2). The total area of other land transformed into shrubland was 4.52 × 104 km2, mainly from the evergreen broad-leaved forest (+50.25%, +2.27 × 104 km2) and deciduous broad-leaved forest (+39.91%, +1.80 × 104 km2). The shrinkage of the shrubland has a small distribution in the central plains of Mandalay, Mandalay, and other areas. Other areas are scattered. Shrubland expansion mainly occurred in the eastern plateau, western and northern mountains, and southern coastal areas (Figure 5c).

3.4. Climate Change and Socio-Economic Development

During 2000–2020, the overall climate in Myanmar showed a slightly warm and dry trend (Figure 6a). Over the past 20 years, the average annual temperature for many years has been 27.04 °C, showing an overall upward trend (+0.09 °C/10a), reaching a peak in 2019 (28.26 °C). The average annual total precipitation for many years is 1799.93 mm, showing a fluctuating decline overall trend (−118.8 mm/10a). The total annual precipitation in 2017 was the least (1333.98 mm).
From 2000 to 2020, the overall population of Myanmar showed a continuous increase trend (Figure 6b). The total population of Myanmar increased from 2000 (4671.97 × 104) to 2020 (5440.98 × 104). These years have increased by 769.01 × 104. It grew at an annual rate of 7.6‰. The rural population increased from 2000 (3409.37 × 104) to 2020 (3746.60 × 104). These years have increased by 337.23 × 104. It grew at an annual rate of 4.7‰. The urban population increased from 2000 (1262.6 × 104) to 2020 (1694.38 × 104). These years have increased by 431.78 × 104. It grew at an annual rate of 1.48%. The rate of urbanization continues to accelerate from 2000 (27.03%) to 2020 (31.14%). These years have increased by 4.12 percentage points.
From 2000 to 2020, the overall economic development of Myanmar showed an upward trend (Figure 6c). GDP increased from 2000 (USD 68.49 billion) to 2020 (USD 789.30 billion). These years have increased by USD 720.81 billion. It grew at an annual rate of 13.00%. The agricultural value added increased from 2000 (USD 35.54 billion) to 2020 (USD 152.65 billion). These years have increased by USD 117.10 billion. It grew at an annual rate of 7.56%. The industrial value added increased from 2000 (USD 6.35 billion) to 2020 (USD 229.14 billion). These years have increased by USD 222.79 billion. It grew at an annual rate of 19.64%.
From 2000 to 2020, the overall forestry development in Myanmar showed an upward trend (Figure 6d). Except for a slight decline in natural rubber production after 2018, it rose steadily in other years and reached its peak in 2018 (27.55 × 104 t). Wood charcoal production fluctuated slightly during the rising process, reaching peaks (43.82 × 104 t) in 2019 and 2020. Wood production increased rapidly before 2005 and from 2010 to 2013, showed a downward trend after 2013, and reached a peak in 2013 (4690.84 × 104 m3).
From 2000 to 2020, the overall development of agriculture in Myanmar showed an increasing trend (Figure 6e). Vegetable production declined slightly in 2012 and then increased year by year. In 2016, it began to decline slowly and fluctuated, maintaining at about 490 × 104 t, and reached a peak in 2015 (522.25 × 104 t). Oil crop production continued to rise before 2010, then fluctuated slowly, and reached a peak in 2010 (330.57 × 104 t). Food production continued to rise before 2010, then declined, and reached a peak in 2010 (4061.62 × 104 t). The overall fruit production is on the rise, from 757.57 × 104 t in 2000 to 1536.36 × 104 t in 2020, an increase of 778.79 × 104 t. It grew at an annual rate of 3.60%.

3.5. Analysis of Driving Factors

In Table 3, the magnitudes of the various driving factors we have selected are different. To address this issue, we adopted a standard deviation normalization method for the socio-economic and climate data of Myanmar.
Through principal component analysis, we obtained its component matrix (Table 4). From this we know that from 2000 to 2020, the driving factors of Myanmar’s economy, society, and climate can be divided into three dimensions. F1 can be classified as the socio-economic dimension. It shows an extremely strong positive correlation with eight driving factors. They are total population (X3), rural population (X4), urban population (X5), urbanization rate (X6), industrial value added (X9), wood charcoal production (X11), vegetable production (X13), and fruit production (X16). F2 can be classified as a climate dimension. It shows an extremely strong positive correlation with the total annual precipitation (X2). F3 is classified as the dimension of agriculture (grain planting), which shows a strong negative correlation with food production (X15). In terms of contribution rate, F1 (56.29%) is far greater than F2 (27.14%) and F3 (13.46%). National development (F1) has more impact on LUCC than climate (F2) and food planting (F3).
In the multivariate linear stepwise regression analysis, we exclude the wetland water type. The independent variables are three principal components (F1, F2, and F3). The dependent variable is the area of the remaining seven LUC types. The independent variable and dependent variable are analyzed, respectively, and the regression equation between them is established (Table 5).
From the regression equation, it can be seen that the area of cropland has an obvious positive relationship with socio-economic (F1) and agriculture (staple food cultivation) (F3). It shows that the cropland area has increased while Myanmar’s socio-economic and agricultural development has taken place. The area of evergreen broad-leaved forests and deciduous broad-leaved forests have an obvious negative correlation with socio-economic (F1). It shows that with the development of the social economy in Myanmar, the area of evergreen broad-leaved forests and deciduous broad-leaved forests decreased. The area of shrubland, grassland, and impervious surfaces have a significant positive correlation with socio-economic (F1), respectively. It shows that with the development of society and economy, the area of shrubland, grassland, and impervious surfaces are all increasing. In summary, the main influence on the LUCC situation in Myanmar is socio-economic. The small amount of impact is on agriculture (staple food cultivation). There is no credible statistical relationship between the evergreen needle-leaved forest area and the main components F1, F2, and F3.
Furthermore, based on principal component analysis, we obtained the 10 key driving factors, which were total annual precipitation (X2), total population (X3), rural population (X4), urban population (X5), urbanization rate (X6), industrial value added (X9), wood charcoal production (X11), vegetable production (X13), food production (X15), and fruit production (X16). We used them as independent variables. We used the area of the seven aforementioned LUC types as dependent variables. Then, we established a regression equation using multiple linear stepwise regression methods (Table 6).
From Table 6, we can know that there is a significant positive correlation between the cropland area and the urbanization rate (X6). It has a significant negative correlation with food production (X15). The accelerated rate of urbanization has promoted population growth, and the demand for food, vegetables, and fruits will increase. Due to these markets’ needs, the cropland area will increase. After 2010, while the area of cropland continued to expand, the output of staple food and vegetables continued to decline, but the output of oil crops and fruits continued to increase (Figure 6). This shows that a considerable part of the arable land has been transferred from traditional staple food and vegetable cultivation to high-value-added oil crops and fruit cultivation. Both the area of evergreen broad-leaved forests and deciduous broad-leaved forests showed a significant negative correlation with the urbanization rate (X6). This shows that with the increase in urbanization rate and the construction and expansion of urban land, the evergreen broad-leaved forest and deciduous broad-leaved forest located in the low-lying, suitable for living and development areas are cut down first, reducing their area. The area of shrubland has a significant positive correlation with the urban population (X5). Urbanization increases the size of shrubland by destroying forests and converting them into shrubland. The grassland area has a significant positive correlation with industrial value added (X9). It also has a significant negative correlation with food production (X15). Industrial development will increase the change from forest to shrubland. This will also allow the change of shrubland to grassland. To meet the growing demand for food, it depends on the expansion of the cropland area. The expansion of cropland area is first to occupy the grassland that is easy to cultivate, and then to occupy the forest. There is an obvious positive correlation between the area of impervious surfaces and the total population (X3). The overall population increases, resulting in the need for more places to live. This requires building more houses, building roads, etc. This will increase the area of impervious surfaces. After the analysis, we did not obtain regression equations for the area of evergreen needle-leaved forest and 10 key drivers.

4. Discussion

4.1. LUCC and Its Influencing Factors

Our research shows that from 2000 to 2020, the area of evergreen broad-leaved forests and deciduous broad-leaved forests in Myanmar decreased. The area of the rest of the land types increased. Among them, the absolute area of shrubland increased the most. This result is the same as the regional-scale research results of Wang et al. [40] in Yangon, Karki et al. [41] in Inle Lake, and McGinn et al. [25] in the Chindwin River Basin. All of their studies also indicate that the forest area decreased, and the areas of shrubland, grassland, impervious surfaces, and wetlands all increased within the study area. The changing trend of LUCC in Myanmar is similar to that of other Indochina countries. Research by Hu et al. [18], Zhang et al. [42], and Niu et al. [43] in Vietnam, Laos, Cambodia, and other countries also found the phenomenon of shrinking forests, shrubland, and the expansion of impervious surfaces.
Studies have shown that 60% of LUCC is related to direct human activities, and the other 40% is indirectly related to climate change [44]. Studies on LULC changes in Vietnam, Cambodia, and Laos also found that national economic development and human activities have a greater impact on LUCC in these countries. Among the factors that have a more pronounced impact are GDP, population, and urbanization rates [42,43,45]. Our research also shows that at the national scale, socio-economic and agricultural development have a greater impact on LUCC in Myanmar. The impact of climate change is smaller.
In terms of forest resources in Myanmar, infrastructure development (such as the construction of China–Myanmar oil and gas pipelines) and logging for commercial purposes have had a significant impact on the country’s forest resources [46,47,48]. Since the 1990s, the Myanmar government has successively promulgated laws and regulations such as “Standards and Indicators for Sustainable Forestry Management”, “Forest Law of the Union of Myanmar”, “Forestry Regulations of the Union of Myanmar”, “Regulations for the Implementation of National Forest Logging”, and other laws and regulations, and from April 1 in 2014, a total ban on the export of logs [49,50]. The protection of forest resources not only depends on the formulation and implementation of laws and regulations but also depends on improving the livelihood of farmers and establishing a sustainable forest industry chain. Therefore, there are two main important tasks for Myanmar in terms of land management and achieving sustainable development. First, Myanmar should make use of the good water and heat conditions provided by the tropical and subtropical climate in which it is located. Secondly, Myanmar should establish a modern forest industry with scientific afforestation, nurturing, logging and high value-added forest products processing industry as its core.
At the same time, the support of the international community for the sustainable development of Myanmar’s national land is also crucial. The IFAD, FAO, and the Myanmar government have reached relevant agreements on these aspects, that is, to strengthen the sustainable management and governance of land, forests, water sources, and ecosystems, and to improve local communities and the ability of farmers to cope with natural and man-made disasters, climate change, and transboundary diseases [51,52]. With the gradual construction and development of the China–Myanmar Economic Corridor planned in the “Belt and Road” initiative, the LUC pattern and evolution process in various parts of Myanmar will be increasingly affected by socio-economic development factors and specific construction projects. How will the land use/cover change in Myanmar and other Indochina Peninsula countries? What impact will the LUCC change have on regional ecological services and farmers’ livelihoods? Can regional sustainability be effectively improved? This will be important content of LUCC research at the national and regional scales in Myanmar in the future.

4.2. Uncertainty in Data and Results

A study has been conducted for the GLC_FCS30 data accuracy. After the study, it achieves an overall accuracy of 82.5% [10]. In this dataset, the forest accuracy is 94%, the cropland accuracy is 88%, the shrubland accuracy is 56.8%, the grassland accuracy is 67.3%, the impervious surfaces accuracy is 79.3%, and the water body accuracy is 83.8%. Actually, in most land cover datasets based on remote sensing interpretation, the accuracy of cropland and shrubland is not too high. Cropland and shrubland are easily confused with other types, such as cropland with grassland, shrubland with forest, and shrubland with grassland [53]. The accuracy analysis carried out by Wang et al. [54] based on four sets of LUC datasets based on different satellite remote sensing information sources and different resolutions in the Indochina Peninsula also shows that the degree of confusion of grassland, shrubs, bare land, and construction land is high, and the overall accuracy is low. The impact of the accuracy of the GLC_FCS30 dataset on the results of this paper is inevitable.
The rainfall in various parts of Myanmar from May to October accounts for about 90–95% of the annual precipitation. In different seasons of the year, the range and boundaries of wetlands and water bodies are affected by annual precipitation, and the changes are very large [55]. Therefore, there are great uncertainties in the area statistics and change analysis of wetlands and water bodies in Myanmar in this study. Therefore, we excluded wetlands and water bodies from the LUCC driving factors analysis.
The two driving factor analysis methods used in this study also have limitations. We simply established a linear regression relationship between the values of the driving factors and the values of the LUC areas. However, it does not correspond to it in space. We also did not take into account human factors such as culture, policy, and people’s ideological changes. Therefore, the results of our analysis of driving factors presented so far are coarse resolution in spatial terms. It also does not present spatial differences in the driving factors of LUCC within different provinces, watersheds, and agricultural zones in Myanmar. Moreover, the driving factors presented in this paper are macro-oriented and fail to reflect the influence of Myanmar’s unique Buddhist culture, military government system, national and provincial development plans, and specific construction projects. These shortcomings will depend on the development of a more detailed spatial analysis model of the driving mechanism in the future and the development of multi-scale and cross-scale LUCC driving mechanism analysis cases.

5. Conclusions

The data for our study is the GLC_FCS30 dataset and official Myanmar socio-economic statistics. We analyzed the dynamic change characteristics of LUCC in Myanmar from 2000 to 2020. The driving mechanisms of LUCC in Myanmar were summarized. It also suggested possible directions for Myanmar’s forest resource conservation and LUCC research. The present study is a regional reference that can be used for global LUCC studies. It is also the first study of long-time series LUCC processes and driving factors at the national scale in Myanmar.
As a result of our study, it is known that the main LUC types in Myanmar are forests and cropland. During 2000–2020 in Myanmar, forests showed a significant shrinkage trend, and the rest of the land types showed a significant expansion trend. Socio-economic development had a greater impact on LUCC change in Myanmar, and climate change had a smaller impact on LUCC change in Myanmar. We emphasize that in order to protect Myanmar’s disappearing forest resources, the Myanmar government and the international community should work together and take effective measures such as returning farmland to forests, improving national literacy, and developing a green economy.
In the study, we established a method to analyze spatial pattern changes in LUCC. The research basis of this method is the GLC-FCS30 dataset, land use dynamic degree, and Sankey map. We also established an analysis method for the driving factors. The research basis of this method is the principal component analysis and multiple linear stepwise regression analysis. In the future, if scholars want to conduct research on similar LULC changes in other countries or regions, they can learn from and apply these methods. However, the current study is still inadequate in the refinement of the driving mechanism and spatial modeling of the driving mechanism analysis, and more in-depth research is needed in the future. In addition, after understanding the basic situation of land change in Myanmar, the hydrological effects, biodiversity changes, and ecological service function changes due to land change in the region are also directions for future research.

Author Contributions

Writing—original draft, Y.W.; software, Y.W. and X.N.; data curation, Y.W. and X.N.; conceptualization, Y.H.; methodology, Y.H. and Y.W.; writing—review and editing, Y.H. and Y.W.; supervision, Y.H. and L.Z.; project administration, Y.H.; funding acquisition, Y.H.; formal analysis, H.Y. and L.Z.; investigation, L.Z. and H.Y.; visualization, X.N. and H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (42130508) and the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA20010202).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their sincere thanks to the anonymous reviewers because the comments and suggestions were of great help to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, J.; Kuang, W.; Zhang, Z.; Xu, X.; Qin, Y.; Ning, J.; Zhou, W.; Zhang, S.; Li, R.; Yan, C.; et al. Spatiotemporal characteristics, patterns and causes of land use changes in China since the late 1980s. Acta Ecol. Sinica 2014, 69, 3–14. [Google Scholar] [CrossRef]
  2. Verburg, P.H.; Crossman, N.; Ellis, E.C.; Heinimann, A.; Hostert, P.; Mertz, O.; Nagendra, H.; Sikor, T.; Erb, K.-H.; Golubiewski, N.; et al. Land system science and sustainable development of the earth system: A global land project perspective. Anthropocene 2015, 12, 29–41. [Google Scholar] [CrossRef] [Green Version]
  3. Zhou, T.-J. Implementation of the Future Earth programme in China: Chinese scientists at the forefront. Adv. Clim. Chang. Res. 2015, 6, 81–83. [Google Scholar] [CrossRef]
  4. Wilson, M.F.; Henderson-Sellers, A. A global archive of land cover and soils data for use in general circulation climate models. J. Climatol. 1985, 5, 119–143. [Google Scholar] [CrossRef]
  5. Loveland, T.; Brown, J.; Ohlen, D.; Reed, B.; Zhu, Z.; Yang, L.; Howard, S.; Hall, F.; Collatz, G.; Meeson, B. ISLSCP II IGBP DISCover and SiB Land Cover, 1992–1993; ORNL DAAC: Oak Ridge, TN, USA, 2009. [Google Scholar]
  6. Bartholomé, E.; Belward, A.S. GLC2000: A new approach to global land cover mapping from Earth observation data. Int. J. Remote. Sens. 2005, 26, 1959–1977. [Google Scholar] [CrossRef]
  7. Defourny, P.; Schouten, L.; Bartalev, S.; Bontemps, S.; Cacetta, P.; De Wit, A.; Di Bella, C.; Gérard, B.; Giri, C.; Gond, V. Accuracy Assessment of a 300 m Global Land Cover Map: The GlobCover Experience; International Center for Remote Sensing of Environment (ICRSE): Tucson, AZ, USA, 2009. [Google Scholar]
  8. Defourny, P.; Vancutsem, C.; Bicheron, P.; Brockmann, C.; Nino, F.; Schouten, L.; Leroy, M. Globcover: A 300 m global land cover product for 2005 using Envisat MERIS time series. In Proceedings of the ISPRS Commission VII Mid-Term Symposium: Remote Sensing: From Pixels to Processes, Enschede, The Netherlands, 8–11 May 2006. [Google Scholar]
  9. Jun, C.; Ban, Y.; Li, S. Open access to Earth land-cover map. Nature 2014, 514, 434. [Google Scholar] [CrossRef] [Green Version]
  10. Zhang, X.; Liu, L.; Chen, X.; Gao, Y.; Xie, S.; Mi, J. GLC_FCS30: Global land-cover product with fine classification system at 30 m using time-series Landsat imagery. Earth Syst. Sci. Data 2021, 13, 2753–2776. [Google Scholar] [CrossRef]
  11. Zhang, X.; Liu, L.Y.; Wu, C.S.; Chen, X.D.; Gao, Y.; Xie, S.; Zhang, B. Development of a global 30 m impervious surface map using multisource and multitemporal remote sensing datasets with the Google Earth Engine platform. Earth Syst. Sci. Data 2020, 12, 1625–1648. [Google Scholar] [CrossRef]
  12. Liu, L.; Zhang, X.; Gao, Y.; Chen, X.; Shuai, X.; Mi, J. Finer-Resolution Mapping of Global Land Cover: Recent Developments, Consistency Analysis, and Prospects. J. Remote. Sens. 2021, 2021, 5289697. [Google Scholar] [CrossRef]
  13. Zanaga, D.; Van De Kerchove, R.; De Keersmaecker, W.; Souverijns, N.; Brockmann, C.; Quast, R.; Wevers, J.; Grosu, A.; Paccini, A.; Vergnaud, S. ESA WorldCover 10 m 2020 v100; The European Space Agency (ESA): Paris, France, 2021. [Google Scholar]
  14. Karra, K.; Kontgis, C.; Statman-Weil, Z.; Mazzariello, J.C.; Mathis, M.; Brumby, S.P. Global land use/land cover with Sentinel 2 and deep learning. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS (2021), Brussels, Belgium, 12–16 July 2021; pp. 4704–4707. [Google Scholar]
  15. Millennium Ecosystem Assessment. Ecosystems and Human Well-Being: Biodiversity Synthesis; World Resources Institute: Washington, DC, USA, 2017. [Google Scholar]
  16. Liu, J.; Zhang, Z.; Zhuang, D.; Wang, Y.; Zhou, W.; Zhang, S.; Li, R.; Jiang, N.; Wu, S. A study on spatial-temporal dynamic changes of land-use and driving forces analyses of China in the 1990s. Geogr. Res. 2003, 22, 1–12. [Google Scholar]
  17. Moisa, M.B.; Dejene, I.N.; Roba, Z.R.; Gemeda, D.O. Impact of urban land use and land cover change on urban heat island and urban thermal comfort level: A case study of Addis Ababa City, Ethiopia. Environ. Monit. Assess. 2022, 194, 123456789. [Google Scholar] [CrossRef] [PubMed]
  18. Hu, W.; Su, F.; Wang, W.; Feng, X. The Characteristics of Land Use Change during the Different Periods in Halong City. J. Geo-Inf. Sci. 2017, 19, 570–579. [Google Scholar]
  19. Hossein, T.K.; Raoof, M.; Ali, A.M.; Khamsa, A.N.S. Temporal land use change and its economic values under competing driving forces in a diverse land use configuration. Arab. J. Geosci. 2022, 15, 20. [Google Scholar]
  20. Qinghe, Z.; Shiliang, L.; Li, D.; Shikui, D.; Zhifeng, Y.; Qi, L. Determining the influencing distance of dam construction and reservoir impoundment on land use: A case study of Manwan Dam, Lancang River. Ecol. Eng. 2013, 53, 235–242. [Google Scholar]
  21. Yuanyuan, X.; Yongdong, W.; Yuan, Y.; Qin, S.; Loje, H.M.; Bora, J.Z.; Huka, J.G. Dynamic Changes of the Bush Encroachment in Low Altitude Area of Ethiopia. J. Resour. Ecol. 2018, 9, 281–289. [Google Scholar]
  22. Xiaohui, J.; Haoqi, G.; Yuxin, L. Changes and driving factors of land use types in Kuye Basin in recent 40 years. South North Water Transf. Water Sci. Technol. 2020, 18, 17–38. [Google Scholar]
  23. Weilin, L.; Hao, H.; Shengnan, Z.; Xiang, L.; Bin, W.; Jingrong, Z.; Qing, Y.; Lina, L. Spatial and Temporal Variability Characteristics and Driving Mechanisms of Land Use in Jiangxi Province from 1980 to 2015. Res. Soil Water Conserv. 2023, 30, 1–8. [Google Scholar]
  24. UNDP. Human Development Report 2021–22; UNDP (United Nations Development Programme): New York, NY, USA, 2022. [Google Scholar]
  25. McGinn, A.J.; Wagner, P.D.; Htike, H.; Kyu, K.K.; Fohrer, N. Twenty years of change: Land and water resources in the Chindwin catchment, Myanmar between 1999 and 2019. Sci. Total Environ. 2021, 798, 148766. [Google Scholar] [CrossRef]
  26. Zaehringer, J.G.; Lundsgaard, H.L.; Thein, T.T.; Llopis, J.C.; Tun, N.N.; Myint, W.; Schneider, F. The cash crop boom in southern Myanmar: Tracing land use regime shifts through participatory mapping. Ecosyst. People 2020, 16, 36–49. [Google Scholar] [CrossRef]
  27. Ding, Y.; Yang, X.; Wang, Z.; Fu, D.; Li, H.; Meng, D.; Zeng, X.; Zhang, J. A Field-Data-Aided Comparison of Three 10 m Land Cover Products in Southeast Asia. Remote. Sens. 2022, 14, 5053. [Google Scholar] [CrossRef]
  28. Helen Pawel, J.M.; Alexandros, G. Land use change, carbon stocks and tree species diversity in green spaces of a secondary city in Myanmar, Pyin Oo Lwin. PLoS ONE 2019, 14, e0225331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Rodell, M.; Kato, H.; Zaitchik, B. Ongoing Development of NASA’s Global Land Data Assimilation System; American Geophysical Union: Washington, DC, USA, 2008. [Google Scholar]
  30. Ashouri, H.; Hsu, K.-L.; Sorooshian, S.; Braithwaite, D.K.; Knapp, K.R.; Cecil, D.L.; Nelson, B.R.; Prat, O.P. PERSIANN-CDR: Daily Precipitation Climate Data Record from Multisatellite Observations for Hydrological and Climate Studies. Bull. Am. Meteorol. Soc. 2015, 96, 69–83. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, G.; Wang, X.; Xiang, A.; Wang, X.; Wang, B.; Xiao, S. Spatial heterogeneity and driving factors of land use change in the middle and upper reaches of Ganjiang River, southern China. Chin. J. Appl. Ecol. 2021, 32, 2545–2554. [Google Scholar]
  32. Li, L.; Yang, J.; Jiang, D.; Li, J.; Hou, X. GIS based study on spatial temporal changes of land use in Wuding River Basin in the 1990s. Geogr. Res. 2005, 24, 527–534+656. [Google Scholar]
  33. Cui, X.; Guo, Y. Analysis on the spatio-temporal dynamic evolution of land use structure of western urban agglomerations in the past 40 years. J. Arid. Land Resour. Environ. 2022, 36, 16–24. [Google Scholar]
  34. Wu, L.; Yang, S.; Liu, X.; Luo, Y.; Zhou, X.; Zhao, H. Response analysis of land use change to the degree of human activities in Beiluo River basin since 1976. Acta Geogr. Sin. 2014, 69, 54–63. [Google Scholar]
  35. Zhang, P. A Comprehensive Evaluation Study Based on Principal Component Analysis; Nanjing University of Science & Technology: Nanjing, China, 2004. [Google Scholar]
  36. Li, Q.; Zhang, H.; Guo, S.; Fu, K.; Liao, L.; Xu, Y.; Cheng, S. Groundwater pollution source apportionment using principal component analysis in a multiple land-use area in southwestern China. Environ. Sci. Pollut. Res. 2020, 27, 9000–9011. [Google Scholar] [CrossRef]
  37. Li, X.; Chen, Z. Correctly Using SPSS Software for Principal Components Analysis. Stat. Res. 2010, 27, 105–108. [Google Scholar]
  38. He, X. Research on Multivariate Linear Model and Ridge Regression; Huazhong University of Science & Technology: Wuhan, China, 2005. [Google Scholar]
  39. You, S.; Yan, Y. Stepwise Regression Analysis and Its Application. Stat. Decis. 2017, 14, 31–35. [Google Scholar]
  40. Wang, Y.-C.; Hu, B.K.H.; Myint, S.W.; Feng, C.-C.; Chow, W.T.L.; Passy, P.F. Patterns of land change and their potential impacts on land surface temperature change in Yangon, Myanmar. Sci. Total Environ. 2018, 643, 738–750. [Google Scholar] [CrossRef]
  41. Karki, S.; Thandar, A.M.; Uddin, K.; Tun, S.; Aye, W.M.; Aryal, K.; Kandel, P.; Chettri, N. Impact of land use land cover change on ecosystem services: A comparative analysis on observed data and people’s perception in Inle Lake, Myanmar. Environ. Syst. Res. 2018, 7, 25. [Google Scholar] [CrossRef] [Green Version]
  42. Zhang, Y.; Niu, X.; Hu, Y.; Yan, H.; Zhen, L. Temporal and Spatial Evolution Characteristics and Its Driving Mechanism of Land Use/Land Cover Change in Laos from 2000 to 2020. Land 2022, 11, 1188. [Google Scholar] [CrossRef]
  43. Niu, X.; Hu, Y.; Lei, Z.; Wang, H.; Zhang, Y.; Yan, H. Spatial and Temporal Evolution Characteristics of Land Use/Cover and Its Driving Factor in Cambodia during 2000–2020. Land 2022, 11, 1556. [Google Scholar] [CrossRef]
  44. Song, X.-P.; Hansen, M.C.; Stehman, S.V.; Potapov, P.V.; Tyukavina, A.; Vermote, E.F.; Townshend, J.R. Global land change from 1982 to 2016. Nature 2018, 560, 639–643. [Google Scholar] [CrossRef] [PubMed]
  45. Niu, X.; Hu, Y.; Lei, Z.; Yan, H.; Ye, J.; Wang, H. Temporal and Spatial Evolution Characteristics and Its Driving Mechanism of Land Use/Cover in Vietnam from 2000 to 2020. Land 2022, 11, 920. [Google Scholar] [CrossRef]
  46. He, Z. The Change of Landscape Pattern and Its Driving Mechanism in the China-Myanmar Economic Corridor from 1995 to 2015; Yunnan University: Kunming, China, 2021. [Google Scholar]
  47. Zhou, J.; Li, S. Spatiotemporal Variations of Land Use/Cover in Important Node Area along the China-Myanmar Oil and Gas Pipeline during 2012–2015. Bull. Surv. Mapp. 2018, 3, 113–116. [Google Scholar]
  48. Zaw, N.T.; Paul, D.; McMoran, D.J.; Clive, M.; Genia, H. Patterns and Drivers of Deforestation and Forest Degradation in Myanmar. Sustainability 2021, 13, 7539. [Google Scholar]
  49. Ouyang, A. Legal System of Forest Protection in Myanmar. World Environ. 2010, 5, 58–59. [Google Scholar]
  50. Wen, L.; Wang, N. Exploration of the Forestry and the Protection Measures in Myanmar. J. Guangxi Univ. Natl. 2014, 36, 99–101. [Google Scholar]
  51. Huang, Q. Myanmar Government signed a framework agreement with the Food and Agriculture Organization of the United Nations. World Trop. Agric. Inf. 2019, 3, 26. [Google Scholar]
  52. Zeng, Y. Development Strategies for Rice Industry in Myanmar under the Perspective of Sino-Myanmar Cooperation. Chin. J. Trop. Agric. 2022, 42, 103–108. [Google Scholar]
  53. Pendrill, F.; Persson, U.M. Combining global land cover datasets to quantify agricultural expansion into forests in Latin America: Limitations and challenges. PLoS ONE 2017, 12, e0181202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wang, H.; Yan, H.; Hu, Y.; Xi, Y.; Yang, Y. Consistency and Accuracy of Four High-Resolution LULC Datasets—Indochina Peninsula Case Study. Land 2022, 11, 758. [Google Scholar] [CrossRef]
  55. Li, C.; Chen, J.; Ge, Y.; Chen, L.; Zhang, W.; Zhang, Y.; Zhao, Y. Land cover changes and their driving forces in Myanmar based on Globe-Land30. World Reg. Stud. 2021, 30, 478–489. [Google Scholar]
Figure 1. Location and topography map of Myanmar. AY: Ayeyarwady; BG: Bago; CS: Chin State; KC: Kachin State; KR: Karen State: KY: Kayah State; MD: Mandalay; MG: Magway; MS: Mon State; RK: Rakhine; SG: Sagaing; SS: Shan State; TT: Tanintharyi; and YG: Yangon.
Figure 1. Location and topography map of Myanmar. AY: Ayeyarwady; BG: Bago; CS: Chin State; KC: Kachin State; KR: Karen State: KY: Kayah State; MD: Mandalay; MG: Magway; MS: Mon State; RK: Rakhine; SG: Sagaing; SS: Shan State; TT: Tanintharyi; and YG: Yangon.
Ijerph 20 02409 g001
Figure 2. LUC map of Myanmar in 2020.
Figure 2. LUC map of Myanmar in 2020.
Ijerph 20 02409 g002
Figure 3. Change tendencies of LUC types of Myanmar during 2000–2020.
Figure 3. Change tendencies of LUC types of Myanmar during 2000–2020.
Ijerph 20 02409 g003
Figure 4. Conversion between different LUC types in Myanmar from 2000 to 2020.
Figure 4. Conversion between different LUC types in Myanmar from 2000 to 2020.
Ijerph 20 02409 g004
Figure 5. Spatial distributions of different LUC gains/losses from 2000 to 2020.
Figure 5. Spatial distributions of different LUC gains/losses from 2000 to 2020.
Ijerph 20 02409 g005
Figure 6. Development status in Myanmar from 2000 to 2020. (a) Climate factors; (b) demographic factors; (c) economic factors; (d) forestry factors; and (e) agricultural factors.
Figure 6. Development status in Myanmar from 2000 to 2020. (a) Climate factors; (b) demographic factors; (c) economic factors; (d) forestry factors; and (e) agricultural factors.
Ijerph 20 02409 g006aIjerph 20 02409 g006b
Table 1. LULC classification system of GLC_FCS30.
Table 1. LULC classification system of GLC_FCS30.
CodeLevel 1 ClassesLC IDLevel 2 Classes
1Cropland10Rainfed cropland
11Herbaceous cover
12Tree or shrub cover (orchard)
20Irrigated cropland
2Forest51Open evergreen broad-leaved forest
52Closed evergreen broad-leaved forest
61Open deciduous broad-leaved forest (0.15 < fc < 0.4)
62Closed deciduous broad-leaved forest (fc > 0.4)
71Open evergreen needle-leaved forest (0.15 < fc < 0.4)
72Closed evergreen needle-leaved forest (fc > 0.4)
81Open deciduous needle-leaved forest (0.15 < fc < 0.4)
82Closed deciduous needle-leaved forest (fc > 0.4)
91Open mixed-leaf forest (broad-leaved and needle-leaved)
92Closed mixed-leaf forest (broad-leaved and needle-leaved)
3Shrubland120Shrubland
121Evergreen shrubland
122deciduous shrubland
4Grassland130Grassland
5Wetlands180Wetlands
6Impervious surfaces190Impervious surfaces
7Bare areas140Lichens and mosses
150Sparse vegetation (fc < 0.15)
152Sparse shrubland (fc < 0.15)
153Sparse herbaceous (fc < 0.15)
200Bare areas
201Consolidated bare areas
202Unconsolidated bare areas
8Water body210Water body
9Permanent ice and snow220Permanent ice and snow
250Filled value
Table 2. LULC classification system in Myanmar.
Table 2. LULC classification system in Myanmar.
CodeLevel 1 ClassesLC IDLevel 2 Classes
1Cropland10Rainfed cropland
11Herbaceous cover
12Tree or shrub cover (orchard)
20Irrigated cropland
2Evergreen broad-leaved forest51Open evergreen broad-leaved forest
52Closed evergreen broad-leaved forest
3Deciduous broad-leaved forest61Open deciduous broad-leaved forest
(0.15 < fc < 0.4)
62Closed deciduous broad-leaved forest (fc > 0.4)
4Evergreen needle-leaved forest71Open evergreen needle-leaved forest
(0.15 < fc < 0.4)
72Closed evergreen needle-leaved forest (fc > 0.4)
5Shrubland120Shrubland
121Evergreen shrubland
122deciduous shrubland
6Grassland130Grassland
150Sparse vegetation (fc < 0.15)
200Bare areas
201Consolidated bare areas
202Unconsolidated bare areas
7Impervious surfaces190Impervious surfaces
8Wetlands & Water body180Wetlands
210Water body
220Permanent ice and snow
Table 3. Driving factors system of LUCC in Myanmar.
Table 3. Driving factors system of LUCC in Myanmar.
CategoryIndexunit
ClimateX1 Average annual temperature°C
X2 Total annual precipitationmm
Social developmentX3 Total population104
X4 Rural population104
X5 Urban population104
X6 Urbanization rate%
Economic developmentX7 Gross Domestic Product (GDP)USD billion (current USD)
X8 Agricultural value addedUSD billion (current USD)
X9 Industrial value addedUSD billion (current USD)
X10 Natural rubber production104 t
X11 Wood charcoal production104 t
X12 Wood production104 m3
X13 Vegetable production104 t
X14 Oil crop production104 t
X15 Food production104 t
X16 Fruit production104 t
Note: X1–X2 are from GLDAS-2.1 and PERSIANN-CDR datasets; X3–X9 are from the World Bank; X10–X16 are from the World Food and Agriculture Organization.
Table 4. Component matrix of the principal component analysis.
Table 4. Component matrix of the principal component analysis.
VariablesDescriptionComponent
F1-Socio-EconomicF2-ClimateF3-Agriculture
X1Average annual temperature0.637−0.356−0.425
X2Total annual precipitation−0.2030.8810.362
X3Total population 0.979 −0.1780.097
X4Rural population 0.986 −0.1440.040
X5Urban population 0.968 −0.2040.142
X6Urbanization rate0.964−0.2060.157
X7Gross Domestic Product (GDP)0.945−0.1080.307
X8Agricultural value added0.8840.4670.007
X9Industrial value added 0.969 0.1040.223
X10Natural rubber production0.949−0.1710.263
X11Wood charcoal production 0.966 0.1550.121
X12Wood production 0.839 −0.066−0.496
X13Vegetable production 0.956 0.132−0.130
X14Oil crop production0.8580.466−0.181
X15Food production0.6640.489 −0.563
X16Fruit production 0.989 −0.0770.111
Variance (%) 56.29%27.14%13.46%
Eigenvalues 9.004.342.15
Table 5. The relationship between the area of different LUC types and the three principal components.
Table 5. The relationship between the area of different LUC types and the three principal components.
Cropland Y 1 = 188 , 745.6 * * * + 4258.8 ×   F 1 * * + 1089.7 × F 3 * R2 = 0.99, p < 0.01
Evergreen broad-leaved forest Y 2 = 202 , 771.6 * * * 7311.9 ×   F 1 * R2 = 0.90, p < 0.05
Deciduous broad-leaved forest Y 3 = 193 , 633.8 * * * 8608.8 × F 1 * R2 = 0.84, p < 0.05
Shrubland Y 4 = 100 , 310.6 * * * + 6547.5 × F 1 * R2 = 0.82, p < 0.05
Grassland Y 5 = 4162.1 * * * + 546.7 × F 1 * * R2 = 0.97, p < 0.01
Impervious surfaces Y 6 = 3842.3 * * * + 722.5 × F 1 *   R2 = 0.79, p < 0.05
Note: * indicates passing the test of p = 0.05, ** indicates passing the test of p = 0.01, and *** indicates passing the test of p = 0.001.
Table 6. Results of multiple stepwise regression analysis of the area and driving factors for different LUC types.
Table 6. Results of multiple stepwise regression analysis of the area and driving factors for different LUC types.
Cropland Y 1 = 109 , 915.1 * * + 3037.2 × X 6 * * 2.8 × X 15 * R2 = 0.99, p < 0.01
Evergreen broad-leaved forest Y 2 = 335 , 789.6 * * * 4591.4 × X 6 * R2 = 0.91, p < 0.05
Deciduous broad-leaved forest Y 3 = 468 , 415.7 * * * 10454.4 × X 6 * * * R2 = 0.98, p < 0.001
Shrubland Y 4 = 37 , 508.0 * * * + 42.7 × X 5 * * * R2 = 0.99, p < 0.001
Grassland Y 5 = 4817.4 * * + 6.3 × X 9 * * 0.4 × X 15 * R2 = 0.99, p < 0.05
Impervious surfaces Y 6 = 9808.8 * * * + 2.7 × X 3 * * *   R2 = 0.99, p < 0.001
Note: * indicates passing the test of p = 0.05, ** indicates passing the test of p = 0.01, and *** indicates passing the test of p = 0.001.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, Y.; Hu, Y.; Niu, X.; Yan, H.; Zhen, L. Myanmar’s Land Cover Change and Its Driving Factors during 2000–2020. Int. J. Environ. Res. Public Health 2023, 20, 2409. https://doi.org/10.3390/ijerph20032409

AMA Style

Wang Y, Hu Y, Niu X, Yan H, Zhen L. Myanmar’s Land Cover Change and Its Driving Factors during 2000–2020. International Journal of Environmental Research and Public Health. 2023; 20(3):2409. https://doi.org/10.3390/ijerph20032409

Chicago/Turabian Style

Wang, Yiming, Yunfeng Hu, Xiaoyu Niu, Huimin Yan, and Lin Zhen. 2023. "Myanmar’s Land Cover Change and Its Driving Factors during 2000–2020" International Journal of Environmental Research and Public Health 20, no. 3: 2409. https://doi.org/10.3390/ijerph20032409

APA Style

Wang, Y., Hu, Y., Niu, X., Yan, H., & Zhen, L. (2023). Myanmar’s Land Cover Change and Its Driving Factors during 2000–2020. International Journal of Environmental Research and Public Health, 20(3), 2409. https://doi.org/10.3390/ijerph20032409

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop