Earth Engine Based Three Decadal Landsat Imagery Analysis for Mapping of Mangrove Forests and Its Surroundings in the Trat Province of Thailand

Monitoring and understanding the changes in mangrove ecosystems and their surroundings are required to determine how mangrove ecosystems are constantly changing while influenced by anthropogenic, and natural drivers. Consistency in high spatial resolution (30 m) satellite and high performance computing facilities are limiting factors to the process, with storage and analysis requirements. With this, we present the Google Earth Engine (GEE) based approach for long term mapping of mangrove forests and their surroundings. In this study, we used a GEE based approach: 1) to create atmospheric contamination free data from 1987-2017 from different Landsat satellite imagery; and 2) evaluating the random forest classifier and post classification change detection method. The obtained overall accuracy for the years 1987 and 2017 was determined to be 0.87 and 0.96, followed by a Kappa coefficient 0.80 and 0.94. The change detection results revealed a significant decrease in the agricultural area, while there was an increase in mangrove forest, shrimp/fish farm, and bareland area. The results suggest that interconversion of land use and land cover is affecting the landscape dynamics within the study area. How to cite this paper: Pimple, U., Simonetti, D., Sitthi, A., Pungkul, S., Leadprathom, K., Skupek, H., Som-ard, J., Gond, V. and Towprayoon, S. (2018) Google Earth Engine Based Three Decadal Landsat Imagery Analysis for Mapping of Mangrove Forests and Its Surroundings in the Trat Province of Thailand. Journal of Computer and Communications, 6, 247-264. https://doi.org/10.4236/jcc.2018.61025 Received: December 8, 2017 Accepted: December 26, 2017 Published: December 29, 2017


Introduction
Mangrove forests are located throughout the tropical and sub-tropical regions of the world, and are claimed to be one of the most vulnerable ecosystems to be affected by natural disturbances and human interference [1] [2] [3] [4].Mangrove forests are unique ecosystems that provide important ecological services for coastal habitat and coastal protection [5] [6].These ecosystems, however, are under high pressure due to over exploitation, and are declining at an alarming rate [2] [6].Despite this rapid decline, and their importance, mangrove forest ecosystems have not received much publicity, particularly in regard to conservation and rehabilitation [6] [7].Quantifying and monitoring the spatial and temporal dynamics of the mangrove ecosystem is essential for a better understanding of the many coastal land and sea processes.Traditionally, mapping a mangrove forest requires intensive field work, which is costly in time and money, as mangroves are inaccessible or difficult to field survey [6] [8].Satellite remote sensing has a great potential for mapping and monitoring changes in mangrove forests, as the space based technology allows for collecting information from the landscape which is otherwise particularly difficult to access [9] [10].
[2] [11] have provided a detailed summary with an overview of remote-sensing research activities, including critical analysis that has been performed in the last few decades.Green et al. (1998) and Kuenzer et al. (2011) highlighted the importance of understanding the local environment when using remote sensing based mapping and monitoring.In recent years, several studies have been published, illustrating hyperspectral airborne and spaceborne data applications, including AISA, CASI, Hymap, AVIRIS, Dedalus, and EO-1 hyperian [6] [11] [12] [10].High resolution imagery applications such as IKONOS, Quickbird, Rapid Eye and WorldView-3, have been very effective in discriminating mangrove forests from other forms of land use [13] [14] [15] [16].However, limited spectral bands, complex data collection, and analysis methods, along with their high cost, are major limiting factors in using hyperspectral and high resolution data [2] [9] [17].Few studies reported the use of RADAR data, with the inclusion of ALOS PALSAR, ERS-1/2 and Radarsat-1 SAR, as a tool in the mangrove classification framework, having found that classification and mapping accuracy requires improvement [2] [18] [19] [20].In previous studies, optical remote sensing imagery, like the Landsat Multispectral Scanner (MSS), Landsat Thematic Mapper (TM), and Indian Remote Sensing satellites (IRS), SPOT XS, and SPOT-5, have been commonly used for mangrove forest mapping as the data is available via free access, or at a low cost [10] [21].The recent developments in the series of Landsat satellites, such as the Operational Land Imager (OLI)-8 and Sentinel-2, have added a new dimension to long term data.The availability of multiple spectral band records over the long term, means the Landsat series can be used for accurate mapping and monitoring of mangrove forests [22].Additionally, the application of the non-parametric or machine learning classifiers are very efficient for land use and land cover (LULC) mapping, even if still not abundant in mangrove mapping studies [23].
In recent years, there has been an increase in high-performance cloud computing platforms, such as the NASA Earth Exchange (NEX), Amazon Web Service (AWS) and Google Earth Engine (GEE).These high performance cloud computing platforms allow free access to the vast and fast growing earth observation data for global, as well as regional studies [24] [25].For example, GEE provides preprocessed Landsat data (1982-present), along with the required disk space and advanced classification machine learning algorithms [25].
In Southeast Asia, large areas of the coastal zones have been occupied by mangrove forests [26].According to the Asian Development Bank Regional Review on the Economics of Climate Change in Southeast Asia [27], the reduction in the size of mangroves resulted in coastal erosion in Thailand, and in neighboring countries.While the extent of these changes remains limited, means of sustainable management and future rehabilitation remains highly uncertain [2] [11] [28].In recent decades, Thailand's mangrove forest area has substantially decreased as a result of human settlements, transport infrastructure, agriculture, and aquaculture [26] [29].According to NESDB, while inconclusive in most regions, 30% of mangrove forest was lost during 1961-1996 due to the conversion of mangrove forest to shrimp farms.
Given the above factors, the main objective of this study is to quantify the presence of mangrove forests in Thailand's Trat Province over the last 30 years (1987-2017) using Landsat imagery and GEE cloud computing, as well as developing an operational wall-to-wall change detection methodology based on long time series analysis.

Study Area
The study area is located in the Trat Province, eastern Thailand, on the border with Cambodia and along the Gulf of Thailand (Figure 1).It covers an area of approximately 240 km 2 of witch 106 km 2 of mangrove forest.Since ever, local communities have benefit from goods and services provided by the forest but in recent years, mangroves have been heavily exploited for timber extraction, charcoal production and shrimp farms.[30] reported that, after the 1980's, thanks to the effort of the Department of Marine and Coastal Resources, Thailand has reaerated this area.Villagers to this day, continue to use mangrove wood for various domestic purposes.

Landsat Imagery
Data consists of Landsat TM, ETM+, and OLI Tier 1 top-of-atmosphere (TOA) reflectance as obtained from GEE image collections.The red, green, blue, NIR, SWIR-1, and SWIR-2 spectral bands of the TM-5, 7 ETM+, and OLI-8 platforms, were considered in the analysis.Jagged pixels at the edges of Landsat TM-5 images were removed using the 450 m inward buffer, which ensures the best available reflectance values for image analysis [31].Annual composites were created by using the median reflectance values of the collection (all images for a target year e.g.1987, from 01/01/1987 to 31/12/1987) [24], after been cleaned from cloudy or no-data pixels following the algorithm proposed by [32], and available in GEE.The algorithm is driven by predefined knowledge-based rules built upon the spectral signature collected on a global scale, and generates a thematic output including a cloud mask.
The TOA imagery was atmospherically corrected using a Dark Object Subtraction (DOS) method [33].Using a forest normalization method, the median value of the mangrove forest pixels was used to apply a linear shift to each spectral band [34] [35].Prior to classification, the Normalized difference vegetation index (NDVI) [11] [36], and Normalized difference infrared index (NDII) [37], has be computed to mask the composite from residual clouds or no data.

Training and Validation Data
One of the major issues of classifying historical images by using training and va-lidation samples, is often the lack of field data or so called ground-truth [38].The design of a systematic training and validation dataset across a specified area is crucial for identifying major changes in mangrove forests over time.The design must be a good representation of major Land Use Land Cover (LULC) classes, with the dataset sufficiently large to provide reliable estimates [35].The unsupervised classification was performed on the Landsat OLI-8 year 2017 images to establish thematic categories of LULC.The most recent image was chosen for stratification to avoid the effect of the LULC change on training and testing design.Later the stratified random sampling approach was used to estimate the total number of samples per class [39] [40].The stratification and selection of independent sample design focused on each five-year interval to ensure a stable change identification [39].In total, 414 sample locations [Figure 2] were selected for LULC classification: class 1, active agriculture; class 2, bare land and urban areas (some agricultural land without vegetation was included here); class 3, mangrove forest; and class 4, shrimp and fish farms.To establish each class, training samples were obtained: 150 for agriculture, 50 for bare land and urban areas, 164 for mangrove forest, and 50 for shrimp and fish farms.Among these sample locations, 109 were randomly selected to be set aside as validation samples.In March, 2015 and October 2016, a field mission was conducted in the study area to collect training and validation data for mangrove mapping.About 60 mangrove forest samples were collected during these missions, and were used along with a combination of Google Earth images, high-resolution satellite imagery, aerial photographs, and prior knowledge, for use as samples of the remaining classes.A distance of 500 meters separated each sample to avoid spatial autocorrelation, while training and validation pixels remained independent of each other.In the process of classification, the reference training and testing data were first developed for each year and then used in the classifier [41].

Pixel Based Random Forest (RF) Classifier
We performed a supervised pixel-based classification using a Random Forest (RF), a tree based classifier that includes K-decision trees [42] [43]).RF overcome the problem of overfitting by constructing an ensemble of decision trees [43] [44].[44] reported that there are accurate and higher performance RF classifiers in land cover classification studies.The RF classifier was used to classify the extent of mangrove vegetation and other LULC in a study area, as shown in Figure 1.We trained the RF classifier (20 trees) in the GEE environment with 305 training samples to then classify the annual composite Landsat images into four LULCs: active agriculture (orchid plants, coconut grows, oil palm and rubber plantation), bare land/urban area/non-active agriculture plots, primary mangrove forest, and water bodies (fish and shrimp farms, and other water resources including water canals within mangrove forest).The categories were based on a detailed analysis of the study area, and reviewing previous studies and field surveys conducted for training and testing data.The RF classifier was performed on red, green, blue, NIR, SWIR-1, and SWIR-2 spectral bands of each annual composite.

Validation
During the classification process, composite from different years are trained and validated individually.About 70% of the sample points are used to train the classifier, while the remaining 30% of samples were used to test the accuracy and validate the RF classifier (Figure 2).The RF classifier accuracy and Kappa statistics is assessed by an error matrix.The final maps were compared with high resolution aerial imagery available in Google Earth for visual refinement.

Post Classification Change Detection
Several methods such as an image overlay, change vector analysis, image rationing, and principal component analysis have been used in LULC mapping studies [44] [45].In this study, the "From-to" change detection algorithm (Post-classification comparison), has been used to provide detailed information about the type of LULC change [46] [47].The main advantage of post classification change detection is indicating the nature and magnitude of the LULC changes that had taken place over a time.

Cloud Free Annual Mosaic of Landsat Series
Clouds and shadows represent one of the main sources of issue while working with optical remote sensed imagery such as Landsat, particularly when working in tropical regions.Figure 3(a) illustrates the influence of clouds, haze and missing pixels on Landsat series imagery, which could be the main limiting factor on the spatial and temporal consistency of long term mangrove ecosystem changes mapping and monitoring.The cloud cover, shadows, availability of haze and missing data, influence many data analysis processes including inaccurate atmospheric correction, biased vegetation indexes, mistakes in land cover classification and false detection of land use and land cover change [48].The cloud free seamless mosaic of Landsat series imagery was created with a predefined knowledge-based rule built upon the spectral signature [32].Annual composites were created by taking median reflectance values of the collection.Seamless and cloud free image mosaics can be important when mapping mangrove forests, because the cloud and seams can affect the visual interpretation of training sample collection, or leading to erroneous classifications [49].Cloud free and seamless mosaic images would likely improve the results of the forest normalization method described in Section 2.2.
The comparable visual results presented in Figure 3

Magnitude of Mangrove and Its Surrounding Change
In order to test the performance of consistent long term imagery from different Landsat sensors for mangrove forests and their surroundings, Landsat imagery was used representing two time periods (1987 and 2017).The stratified buffer was generated around the mangrove forest.A 850 meter buffer was generated to delineate the potential loss or gain of mangrove forest [Figure 1] within the surrounding landscape.The buffer was designed and based on the change in mangrove pixels within the study area using Normalized Differential Vegetation Index (NDVI), Normalized Difference Infrared Index (NDII), Digital Elevation Model(DEM), and Automatic Classification [10] [11] [32] [36] [50] [51].A visual refinement has been carried out for year 1987 and 2017 (Landsat TM-5, ETM+ 7, and OLI-8 imagery respectively), to ensure the high quality of the mask over the time interval.
Independent training and validation data was used for the years 1987 and 2017.The LULC maps for the years 1987 and 2017 were based on Landsat TM-5 and OLI-8 satellite imagery that were prepared with four LULC types: class 1, active agriculture, class 2, bare land and urban area (some agricultural land without vegetation was included here), class 3, mangrove forest and class 4, shrimp and fish farms.Figure 4 shows the final classification of the RF classifier, which consists of classified maps of the study area, for the year 1987 and 2017.
Table 1 shows the results obtained from the classified map of 1987 and 2017.
An overall accuracy of 0.87 and 0.96, followed by a Kappa coefficient 0.80 and 0.94, were obtained.As a result, the performance of the RF classifier for the year 2017 produced a higher accuracy classification map with an overall accuracy of 0.96.On the other hand, the classification performance in 1987 was less than  Table 2 summarizes the results of the LULC change in the study area of each LULC class.The agriculture and mangrove forest was the main LULC in 1987 with 49.18% and 34.20%, followed by bareland/urban areas and water bodies/shrimp farms with 9.80% and 6.82% respectively.Agricultural area decreased from 49.18% (111.36 km 2 ) in 1987 to 41.25% (93.40 km 2 ) in 2017, while, mangrove forest area increased from 34.20% (77.43 km 2 ) in 1987 to 36.17% (81.90 km 2 ) in 2017.Additionally, some of the disturbed mangrove forest area has shown significant recovery (Figures 5(a)-(d)).The bareland/urban areas increased from 9.80% (22.17 km 2 ) in 1987 to 11.08% (25.10 km 2 ) in 2017.Shrimp/fish farms progressively increased from 6.82% (15.16 km 2 ) in 1987 to 11.49% (26.02 km 2 ).
The decline in active agriculture was observed to be 7.93% between 1987 and 2017.Bareland/urban areas and shrimp/fish farms did experience an expansion during the period of 1987 to 2017.Another increased rate of change was observed in the mangrove forest.

LULC Change Transition from 1987 to 2017
The transition matrix was used to analyse the rates of LULC conversion from one LULC to another for the years 1987 and 2017.The corresponding probablities of change are shown in Table 3.There was a major conversion of agriculture to bareland (15.15 km 2 ) and from agricuture to mangrove forest (8.40 km 2 ) during this period.At the same time, mangrove forest changed to shrimp/fish farm (4.72 km 2 ), and bareland was converted to shrimp/fish farms.

Change in LULC
The results show that a significant change occurred in land cover, particularly in mangrove forests, in the surrounding areas of Trat, between 1987 and 2017.It appears that agriculture, bareland and shrimp farmshad undergone major changes.However, there is very little research has been conducted study of these changes.The mangrove forests made a significant recovery over time.This trend is indicative of the local community's awareness for mangrove forest conservation and in the detrimental effect that shrimp farming can have on mangrove forest conservation.[56] reported on the expansion of shrimp farming , and low-salinity shrimp farming (1996)(1997)(1998)(1999)(2000)(2001)(2002).He also reported that recently, (2003)(2004), the Thai government's policies restricted the expansion of low-salinity shrimp farming within the freshwater regions of the country.The result of this study suggests that the rapid changes in agriculture, bareland, and shrimp farms, and their interconversion, are a major driver of the change.

Conclusions
This paper presents a new strategy in attempting to achieve error free 30 year annual composites of Landsat satellites imagery for mapping mangroves and their surrounding LULC changes on the GEE cloud computing platform.This strategy uses pre-defined knowledge-based rules to remove contaminated pixels from all available imagery and uses annual median reflectance values in the collection.A fast, accurate and stable detection of change in agriculture, bare lands, mangrove forests, and shrimp/fish farms generated from consistent seamless mosaic and the RF classifier, demonstrates these results.The study area experienced drastic interchange between agriculture, bare land and shrimp/fish farms, while mangrove forests had made a recovery over a period of time.
The study contributes to the application of cloud computing GEE and its potential for costal ecosystem mapping and monitoring.The provided reliable and consistent long term satellite data and high performance classification approach could be beneficial for finding changes in mangrove ecosystems and their surroundings to fill the gaps necessary for forest management, conservation, as well as in understanding their carbon sequestration potential.

Figure 1 .
Figure 1.Location and extent of the study area in the Trat Province of Thailand including 850 mt buffer around mangrove forest (source: Google Earth 2017).

Figure 2 .
Figure 2. Locations of training and testing samples in the study area.A total of 414 points were created using the stratified random sampling approach, of which 109 were used as validation points (a) Training and testing samples for the year 1987; (b) Training and testing samples for the year 2017 (The figure contains back ground Google Earth engine map background and Landsat image for individual year).
(a) and Figure 3(b), were acquired by Landsat TM-5, ETM+ 7, and OLI-8.Their corresponding pseudo color composite (SWIR1, NIR, Red) before and after pre-processing for Landsat TM-5, ETM+ 7 OLI-8 are shown in Figure 3(a) and Figure 3(b).Several images in Figure 1(a) contains cloud, shadows, haze and missing data.The automatic rule-based algorithm was used to remove contamination in individual Landsat imagery.Figure 3(b) illustrates the cloud, haze and missing data free annual composite for the period of 1987-2017.In addition, the obtained composite revealed more vivid tone when compared with the original imagery.The selected contaminated pixels were tested using NDVI and NDII to mask the composite from residual clouds or no data.Only a selected annual composite image was classified, and stratified random samples generated.Journal of Computer and Communications

Table 1 .
Error matrix and accuracy statistics for classification for the years 1987 and 2017 (class 1, active agriculture; class 2, bare land and urban area (some agricultural land without vegetation was included here); class 3, mangrove forest; and class 4, shrimp and fish farms).
Overall accuracy = 0.96, Kappa statistics = 0.94 Journal of Computer and Communications

Table 2 .
Change in LULC area and percentage in study area.

Table 3 .
Land Use Land Cover (LULC) change matrix between 1987 and 2017.: The bold numbers indicate that there is no change in LULC over the study period. Note

Overcoming the Limitations with Landsat Series Imagery for Mangrove Forest Mapping
 easy and user friendly programming environment  virtually unlimited processing power with high computationally efficiency  raw satellite data are already available on GEE servers hence no need to retrieve and download huge amount of data (less local resources) presence of atmospheric contamination and sensor artifacts such as the SLC-Off (Landsat-7) and missing pixel data.Benefits introduced by implementing this methodological framework in GEE are threefold: