Land Cover Changes from 1990 to 2019 in Papua, Indonesia: Results of the Remote Sensing Imagery

Long-term land cover changes play a significant driver of ecosystem and function of natural biodiversity. Hence, their analysis can be used for evaluating and supporting government plans, especially conservation and management of natural habitats such as sago palm. In Papua Province of Indonesia, sago palm has been stated as one of the priority plants in the Medium-Term Development Plan (R.P.J.M.). However, limited studies have examined this palm in one of the Regencies of Papua Province, namely, Merauke Regency. In this study, we performed remotely sensed data imagery and supervised classification to produce land cover maps from 1990 to 2019. During the study period, twenty-one land cover classes were identified. The six classes of the natural forest consist of primary dryland forest, secondary dryland forest, primary mangrove forest, secondary mangrove forest, primary swamp forest, and secondary swamp forest; thus, fifteen classes of non-forested area. Concerning the sago palm habitat, our study evaluated two different categories (1) based on the land cover scheme from the Ministry of Environment and Forestry and (2) according to the peatland land cover ecosystem in Papua. Based on paired samples t-test, the result indicated statistically significant changes specifically at primary dryland (p-value = 0.015), grassland (p-value = 0.002) and swamp (p-value = 0.007). Twelve from 20 districts of Merauke Regency tend to lose the forecasted natural habitat of the sago palm. Therefore, this study suggests the further need to recognize and estimate the yield of sago palm area in these various ecosystems.


Introduction
The land cover indicates the physical land class covered by swamp forests, mining areas, and other land cover classes. In contrast, land use refers to the purpose land serves, for example, recreation and wildlife habitat. Land cover and land use are often used reciprocally, but both of them can be performed to support various purposes, for instance, identification and change detection [1,2]. Land cover changes information is useful to achieve a better perspective of landscape dynamics and is also proper for evaluating the sustainability of natural resources [3,4]. Thus, ground cover monitoring and mapping are required to investigate spatial planning and environmental examination [5,6]. Additionally, land Thus, we calculated the rate of land cover change throughout the study. Regarding the sago palm habitat, this study examined the natural ecosystem of sago palm habitat based on two categories, i.e., (1) based on MoEF land cover scheme [44], and (2) peatland land cover in Papua [45]. We also applied statistical analysis to investigate the significant change of sago palm habitat based on their land cover. Our findings provide current land cover maps and the evaluation of land cover map changes over time, which can be used to support decision-making of the local Government for managing and conserving natural resources.

Study Area
The study was carried out in Merauke Regency (137 • 38 52.9692 E-141 • 0 13.3233 E and 6 • 27 50.1456 S-9 • 10 1.2253 S; Figure 1), which is located on the southern part of Papua Province. Papua Province is one of the thirty-four provinces of Indonesia, with a total area of about 31,509.162 ha and is noticed as a province with the most significant area. Papua and West Papua contribute to approximately 10% of the world mangrove area within various ecosystems [46]. In this study, the selected sago palm tree is scientifically known as Metroxylon sagu Rottb. The palm tree is one of the ecological tree species that grows in Indonesia, particularly in Papua and West Papua [21]. The palm has tremendous advantages to support the food sector, bioethanol [47,48], as a raw material for the agro-industry, and other aspects of sago [20]. This palm rises well in freshwater such as in Jayapura and swamps or salty areas of tropical lowlands, for example, in Merauke Regency.
Merauke Regency encompasses twenty districts, namely, Ulilin, Muting, Kaptel, Ngguti, Ilwayab, Tabonji, Waan, Kimaam, Tubang, Okaba, Malind, Kurik, Animha, Elikobel, Jagebob, Tanah miring, Semangga, Sota, Naukenjerai, and Merauke). The capital of Regency is located in Merauke. This Regency occupies an area of around 4,851,715 ha and is well known as the largest Regency in Papua Province. Merauke Regency is also a leader of the top three paddies providers, with about 91.47% of paddy production over Papua Province [49].The site Regency is bordered directly with Papua New Guinea and Australia, with an average temperature of around 27 • C, while the humidity is about 81%. The type of forest in this regency consists of primary forest and secondary forest, including mangrove, swamp, and dryland [49,50].
Sustainability 2020, 12, x FOR PEER REVIEW 3 of 20 sago palm habitat, this study examined the natural ecosystem of sago palm habitat based on two categories, i.e., (1) based on MoEF land cover scheme [44], and (2) peatland land cover in Papua [45].
We also applied statistical analysis to investigate the significant change of sago palm habitat based on their land cover. Our findings provide current land cover maps and the evaluation of land cover map changes over time, which can be used to support decision-making of the local Government for managing and conserving natural resources.

Study Area
The study was carried out in Merauke Regency (137°38'52.9692"E-141°0'13.3233"E and 6°27'50.1456"S-9°10'1.2253"S; Figure 1), which is located on the southern part of Papua Province. Papua Province is one of the thirty-four provinces of Indonesia, with a total area of about 31.509.162 ha and is noticed as a province with the most significant area. Papua and West Papua contribute to approximately 10% of the world mangrove area within various ecosystems [46]. In this study, the selected sago palm tree is scientifically known as Metroxylon sagu Rottb. The palm tree is one of the ecological tree species that grows in Indonesia, particularly in Papua and West Papua [21]. The palm has tremendous advantages to support the food sector, bioethanol [47,48], as a raw material for the agro-industry, and other aspects of sago [20]. This palm rises well in freshwater such as in Jayapura and swamps or salty areas of tropical lowlands, for example, in Merauke Regency.
Merauke Regency encompasses twenty districts, namely, Ulilin, Muting, Kaptel, Ngguti, Ilwayab, Tabonji, Waan, Kimaam, Tubang, Okaba, Malind, Kurik, Animha, Elikobel, Jagebob, Tanah miring, Semangga, Sota, Naukenjerai, and Merauke). The capital of Regency is located in Merauke. This Regency occupies an area of around 4,851,715 ha and is well known as the largest Regency in Papua Province. Merauke Regency is also a leader of the top three paddies providers, with about 91.47% of paddy production over Papua Province [49].The site Regency is bordered directly with Papua New Guinea and Australia, with an average temperature of around 27 °C, while the humidity is about 81%. The type of forest in this regency consists of primary forest and secondary forest, including mangrove, swamp, and dryland [49,50].

Data Acquisition and Preprocessing
In this study, two types of data from remotely sensed and secondary data were analyzed. Several secondary data such as the type of forest, area of forest by function, were contributed by Government Sustainability 2020, 12, 6623 4 of 18 agencies, i.e., Plantation and Forestry, and statistics agencies through their catalogues: Papua Province in figure 2020 (Number of catalogue: 1102002.94) [49], and Merauke in figure 2020 (Number of catalogue: 1102002.9401) [50]. Existing land cover maps for 1990, 1996, 2003, 2006, 2011, and 2014 were published on www.webgis.menlhk.go.id. The provincial boundary spatial data were acquired from the Regional Development Planning Agency or Badan Perencanaan Pembangunan Daerah (BAPPEDA) of Papua Province that we used as supporting data. Land cover classes of Indonesia and the description are referred to the Ministry of Environment and Forestry (MoEF), which includes the Standardization Agency of Indonesia or Badan Standardisasi Nasional (B.S.N.), specifically, the Indonesian National Standard or Standard Nasional Indonesia (S.N.I. 8033:2014). The land cover is classified into twenty-three classes that consist of six classes of forests, one plantation of the forest, 16 classes of non-forests (Table 1). To achieve remotely sensed data (  (2013). In this study, we applied Landsat 7 ETM + and Landsat 8 OLI/TIRS (L1T) multispectral images covering the study area. The L-1Terrain (L1T) product will automatically correct their geometrics and radiometric based on inputs from the sensors as well as the Ground Control Point (G.C.P.) and Digital Elevation Models (D.E.M.). We obtained data freely for the years 2000,2009,2015,2016,2017,2018, and 2019 (E.T.M. + and OLI/TIRS), with 30 m of the resolution, 705 km of altitude, and less than 50% cloud cover. To reduce this cloud cover, we combined multitemporal Landsat images from two seasons of Papua, which were appropriately selected from Figure 2.  All remote sensing data were downloaded from the archives of the United States Geological Survey (USGS), available on https://earthexplorer.usgs.gov/. We also recorded field data location using a hand-held global positioning system (G.P.S.).
After extracting the images, we carried out preprocessing steps consisting of radiometric calibration to remove the effect of atmospheric, illumination angles, and seasonal variation across the image [51]. Then, we projected the image according to our study area (European Petroleum Survey Group/EPSG: 23894), which is Universal Transverse Mercator (U.T.M.) World Geodetic System (W.G.S.) 84, in zone 54S local projection type. Preliminary image based on the area of interest was conducted using false color composites of red, green, and blue. To add these, we also did image enhancement and mosaicking to merge the image frame; all the preprocessing tasks here were carried out using the System for Automated Geoscientific Analyses (SAGA) Geographic Information System (GIS) that was included in the Quantum GIS (QGIS) software.
Survey Group/EPSG: 23894), which is Universal Transverse Mercator (U.T.M.) World Geodetic System (W.G.S.) 84, in zone 54S local projection type. Preliminary image based on the area of interest was conducted using false color composites of red, green, and blue. To add these, we also did image enhancement and mosaicking to merge the image frame; all the preprocessing tasks here were carried out using the System for Automated Geoscientific Analyses (SAGA) Geographic Information System (GIS) that was included in the Quantum GIS (QGIS) software.

Data Processing
For image classification, firstly, the actual land cover maps from MoEF (1990MoEF ( , 1996MoEF ( , 2003MoEF ( , 2006MoEF ( , 2011MoEF ( , 2014, and 2017) were clipped and overlaid on Landsat imagery in the same year. During this process, we looked at the entire map and made corrections wherever it was needed. Afterwards, supervised classification was applied [52]; training samples were selected by delineating polygons at characteristics sites. We chose 15 of each class as the training data; at this step of the procedure, we also compared each class that was collected in the fields using G.P.S. The accuracy assessments of the image classification were done as an integral part of the image classification process using QGIS. Next, we developed land cover maps for other years (2000,2009,2015,2016,2018, 2019) by analyzing the supervised classification result and the existing land cover maps. Since the visual interpretation should be more standardized, we therefore complied with the Ministry of Environment and Forestry scheme [44], which included the Indonesian National Standard (S.N.I.), particularly SNI 8033:2014, that allows us to recognize the image through structure, texture, shape, pattern and color; at the same time, grants land cover classification in Indonesia (Table 1). According to that scheme, there were twenty-three land cover classes in total, consisting of 6 classes of natural forests, 1 class of plantation forest, 16 classes of non-forest. Following the natural habitat of sago palm, there were two

Data Processing
For image classification, firstly, the actual land cover maps from MoEF (1990MoEF ( , 1996MoEF ( , 2003MoEF ( , 2006MoEF ( , 2011MoEF ( , 2014, and 2017) were clipped and overlaid on Landsat imagery in the same year. During this process, we looked at the entire map and made corrections wherever it was needed. Afterwards, supervised classification was applied [52]; training samples were selected by delineating polygons at characteristics sites. We chose 15 of each class as the training data; at this step of the procedure, we also compared each class that was collected in the fields using G.P.S. The accuracy assessments of the image classification were done as an integral part of the image classification process using QGIS. Next, we developed land cover maps for other years (2000,2009,2015,2016,2018, 2019) by analyzing the supervised classification result and the existing land cover maps. Since the visual interpretation should be more standardized, we therefore complied with the Ministry of Environment and Forestry scheme [44], which included the Indonesian National Standard (S.N.I.), particularly SNI 8033:2014, that allows us to recognize the image through structure, texture, shape, pattern and color; at the same time, grants land cover classification in Indonesia (Table 1). According to that scheme, there were twenty-three land cover classes in total, consisting of 6 classes of natural forests, 1 class of plantation forest, 16 classes of non-forest. Following the natural habitat of sago palm, there were two sites, i.e., primary swamp forest and secondary swamp forest. At the final step, we validated the land cover maps by using a geographic browser, which is Google Earth Pro [53], that provides a higher resolution of satellite imagery. We ensured that checkpoints were spread throughout the study area and expressed all land cover classes; nonetheless, these numbers were selected in a different amount due to the possibility of high resolution images.

Data Analysis
The land area for general characteristics of the study site were presented in mean and standard deviation (minimum, maximum). Land area changes over time were analyzed using the paired t-test [54,55] for the years 1990 and 2019. To designate the statistical significance in all analyses, a p-value of less than 0.05 was used. Statistical analysis was performed using IBM SPSS statistics version 25 (IBM Corp., Armonk, NY, USA). In this study, a t-test was performed to examine the means of land cover for the general characteristics of the site and whether or not the natural habitat of sago palm has significantly changed throughout the study.
Moreover, we also calculated gain and losses [56], the rate of land cover change, in the following way (Equation (1)) [57,58].
Area(f ) and Area(i) are the areas of a certain land type at the final area or at the end, while Area(i) represents an initial or at the beginning of the research period, respectively. The negative rate number of the land cover shows a decreasing trend during the period of study, while in contrast, the positive number indicates an increasing area of each class or category.

Land Cover Changes in Merauke Regency
Recent land cover maps from 1990 to 2019 are presented in Figure 3; we made 13 years of land cover maps in this Regency, consisting of land cover maps in 1990,1996,2000,2003,2006,2009,2011, and from 2014 to 2019. Twenty-one land cover categories were identified: (1) 6 classes of the natural forest included primary dryland forest, secondary dryland forest, primary mangrove forest, primary swamp forest, secondary mangroves forest, secondary swamp forest; (2) 15 classes of non-forest consisting of swamp shrub, swamp, bush/shrub, estate crop plantation, settlement area, barren land, clouds, grassland, water body, dryland agriculture, shrub-mixed dryland, paddy field, fishpond, airport, transmigration area ( Figure 3). By using supervised classification [1] of remotely sensed imagery, it is possible to get the resultant area estimates that occurred in this Regency, as shown by Tables 3-5.     In 1990, Merauke Regency was covered by the natural forest of around 2,440,396 ha, or approximately 50.3% (Figure 4), compared to the non-forest area of about 2,411,319 ha, or 49.70% of total area (Table 3). Although this Regency was dominated by natural forest, the area cropped highest by swamp shrub and was followed by three classes of natural forest area, namely, primary dryland, secondary dryland and secondary swamp forest, around 19.17%, 14.32%, 13.15%, and 10.95%, with areas of 930,069 ha, 694,737 ha, 638,049 ha, and 531,109 ha, respectively. Also in 1990, the land cover map did not include estate crop plantation and fishpond. Conversely in 2003, this site was mostly covered by non-forest area for about 2,780,824 ha (Table 3), or approximately 57.32% of the total area, while the natural forest had an area of 2,070,891 ha, or about 42.68%.

Land Cover Losses and Gains from 1990 to 2019
Land cover category losses and gains were also examined in this study. The total changed areas for each class are shown in Table6. The result indicates that five classes of natural forests were generally decreasing; only one class of it, that has a small number, was increasing. Primary swamp forest was reduced by around 40.95%, followed by primary dryland forest, that also decreased by about 27.98% during this period. Likewise, primary mangrove, secondary mangrove, and thus, the secondary dryland forest were also reduced by around 6.39%, 5.07%, and 1.06% respectively, while only secondary swamp forest had gained not up to 0.03% in terms of the total area and the percentage of change. In contrast, in the non-forest category, the paddy field was the significantly increased area with anoverall percentage of change of about 316.26%. Thus, an expansion of about 124.38% in the settlement area was observed. On the other hand, we could not examine the percentage of the changed estate crop plantation and fishpond due to the lack of initial value in 1990. Furthermore, the transmigration area and water body lost their aerial extent around 30.20% and 0.61%. Other non-forest categories gained less than 80%, i.e., bush/shrub, shrub-mixed dryland, dryland agriculture, grassland swamp, airport, barren land, and swamp shrub subsequently around 79.95%, 57.84%, 50.73%, 17.72%, 10.94%, 10.27%, 8.85%, and 1.39%, respectively. Although the transmigration and water body lost their areas, the overall non-forest gained an increasing trend in the 1990-2019 change period. In 2006, the Regency was also coated by a non-forested area of 2,828,047 ha, or 58.29%, while the natural forest was decreased to 41.71% (Table 4). The most typical land cover of this Regency was swamp shrub with the area of 949,786 ha, followed by two other land cover types, i.e., grassland about 655,175 ha, and secondary dryland forest with the area around 627,494 ha. The pattern of land cover changes continued in 2014; this Regency was also mostly covered by the non-forest area, which is swamp shrub, grassland, and another class of natural forest, namely secondary dryland forest. The total area change in 2014 was about 906,111 ha, 708,703 ha, and 678,803 ha, respectively.
In 2015 (Table 5), it is shown that the natural forests were 1,977,080 ha, or approximately 40.75%, of the total area. Thus, the non-forest occupied 2,874,635 ha or near to 59.25%. This percentage explained why, in 2015, the non-forest category revealed an increasing trend which was about 18.5% higher than the natural forest. The area was still wrapped highest by swamp shrub for about 860,813 ha and was followed by the secondary dryland forest (664,888 ha). The grassland seems to decrease compared to the previous year. In 2018, the non-forest area reached the highest percentage of change: a loss of about 60.54% more than the natural forest, which tended to lose space for around 39.46% of the total Regency area. The dominant land cover in 2018 was swamp shrub, secondary dryland forest, and then grassland. Even though in 2019 the Regency was still covered by non-forested area (for about 2,767,158 ha or 57.03% of the entire site), the area of natural forest slightly increased about 3.51% compared to 2018 (Figure 4).

Land Cover Losses and Gains from 1990 to 2019
Land cover category losses and gains were also examined in this study. The total changed areas for each class are shown in Table 6. The result indicates that five classes of natural forests were generally decreasing; only one class of it, that has a small number, was increasing. Primary swamp forest was reduced by around 40.95%, followed by primary dryland forest, that also decreased by about 27.98% during this period. Likewise, primary mangrove, secondary mangrove, and thus, the secondary dryland forest were also reduced by around 6.39%, 5.07%, and 1.06% respectively, while only secondary swamp forest had gained not up to 0.03% in terms of the total area and the percentage of change. In contrast, in the non-forest category, the paddy field was the significantly increased area with an overall percentage of change of about 316.26%. Thus, an expansion of about 124.38% in the settlement area was observed. On the other hand, we could not examine the percentage of the changed estate crop plantation and fishpond due to the lack of initial value in 1990. Furthermore, the transmigration area and water body lost their aerial extent around 30.20% and 0.61%. Other non-forest categories gained less than 80%, i.e., bush/shrub, shrub-mixed dryland, dryland agriculture, grassland swamp, airport, barren land, and swamp shrub subsequently around 79.95%, 57.84%, 50.73%, 17.72%, 10.94%, 10.27%, 8.85%, and 1.39%, respectively. Although the transmigration and water body lost their areas, the overall non-forest gained an increasing trend in the 1990-2019 change period.

Land Cover Changes and Their Impact on the Prediction of the Natural Habitat of the Sago Palm
We evaluated the natural habitat of sago palm based on the land cover changes of the Regency. We used two different categories, i.e., (1) based on the Ministry of Environment and Forestry land cover schemes (Table 1) and (2) the peatland land cover ecosystem of Papua [45]. In [44] (Table 1), as mentioned previously, it predicts two typical habitats of the sago palm, namely, primary swamp forest and secondary swamp forest. Other natural ecosystems were forecasted, i.e., dryland forest, bush/shrub, swamp, swamp shrub, and savanna/grassland [45]; the local community usually utilizes this peat ecosystem to find fish, or to plant and harvest sago palm [45]. Here, we assessed statistically twenty districts of Merauke Regency, by including eight classes of the potential habitat of the sago palm(primary dryland forest, secondary dryland forest, primary swamp forest, secondary swamp forest, bush/shrub, grassland, swamp shrub, swamp), as shown by Tables 7 and 8 using the paired samples t-test for the years 1990 and 2019. We used this result to examine the significant changes of the various natural sago palm ecosystems in Merauke Regency as well as to support our third aim of this study. Table 7 presents that, in 1990, Semangga, Malind, Kurik, and Naukenjerai did not have a possible area of sago of more than 11 ha, while Waan, Ulilin, Kimaam, and Ngguti had a larger potential area of sago habitat compared to other districts, with the mean value of the land area of 60,482.52 ha, 56,469.20 ha, 45,414.98 ha, and 39,547.46 ha, respectively. In the year 2019, it is observed that those districts still had a large area in comparison to other districts based on their mean values, i.e., Waan, Ulilin, Kimaam, and Muting. However, Ulilin, Kimaam, and Muting gradually decreased the land area by around 1.69 ha to 451.68 ha throughout the study. In contrast, Waan district extended the area by approximately 899.53 ha. Semangga, Malind, Kurik, and Naukenjerai still had an inadequate site according to their mean value, and overall, the area progressively declined during the study period. Table 7. General characteristics of the prediction of sago palm habitat in Merauke Regency (N = 8).

District 1990 2019
Statistical analyses using paired t-test were carried out on data generated from supervised classification-based restricted by remote sensing images of the study area. Specifically, the potential habitat of the sago palm is shown in Table 8. The result of the paired samples t-test (Table 8) demonstrates that the p-value is less than 0.05 for primary dryland (p-value = 0.015), grassland (p-value = 0.002) and swamp (p-value = 0.007). Other categories are greater than 0.05, namely, secondary dryland, primary swamp forest, secondary swamp forest, bush/shrub, and swamp shrub. Therefore, the result reveals that primary dryland forest, grassland, and swamp experienced a statistically significant change from 1990 to 2019.

Discussion
We used the Landsat satellite imagery and supervised classification to develop recent land cover changes in Merauke Regency over the past 29-year period ( Figure 3). The land cover of this Regency today consists of six classes of the natural forest and fifteen classes of non-forested areas. The result from this study showed that a land cover map of Merauke Regency is now generally dominated by non-forested regions, followed by natural forests. Concerning the proportion of natural versus non-forest regions (Figure 4), it is shown that, in 1990, the percentage of the natural forested landscape was 49.70%; however, in 2019, it declined by 7.33%. In the opposite category, the non-forested area was about 49.70% in 1990 and then expanded in 2019, approximately by 57.03%. During this study period, we found two classes of non-forest land cover that extended to up to 100% in terms of total area and the percentage of change, i.e., settlement area and paddy field; 124%, 316.26%, respectively ( Table 6). The land area changes also explain why Merauke Regency nowadays is the largest producer of paddy commodities by around 91.47% over Papua Province [49,50]. The settlement area was naturally expanding, expressing population growth in this Regency. Inversely, two other classes of non-forest, for instance water body and transmigration area, had slightly decreased around 0.61% and 30.20% ( Table 6). The increment of built-up areas can bea result of migration or urban growth and may lead to environmental problems such as degradation in vegetation and water bodies [59].
In this study, the natural forest cover consists of primary dryland, secondary dryland, primary mangrove, secondary mangrove, primary swamp, and secondary swamp forest. Our evidence confirmed that primary swamp forest and primary dryland forest had higher area losses compared to other forest types, such as primary mangrove, or secondary mangrove. As shown in Table 6, the primary dryland forest lost about 27.98%. Meanwhile, secondary dryland forest, primary mangrove forest, and secondary mangrove forest lost around 1.06%, 6.39%, and 5.07%, respectively. Only the secondary swamp forest gained about 0.03% during the study period. Decrement of the forested landscape was also experienced by some European countries due to increased temperatures caused by climate change. In the Czech Republic, the forest area is currently diminished because of the bark beetle attack [60,61]. Most of the Czech forests are productive forests with less-diverse trees [62]. Therefore, to restore the forested landscape, it is suggested to cultivate various plants with different ages of trees and appropriate width stands of the trees in the Czech forests to provide a natural interconnection in the food chain or food web, to reduce the damage caused by insect attacks [63,64]. The challenge can be taken as a lesson learned, and the recommendations can also be adopted for forest management and restoration in Merauke Regency.
In contrast, the agricultural sector, such as dryland agriculture, mix dryland farms, gained from around 50% to 57% of the land cover area in Merauke Regency (Table 6), which is in agreement with the previous notion that agricultural development can be a threat to the sustainability of mangrove forest in Papua [46].Nevertheless, the mangrove environment can support ecological and economic services to neighborhood, society, and also in the country [65]. Another class of natural forest that suffered in terms of the percentage of change was primary swamp forest; this class lost area by around 40.95% during the study period. This phenomenon can lead directly to a decreased forest quality such as loss of plants that live ecologically in this area, for instance, sago palm (Metroxylon sago Rottb) ( Table 1) [45]. It also refers to previous studies that found sago palm raising well in particular ecosystems, for example, swamp forest, peat soil, saline [19].
Moreover, sago forests have the opportunity of excellent carbon sinks for carbon absorption, reducing the greenhouse impact and global warming [66]. The wetlands changes and the effect on biodiversity are also acquainted by East African countries such as Tanzania [12]; the number of native species, for instance, woody plant species, bird species, have declined due to some reasons such as deforestation, agricultural expansion, and the extension of the settlement area. Deforestation in peat swamp areas is also experienced by Malaysia [67]; it mentioned that swamp forests declined mostly as the result of the transformation to other non-forested regions such as the agriculture sector. In Papua, sago is believed to be one of the natural resources which can enhance the needs of communities such as to support households' low energy consumption [47] and to strengthen food security as a staple food [68]. Therefore, more preventive approaches, such as regulation or desertification programs from the local government or relevant stakeholders, are highly required.
Other natural ecosystems of sago palm, as estimated in the peatland ecosystem [45], were dryland, bush/shrub, grassland, swamp shrub, and swamp. We analyzed these areas in twenty districts of Merauke Regency by using paired t-test (Tables 7 and 8 To add to this, we analyzed the natural habitat of sago palm in eight possible ecosystems from land cover changes in 20 districts of Merauke Regency, especially primary dryland, secondary dryland, primary swamp forest, secondary swamp forest, bush/shrub, grassland, swam shrub, and swamp. Based on land cover losses and gains (Table 6), the peatland ecosystem gained in the area of bush/shrub (+79.95%), swamp shrub (+1.39%), swamp (+10.94%), savanna/grassland (+17.72%), and secondary swamp forest (+0.03%). However, this regency lost the peatland area mostly at primary swamp (−40.95%), primary dryland (−27.98%), and secondary dryland (−1.06%) from 1990 to 2019. Nevertheless, the decreasing trend of the natural forest seems to be broke out in 2019. It can be seen in Figure 4 that the area of natural forest has increased by around 3.5% compared to the previous year. This circumstance could be achieved through support by some strategic plans from the Indonesian Government such as the Peat Restoration Agency or Badan Restorasi Gambut (B.R.G.), the mitigation of forest, land fires, and forest evaluation. The Peat Restoration Agency takes a responsibility to restore peat lands ecosystem in several provinces including Papua [16,17]. The forest restoration is expected to be continued in the coming years. In this respect, the annual land cover maps on this site are required by the region, as well as on a national scale, to monitor and to evaluate particular habitat; since the habitat is damaged overtime, Indonesia will suffer from the loss of various natural plants [43,69]. Moreover, the outcomes of the paired sample t-test as presented in Table 8 reports that primary dryland (p-value = 0.015), grassland (p-value = 0.002) and swamp (p-value = 0.007) significantly changed statistically. The area of primary dryland seems to have decreased; however, swamp and grassland have increased in the statistical term.
On the other hand, this study has not compared the trend of sago palm over time, due to the lack of supporting data, especially the area of sago palm in this site. Therefore, the yield estimation or calculation of this plant in this particular area should be considered annually produced. This study recommends further research in terms of detecting and recognizing the sago palm in this specific area of Merauke Regency.

Conclusions
The study has produced recent land cover maps and examined the spatial temporal patterns and rate of land cover changes in Merauke Regency from 1990 to 2019 using remote sensing techniques and supervised classification. Merauke land cover consists of twenty-one classes, which are classified into six classes of natural forest and fifteen classes of non-forested areas. Our analysis of the land cover map presented that the largest declines mostly occurred in the natural forest, namely primary dryland forest, secondary dryland forest, primary mangrove, secondary mangrove, primary swamp forest, for about −81.45% altogether, and only secondary swamp forest has slightly increased by about 0.03% over time. According to the natural habitat of the sago palm, we evaluated eight possible ecosystems, namely, dryland, bush/shrub, grassland, swamp shrub, swamp, and swamp forest, using paired sample t-test. The result indicated statistically significant changes specifically at primary dryland (p-value = 0.015), grassland (p-value = 0.002), and swamp (p-value = 0.007). We also analyzed these areas in twenty districts of Merauke Regency. Our findings confirmed that 12 districts from 20 districts of Merauke Regency tend to lose the natural habitat of the sago palm, while only one district remains unchanged. Nonetheless, these particular ecosystems are beneficial to support local community life; for example, planting and harvesting sago palm. Therefore, these outcomes could be integrated within decision-makers and stakeholders to evaluate and to establish government development plans.