Next Article in Journal
Application of the DEA Double Bootstrap to Analyze Efficiency in Galician Sheltered Workshops
Previous Article in Journal
How does the Traditional Heavy Industry Use Ecotechnology to Achieve the Ecological Innovation Goal? Analysis Based on Expert Survey in China’s Shipbuilding Industry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Sri Murniani Angelina Letsoin
1,2,*,
David Herak
1,
Fajar Rahmawan
3 and
Ratna Chrismiari Purwestri
4,5,*
1
Department of Mechanical Engineering, Faculty of Engineering, Czech University of Life Sciences Prague, Kamycka 129, 16500 Praha-Suchdol, Czech Republic
2
Faculty of Engineering, University of Musamus, Merauke Regency, Papua 99611, Indonesia
3
INTSIA Foundation of Papua Province, Furia 3 Number 116 Abepura, Jayapura City, Papua 99225, Indonesia
4
Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, Kamycka 129, 16500 Praha-Suchdol, Czech Republic
5
Institute of Nutritional Sciences, University of Hohenheim, Garbenstrasse 28, 70599 Stuttgart, Germany
*
Authors to whom correspondence should be addressed.
Sustainability 2020, 12(16), 6623; https://0-doi-org.brum.beds.ac.uk/10.3390/su12166623
Submission received: 16 July 2020 / Revised: 10 August 2020 / Accepted: 11 August 2020 / Published: 16 August 2020

Abstract

:
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.

1. 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 cover and land use analysis will help in the reliable prediction of future circumstances. For example, future changes in forest cover can be predicted using the substance gained from historical datasets and remote sensing observations [7].
Numerous studies have been conducted on the measurement of land cover change at a national as well as a global scale because of its tremendous impact on various aspects such as urban development, water supplies, or environmental studies [8,9]. Long-term monitoring of land cover changes is significantly needed to estimate carbon stock, ecosystem service, and biodiversity losses [10,11,12]. Generally, the main focus on developing countries such as Indonesia is related to the conversion of forest areas for supporting agriculture expansion, urban or infrastructure planning [13,14,15]. However, the loss of forest cover changes can rigidly decrease the natural sustainability of the site, including the loss of biodiversity, which leads to a reduction in ecosystem functioning and an increase in the climate change factor. To handle this, some strategic ways have been set up, for instance, the Indonesian Government in Presidential Instruction Number 3 of 2020 has mentioned about the mitigation of forest and land fires prevention, monitoring, and forest evaluation. Also, Presidential Decree Number 1 of 2016 established The Peat Restoration Agency or Badan Restorasi Gambut (B.R.G.) to coordinate and facilitate efforts to restore peat lands in several provinces, including in Papua [16,17].
In Papua, forest areas play an essential role as a natural habitat of sago palm that grows in the moist upland rainforest, freshwater, peat lands, swamps, or salty areas of tropical lowlands up to 700 m above sea level. The favorable preconditions date temperatures of above 25 °C and relative air humidity of around 70%; at the time of vegetative just before flowering, the plant transforms its saved nutrition towards starch, which fills the trunk. Thus, at the mature phase, it occupies a huge trunk and may reach a height of 6–10 m. Sago palm reaches commercial maturity at 9–12 years of age, when fruits start to develop and starch growth in the trunk reaches its highest level [18]. In comparison to other starches, sago produces an amount of approximately 400kg dry starch for each tree, while cassava or potato produces just around 40 kg [19]. Sago palm is also a part of indigenous costumes and, as a staple food, it will be cooked traditionally using burnt stones while the sago waste is used as livestock feed. Today, sago palm is said to be one of the priority plants in the Medium-Term Development Plan, or Rencana Pembangunan Jangka Menengah (R.P.J.M.), of Papua Province, because it provides many benefits in various sectors such as food security, the agro-industry, and environmental issues [20,21].
At the moment, current satellite imagery has been explored extensively for mapping and monitoring land cover changes using Landsat [22], Sentinel [23], Moderate Resolution Imaging Spectroradiometer (MODIS) [24], Synthetic Aperture Radar (SAR) [25], Satellite Pour I ‘Observation de la Terre (SPOT) [26], Pleiades [27], and other openness of data satellites. Furthermore, a broad range of suitable spectral bands with a very high resolution and accessibility with various computer-aided software has improved the use of remotely sensed data. For instance, the System for Automated Geoscientific Analyses (SAGA) [28], Quantum Geographic Information System (QGIS) [29], Sentinel Application Platform (SNAP) [30], eCognition [31], Earth Resources Data Analysis System (ERDAS) [32], and other platforms such as Google Earth [33] have improved the use of remotely sensed data. The long-term series of Landsat has been demonstrated successively to monitor land cover changes because of land degradation [34], to investigate land cover changes as a result of some variables such as urban expansion [35], flooding [36], deforestation [37], coastal abrasion [38], and vegetation spreads [39]. It continues to measure forest carbon stock [40], climate [41], biodiversity [42], and other land cover changes.
Although some studies have focused on monitoring land cover changes of Indonesia in general, we found studies applying remote sensing to evaluate changes in swamp forest land cover, and its impact on the natural habitat in this location are still unreported [43]. Therefore, this paper aims (1) to provide the current land cover maps for the Merauke Regency in Papua Province of Indonesia, (2) to evaluate land cover changes, and (3) to estimate the potential area forthe natural habitat ofthe sago palm. To deal with these objectives, we combined existing land cover maps that we obtained from the Ministry of Environment and Forestry (MoEF) and Landsat imagery to create new land cover maps. 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.

2. Materials and Methods

2.1. 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].

2.2. 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 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 (Table 2), we used Landsat multispectral images that, since 1972, have provided sensor Multispectral Scanner Systems (MSS)/Return Beam Vidicon (RBV) (1972/1978), MSS/Thematic Mapper (TM) (1982/1984), Enhanced Thematic Mapper (ETM+) (1993/1999) and Operational Land Imager/Thermal Infrared Scanner (OLI/TIRS) (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.

2.3. Data Processing

For image classification, firstly, the actual land cover maps from MoEF (1990, 1996, 2003, 2006, 2011, 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.

2.4. 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].
t h e   c h a n g e   r a t e   o f   l a n d   ( C ) = A r e a ( f ) A r e a ( i ) A r e a ( i )
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.

3. Results

3.1. 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 Table 3, Table 4 and Table 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%.
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).

3.2. 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.

3.3. 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 Table 7 and Table 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.
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.

4. 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 (Table 7 and Table 8). Based on the mean values, our findings confirmed that 12 districts from 20 districts of Merauke Regency lost the natural habitat of sago palm, namely, Elikobel (mean = 17,495.36), Jagebob (mean = 15,082.23), Kurik (mean = 7971.55), Kaptel (mean = 27,699.84), Kimaam (mean = 45,413.29), Malind (mean = 4883.73), Muting (mean = 39,489.52), Naukenjerai (mean = 9737.45), Ngguti (mean = 38,816.99), Tanah Miring (mean = 16,359.38), Tabonji (mean = 33,030.82), and Ulilin (mean = 56,017.52) (Table 7). Only one district remained unchanged, i.e., Semangga, thus, seven districts increased slightly, specifically Animha (mean = 16,983.27), Ilwayab (mean = 21,900.72), Merauke (mean = 15,126.20), Okaba (mean = 18,900.30), Sota (mean = 30,931.49), Tubang (mean = 32,655.19) and Waan (mean = 61,382.05). Waan had the largest potential of sago palm area in Merauke Regency.
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.

5. 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.

Author Contributions

Conceptualisation, methodology, validation, S.M.A.L., R.C.P., and D.H.; software, S.M.A.L., F.R., and R.C.P.; formal analysis, S.M.A.L. and R.C.P.; investigation, S.M.A.L. and F.R.; resources, D.H.; writing—original draft preparation, S.M.A.L., R.C.P., and D.H.; writing—review and editing, S.M.A.L., R.C.P., and D.H.; visualization, S.M.A.L. and F.R.; supervision, D.H. and R.C.P.; project administration, D.H.; funding acquisition, D.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study is part of on-going research being undertaken for PhD work of the first author. The first author is deeply grateful to the Indonesia Endowment Fund for Education (LPDP-Indonesia) for funding and supporting her PhD study. R.C.P. is also supported by grant EVA 4.0 No. CZ.02.1.01/0.0/0.0/16_019/0000803, financed by the OPRDE-Ministry of Education of the Czech Republic.

Acknowledgments

All the authors would like to thank the Ministry of Environment and Forestry (MoEF) of Indonesiafor granting existing land cover maps of Merauke Regency as well as the boundary map. We are also thankful tothe Provincial Development Planning Agency of Papua for providing an additional provincial boundary map. We appreciate the assistance of Meryem Tahri, during idea development, Julia V. Latul and Bambang A. Nugroho (Forestry Office of Merauke Regency) for sharing the peatland LULC map data, and Jitka Kumhalova, for peer reviewing the earlier draft of the manuscript. We are also grateful to our anonymous reviewers, editor, whose comments and suggestions improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, Y.; Fang, G.; Xu, Y.-P.; Tian, X.; Xie, J. Identifying how future climate and land use/cover changes impact streamflow in Xinanjiang Basin, East China. Sci. Total Environ. 2020, 710, 136275. [Google Scholar] [CrossRef]
  2. Halmy, M.W.A.; Gessler, P.E.; Hicke, J.A.; Salem, B.B. Land use/land cover change detection and prediction in the north-western coastal desert of Egypt using Markov-CA. Appl. Geogr. 2015, 63, 101–112. [Google Scholar] [CrossRef]
  3. Whittle, M.; Quegan, S.; Uryu, Y.; Stüewe, M.; Yulianto, K. Detection of tropical deforestation using ALOS-Palsar: A Sumatran case study. Remote. Sens. Environ. 2012, 124, 83–98. [Google Scholar] [CrossRef]
  4. Aliani, H.; Malmir, M.; Sourodi, M.; Kafaky, S.B. Change detection and prediction of urban land use changes by CA–Markov model (case study: Talesh County). Environ. Earth Sci. 2019, 78, 546. [Google Scholar] [CrossRef]
  5. Cheng, K.; Wang, J. Forest type classification based on integrated spectral-spatial-temporal features and random forest algorithm—A case study in the qinling mountains. Forests 2019, 10, 559. [Google Scholar] [CrossRef]
  6. Tripathy, P.; Kumar, A. Monitoring and modelling spatio-temporal urban growth of Delhi using cellular automata and geoinformatics. Cities 2019, 90, 52–63. [Google Scholar] [CrossRef]
  7. Hamad, R.; Balzter, H.; Kolo, K. Predicting land use/land cover changes using a CA-markov model under two different scenarios. Sustainability 2018, 10, 3421. [Google Scholar] [CrossRef] [Green Version]
  8. Anand, V.; Oinam, B. Future land use land cover prediction with special emphasis on urbanization and wetlands. Remote Sens. Lett. 2020, 11, 225–234. [Google Scholar] [CrossRef]
  9. Wu, Z.; Dai, E.; Wu, Z.; Lin, M. Future forest dynamics under climate change, land use change, and harvest in subtropical forests in Southern China. Landsc. Ecol. 2019, 34, 843–863. [Google Scholar] [CrossRef]
  10. Saragi-Sasmito, M.F.; Murdiyarso, D.; June, T.; Sasmito, S.D. Carbon stocks, emissions, and aboveground productivity in restored secondary tropical peat swamp forests. Mitig. Adapt. Strateg. Glob. Chang. 2019, 24, 521–533. [Google Scholar] [CrossRef] [Green Version]
  11. Zhao, M.; He, Z. Evaluation of the effects of land cover change on ecosystem service values in the upper reaches of the heihe river basin, northwestern China. Sustainability 2018, 10, 4700. [Google Scholar] [CrossRef] [Green Version]
  12. Seki, H.A.; Shirima, D.D.; Courtney Mustaphi, C.J.; Marchant, R.; Munishi, P.K.T. The impact of land use and land cover change on biodiversity within and adjacent to Kibasira Swamp in Kilombero Valley, Tanzania. Afr. J. Ecol. 2018, 56, 518–527. [Google Scholar] [CrossRef]
  13. Obidzinski, K.; Takahashi, I.; Dermawan, A.; Komarudin, H.; Andrianto, A. Can large scale land acquisition for agro-development in Indonesia be managed sustainably? Land Use Policy 2013, 30, 952–965. [Google Scholar] [CrossRef]
  14. Acosta, P.; Curt, M.D. Understanding the expansion of oil palm cultivation: A case-study in Papua. J. Clean. Prod. 2019, 219, 199–216. [Google Scholar] [CrossRef]
  15. Sloan, S.; Campbell, M.J.; Alamgir, M.; Engert, J.; Ishida, F.Y.; Senn, N.; Huther, J.; Laurance, W.F. Hidden challenges for conservation and development along the Trans-Papuan economic corridor. Environ. Sci. Policy 2019, 92, 98–106. [Google Scholar] [CrossRef]
  16. Cabinet Secretariat of the Republic of Indonesia. Available online: https://setkab.go.id/en/president-jokowi-establishes-peat-land-restoration-agency-brg/ (accessed on 9 June 2020).
  17. Cabinet Secretariat of the Republic of Indonesia. Available online: https://setkab.go.id/en/president-jokowi-signs-presidential-instruction-on-forest-land-fires-mitigation/ (accessed on 8 June 2020).
  18. Ehara, H.; Toyoda, Y.; Johnson, D.V. (Eds.) Sago Palm: Multiple Contributions to Food Security and Sustainable Livelihoods; Springer: Singapore, 2018; ISBN 978-981-10-5268-2. [Google Scholar]
  19. Zhu, F. Recent advances in modifications and applications of sago starch. Food Hydrocoll. 2019, 96, 412–423. [Google Scholar] [CrossRef]
  20. Singhal, R.S.; Kennedy, J.F.; Gopalakrishnan, S.M.; Kaczmarek, A.; Knill, C.J.; Akmar, P.F. Industrial production, processing, and utilization of sago palm-derived products. Carbohydr. Polym. 2008, 72, 1–20. [Google Scholar] [CrossRef]
  21. Karim, A.A.; Tie, A.P.-L.; Manan, D.M.A.; Zaidul, I.S.M. Starch from the Sago (Metroxylonsagu) palm treeproperties, prospects, and challenges as a new industrial source for food and other uses. Compr. Rev. Food Sci. Food Saf. 2008, 7, 215–228. [Google Scholar] [CrossRef]
  22. Yan, E.; Zhao, Y.; Lin, H.; Wang, G.; Mo, D. Improving the estimation of forest carbon density in mountainous regions using topographic correction and landsat 8 images. Remote Sens. 2019, 11, 2619. [Google Scholar] [CrossRef] [Green Version]
  23. Poortinga, A.; Tenneson, K.; Shapiro, A.; Nquyen, Q.; San Aung, K.; Chishtie, F.; Saah, D. Mapping plantations in myanmar by fusing landsat-8, sentinel-2 and sentinel-1 data along with systematic error quantification. Remote Sens. 2019, 11, 831. [Google Scholar] [CrossRef] [Green Version]
  24. Xing, L.; Niu, Z. Mapping and analyzing China’s wetlands using MODIS time series data. Wetl. Ecol. Manag. 2019, 27, 693–710. [Google Scholar] [CrossRef]
  25. Richard Allen, T.; Wang, Y.; Gore, B. Coastal wetland mapping combining multi-date SAR and LiDAR. Geocarto Int. 2013, 28, 616–631. [Google Scholar] [CrossRef]
  26. Flores-Casas, R.; Ortega-Huerta, M.A. Modelling land cover changes in the tropical dry forest surrounding the Chamela-Cuixmala biosphere reserve, Mexico. Int. J. Remote Sens. 2019, 40, 6948–6974. [Google Scholar] [CrossRef]
  27. Wang, D.; Wan, B.; Qiu, P.; Su, Y.; Guo, Q.; Wu, X. Artificial mangrove species mapping using pléiades-1: An evaluation of pixel-based and object-based classifications with selected machine learning algorithms. Remote Sens. 2018, 10, 294. [Google Scholar] [CrossRef] [Green Version]
  28. Cando-Jácome, M.; Martínez-Graña, A. Determination of primary and secondary lahar flow paths of the fuego volcano (Guatemala) using morphometric parameters. Remote Sens. 2019, 11, 727. [Google Scholar] [CrossRef] [Green Version]
  29. Jakimow, B.; van der Linden, S.; Thiel, F.; Frantz, D.; Hostert, P. Visualizing and labeling dense multi-sensor earth observation time series: The EO time series viewer. Environ. Model. Softw. 2020, 125, 104631. [Google Scholar] [CrossRef]
  30. Elfadaly, A.; Abouarab, M.A.R.; Shabrawy, R.R.M.E.; Mostafa, W.; Wilson, P.; Morhange, C.; Silverstein, J.; Lasaponara, R. Discovering potential settlement areas around archaeological tells using the integration between historic topographic maps, optical, and radar data in the northern Nile Delta, Egypt. Remote Sens. 2019, 11, 3039. [Google Scholar] [CrossRef] [Green Version]
  31. Gurjar, S.K.; Tare, V. Estimating long-term LULC changes in an agriculture-dominated basin using CORONA (1970) and LISS IV (2013–14) satellite images: A case study of Ramganga River, India. Environ. Monit. Assess. 2019, 191, 217. [Google Scholar] [CrossRef]
  32. Patel, S.; Shah, E.; Jayaprasad, P.; James, M.E. Changes in Antarctic coastline between 1997 and 2016 using RADARSAT and MODIS data. Int. J. Remote Sens. 2020, 41, 1389–1414. [Google Scholar] [CrossRef]
  33. Stromann, O.; Nascetti, A.; Yousif, O.; Ban, Y. Dimensionality reduction and feature selection for object-based land cover classification based on sentinel-1 and sentinel-2 time series using Google earth engine. Remote Sens. 2019, 12, 76. [Google Scholar] [CrossRef] [Green Version]
  34. Venter, Z.S.; Scott, S.L.; Desmet, P.G.; Hoffman, M.T. Application of Landsat-derived vegetation trends over South Africa: Potential for monitoring land degradation and restoration. Ecol. Indic. 2020, 113, 106206. [Google Scholar] [CrossRef]
  35. Zhang, Y.; Zhao, H. Land–use and land-cover change detection using dynamic time warping–based time series clustering method. Can. J. Remote Sens. 2020, 46, 67–83. [Google Scholar] [CrossRef]
  36. Sajjad, A.; Lu, J.; Chen, X.; Chisenga, C.; Saleem, N.; Hassan, H. Operational monitoring and damage assessment of riverine flood-2014 in the lower chenab plain, punjab, pakistan, using remote sensing and GIS techniques. Remote Sens. 2020, 12, 714. [Google Scholar] [CrossRef] [Green Version]
  37. Fortin, J.A.; Cardille, J.A.; Perez, E. Multi-sensor detection of forest-cover change across 45 years in Mato Grosso, Brazil. Remote Sens. Environ. 2020, 238, 111266. [Google Scholar] [CrossRef]
  38. Kalther, J.; Itaya, A. Coastline changes and their effects on land use and cover in Subang, Indonesia. J. Coast. Conserv. 2020, 24, 16. [Google Scholar] [CrossRef]
  39. Anderson, K.; Fawcett, D.; Cugulliere, A.; Benford, S.; Jones, D.; Leng, R. Vegetation expansion in the subnival Hindu Kush Himalaya. Glob. Chang. Biol. 2020, 26, 1608–1625. [Google Scholar] [CrossRef] [Green Version]
  40. Ma, W.; Domke, G.M.; Woodall, C.W.; D’Amato, A.W. Contemporary forest carbon dynamics in the northern U.S. associated with land cover changes. Ecol. Indic. 2020, 110, 105901. [Google Scholar] [CrossRef]
  41. Tao, S.; Fang, J.; Ma, S.; Cai, Q.; Xiong, X.; Tian, D.; Zhao, X.; Fang, L.; Zhang, H.; Zhu, J.; et al. Changes in China’s lakes: Climate and human impacts. Natl. Sci. Rev. 2019, 7, 132–140. [Google Scholar] [CrossRef] [Green Version]
  42. Rigge, M.; Homer, C.; Cleeves, L.; Meyer, D.K.; Bunde, B.; Shi, H.; Xian, G.; Schell, S.; Bobo, M. Quantifying Western U.S. rangelands as fractional components with multi-resolution remote sensing and in situ data. Remote Sens. 2020, 12, 412. [Google Scholar] [CrossRef] [Green Version]
  43. Uda, S.K.; Hein, L.; Sumarga, E. Towards sustainable management of Indonesian tropical peatlands. Wetl. Ecol. Manag. 2017, 25, 683–701. [Google Scholar] [CrossRef] [Green Version]
  44. MoEF. National Forest Reference Emission Level for Deforestation and Forest Degradation: In the Context of Decision 1/CP.16 para 70 UNFCCC (Encourages Developing Country Parties to Contribute to Mitigation Actions in the Forest Sector); DG-PPI MoEF: New Delhi, India, 2015. [Google Scholar]
  45. Papua Province. RREG Papua Province 2017–2020; Peatland Restoration Agency: Jakarta, Indonesia, 2017; pp. 1–2. [Google Scholar]
  46. Sasmito, S.D.; Sillanpää, M.; Hayes, M.A.; Bachri, S.; Saragi-Sasmito, M.F.; Sidik, F.; Hanggara, B.B.; Mofu, W.Y.; Rumbiak, V.I.; Hendri; et al. Mangrove blue carbon stocks and dynamics are controlled by hydrogeomorphic settings and land-use change. Glob. Chang. Biol. 2020, 26, 3028–3039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Jonatan, N.J.; Ekayuliana, A.; Dhiputra, I.M.K.; Nugroho, Y.S. The Utilization of Metroxylon Sago (Rottb.) dregs for low bioethanol as fuel households needs in Papua Province Indonesia. KLS 2017, 3, 150. [Google Scholar] [CrossRef] [Green Version]
  48. Thangavelu, S.K.; Rajkumar, T.; Pandi, D.K.; Ahmed, A.S.; Ani, F.N. Microwave assisted acid hydrolysis for bioethanol fuel production from sago pith waste. Waste Manag. 2019, 86, 80–86. [Google Scholar] [CrossRef] [PubMed]
  49. BPS. Papua Province in Figures 2020; BPS-Statistics of Papua Province, BPS: Jakarta, Indonesia, 2020. [Google Scholar]
  50. BPS. Merauke Regency in Figures 2020; BPS-Statistics of Merauke Regency, BPS: Jakarta, Indonesia, 2020. [Google Scholar]
  51. Leinenkugel, P.; Deck, R.; Huth, J.; Ottinger, M.; Mack, B. The potential of open geodata for automated large-scale land use and land cover classification. Remote Sens. 2019, 11, 2249. [Google Scholar] [CrossRef] [Green Version]
  52. Alonso-Sarria, F.; Valdivieso-Ros, C.; Gomariz-Castillo, F. Isolation forests to evaluate class separability and the representativeness of training and validation areas in land cover classification. Remote Sens. 2019, 11, 3000. [Google Scholar] [CrossRef] [Green Version]
  53. Gilani, H.; Shrestha, H.L.; Murthy, M.S.R.; Phuntso, P.; Pradhan, S.; Bajracharya, B.; Shrestha, B. Decadal land cover change dynamics in Bhutan. J. Environ. Manag. 2015, 148, 91–100. [Google Scholar] [CrossRef]
  54. Jeong, S.; Ko, J.; Yeom, J.-M. Nationwide projection of rice yield using a crop model integrated with geostationary satellite imagery: A case study in South Korea. Remote Sens. 2018, 10, 1665. [Google Scholar] [CrossRef] [Green Version]
  55. Ye, S.; Rogan, J.; Sangermano, F. Monitoring rubber plantation expansion using Landsat data time series and a Shapelet-based approach. ISPRS J. Photogramm. Remote Sens. 2018, 136, 134–143. [Google Scholar] [CrossRef]
  56. Gondwe, M.F.; Cho, M.A.; Chirwa, P.W.; Geldenhuys, C.J. Land use land cover change and the comparative impact of co-management and government-management on the forest cover in Malawi (1999–2018). J. Land Use Sci. 2019, 14, 281–305. [Google Scholar] [CrossRef]
  57. Entwistle, C.; Mora, M.A.; Knight, R. Estimating coastal wetland gain and losses in Galveston County and Cameron County, Texas, USA: Wetland Gains and Losses in 2 Texas Coastal Counties. Integr. Environ. Assess. Manag. 2018, 14, 120–129. [Google Scholar] [CrossRef] [Green Version]
  58. Martínez, M.L.; Pérez-Maqueo, O.; Vázquez, G.; Castillo-Campos, G.; García-Franco, J.; Mehltreter, K.; Equihua, M.; Landgrave, R. Effects of land use change on biodiversity and ecosystem services in tropical montane cloud forests of Mexico. Forest Ecol. Manag. 2009, 258, 1856–1863. [Google Scholar] [CrossRef]
  59. Mathan, M.; Krishnaveni, M. Monitoring spatio-temporal dynamics of urban and peri-urban land transitions using ensemble of remote sensing spectral indices—A case study of Chennai Metropolitan Area, India. Environ. Monit. Assess. 2020, 192, 15. [Google Scholar] [CrossRef]
  60. Holusa, J.; Lubojacky, J.; Knizek, M. Distribution of the double-spined spruce bark beetle Ipsduplicatus in the Czech Republic: Spreading in 1997–2009. Phytoparasitica 2010, 38, 435–443. [Google Scholar] [CrossRef] [Green Version]
  61. Purwestri, R.C.; Hájek, M.; Šodková, M.; Sane, M.; Kašpar, J. Bioeconomy in the National Forest Strategy: A comparison study in Germany and the Czech Republic. Forests 2020, 11, 608. [Google Scholar] [CrossRef]
  62. Food and Agriculture Organization of the United Nations (FAO) Global Forest Resources Assessment 2015: Desk Reference. 2017. Available online: http://www.fao.org/forest-resources-assessment/past-assessments/fra-2015/en/ (accessed on 9 August 2019).
  63. Martin, K.; Norris, A.; Drever, M. Effects of bark beetle outbreaks on avian biodiversity in the British Columbia interior: Implications for critical habitat management. JEM 2006, 7, 10–24. [Google Scholar]
  64. Bentz, B.; Logan, J.; MacMahon, J.; Allen, C.D.; Ayres, M.; Berg, E.; Carroll, A.; Hansen, M.; Hicke, J.; Joyce, L.; et al. Bark beetle outbreaks in western North America: Causes and consequences. In Bark Beetle Symposium; Snowbird, Utah; November, 2005; University of Utah Press: Salt Lake City, UT, USA, 2009. [Google Scholar]
  65. Zhang, X.; Treitz, P.M.; Chen, D.; Quan, C.; Shi, L.; Li, X. Mapping mangrove forests using multi-tidal remotely-sensed data and a decision-tree-based procedure. Int. J. Appl. Earth Obs. Geoinf. 2017, 62, 201–214. [Google Scholar] [CrossRef]
  66. Chew, T.A.; Isa AH, B.M.; Mohayidin, M.G.B. Sago (MetroxylonsaguRottboll, the forgotten palm. J. Sustain. Agric. 1999, 14, 5–17. [Google Scholar] [CrossRef]
  67. Kamlun, K.U.; Bürger Arndt, R.; Phua, M.-H. Monitoring deforestation in Malaysia between 1985 and 2013: Insight from South-Western Sabah and its protected peat swamp area. Land Use Policy 2016, 57, 418–430. [Google Scholar] [CrossRef]
  68. Salosa, S.T. Challenge of Sago (MetroxylonSp) as Papua Food Identity. KSS 2016, 1, 62–65. [Google Scholar] [CrossRef] [Green Version]
  69. Vetrita, Y.; Cochrane, M.A. Fire frequency and related land-use and land-cover changes in indonesia’s peatlands. Remote Sens. 2019, 12, 5. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geographical position of study area. The study area has boundaries with Mappi and BovenDigoel Regency to the north, the Arafuru Ocean to the south and west, Papua New Guinea to the east.
Figure 1. Geographical position of study area. The study area has boundaries with Mappi and BovenDigoel Regency to the north, the Arafuru Ocean to the south and west, Papua New Guinea to the east.
Sustainability 12 06623 g001
Figure 2. The study area: Landsat imagery path and row scenes.
Figure 2. The study area: Landsat imagery path and row scenes.
Sustainability 12 06623 g002
Figure 3. Merauke Regency land cover maps in 1990, 2019, and land use changes map.
Figure 3. Merauke Regency land cover maps in 1990, 2019, and land use changes map.
Sustainability 12 06623 g003
Figure 4. Percentage of non-forest and natural forest land cover changes during 29 years.
Figure 4. Percentage of non-forest and natural forest land cover changes during 29 years.
Sustainability 12 06623 g004
Table 1. Land cover classes of Indonesia 1.
Table 1. Land cover classes of Indonesia 1.
No.ClassAbbreviationDefinitionCategoryIPCC 1
1.Primary dryland forestP.D.F.The natural tropical forest grows on dryland habitat including lowland, upland, and mountain forests, with no signs of human or logging activities.Natural forestForest
2.Secondary dryland forestS.D.F.The natural tropical forest grows on dryland habitat including lowland, upland, and spotting of logging.Natural forestForest
3.Primary mangrove forestP.M.F.Inundated forest with access to sea/brackish water and dominated by species of mangrove and nipa that has no signs of logging activities.Natural forest Forest
4.Secondary mangrove forestS.M.F.Inundated forest with access to sea/brackish water and dominated by species of mangrove and nipa that exhibit signs of logging activities, indicated by patterns and spotting of logging.Natural forestForest
5.Primary swamp forestP.S.F.The natural tropical forest grows in wet habitat including brackish swamp, sago, and peat swamp.Natural forestForest
6.Secondary swamp forestS.S.F.The natural tropical forest grows in wet habitat including brackish swamp, sago, and peat swamp, with signs of human intervention or logging activities.Natural forestForest
7.Plantation forestP.F.Planted forest including areas of reforestation, industrial plantation forest, and community plantation forest.Plantation forestForest
8.Estate croplandE.C.Estate areas that have been planted, mostly with perennial crops or other agriculture trees commodities.Non-forestCropland
9.Pure dry agricultureP.D.A.All land cover associated with agriculture activities on dry/nonwet land such as moors, mixed garden, and agriculture fields.Non-forestCropland
10Mixed dry agricultureM.D.A.All land cover associated with agriculture activities on dry land that is mixed with shrubs, thickets, and logs over the forest. This cover type often results in shifting cultivation and its rotation, including on carts.Non-forestCropland
11Dry shrubD.S.Highly degraded log over areas on dryland habitat that is undergoing succession but have not yet hit a stable forest ecosystem, having natural scattered trees or shrubs.Non-forestGrassland
12Paddy fieldP.F.Agriculture areas on nondry habitats, especially for paddies that typically exhibit dyke patterns, rained, seasonal, and irrigated paddy fields.Non-forestCropland
13Wet shrubW.S.Highly degraded log over areas on nondryland habitat of wet habitat that have not yet reached a stable forest ecosystem, having natural separated trees or shrubs.Non-forestGrassland
14Savanna and grassesS.G.Areas with grasses and scattered trees and shrubs, could be in wet or nonwet habitats. Typical of the natural ecosystem on Sulawesi Tenggara, NTT, and the south part of Papua.Non-forestGrassland
15Open swampO.SOpen swamp with less vegetation.Non-forestWetland
16Open waterO.W.Water body including ocean, rivers, lakes, ponds.Non-forestWetland
17FishpondF.P.Areas exhibiting aquaculture activities such as fish, shrimp, or salt ponds.Non-forestWetland
18Port and harborP.H.Ports and harbors big enough to be independently delineated as independent objects.Non-forestOther land
19Transmigration areasT.AUnique settlement, association with houses, agroforestry, and garden in and around.Non-forestSettlement
20Settlement areasS.A.Including rural, urban, industrial, and other typical appearances of settlement.Non-forestSettlement
21Mining areasM.AMining areas such as open-pit mining, tailing ground.Non-forestOther land
22Bare groundB.GBarren land with no vegetation cover yet, open exposure areas, craters, sandbanks, sediments, and area post-fire that has not shown regrowth yet.Non-forestOther land
23Clouds and no-dataC.OoClouds, cloud shadows with size up to 4 cm 2 at 100.000 scales.Non-forestNo data
1 [44] (pp. 5–6 and 49–51); 2 Intergovernmental Panel on Climate Change (IPCC).
Table 2. Characteristics of the satellite data used in this study.
Table 2. Characteristics of the satellite data used in this study.
PropertyLandsat 5Landsat 7 Landsat 8
Spatial resolution30 m for visible and I.R.,30 m for visible and Infrared (I.R.)30 m for visible and I.R.
120 m for thermal15 m for Panchromatic (Pan) and 60 m for thermal15 m for (Pan) and 100 m for thermal
Spectral resolution7 bands (visible, I.R., thermal band)8 bands (visible, I.R., Pan, and thermal band)11 bands (visible, I.R., Pan, and thermal)
Radiometric resolution8 bit8 bit16 bit
Temporal resolution16 day16 day16 day
Details of spectral resolutions (μm)Band 1: (blue) 0.450–0.515Band 1: (blue) 0.450–0.515Band 1: (blue) 0.43–0.45
Band 2: (green) 0.525-0.605Band 2: (green) 0.525–0.605Band 2: (blue-green) 0.45–0.51
Band 3: (red) 0.63-0.69Band 3: (red) 0.63–0.69Band 3: (green) 0.53–0.59
Band 4: Near-Infrared (N.I.R.) 0.76–0.90Band 4: (N.I.R.) 0.76–0.90Band 4: (red) 0.64–0.67
Band 5: Short-Wave Infrared (SWIR-1) 1.55–1.75Band 5: (SWIR-1) 1.55–1.75Band 5: (N.I.R.) 0.85–0.88
Band 6: (thermal) 10.4–12.5Band 6: (thermal) 10.4–12.5Band 6: (SWIR-1) 1.57–1.65
Band 7: (SWIR-2) 2.09–2.35Band 7: (SWIR-2) 2.09–2.35Band 7: (SWIR-2) 2.11–2.29
Band 8: (Pan) 0.52–0.92Band 8: (Pan) 0.50–0.68
Band 9: (Cirrus) 1.36–1.38
Band 10: (Thermal I.R.) 10.60–11.19
Band 11: (Thermal I.R.) 11.50–12.51
Table 3. Land cover area (ha) and the percentage of change from 1990 to 2003.
Table 3. Land cover area (ha) and the percentage of change from 1990 to 2003.
LC Class1990 (ha)1996 (ha)2000 (ha)2003 (ha)
Natural Forest
Primary dryland forest694,737664,757634,776619,004
Secondary dryland forest638,049620,773603,496618,381
Primary mangrove forest208,727207,345205,963201,768
Secondary mangrove forest25,34524,20923,07325,776
Primary swamp forest342,429329,304316,179292,789
Secondary swamp forest531,109419,213307,317313,173
Total area (ha)2,440,3962,265,6002,090,8042,070,891
Percentage of change (%)50.3046.7043.0942.68
Change rate (ha/yr)=−7.1626−7.715−0.952
Non-Forest
Bush/shrub71,94624,194176,443177,229
Estate crop plantation---101
Settlement area3160336635713667
Barren land81,71451,75921,80521,805
Cloud covered764764764764
Savanna/grassland471,693549,087626,480646,258
Water body352,031352,012351.993351,992
Swamp shrub930,069931,438932,806929,360
Dryland agriculture14,37715,36816,35816,772
Shrub-mixed dryland farm43,46249,01354,56354,563
Paddy field10,93210,93210,93210,974
Fishpond----
Airport/harbor159159159159
Transmigration area36,63841,43046.22146,221
Swamp394,375456,596518,816521,051
Total area (ha)2,411,3192,586,1152,760,9122,780,824
Percentage of change (%)49.7053.3056.9157.32
Change rate (ha/yr)=724967590.721
Table 4. Land cover area (ha) and the percentage of change in two categories from 2006 to 2014.
Table 4. Land cover area (ha) and the percentage of change in two categories from 2006 to 2014.
LC Class2006 (ha)2009 (ha)2011 (ha)2014 (ha)
Natural Forest
Primary dryland forest598,828553,728553,098543,670
Secondary dryland forest627,494672,086672,425678,803
Primary mangrove forest196,510196,510196,510197,808
Secondary mangrove forest23,67823,57423,57423,675
Primary swamp forest238,249205,343205,343206,530
Secondary swamp forest338,909371,810371,810374,446
Total area (ha)2,023,6682,023,0512,022,7602,024,932
Percentage of change (%)41.7141.7041.6941.74
Change rate (ha/yr)−2280−0.030−0.014−0.107
Non-Forest
Bush/shrub178,032178,463177,262174,273
Estate crop plantation101101153316,535
Settlement area3891389138913917
Barren land21,85321,85321,91323,501
Cloud covered764764764-
Savanna/grassland655,175704,034704,044708,703
Water body351,995351,994351,994322,264
Swamp shrub949,786900,908900,838906,111
Dryland agriculture16,80316,88016,88017,184
Shrub-mixed dryland farm65,25065,37965,37965,760
Paddy field10,97410,97411,04411,463
Fishpond----
Airport/harbor159159159159
Transmigration area46,22146,22146,22146,440
Swamp527,044527,044527,034530,472
Total area (ha)2,828,0472,828,6642,828,9552,826,783
Percentage of change (%)58.2958.3058.3158.26
Change rate (ha/yr)16980.0220.010−0.077
Table 5. Land cover area (ha) and the percentage of change in two categories from 2015 to 2019.
Table 5. Land cover area (ha) and the percentage of change in two categories from 2015 to 2019.
LC Class2015 (ha)2016 (ha)2017 (ha)2018 (ha)2019 (ha)
Natural Forest
Primary dryland forest529,715522,977519,144401,879500,359
Secondary dryland forest664,888654,663652,518732,934631,295
Primary mangrove forest196,758195,162195,007195,660195,384
Secondary mangrove forest23,52123,87623,82923,93224,060
Primary swamp forest202,799200,958200,400202,694202,193
Secondary swamp forest359,399356,270358,089357,151531,266
Total changed area (ha)1,977,0801,953,9061,948,9871,914,2502,084,557
Percentage of change (%)40.7540.2740.1739.4642.97
Change rate (ha/yr)−2363−1172−0.252−17828897
Non-Forest
Bush/shrub169,262166,111170,801169,65629,465
Estate crop plantation19,88527,39753,85780,2314359
Settlement area36533878348072167090
Barren land263,85975,08156,53977,99488,946
Cloud covered-----
Savanna/grassland568,723700,156603,422576,528555,274
Water body322,282351,749351,734349,816349,884
Swamp shrub860,813917,482969,770978,818942,998
Dryland agriculture16,39617,07216,37718,27821,671
Shrub-mixed dryland farm62,13965,07165,34470,6928600
Paddy field11,45911,38811,38848,79545,505
Fishpond---44880
Airport/harbor159159159175175
Transmigration area46,44046,15245,50426,52625,575
Swamp529,565516,113554,354532,29137,538
Total area (ha)2,874,6352,897,8092,902,7282,937,4652,767,158
Percentage of change (%)59.2559.7359.8360.5457.03
Change rate (ha/yr)16930.8060.1701197-5798
Table 6. Land cover changes of each class in Merauke Regency, Papua, Indonesia, from 1990 to 2019.
Table 6. Land cover changes of each class in Merauke Regency, Papua, Indonesia, from 1990 to 2019.
L.C. ClassChanged Rate (ha/yr)Total Changed Area
Gain(+)Loss(−)Net(±)Ha%
Primary dryland forest8206.6724,404.83−16,198.17−194,378.00−27.98
Secondary dryland forest12,976.9213,539.75−562,83−6754.00−1.06
Primary mangrove forest162.581274.50−1111.92−13,343.00−6.39
Secondary mangrove forest282.60389.74−107,14−1285.70−5.07
Primary swamp forest290.0811,976.42−11,686.33−140,236.00−40.95
Secondary swamp forest20,255.2520,242.1713.08157.000.03
Bush/shrub9267.264474.004793.2657,519.1079.95
Estate crop plantation7863.22-7863.2294,358.60-
Settlement area393.1965.68327.523930.23124.38
Barren land22,871.7522,269.04602.717232.508.85
Cloud covered-63.64−63.64−763.65−100
Grassland30,703.5823,738.506965.0883,581.0017.72
Water body2462.992641.91−178.93−2147.15−0.61
Swamp shrub12,203.4211,126.001077.4212,929.001.39
Dryland agriculture731.42123.60607.827293.8050.73
Shrub-mixed dryland2570.88476.092094.7825,137.4057.84
Paddy field3161.55280.482881.0834,572.90316.26
Fishpond37.3530.716.6479.67-
Airport1.380.021.3616.3010.27
Transmigration area816.871738.78−921.92−11,063.00−30.20
Swamp14,529.0010,932.083596.9243,163.0010.94
Table 7. General characteristics of the prediction of sago palm habitat in Merauke Regency (N = 8).
Table 7. General characteristics of the prediction of sago palm habitat in Merauke Regency (N = 8).
District19902019
Animha16,975.06 ± 12,669.66 (125.05, 36,055.10)16,983.27 ± 13,440.54 (437.21, 43,225.50)
Elikobel18,216.17 ± 33,037.63 (0.00, 96,199.60)17,495.36 ± 34,366.44 (0.00, 96,199.60)
Ilwayab21,715.66 ± 27,109.94 (0,00, 80,960.40)21,900.72 ± 21,947.41 (0.00, 80,960.40)
Jagebob15,381.23 ± 21,832.78 (841.02, 6,459.20)15,082.23 ± 23,852.07 (841.02, 64,359.20)
Kurik8312.54 ± 7576.85 (482.92, 20,296.30)7971.55 ± 6519.25 (482.92, 20,296.30)
Kaptel27,868.15 ± 22,683.46 (2305.88, 63,204.40)27,699.84 ± 20,332.62 (5649.91, 63,317.30)
Kimaam45,414.98 ± 63,137.85 (434.43, 181,539.00)45,413,29 ± 64,595.09 (20.49, 1,9042.00)
Malind4975.18 ± 3897.71 (0.00, 10,343.10)4883.73 ± 4361.10 (0.00, 11,261.90)
Merauke15,106.11 ± 22,450.10 (0.00, 60,288.60)15,126.20 ± 18,966.08 (0.00, 60,228.60)
Muting39,547.46 ± 46,345.29 (3699.06, 118,800.00)39,489.52 ± 41,776.44 (3954.98, 112,000.00)
Naukenjerai10,682.42 ± 16,181.12 (0.00, 46,611.40)9737.45 ± 10,263.51 (0.00, 50,872.70)
Ngguti40,029 ± 24,706.78 (11,205.70, 70,419.90)38,816.99 ± 28,004.91 (11,293.40, 89,504.20)
Okaba17,481.54 ± 32,317.97 (37.903, 94,925.10)18,900.30 ± 25,688,46 (505.92, 76,309.70)
Semangga2702.08 ± 3717.96 (0.00, 10,013.90)2702.08 ± 3978.13 (0.00, 11,292.90)
Sota31,005.28 ± 35,952.15 (1608.15, 110,369.00)30,931.49 ± 29,605.80 (6158.29, 99,919.30)
Tanah Miring16,414.01 ± 9013.24 (4923.28, 28,562.60)16,359.38 ± 11,134.00 (347.94, 30,763.80)
Tabonji33,038.64 ± 40,211.73 (0.00, 111,643.00)33,030.82 ± 42,503.87 (0.00, 120,527.00)
Tubang26,192.51 ± 108,655.74 (1027.90, 73,437.70)32,655.19 ± 42,503.87 (8534.63, 89,472.00)
Ulilin56,469.20 ± 108,655.74 (1409.11, 315,111.00)56,017.52 ± 103,832.79 (742.51, 299,073.00)
Waan60,482.52 ± 60,071.43 (0.00, 165,196.00)61,382.05 ± 64,152.86 (742.51, 29,9073.00)
Table 8. Land cover changes from the natural habitat of sago in 1990 and 2019.
Table 8. Land cover changes from the natural habitat of sago in 1990 and 2019.
LC19902019p-Value
Primary dryland34,736.82 ± 71,532.46 (0.00, 315,111.00)27,686.42 ± 67,227.85 (0, 299,073.00)0.015
Secondary dryland31,902.33 ± 38,007.26 (1.02, 118,800.00)33,604.22 ± 39,934.11 (0, 112,000.00)0.313
Primary swamp forest17,126.28 ± 23,169.16 (1276.23, 107,615)10,271.99 ± 8519.85 (531.72, 24,711.10)0.107
Secondary swamp forest26,555.19 ± 24,072.41 (4668.14, 94,925.10)18,590.47 ± 23,439.27 (949.07, 105.92)0.152
Bush/shrub3597.31 ± 6055.62 (0, 24,048.80)8923.07 ± 16,655.05 (0, 63,317.30)0.081
Grassland23,585.31 ± 36,748.43 (0, 111.643.00)35,202.67 ± 42,540.96 (0, 152,745.00)0.002
Swamp shrub46,503 ± 52,913.31 (51.08, 181,539.00)45,045.15 ± 50,975.60 (51.08, 190.427)0.723
Swamp19,197.62 ± 16,473.24 (79.92, 62,207.50)25,707.58 ± 17,481.00 (34.41, 68,235.40)0.007

Share and Cite

MDPI and ACS Style

Letsoin, S.M.A.; Herak, D.; Rahmawan, F.; Purwestri, R.C. Land Cover Changes from 1990 to 2019 in Papua, Indonesia: Results of the Remote Sensing Imagery. Sustainability 2020, 12, 6623. https://0-doi-org.brum.beds.ac.uk/10.3390/su12166623

AMA Style

Letsoin SMA, Herak D, Rahmawan F, Purwestri RC. Land Cover Changes from 1990 to 2019 in Papua, Indonesia: Results of the Remote Sensing Imagery. Sustainability. 2020; 12(16):6623. https://0-doi-org.brum.beds.ac.uk/10.3390/su12166623

Chicago/Turabian Style

Letsoin, Sri Murniani Angelina, David Herak, Fajar Rahmawan, and Ratna Chrismiari Purwestri. 2020. "Land Cover Changes from 1990 to 2019 in Papua, Indonesia: Results of the Remote Sensing Imagery" Sustainability 12, no. 16: 6623. https://0-doi-org.brum.beds.ac.uk/10.3390/su12166623

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