Dynamic Analysis of Mangrove Forests Based on an Optimal Segmentation Scale Model and Multi-Seasonal Images in Quanzhou Bay , China

Mangrove forests are important coastal ecosystems and are crucial for the equilibrium of the global carbon cycle. Monitoring and mapping of mangrove forests are essential for framing knowledge-based conservation policies and funding decisions by governments and managers. The purpose of this study was to monitor mangrove forest dynamics in the Quanzhou Bay Estuary Wetland Nature Reserve. To achieve this goal, we compared and analyzed the spectral discrimination among mangrove forests, mudflats and Spartina using multi-seasonal Landsat images from 1990, 1997, 2005, 2010, and 2017. We identified the spatio-temporal distribution of mangrove forests by combining an optimal segmentation scale model based on object-oriented classification, decision tree and visual interpretation. In addition, mangrove forest dynamics were determined by combining the annual land change area, centroid migration and overlay analysis. The results showed that there were advantages in the approaches used in this study for monitoring mangrove forests. From 1990 to 2017, the extent of mangrove forests increased by 2.48 km2, which was mostly converted from mudflats and Spartina. Environmental threats including climate change and sea-level rise, aquaculture development and Spartina invasion, pose potential and direct threats to the existence and expansion of mangrove forests. However, the implementation of reforestation projects and Spartina control plays a substantial role in the expansion of mangrove forests. It has been demonstrated that conservation activities can be beneficial for the restoration and succession of mangrove forests. This study provides an example of how the application of an optimal segmentation scale model and multi-seasonal images to mangrove forest monitoring can facilitate government policies that ensure the effective protection of mangrove forests.


Introduction
Mangrove forests, situated in the intertidal zone of tropical and sub-tropical coastal regions, are one of the most productive ecosystems on Earth.Their distinct marine and terrestrial characteristics [1][2][3] provide significant ecological services and functions in terms of coastal water purification, biodiversity conservation, shoreline stabilization, storm protection and fishery harvest [4][5][6].Moreover, as an important component of blue carbon [7], mangrove forests can sequester carbon in aboveground biomass [8], below-ground biomass [9,10], and in sediments [11,12], which is vital for the equilibrium of the global carbon cycle.Despite an abundant range of economic and ecological values, a third of mangrove forests worldwide have been lost in the last fifty years because of rapid urban growth, increasing population pressure, aquaculture expansion, and other impacts caused by anthropogenic disturbances and climate change [13][14][15][16].However, there is an emerging demand for conservation and restoration efforts.To make appropriate decisions and polices, the spatio-temporal extent of mangrove forests needs to be inventoried and monitored, and the driving factors of change need to be identified [17,18].
Mangrove forests are difficult to monitor because of their inaccessibility and large area, which makes field observation problematic.A considerable number of remote sensing monitoring and mapping studies have been conducted, from the local to the global scale, over the past two decades [19,20].Various remote sensing-based methodologies, either exclusively or in combination, have been applied to monitor the extent of change in mangrove forests.Traditional mapping approaches, including visual interpretation and on-screen digitization, have been used to map mangrove forests.Applying visual interpretation as an auxiliary support improves the accuracy of classification results obtained by supervised and unsupervised classification methodologies [21].A variety of classification methods have been investigated and compared to enhance the effectiveness of spectral discrimination for mangrove forests, such as the support vector machine [22], artificial neural networks [23], maximum likelihood classifier [24,25], machine learning [26], and the iterative self-organizing data analysis technique (ISODATA) algorithm [27].Nevertheless, due to "salt-and-pepper" effects, pixel-based classification methods frequently generate erroneous classifications of pixels [28].In contrast, the object-oriented classification method segments an image into groups of contiguous and homogeneous pixels (image objects) as the mapping unit [29], which can reduce the "within-class" spectral variation and effectively overcome the salt-and-pepper effect.In addition, these classification methods consider not only the spectral properties of the objects, but also their texture, shape, and geometric features in the classification process, and as a result, more accurate and effective performances are obtained than with pixel-based approaches [30][31][32].
Image segmentation is deemed to be a critical prerequisite for object-oriented classification, because its quality largely affects the final performance of geo-object recognition.Understanding how to effectively determine the optimal segmentation scale is crucial to the improvement of segmentation quality [33].To date, visual inspection has mainly been used to assess the accuracy of segmentation results.For example, Liu et al. (2017) determined the advisable segmentation scale of Google Earth images based on multiple tests and visual analysis [31].Yet, visual inspection, as a qualitative approach, cannot provide a quantitative evaluation and may suffer from subjectivity, since different people are likely to have different opinions on which segmentation result is better [33,34].Quantitative evaluation approaches can score a range of segmentation results.However, there are few studies focused on optimal segmentation scale selection in terms of monitoring and mapping for mangrove forests.
Given the importance of mangrove forests, as well as the limitations of segmentation scale evaluation and selection methods in mangrove forest mapping, the objectives of this study were to: (1) Apply a systematic optimal segmentation scale selection method to assist in object-oriented image classification; (2) monitor the dynamics of mangrove forest extent during the period 1990-2017, and transitions between mangrove forests and other land cover types; (3) analyze the influence of anthropogenic activities, climate change, and plant invasion on the spatio-temporal changes of mangrove forests; and (4) propose more feasible conservation policies and local management plans for mangrove forests.We used the mangrove forests in Quanzhou Bay, Fujian Province, China, as the study site.An optimal segmentation scale model based on object-oriented classification was applied to determine the best segmentation scale.We combined image segmentation with a decision tree and visual interpretation to produce land cover maps.The mangrove forest dynamics and conversion to and from other land cover types were described using annual land change rate (ALCR), centroid migration, and overlay analysis.

Study Area
Quanzhou Bay Estuary Wetland Nature Reserve (QBEWNR) is located in the southeast of Fujian Province, China.The latitude ranges from E118.63 • to E118.71 • and the longitude ranges from N24.79 • to N24.99 • (Figure 1).There are two rivers (i.e., Luoyang River and Jinjiang River) in the area, and the Bay has a regular semi-diurnal tide.The nature reserve runs across five districts and cities (Huian, Luojiang, Fengze, Jinjiang, Shishi) from north to west, with a total area of 7,130 ha.The annual mean temperature is 20.4 • C and the average annual precipitation is 1095.4mm.The prevalent climate of this region is an oceanic monsoon climate, characterized by a warm and wet winter and a hot and rainy summer.The leaf-on season is from April to November, and the leaf-off season is from December to March.
In the nature reserve, the mangrove swamp ecosystem is the fundamental object for protection.The main mangrove forest species include Kandelia candel, Aegiceras corniculatum, Avicennia marina, and Acanthus ilicifolius.In addition, there are two species of Spartina-Spartina anglica and Spartina alterniflora.The mangrove swamps provide habitat for endangered and rare species of aquatic animals and waterfowl (Sousa chinensis, Acipenser sinensis, Larus saundersi, and Egretta eulophotes).They also support extensive fisheries and serve as a natural reservoir that provides vital flood control.Province, China, as the study site.An optimal segmentation scale model based on object-oriented classification was applied to determine the best segmentation scale.We combined image segmentation with a decision tree and visual interpretation to produce land cover maps.The mangrove forest dynamics and conversion to and from other land cover types were described using annual land change rate (ALCR), centroid migration, and overlay analysis.

Study Area
Quanzhou Bay Estuary Wetland Nature Reserve (QBEWNR) is located in the southeast of Fujian Province, China.The latitude ranges from E118.63° to E118.71° and the longitude ranges from N24.79° to N24.99° (Figure 1).There are two rivers (i.e., Luoyang River and Jinjiang River) in the area, and the Bay has a regular semi-diurnal tide.The nature reserve runs across five districts and cities (Huian, Luojiang, Fengze, Jinjiang, Shishi) from north to west, with a total area of 7,130 ha.The annual mean temperature is 20.4 °C and the average annual precipitation is 1095.4mm.The prevalent climate of this region is an oceanic monsoon climate, characterized by a warm and wet winter and a hot and rainy summer.The leaf-on season is from April to November, and the leaf-off season is from December to March.
In the nature reserve, the mangrove swamp ecosystem is the fundamental object for protection.The main mangrove forest species include Kandelia candel, Aegiceras corniculatum, Avicennia marina, and Acanthus ilicifolius.In addition, there are two species of Spartina-Spartina anglica and Spartina alterniflora.The mangrove swamps provide habitat for endangered and rare species of aquatic animals and waterfowl (Sousa chinensis, Acipenser sinensis, Larus saundersi, and Egretta eulophotes).They also support extensive fisheries and serve as a natural reservoir that provides vital flood control.

Data Preparation and Fieldwork
Considering the data consistency, spectral discrimination ability and specific distribution characteristics of mangrove forests in the study area, Landsat Thematic Mapper (TM) and Operational Land Imager (OLI) images (30 m spatial resolution) were selected as the basic data sources with which to analyze the temporal and spatial dynamics of mangrove forests for 1990, 1997, 2005, 2010 and 2017.In addition to mangrove forests, Spartina and mudflats are found within the study area.Previous studies have reported that Spartina expansion has a significant impact on mangrove forest distribution and expansion [31,35].Consequently, mapping the distribution of Spartina enhances our understanding of the factors influencing mangrove forest change.
In the leaf-on season, there are few differences in spectral characteristics between mangrove forests and Spartina.In the leaf-off season, the spectral characteristics of Spartina are similar to those of mudflats.Therefore, it is difficult to distinguish land cover accurately using images from one season.Differentiating mangrove forests, Spartina and mudflats can be improved by using their spectral characteristics in different seasons (Figure 2).

Data Preparation and Fieldwork
Considering the data consistency, spectral discrimination ability and specific distribution characteristics of mangrove forests in the study area, Landsat Thematic Mapper (TM) and Operational Land Imager (OLI) images (30 m spatial resolution) were selected as the basic data sources with which to analyze the temporal and spatial dynamics of mangrove forests for 1990, 1997, 2005, 2010 and 2017.In addition to mangrove forests, Spartina and mudflats are found within the study area.Previous studies have reported that Spartina expansion has a significant impact on mangrove forest distribution and expansion [31,35].Consequently, mapping the distribution of Spartina enhances our understanding of the factors influencing mangrove forest change.
In the leaf-on season, there are few differences in spectral characteristics between mangrove forests and Spartina.In the leaf-off season, the spectral characteristics of Spartina are similar to those of mudflats.Therefore, it is difficult to distinguish land cover accurately using images from one season.Differentiating mangrove forests, Spartina and mudflats can be improved by using their spectral characteristics in different seasons (Figure 2).To account for the influence of tidal phenomena on the results of coastal land cover identification, two seasonal (leaf-on season and leaf-off season) images were selected in each base year, and these were all acquired during low tide.Detailed tidal information for the images is shown in Table 1.All images were geo-rectified with the registration error being less than half a pixel, and they were all atmospherically corrected using Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) [36].From 2015 to 2017, 171 sampling points were obtained via four field surveys (Figure 1), which were used to evaluate the accuracy of the land cover classification results in 2017.Owing to the lack of field survey data in 1990, 1997, 2005 and 2010, 150 independent points for each year were created by a random sampling scheme.These random points were classified into different land cover types (described in Section 2.3) by consulting with experienced interpreters and local experts, and were then used as validation points.The overall accuracy, user accuracy, producer accuracy, and Kappa coefficient were used to assess the accuracy of each land cover map [37].
Data for annual mean temperature and precipitation from 1990 to 2017 were collected at 30 meteorological stations in and around the study area.This allowed the comprehensive and objective analysis of climate change factors driving mangrove forest change.A spatial interpolation method was applied to these meteorological data to obtain a spatially continuous surface.

Land Cover Classification System
A landscape classification system was created, including ten landscape types (Mangrove, Spartina, Mudflat, Water body, Aquaculture pond, Built-up area, Woodland, Cropland, Grassland and Barren land).The details of the classification system and image features for these ten land cover types are given in Table 2.

Optimal Segmentation Scale Model Based on Object-Oriented Classification
To improve the boundary consistency between the segmented image objects and real land cover types, an optimal segmentation scale was applied to process the scale segmentation and selection.It is assumed that when the internal homogeneity (i.e., no mixed objects) and neighboring heterogeneity (i.e., high differentiation among adjacent objects) of segmented objects reaches their optimum, the segmentation scale is optimal.The basic concept of the optimal segmentation scale model was referenced by Chabrier et al. [38] and Espindola et al. [39].Using the model, the relationship between the internal homogeneity and the heterogeneity of segmented objects is expressed by the segmentation quality function, which can objectively evaluate the quality of segmentation effects.In order to obtain the optimal segmentation scale, a series of pre-segmentation experiments for different segmentation scale was performed.In the study, a segmented image series was obtained using 15 different scales, which started with a scale parameter of 10 and ended at 80, increasing by increments of 5.
Apart from the segmentation scale, the image segmentation parameters also included shape and compactness factor, which were determined by referencing previous studies and from practical experience [31].After segmentation, the image objects were classified into specific land cover types using a decision tree and visual interpretation.

Mangrove Forest Change
To analyze and assess the dynamic degree of mangrove forests objectively, the annual land change rate (ALCR) of each land cover type was calculated, which is expressed as follows [40]: where U a and U b represent the area of each land cover type at the beginning and the end of the study period, respectively.T is the number of years.In the study, the time interval was divided into 4 stages: 1990-1997, 1997-2005, 2005-2010 and 2010-2017.
The invasive species of Spartina (Spartina anglica and Spartina alterniflora) pose a significant threat to the existence and expansion of mangrove forests [41].From the perspective of species competition, the adverse effects of Spartina on mangrove forests can be reflected by analyzing the spatial and temporal reproduction intensity of Spartina.An exponential growth model was applied to describe the expansion of Spartina, which is calculated as follows [31]: where S a and S b are defined as the area of Spartina in the beginning year and ending year in each stage, respectively; n is the number of years; and p is the annual expansion rate of Spartina.In addition, to analyze the characteristics of spatial change of mangrove forests more explicitly, an overlay analysis in ArcGIS 10 [42] was used to create a conversion matrix between mangrove forests and other land cover types for the time periods 1990-1997, 1997-2005, 2005-2010 and 2010-2017.Meanwhile, a Sankey diagram [43] was used to illustrate the conversion results of each land cover type, which can visualize the dynamics of total land cover types.

Calculation of Centroid Migration
Centroid migration of different landscapes can reflect their development characteristics in terms of spatial and temporal scales [44].Analysis of the centroid change of mangrove forests and other land cover types not only allows for the quantitative investigation of the spatial and temporal evolution (direction and distance) of mangrove forests, but also the relationship between mangrove forests and other land cover types.For this purpose, the centroid migrations of mangrove forests and Spartina were analyzed.The equation of centroid is expressed as follows [45]: where X t and Y t represent the latitude and longitude coordinates of the centroid for one land cover type in the t year, respectively, X i and Y i are the latitude and longitude coordinates for the centroid of the ith patch for one land cover type, C ti is the area of the ith patch, and n is the number of patches for one land cover type.

Optimal Segmentation Scale
We determined the optimal segmentation scale using the optimal segmentation scale model.The segmentation effect and the segmentation quality of different segmentation scales are shown in Figure 3.According to Figure 3, when the segmentation scale was less than 30, the image was over-segmented and one real geo-object was fragmented by more than one resulting segment.When the segmentation scale was greater than 60, the image was under-segmented and one resulting segment covered an area that did not belong to the target geo-object.Moreover, as the segmentation scale increased, additional time was required to segment.Using qualitative visual analysis, we found that the best segmentation effect occurred when the segmentation scale equaled 40, and the corresponding segmentation quality was 0.95.The optimal segmentation scale model calculated the optimal segmentation scale to be 38 and the corresponding segmentation quality was 0.99.The segmentation quality of the optimal segmentation scale obtained by optimal segmentation scale model was better than that of visual judgment.

Temporal and Spatial Changes of Mangrove Forests
Combining the object-oriented classification, decision tree and visual interpretation, we produced land cover maps of the study area (Figure 4).Table 3 presents the accuracy assessment results of the land cover types in each study year.The overall accuracies of all the classification results were more than 0.91 and all the Kappa coefficients were more than 0.90, which means that our classification results were consistent with those obtained from the validation points.Table 4 shows the change area and ALCR of each land cover type.The comparison in extent of mangrove, Spartina and aquaculture pond are illustrated in Figure 5.The results indicate that the mangrove forests experienced a general increasing trend from 0.09 km 2 in 1990 to 2.57 km 2 in 2017, with a change rate of 0.09 km 2 /y and an ALCR of 102.06%.The mangrove forest area decreased only in the period 1997-2005, while from 2010 to 2017, the mangrove forest area increased rapidly.Spartina exhibited a dramatic increase from 1990 to 2017, implying an average rate of gain of 0.29 km 2 /y.During 1997-2005, Spartina expanded markedly, in contrast to the mangrove forest change trend during this period.The distribution range of Spartina expanded from sporadic patches on the West Luoyang River in 1990 to large-scale continuous stands along both sides of the Luoyang River in 2017.During the period 1990-2017, aquaculture ponds more than tripled in area (from 1.54 to 5.14 km 2 ).Built-up areas increased by 1.76 km 2 from 1990 to 2017.Although there were small changes in area, the ALCR of grassland and barren land were 24.07%and 11.11%, respectively, which were relatively higher than other land cover types from 1990 to 2017.In contrast, mudflats decreased sharply with a loss of 15.75 km 2 during 1990-2017.).Built-up areas increased by 1.76 km 2 from 1990 to 2017.Although there were small changes in area, the ALCR of grassland and barren land were 24.07%and 11.11%, respectively, which were relatively higher than other land cover types from 1990 to 2017.In contrast, mudflats decreased sharply with a loss of 15.75 km 2 during 1990-2017.

Conversions between Mangrove Forests and Other Land Cover Types
Figure 6b and Table 5 illustrate the spatial distribution and area of conversion between mangrove forests and other land cover types.It was found that the conversion region gradually expanded, and the most marked mangrove forest recession occurred from 1997 to 2005, attributed to a large area conversion to mudflats and Spartina, which accounted for 33.51% and 29.74% of mangrove forest reduction in this stage, respectively.During the period 2010-2017, nearly 181 ha of other land cover types were transformed into mangrove forests and resulted in a significant expansion of mangrove forests in this stage.The proportion of mudflats and Spartina converted into mangrove forests were the highest of the different land cover types, accounting for 72.97% and

Conversions between Mangrove Forests and Other Land Cover Types
Figure 6b and Table 5 illustrate the spatial distribution and area of conversion between mangrove forests and other land cover types.It was found that the conversion region gradually expanded, and the most marked mangrove forest recession occurred from 1997 to 2005, attributed to a large area conversion to mudflats and Spartina, which accounted for 33.51% and 29.74% of mangrove forest reduction in this stage, respectively.During the period 2010-2017, nearly 181 ha of other land cover types were transformed into mangrove forests and resulted in a significant expansion of mangrove forests in this stage.The proportion of mudflats and Spartina converted into mangrove forests were the highest of the different land cover types, accounting for 72.97% and 25.48%, respectively.There was little conversion among mangrove forests and cropland, built-up areas or aquaculture ponds during the study period.These results suggest that the change in mangrove forest area can be attributed to their expansion in mudflats and the ecological competition with Spartina.The centroid migration routes of mangrove forests and Spartina are displayed in Figure 7. From 1990 to 2017, the centroid of mangrove forests experienced a trend of shifting landward before 1997 and a process of shifting seaward after 1997.The centroid of Spartina showed a persistent seaward migration throughout the study period.It is noteworthy that the distance between the mangrove forest centroid and the centroid of Spartina was increasing, changing from 159.57m in 1990 to 7308.92 m in 2017.The centroid migration routes of mangrove forests and Spartina are displayed in Figure 7. From 1990 to 2017, the centroid of mangrove forests experienced a trend of shifting landward before 1997 and a process of shifting seaward after 1997.The centroid of Spartina showed a persistent seaward migration throughout the study period.It is noteworthy that the distance between the mangrove forest centroid and the centroid of Spartina was increasing, changing from 159.57m in 1990 to 7308.92 m in 2017.

Mangrove Forest Losses Associated by Human Activities and Possible Environmental Threats
Spartina, as invasive plant species, were first introduced to China's coastal provinces from the Atlantic Coast of the U.S. for coastal environment promotion in 1976 [46].However, due to their rapid growth, high productivity and strong adaptability, introduced Spartina expanded rapidly in tidal zones and now pose a tremendous threat to native biodiversity.Once Spartina have invaded the mangrove forest gaps, they inhibit the germination and growth of mangrove seedlings by competing for space, light and nutrients [47,48].Our results found that the propagation of mangrove forests was suppressed where they grew in or around Spartina (Figure 4).Conversion to Spartina was the primary contributor of mangrove forest losses, especially from 1997 to 2005 (Table 5).Previous studies have also shown Spartina expansion increased the pressure on survival and breeding of mangrove forests [48,49], which is consistent with our findings.Moreover, in other coastal zones of China, for instance, Shankou of Guangxi province, Qi'ao island, and Shenzhen Bay, Spartina invasion has become one of the most serious threats to mangrove forest growth [41,45,50].Spartina can become established in mudflats, easily resulting in the loss of suitable habitats which are available for mangrove forest reproduction [31,51].Our results also showed that the Spartina area increased by 7.84 km 2 with an ALCR of 100.13% from 1990 to 2017, with mudflats as the priority areas of encroachment (Table 4, Figures 4-6).Without effective systematic intervention, we posit that Spartina would gradually cover the entire intertidal mudflat and become an irrevocable threat to mangrove ecosystems.Figure 8a,b illustrates the inhibition of mangrove seedling growth by Spartina and its expansion in mudflats, respectively.
available for mangrove forest reproduction [31,51].Our results also showed that the Spartina area increased by 7.84 km 2 with an ALCR of 100.13% from 1990 to 2017, with mudflats as the priority areas of encroachment (Table 4, Figures 4-6).Without effective systematic intervention, we posit that Spartina would gradually cover the entire intertidal mudflat and become an irrevocable threat to mangrove ecosystems.Figure 8a,b illustrates the inhibition of mangrove seedling growth by Spartina and its expansion in mudflats, respectively.Aquaculture development is another factor that has contributed to the losses of mangrove forest in the QBEWNR.Over the period 1990-2017, aquaculture ponds increased significantly by 3.57 km 2 (Figure 5) and most were established in mudflat areas (Figures 4 and 6).Diseases caused by aquaculture have a negative effect on mangrove forests [52].The untreated water discharged by aquaculture can bring more pressures against the nutrient load in water bodies.Once the self-purification capacity of water bodies has been exceeded, the health of mangrove forests faces Aquaculture development is another factor that has contributed to the losses of mangrove forest in the QBEWNR.Over the period 1990-2017, aquaculture ponds increased significantly by 3.57 km 2 (Figure 5) and most were established in mudflat areas (Figures 4 and 6).Diseases caused by aquaculture have a negative effect on mangrove forests [52].The untreated water discharged by aquaculture can bring more pressures against the nutrient load in water bodies.Once the self-purification capacity of water bodies has been exceeded, the health of mangrove forests faces serious threats [53].In addition, suitable habitat for mangrove forests was encroached upon by constructing aquaculture ponds in mudflats.Consequently, the losses of mangrove forests were due to aquaculture expansion to some extent.In previous studies, researchers have indicated that of the one-third of mangrove forests that have disappeared worldwide in the last 30 years, 35% was lost to aquaculture; this figure may reach 60% by 2030 [54,55].Figure 8c,d displays aquaculture ponds located adjacent to mangrove forests and aquaculture development in mudflats, respectively.
Due to their location in low coastal elevation areas, mangrove forests are particularly vulnerable to sea-level rise [56][57][58].Previous research on sea-level change has reported that the sea-level of the Fujian coastal zone increased at an average rate of 2 mm/y during the last half century [59].In response to rising sea-levels, mangrove forests are more likely to migrate landward [30,60].However, such a landward migration and establishment has been obstructed by artificial seawalls in the study area (Figure 9), which meant that mangrove forests could only expand seaward.In addition, the roots of mangrove forests should be exposed to the air at a certain time, otherwise they would not be able to complete their own physiological processes and propagation [58,61,62].As a result of sea-level rise, lower elevation areas would be inundated and waterlogged longer, which means mangrove forests would not meet the time required for mangrove root respiration, so mangrove forests cannot easily establish and propagate seaward.As shown in Figure 10, the annual average temperature and annual precipitation both increased in the QBEWNR from 1990 to 2017.Although a warmer and wetter climate would favor the establishment of mangrove forests [63], such climate change would promote sea-level rise, leading directly to the space loss of mangrove forest expansion seaward.Therefore, it is evident that, in the context of global warming, the potential for major coastal change aggravates the survival risk and pressures of seaward succession for mangrove forests.Overall, as for the study area, climate change and sea-level rise would potentially have adverse effects on mangrove forests.

Positive Effects of Reforestation Projects and Spartina Control
One notable caveat to these results is that the area of mangrove forests in the QBEWNR was only 0.09 km 2 in 1990, had doubled to 0.18 km 2 in 1997, and then declined to 0.16 km 2 in 2005.These areas represent a limited distribution.Since 2005, mangrove forests have increased rapidly (Figures 5 and 6).Under the ecological restoration schemes for coastal wetlands in Quanzhou Bay, various national and local ecological conservation projects have been carried out since 2000, including "National Wetland Conservation and Restoration", "Marine Environment Ecological Restoration Project in Quanzhou Bay", and "Mangrove Forest Restoration Project in Luoyang River" [64].One of the common objectives for these ecological engineering projects is to cultivate and replant mangrove

Positive Effects of Reforestation Projects and Spartina Control
One notable caveat to these results is that the area of mangrove forests in the QBEWNR was only 0.09 km 2 in 1990, had doubled to 0.18 km 2 in 1997, and then declined to 0.16 km 2 in 2005.These areas represent a limited distribution.Since 2005, mangrove forests have increased rapidly (Figures 5 and 6).Under the ecological restoration schemes for coastal wetlands in Quanzhou Bay, various national and local ecological conservation projects have been carried out since 2000, including "National Wetland Conservation and Restoration", "Marine Environment Ecological Restoration Project in Quanzhou Bay", and "Mangrove Forest Restoration Project in Luoyang River" [64].One of the common objectives for these ecological engineering projects is to cultivate and replant mangrove forests (i.e., reforestation projects) to increase mangrove forest area.According to the previous report [65], the afforestation area of mangrove forests was approximately 200 ha until 2015, indicating an 80%-85% increase in mangrove forests because of the reforestation projects.Studies focusing on the mangrove forest dynamics in China from 1973 to 2015 had indicated that protection and reforestation actions had supported mangrove forest restoration greatly, which is consistent with our results [30]. Figure 8e,f demonstrates a monument of mangrove forest reforestation project and the mangrove seedling cultivation base, respectively.
Although awareness of the need to protect mangrove forests has increased, the threats of Spartina invasion and the importance of Spartina control were recognized relatively late in China [30,66].This probably explains why reforestation engineering was implemented in 2000, but the mangrove forests continued to decrease between 1997 and 2005 (Table 4).Since 2007, both national and local agencies have promoted the implementation and enforcement of Spartina control projects in the QBEWNR [67].Up to 2017, the area of Spartina converted to mangrove forests had increased to 53.42 ha, including 7.20 ha from 2005 to 2010 and 46.22 ha from 2010 to 2017 (Table 5).Meanwhile, the annual expansion rate of Spartina had declined from 20.66% at the beginning of the same period, to 3.64% at the end.Moreover, the centroid distance between Spartina and mangrove forests has increased (Figure 7).Therefore, it can be inferred that, due to the management and control of Spartina, the disturbance caused by Spartina to mangrove forests has been mitigated to some degree.Figure 8g,h show the specific implementation scenario and subsequent effects of Spartina control, respectively.
As mentioned above, whether through mangrove forest reforestation projects or Spartina control, the performance of conservation activities plays a substantial role in the existence and expansion of mangrove forests, which to a large extent neutralizes or offsets the losses of mangrove forests caused by negative factors.

Suggestions for Conserving and Managing Mangrove Forests
Monitoring and dynamic change analysis for mangrove forests are critical for their conservation and management, and to help with formulating and implementing government policies [68,69].The results from this study on the spatio-temporal change of mangrove forests and conversion between mangrove forests and other land cover types (Figures 4-6, Tables 4 and 5) can be used as a guideline for conserving and managing mangrove forests.
First, the implementation of artificial afforestation projects should be continued to increase the area of mangrove forests.Establishing a mangrove forest monitoring system is also indispensable and would allow for the elimination and control of insect pests and effective feedback on all aspects of mangrove forests.Meanwhile, cutting mangrove forests indiscriminately and disposing of rubbish within mangrove forests must be prohibited.Second, prevention and control of Spartina should be further studied, and more feasible measures should be taken to impede the expansion of Spartina in mudflats, with the aim of restoring suitable habitat for mangrove forests.For the zones in which Spartina has been eradicated, breeding conditions should be monitored regularly.Third, strict limitations on aquaculture development should be implemented in the reserve, and the pollution problems caused by aquaculture must be contained.Fourth, in focusing on the important ecological value of mangrove forests, local managers and conservationists should reinforce mangrove forest conservation by teaching local residents about the conservation implications and by disseminating information about the severe threats faced by mangrove forests.

Advantages and Uncertainties of the Methods for Mangrove Forest Monitoring
So far, the primary segmentation scale determination approach is based on visual analysis, which is subjective and cannot effectively avoid over-segmentation and under-segmentation errors [70].The optimal segmentation scale model used in our study not only gave a quantitative result, which is based on objective evaluation, but also kept better boundary consistency between the segmented image objects and real land cover types (Figure 3).Moreover, in the visual interpretation process, the manual modification of objects obtained using the optimal segmentation scale model was less than the objects generated through visual analysis determining the best scale, which improved classification efficiency and reduced the workload of interpreters.
Classification accuracy can be improved by using the spectral difference among different land cover types effectively [71].In this study, the spectral discrimination of mangrove forests, Spartina, and mudflats derived from seasonal change was given a full comparison and analysis, and we found that they were more accurately distinguished than in previous research [21,30,31].Therefore, multi-seasonal images appear to be feasible and desirable data sources for monitoring mangrove forests located in subtropical areas.In tropical zones, owing to very small phenological differences, it may be difficult to discern spectral distinction between mangrove forests and Spartina.Accordingly, multi-seasonal images may not be applicable for mapping mangrove forests growing in tropical regions.
Our analysis focused on mangrove forest dynamics over a long time series (from 1990 to 2017).This is significant for maintaining the spatial coherence of classification maps derived from various sensor images.Due to the spatial uncertainties that arose from the scanning system, wavelength setting, over-pass time, and angular effect, unbiased maps were difficult to produce.In addition, several uncertainties and limitations also existed in the segmentation parameter selection.Besides the segmentation scale, the shape and compactness variables, which respectively balance spectral homogeneity with the shape of the objects and the compactness with smoothness, are essential for controlling the clustering decision process of image objects.Although the optimal segmentation scale was determined objectively, the selection of other segmentation factors chosen by referencing previous research was subjective.

Conclusions
Considering the importance of tidal information for mangrove forest monitoring and the spectral differences among mangrove forests, mudflats, and Spartina, several multi-seasonal Landsat images acquired during low tide were selected as the basic data sources.Combining an optimal segmentation scale model based on object-oriented classification, centroid migration calculations and spatial analysis, we discerned dynamic changes in the mangrove forest and their influencing factors in the QBEWNR from 1990 to 2017.Our results showed that there were some advantages for the approaches used in this study for mangrove forest monitoring, since the classification accuracy and efficiency of land cover map improved.Mangrove forests significantly expanded during the period of 1990-2017: the total area increased by 2.48 km 2 , with a dynamic degree of 102.06%.Most of the expanding mangrove forests transitioned from mudflats and Spartina.Mangrove forest changes were influenced by many factors.Environmental threats, including climate change and sea-level rise, Spartina invasion and aquaculture development exerted negative effects, while reforestation projects and Spartina control had a positive effect.In this study, the role of the latter was greater than the former.We demonstrated that conservation activities benefit the existence and expansion of mangrove forests.These conclusions can be used as a guide for governments and conservationists to make policies and effectively protect and monitor mangrove forests.

Figure 1 .
Figure 1.Location of the study area (QBEWNR represents the Quanzhou Bay Estuary Wetland Nature Reserve).

Figure 1 .
Figure 1.Location of the study area (QBEWNR represents the Quanzhou Bay Estuary Wetland Nature Reserve).

Figure 2 .
Figure 2. Comparisons of image characteristics and spectral reflectance for mangrove forests, Spartina, and mudflats in different seasons.(a) Image of leaf-on season; (b) image of leaf-off season; (c) spectral reflectance curve of leaf-on season for mangrove forests, Spartina, and mudflats; (d) spectral reflectance curve of leaf-off season for mangrove forests, Spartina, and mudflats (A: mangrove forests, B: Spartina, C: mudflats; Band combination Landsat OLI: R: G: B = Band 5: Band 4: Band 3).

Figure 2 .
Figure 2. Comparisons of image characteristics and spectral reflectance for mangrove forests, Spartina, and mudflats in different seasons.(a) Image of leaf-on season; (b) image of leaf-off season; (c) spectral reflectance curve of leaf-on season for mangrove forests, Spartina, and mudflats; (d) spectral reflectance curve of leaf-off season for mangrove forests, Spartina, and mudflats (A: mangrove forests, B: Spartina, C: mudflats; Band combination Landsat OLI: R: G: B = Band 5: Band 4: Band 3).
With bright white, yellowish white or white brown; irregular shape, uniform texture * The example images are the color composite image from Landsat OLI: R: G: B = Band 5: Band 4: Band 3. Remote Sens. 2018, 10, 2020 8 of 23

Figure 3 .
Figure 3. Segmentation effects and segmentation quality evaluation with different scales.Figure 3. Segmentation effects and segmentation quality evaluation with different scales.

Figure 3 .
Figure 3. Segmentation effects and segmentation quality evaluation with different scales.Figure 3. Segmentation effects and segmentation quality evaluation with different scales.

Spartina
exhibited a dramatic increase from 1990 to 2017, implying an average rate of gain of 0.29 km 2 /y.During 1997-2005, Spartina expanded markedly, in contrast to the mangrove forest change trend during this period.The distribution range of Spartina expanded from sporadic patches on the West Luoyang River in 1990 to large-scale continuous stands along both sides of the Luoyang River in 2017.During the period 1990-2017, aquaculture ponds more than tripled in area (from 1.54 to 5.14 km 2

Figure 6 .
Figure 6.Comparisons of land cover dynamics.(a) Sankey diagram for comparison of total land cover dynamics from 1990 to 2017; (b) spatial distribution of conversion between mangrove forests and other land cover types in four time intervals from the years 1990, 1997, 2005, 2010, and 2017.

Figure 7 .
Figure 7. Spatial distribution map of centroid migration for mangrove forests and Spartina from 1990 to 2017.

Figure 7 .
Figure 7. Spatial distribution map of centroid migration for mangrove forests and Spartina from 1990 to 2017.

Figure 8 .
Figure 8. Positive and negative effects for mangrove forests.(a) Spartina invading the mangrove seedling area; (b) Spartina expansion into mudflats; (c) aquaculture ponds as a threat to mangrove forests; (d) aquaculture development in mudflats; (e) mangrove seedling milestone recorded during the mangrove afforestation project carried out initially in 2001; (f) mangrove seedling cultivation base; (g) artificial control of Spartina by mowing; (h) Spartina were cut off.

Figure 8 .
Figure 8. Positive and negative effects for mangrove forests.(a) Spartina invading the mangrove seedling area; (b) Spartina expansion into mudflats; (c) aquaculture ponds as a threat to mangrove forests; (d) aquaculture development in mudflats; (e) mangrove seedling milestone recorded during the mangrove afforestation project carried out initially in 2001; (f) mangrove seedling cultivation base; (g) artificial control of Spartina by mowing; (h) Spartina were cut off.

Figure 9 .
Figure 9. Artificial seawalls with landward mangrove forests.(a) Artificial seawall located in the western boundary of the study area; (b) artificial seawall located in the eastern boundary of the study area.

Figure 9 .
Figure 9. Artificial seawalls with landward mangrove forests.(a) Artificial seawall located in the western boundary of the study area; (b) artificial seawall located in the eastern boundary of the study area.As shown in Figure10, the annual average temperature and annual precipitation both increased in the QBEWNR from 1990 to 2017.Although a warmer and wetter climate would favor the establishment of mangrove forests [63], such climate change would promote sea-level rise, leading directly to the space loss of mangrove forest expansion seaward.Therefore, it is evident that, in the context of global warming, the potential for major coastal change aggravates the survival risk and pressures of seaward succession for mangrove forests.Overall, as for the study area, climate change and sea-level rise would potentially have adverse effects on mangrove forests.

Figure 10 .
Figure 10.Change in annual average temperature and annual precipitation in the Quanzhou Bay Estuary Wetland Nature Reserve (QBEWNR) from 1990 to 2017.(a) Annual average temperature in QBEWNR from 1990 to 2017; (b) annual precipitation in QBEWNR from 1990 to 2017.

Figure 10 .
Figure 10.Change in annual average temperature and annual precipitation in the Quanzhou Bay Estuary Wetland Nature Reserve (QBEWNR) from 1990 to 2017.(a) Annual average temperature in QBEWNR from 1990 to 2017; (b) annual precipitation in QBEWNR from 1990 to 2017.

Table 1 .
Tidal information of selected satellite images.

Tidal Level Tidal Level of High Tide * Tidal Level of Low Tide *
The study area has a regular semi-diurnal tide.The high and low tide in the table represent the high and low tide, which is the smallest interval with the image transit time, respectively. *

Table 2 .
Description of the landscape classification system and image feature.

Table 2 .
Description of the landscape classification system and image feature.
MangroveAreas covered by mangrove forests With dark red or red, obvious boundary, irregular shape, smooth texture Spartina Salt marshes covered by Spartina With red or light red, fan-shaped or dot shape, smooth texture Mudflat Muddy beaches in the intertidal zone With grey or dark grey, irregular shape, fine and uniform texture Water body Land covered by rivers and shallow sea areas With blue or dark blue, obvious geometric shape, fine and uniform texture Aquaculture pond Man-made farming of aquatic plants and animals in enclosures With dark blue, blue or light blue, obvious configuration, small rectangle shape, smooth texture With dark red or red, obvious boundary, irregular shape, smooth texture Spartina Salt marshes covered by Spartina Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 22

Table 2 .
Description of the landscape classification system and image feature.
MangroveAreas covered by mangrove forests With dark red or red, obvious boundary, irregular shape, smooth texture Spartina Salt marshes covered by Spartina With red or light red, fan-shaped or dot shape, smooth texture Mudflat Muddy beaches in the intertidal zone With grey or dark grey, irregular shape, fine and uniform texture Water body Land covered by rivers and shallow sea areas With blue or dark blue, obvious geometric shape, fine and uniform texture Aquaculture pond Man-made farming of aquatic plants and animals in enclosures With dark blue, blue or light blue, obvious configuration, small rectangle shape, smooth texture With red or light red, fan-shaped or dot shape, smooth texture Mudflat Muddy beaches in the intertidal zone Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 22

Table 2 .
Description of the landscape classification system and image feature.
MangroveAreas covered by mangrove forests With dark red or red, obvious boundary, irregular shape, smooth texture Spartina Salt marshes covered by Spartina With red or light red, fan-shaped or dot shape, smooth texture Mudflat Muddy beaches in the intertidal zone With grey or dark grey, irregular shape, fine and uniform texture Water body Land covered by rivers and shallow sea areas With blue or dark blue, obvious geometric shape, fine and uniform texture Aquaculture pond Man-made farming of aquatic plants and animals in enclosures With dark blue, blue or light blue, obvious configuration, small rectangle shape, smooth texture With grey or dark grey, irregular shape, fine and uniform texture Water body Land covered by rivers and shallow sea areas Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 22

Table 2 .
Description of the landscape classification system and image feature.
Mangrove Areas covered by mangrove forests With dark red or red, obvious boundary, irregular shape, smooth texture Spartina Salt marshes covered by Spartina With red or light red, fan-shaped or dot shape, smooth texture Mudflat Muddy beaches in the intertidal zone With grey or dark grey, irregular shape, fine and uniform texture Water body Land covered by rivers and shallow sea areas With blue or dark blue, obvious geometric shape, fine and uniform texture Aquaculture pond Man-made farming of aquatic plants and animals in enclosures With dark blue, blue or light blue, obvious configuration, small rectangle shape, smooth texture With blue or dark blue, obvious geometric shape, fine and uniform texture Aquaculture pond Man-made farming of aquatic plants and animals in enclosures Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 22

Table 2 .
Description of the landscape classification system and image feature.
Land Cover Type Description Image Feature Mangrove Areas covered by mangrove forests With dark red or red, obvious boundary, irregular shape, smooth texture Spartina Salt marshes covered by Spartina With red or light red, fan-shaped or dot shape, smooth texture Mudflat Muddy beaches in the intertidal zone With grey or dark grey, irregular shape, fine and uniform texture Water body Land covered by rivers and shallow sea areas With blue or dark blue, obvious geometric shape, fine and uniform texture Aquaculture pond Man-made farming of aquatic plants and animals in enclosures With dark blue, blue or light blue, obvious configuration, small rectangle shape, smooth texture With dark blue, blue or light blue, obvious configuration, small rectangle shape, smooth texture
Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 22 Built-up area Lands used for urban and rural settlements, factories or transportation facilities With obvious geometric configuration, cyan or grey, coarse structure Woodland Woody plants grew in terrene greater than 30% With scarlet or dark red, irregular shape, fine and smooth texture Cropland Cultivated land for crops, including paddy field and dry land With dark red or grayish yellow, clear boundary and larger rectangle shape, smooth texture Grassland Natural areas with herbaceous vegetation greater than 30% With red or brown, unclear boundary, irregular shape, smooth structure Barren land Sandy land and areas with less than 5% vegetation cover With bright white, yellowish white or white brown; irregular shape, uniform texture With obvious geometric configuration, cyan or grey, coarse structure Woodland Woody plants grew in terrene greater than 30% Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 22 Built-up area Lands used for urban and rural settlements, factories or transportation facilities With obvious geometric configuration, cyan or grey, coarse structure Woodland Woody plants grew in terrene greater than 30% With scarlet or dark red, irregular shape, fine and smooth texture Cropland Cultivated land for crops, including paddy field and dry land With dark red or grayish yellow, clear boundary and larger rectangle shape, smooth texture Grassland Natural areas with herbaceous vegetation greater than 30% With red or brown, unclear boundary, irregular shape, smooth structure Barren land Sandy land and areas with less than 5% vegetation cover With bright white, yellowish white or white brown; irregular shape, uniform texture With scarlet or dark red, irregular shape, fine and smooth texture Cropland Cultivated land for crops, including paddy field and dry land Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 22 Built-up area Lands used for urban and rural settlements, factories or transportation facilities With obvious geometric configuration, cyan or grey, coarse structure Woodland Woody plants grew in terrene greater than 30% With scarlet or dark red, irregular shape, fine and smooth texture Cropland Cultivated land for crops, including paddy field and dry land With dark red or grayish yellow, clear boundary and larger rectangle shape, smooth texture Grassland Natural areas with herbaceous vegetation greater than 30% With red or brown, unclear boundary, irregular shape, smooth structure Barren land Sandy land and areas with less than 5% vegetation cover With bright white, yellowish white or white brown; irregular shape, uniform texture With dark red or grayish yellow, clear boundary and larger rectangle shape, smooth texture Grassland Natural areas with herbaceous vegetation greater than 30% Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 22 Built-up area Lands used for urban and rural settlements, factories or transportation facilities With obvious geometric configuration, cyan or grey, coarse structure Woodland Woody plants grew in terrene greater than 30% With scarlet or dark red, irregular shape, fine and smooth texture Cropland Cultivated land for crops, including paddy field and dry land With dark red or grayish yellow, clear boundary and larger rectangle shape, smooth texture Grassland Natural areas with herbaceous vegetation greater than 30% With red or brown, unclear boundary, irregular shape, smooth structure Barren land Sandy land and areas with less than 5% vegetation cover With bright white, yellowish white or white brown; irregular shape, uniform texture With red or brown, unclear boundary, irregular shape, smooth structure Barren land Sandy land and areas with less than 5% vegetation cover Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 22 Built-up area Lands used for urban and rural settlements, factories or transportation facilities With obvious geometric configuration, cyan or grey, coarse structure Woodland Woody plants grew in terrene greater than 30% With scarlet or dark red, irregular shape, fine and smooth texture Cropland Cultivated land for crops, including paddy field and dry land With dark red or grayish yellow, clear boundary and larger rectangle shape, smooth texture Grassland Natural areas with herbaceous vegetation greater than 30% With red or brown, unclear boundary, irregular shape, smooth structure Barren land Sandy land and areas with less than 5% vegetation cover With bright white, yellowish white or white brown; irregular shape, uniform texture

Table 3 .
Summary of land cover classification accuracies from 1990 to 2017.
Pro denotes producer accuracy; Use denotes user accuracy.

Table 3 .
Summary of land cover classification accuracies from 1990 to 2017.
Pro denotes producer accuracy; Use denotes user accuracy.

Table 4 .
Change area and annual land change rate (ALCR) of land cover types in study area.

Table 5 .
Conversion comparison between mangrove forests and other land cover types in the Quanzhou Bay Estuary Wetland Nature Reserve (QBEWNR).