Joint Effect of Spartina alterniflora Invasion and Reclamation on the Spatial and Temporal Dynamics of Tidal Flats in Yangtze River Estuary

Tidal flats, which are non-vegetated land–sea transition areas, have an important ecological function in the global ecosystem. However, they have been shrinking in recent years due to natural and anthropogenic activities. Although many studies focus on tidal flats in the Yangtze River estuary (YRE) in China, how reclamation and plant invasion affect the expansion and erosion of tidal flats are still unclear. In this study, we analyzed all of the available Landsat TM/ETM+ /OLI imagery from the period 1996 to 2018 using the Google Earth Engine (GEE) cloud computing platform to obtain annual maps of coastal tidal flats of YRE at 30 m spatial resolution. We chose three sample tidal flats where severe Spartina alterniflora (S. alterniflora) invasion, reclamation, and control areas existed to explore the joint impact of plant invasion and reclamation on tidal flats. We also point out the main driving factor of tidal flat expansion of each island in YRE by multiple linear regression. Our results suggest that the tidal flats of YRE had obvious expansion from 1996 to 2018, and the speed of expansion is getting slower because of the decreasing deposits in the Yangtze River. Invasive S. alterniflora is effective at promoting silting, and tidal flats with S. alterniflora invasion expanded 2.54 times faster than the control group. Chongming and Hengsha Islands were mainly affected by sediment concentration, while Changxing and Jiuduansha Islands were affected by reclamation and S. alterniflora invasion, respectively. The results could be used to support coastal zone management and biodiversity conservation of the YRE.


Introduction
Coastal tidal flats are land-sea transition zones neighboring coastal vegetation areas, and are usually referred to as non-vegetated coastal areas. Coastal tidal flats often serve as important habitats for zoo-benthos (crabs, nematodes, mollusks, etc.), fish, waterfowl, and migrating birds [1][2][3][4][5]. Moreover, they help prevent coastal erosion and protect terrestrial ecosystems from storms and other natural disasters, due to their location as a buffer zone [6]. Although tidal flats have crucial ecological functions, around 16% of global tidal flats were lost between 1984 and 2016, and 50-80% were lost in the Yellow River Delta between the 1950s and the 2000s, arousing global concern [7,8].  [36] 1978,1990,2000,2008 China Niu et al. (2012) [37] 1985,1990,1995,2000, 2005,2010,2014 Yangtze River estuary Chen et al. (2016) [39] 2010 Global Gong et al. (2013) [ The Yangtze River Estuary (YRE), with severe S. alterniflora invasion, reclamation, and control areas existing together, was selected as our study area to compare the impact of plant invasion and reclamation on tidal flat development through analyzing time-series Landsat TM/ETM+ /OLI images from 1996 to 2018 and using the GEE cloud computing platform. The main objectives of this study include the following: 1.
To generate tidal flat mapping in YRE during the period 1996-2018; 2.
To analyze the spatial and temporal dynamics of tidal flats in YRE; 3.
To explore the influence of reclamation and plant invasion on the development of tidal flats in YRE.

Study Area
The study area is located at the mouth of the Yangtze River in eastern China, between 31 • 05 -31 • 40 N and 121 • 40 -122 • 10 E, including Chongming Island (CMI), Changxing Island (CXI), Hengsha Island (HSI), and Jiuduansha Island (JDSI). The development of CXI and HSI is closely associated with reclamation, and their shorelines are smooth and straight. Eastern CMI and JDSI have natural developing tidal flats and rough shorelines. Therefore, our study focused on eastern CMI and middle and lower JDSI with natural shorelines (Figure 1). There is a simple plant zonation pattern in YRE. The three dominant species are P. australis, S. alterniflora, and S. mariqueter, which have zonal distribution. S. mariqueter occupies the low tide zone and P. australis occupies the high tide zone [46]. S. alterniflora competes with P. australis and S. mariqueter in the middle tide zone. P. australis and S. alterniflora usually grow to 1.5-2.0 m tall, while the height of mature S. mariqueter is only 0.3-0.7 m [47]. S. alterniflora was introduced to JDSI in 1997 and Dongtan in 2001 [3,17], so our study period started at 1996, one year before the introduction of S. alterniflora.
Three sample bands 2 km wide were selected to compare the effect of plant invasion and reclamation on tidal flats: Dongtan (DT) on eastern CMI, with both plant invasion with disturbance and reclamation; Jiuduansha (JDS), with only plant invasion and no disturbance; and the control group (CG) on southeastern CMI, with neither plant invasion nor reclamation ( Figure 2). The control group and two sample areas are at the same estuary, which reduces the background interference and allows a comparison of comprehensive disturbances at the same time. At the point of 31 •  Three sample bands 2 km wide were selected to compare the effect of plant invasion and reclamation on tidal flats: Dongtan (DT) on eastern CMI, with both plant invasion with disturbance and reclamation; Jiuduansha (JDS), with only plant invasion and no disturbance; and the control group (CG) on southeastern CMI, with neither plant invasion nor reclamation ( Figure 2). The control group and two sample areas are at the same estuary, which reduces the background interference and allows a comparison of comprehensive disturbances at the same time. At the point of 31°32'13"N, 121°52'53"E, two lines were radiated with a depression of 15° and 48° to traverse the Dongtan Nature Reserve and farm on CMI. The DT line was parallel with the outline of dike 2013, and the CG line was through the center of the farm. The other point was the center of the upper JDSI wetland, at 31°11'38"N, 121°58'48"E, then the same operation was conducted with a depression of 12°, in the same direction as the JDSI tidal flat expansion. For each central line of the three samples, 1 km rectangular buffer samples were generated on both sides, leading to a 2 km wide rectangle with a length varying according to the spatial and temporal development of each sample.

Landsat Data and Pre-processing
Compared with the conventional approach using only a single image per year, all available Landsat 5/7/8 images from 1996 to 2018 were included in this research on the GEE platform, where all image pre-processing tasks could be carried out efficiently using programming and imageprocessing techniques [45]. We counted the number of good-quality observations of Landsat images (Figure 3a), the total number of observations of individual pixels by Landsat 5/7/8 (Figure 3b), and the distribution of annual good-quality observations of all pixels from 1996 to 2018 (Figure 3c). We used a total of 1763 images to generate water frequency maps after removing poor-quality observations, including clouds, shadows, and scan-line corrector (SLC) off gaps by quality assurance

Landsat Data and Pre-Processing
Compared with the conventional approach using only a single image per year, all available Landsat 5/7/8 images from 1996 to 2018 were included in this research on the GEE platform, Remote Sens. 2020, 12, 1725 5 of 20 where all image pre-processing tasks could be carried out efficiently using programming and image-processing techniques [45]. We counted the number of good-quality observations of Landsat images (Figure 3a), the total number of observations of individual pixels by Landsat 5/7/8 (Figure 3b), and the distribution of annual good-quality observations of all pixels from 1996 to 2018 (Figure 3c). We used a total of 1763 images to generate water frequency maps after removing poor-quality observations, including clouds, shadows, and scan-line corrector (SLC) off gaps by quality assurance (QA) band.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 23  Nominalized difference vegetation index (NDVI) [48], enhanced vegetation index (EVI) [49,50], modified normalized difference water index (mNDWI) [51], and land surface water index (LSWI) [52,53] were calculated from the surface reflectance data. NDVI and EVI indicate the greenness of vegetation, while mNDWI is related to open surface water bodies. These four indices were combined to reduce the disturbances of vegetation on water mask and water on vegetation mask [54][55][56].

Open Surface Water Body and Vegetation Annual Frequency Mapping
We used the mNDWI and VIs algorithm to detect water bodies with mixed pixels of vegetation and water body [54]. If the water signal is stronger than the vegetation signal (mNDWI > EVI or mNDWI > NDVI) with green vegetation excluded (EVI < 0.1), we considered this pixel as a water pixel. Thus, the final criterion of the surface water detection algorithm is ((mNDWI > EVI or mNDWI > NDVI) and (EVI < 0.1)) [54,57].
We identified water frequency as the proportion of water observations in all good-quality observations, ranging from 0 to 1. The frequency of water was calculated by Equation (5) below, where SUM water is the number of water observations each year recognized using the above-mentioned mNDWI and VIs algorithm, SUM total is the total observations each year, and SUM bad is the number of bad observations each year. WaterFreq = SUM water SUM total − SUM bad (5) Pixels with water frequency larger than 0.05 of the period were classified as open surface water [54]. The adaptive-thresholding algorithm was proved to map open surface water bodies in many regions successfully [42,54,57], which turned out to be more robust than constant thresholds for further analysis.

Identification of Annual Tidal Flats
Based on our previous study, we considered pixels with vegetation frequency < 0.05 and 0.05 < water frequency < 0.95 as non-vegetation tidal flats [42]. After separating tidal flats by the above-mentioned algorithm, annual maps were filtered to remove all tidal flats inside the dikes, which had an orderly structure with smooth edges. The remaining tidal flats outside the dikes were regarded as our research targets.

Identification of Dike and S. alterniflora
Dikes were identified by manual interpretation according to high-resolution satellite images in Google Earth for each year. We used standard false color images of Landsat TM/ETM+/OLI and phenology signals of dominant species to identify S. alterniflora distribution areas from 1997 to 2018 according to the methods in [47]. The intertidal zone was extracted and classified into seven classes: S. alterniflora, S. mariqueter, P. australis, mudflat, water and diked area. High-spatial-resolution images, historical data in the Resource Monitoring Report from Shanghai Chongming Dongtan National Nature Reserve, and distribution data from previous research were collected to assist with feature recognition and validation. Mudflat mask was covered and large areas of certain plants in the historical data were used to train samples [3,59,60]. S. mariqueter can be easily distinguished from P. australis and S. alterniflora by significant texture, lower zonal distribution, and shorter growing season. P. australis sprouts in early April and grows rapidly in early May. S. alterniflora sprouts in early May and grows rapidly in late June and early July [47]. Thus, green P. australis showing bright red, in contrast to withered S. alterniflora showing dark brown, in standard false color images of early May, can be identified by maximum likelihood classification ( Figure 4) [21]. Interpretation results were converted to vector maps and calculated the areas of S. alterniflora.

Identification of Dike and S. alterniflora
Dikes were identified by manual interpretation according to high-resolution satellite images in Google Earth for each year. We used standard false color images of Landsat TM/ETM+/OLI and phenology signals of dominant species to identify S. alterniflora distribution areas from 1997 to 2018 according to the methods in [47]. The intertidal zone was extracted and classified into seven classes: S. alterniflora, S. mariqueter, P. australis, mudflat, water and diked area. High-spatial-resolution images, historical data in the Resource Monitoring Report from Shanghai Chongming Dongtan National Nature Reserve, and distribution data from previous research were collected to assist with feature recognition and validation. Mudflat mask was covered and large areas of certain plants in the historical data were used to train samples [3,59,60]. S. mariqueter can be easily distinguished from P. australis and S. alterniflora by significant texture, lower zonal distribution, and shorter growing season. P. australis sprouts in early April and grows rapidly in early May. S. alterniflora sprouts in early May and grows rapidly in late June and early July [47]. Thus, green P. australis showing bright red, in contrast to withered S. alterniflora showing dark brown, in standard false color images of early May, can be identified by maximum likelihood classification ( Figure 4) [21]. Interpretation results were converted to vector maps and calculated the areas of S. alterniflora.

Calculating Accumulated and Expanded Tidal Flats Area
The same measures of calculating accumulated and expanded tidal flats area from 1996 to 2018 were adopted in the three sample areas of DT, CG, and JDS ( Figure 5). We used accumulated area to demonstrate the distance that tidal flats have extended toward the sea since 1996. The accumulated tidal flats area of year N was the average value of mosaic-extracted sample tidal flats from 1996 to year N. To separate the expanded tidal flats for each year and avoid the issue of calculating the

Calculating Accumulated and Expanded Tidal Flats Area
The same measures of calculating accumulated and expanded tidal flats area from 1996 to 2018 were adopted in the three sample areas of DT, CG, and JDS ( Figure 5). We used accumulated area to demonstrate the distance that tidal flats have extended toward the sea since 1996. The accumulated tidal flats area of year N was the average value of mosaic-extracted sample tidal flats from 1996 to year N. To separate the expanded tidal flats for each year and avoid the issue of calculating the overlapping Remote Sens. 2020, 12, 1725 8 of 20 area twice, we merged the tidal flats of this year and last year and then extracted those of last year. S in Equations (7) and (8) represents the area of samples.
Remote Sens. 2020, 12, x FOR PEER REVIEW 9 of 23 overlapping area twice, we merged the tidal flats of this year and last year and then extracted those of last year. in Equation (7) and (8) represents the area of samples.

Stastical Analysis of Relationship between Tidal Flats Dynamics and Driving Factors
We took river sediment deposits, S. alterniflora and reclamation area into consideration as three driving factors for tidal flat development. For the time frame of 2000-2018, we conducted multiple linear regressions to quantify and rank the contributions of three features (annual sediment concentration, S. alterniflora area of each year, and reclamation area of each year) among the four islands (CMI, CXI, HSI, and JDSI) and three samples (DT, CG, and JDS) to tidal flat areas. For sample DT, we conducted two models for 2000-2013 and 2000-2018 separately, as the S. alterniflora control project changed the pattern in DT. The S. alterniflora area, reclamation area, and tidal flats area calculated at the same place were correspondingly at the same scale. All subsets regression models were applied to take all independent features into account and inspect all possible multiple linear regression models. The best model was selected according to the adjusted R 2 of the possible models. D-W tests were conducted to guarantee that there was no multi-collinearity among features. The coefficient in the formulas indicates the contribution of each feature.

Accuracy Assessment for Image Classification
We verified the supervised classification results with QuickBird (Date: 2011.09.24) and SuperView-1 (Date: 2018.08.05) high-resolution imagery, assisted by historical field data. A confusion matrix was calculated to compute classification accuracy (Table 2), which presented a high

Stastical Analysis of Relationship between Tidal Flats Dynamics and Driving Factors
We took river sediment deposits, S. alterniflora and reclamation area into consideration as three driving factors for tidal flat development. For the time frame of 2000-2018, we conducted multiple linear regressions to quantify and rank the contributions of three features (annual sediment concentration, S. alterniflora area of each year, and reclamation area of each year) among the four islands (CMI, CXI, HSI, and JDSI) and three samples (DT, CG, and JDS) to tidal flat areas. For sample DT, we conducted two models for 2000-2013 and 2000-2018 separately, as the S. alterniflora control project changed the pattern in DT. The S. alterniflora area, reclamation area, and tidal flats area calculated at the same place were correspondingly at the same scale. All subsets regression models were applied to take all independent features into account and inspect all possible multiple linear regression models. The best model was selected according to the adjusted R 2 of the possible models. D-W tests were conducted to guarantee that there was no multi-collinearity among features. The coefficient in the formulas indicates the contribution of each feature.

Accuracy Assessment for Image Classification
We verified the supervised classification results with QuickBird (Date: 2011.09.24) and SuperView-1 (Date: 2018.08.05) high-resolution imagery, assisted by historical field data. A confusion matrix was calculated to compute classification accuracy (Table 2), which presented a high consistency between ground reference maps and classification maps. The overall accuracy of the classification result was 94.60%, and the kappa coefficient was 0.9293.

Spatiotemporal Dynamics of Tidal Flats Area in YRE
The tidal flats of YRE showed an obvious trend of expansion from 1996 to 2018, following the flow direction of the Yangtze River towards the East China Sea (Figure 6). The northern CMI expanded toward the north, which was nearly in union with southern Jiangsu Province, while the eastern CMI  The total tidal flat areas of CMI and JDSI had large inter-annual variations from 1996 to 2018 (Figure 7, Table 3), indicating natural growing tidal flats of the two islands, and they both had increasing island area from 1996 to 2018. Contrary to the fluctuant tidal flats in CMI and JDSI, CXI had a steady tidal flat area. The sudden increase in island area in 2009 was due to the enclosed construction of the Qingcaosha reservoir, which started in 2007 and was completed in 2010, after which the island area kept steady. The tidal flats area in CXI was also steady before 2006; the consequent increase was due to reclamation along the eastern shoreline. For the same reason, the island area of HSI was steady before 2006, and kept increasing until the reclamation completed in 2017. After the reclamation, HSI was larger than CXI. The total tidal flat areas of CMI and JDSI had large inter-annual variations from 1996 to 2018 (Figure 7, Table 3), indicating natural growing tidal flats of the two islands, and they both had increasing island area from 1996 to 2018. Contrary to the fluctuant tidal flats in CMI and JDSI, CXI had a steady tidal flat area. The sudden increase in island area in 2009 was due to the enclosed construction of the Qingcaosha reservoir, which started in 2007 and was completed in 2010, after which the island area kept steady. The tidal flats area in CXI was also steady before 2006; the consequent increase was due to reclamation along the eastern shoreline. For the same reason, the island area of HSI was steady before 2006, and kept increasing until the reclamation completed in 2017. After the reclamation, HSI was larger than CXI.   1 Uncertainties within the 95% confidence interval. 2 Year of maximum and minimum area in brackets.
Regarding the accumulated tidal flats areas of the three sample areas, DT had the largest accumulated tidal flats area in 2018, which was about 15.66 km 2 . JDS had a relatively larger accumulated tidal flats area, covering about 13.74 km 2 , while CG had the smallest, at 10.71 km 2 . The accumulated area of DT increased from 3.90 km 2 in 1996 to 13.78 km 2 in 2006, and gradually increased afterward, then a sudden rapid increase of 0.91 km 2 was found in 2013. The accumulated tidal flats area of CG increased 5.77 km 2 before 2002, followed by a steady slight increase of 0.65 km 2 for 7 years, then was maintained at 10.7 km 2 since 2009. The accumulated tidal flats area of JDS increased over all the years, from 1.07 km 2 in 1996 to 13.74 km 2 in 2018, as the rate of increase became faster (Figure 8). The accumulated tidal flats area in DT was 0.36 km 2 smaller than CG in 1996. However, the rate of increase of DT was much larger than that of CG in the following years. The accumulated tidal flats area of DT exceeded that of CG in 2000, and JDS exceed CG in 2016, but was always smaller than DT. The decline in 1998, followed by subsequent increases in DT (red up arrow) and CG (blue arrow), indicated the dike constructed in 1998 in CMI, which was across our sample area DT and CG. The sudden increase in 2013 in DT (red down arrow) indicated the S. alterniflora control project across DT that year. accumulated tidal flats area in 2018, which was about 15.66 km 2 . JDS had a relatively larger accumulated tidal flats area, covering about 13.74 km 2 , while CG had the smallest, at 10.71 km 2 . The accumulated area of DT increased from 3.90 km 2 in 1996 to 13.78 km 2 in 2006, and gradually increased afterward, then a sudden rapid increase of 0.91 km 2 was found in 2013. The accumulated tidal flats area of CG increased 5.77 km 2 before 2002, followed by a steady slight increase of 0.65 km 2 for 7 years, then was maintained at 10.7 km 2 since 2009. The accumulated tidal flats area of JDS increased over all the years, from 1.07 km 2 in 1996 to 13.74 km 2 in 2018, as the rate of increase became faster ( Figure  8). The accumulated tidal flats area in DT was 0.36 km 2 smaller than CG in 1996. However, the rate of increase of DT was much larger than that of CG in the following years. The accumulated tidal flats area of DT exceeded that of CG in 2000, and JDS exceed CG in 2016, but was always smaller than DT. The decline in 1998, followed by subsequent increases in DT (red up arrow) and CG (blue arrow), indicated the dike constructed in 1998 in CMI, which was across our sample area DT and CG. The sudden increase in 2013 in DT (red down arrow) indicated the S. alterniflora control project across DT that year.

Major Driving Factors of Tidal Flats Expansion Patterns
Chongming Dongtan experienced rapid growth due to a large quantity of deposits by the Yangtze River over the last century. However, the deposits decreased since the beginning of the 21st century. According to data from the Yangtze River Datong Hydrological Station, the annual sediment

Major Driving Factors of Tidal Flats Expansion Patterns
Chongming Dongtan experienced rapid growth due to a large quantity of deposits by the Yangtze River over the last century. However, the deposits decreased since the beginning of the 21st century. According to data from the Yangtze River Datong Hydrological Station, the annual sediment load decreased from 3.39 × 10 8 (Figure 4).
To quantify the contribution of river sediment, S. alterniflora invasion, and reclamation in each island, a multiple linear regression was conducted (Table 4). Among the four islands, we found that river sediment concentration was the factor contributing the most in CMI and HSI. CXI was largely affected by reclamation, while JDSI was totally promoted by S. alterniflora. The most influential factor in sample CG was sediment concentration. Sample JDS had the same result as JDSI. For sample DT, S. alterniflora was the factor contributing the most before 2013, but, for the whole period, river sediment was still the major factor as CMI. Table 4. Statistical summary of the multiple linear regression formula in four islands and two samples. S. alterniflora is mainly distributed in the northern and northeast CMI. Its area increased from 530.12 ha in 2001 to 2067.46 ha in 2010, and then decreased to 1714.52 ha in 2013. After ecological engineering covered most of the S. alterniflora population, there were still 718.87 ha of patches remaining outside the dike. The area of S. alterniflora in JDSI increased by 214% in the past 15 years. S. alterniflora in JDSI experienced less disturbance and demonstrated clear zonal distribution of natural succession, of which P. australis is in the highest intertidal zone, S. alterniflora is in the middle and lower zones, and S. mariqueter is in the lowest intertidal zone (Figure 4).

Islands and Samples
To quantify the contribution of river sediment, S. alterniflora invasion, and reclamation in each island, a multiple linear regression was conducted (Table 4). Among the four islands, we found that river sediment concentration was the factor contributing the most in CMI and HSI. CXI was largely affected by reclamation, while JDSI was totally promoted by S. alterniflora. The most influential factor in sample CG was sediment concentration. Sample JDS had the same result as JDSI. For sample DT, S. alterniflora was the factor contributing the most before 2013, but, for the whole period, river sediment was still the major factor as CMI. Table 4. Statistical summary of the multiple linear regression formula in four islands and two samples.

Reliability and Uncertainty of Tidal Flats Mapping
In this study, we demonstrated the reliability of generating annual tidal flats mapping at 30 m spatial resolution in YRE. The user's accuracy, producer's accuracy and overall accuracy of classification results were all higher than 90% (Table 2). We used an average of 75 Landsat images each year in the study area (Figure 3b). More than 96.6% of the pixels had four or more good-quality observations, and each pixel had an average of 431 good-quality observations from 1996 to 2018 (Figure 3c). The study was conducted successfully with plenty of good-quality data and accurate classification.
Many global-scale land cover maps including a wetland layer were already produced [29][30][31][32]34], and several researchers have mapped coastal wetlands in China [42,62,63]. In 2019, Murray [8] mapped global tidal flats during 1984-2016 using a random forest algorithm within 3-year windows on the GEE platform. Our results for the trend of tidal flats area in YRE from 1996 to 2016 are consistent with previous results (Figure 11a-d). However, Murray's values were obviously larger due to the inclusion of fish ponds and water channels (outline of reclamation projects) as tidal flats (Figure 11e). The water frequency of the orange area in JDSI (Figure 11g) was 0.99-1, classified as tidal flats in Murray's study, would be classified as permanent water in this study. Considering the global product, a detailed revision for more representative samples should be conducted in YRE.
Although we use adequate Landsat images, annual good-quality observations of pixels showed great spatial and inter-annual variation, and 66% of the pixels had fewer than 15 good-quality observations annually (Figure 3). The variation in Landsat images may lead to uncertainty in annual tidal flats mapping [38], but it is still useful to track inter-annual dynamics.
Besides image data and algorithms, the background environment of three sample areas in the past 30 years were not exactly the same, especially the coastal surface water salinity, which may result in the different growth rate of wetland vegetation and affect tidal flat development. According to the data of YRE in the comprehensive investigation of the coastal zone and tidal resources and the project of sediment resources and utilization in 2012, the surface salinity of both CG and JDS was around 10 PSU in winter [64]. However, after construction of the YRE deep water channel, a higher level of saltwater encroachment occurred in the south passage (south of JDSI) and north branch (between CMI and Jiangsu Province), leading to higher surface salinity of JDS and DT. The surface salinity of CG was around 2PSU lower than JDS during the neap tide. The salinity differences of wetland cause divergent S. alterniflora competitive capacity, as it grows better in higher salinity wetlands [3], which is also why it leads to faster sediment accumulation in JDS and DT.

Impact of S. alterniflora Invasion on Tidal Flats Expansion Patterns
With a great capacity for decreasing tidal velocity, mitigating erosion, and trapping sediments, S. alterniflora has been introduced to many coastal regions in the world to protect dikes and promote silting for natural reclamation [65], which turned out to be successful to some extent. The plant did reduce the influence of typhoons and protect the coastal environment of Zhejiang against typhoons in 1990 and 1994 [66]. Also, S. alterniflora is responsible for the rapid sediment accretion in CMI and JDSI besides anthropogenic reclamation. Even though deposits from the Yangtze River are decreasing, YRE is still expanding overall. S. alterniflora has larger stem and rhizome density compared with S. mariqueter and P. australis, so it intercepts more suspended particles with smaller sizes, which helps make it an ecosystem engineer [3,15]. Enhanced net primary production by S. alterniflora strengthened the carbon sink of the wetland [67]. In addition, S. alterniflora provided compatible habitats for native crab S.dehaani in YRE [68], and crab increased the competition ability of S. alterniflora by increasing soil nitrogen [69]. However, S. alterniflora replaced S. mariqueter in YRE and occupied mangrove zones in Southern China [22,23], reducing the food resources and habitat of local wetland animals. Invasive S. alterniflora homogenized and lowered the species diversity of the

Impact of S. alterniflora Invasion on Tidal Flats Expansion Patterns
With a great capacity for decreasing tidal velocity, mitigating erosion, and trapping sediments, S. alterniflora has been introduced to many coastal regions in the world to protect dikes and promote silting for natural reclamation [65], which turned out to be successful to some extent. The plant did reduce the influence of typhoons and protect the coastal environment of Zhejiang against typhoons in 1990 and 1994 [66]. Also, S. alterniflora is responsible for the rapid sediment accretion in CMI and JDSI besides anthropogenic reclamation. Even though deposits from the Yangtze River are decreasing, YRE is still expanding overall. S. alterniflora has larger stem and rhizome density compared with S. mariqueter and P. australis, so it intercepts more suspended particles with smaller sizes, which helps make it an ecosystem engineer [3,15]. Enhanced net primary production by S. alterniflora strengthened the carbon sink of the wetland [67]. In addition, S. alterniflora provided compatible habitats for native crab S.dehaani in YRE [68], and crab increased the competition ability of S. alterniflora by increasing soil nitrogen [69]. However, S. alterniflora replaced S. mariqueter in YRE and occupied mangrove zones in Southern China [22,23], reducing the food resources and habitat of local wetland animals. Invasive S. alterniflora homogenized and lowered the species diversity of the nematode community and altered the diets and community structure of local arthropods and macro-benthonic invertebrates [3,[70][71][72][73], leading to further impact on migratory birds [19]. According to the growth features of S. alterniflora and clustered plantation by human introduction, small patches first formed on the low-tide areas of tidal flats, and then merged and formed larger patches sprawling toward higher tidal flats and competing with native vegetation. Thus, tidal flats with S. mariqueter at the front experienced much slower S. alterniflora expansion than tidal flats [74]. Moreover, human disturbances to native vegetation left gaps for S. alterniflora to colonize and suppress native vegetation recruitment [22], which provided the environmental conditions for S. alterniflora expansion [75]. On the other hand, larger patches of S. alterniflora filled in the tidal creeks toward the sea and occupied the area of tidal flats, which resulted in the smaller expansion area of tidal flats in one year and larger expansion area in previous years.

Impact of Ecological Engineering on Tidal Flats Expansion Patterns
A dike with elevation above sea level can control the ecological consequences of previous S. alterniflora invasions [25], and meet the needs of the land in cities and keep the expansion of tidal flats. The construction of dikes on Dongtan promoted the outward expansion of tidal flats, in that the dikes intercepted the route of the flood tide, leading to decreased tidal velocity and causing a large amount of sediment deposit on the tidal flats outside the dike [76]. The promotion effect is obvious in the first year and is weakened over time until the sediment flux of upper and lower tidal flats reaches a new balance. During the construction of dikes, the sediment environment of neighboring tidal flats was disturbed, the supply of which reduced tremendously in a short time due to large-scale sand borrowing, leading to erosion that year. However, the sediment concentration after construction was higher than pre-construction [77], which resulted in the expansion of tidal flats in the following years. However, the dike constructed below sea level in northern CMI block the tide from entering tidal flats inside the dike, leading to a tremendous reduction in and even disappearance of Polychaeta species and increased insect larvae and mollusks species [78]. The disappearance of important benthic animals because of habitat loss is direct evidence of excessive reclamation [79].

Joint Effect of S. alterniflora Invasion and Reclamation on Tidal Flats
Land reclamation played an important role on CMI. Anthropogenic reclamation was dispersed and small scale in [1985][1986][1987][1988][1989][1990]. Two large dikes, dike 1992 and dike 1998, were constructed in the 1990s [80], leading to a rapid decline in tidal flats at that time. The reclamation dike of the S. alterniflora invasion control project, which replaced S. alterniflora with S. mariqueter, was carried out in 2013. A large area of bare land was identified for tidal flats, which was covered by S. alterniflora during the gap between S. alterniflora cutting and S. mariqueter replantation, leading to a tremendous increase in tidal flats in 2013 and a low growth rate in the following years after replantation [42]. The control project is the vital reason for S. alterniflora shrinkage. Liu et al. [22] provided a similar result, proving that land reclamation on tidal flats or S. alterniflora vegetation to aquaculture ponds and mangrove replantation led to S. alterniflora shrinkage.
The control project was effective in the controlled region, but there were still S. alterniflora patches dispersed outside the dike and quite a large area observed along the northern CMI shoreline [23]. Since S. alterniflora was proved to be robust in promoting silting, it would be time consuming and labor intensive to completely remove them without large construction. It is strongly recommended that the S. alterniflora patches on the eastern part of island be removed, otherwise, the tragedy would be repeated. However, for the land ascertained as a reclamation region for fields, ponds, or building areas in the near future, we suggest retaining S. alterniflora to promote silting in the diked region. The outline dikes were important because S. alterniflora must be under control. After the reclamation project is complete, it should be removed if the region will no longer be wetland, in case of its fast spread in the surrounding area.

Future Tidal Flats in YRE
Besides river sediment deposits, the main factor that affected CXI and HSI tidal flat areas was reclamation, and their shorelines were both solidified. The Yangtze estuary deep-water channel was located along the northern shoreline of the two islands. There is a nearly 20 km deep shoreline with little sediment deposit along southern CXI, which is an ideal place for a port. The largest shipyard in China, Shanghai South of Yangzi River Changxing Shipyard, is located at the southern shoreline, while Qingcaosha reservoir was located at the northern shoreline. There was a large area of reclamation along eastern HSI during 2006-2017, which is nearly the same area as CXI. We can predict that both the tidal flats and island areas of CXI and HSI will be stable in the future, since their shorelines have been smoothed and solidified by man-made shorelines.
S. alterniflora was the main factor that promoted tidal flats in JDSI. One reason to introduce S. alterniflora to JDSI was to attract migratory birds and draw them away from the nearby Pudong International Airport [10]. Shanghai Jiuduansha Wetland National Nature Reserve covers all of JDSI, so no reclamation project could be easily conducted there. Although S. alterniflora has been proved to be harmful to migratory bird habitats, the focus of the JDS nature reserve is to protect the wetlands rather than the bird habitats. To a large extent, S. alterniflora would be kept in the future. The situation in CMI would be more complicated. Sediment deposits had a much larger contribution than S. alterniflora and reclamation. Contrary to the JDS nature reserve, Chongming Dongtan National Nature Reserve only covers eastern CMI, not the whole island, so the control project discussed above was only conducted there for bird habitat protection. However, there was still S. alterniflora left outside the dike. If the management suggestions we provided above had been accepted, the tidal flats and island areas of CMI and JDSI would still increase. The tidal flat area of JDSI would keep its present speed with the help of S. alterniflora, and eastern CMI would increase rather more slowly with several sudden increases as most S. alterniflora is removed, but more reclamation projects may be conducted in the northern island.

Conclusions
This paper demonstrated the spatial and temporal dynamics of tidal flats in YRE during 1996-2018 using the GEE cloud computing platform, analyzing the effects of invasion by S. alterniflora and reclamation on the expansion and erosion of tidal flats. Three typical sample areas in YRE, CG, JDS, and DT, represent the control group, S. alterniflora invasion, and joint controls, respectively, which reduce the background interference and allow the comparison of comprehensive disturbance at the same time. The results illustrate the joint effect of S. alterniflora invasion and reclamation by comparing the three sample areas and point out the main driving factors of tidal flats in four islands. Due to the decreasing river deposits in the Yangtze River, the expansion speed of tidal flats in YRE is getting slower, but the trend of expansion is still obvious. River sediment concentration had the largest contribution in CMI and HSI, while reclamation and S. alterniflora invasion were the most vital factor in CXI and JDSI, respectively. S. alterniflora has a robust capability of silting and can accelerate tidal flat expansion by 2.54 times. Further monitoring of S. alterniflora on Dongtan is necessary for migratory birds, and retaining the remaining S. alterniflora in northern CMI would be an effective approach to accelerate reclamation if carefully enclosed in a limited region. Therefore, our findings could be used to support coastal zone management and biodiversity conservation of the YRE in the future.