Mapping the Spatiotemporal Pattern of Sandy Island Ecosystem Health during the Last Decades Based on Remote Sensing

: Sandy islands suffer multiple external disturbances and thus experience drastic temporal ecological changes, and the two parts, that is, the states of essential components (Part 1) and the resilience under multiple disturbances (Part 2), are both indispensable to the sandy island ecosystem health. In this study, a model for the sandy island ecosystem health was established by integrating the two parts. In Part 1, the states were measured following the framework of vegetation, soil, and landscape, and a total of 12 factors in the three components were adopted. In Part 2, three typical disturbances, namely, geomorphological change, soil salinization, and human inﬂuence, were identiﬁed, and the resilience across different time intervals was measured by clarifying the intrinsic correlations between the components and disturbances. A sandy island ecosystem health index (SIEHI) was proposed based on the two parts, and Chongming Island was selected as the study area to demonstrate the model. The results indicated that the SIEHI continuously increased from 1988 to 2017, denoting the good effects of “Eco-Island construction” on Chongming Island. In different components and factors, the vegetation and soil components contributed more than the landscape component to the sandy island ecosystem health, and vegetation quality and soil carbon/nitrogen were the factors that made the most contributions. In different disturbances, the human inﬂuence played a major role in driving the spatiotemporal variations of the sandy island ecosystem health. Farming and building construction contributed the most and accounted for 37.12% and 35.59% of the total human inﬂuence, respectively, while trafﬁc development exerted the highest inﬂuence per area. Then, inﬂuence coefﬁcients of different human activities on the sandy island were determined, and measures for different functional zones were proposed for balancing the protection and development and achieving the sandy island ecosystem-based management.


Introduction
The sandy island is one of the three main types of islands according to the material constitution (the other two are rocky and coral reef islands), and by definition, this type of island is mainly constituted by sand and mud. It is always located in the estuarine areas and formed and changed by the sediment input via the river and deposition in the estuary. Sandy islands are in a low number compared with the total number of islands in China (more than 10,000) yet occupy a considerable proportion of the island areas and become an integral part of the coastal areas for human inhabitation, economic development, and ecological conservation [1]. Saltmarshes and mangroves are always developed on sandy islands and exhibit enormous ecological functions, including biodiversity maintenance, habitat provision, pollutant purification, erosion prevention, and carbon sequestration [2][3][4][5][6][7]. Meanwhile, flat terrain, fertile land, and superior geographic location render the sandy island a hotspot for human activities. Humans live and conduct economic activities by utilizing various resources of sandy islands. In China, urban and rural constructions, coastal reclamation, farming, pond aquaculture, and tourism activities are frequent and Figure 1. Framework for assessing the sandy island ecosystem health: VA: vegetation abundance (dimensionless); NDVI: normalized difference vegetation index (dimensionless); NPP: net primary productivity (g C m −2 a −1 ); H': Shannon-Wiener index (dimensionless); BD: bulk density (g/cm 3 ); TOC: total organic carbon (g/kg); TN: total nitrogen (g/kg); AK: available potassium (mg/kg); ILC: important landscape coverage (dimensionless); NP: number of patches (dimensionless); TE: total edges (m); LII: landscape isolation index (dimensionless).

Study Area
Chongming Island is in the estuary of the Yangtze River, adjacent to the East China Sea, and close to Shanghai City ( Figure 2). It is the third largest island in China and the largest estuarine sandy island in the world, and exhibits typical features of estuarine sandy islands, that is, changeable geomorphology, abundant ecological functions, intensive land-sea interaction, and increasing human activities [5,23,67]. Chongming Island is the location of important ecological function areas, including Dongtan Wetland, Beihu Wetland, Xisha Wetland, and Dongping Forest Park (Figure 2), of which Dongtan Wetland was designated as a national nature reserve in 2005. Extensive studies on Chongming Island's ecosystem have been conducted. Current studies covered different aspects, including geomorphology evolution [13,68], biological resources [7,69], soil quality [70][71][72], landscape character [5,73], as well as ideas, policies, and measures for constructing an ecoisland [12,74,75]. Furthermore, evaluations on the island ecosystem were also implemented using different methods [24,76]. However, the spatiotemporal characteristics of ecosystem health that integrate the states of essential components and the resilience under multiple disturbances are still unclear. The spatial characteristics refer to the spatial heterogeneity of ecosystem health within the island based on the assessment unit, which is the minimum unit to conduct the assessment; the temporal characteristics indicate the temporal changes of ecosystem health across different years during the last decades.
The island has a subtropical, oceanic, monsoon climate with an average annual temperature and rainfall of 15.2 • C and 1025 mm, respectively [5]. Most of the island belongs to Chongming District in Shanghai City. In 2021, the gross domestic product (GDP) of the district summed up to CNY 40.97 billion (USD 5.76 billion), and the GDP per capita was CNY 60,889 (USD 8555); at the end of 2021, the total population summed up to approximately 672,900 [78]. A small part of the island in the northwest belongs to Nantong City in Jiangsu Province because of the historical issue. The spatial scope of this study is the entire Chongming Island, and the temporal scope is the last 30 years. Specifically, the entire Chongming Island in the 1988,1995,2007, and 2017 were selected as the objects, and the spatiotemporal changes over the island and across the years were studied.

Data Source
The assessment of the sandy island ecosystem health needs the area data that cover the entire island and range from 1988 to 2017. The needed data for the assessment were sourced by remote sensing and field investigation, of which remote sensing served as the major data source.
The remote sensing data from Landsat 5 in 1988Landsat 5 in , 1995Landsat 5 in , and 2007 and that from Landsat 8 in 2017 were collected on the premises of the consistency in season, the cleanliness of images with no or little clouds, and the time interval of approximately 10 years. The spatial resolution of the multispectral bands is 30 × 30 m. The spectral reflectance of each band in different years was extracted through radiometric calibration and atmospheric correction, and reflectance-based ecological indices were then calculated. Furthermore, island outlines and land cover types in different years were obtained through visual interpretation and supervised classification methods, and the land covers were divided into wetland vegetation, mudflat, woodland, farmland, water area, pond, building land, traffic land, and industrial land ( Figure 2). The classification accuracies were assessed using error matrices, and high accuracies were validated in a previous study of the authors [5]. The land cover data serve as the basic data for calculating the factors in vegetation and landscape components. Location, sampling sites, island outlines, and land cover types of the study area: the sampling sites, island outlines, and land cover types were quoted from previous relevant studies by the authors [5,71]. The protected areas were sourced from the People's Government of Chongming District in Shanghai [77].
The field investigation was conducted in May 2018. A total of 110 sampling sites were set based on the even distribution in the land surface, the representativeness of plant communities, and the accessibility to the sites. The plant community in each site was surveyed, and plant species data, including abundance, height, and coverage of each species in tree, shrub, and herb layers, were recorded. The surface (0-20 cm) soil was sampled in each site, and soil physicochemical properties, including bulk density (BD), total organic carbon (TOC), total nitrogen (TN), available phosphorus (AP), and available potassium (AK), were measured in the laboratory. The Part 1 of the sandy island ecosystem health is the states of essential components of the ecosystem. Vegetation, soil, and landscape are the three essential components and refer to the overall characteristics, vitality, and base of the island ecosystem, respectively [41,[79][80][81]. The three essential components involve different aspects of the island ecosystem from macro-to micro perspectives and from upper to bottom layers. In terms of the perspectives, the vegetation and soil indicate the specific ecological functions from the micro perspective, including biodiversity maintenance, habitat provision, carbon sequestration, and water conservation [79,[81][82][83]; the landscape refers to the compositional structure and spatial configuration of land surface elements from a macro perspective [5,84,85]. In terms of the spatial layers, landscape, vegetation, and soil are situated at the upper, intermediate, and bottom layers over the island surface, and can comprehensively represent the essential components and overall conditions of the island ecosystem. Some other natural factors such as geological background, climate condition, water resource, and animals were not considered. Geological background and climate condition are always spatially homogeneous over an island, and they are the background conditions rather than the indicators of the ecosystem health. Naturally developed surface runoffs are rare on islands because of the small catchment areas. The existing water areas are mostly artificial and have been considered in the landscape component; thus, the water resource was not selected as an independent component. For biological resources, only vegetation was selected, and other biological resources such as animals and microorganisms were not included in the framework. The vegetation is the main producer of the ecosystem; it refers to various ecological functions and largely represents the habitat quality of islands. Moreover, the investigation of vegetation is relatively simple with low uncertainty, and some parameters can be easily monitored through remote sensing in the long term. By contrast, difficulty and cost for investigating and monitoring animals and microorganisms are high, with distinct uncertainty. Thus, they were not considered to ensure the applicability of the framework. In previous studies, the three components have been successfully used in the assessments of island ecosystem health on rocky islands [10,41]. In this study, the framework of these three components was adopted to validate its effectiveness in the sandy island. In the factor selection for each component, the factors should be typical for specific aspects of the component and different factors should be mutually independent, and a total of 12 factors were selected under the framework (Figure 1).
The vegetation component was assessed in two aspects, namely, quantity and quality. Vegetation abundance (VA) refers to the quantity of the vegetation in the assessment unit, and normalized difference vegetation index (NDVI), net primary productivity (NPP), and Shannon-Wiener index (H') denote the quality. The NDVI directly represents the vegetation growth condition and neglects the difference among different vegetation types, while NPP considers the growth condition and the difference in light use efficiencies among vegetation types [86,87]. The H' refers to the species composition and diversity and indicates the stability of the island ecosystem [82,88]. The VA was calculated using the following equation: where VA i is the area of vegetation type i, and i = 1, 2, and 3 refer to wetland vegetation, woodland, and cropland, respectively; w i is the weighted value of vegetation type i and indicates the existence time of vegetation type i over the year; TA is the total area of the assessment unit. Wetland and woodland are covered by vegetation throughout the year, and thus w 1 and w 2 are given as 1; by contrast, cropland is uncovered after harvest, and thus w 3 is given as 0.5. NDVI was obtained through band calculation of the spectral reflectance using the following equation: where Re x is the spectral reflectance of band x. The serial number (x) of bands in the equation is for Landsat 8 data and should be (x − 1) for Landsat 5 data. NPP was estimated using a Carnegie-Ames-Stanford approach [89]. H' was calculated based on shrub and herb species because the tree species are mostly artificial and in a low number, and the method has been reported in [69]. Of the four factors, H' was sourced by the field investigation and the remaining three factors were sourced by remote sensing. The soil component was measured based on the physicochemical properties. Four soil factors, that is, BD, TOC, TN, and AK, were selected. The BD is an important physical factor that reflects the structural support function [90]. The TOC indicates the soil carbon storage; TN and AK refer to different nutrients; and all the three factors involve biogeochemical cycles [41,91]. All of the four factors are sourced by field investigation.
The landscape component consists of composition and configuration. The former refers to the land cover types and their areas in an assessment unit, while the latter focuses on the spatial distributions and inter-relations of different landscape patches in the assessment unit [5]. Four factors, namely, important landscape coverage (ILC), number of patches (NP), total edges (TE), and landscape isolation index (LII), were used. The ILC denotes the landscape composition, and the remaining three factors indicate the landscape configuration [5]. They are frequently used factors for representing the landscape patterns on different islands [5,10,92]. The ILC was calculated using the equation in [10]. The wetland vegetation was selected as the very important landscape because of its enormous ecological values; mudflat, woodland, and water area also provide specific ecological functions, yet to a lesser degree; thus, they were considered importance landscapes. The NP, TE, and LII were calculated using the methods reported in [5]. All of the four factors are obtained based on the land cover data sourced by remote sensing.
The assessment unit was set as 100 m × 100 m to balance the demands of spatial heterogeneity exhibition, landscape factor calculations, and appropriate data amount [11]. The assessment units in the four years were generated using the Fishnet tool in ArcGIS 10.0 (ESRI, Redlands, CA, USA).

Spatial Simulations of Field Data
In this study, the factors that are sourced by filed investigation include H', BD, TOC, TN, and AK. They are point data at a time point. The historical data from field investigation are always discontinuous and unavailable. Thus, these factors should be converted from point data at a time point to multitemporal area data to meet the assessment requirement. The method of spatial simulation can help realize the conversion. It is always conducted by establishing the relationships between field and remote sensing data and then simulating the spatial distributions of the field data based on the relationships [71,93,94]. In previous studies, the spatial simulations that convert the point data to area data at the same time point have been widely conducted [71,[95][96][97], and the spatial simulations that aim to convert the data "from present to past" were also explored [72]. The previous studies achieved good results and provided necessary methods for the conversion in the present study. The reliable and continuous data provided by Landsat series satellites lay the foundation for the cross-time simulation. The spectral reflectance and reflectance-based ecological indexes at a time point can correspond to those at another time point after radiometric calibration and atmospheric correction, and thus the simulations can be conducted from point to area and from present to past based on the remote sensing data. The detailed processes are as follows: First, a total of 15 predictors were selected following a predictor system that covers different aspects of the land surface characteristics [71]. The predictor system consists of two aspects, namely, spectral values and ecological indexes. The spectral values refer to the original spectral reflectance of each band. According to the spectral range of each band in Landsat 8 and Landsat 5 data, Re x (for Re 2 -Re 6 ) in Landsat 8 data corresponds to Re (x−1) in Landsat 5 data, while Re 7 in Landsat 8 data corresponds to Re 7 in Landsat 5 data. The ecological indexes are indexes that possess specific ecological significance and are calculated based on the spectral reflectance. Various types of ecological indexes exist and represent different aspects of the ecological conditions [71]. Vegetation, salinity, humidity, and aridity indexes indicate the vegetation growth condition, soil salinity, moisture degree, and dry degree, respectively [98][99][100]. Vegetation indexes consist of difference vegetation index (DVI), NDVI, and soil-adjusted vegetation index (SAVI); salinity indexes include salinity index 1 (SI1), salinity index 2 (SI2), and salinity index 3 (SI3); humidity index indicates land surface wetness (LSW); and aridity indexes are composed of index-based built-up index (IBI) and bare soil index (BSI). NDVI was calculated using Equation (2), and the remaining indexes were calculated using the following equations: where Re x indicates the spectral reflectance of Band x and Band (x − 1) in Landsat 8 and Landsat 5 data, respectively, for DVI, SAVI, SI1, SI2, SI3, IBI, and BSI; L is given as 0.5; LSW Lt8 and LSW Lt5 are the LSWs using Landsat 8 and Landsat 5 data, respectively; C x are coefficients of Re x for the calculation of LSW and are sourced from [101] and [102] for LSM Lt8 and LSM Lt5 , respectively. The predictors in all the assessment units in different years were obtained through ArcGIS 10.0. The correlations of the field data with the 15 predictors were analyzed using IBM SPSS 18 (IBM, Armonk, NY, USA), and certain significant correlations were observed (Table 1). All the simulated factors, that is, H', BD, TOC, TN, and TK, possessed significant correlations with different predictors, indicating the feasibility of their spatial simulations. In addition, another two factors, namely, Pielou index in the vegetation component and AP in the soil component, were also initially selected as the potential factors of the sandy island ecosystem health. However, no significant correlations of the two factors with the predictors existed, and thus the two factors were eliminated from the factors. The soil salinity (Sa), which would be used for indicating one of the disturbances in the next step, was originally point data sourced by filed investigation. It was also included in the simulated objects and its significant correlations with the predictors were found.
Second, 20% of the 110 sampling sites were randomly selected as the validating samples, and the other 80% were considered training samples. Before conducting the simulations, the outliers for each simulated factor were identified and removed. The partial least square regression (PLSR), which is a commonly used and highly repeatable approach [71,103,104], was used to implement the simulations using Minitab 17. The training and validating samples were entered, the maximum component number was given as 10, and a leave-one-out cross validation method was adopted. The results indicated the acceptable simulation accuracies for all the simulated factors. Third, the predictor values of all the assessment units in different years were entered, and the area data of the simulated factors were obtained, i.e., the simulations from point to area and from present to past were accomplished.
The spatial distributions of the 12 factors in the three components are shown in Figure 3.

Calculation of SIEHI-1
For measuring the states of different components and calculating the SIEHI-1, all the factor values should be standardized for the comparability of different factors in different years. For each of the factors, the 95th and 5th percentiles of all the factor values in the four years were considered the upper and lower limits, respectively, to avoid the extreme value effects. Of the 12 factors, BD, NP, TE, and LII are "less is better" factors, and the remaining 8 factors are "more is better" factors. The standardization was then conducted using the equation in [5].
The states of the three components were measured using the following equation: where C is the state of one component; F i is the standardized value of factor i contained in the component; and w i is the weighed value of factor i. Then, the SIEHI-1 was calculated using the following equation: where C i and w i are the state and weighed value of component i, respectively. An equal weighted method was adopted to represent the equal importance of the factors and components [5,10]. The f s is the function to standardize the value of SIEHI-1 using the 99th and 1st percentiles of all the index values in the four years as the upper and lower limits, respectively, for clearly exhibiting the spatiotemporal heterogeneity of the sandy island ecosystem health.

Part 2 of the Sandy Island Ecosystem Health: The Resilience
The framework for assessing the resilience is presented in Figure 4. The method was first proposed in the present study to assess the resilience of sandy ecosystem health by integrating the results of previous studies. In the method, three main external disturbances of the sandy island ecosystem were identified, the intrinsic correlation of the variation of each component with that of each disturbance was clarified, and the resilience was then measured based on the variations of component and disturbance, as well as their intrinsic correlations.

Three Disturbances
The sandy island ecosystem is deeply influenced by the intensive interactions of river, land, and sea because the formation and change of the sandy island are closely related to the water and sediment discharge via the river, the hydrodynamic condition in the estuary, and the continuous force from the seawater [13,71,105]. Furthermore, intensive human activities exert increasing influences on the sandy island ecosystem [5,14,69,106]. Therefore, multiple natural and anthropogenic disturbances jointly drive the change of sandy island ecosystem health. The main disturbances were identified in the following three types.
(1) Geomorphological change The geomorphological change mainly refers to the variations of area and shape of the island. It is the most conspicuous sign of island evolution and affects the basic parameter of an island, i.e., the island area, which determines the carrying capacity for biota and human activities [92,107]. It is also a fundamental parameter of coastal wetlands [52]. In terms of a sandy island, the increase in the island area indicates the increase in land resources and spaces for performing more functions and thus enhances the ecosystem health, and vice versa. The disturbance of geomorphological change across a time interval was measured using the following equation: (14) where D1 is the disturbance of geomorphological change; ∆D1 a-b is the variation of this disturbance across years a and b, and a-b indicates 1988-1995, 1995-2007, and 2007-2017 in this study; and A b and A a are island areas in an assessment unit in years b and a, respectively. A higher ∆D1 a-b indicates that the island area is increasing, the disturbance is decreasing, and the influence of this disturbance is getting more positive across years a and b.
(2) Soil salinization As a common ecological problem in coastal lowlands, the soil salinization is controlled by different factors. Specifically, seawater intrusion provides the salt source, and high evaporation rate, frequent droughts, irregular precipitations, and sea level changes accelerate the process [108][109][110][111]. Meanwhile, human activities have played an important role in influencing the soil salinization. They aggravate it by saline water irrigation and groundwater use and mitigate it by ecological measures and reasonable tillage practice [15,16,112]. The soil salinization was judged by the salinity content, and the disturbance of soil salinization across a time interval was measured using the following equation: where D2 is the disturbance of soil salinization; SS b and SS a are standardized salinity contents in an assessment unit in years b and a, respectively. The standardization method for the salinity is the same as for the aforementioned negative factors, and a higher ∆D2 a-b indicates that the degree of soil salinization is decreasing, the disturbance is decreasing, and the influence of this disturbance is getting more positive across years a and b.
(3) Human influence In recent decades, human influence has been the main driving factor of the ecological evolution on many islands [5,69,113,114]. Human activities not only participate in processes of the aforementioned disturbances, but also generate considerable influences on the sandy island ecosystem in different aspects. The classification and quantification of the human influence are complex and difficult [115,116]. A composite index, namely, human net influence index (HNII), was proposed by the authors [11] to identify the integrative influence of different types of human activities on the sandy island ecosystem. The HNII was calculated based on remote sensing by measuring the negative and positive influences of human activities on five aspects of the sandy island ecosystem, including island geomorphology, landscape pattern, plant community, physical quality, and chemical environment. The HNII has been applied in Chongming Island and proven to be effective in representing the human integrative influence on the sandy island ecosystem, and a higher HNII indicates a more positive influence. The disturbance of human influence across a time interval was measured using the following equation: where D3 is the disturbance of human influence; HNII b and HNII a are HNIIs in an assessment unit in years b and a, respectively. A higher ∆D3 a-b indicates that the disturbance is decreasing, and the human influence is getting more positive across years a and b.
The spatial distributions of variations of the three disturbances across different years are presented in Figure 5.

Resilience of Different Components under the Multiple Disturbances
The intrinsic correlation between components and disturbances at a single time point has been studied on Dongtou Archipelago [10]. In that study, the partial correlation analysis was employed to identify the factual influence of the disturbance on the component and quantify the influence degree based on the correlation coefficient; the results revealed the intrinsic correlations between components and disturbances [10]. However, the practice in the previous study was in a static state, i.e., at a single time point. The intrinsic correlations between the variations of components and disturbances of the sandy island ecosystem across different time points remain to be clarified in the present study for assessing the resilience. The resilience refers to the responses of different components to the three disturbances across a time interval. The variations of components (∆C) and disturbances (∆D) occur across the same time interval, that is, 1988-1995, 1995-2007, and 2007-2017.
A positive ∆C indicates an increase in the state value of a component across the time interval, and the higher ∆C, the higher the increase in the state value. A positive ∆D denotes a decrease in the degree of disturbance across the time interval, and the higher ∆D, the higher the decrease in the degree of disturbance. Meanwhile, ∆C corresponds to ∆D in each assessment unit. Thus, the components and disturbances are comparable after quantification and standardization, the intrinsic correlations between components and disturbances can be measured, and the resilience can be then obtained based on the intrinsic correlations. Two circumstances indicate the high resilience across a time interval: (1) when the ∆D is positive, the increase in the states of components is higher than the decrease in the degree of disturbances; (2) when the ∆D is negative, the decrease in the state of components is lower than the increase in the degree of disturbances. To be brief, higher values of (∆C − ∆D) denote higher resilience.
Furthermore, the intrinsic correlations between components and disturbances are complicated, that is, different components respond variously to the disturbance. Moreover, correlations may also differ among different disturbances. Thus, the partial correlation of each of the ∆Cs with each of the ∆Ds was analyzed through IBM SPSS 18 to explore the factual response of each component to each of the disturbances after eliminating the joint effects of the two other disturbances across different time intervals (Table 2). The positive partial correlations were considered the factual responses, and the negative ones were ruled out. Then, the resilience of each component under the multiple disturbances across different time intervals was calculated using the following equations: where RC1, RC2, and RC3 are the resilience of vegetation, soil, and landscape components, respectively, across a time interval; ∆C1, ∆C2, and ∆C3 are the variations of vegetation, soil, and landscape components, respectively, across a time interval; ∆D i is the variation of disturbance i across a time interval; and PCC i is the partial correlation coefficient of ∆C with ∆D i across a time interval, and only positive correlations are considered. PCC indicates the factual influence degree of a disturbance on a component, and (∆C − ∆D) denotes the change of the component under the disturbance. Then, the resilience of all the three components (RC) is calculated using the following equation: where RC i is the resilience of component i across a time interval.

Calculation of SIEHI-2
The SIEHI-2 should be measured for each year (1988, 1995, 2007, and 2017) to ensure the consistency with SIEHI-1 and to conduct the calculation of the following SIEHI. The resilience in one year is the average of the resilience in the time interval before the year and that after the year. Thus, the SIEHI-2 was calculated using the following equation: where RC before and RC after are the RCs in the time interval before one year and that after the year, respectively. Specifically, SIEHI-2 in 1995 is calculated using the average value of The SIEHI is composed of two parts, that is, SIEHI-1 and SIEHI-2. The two parts are closely inter-related and interacted, because the resilience is obtained based on the states of different components and their responses to the multiple disturbances. Furthermore, each of the two parts is essential to the sandy island ecosystem health, and the deterioration of each of the two parts will result in the distinct decrease in SIEHI. At the extreme, if each of SIEHI-1 and SIEHI-2 = 0, then SIEHI = 0. Thus, the SIEHI was calculated using the following equation:

Analyses of SIEHI-1, SIEHI-2, and SIEHI in Chongming Island
The maps of SIEHI-1, SIEHI-2, and SIEHI were generated through ArcGIS 10.0. Based on the results of the assessment units, the temporal changes of SIEHI-1, SIEHI-2, and SIEHI in the entire study area, in different land cover types, and for different protected areas from 1988 to 2017, were analyzed. Each of SIEHI-1, SIEHI-2, and SIEHI was divided into five levels, that is, Levels 1-5, when the index values are 0-0.2, 0.2-0.4, 0.4-0.6, 0.6-0.8, and 0.8-1, respectively, and area proportions of different levels were determined from 1988 to 2017. To verify the gradient effect of the distance to the shoreline (DTS) on the sandy island ecosystem health, changes of SIEHI-1, SIEHI-2, and SIEHI along the DTS gradient were analyzed. The spatiotemporal characteristics of the three components in the four years are presented in Figure 6. The state of vegetation component was better in the alongshore wetland vegetation than in the inner island, and in woodland and cropland than in built-up areas, water areas, and ponds. From 1988 to 2017, the state was becoming better over the entire island, which was largely manifested in the improvement of vegetation conditions in cropland and woodland in the inner island. In the alongshore areas, the state of the newly formed wetland vegetation was getting better, whereas that of the newly reclaimed ponds was getting worse. In terms of the soil component, most areas of the island were in the intermediate or bad state in 1988, and only a small part in the alongshore areas exhibited a good state. In 1995, the areas in an intermediate state increased compared with those in 1988. In 2007, the areas in a good state slightly increased compared with those in 1995. The soil component was generally in a good state in 2017, and the state over the island was distinctly better than that in the three other years. From 1988 to 2017, the state was generally becoming better. For the landscape component, the state was very good in the alongshore areas. By contrast, the areas in the inner island, particularly those along the roads or waterways and in the urbanized areas, were in a bad state. The state was generally becoming worse from 1988 to 2017.

Maps of SIEHI-1 in the Four Years
The maps of SIEHI-1 in the four years are presented in Figure 7. In 1988, the SIEHI-1 was high in the alongshore vegetation areas and generally low over the inner areas of the island. In 1995, the SIEHI-1 slightly decreased in the alongshore areas and the expanded built-up areas, while it increased in cropland and woodland compared with that in 1988. In 2007, the SIEHI-1 generally increased compared with that in 1995. Specifically, the SIEHI-1 in the northern and western alongshore areas distinctly increased, and that in the surrounding areas of built-up areas in the inner island slightly increased. In 2017, the SIEHI-1 increased over the island, except for in the eastern alongshore areas.

Assessment Results of Part 2 3.2.1. Resilience under Multiple Disturbances across Different Years
The spatiotemporal characteristics of the resilience of components under multiple disturbances are shown in Figure 8. The resilience of vegetation component was generally low across 1988-1995, and only a small part in the east of the island exhibited high resilience. The resilience was much higher across 1995-2007 than that across 1988-1995, and high resilience was observed over nearly the entire island. Across 2007-2017, this resilience showed clear spatial heterogeneity, that is, the assessment units with different values of the resilience were mixed. Furthermore, the western part of the island showed higher resilience than the eastern part. Overall, the three time intervals were in the following descending order of the resilience of vegetation component: 1995-2007, 2007-2017, and 1988-1995. In terms of the soil component, the resilience was low in the eastern alongshore areas and the northwestern part in the inner island across 1988-1995. Across 1995-2007, the western, northern, and eastern alongshore areas exhibited low resilience, yet the resilience in inner areas of the island was generally higher than that across 1988-1995. Across 2007-2017, the resilience in the alongshore areas was still low, while the inner island generally showed higher resilience compared with that across 1995-2007. Overall, the resilience in the alongshore areas remained low, yet that in the inner island increased along the three time intervals. The resilience of landscape component stayed low in the alongshore areas across the three time intervals. Meanwhile, the resilience in the inner island was distinctly higher across 1995-2007 than that across 1988-1995 and 2007-2017. For the resilience of all the three components, it combined the spatiotemporal characteristics of vegetation, soil, and landscape components. Specifically, the alongshore areas exhibited low resilience across the three time intervals, which was in accordance with the characteristics of the resilience of soil and landscape components; the resilience in the inner areas of the island was generally low and high across 1998-1995 and across 1995-2007, respectively, and showed clear spatial heterogeneity across 2007-2017, which was in accordance with the characteristics of the resilience of vegetation component.

Spatiotemporal Characteristics of the Sandy Island Ecosystem Health
The overall temporal changes of the sandy island ecosystem health from 1988 to 2017 are shown in Figure 11. The SIEHI-1 slightly decreased from 1988 to 1995 and then distinctly increased from 1995 to 2017. The SIEHI-2 initially increased from 1998 to 2007 and then decreased from 2007 to 2017. The SIEHI continuously increased from 1988 to 2017.
The SIEHI-1, SIEHI-2, and SIEHI in different land cover types and for different protected areas are presented in Table 3. In different land cover types, SIEHI-1 was high in wetland vegetation, intermediate in cropland and woodland, and low in the remaining land cover types. From 1988 to 2017, SIEHI-1 in wetland vegetation, mudflat, and traffic land initially decreased and then increased; that in woodland showed fluctuating characteristics; that in cropland, water areas, building land, and industrial land continuously increased; and that in ponds did not present distinct changes. For the SIEHI-2 from 1988 to 2017, it remained low in wetland vegetation, showed decreasing characteristics in mudflat, and initially distinctly increased and then slightly decreased in the remaining land cover types. The SIEHI showed the similar characteristics to the SIEHI-1 for different land cover types. In terms of the protected areas, SIEHI-1 was higher in core and ordinary protected areas than in nonprotected areas in all the four years. In 1988, SIEHI-1 was much higher in core protected areas than in ordinary protected areas. In 1995, 2007, and 2017, ordinary protected areas achieved the highest SIEHI-1. The SIEHI-2 generally showed decreasing and increasing characteristics in core and ordinary protected areas, respectively, from 1988 to 2017. It first increased from 1988 to 2007 and then slightly decreased from 2007 to 2017 in nonprotected areas. The SIEHI was highest in core protected areas in 1988, in ordinary protected areas in 1995 and 2007, and in nonprotected areas in 2017.  The proportions of different levels of the sandy island ecosystem health are presented in Table 4.  The changes of the sandy island ecosystem health along the DTS gradient are shown in Figure 12. In 1988, SIEHI-1, SIEHI-2, and SIEHI all showed decreasing characteristics along the gradient, and the characteristics for SIEHI-1 and SIEHI were more distinct than those for SIEHI-2. The lines were gentle when DTS > 3 km. Along the gradient in 1995, the SIEHI-1 decreased when DTS < 3 km and then gently fluctuated, the SIEHI-2 increased when DTS < 5 km and then gently fluctuated, and the SIEHI showed gently fluctuating characteristics. In 2007, SIEHI-1 and SIEHI decreased when DTS < 3 km and then gently fluctuated, and the SIEHI-2 slightly increased along the gradient. In 2017, the SIEHI-1 distinctly decreased when DTS < 3 km and then slightly increased, the SIEHI-2 continuously increased and the increasing characteristics were more distinct than those in 2007, and the SIEHI initially increased when DTS < 1 km, then decreased when 1 km < DTS < 4 km, and finally increased when DTS > 4 km.

. Two Parts of the Sandy Island Ecosystem Health
Two parts of the sandy island ecosystem health were identified, analyzed, and assessed in this study.
Part 1 was assessed based on the recognized framework of the three components, namely, vegetation, soil, and landscape. In previous studies, this framework has been successfully applied in Dongtou Archipelago in southern China and Miaodao Archipelago in northern China, which are constituted mainly by rocky islands [10,41]. This study demonstrated the effectiveness of the framework in assessing the states of essential components in the sandy island ecosystem. Rocky and sandy islands are two major types of islands and occupied most of the world's islands and nearly all the islands in middle and high latitudes [5,92,117]. Previous and present studies proved the good applicability of the framework in these two types of islands, and its applicability in coral reef islands remains to be studied in future work. Part 2 was assessed by identifying the multiple disturbances, tracking the rapid temporal variations, and clarifying the intrinsic correlations between the components and disturbances. Various natural and anthropogenic disturbances, including sea level rise [47], saltwater intrusion [48], sea reclamation [49], onshore industries [50], and aquaculture [51], existed and set a special and complex context of the coastal ecosystem. In the present study, three main disturbances were identified by considering the above potential disturbances that generate distinct influences. All of the three disturbances were measured using quantifiable indicators, parts of which were sourced from previous reliable results [11]. The temporal variations were analyzed on the premise of the comparability of the components and disturbances, and the three time intervals were exhibited based on the existing data in the 4 years during the last 30 years. The intrinsic correlations could be considered causal links between the disturbances (causes) and components (effects) after the aforementioned identification of multiple disturbances. A partial correlation analysis was used to clarify the intrinsic correlations. This analysis could reveal the practical effects generated by one type of disturbance by eliminating the joint effects from the two other disturbances [10]. In this study, nearly all of the partial correlation coefficients were positive, indicating the distinct causal links between the disturbances and components. The partial correlation coefficients greatly varied in different disturbances and components, thereby revealing the different influences from different disturbances on different components of the sandy island ecosystem. Based on the aforementioned work, the resilience was measured by connecting the components with the disturbances using their comparable values and intrinsic correlations in different time intervals. Then, the resilience in different years was obtained by resampling the results across the three time intervals to those in the four years.
The two parts can be measured and integrated to comprehensively represent sandy island ecosystem health. Part 1 contributed to the wide application of the framework of vegetation, soil, and landscape in different types of islands. Part 2 provided a new method for measuring the resilience of the sandy island ecosystem responding to multiple disturbances across different time points, and the applicability of the method also remains to be validated on different islands in future studies. The method proposed in the present study corresponds to the three points existing in the previous studies, which were mentioned in the Introduction section. Specifically, the method comprehensively considers the essential components of the island ecosystem and the main external disturbances by integrating the field and remote sensing data, clarifies the intrinsic correlations between the variations of components and disturbances across different time points, and realizes the spatial visualization of the state and resilience in the assessment results, thereby effectively solving the three points and mapping the spatiotemporal pattern of sandy island ecosystem health during the last decades.

Adequate Utilization of Remote Sensing Data in the Assessment
The other key issue for realizing the assessment was the adequate utilization of remote sensing data. Remote sensing as a rapidly developing and highly applicable technique has been widely applied in ecosystem health research for the abundant spatiotemporal information and accessible data source [4,27,52]. In the present study, the spatiotemporal information of remote sensing data was adequately utilized in two aspects.
First, the remote sensing data were directly used for the extractions of island outlines and land cover types, as well as the calculations of ecological indexes. The island outlines set the extent of the study area in different years and were used to analyze the geomorphological change. The land cover types served as the basic data for calculating VA in the vegetation component and all the factors (ILC, NP, TE, and LII) in the landscape component. Two of the ecological indexes, namely, NDVI and NPP, were factors in the vegetation component. Moreover, HNII as the indicator of human interference was calculated by combining the land cover and ecological indexes [11]. Therefore, seven factors in the components, including VA, NDVI, NPP, ILC, NP, TE, and LII, and two disturbances were all obtained directly using the remote sensing data. Their results were area data that covered the entire study area across different years during the last decades, thereby meeting the requirement of assessment of the sandy island ecosystem health.
Second, the remote sensing data functioned as the predictors to conduct the spatial simulations of field data from point to area and from present to past. Some factors in the components and disturbances, including H', BD, TOC, TN, AK, and Sa, were sourced by field work. Unlike the above factors that were directly sourced by remote sensing and covered the entire island in different years, these field data were point data only at a time point. In previous studies, the spatial simulations of the field data from point to area and from present to past have been implemented, and good effects have been achieved, thereby proving the feasibility to conduct such simulations [72,97,118]. In this study, we conducted the spatial simulations based on the multitemporal remote sensing data (Landsat series satellites), various predictors that adequately utilized the ecological information of the remote sensing data (spectral reflectance and reflectance-based ecological indexes), and reliable algorithm that had high maneuverability and repeatability (PLSR). To avoid the potential errors generated by the temporal variations of the remote sensing data, controlling measures such as radiometric calibration and atmospheric correction were conducted, and the changing bands and indexes such as thermal infrared bands and brightness temperatures were not selected as the predictors. The simulation results showed acceptable accuracies for different simulated factors. Thus, the factors that ranged from 1988 to 2017 and covered the entire study areas were obtained, thereby meeting the requirement of the assessment. However, two points remain to be promoted, that is, the collection of measured data to validate the simulation results in different years and the improvement of the simulation accuracy using more sophisticated procedures.
Therefore, remote sensing provided the major data source of the present study, and the utilization of remote sensing data in the two aspects adequately tapped their ecological significance and effectively supported the assessments of the two parts of the sandy island ecosystem health.

Contributions of Components and Factors to the Spatiotemporal Variations
To identify the main driving factors of spatiotemporal variations of the ecosystem health of Chongming Island, the contributions of components and factors to the spatiotemporal variations were measured using correlation analyses. The correlation coefficients of the sandy island ecosystem health with different components and factors in the four years were presented in Table 5. All of the correlations were significant because of the enormous analyzed samples, yet the correlation coefficients varied considerably. Among SIEHI, SIEHI-1, and SIEHI-2, the correlation coefficient between SIEHI and SIEHI-1 was constantly higher than that between SIEHI and SIEHI-2 in all the four years; the difference was small in 1988 and 2017 and large in 1995 and 2007, which indicated that the spatial characteristics of SIEHI were influenced by both of SIEHI-1 and SIEHI-2 in 1988 and 2017, while mainly controlled by SIEHI-1 in 1995 and 2007. The correlation coefficient between SIEHI-1 and SIEHI-2 was in a very low level in the four years, exhibiting their weak spatial correlation. It indicated that the two parts of the sandy island ecosystem health were relatively independent to each other in the temporal scope of the study, and each of them was indispensable. That is, the high value of only one of the two parts did not necessarily denote good sandy island ecosystem health. In terms of the three components and 12 factors, SIEHI and SIEHI-1 showed certain high correlation coefficients with them, whereas SIEHI-2 exhibited weak correlations, which was in accordance with the aforementioned correlation between SIEHI-1 and SIEHI-2. Generally, correlation coefficients for vegetation and soil components were higher than those for the landscape component in the four years, of which the vegetation component achieved the highest one, and those for landscape component were particularly lower in 2007 and 2017. The vegetation and soil components, which mainly refer to the microscale state, contributed much more to the spatial variance of sandy island ecosystem health than the landscape component, which mainly represents the macroscale state [5,119]. This result was not in accordance with the previous studies on the rocky islands, which showed that the soil and landscape components were the major components on Dongtou Archipelago and Miaodao Archipelago, respectively [10,41]. Though the three components are all essential to the island ecosystem, their importance degrees may change across different islands. In the vegetation component, NDVI and NPP showed high correlation coefficients in all the four years, indicating the large contribution of vegetation growth condition to the sandy island ecosystem health, which was in accordance with the previous studies on the rocky islands [10,41]. The correlation coefficient for VA was high in 1988, 1995, and 2007 and distinctly low in 2017. By contrast, that for H' continuously increased from 1988 with a low value to 2017, when a high one was achieved. It indicated that the plant diversity has played an increasingly important role, and the vegetation quality rather than the quantity greatly contributed to sandy island ecosystem health in recent years. For the four factors in the soil component, TOC and TN exhibited high and increasing correlation coefficients from 1988 to 2017, whereas correlation coefficients for the remaining two factors were low. TOC and TN, which refer to the two basic elements in the soil [90], made large contributions to the sandy island ecosystem. In the landscape component, the correlation coefficient for ILC was high only in 1988 and low in 1995, 2007, and 2017. The remaining three factors showed intermediate correlation coefficients in all the four years, which indicated that landscape fragmentation was not the decisive factor influencing the sandy island ecosystem health and explained the relatively lower contribution of landscape component than the two other components.

Contributions of Disturbances and Resilience to the Spatiotemporal Variations
The influences of different disturbances on the states of components have been shown in Table 2. The partial correlation coefficient between ∆C and ∆D denoted the independent influence of the variation of one disturbance on that of one component. The ∆D1 showed positive correlations with ∆C1 and ∆C3 and a negative correlation with ∆C2 in all the three time intervals. The increase in the island area brought about the improvement in states of vegetation and landscape components because the newly formed areas were always covered by wetland vegetation, which showed good growth condition and high landscape connectivity. Meanwhile, the soil in these newly formed areas was severely influenced by seawater and thus had a low state. Distinct positive correlation was observed between ∆D2 and ∆C1, indicating that the decrease in the degree of soil salinization remarkably helped improve the state of the vegetation component, particularly in the inner areas of the island covered by crops and forests. This result was in accordance with a previous study on coastal forested wetlands [43]. Furthermore, ∆D2 was positively and negatively correlated with ∆C2 and ∆C3, respectively, though in a limited degree, showing the influences of soil salinization on the soil and landscape components. As shown in Figure 5, areas with positive ∆D2 covered most of the island across 1988-1995 and 2007-2017, while those with negative ∆D2 occupied the most across 1995-2007. The degree of soil salinization showed a fluctuating temporal change under the natural and anthropogenic factors and was generally getting better after the long-term and island-wide tillage practice [71]. High partial correlation coefficients of ∆D3 with ∆C1, ∆C2, and ∆C3 revealed the distinct human influences on the three components. The results denoted the high sensitivity of all the three components to human influences. That is, positive human influences contributed to increasing the states of the three components, while negative human influences resulted in the decreases in the states.
The correlation coefficients among the resilience of different components are presented in Table 6. Across the three time intervals, RC2, RC1, and RC3 were in the descending order of the correlation coefficients with RC, indicating their different contributions to the spatiotemporal variations of RC. Among RC1, RC2, and RC3, only the correlation between RC2 and RC3 was distinct, revealing the intrinsic link between the resilience of soil and landscape components. In addition, RC2 also showed a certain positive correlation with RC1. The soil served as the base of the ecosystem, and high resilience of the soil component may improve the resilience of the two other components [10,120].

Driving Factors of the Spatiotemporal Variations
Based on the aforementioned analyses, the spatiotemporal variations of sandy island ecosystem health were controlled by the two parts, namely, the state and resilience. In different components, the vegetation and soil components contributed more than the landscape component to the sandy island ecosystem health. In the specific factors, vegetation quality and soil carbon/nitrogen were the ones that showed the most contributions.
Of the three disturbances, the human influence on the three components was generally higher than the two other disturbances, indicating the major role of human activities in driving the spatiotemporal variations of the sandy island ecosystem health. This result was in accordance with the previous studies on rocky islands [41,113,114] and coastal wetlands [4,52]. The states of the three components were quickly getting better when the human influence was getting more positive, and vice versa. However, the influences generated by human activities were not overwhelmingly dominant. The geomorphological change and soil salinization also exerted practical influences on different components to certain degrees. The vegetation component received an increasing influence from geomorphological change yet a decreasing influence from soil salinization from 1988 to 2017, and the increase in the island area and the decrease in the degree of soil salinization brought about the improvement of the vegetation component. In addition, the increase in the island area resulted in the deterioration and improvement of the soil and landscape components, respectively. These results confirmed that the spatiotemporal variations of the sandy island ecosystem health were controlled by multiple natural and anthropogenic disturbances, and human influence played a major role during the last decades.

Influence Coefficients of Different Human Activities on the Sandy Island Ecosystem
The human activity has been proven to be the major factor influencing sandy island ecosystem health. Accurately measuring the influences of different human activities on the island ecosystem is the premise to balance the island protection and development and achieve island ecosystem-based management. The influence coefficients of different human activities can provide a universal standard for measuring the anthropogenic influences through a simple process. Specifically, the single influence from one type of human activity and the comprehensive influences from multiple types of human activities can be easily quantified based on the areas and the influence coefficients of different human activities. In most of the previous studies, the influence coefficients were empirically assigned to ensure the simplicity, yet the accuracy was uncertain [92,115,121]. To meet the demands of high simplicity and accuracy, a previous study on the ecosystem health of rocky islands has provided the influence coefficients of typical human activities on the rocky island [41]. In this study, the influence coefficients for the sandy island could also be obtained based on the results. Seven typical types of human activities were extracted based on the land cover types, that is, plantation, farming, waterway construction, pond aquaculture, building construction, traffic development, and factory construction. The wetland vegetation and mudflat, which received little or low anthropogenic influence, were considered natural areas of the sandy island, and the area-weighted SIEHI of wetland vegetation and mudflat was used as the baseline to measure the changes of SIEHI under different types of human activities. The relative SIEHI of one type of human activity was obtained using the ratio of the SIEHI of the human activity to the aforementioned area-weighted SIEHI. Then, the influence coefficients of different human activities were calculated using the following equation [41]: where IC is the influence coefficient of one type of human activity; rSIEHI is the relative SIEHI of one type of human activity; T is the threshold. Two schemes exist for the threshold selection. Scheme 1: A threshold of 0.7 was used based on the island carrying capacity [122], and the influence coefficients are comparable to those for the rocky islands. Scheme 2: The lowest relative SIEHI of different types of human activities (i.e., 0.4384 for the traffic development) was used as the threshold, and the influence coefficients are applicable for only the sandy island and can better represent the differences among different human activities. The calculated influence coefficients of different types of human activities on the sandy island ecosystem are shown in Table 7. In terms of Scheme 2, traffic development, waterway construction, and factory construction were the three types that exerted the highest influences per area. In consideration of the area of each type, farming and building construction were the two types that showed the highest influences, and accounted for 37.12% and 35.59%, respectively, of the anthropogenic influence.

Measures for the Sandy Island Ecosystem-Based Management
Chongming Island is the hotspot of ecosystem studies, and various measures for conserving the island ecosystem have been proposed by scholars and administrators. As a time-honored site for human inhabitation and agricultural activities and an indispensible part of Shanghai City, the largest city in China, the economic development is one of the main goals of the island [77]. Meanwhile, the ecological functions of the island received increasing attention in recent decades and were unprecedented valued in the context of the China's theory of "Lucid Water and Lush Mountains are Invaluable Assets" [123]. The development and protection are both essential to the island and should be well-balanced to pursue the sandy island ecosystem-based management. Different functional zones can be divided based on the results of the state and resilience of the sandy island ecosystem health for balancing the development and protection through conducting different measures in different zones. According to the five levels for SIEHI-1 and SIEHI-2 (Table 4), three functional zones, namely, Zone A, Zone B, and Zone C, were assigned for different purposes in 2017. Along the zones from Zone A to Zone C, the protection and development functions were decreasing and increasing, respectively. The division scheme for different zones is as follows: Zone A: areas in Levels 4 or 5 of SIEHI-1 and in Levels 1, 2, or 3 of SIEHI-2, as well as areas in Level 3 of SIEHI-1 and in Levels 1 or 2 of SIEHI-2; Zone C: areas in Levels 1 or 2 of SIEHI-1 and in Levels 3, 4, or 5 of SIEHI-2, as well as areas in Level 3 of SIEHI-1 and in Levels 4 or 5 of SIEHI-2; Zone B: the remaining areas. The scheme was proposed by referencing the previous study [10]. Thus, the three functional zones in 2017 were obtained and are shown in Figure 13. The area proportions of Zones A, B, and C were 24.01%, 36.88%, and 39.10%, respectively. Zone A was mainly along the shoreline with continuous distribution and also in areas with good state and bad resilience in the inner island with scattered distribution, and it covered most of the protected areas of the island. Zone C was in different types of development areas, including urban constructions, pond aquaculture, waterways, and parts of croplands. Zone B was in the intermediate areas between Zones A and C. The areas and spatial distributions of different zones ensured the core functions of protection (Zone A) and development (Zone C), as well as the multiple functions (Zone B) on the basis of factual conditions. Different measures should be formulated for different zones. In Zone A, strict protection measures should be implemented. In this zone, the SIEHI-1 was generally high, while the SIEHI-2 was generally low. For the three components, the damaged vegetation and soil should be ecologically restored; the current vegetation should be protected from human occupation and biological invasion; and the landscape connectivity should be improved by removing the obstacles and constructing the ecological corridors. In terms of the three disturbances, most of the areas in Zone A were in the alongshore areas and were susceptible to geomorphological change and soil salinization. In areas in which coastal erosions often occurred, revetment measures with the ecological revetment as the primary means should be conducted. Specific measures to control the soil salinization are not needed in these areas because the wetland vegetation species have certain salt-tolerance. The human influences also need attention, and human activities that are not for protection purposes should be prohibited. In Zone B, the harmonious coexistence of protection and development is pursued. This zone was mainly composed of cropland, woodland, and parts of wetland vegetation and built-up areas. Within this zone, areas with different levels of SIEHI-1 or SIEHI-2 could be observed. The areas with high SIEHI-1 had high SIEHI-2, and vice versa. Soil salinization and human influence were two major disturbances to this zone. The former may severely threaten the vegetation growth in cropland and woodland. Different measures for controlling the soil salinization, including physical, chemical, and biological modifications, should be conducted to mitigate or even eliminate the negative influence from soil salinization [112,124]. For the human influence, the spatial patterns of human activities should be optimized based on the SIEHI-1 and SIEHI-2. In areas with high SIEHI-1 and SIEHI-2 (e.g., parts of wetland vegetation), moderate human activities can be conducted in forms of ecological tourism. In areas with intermediate SIEHI-1 and SIEHI-2 (e.g., cropland and woodland), high-quality agricultural activities should be conducted through rational tillage practice. In areas with low SIEHI-1 and SIEHI-2 (e.g., parts of built-up areas), the constructions that have low density and intensity are encouraged and ecological restoration measures should be simultaneously carried out to improve SIEHI-1 and SIEHI-2. In Zone C where the SIEHI-1 was generally low while the SIEHI-2 was generally high, measures for reasonable and efficient development should be conducted for maximizing the economic and social benefits and minimizing the negative effects on the sandy island ecosystem. The urban construction should be organized and compact with high intensities for reserving space for future generations [125,126]. Meanwhile, the greenway network should be continuously constructed to improve the ecological connectivity of built-up areas, which aims to mitigate the influences from landscape fragmentation on the island biodiversity [11]. From the perspective of a larger scale, the pressure from the surrounding and relevant areas, such as the decline in water and sediment discharge along the Yangtze River [13], the high threat degree of wastewater pollutants in the Yangtze River Delta [50], and the tourist visits from the urban agglomeration [127], should also be noticed and the corresponding measures should be made in the way of regional linkage and cooperation.

Conclusions
This study assessed the sandy island ecosystem health by integrating the two parts, that is, the states of essential components and the resilience under multiple disturbances. The states were measured following the framework of vegetation, soil, and landscape, and a total of 12 factors in the three components were adopted. The results demonstrated the effectiveness of the framework in the sandy island ecosystem, which corresponded to the first scientific question proposed at the end of the Introduction section. Three typical disturbances on the sandy island, namely, geomorphological change, soil salinization, and human influence, were identified, and the resilience was measured by clarifying the intrinsic correlations between the variations of components and disturbances across different time intervals, thereby solving the second scientific question of the study. Then, the contributions of components, factors, and disturbances to the spatiotemporal variations of sandy island ecosystem health were analyzed, and the main driving factors were identified, which corresponded to the third scientific question. Based on the assessment results, influence coefficients of different human activities were determined for providing a universal standard and a simple method for measuring the anthropogenic influences on the sandy island ecosystem, and different measures for the sandy island ecosystem-based management were proposed for different functional zones. The novel suggestions answered the fourth scientific question.
The results in Chongming Island indicated that the SIEHI-1 in the alongshore areas was generally high in the four years, and that in the inner island was low in 1988 yet showed increasing characteristics from 1988 to 2017. By contrast, the SIEHI-2 was generally low in the alongshore areas in the four years, while that in the inner island was low, intermediate, and high in 1988, 1995, and 2007, respectively, and showed distinct heterogeneity in 2017. The SIEHI combined the spatiotemporal characteristics of SIEHI-1 and SIEHI-2. In the entire study area from 1988 to 2017, the SIEHI-1 initially slightly decreased and then distinctly increased, the SIEHI-2 first distinctly increased and then slightly decreased, and the SIEHI continuously increased. In different land cover types, the SIEHI and SIEHI-1 were high in wetland vegetation, intermediate in cropland and woodland, and low in the remaining land cover types in the four years. The SIEHI-2 remained low in wetland vegetation, showed decreasing characteristics in mudflat, and initially distinctly increased and then slightly decreased in the remaining land cover types from 1988 to 2017. Along the DTS gradient, the SIEHI showed decreasing and gently fluctuating characteristics in 1988 and 1995, respectively; it first decreased and then gently fluctuated in 2007; and it initially increased, then decreased, and finally increased in 2017. In different components and factors, the vegetation and soil components contributed more than the landscape component to the sandy island ecosystem health, and vegetation quality and soil carbon/nitrogen were the factors that showed the most contributions; in different disturbances, the human influence played a major role in driving the spatiotemporal variations of the sandy island ecosystem health, and the geomorphological change and soil salinization also exerted practical influences on different components to certain degrees.

Data Availability Statement:
The data presented in this study are available upon reasonable request.