Change Detection of Mangrove Forests in Coastal Guangdong during the Past Three Decades Based on Remote Sensing Data

: Mangrove forests are among the most productive ecosystems on Earth and mainly grow at tropical and subtropical latitudes. They provide many important ecological and societal functions. However, rapid spatiotemporal variations in mangroves have been observed worldwide, especially in the coastal zones of developing areas, and the integrity of mangroves has been signiﬁcantly a ﬀ ected by anthropogenic activities in recent decades. The goal of this study was to determine the spatiotemporal characteristics of mangrove distribution over the past 30 years in Guangdong Province. This goal was achieved by classifying multi-temporal Landsat images using a decision tree method based on Classiﬁcation and Regression Tree (CART) algorithm. The driving forces resulting in these spatiotemporal variations of mangroves were then discussed. Our analysis revealed that the classiﬁcation method used in this study yielded good accuracy, with an overall accuracy and kappa coe ﬃ cient of higher than 90% and 0.8, respectively. In Guangdong province, the mangrove forests covered areas of 9305, 9556, 6793, and 9700 ha in 1985, 1995, 2005, and 2015, respectively, with remarkable inter-annual changes. Mangrove forests are mainly located in Western Guangdong, and few are located in Eastern Guangdong. The distribution of mangrove patches became more fragmented from 1985 to 2005 and less fragmented from 2005 to 2015, and the distribution pattern in 2015 showed stronger connectivity than that in 1985. Natural factors, such as temperature, sea level rise, extreme weather events, and the length of the coastline, have macroscopic e ﬀ ects on the distribution of mangrove forests. Anthropogenic activities, such as deforestation, urbanization, and aquaculture development, have negative e ﬀ ects on the distribution of mangroves. On the other hand, the establishment of nature reserves has positive e ﬀ ects on the distribution of mangroves. The ﬁndings of this study provide a reference for the management and protection of mangroves, which is of great practical signiﬁcance.


Introduction
Mangrove forests are among the most productive ecosystems on Earth and can be found at tropical and subtropical latitudes [1][2][3]. This ecosystem has received increasing attention from not only ocean and land managers but also academia and conservation communities [4,5] due to its important ecological and societal functions [6,7], such as providing food and breeding areas as well as nursing grounds for many fauna, reducing pollution, precipitating fine sediments, protecting coastlines, and storing carbon, among other functions [8][9][10][11][12][13]. However, mangroves have undergone rapid spatiotemporal As it is often cloudy in the studied coastal area, which is vast and required eight mosaicked Landsat scenes, it was difficult to obtain cloud-free images of all scenes in the same year.

Image Preprocessing
First, radiometric calibration was performed for all images to eliminate the errors produced by variations in the sensor performance and characteristics over time as well as to generate the exoatmospheric reflectance images [40]. Second, atmospheric correction was implemented using the Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) tool to eliminate atmospheric and light effects and to obtain the actual surface reflectance. Third, evenly distributed ground control points (GCPs) were selected to geometrically correct the images collected in 1985, 1995, and 2015 with the reference image collected in 2005. The fourth step was image mosaicking. For this process, the Seamless Mosaic tool was used to mosaic eight scenes for each study year, and the cubic convolution resampling method was selected because it provided a higher spatial accuracy than the nearest-neighbor resampling method [16]. The mosaic image from 1985 is shown in Figure 1b. Finally, we cropped the mosaic images to extract the subset images of coastal zone where mangrove forests grow. The result of cropped mosaic image covering the coastal zone of Guangdong province from 1985 is shown in Figure 1c.
This process improved the efficiency of the analysis by avoiding unnecessary calculations, and it improved the classification accuracy by reducing interference from inland areas. In addition, the final study region included Macao and Hong Kong. The mangroves in Macao and Hong Kong were counted with those in Zhuhai City and Shenzhen City, respectively.

Image Classification
To distinguish mangroves from other types of land use and to further discuss the factors influencing the mangrove distribution, six land cover types were defined: (1) mangroves, (2) water, (3) farmland, (4) built-up land, (5) other vegetation, and (6) other types (e.g., shadow areas, tidal flats, and bare land). These were mainly based on the consideration of the actual conditions, the characteristics identified through high spatial resolution images available on Google Earth, and other related studies [4,38]. Then, the decision tree method was used to retrieve the information about mangrove distribution. The decision tree method has been proven to be effective for mapping mangrove forests, as shown by [27]. The decision tree method based on CART algorithm is known as a machine learning or expert system and provides a non-parametric discriminating statistical relationship between multiple data layers to produce a binary decision tree. It splits the training data pixels into more and more homogeneous subsets through recursion and measures homogeneity according to classes defined by training data [24]. One advantage of this method is that it can incorporate multiple sources of data to define classification rules.
The main steps of the classification process, which are shown in Figure 2, include the following: (1) a multi-band image was composited to retrieve mangrove information; (2) the datasets used for classification were incorporated; (3) the training samples for each type of land use were selected with the Region of Interest (ROI) selection tool in ENVI (The Environment for Visualizing Images) 5.3. The false color composites combined with bands 3, 4, and 5 for TM and 4, 5, and 6 for OLI were used to assist the selection of training samples by referring to an established interpretation key ( Table 1). The samples shown as polygons that were selected for training the classification rules are listed in Table 2. (4) The features of training samples were obtained and the CART algorithm was used to discover the relevant classification rules for each defined type of land use with the extension tool RuleGen of ENVI, which can analyze the ROI samples and generate the ENVI decision tree project file automatically; and (5) the rules were applied to stacked images to retrieve the classification results.
In this study, the original multi-spectral images; Normalized Difference Vegetation Index (NDVI) data, which is expressed as Formula (1); Inundated Mangrove Forest Index (IMFI) data, as proposed by Jia [41] and calculated using Formula (2); and the initial classification imagery obtained from ISODATA, were added as bands to improve the ROI separability. The NDVI, which ranges from −1 to 1, is an indicator of the health and density of vegetation, and a higher density of green leaves results in a value close to 1 [42]. The IMFI can help to distinguish between submerged mangroves and sea water. Evenly distributed samples were selected to create decision trees, and the ROI separability values between different samples were greater than 1.8, which indicated that the samples were suitable.
where R red , R blue , R green and R nir represent the reflectance of the red band, blue band, green band, and near infrared band, respectively.  Figure 2. Flow chart of the classification process by using the Classification and Regression Tree (CART) algorithm.  Figure 2. Flow chart of the classification process by using the Classification and Regression Tree (CART) algorithm.  Figure 2. Flow chart of the classification process by using the Classification and Regression Tree (CART) algorithm.      Figure 2. Flow chart of the classification process by using the Classification and Regression Tree (CART) algorithm.

Accuracy Validation
After classification, majority analysis was implemented, and further editing was performed to eliminate obvious errors. In addition, the accuracy of the classification results was further validated. To measure the quality of the classification, accuracy measures were used in this study. Generally, these measures were based on comparing a classification result with a reference classification obtained via a different method, e.g., on-site ground measurements, which are considered reliable and true. A confusion matrix that contains all the information about the relation between the classification and reference classification can be derived. Then, the overall accuracy (OA), producer's accuracy (PA), and user's accuracy (UA) can be calculated from the confusion matrix. The OA is defined as the proportion of all reference pixels that are classified correctly (in the sense that the class assignment of the classification and of the reference classification agree). As the OA, PA, and UA give no information about what classes are classified with high accuracy and cannot reflect the agreement between the two classification sets for each class, the kappa coefficient was further used to reflect this agreement.
Typically, high spatial resolution images from Google Earth can be regarded as ground truth data and have frequently been used to evaluate classification accuracy [43]. In this study, due to the difficulty in obtaining the ground truth data of mangrove distribution over such a large scope, high spatial resolution images derived from Google Earth were also regarded as the true references to calculate the OA, PA, UA, and kappa coefficient. The mangrove appearance on high spatial resolution images was confirmed by referring to the locations of mangrove reserves. In addition, accuracy validation was performed only for the classification result obtained in 2015 because of the difficulty in acquiring high-resolution images in previous study years. Evenly distributed samples were randomly selected from the Google Earth images collected in 2015 to calculate the indices. One hundred samples were selected for every land use type, except for the mangrove samples in Eastern Guangdong, where only 53 samples were selected, due to the small distribution of mangrove forest, and in Western Guangdong, Pearl River Delta, and Eastern Guangdong, where only 30 samples were selected for other types, due to the same reason. The distribution of these validation samples is shown in Figure 3. Image preprocessing and accuracy validation were performed in ENVI 5.3, and image classification was mainly implemented in ENVI 5.3, except for some editing that was done in ArcGIS 10.3.

Landscape Pattern of Mangrove Distribution
Furthermore, this study also focused on the spatiotemporal patterns of the mangrove distribution during the study period. Measures such as landscape metrics, which can reflect some ecological significance for land use change, have been adopted to discuss the spatiotemporal patterns of mangrove distribution based on the classification results. A total of six landscape metrics at the class level were selected in this study. These metrics are: (1) the mean patch area (AREA_MN); (2) the mean Perimeter-Area ratio (PARA_MN); (3) the mean shape index (SHAPE_MN); (4) the number of patches (NP); (5) the patch density (PD); and (6) the patch cohesion index (COHESION). Therein, AREA_MN and NP are size metrics, and the patch area is the basis for the calculation of many other metrics from the patch level to the class and landscape levels. PARA_MN and SHAPE_MN can measure the shape complexity of the patches. The higher the value of PARA_MN is, the more complex the patch shape is. The values of SHAPE_MN are usually larger than or equal to 1. When almost all the patch shapes are close to square-shaped, the value of SHAPE_MN is equal to 1, and when the patch shapes are more irregular, the value of SHAPE_MN will correspondingly increase. PD can reflect the spatial uniformity of mangrove distribution. Under the condition of keeping the total landscape area unchanged, the information reflected by NP and PD is consistent. For the COHESION, the values of it range from 0 to 100, and it can measure the connectivity between patches of the same landscape type, and larger values symbolize stronger connectivity between patches. All these metrics were calculated using Fragstats v4.2.1.

Results
The classifications of mangrove forests and other types of land use were derived from Landsat images based on the CART algorithm. Based on the selected training samples, classification rules were retrieved from the RuleGen.

Accuracy of the Classification Results
The classification results of the above defined land use types in 2015 are shown in Figure 4. The accuracy indices in 2015 are shown in Table 3. The producer's accuracy of the mangrove classification was around 93%, 83%, and 81% in Western Guangdong (WG), Pearl River Delta (PRD), and Eastern Guangdong (EG), respectively, and the user's accuracy was 88%, 97%, and 90%, respectively. Thus, the mangrove forests were extracted with acceptable accuracy. The OA was greater than 90% in the three regions, and the kappa coefficients of the classification results were between 0.81 and 1, indicating excellent to perfect agreement [44]. Users' accuracies and producers' accuracies of other types obtained almost acceptable values. According to the confusion matrices, the accuracies of classification were acceptable, suggesting that further analysis could be performed based on the classification results.

Spatial Pattern of Mangrove Distribution
The distribution of mangroves over the past thirty years is shown in Figure 5 in detail. The classification results indicate that the distribution of mangrove forests has displayed similar patterns over the past thirty years. Mangrove forests are mainly distributed in Western Guangdong, especially in Zhanjiang City. Additionally, there are a few mangrove forests in Eastern Guangdong.

Spatial Pattern of Mangrove Distribution
The distribution of mangroves over the past thirty years is shown in Figure 5 in detail. The classification results indicate that the distribution of mangrove forests has displayed similar patterns over the past thirty years. Mangrove forests are mainly distributed in Western Guangdong, especially in Zhanjiang City. Additionally, there are a few mangrove forests in Eastern Guangdong. From Figure 5, it can be seen that in 1985, mangroves were mainly distributed in Zhanjiang city, Leizhou city, and Zhuhai city, which represented 45% of mangroves in Guangdong province. The areas of the mangroves in these three cities were 2199, 1054, and 895 ha, respectively, which represented about 24%, 11%, and 10% of the mangroves in Guangdong province, respectively. The mangrove forests in Zhanjiang city were mainly located in Leizhou bay, especially in the marine From Figure 5, it can be seen that in 1985, mangroves were mainly distributed in Zhanjiang city, Leizhou city, and Zhuhai city, which represented 45% of mangroves in Guangdong province. The areas of the mangroves in these three cities were 2199, 1054, and 895 ha, respectively, which represented about 24%, 11%, and 10% of the mangroves in Guangdong province, respectively. The mangrove forests in Zhanjiang city were mainly located in Leizhou bay, especially in the marine ecological nature reserve of Leizhou bay, where many silty coastlines and wide mudflats are distributed, and the land around this area is dominated by aquaculture. Mangroves in Zhuhai city were mainly distributed in the bay area and around islands, for example, Hengqin Island and Qi'ao Island. The main types of land use around mangrove growth areas were forest, aquaculture, and agricultural land.
In 1995, mangroves were mainly distributed in Zhanjiang city, Leizhou city, and Suixi County, which represented about 57% of the mangroves in Guangdong province. The mangrove areas in these three districts were about 3313, 1193, and 986 ha, representing about 35%, 12%, and 10% of the mangrove area in Guangdong province, respectively. Mangroves in Suixi County were mainly distributed in estuarine areas where the main types of land use were aquaculture and agricultural land. The situations in Zhanjiang city and Leizhou city were similar to those in 1985.
In 2005, mangroves were mainly distributed in Zhanjiang city, Leizhou city, and Lianjiang city, which represented about 44% of the mangroves in Guangdong province. The mangrove areas in these three districts were about 1407, 898, and 701 ha, representing about 21%, 13%, and 10% of the mangrove area in Guangdong province, respectively. Mangroves in Lianjiang city were mainly distributed in Zhanjiang mangrove national nature reserve, which was established in 1990; other estuarine areas also had large distributions. The main types of land use around mangrove growth areas were aquaculture and agricultural land. The situations in Zhanjiang city and Leizhou city were similar to those in 1985, except that more aquaculture existed around mangroves.
In 2015, mangroves were mainly distributed in Leizhou city, Zhanjiang city, and Lianjiang city, which represented about 50% of the mangroves in Guangdong province. The mangrove areas in these three districts were about 1840, 1502, and 1501 ha, representing about 19%, 15%, and 15% of the mangrove area in Guangdong province, respectively. The situations in Zhanjiang city, Leizhou city, and Lianjiang city were similar to those in 2005, except for the existence of wider mudflats around the mangroves.

Spatiotemporal Variations in the Mangrove Distribution
The variation in the total area of mangrove forest over the past thirty years is shown in Figure 6, and the spatial variation in mangroves from 1985 to 2015 is shown in Figure 7.
Remote Sens. 2019, 11, x FOR PEER REVIEW 10 of 20 ecological nature reserve of Leizhou bay, where many silty coastlines and wide mudflats are distributed, and the land around this area is dominated by aquaculture. Mangroves in Zhuhai city were mainly distributed in the bay area and around islands, for example, Hengqin Island and Qi'ao Island. The main types of land use around mangrove growth areas were forest, aquaculture, and agricultural land.
In 1995, mangroves were mainly distributed in Zhanjiang city, Leizhou city, and Suixi County, which represented about 57% of the mangroves in Guangdong province. The mangrove areas in these three districts were about 3313, 1193, and 986 ha, representing about 35%, 12%, and 10% of the mangrove area in Guangdong province, respectively. Mangroves in Suixi County were mainly distributed in estuarine areas where the main types of land use were aquaculture and agricultural land. The situations in Zhanjiang city and Leizhou city were similar to those in 1985.
In 2005, mangroves were mainly distributed in Zhanjiang city, Leizhou city, and Lianjiang city, which represented about 44% of the mangroves in Guangdong province. The mangrove areas in these three districts were about 1407, 898, and 701 ha, representing about 21%, 13%, and 10% of the mangrove area in Guangdong province, respectively. Mangroves in Lianjiang city were mainly distributed in Zhanjiang mangrove national nature reserve, which was established in 1990; other estuarine areas also had large distributions. The main types of land use around mangrove growth areas were aquaculture and agricultural land. The situations in Zhanjiang city and Leizhou city were similar to those in 1985, except that more aquaculture existed around mangroves.
In 2015, mangroves were mainly distributed in Leizhou city, Zhanjiang city, and Lianjiang city, which represented about 50% of the mangroves in Guangdong province. The mangrove areas in these three districts were about 1840, 1502, and 1501 ha, representing about 19%, 15%, and 15% of the mangrove area in Guangdong province, respectively. The situations in Zhanjiang city, Leizhou city, and Lianjiang city were similar to those in 2005, except for the existence of wider mudflats around the mangroves.

Spatiotemporal Variations in the Mangrove Distribution
The variation in the total area of mangrove forest over the past thirty years is shown in Figure 6, and the spatial variation in mangroves from 1985 to 2015 is shown in Figure 7.  areas exhibited mangrove loss because of sea reclamation (subgraphs D and F) and aquatic breeding (subgraph B). In PRD, the mangrove area also increased in many zones, for example, substantial mangrove increases occurred in the nature reserve on Qi'ao island (subgraph C) and at Mai Po Inner Deep Bay Ramsar Site (subgraph G). In EG, increases and decreases in mangroves were both rare due to the limited distribution of mangroves in the region. Mangrove loss always occurred landward, and increases often occurred seaward (Figure 7). In addition, the spatial pattern of mangrove distribution transformed from fragmented to aggregated. In detail, the mangrove estimates in Western Guangdong in 1985Guangdong in , 1995Guangdong in , 2005, and 2015 were about 5564, 7630, 5004, and 7380 ha, respectively, accounting for about 60%, 80%, 74%, and 76% of the total provincial mangrove area in these years. The mangrove area was extremely small in Eastern Guangdong in the four study years, with values of about 989, 515, 512, and 741 ha in 1985, 1995, 2005, and 2015, respectively, accounting for only about 11%, 5%, 8%, and 8% of the total provincial mangrove area in those years. The mangrove estimates for the Pearl River Delta were between those for Western Guangdong and Eastern Guangdong, with values of about 30%, 15%, 19%, and 16% in 1985, 1995, 2005, and 2015, respectively. Figure 8 shows the variations in the mangrove distribution at the county and city levels. Although there was a sharp decrease in the mangrove area from 1995 to 2005 in the province, remarkable increases of about 160, 86, and 59 ha were observed in Wuchuan City, Shenzhen City, and Yangxi County, respectively. From 2005 to 2015, a significant increase in mangroves occurred in the province, but notable decreases of about 83, 46, and 45 ha were found in Taishan City, Jiangmen City, and Wuchuan City, respectively. Thus, changes in the province did not represent regional trends. Similarly, long-term changes differed from short-term changes. For example, the mangrove  Figure 7 shows, a substantial increase in mangroves occurred in the national nature reserve in Zhanjiang City, which was established in 1990. Different from Western Guangdong, the Pearl River Delta and Eastern Guangdong showed net decreases of about 1173 and 248 ha, respectively, from 1985 to 2015. Most areas exhibited mangrove loss because of sea reclamation (subgraphs D and F) and aquatic breeding (subgraph B). In PRD, the mangrove area also increased in many zones, for example, substantial mangrove increases occurred in the nature reserve on Qi'ao island (subgraph C) and at Mai Po Inner Deep Bay Ramsar Site (subgraph G). In EG, increases and decreases in mangroves were both rare due to the limited distribution of mangroves in the region. Mangrove loss always occurred landward, and increases often occurred seaward (Figure 7). In addition, the spatial pattern of mangrove distribution transformed from fragmented to aggregated.
In detail, the mangrove estimates in Western Guangdong in 1985Guangdong in , 1995Guangdong in , 2005, and 2015 were about 5564, 7630, 5004, and 7380 ha, respectively, accounting for about 60%, 80%, 74%, and 76% of the total provincial mangrove area in these years. The mangrove area was extremely small in Eastern Guangdong in the four study years, with values of about 989, 515, 512, and 741 ha in 1985, 1995, 2005, and 2015, respectively, accounting for only about 11%, 5%, 8%, and 8% of the total provincial mangrove area in those years. The mangrove estimates for the Pearl River Delta were between those for Western Guangdong and Eastern Guangdong, with values of about 30%, 15%, 19%, and 16% in 1985, 1995, 2005, and 2015, respectively. Figure 8 shows the variations in the mangrove distribution at the county and city levels. Although there was a sharp decrease in the mangrove area from 1995 to 2005 in the province, remarkable increases of about 160, 86, and 59 ha were observed in Wuchuan City, Shenzhen City, and Yangxi County, respectively. From 2005 to 2015, a significant increase in mangroves occurred in the province, but notable decreases of about 83, 46, and 45 ha were found in Taishan City, Jiangmen City, and Wuchuan City, respectively. Thus, changes in the province did not represent regional trends. Similarly, long-term changes differed from short-term changes. For example, the mangrove area in Suixi County increased by about 24 ha from 1985 to 2015 but decreased by about 475 ha from 1995 to 2005. Additionally, the mangrove area decreased by approximately 697 ha in Zhanjiang City from 1985 to 2015 but increased by 1115 ha from 1985 to 1995. In addition, a decrease of 328 ha occurred in Zhuhai City from 1985 to 2015, but an increase of 222 ha was observed from 2005 to 2015. As a result, appropriate temporal and spatial scales should be selected according to the study objectives. As a result, appropriate temporal and spatial scales should be selected according to the study objectives.

Discussion
The total area and overall variation trend of mangrove distribution were close to those reported in other studies [38]. For example, the mangrove area in 2015 estimated in this study is similar to the estimation of 8136 ha provided in publication [35], and the mature mangrove area of Guangdong Province in 2001 was reported to be 9084 ha [45], which is between the areas of 9556 ha in 1995 and 6793 ha in 2005 derived in this study. As there are differences in the resolution of remote sensing imageries, classification methods, and the collection times of the images, there is a certain discrepancy in the classification results of mangrove forests. In the study area, the length of the coastline in Western Guangdong is the longest, followed by those in the Pearl River Delta and Eastern Guangdong. These differences may be why Western Guangdong is home to the majority of mangroves in Guangdong Province and why Eastern Guangdong had few mangrove forests compared to the other regions during the study period.

Distribution Characteristics of Mangrove Patches
Further analysis was carried out to discuss the distribution characteristics of mangrove patches; the results are shown in Figure 9. The Mean Patch Area showed a decreasing tendency from 1985 to 2005, while an increasing tendency was shown from 2005 to 2015. The change trend in the Mean Shape Index is consistent with that of Mean Patch Area. This means that the shapes of mangrove

Discussion
The total area and overall variation trend of mangrove distribution were close to those reported in other studies [38]. For example, the mangrove area in 2015 estimated in this study is similar to the estimation of 8136 ha provided in publication [35], and the mature mangrove area of Guangdong Province in 2001 was reported to be 9084 ha [45], which is between the areas of 9556 ha in 1995 and 6793 ha in 2005 derived in this study. As there are differences in the resolution of remote sensing imageries, classification methods, and the collection times of the images, there is a certain discrepancy in the classification results of mangrove forests. In the study area, the length of the coastline in Western Guangdong is the longest, followed by those in the Pearl River Delta and Eastern Guangdong. These differences may be why Western Guangdong is home to the majority of mangroves in Guangdong Province and why Eastern Guangdong had few mangrove forests compared to the other regions during the study period.

Distribution Characteristics of Mangrove Patches
Further analysis was carried out to discuss the distribution characteristics of mangrove patches; the results are shown in Figure 9.  This spatial pattern of mangroves is actually dominated by characteristics of the distribution locations and coastal lines, and the increased area of mangroves was mainly due to artificial planting. It has been highlighted that biodiversity conservation is closely related to mangrove fragmentation [4]. This is of great significance for the formation of a compact pattern to expand the size of mangrove patches. This requires to further understand the difference between natural growing mangroves and artificial planting mangroves well, which can help improve the landscape pattern of mangrove patches to provide better ecological services. High spatial resolution images can help retrieve such information.

Factors Influencing Spatiotemporal Variations in Mangrove Forests
The variation in the area of mangrove forest is mainly influenced by natural factors and This spatial pattern of mangroves is actually dominated by characteristics of the distribution locations and coastal lines, and the increased area of mangroves was mainly due to artificial planting. It has been highlighted that biodiversity conservation is closely related to mangrove fragmentation [4]. This is of great significance for the formation of a compact pattern to expand the size of mangrove patches. This requires to further understand the difference between natural growing mangroves and artificial planting mangroves well, which can help improve the landscape pattern of mangrove patches to provide better ecological services. High spatial resolution images can help retrieve such information.

Factors Influencing Spatiotemporal Variations in Mangrove Forests
The variation in the area of mangrove forest is mainly influenced by natural factors and anthropogenic activities [28]. From the perspective of natural influence, temperature variation and sea level rise, which are also important indicators for measuring climate change, have been generally discussed. The related research has reported that the temperature has a macroscopic effect on the distribution of mangrove forests. In the northern hemisphere, the mangrove forest area, number of tree species, and tree height decrease with increasing latitude [46]. The results about the distribution pattern of mangroves among Eastern Guangdong, Pearl River Delta, and Western Guangdong also validated this evidence. Meanwhile, according to the results of the national wetland survey in 2001, 90% of the mangroves in Guangdong province are sensitive to sea level rise since they are in front of the seawalls [47]. It was reported that the relative sea level rise rate along the coast of Guangdong was 2-3 mm/a from 1986 to 2008, and the sea level will be 38-45 cm higher at the end of the 21st century than in 2000 [48]. The continuing variation trend will result in the serious disturbances to the mangrove ecosystem including a decrease in the total area and degeneration of the ecological quality over a very long period of time. Actually, the influence of the sea level rise requires further validation using the data materials collected over a very long period of time. On the other hand, Guangdong province, located along the South China Sea, suffers frequent typhoon attacks. Extreme weather events also play an important role in the distribution of mangroves; although mangroves can weaken wind and waves, they also experience severe disturbances [5]. In a study about Shenzhen city, it was reported that the total litter fall from mangroves increases after the passage of a typhoon. Typhoons can damage the mangrove ecosystem when the wind force is higher than the level of category 11, and the destruction depends on the origin, age, and density of the mangrove forests [49,50].
However, anthropogenic activities have had drastic influences on the spatiotemporal variations of the mangrove area in coastal Guangdong since the reform and opening-up, and this influence on the mangrove ecosystem may exceed that of natural factors. It was reported that, during the 1980s and 1990s, plenty of mangroves in Guangdong were degenerated into aquaculture area [51]. At the provincial level, a sharp decrease in the mangrove area was observed from 1995 to 2005, especially in Zhanjiang City, Suixi County, Leizhou City, Zhuhai City, and Dianbai County, with decreases of 1906, 475, 295, 102, and 69 ha, respectively. The majority of these decreases were due to urbanization and aquaculture development. Table 4 shows some examples of mangrove erosion in Zhuhai, Macao, and Zhanjiang from 1995 to 2005. Some causes of erosion are unclear, but we speculate that erosion may be caused by deforestation, urbanization or sea level rise in the study area. Further research on this topic is needed.
Anthropogenic activities may have a positive effect on the mangrove distribution, such as by establishing mangrove nature reserves. The first reserve of 815 ha was established in Shenzhen City in 1984, and the total area of mangrove nature reserves sharply increased to 27,489 ha in 1995 before steadily but continuously increasing from 1995 to 2015, as shown in Figure 10. Detailed information corresponding to the mangrove nature reserve in Guangdong is listed in Table 5.
Since the establishment of those nature reserves, some of them have achieved enormous success in protecting mangroves [29]. For example, the mangrove area increased by 1115 ha in Zhanjiang City from 1985 to 1995 because of the establishment of national nature reserves in 1990, and that in Shenzhen City increased by 86 ha from 2005 to 2015 because of the establishment of city nature reserves in 2010. The effect of some nature reserves may be neutralized by other factors, such as urbanization. For example, a decrease of 448 ha occurred in Zhuhai City from 1985 to 1995, even though a provincial nature reserve was established in 1989. However, some nature reserves had no obvious effect on the regional mangrove area; potentially, these reserves were lower-level types so that they did not have sufficient budget to support the conservation [38]. Further studies are necessary to confirm these evidences.

Region 1995 2005 Erosion Reason
Zhuhai City of the mangrove area in coastal Guangdong since the reform and opening-up, and this influence on the mangrove ecosystem may exceed that of natural factors. It was reported that, during the 1980s and 1990s, plenty of mangroves in Guangdong were degenerated into aquaculture area [51]. At the provincial level, a sharp decrease in the mangrove area was observed from 1995 to 2005, especially in Zhanjiang City, Suixi County, Leizhou City, Zhuhai City, and Dianbai County, with decreases of 1906, 475, 295, 102, and 69 ha, respectively. The majority of these decreases were due to urbanization and aquaculture development. Table 4 shows some examples of mangrove erosion in Zhuhai, Macao, and Zhanjiang from 1995 to 2005. Some causes of erosion are unclear, but we speculate that erosion may be caused by deforestation, urbanization or sea level rise in the study area. Further research on this topic is needed.

Region 1995 2005 Erosion Reason
Zhuhai City aquaculture unclear of the mangrove area in coastal Guangdong since the reform and opening-up, and this influence on the mangrove ecosystem may exceed that of natural factors. It was reported that, during the 1980s and 1990s, plenty of mangroves in Guangdong were degenerated into aquaculture area [51]. At the provincial level, a sharp decrease in the mangrove area was observed from 1995 to 2005, especially in Zhanjiang City, Suixi County, Leizhou City, Zhuhai City, and Dianbai County, with decreases of 1906, 475, 295, 102, and 69 ha, respectively. The majority of these decreases were due to urbanization and aquaculture development. Table 4 shows some examples of mangrove erosion in Zhuhai, Macao, and Zhanjiang from 1995 to 2005. Some causes of erosion are unclear, but we speculate that erosion may be caused by deforestation, urbanization or sea level rise in the study area. Further research on this topic is needed.

Region 1995 2005 Erosion Reason
Zhuhai City aquaculture unclear aquaculture the mangrove ecosystem may exceed that of natural factors. It was reported that, during the 1980s and 1990s, plenty of mangroves in Guangdong were degenerated into aquaculture area [51]. At the provincial level, a sharp decrease in the mangrove area was observed from 1995 to 2005, especially in Zhanjiang City, Suixi County, Leizhou City, Zhuhai City, and Dianbai County, with decreases of 1906, 475, 295, 102, and 69 ha, respectively. The majority of these decreases were due to urbanization and aquaculture development. Table 4 shows some examples of mangrove erosion in Zhuhai, Macao, and Zhanjiang from 1995 to 2005. Some causes of erosion are unclear, but we speculate that erosion may be caused by deforestation, urbanization or sea level rise in the study area. Further research on this topic is needed.

Region 1995 2005 Erosion Reason
Zhuhai City aquaculture unclear the mangrove ecosystem may exceed that of natural factors. It was reported that, during the 1980s and 1990s, plenty of mangroves in Guangdong were degenerated into aquaculture area [51]. At the provincial level, a sharp decrease in the mangrove area was observed from 1995 to 2005, especially in Zhanjiang City, Suixi County, Leizhou City, Zhuhai City, and Dianbai County, with decreases of 1906, 475, 295, 102, and 69 ha, respectively. The majority of these decreases were due to urbanization and aquaculture development. Table 4 shows some examples of mangrove erosion in Zhuhai, Macao, and Zhanjiang from 1995 to 2005. Some causes of erosion are unclear, but we speculate that erosion may be caused by deforestation, urbanization or sea level rise in the study area. Further research on this topic is needed. Anthropogenic activities may have a positive effect on the mangrove distribution, such as by establishing mangrove nature reserves. The first reserve of 815 ha was established in Shenzhen City in 1984, and the total area of mangrove nature reserves sharply increased to 27,489 ha in 1995 before steadily but continuously increasing from 1995 to 2015, as shown in Figure 10. Detailed information corresponding to the mangrove nature reserve in Guangdong is listed in Table 5.

Macao urbanization
Zhanjiang City aquaculture Anthropogenic activities may have a positive effect on the mangrove distribution, such as by establishing mangrove nature reserves. The first reserve of 815 ha was established in Shenzhen City in 1984, and the total area of mangrove nature reserves sharply increased to 27,489 ha in 1995 before steadily but continuously increasing from 1995 to 2015, as shown in Figure 10. Detailed information corresponding to the mangrove nature reserve in Guangdong is listed in Table 5.   Anthropogenic activities may have a positive effect on the mangrove distribution, such as by establishing mangrove nature reserves. The first reserve of 815 ha was established in Shenzhen City in 1984, and the total area of mangrove nature reserves sharply increased to 27,489 ha in 1995 before steadily but continuously increasing from 1995 to 2015, as shown in Figure 10. Detailed information corresponding to the mangrove nature reserve in Guangdong is listed in Table 5. Since the establishment of those nature reserves, some of them have achieved enormous success in protecting mangroves [29]. For example, the mangrove area increased by 1115 ha in Zhanjiang City from 1985 to 1995 because of the establishment of national nature reserves in 1990, and that in Shenzhen City increased by 86 ha from 2005 to 2015 because of the establishment of city nature reserves in 2010. The effect of some nature reserves may be neutralized by other factors, such as urbanization. For example, a decrease of 448 ha occurred in Zhuhai City from 1985 to 1995, even  Actually, there are a variety of factors that affect mangrove distribution. Future studies will explore the dominant factors to provide referable information for mangrove preservation.

Conclusions
In this study, the spatiotemporal patterns of mangrove distribution as well as its variation were detected from multi-temporal remote sensing images. The classification results indicated that the incorporation of NDVI, IMFI, and ISODATA data can improve the discrimination between mangrove forests and other land use types. The accuracy validation with high spatial resolution images indicated that the classification method used in this study yielded good accuracy, with overall accuracy and kappa coefficient higher than 90% and 0.8, respectively.
It can be concluded that mangrove distribution in the coastal zone of Guangdong province was shown to have a significant spatial difference and periodical evolution pattern, the total area of mangrove forests has had a net increase over the last 30 years. Additionally, the majority of mangroves in Guangdong province are distributed in Western Guangdong, while Eastern Guangdong has the smallest area of mangroves. Mangrove patches had stronger connectivity in 2015 than in 1985, the distribution pattern of mangrove patches still presented a fragmented trend from 1985 to 2015. It is necessary to recover the spatial distribution of mangroves with appropriate patterns during the recoverable period. The dynamic changes detected in this study are of great significance to mangrove conservation.
Furthermore, natural and anthropogenic factors influencing the spatial distribution of mangroves were discussed in this study. Natural factors have macroscopic effects on the distribution of mangrove forests, whereas anthropogenic activity may have a negative effect on the distribution of mangroves due to processes such as deforestation, urbanization, and aquaculture development. Some anthropogenic activities, such as the establishment of nature reserves, can promote an increase in the mangrove area. This study showed that enormous success in protecting mangroves has been achieved since the establishment of nature reserves. The findings of this study can provide referable information for coastal management particularly ecological protection.