Remote Sensing of Spatiotemporal Changes in Wetland Geomorphology Based on Type 2 Fuzzy Sets : A Case Study of Beidagang Wetland from 1975 to 2015

Few studies have considered the spatiotemporal changes in wetland land cover based on type 2 fuzzy sets using long-term series of remotely sensed data. This paper presents an improved interval type 2 fuzzy c-means (IT2FCM*) approach to analyse the spatial and temporal changes in the geomorphology of the Beidagang wetland in North China from 1975 to 2015 based on long-term Landsat data. Unlike traditional type 1 fuzzy c-means methods, the IT2FCM* algorithm based on interval type-2 fuzzy set has an ability to better handle the spectral uncertainty. Four indexes were adopted to validate the separability of classes with the IT2FCM* algorithm. These four validity indexes showed that IT2FCM* obtained better results than traditional methods. Additionally, the accuracy of the classification results was assessed based on the confusion matrix and kappa coefficient, which were high for the analysis of wetland landscape changes. Based on the analysis of separability of classes with the IT2FCM* algorithm using four validity indexes, the classification results, and the membership value images, the long-term series of satellite datasets were processed using the IT2FCM* method, and the study area was classified into six classes. Because water resources and vegetation are two key wetland components, the water resource dynamics and vegetation dynamics, based on the normalized difference vegetation index (NDVI), were analysed in detail according to the spatiotemporal classification results. The results show that the changes in vegetation types have historically been associated with water resource variations and that water resources play an important role in the evolution of vegetation types.


Introduction
The official federal definition of wetlands was provided by Cowardin et al. [1] as follows: "Wetlands are transitional lands between terrestrial and aquatic systems where the water table is usually at or near the surface or the land is covered by shallow water".Wetlands are considered an immensely important part of the global ecosystem and support high levels of biodiversity, including significant numbers of rare and endangered plant and animal species [2,3], although they only cover approximately 6% of the Earth's land surface [4,5].Wetlands also provide flood protection, protection from storm and wave damage, water quality improvements through the filtering of agricultural and industrial waste, and aquifer recharge [6][7][8].Although wetlands play an important role in environmental processes and provide many ecological services for mankind, they are under increasing threats from climate change, dredge and fill operations, hydrologic variations, pollutant runoff, eutrophication, impoundments, and land modification, as well as the growing load placed on wetlands by agriculture, aquaculture, forestry, and development [9][10][11][12][13].
To better conserve and manage wetland resources, it is important to inventory and monitor wetlands and their adjacent uplands.Photographic evidence of hydrological conditions, in combination with collateral data, is sufficient to accurately document wetland existence.Recently, with the rapid development of image processing techniques, wetland resources monitoring and corresponding land cover changes based on satellite remote sensing has displayed considerable advantages over other traditional methods, especially over large geographic areas.Many clustering algorithms and methods for land cover classification have been proposed based on remotely sensed datasets with different spatial and spectral resolutions, such as NOAA AVHRR, Terra MODIS and ASTER, Landsat TM/ETM+, and EO-1 Hyperion and ALI, as well as airborne remote sensing datasets, including AVIRIS, HyMap, and others [14][15][16][17][18][19][20].Compared with traditional Boolean classification methods, clustering algorithms based on fuzzy sets are more appropriate for land cover mapping [21][22][23].
Fuzzy c-means (FCM) clustering is an important image classification algorithm and is widely used in remote sensing image analysis.In image segmentation, a record in a fuzzy classification reflects the degree to which any object belongs to all candidate classes [24].Thus, the outcome of a fuzzy classification is a record for every object being analysed of the degree to which that object belongs to every class being considered.In contrast, a Boolean classification simply records that an object either is (1) or is not (0) part of the set.Typically, this is done by recording only the set to which each object belongs.Therefore, in land cover classification, for any pixel or object under consideration, a Boolean classification might simply record that the pixel or object belongs to the class grassland among a set of five classes [grassland, woodland, grass land, bare soil, water], while a fuzzy classification might record a list of the degree of membership of each class, e.g., [0.5, 0.2, 0.2, 0.1, 0].It is not required that the membership degrees should sum to 1, but this is a common restriction of fuzzy classification [23].
The original clustering algorithms, such as k-means, fuzzy c-means (FCM), and kernel fuzzy c-means, are all type 1 FCM methods, and most use Euclidean distance to quantify the degree of similarity between data points and the corresponding membership degree [23,25,26].FCM with spatial information has been applied to enhance segmentation algorithms for use with remotely sensed imagery containing noise [27,28].Shao et al. proposed a novel and robust semi-supervised FCM (RSFCM) clustering method to detect more changes than traditional methods and provide noise immunity by synergistic exploitation of pseudo-labels based on a difference image and spatial information [26].Wang et al. proposed an adaptive spatial-information-based fuzzy clustering algorithm for image segmentation that overcomes the sensitivity to noisy information and the lack of spatial information, which improved the robustness of traditional FCM methods [29].
Since type-1 fuzzy set theory was first introduced in [29], conventional type 1 FCM methods such as k-means and FCM have been developed and applied to many domains due to their ability to model fuzziness.However, the classical FCM clustering methods are based on type-1 fuzzy set theory, which cannot address uncertainties associated with membership grade; they do not perform well on data with complicated geometry, and they are unable to handle and quantify uncertainty when determining the associated membership function [30].This fault can be overcome using the type-2 fuzzy set and the type-N fuzzy set introduced by Zadeh (1975) [31].The type-2 fuzzy set is an extension of the type-1 fuzzy set that makes computation extremely difficult, and it is also difficult to handle the secondary membership value or function.The interval type-2 fuzzy set (IT2 FS) is a special case of the type-2 fuzzy set for which the secondary membership grade is a constant equal to one (Mendel et al., 2002) [32]; consequently, it requires less computation than the type-1 fuzzy set and has been adopted in many domains [33].In contrast, type 2 fuzzy sets have unique advantages and can be used to assess the uncertainty that exists in image data due to environmental factors such as weather conditions and sensor issues.Thus, these methods are used in data clustering, such as via the interval type 2 FCM clustering algorithm.A boosted genetic fuzzy classifier (BGFC) was developed by Stavrakoudis et al. for land cover classification using multispectral remote sensing data [34].In this method, the fuzzy rules are generated in an iterative fashion, and a genetic tuning stage is also adopted to improve the cooperation among the fuzzy rules.The classification results of this method are generally good.Ghaffarian developed an automatic FCM algorithm based on a histogram of land cover classification and multispectral data [35].Ghosh performed supervised change detection using fuzzy clustering combined with the spatial correlation between neighbouring pixels in a difference image, which was calculated from two images of the same location at different times [36].In the standard FCM method, two parameters that can be set by users are the class number, C, and the fuzzifier, m [36].The cluster centres are expressed by real number vectors, and the distance between a sample and the cluster centres is used to determine the membership degree of a sample belonging to one class.As mentioned before, classical FCM clustering methods cannot handle the uncertainty of membership degree.An interval type 2 fuzzy c means (IT2FCM) was proposed in [33] based on the IT2 FS to handle fuzziness uncertainty in fuzzy clustering tests.
Until recently, however, very few studies have focused on land cover and land use based on IT2FS methods, and there are almost no studies regarding the use of such methods in spatiotemporal analyses of wetland dynamics.Thus, the first objective of this paper is to develop a study of wetland land cover change analysis based on the interval type-2 fuzzy set.Notably, the spectrum of a geographic feature does not refer to a single spectral curve but rather to a connected set of possible spectral curves, namely, a spectral curve with a certain width.The existing IT2FCM method does not take the width of a band into account; instead, it typically uses the ranks of the average values of lower and upper membership grades to determine whether the considered pixel belongs to a specific class.Moreover, the interval number ranking based on probability is never considered.Hence, the second objective of is to propose an approach that incorporates the interval number distance and ranking methods into IT2FCM.The method is called IT2FCM*, and it reflects the uncertainty in remotely sensed data [37].Additionally, although many researchers have used the clustering validity indexes of the classical FCM algorithm, four extended cluster validity indexes are adopted in this paper: the partition coefficient (PC-), the partition entropy (PE-), the Fukuyama and Sugeno index (FS-), and the Xie and Beni index (XB-).Finally, based on IT2FCM*, the spatial and temporal changes (from 1975 to 2015) in the Beidagang wetland, North China, are analysed based on a long-term series of Landsat satellite data.Three important factors indicating wetland geomorphology changes are introduced to the spatiotemporal analysis: water area, vegetation index, and vegetation types.Water plays a key role in wetlands.Thus, monitoring the changes in water areas is very important for the wetland.Aquatic vegetation in a wetland not only plays an important role in water purification by directly uptaking and utilizing nutrients in wastewater and absorbing heavy metals and other hazardous substances that accumulate in sewage but also enhances the capacity for the absorption and storage of carbon dioxide in the ecosystem through carbon sequestration.Thus, it is also very important for environmental wetland research to monitor the growth conditions of aquatic vegetation and their changes and evolution.

Study Area
This research was performed in the study area of Tianjin Municipality, a metropolis located in northern coastal China (see Figure 1a) with a latitude ranging from 38 • 34 to 40 • 15 N and a longitude ranging from 116 • 43 to 118 • 04 E, covering an area of 11,860.63square km (4579.41square miles).Tianjin is located along the west coast of the Bohai Sea (the West Pacific), and it is the largest coastal city in northern China.There are 153 km (95 miles) of coastline and 1137.48 km (706.80 miles) of land border.Tianjin is divided into 16 county-level divisions, all of which are districts (see Figure 1b) with a typical East Asian monsoon-influenced climate with windy and dry winters and hot and humid summers.The monthly 24-h average temperature ranges from −3.4 • C in January to 26.8 • C in July, with an annual mean of 12.90 • C. Annual precipitation averages only 511 mm, of which nearly three-fifths occurs in July and August, and annual evaporation totals 1120.5 mm.Thus, the city lies within a semi-arid zone [38].
Remote Sens. 2017, 9, 683 4 of 24 with an annual mean of 12.90 °C.Annual precipitation averages only 511 mm, of which nearly three-fifths occurs in July and August, and annual evaporation totals 1120.5 mm.Thus, the city lies within a semi-arid zone [38].The wetland area covered 5471 square km in Tianjin in the 1820s, accounting for 45.9% of the total area of the city.However, the area of natural wetland had decreased to 3.6% in 2000.This represents a 54.7% reduction in natural wetland area compared with that in the 1950s, and the urban wetland area decreased by 80% [39].Even worse, this decreasing trend continues.In this context, the variation in the wetland area in Tianjin District has been studied.Dahuangpu, Beidagang and Qilihai are the main parts of the wetland in Tianjin.Among them, Dahuangpu has been the most The wetland area covered 5471 square km in Tianjin in the 1820s, accounting for 45.9% of the total area of the city.However, the area of natural wetland had decreased to 3.6% in 2000.This represents a 54.7% reduction in natural wetland area compared with that in the 1950s, and the urban wetland area decreased by 80% [39].Even worse, this decreasing trend continues.In this context, the variation in the wetland area in Tianjin District has been studied.Dahuangpu, Beidagang and Qilihai are the main parts of the wetland in Tianjin.Among them, Dahuangpu has been the most severely affected by humans; its ecological value has decreased significantly as a result of the excavation of large areas of fishponds.In this research, the Beidagang wetland is selected as the study area due to its provincial significance regarding local climate conditions, its suitable size and the availability of moderate-resolution Landsat satellite data in the area.This area includes three or more wetland types (see Figure 1c).There are no clear boundaries between marshes and reservoirs, and wet areas frequently vary between each classification in the Beidagang wetland.Thus, both the Beidagang wetland and the Beidagang Reservoir were studied.
The soil of Beidagang is mainly composed of low-lying silt produced due to epeirogenesis (the formation and submergence of continents by broad, relatively slow displacements of the earth's crust) of the seacoast.With many crisscrossing rivers and ponds and abundant groundwater resources, the Beidagang wetland mainly includes five types of land cover types: lakes, beaches, marshes, rivers, and shallow sea areas.These land cover types encompass numerous ecosystem and biological resources.Notably, six types of first-class nationally protected birds and seventeen types of second-class nationally protected birds inhabit these areas.The wading bird population, which includes some endangered birds, can reach nearly eight million each year [39].

Datasets and Preprocessing
In this paper, three types of remotely sensed datasets are used, including the Advanced Land Imager (ALI) onboard the Earth Observation (EO-1) Mission, the multi-spectral scanner onboard the Sentinel-2 and the long-time series of Landsat datasets.The climate datasets for 22 August 2015 and 28 August 2016 in the study area were freely collected through the China Meteorological Administrative official website [40].This climate dataset was used in combination with the Second Simulation of the Satellite Signal in the Solar Spectrum (6S) software to remove the atmospheric effects for the ALI and the Sentinel-2 datasets.In addition, the field survey dataset was collected on 27 June 2016 in the study area and its surrounding areas.The EO-1 ALI, Sentinel-2, and Landsat 8 satellite datasets were used to carry out an experiment to test the separability of the improved IT2FCM* algorithm.The long-time series of Landsat satellite datasets were used to make the spatiotemporal analysis of wetland cover changes.The acquisition and preprocessing of these datasets are depicted in the following section.

EO-1 ALI Dataset
The multi-spectral dataset acquired on 22 August 2015 collected from the ALI instrument onboard the EO-1 satellite was used in this paper.The ALI instrument onboard the EO-1 satellite is a technology demonstration that serves as a prototype for the Landsat Data Continuity Mission.The ALI observes the Earth in 10 spectral bands whose spectral region ranges from 0.43 to 2.35 µm, and the pan band has a spatial resolution of 10 m, with the spectral region ranging from 0.48 to 0.69 µm.The ALI dataset is calibrated from digital number to radiance value.With the help of the 6S radiometric transfer model and a climate dataset, the atmospheric effects were removed to build the reflectance imagery [41,42].

Sentinel-2 Dataset
Sentinel-2A was launched on 23 June 2015 as a polar-orbiting, multispectral high-resolution imaging mission for land monitoring [43].In this paper, a cloud-free multispectral high-resolution image from the Sentinel-2 satellite was acquired on 30 May 2016, consisting of 13 spectral channels in the visible, near infrared, and short wave infrared part of the spectrum.Among these 13 bands, 3 channels are designed to detect the coastal aerosol (0.443 µm), the water vapour (0.945 µm), and the cirrus (1.375 µm), in which the spatial resolution is 60 m.The spatial resolution of the visible bands 2-4 and NIR band 8 (whose spectral region ranges from 0.705 to) is 10 m, and the spatial resolution of bands 5-7 (ranging from 0.705 to 0.783 µm), band 8 (0.865 µm), and bands 11 and 12 (whose wavelengths are 1.610 µm and 2.190 µm, respectively) is 20 m.The satellite dataset used in this study is a Level 1C product.Before classification, the Level 1C multi-spectral images were processed into Level 2A reflectance images.This step was accomplished using the Sentinel 2 Level 2A Processor (Sen2Cor) plugged into the Sentinel Application Platform (SNAP) software provided by the European Space Agency (ESA) official website [44] with the anniversary dataset, water vapor and cloud detecting bands (for more details, please refer to the official website [45]).

Landsat Dataset
Landsat 2 multispectral scanner (MSS), Landsat 5 thematic mapper (TM) and enhanced thematic mapper (ETM+), and Landsat 8 operational land imager (OLI) images were selected to conduct a spatiotemporal analysis of the Beidagang wetland.The spatial resolution of the MSS sensor on board Landsat 2 is approximately 79 m and is often processed to the pixel size of 60 m [46].The TM sensors on board the Landsat 4/5 satellites have a spatial resolution of 30 m, 6 spectral channels ranging from 0.4 µm to 2.25 µm, a sixth band, and a thermal-spectral channel (10.45~12.42µm) with a spatial resolution of 120 m.The Landsat ETM+ sensor has spatial resolutions of 30 m for the six reflective bands, 60 m for the thermal band, and 15 m for the panchromatic band.Compared with Landsat ETM+ images, Landsat 8 OLI images consist of nine spectral bands in the VNIR and SWIR spectral regions.Bands 1 to 7 and 9 have a spatial resolution of 30 m.The resolution of band 8 (panchromatic) is 15 m, and bands 10 and 11 are thermal-spectral channels with a spatial resolution of 100 m, which are used to provide accurate surface temperatures.
In this study, long-term series of cloud-free summer Landsat satellite data from 1975 to 2015 were collected at intervals of 5 years and used to evaluate the dynamics of the wetland and vegetation fractions in the Beidagang wetland.Landsat 8 OLI dataset of 2016 was also collected and used for classification testing with the combination of the abovementioned ALI and Sentinel-2 dataset.All the Landsat satellite datasets were freely obtained from USGS website [47] and are reflectance images that were corrected by removing atmospheric effects (See Table 1).The preprocessing of the Landsat 4-5 and 7 images used the radiative transfer model 6S , while the pressure, water, and air temperature measurements were obtained from the National Centres for Environmental Prediction (NCEP) grid dataset.The preprocessing of the Landsat 8 images used the radiative transfer model internal algorithm, the water vapor and air temperature measurements were obtained from Moderate Resolution Imaging Spectroradiometer (MODIS) Climate Modelling Grid-Aerosol (CMA), and the pressure was calculated internally based on the elevation.For detailed information we refer readers to [48][49][50].

Field Survey Dataset
In this study, the field survey was carried out on 27 June 2016.Numerous photographs were collected using a digital camera (a Cannon EOS 60D, whose resolution is 18.1 megapixels and uses an APS-C CMOS Sensor, produced by the Cannon Company, made in Japan) (see Figure 2, these photos correspond to the location in the red circle shown in Figure 1c).Additionally, ground truth data was collected through a global positioning system (GPS) terminal at approximately 238 points within the study area (approximately 82 polygons were delineated) and another 105 points around the study area to validity the classification results from the ALI dataset (approximately 37 polygons were delineated).Due to the difficulty of accessing the marsh areas, most of the polygons were collected on both sides of roads.These field survey datasets were used to validate the accuracy of the improved IT2FCM* algorithm.The entire flowchart of the study is given in Figure 3.The entire flowchart of the study is given in Figure 3.

Introduction of IT2FCM
The class number C and the fuzzifier m are two parameters that can be set by users [51] in FCM.The cluster centres are expressed by real number vectors, and the distance between a sample and the cluster centres is used to determine the membership grade of a sample belonging to one class.The fuzzy membership grades of a sample of a specific class can differ for different values of the fuzzifier m for a given classification number, C [23,52].As mentioned above, Classical FCM clustering methods cannot handle uncertainty of membership degree.The concept of a type-2 fuzzy set was introduced by [31] as an extension of the concept of an ordinary fuzzy set (henceforth called a type-1 fuzzy set).A type-1 fuzzy set has a crisp membership degree, whereas a type-2 fuzzy set has fuzzy membership degrees-it can be considered as a "fuzzy-fuzzy set".Such sets are particularly useful in circumstances when it is difficult to determine an exact membership function for a fuzzy set; hence, type-2 fuzzy sets are useful for incorporating uncertainties.Hwang and Rhee proposed IT2FCM based on the interval type II fuzzy set (IT2FS) to address uncertainty in the choice of m, which uses lower and upper membership grades or functions to address the uncertainty in the fuzzy membership value [52].The lower and upper membership functions are constructed using two fuzzifiers, m1 and m2, in IT2FCM.Then, two objective functions can be expressed as follows: where =∥ − ∥ is the distance metric between the sample and the cluster centroid , C is the number of clusters, N is the number of samples, U is the matrix of membership values, k is the sample, i is the clusters, and and are the lower and upper membership grades, respectively, which can be expressed as follows:

Introduction of IT2FCM
The class number C and the fuzzifier m are two parameters that can be set by users [51] in FCM.The cluster centres are expressed by real number vectors, and the distance between a sample and the cluster centres is used to determine the membership grade of a sample belonging to one class.The fuzzy membership grades of a sample of a specific class can differ for different values of the fuzzifier m for a given classification number, C [23,52].As mentioned above, Classical FCM clustering methods cannot handle uncertainty of membership degree.The concept of a type-2 fuzzy set was introduced by [31] as an extension of the concept of an ordinary fuzzy set (henceforth called a type-1 fuzzy set).A type-1 fuzzy set has a crisp membership degree, whereas a type-2 fuzzy set has fuzzy membership degrees-it can be considered as a "fuzzy-fuzzy set".Such sets are particularly useful in circumstances when it is difficult to determine an exact membership function for a fuzzy set; hence, type-2 fuzzy sets are useful for incorporating uncertainties.Hwang and Rhee proposed IT2FCM based on the interval type II fuzzy set (IT2FS) to address uncertainty in the choice of m, which uses lower and upper membership grades or functions to address the uncertainty in the fuzzy membership value [52].The lower and upper membership functions are constructed using two fuzzifiers, m 1 and m 2 , in IT2FCM.Then, two objective functions can be expressed as follows: where d ik = x k − v i is the distance metric between the sample x k and the cluster centroid v i , C is the number of clusters, N is the number of samples, U is the matrix of membership values, k is the sample, i is the clusters, and u ik and u ik are the lower and upper membership grades, respectively, which can be expressed as follows: where i = 1, 2, . . ., C and k = 1, 2, . . ., N.
The cluster centroid v i is represented by the interval between v L and v R and computed in the same way as in FCM methods: The KM algorithm [53] is adopted to determine v L and v R , and v i can be obtained via type reduction: The membership grades can be determined as follows: In addition, u L ij and u R ij can be determined as follows: where M is the number of sample features.Then, the class to which a sample belongs can be determined by the following rule: if u ik > u jk for j = 1, 2, . . ., C and i = j, then xk is assigned to cluster i.
Cluster centres are expressed as real vectors in existing FCMs, so errors of these centres cannot be dealt with effectively.The interval centroids of all clusters must be type reduced to singleton values.Moreover, the type 2 fuzzy membership functions are defuzzified into type 1 fuzzy membership functions during each iteration; therefore, some information is lost [54].All the derivative methods of IT2FCM have the same deficiency, including the interval type 2 fuzzy possibilistic c-means (IFPCM) [55], interval-valued possibilistic FCM (IPFCM) [56], general type 2 FCM (GT2 FCM) [57], interval type 2 FCM clustering with spatial information (IIT2-FCM) [24], and the kernel interval-valued FCM (KIFCM) [30] clustering algorithms.
The spectral characteristic curves of almost all geographical features should be bands with a certain range (i.e., not curves) [58,59]; therefore, the same spectrum in a satellite image may correspond to a different classification in the natural or real world [60].The sensors, data-acquisition process, processing steps, and conversion and transmission processes may produce uncertainty in remote sensing data [59,[61][62][63].

IT2FCM* Algorithm
The lower and upper membership functions are constructed using the two fuzzifiers: m 1 and m 2 , but IT2FCM and its extended algorithm such as IIT2-FCM still have some faults which are discussed in the preceding section.Therefore, we proposed the IT2FCM* based on based on the interval number distance and ranking.There are some distance definitions between interval numbers.For example, the Euclidean distance between interval numbers is commonly used, but its obvious default is that this definition considers only the endpoints of the interval numbers.We tested these distance definitions in IT2FCM* and proved that the definition proposed by Li et al. which is an extension of the definition proposed by Bertoluzza et al. is more effective than others [64].Hence, the interval number distance proposed in [64] was adopted in the improved IT2FCM* algorithm.Let a = [a − , a + ] and b = [b − , b + ] be two interval numbers.Then, the interval distance between a and b is calculated as follows: Just as in IT2FCM.The lower and upper membership functions are constructed using two fuzzifiers: m 1 and m 2 .Then, two objective functions can be established by Equation (1).
All the cluster centroids are interval number vectors in IT2FCM*.Thus, the interval vector distance between a sample and an interval cluster centroid can be expressed as follows: where x is a sample, v is an interval number vector, i = 1, 2, . . .,M, and M is the number of features.Each v i in v has a lower and upper bound (e.g., v − i and v + i ).Based on the interval number distance method and the two different fuzziness parameters, m 1 and m 2 , the lower and upper membership grades of each sample can be expressed as follows: where i = 1, 2, . . ., C, and k = 1, 2, . . ., N.
The objective functions are similar to Equation (2): v L i and v R i are determined by the KM algorithm.Additionally, the iteration can be terminated when J t+1 m (U, v) − J t m (U, v) ≤ ε is satisfied.The possibility ranking method between interval numbers proposed by [65] was adopted in IT2FCM*.
The lower and upper membership grades of each sample in each class are expressed by an interval number vector { u 1k , u 2k , . . . ,u Ck } = u 1k , u 1k , u 2k , u 2k , . . ., u Ck , u Ck .The probability for any two intervals in the vector can be calculated as follows: where L( u ik ) = u ik − u ik and L u jk = u jk − u jk are the widths of the interval numbers u ik and u ik , respectively, for i, j = 1, 2, . . ., C, and k = 1, 2, . . ., N.
We can then obtain a possibility matrix P = (p ij , k).Moreover, the ranking vector w k = (w 1k , w 2k , . . . ,w Ck ) T can be calculated by and the index of the maximum value in w k is the class index of the sample.
In general, the improved IT2FCM* algorithm includes essentially four steps.Initially, the two fuzzifiers m 1 and m 2 and the termination criterion value ε are provided to initialize the lower and upper membership grade matrix u = [u, u] using a random method.Second, all the centroids v are computed and the membership grade matrix is updated.Here, the v L i and v R i relating to the centroids v are determined by the KM algorithm, the Euclidean distance between interval vectors is calculated using Equation ( 8), the lower and upper membership grade matrix u = [u, u] is calculated using Equations (9)  and (10), and then the objective function is calculated via Equation (1).
is satisfied, the algorithm moves to the next step, otherwise, it iteratively repeats the second step; Third, classify each sample using the interval number ranking method, here, the possibility matrix is calculated using Equation ( 11) and-based on the possibility matrix-the ranking vector w k = (w 1k , w 2k , . . . ,w Ck ) T is calculated, and a sample is assigned to a cluster according to the index of the maximum value in the ranking vector.Finally, the algorithm reports the clustering results.The detailed steps in the modified IT2FCM* algorithm are illustrated in Figure 4.
is calculated using Equation ( 8), the lower and upper membership grade matrix = , is calculated using Equations ( 9) and (10), and then the objective function is calculated via Equation (1).If ( , ) − ( , ) ≤ is satisfied, the algorithm moves to the next step, otherwise, it iteratively repeats the second step; Third, classify each sample using the interval number ranking method, here, the possibility matrix is calculated using Equation ( 11) and-based on the possibility matrix-the ranking vector = ( , , … , ) is calculated, and a sample is assigned to a cluster according to the index of the maximum value in the ranking vector.Finally, the algorithm reports the clustering results.The detailed steps in the modified IT2FCM* algorithm are illustrated in Figure 4.

Introduction of Four Validity Indexes
Four cluster validity indexes are selected in this paper to validate the separability of the IT2FCM* method.The PC-method proposed by Bezdek indicates the average relative amount of membership sharing between pairs of fuzzy subsets, in which higher values correspond to better classification results.More details can be found in [51,66].The PE-method, also proposed by Bezdek, is a scalar measure of the amount of fuzziness in a set of results, more details can be found in [51].The validity index FS, proposed by Fukuyama and Sugeno, is designed to measure the discrepancy between fuzzy compactness and fuzzy separation.More detailed information can be found in [67].The validity index XB-, proposed by Xie and Beni, measures the average within-cluster fuzzy compactness against the minimum between-cluster separation.More detailed information can be found in [68,69].Smaller values of PE-, FS-or XB-indicate better clustering performance.

Vegetation Index
A vegetation index can reflect the biophysical parameters of wetland vegetation, e.g., the vegetation cover and chlorophyll content, based on spectral reflectance characteristics in the red and NIR bands.This study adopted the NDVI to address the dynamics of vegetation and its impact on the wetland geomorphology.The equation for calculating the NDVI is as follows [70]:

Introduction of Four Validity Indexes
Four cluster validity indexes are selected in this paper to validate the separability of the IT2FCM* method.The PC-method proposed by Bezdek indicates the average relative amount of membership sharing between pairs of fuzzy subsets, in which higher values correspond to better classification results.More details can be found in [51,66].The PE-method, also proposed by Bezdek, is a scalar measure of the amount of fuzziness in a set of results, more details can be found in [51].The validity index FS, proposed by Fukuyama and Sugeno, is designed to measure the discrepancy between fuzzy compactness and fuzzy separation.More detailed information can be found in [67].The validity index XB-, proposed by Xie and Beni, measures the average within-cluster fuzzy compactness against the minimum between-cluster separation.More detailed information can be found in [68,69].Smaller values of PE-, FS-or XB-indicate better clustering performance.

Vegetation Index
A vegetation index can reflect the biophysical parameters of wetland vegetation, e.g., the vegetation cover and chlorophyll content, based on spectral reflectance characteristics in the red and NIR bands.This study adopted the NDVI to address the dynamics of vegetation and its impact on the wetland geomorphology.The equation for calculating the NDVI is as follows [70]: where ρ nir is the reflectance in the NIR band, and ρ red is the reflectance in the red band.

Experimental Results and Discussion
In this section, the main objectives are to compare the IT2FCM approach and the modified IT2FCM* approach based on four validity indexes to analyse the consistency between the results of these approaches and to further evaluate the separability of the modified IT2FCM* method.In addition, the dynamics of the Beidagang wetland over the last 40 years are analysed based on a long-term series of Landsat satellite data.

Separability Assessment Based on Four Validity Indexes for the Modified IT2FCM* Algorithm
To evaluate the separability of classes with the modified IT2FCM* algorithm, the main goal is to measure the consistency between IT2FCM* and the IT2FCM method.Thus, in this paper, four cluster validity indexes are adopted to assess the consistency: PC-, PE-, FS-, and XB-.As described in the previous section, the PC-value indicates the average relative amount of membership sharing between pairs of fuzzy subsets, thus, a higher PC-value corresponds to a better classification result, and smaller values of PE-, FS-or XB-indicate a better clustering performance.In this paper, to show the suitability and evaluate the performance of the improved IT2FCM* algorithm, different satellite remotely sensed datasets are used, including the EO-1 ALI dataset of 22 August 2015, the Landsat 8 OLS dataset of 29 May 2016, and the Sentinel-2 satellite dataset from 30 May 2016.The associated classification results were used as examples to demonstrate the changes in the four indexes (see Table 2).
Table 2 shows that the four indexes of the IT2FCM* method exhibit the same trends as the four indexes of the IT2FCM method.The PC-and PE-values exhibit slight changes because they account only for the membership degree (i.e., they exclude the clustering centre and source data).From Table 2, we can see that the PC-values of IT2FCM* are generally higher than those of IT2FCM.The most prominent change in these four indexes is the value of FS-, the FS-value of IT2FCM* is typically far smaller than that of IT2FCM.The values of XB-and PE-exhibit slight changes when comparing the two algorithms.According to previous results in the literature [67,71,72], the larger the PC-value is and the smaller the PE-, FS-and XB-values are, the better the classification results are.The abovementioned analysis also shows that the classification results of IT2FCM* are better than those of IT2FCM.The improved IT2FCM* algorithm achieves a better performance on different multi-spectral satellite remotely sensed data.

Separability Assessment Based on Accuracy of the Classification Results
In this section, the separability of classes with the IT2FCM* method was further tested based on the classification results of the different multi-spectral satellite remote sensing datasets (see Figure 5).These different satellite remotely sensed datasets are described in Section 3.For the ALI dataset, a different location around the study area of the Beidagang wetland was selected because the dataset lacked coverage (See Figure 5a,b).In contrast, the Landsat 8 and Sentinel-2 imagery corresponded to the study area (See Figure 5c,d).To evaluate the accuracy of the classification results, ground truth dataset was collected (See Figure 5b,e).In the study area of the Beidagang wetland, because of difficulty accessing the marshes, most of the field survey data were collected and delineated on both sides of the roads within the wetland (See Figure 5e).Then, the confusion matrix coefficients of these different satellite datasets and the accuracy of each class were evaluated based on the kappa coefficients and the overall accuracy.In the area displayed in Figure 5a, the main land cover types include water bodies, human communities and concrete roads, bare soil and small country trails, croplands and thick vegetation, and grasslands or sparse vegetation.In this study, the land cover types in this area were classified into 5 types: water, sparse vegetation, thick vegetation and cropland, bare soil and impervious surface (See Figure 5a).The overall accuracy of the results based on the ALI dataset was 89.55% with a kappa coefficient of 0.86.In the study area of the Beidagang wetland, the land cover types were classified into 6 classes based on the Sentinel-2 and Landsat 8 datasets: water, water with submerged vegetation, emergent vegetation over water, thick emergent vegetation over wetland, vegetation on land, and human communities and bare soil (See Figure 5c,d).The overall accuracy of the results based on the Sentinel-2 dataset was 89.40% with a kappa coefficient of 0.88, and the overall accuracy of the results based on the Landsat 8 OLS dataset was 84.55% with a kappa coefficient of 0.81.As shown in in Figure 5c,d, despite the differences between the classification results due to the different spatial resolution, but the results from Sentine-2 images maintain a good agreement with the results from Landsat 8 images.In general, the clustering results were good, which further demonstrates the separability of classes with the improved IT2FCM* algorithm.

Separability Assessment Based on Accuracy of the Classification Results
In this section, the separability of classes with the IT2FCM* method was further tested based on the classification results of the different multi-spectral satellite remote sensing datasets (see Figure 5).These different satellite remotely sensed datasets are described in Section 3.For the ALI dataset, a different location around the study area of the Beidagang wetland was selected because the dataset lacked coverage (See Figure 5a,b).In contrast, the Landsat 8 and Sentinel-2 imagery corresponded to the study area (See Figure 5c,d).To evaluate the accuracy of the classification results, ground truth dataset was collected (See Figure 5b,e).In the study area of the Beidagang wetland, because of difficulty accessing the marshes, most of the field survey data were collected and delineated on both sides of the roads within the wetland (See Figure 5e).Then, the confusion matrix coefficients of these different satellite datasets and the accuracy of each class were evaluated based on the kappa coefficients and the overall accuracy.In the area displayed in Figure 5a, the main land cover types include water bodies, human communities and concrete roads, bare soil and small country trails, croplands and thick vegetation, and grasslands or sparse vegetation.In this study, the land cover types in this area were classified into 5 types: water, sparse vegetation, thick vegetation and cropland, bare soil and impervious surface (See Figure 5a).The overall accuracy of the results based on the ALI dataset was 89.55% with a kappa coefficient of 0.86.In the study area of the Beidagang wetland, the land cover types were classified into 6 classes based on the Sentinel-2 and Landsat 8 datasets: water, water with submerged vegetation, emergent vegetation over water, thick emergent vegetation over wetland, vegetation on land, and human communities and bare soil (See Figure 5c,d).The overall accuracy of the results based on the Sentinel-2 dataset was 89.40% with a kappa coefficient of 0.88, and the overall accuracy of the results based on the Landsat 8 OLS dataset was 84.55% with a kappa coefficient of 0.81.As shown in in Figure 5c,d, despite the differences between the classification results due to the different spatial resolution, but the results from Sentine-2 images maintain a good agreement with the results from Landsat 8 images.In general, the clustering results were good, which further demonstrates the separability of classes with the improved IT2FCM* algorithm.

Separability Assessment Based on Analysis of the Membership of Classified Land Cover Types
In this section, the separability of the IT2FCM* algorithm was further tested based on the minimum and maximum membership values for each class by using the 2005 Landsat TM imagery as an example.Maps of the minimum and maximum membership values of the different land cover types based on the 2005 Landsat TM date set are shown in Figure 6.Notably, water and water mixed with emergent vegetation, for example, are rarely well classified, confirming the observations based on the fuzzy areas of the two classes.Even the lake does not exhibit a minimum fuzzy membership close to 0. In fact, Figure 5 shows that the maximum membership of each land cover is close to 1 in all the maps.This result indicates that all these classes were clearly differentiated, although the ranges of membership values decreased from the maximum to the minimum maps.Thus, in general, the membership values for each class also show the good separability of the improved IT2FCM* algorithm.

Separability Assessment Based on Analysis of the Membership of Classified Land Cover Types
In this section, the separability of the IT2FCM* algorithm was further tested based on the minimum and maximum membership values for each class by using the 2005 Landsat TM imagery as an example.Maps of the minimum and maximum membership values of the different land cover types based on the 2005 Landsat TM date set are shown in Figure 6.Notably, water and water mixed with emergent vegetation, for example, are rarely well classified, confirming the observations based on the fuzzy areas of the two classes.Even the lake does not exhibit a minimum fuzzy membership close to 0. In fact, Figure 5 shows that the maximum membership of each land cover is close to 1 in all the maps.This result indicates that all these classes were clearly differentiated, although the ranges of membership values decreased from the maximum to the minimum maps.Thus, in general, the membership values for each class also show the good separability of the improved IT2FCM* algorithm.

Forty-Year Spatial and Temporal Dynamics of the Beidagang Wetland Based on a Long-Term Series of Landsat Data
Based on the above analysis from different aspects, the improved IT2FCM* approach shows a good separability in multi-spectral satellite remote sensing image clustering.Using the IT2FCM* method, the spatiotemporal changes of Beidagang wetland geomorphology were detected based on long-term series of multi-spectral Landsat datasets.Based on the classification results, the study area is classified into six classes: water, water with submerged or floating-leaved vegetation, water mixed with emergent vegetation, thick emergent vegetation, thick vegetation on land, and bare soil or impervious area (see Figure 7).The normalized difference vegetation index (NDVI) maps are introduced in this section to address the influence of the dynamic of NDVI on the Beidagang wetland geomorphology (See Figure 8).To address the dynamics of the Beidagang wetland and the important factors that influence these dynamics, altogether, three factors are assessed: the dynamics of water, the dynamics of the vegetation index, and the dynamics of vegetation types.The three factors are introduced to further analyse their influences on the wetland and its degradation.

Forty-Year Spatial and Temporal Dynamics of the Beidagang Wetland Based on a Long-Term Series of Landsat Data
Based on the above analysis from different aspects, the improved IT2FCM* approach shows a good separability in multi-spectral satellite remote sensing image clustering.Using the IT2FCM* method, the spatiotemporal changes of Beidagang wetland geomorphology were detected based on long-term series of multi-spectral Landsat datasets.Based on the classification results, the study area is classified into six classes: water, water with submerged or floating-leaved vegetation, water mixed with emergent vegetation, thick emergent vegetation, thick vegetation on land, and bare soil or impervious area (see Figure 7).The normalized difference vegetation index (NDVI) maps are introduced in this section to address the influence of the dynamic of NDVI on the Beidagang wetland geomorphology (See Figure 8).To address the dynamics of the Beidagang wetland and the important factors that influence these dynamics, altogether, three factors are assessed: the dynamics of water, the dynamics of the vegetation index, and the dynamics of vegetation types.The three factors are introduced to further analyse their influences on the wetland and its degradation.

The Dynamics of Water and Analysis of the Associated Driving Factors
Water is a key component of wetlands.Hence, to assess the spatial dynamics of water and address the factors that influence these dynamics, we quantitatively analysed the area of water bodies in the classification results.Additionally, precipitation data were used in the analysis (see Figures 8 and 9).Figures 7-9, show that the water area in the Beidagang wetland decreased overall from 1975 to 2015.The area of the water bodies in 1976 and 1995 were greater than in other years, especially after 1995, when the area of the water bodies obviously decreased.In 2010, and 2015, the area of water bodies occupied less than one-third of the study area.There are multiple reasons for this trend.Since 1971-and especially since 1980-the annual precipitation in the study area has averaged approximately 550 mm, which represents a decrease of 8.9 mm per decade, while annual evaporation has averaged 1777 mm.Thus, based on Figures 7-9, we can conclude that the water area trend is similar to the trends in annual precipitation and summer precipitation.In addition, the quantity of water resources has decreased since 2000; for example, the Duliujian River (see Figure 1c) had an average annual flow of 1 billion m 3 into the Bohai Sea (Western Pacific) in the 1970s but only 50 million m 3 in the 1990s and less than 10 million m 3 since 2005 (from the Tianjin Bureau of Statistics).In addition, the increasing demands for water from agriculture and industry have influenced the overall quantity of water resources.For example, the Beidagang Oil factory, which is located to the southeast of Beidagang Reservoir, and other large oil refineries and chemical factories (i)

The Dynamics of Water and Analysis of the Associated Driving Factors
Water is a key component of wetlands.Hence, to assess the spatial dynamics of water and address the factors that influence these dynamics, we quantitatively analysed the area of water bodies in the classification results.Additionally, precipitation data were used in the analysis (see Figures 8 and 9).Figures 7-9, show that the water area in the Beidagang wetland decreased overall from 1975 to 2015.The area of the water bodies in 1976 and 1995 were greater than in other years, especially after 1995, when the area of the water bodies obviously decreased.In 2010, and 2015, the area of water bodies occupied less than one-third of the study area.There are multiple reasons for this trend.Since 1971-and especially since 1980-the annual precipitation in the study area has averaged approximately 550 mm, which represents a decrease of 8.9 mm per decade, while annual evaporation has averaged 1777 mm.Thus, based on Figures 7-9, we can conclude that the water area trend is similar to the trends in annual precipitation and summer precipitation.In addition, the quantity of water resources has decreased since 2000; for example, the Duliujian River (see Figure 1c) had an average annual flow of 1 billion m 3 into the Bohai Sea (Western Pacific) in the 1970s but only 50 million m 3 in the 1990s and less than 10 million m 3 since 2005 (from the Tianjin Bureau of Statistics).In addition, the increasing demands for water from agriculture and industry have influenced the overall quantity of water resources.For example, the Beidagang Oil factory, which is located to the southeast of Beidagang Reservoir, and other large oil refineries and chemical factories

The Dynamics of Water and Analysis of the Associated Driving Factors
Water is a key component of wetlands.Hence, to assess the spatial dynamics of water and address the factors that influence these dynamics, we quantitatively analysed the area of water bodies in the classification results.Additionally, precipitation data were used in the analysis (see Figures 8  and 9).Figures 7-9, show that the water area in the Beidagang wetland decreased overall from 1975 to 2015.The area of the water bodies in 1976 and 1995 were greater than in other years, especially after 1995, when the area of the water bodies obviously decreased.In 2010, and 2015, the area of water bodies occupied less than one-third of the study area.There are multiple reasons for this trend.Since 1971-and especially since 1980-the annual precipitation in the study area has averaged approximately 550 mm, which represents a decrease of 8.9 mm per decade, while annual evaporation has averaged 1777 mm.Thus, based on Figures 7-9, we can conclude that the water area trend is similar to the trends in annual precipitation and summer precipitation.In addition, the quantity of water resources has decreased since 2000; for example, the Duliujian River (see Figure 1c) had an average annual flow of 1 billion m 3 into the Bohai Sea (Western Pacific) in the 1970s but only 50 million m 3 in the 1990s and less than 10 million m 3 since 2005 (from the Tianjin Bureau of Statistics).In addition, the increasing demands for water from agriculture and industry have influenced the overall quantity of water resources.For example, the Beidagang Oil factory, which is located to the southeast of Beidagang Reservoir, and other large oil refineries and chemical factories located to the northwest of the Beidagang wetland (see Figure 7f corresponding to Figure 1c) require large quantities of water.Finally, rapid urbanization and farmland reclamation from the wetland have decreased the area of water bodies, and the construction of retaining walls has caused a decrease in intertidal zones.The spatial differences in the water area dynamics are obvious and clearly illustrated in Figures 7 and 8.The natural temporal changes in the water area of Beidagang Reservoir from 1975 to 2015 are notable: the reasons for these changes were addressed above.In addition, large spatial changes in the water area of the Beidagang wetland have occurred since 1995 (see Figure 7f-i).Figure 7f-i illustrate the areas of water in the north of the study area.These areas exhibit unnatural patterns.According to the previous literature [39] and field observations, these water bodies have been manually developed into fish farms.In addition, the conversion rate of natural waters to fishponds has increased in recent decades: before 1995, almost no fishponds existed in this region.The development of fishponds has endangered not only the local water resources, which are key wetland components but also the vegetation in the study area.

The Dynamics of the Vegetation Index
The long-term changes in the NDVI series are shown in Figure 10a-i.These figures show that there is a close relationship between vegetation coverage and the water area and depth.The value of NDVI in the wetland is relatively high compared to the value of NDVI on land except in 1995.This trend occurred because, in 1995, the water area was so large that it nearly inundated the study area.Another interesting phenomenon illustrated by these figures is that the NDVI of some water areas is greater than 0.2 because submerged plants grow in the shallow water areas.Comparing the NDVI with the classification results from 1975 to 2015, especially those in 1985, 1990 and 2000, this phenomenon is obvious.The spatial differences in NDVI are also illustrated in these figures, especially in 2015, which varies considerably in the western part of the study area (corresponding to area D delineated by a black circle in Figure 1c) due to the agricultural land development.However, the NDVI values in 1976 and 1980 differ from those in other years.Notably, many negative values can be observed in these years, especially in the bare soil and built-up area classes.However, these results may be associated with the dataset used in this study.The NDVIs in 1976 and 1980 were calculated from Landsat 2 MSS data, which are associated with different spectral channel settings in the NIR and red bands compared to other satellite-based datasets.The spatial differences in the water area dynamics are obvious and clearly illustrated in Figures 7  and 8.The natural temporal changes in the water area of Beidagang Reservoir from 1975 to 2015 are notable: the reasons for these changes were addressed above.In addition, large spatial changes in the water area of the Beidagang wetland have occurred since 1995 (see Figure 7f-i).Figure 7f-i illustrate the areas of water in the north of the study area.These areas exhibit unnatural patterns.According to the previous literature [39] and field observations, these water bodies have been manually developed into fish farms.In addition, the conversion rate of natural waters to fishponds has increased in recent decades: before 1995, almost no fishponds existed in this region.The development of fishponds has endangered not only the local water resources, which are key wetland components but also the vegetation in the study area.

The Dynamics of the Vegetation Index
The long-term changes in the NDVI series are shown in Figure 10a-i.These figures show that there is a close relationship between vegetation coverage and the water area and depth.The value of NDVI in the wetland is relatively high compared to the value of NDVI on land except in 1995.This trend occurred because, in 1995, the water area was so large that it nearly inundated the study area.Another interesting phenomenon illustrated by these figures is that the NDVI of some water areas is greater than 0.2 because submerged plants grow in the shallow water areas.Comparing the NDVI with the classification results from 1975 to 2015, especially those in 1985, 1990 and 2000, this phenomenon is obvious.The spatial differences in NDVI are also illustrated in these figures, especially in 2015, which varies considerably in the western part of the study area (corresponding to area D delineated by a black circle in Figure 1c) due to the agricultural land development.However, the NDVI values in 1976 and 1980 differ from those in other years.Notably, many negative values can be observed in these years, especially in the bare soil and built-up area classes.However, these results may be associated with the dataset used in this study.The NDVIs in 1976 and 1980 were calculated from Landsat 2 MSS data, which are associated with different spectral channel settings in the NIR and red bands compared to other satellite-based datasets.(a-i) correspond to 1976,1980,1985,1990,1995,2000,2005,2010, and 2015 NDVI results.

Changes in Vegetation Types Based on Water Resource Dynamics
Figure 6 shows that the quantity of natural water resources is decreasing rapidly and that the salt content is increasing, which is important for the evolution of vegetation [39].The number of manmade pools is increasing slowly; and the area of bare land is increasing due to urban development.To assess the evolution of vegetation types in detail, a field survey and observations were performed on 27 June 2016.The geomorphology of the Beidagang wetland is mainly composed of clay, with poor soil and a high salt content.The vegetation in this region is dominated by salt-tolerant vegetation [39].Through the field survey and observations on 27 June 2016 (see Figure 9), we found that the main vegetation is Phragmites australis (P.australis).There are patchy distributions of other vegetation types, including Suaeda salsa and Scirpus planiculmis communities.Almost no woody vegetation exists; only scattered Tamarix chinensis and Nitraria sibirica were observed.Figure 7a-i show that the vegetation is zonally distributed, which is reflective of the spatial evolution.Figure 7 shows that wet vegetation types have large NDVI values.The moderately wet vegetation mainly includes P. australis-Suaeda and Suaeda salsa.As a typical halophyte, Suaeda salsa can grow in high salt content areas, especially in the intertidal zones of salt wetlands.Additionally, it can transform saline ecological environments due to its ability to absorb the salt from the soil.Such processes can improve the soil conditions and allow other vegetation types to grow in the local environment.The xeromorphic vegetation mainly consists of P. australis-Tamarix chinensis communities, both of which are salt-tolerant plants, and some Artemisia annua-Artemisia scoparia communities (which are salt-rejecting plants) in certain land areas, especially those far from water areas.According to the field observations, the submerged vegetation class mainly includes Potamogeton pectinatus, Potamogeton crispus, Najas marina, Myriophyllum verticillatum, and some Ceratophyllum demersum communities.Aquatic vegetation grows close to or in the water and mainly includes P. australis and Scirpus planiculmis communities.Both types of vegetation, which are typical emergent vegetation communities that thrive in water environments, are water resistant and salt tolerant.The Phragmites australis-Suaeda communities are the main types of wet vegetation.They grow in and close to saline wetland and lake areas, such as in saline marshes and other areas with wet conditions and environments and.
The spatiotemporal changes in different types of vegetation suggest that the soil properties underwent significant changes, and the salt content decreased with the water content.With the vegetation types varying from aquatic vegetation to wet vegetation to xeromorphic vegetation, the salt-tolerant plants slowly shifted to salt-secreting and salt-rejecting plants.The spatiotemporal changes in vegetation types reflected not only the temporal gradient of the evolution of different vegetation but also the changes in the water content.Additionally, the salt-tolerance capability of vegetation decreased, resulting in some communities with mainly salt-rejecting vegetation, such as Artemisia annua-Artemisia scoparia.Based on the field observations, among these different types of vegetation, P. australis was widely distributed, implying that P. australis plants are adaptive to a wide range of environments.As the dominant species, P. australis was widely distributed in water areas with different salt contents, and the species accounted for as much as 90 percent of all plant species, forming a "sea of reeds".

Conclusions
In this study, spatiotemporal wetland cover changes were analysed using IT2FCM*, an improved approach proposed by this paper.The IT2FCM* based on interval type 2 fuzzy set was used for multi-spectral image analysis.Compared with other traditional fuzzy methods used in image clustering, the interval number ranking was adopted into this improved IT2FCM* algorithm.This addition improves its ability to address spectral uncertainty.Then, to evaluate the separability of the improved IT2FCM* algorithm, multi-spectral EO-1 ALI, Sentinel-2 and Landsat 8 OLS satellite remotely sensed datasets were used, and four validity indexes, including PC-, PE-, FS-and XB-, were introduced.The values of these four indexes of three satellite remote sensing data perform well in image clustering.To further assess the separability of the improved IT2FCM* approach, the classification results of the abovementioned satellite images were validated with ground truth data collected by a field survey.High overall accuracy and the kappa coefficients show that the IT2FCM* has good performance, which also demonstrated the separability of the algorithm.In addition, the maximum and minimum membership values were used to evaluate the separability of the algorithm; these results also indicate the algorithm's good separability.
Based on the above analysis of the separability of the improved IT2FCM* algorithm, the spatiotemporal changes of the Beidagang wetland geomorphology were analysed.According to spatial and temporal analyses of different land cover types in the study area, the water area has steadily decreased for many reasons including climate conditions (especially precipitation, evaporation, air and land surface temperatures, and sunshine duration), environmental factors (such as the salt content of water and the effects of and anthropogenic effects (including the construction of buildings, land use-land cover changes, and industrial and agricultural development).All these factors have influenced the studied wetland.The numbers of vegetation types, fish species, and birds, including wading birds, are decreasing.These trends will continue if the conditions do not improve in the future.For example, Phragmites and Typha angustifolia were the two dominant plant species before 2000, but the 2016 survey showed that Phragmites and Suaeda salsa have become dominant.Vegetation types underwent dramatic changes.As the volume of water resources decreased, high soil salt content has increased due to the reduction in surface runoff to the Duliujian River and more evaporation than precipitation.In turn, these changes will affect the growth of vegetation and other living creatures.The number of benthic species and the algal area have decreased notably.There were 14 benthic species in 1980, and their average density was 598 ind•m −2 in 1980; however, in 2014, there were only 3 benthic species, and their average density was just 56 ind•m −2 .Moreover, migratory ocean fish have completely disappeared from the Beidagang wetland.In contrast, farmed fish now account for the majority of fish species observed.The other effects produced by these changes require further analysis.
However, we should be aware of the limitations and advantages of satellite remote sensing techniques and multispectral and moderate spatial resolution satellite data in detecting wetland land cover changes.Although satellite remote sensing cannot provide detailed information that could be obtained via a field survey, remotely sensed data can be used to continuously monitor wetlands over large areas and provide researchers with needed complementary information.To better use remote sensing techniques to monitor wetland changes, additional data should be considered.Aerial photographs or field survey data from different years should be collected to strengthen such analyses and help evaluate the results.Unfortunately, such data are rarely collected because of economic limitations and the lack of the necessary equipment.

Figure 1 .
Figure 1.Study area of Tianjin Municipality, North China: (a) the location of Tianjin in China; (b) the location of the study area in Tianjin; and (c) the remotely sensed image obtained by the Landsat 8 satellite (acquired on 2 October 2015) with 7/5/3 bands in the RGB composite colour image.

Figure 1 .
Figure 1.Study area of Tianjin Municipality, North China: (a) the location of Tianjin in China; (b) the location of the study area in Tianjin; and (c) the remotely sensed image obtained by the Landsat 8 satellite (acquired on 2 October 2015) with 7/5/3 bands in the RGB composite colour image.

Figure 2 .
Figure 2. (a-f) The field observations on 27 June 2016; a through e correspond to some photos of different vegetation types and f is the oil finery which corresponds to region C (in the red circle) in Figure 1c.

Figure 2 .
Figure 2. (a-f) The field observations on 27 June 2016; a through e correspond to some photos of different vegetation types and f is the oil finery which corresponds to region C (in the red circle) in Figure 1c.

Figure 3 .
Figure 3. Flowchart of the spatiotemporal analysis of the wetland based on long-term series of Landsat satellite data.

Figure 3 .
Figure 3. Flowchart of the spatiotemporal analysis of the wetland based on long-term series of Landsat satellite data.

Figure 4 .
Figure 4.The flowchart of the modified IT2FCM* algorithm.

Figure 4 .
Figure 4.The flowchart of the modified IT2FCM* algorithm.

Figure 5 .
Figure 5. Classification results from different satellite multi-spectral datasets to further demonstrate the separabilty of the improved IT2FCM*: (a) the classification result from the ALI image of 22 August 2015; (b) the ground truth data collected by the field survey on 27 June 2016; (c) the classification result of the Sentinel-2 dataset acquired on 30 May 2016; (d) the classification result of the Landsat 8 remote sensing imagery acquired on 29 May 2016; and (e) the ground truth data collected by the field survey on 27 June 2016.

Figure 6 .
Figure 6.Maps of the type 2 fuzzy land cover types of the study area based on long-term series of Landsat satellite data starting from 4 September 2005.The land cover types in each year are divided into two columns: one is the cover type with the maximum membership, and the other is the cover type with the minimum membership.

Figure 6 .
Figure 6.Maps of the type 2 fuzzy land cover types of the study area based on long-term series of Landsat satellite data starting from 4 September 2005.The land cover types in each year are divided into two columns: one is the cover type with the maximum membership, and the other is the cover type with the minimum membership.

Figure 8 .
Figure 8.The relationship between the water area and precipitation (annual and summer precipitation).

Figure 8 .
Figure 8.The relationship between the water area and precipitation (annual and summer precipitation).

Figure 8 .
Figure 8.The relationship between the water area and precipitation (annual and summer precipitation).
Remote Sens. 2017, 9, 683 18 of 24 located to the northwest of the Beidagang wetland (see Figure7fcorresponding to Figure1c) require large quantities of water.Finally, rapid urbanization and farmland reclamation from the wetland have decreased the area of water bodies, and the construction of retaining walls has caused a decrease in intertidal zones.

Figure 9 .
Figure 9.The annual precipitation and summer precipitation in the study area.

Figure 9 .
Figure 9.The annual precipitation and summer precipitation in the study area.

Table 1 .
Summary of long-term series of Landsat satellite dataset used in this study.Percent cloud cover corresponds with the date/time of remotely sensed data collection.

Table 2 .
Comparison of validity indexes with from different satellite remote sensing data.