Next Article in Journal
The DPSIR Model for Environmental Risk Assessment of Municipal Solid Waste in Dar es Salaam City, Tanzania
Next Article in Special Issue
Not-in-My-Backyard: Legislation Requirements and Economic Analysis for Developing Underground Wastewater Treatment Plant in China
Previous Article in Journal
Human Papillomavirus-16 DNA Quantitation Differentiates High-Grade Anal Neoplasia
Previous Article in Special Issue
Changes of Ecosystem Service Value in a Coastal Zone of Zhejiang Province, China, during Rapid Urbanization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution of Landscape Ecological Risk at the Optimal Scale: A Case Study of the Open Coastal Wetlands in Jiangsu, China

1
Department of Geographic Information Science, Nanjing University, Nanjing 210023, China
2
Key Laboratory of Coastal Zone Exploitation and Protection, Ministry of Land and Resource, Nanjing 210024, China
3
Collaborative Innovation Center of South China Sea Studies, Nanjing University, Nanjing 210093, China
4
Jiangsu Provincial Key Laboratory of Geographic Information Science and Technology, Nanjing University, Nanjing 210023, China
5
Department of Geography & Spatial Information Techniques, Ningbo University, Ningbo 315211, China
6
Institute for East China Sea, Ningbo University, Ningbo 315211, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(8), 1691; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081691
Submission received: 11 June 2018 / Revised: 19 July 2018 / Accepted: 27 July 2018 / Published: 8 August 2018

Abstract

:
Detailed analysis of the evolution characteristics of landscape ecological risk is crucial for coastal sustainable management and for understanding the potential environmental impacts of a man-made landform landscapes (MMLL). As a typical open coastal wetland, large-scale human activities (e.g., tidal reclamation, fishery activities, wind farm construction, and port construction) have substantially affected the evolution of the coastal ecological environment. Previous landscape ecological risk assessment studies have documented the effectiveness of assessing the quality of ecological environment processes. However, these studies have either focused on the noncoastal zone, or they have not considered the evolution of the spatial characteristics and ecological risk evolution of the landscape at an optimal scale. Here, we present a landscape ecological risk pattern (LERP) evolution model, based on two successive steps: first, we constructed an optimal scale method with an appropriate extent and grain using multi–temporal Landsat TM/OLI images acquired in the years 2000, 2004, 2008, 2013 and 2017, and then we calculated landscape ecological risk indices. Based on this model, the entire process of the spatiotemporal evolution of ecological risk patterns of the open coastal wetlands in Jiangsu, China, was determined. The principal findings are as follows: (1) The main landscape types in the study area are tidal flats and farmland, and the main features of the landscape evolution are a significant increase in aquafarming and a substantial decrease in the tidal flat area, while the landscape heterogeneity increased; (2) In the past 20 years, the areas of low and relatively low ecological risk in the study region were greatly reduced, while the areas of medium, relatively high, and high ecological risk greatly increased; the areas of high-grade ecological risk areas are mainly around Dongtai and Dafeng; (3) The area of ecological risk from low-grade to high-grade occupied 71.75% of the study area during 2000–2017. During the previous periods (2000–2004 and 2004–2008), the areas of low-grade ecological risk were transformed to areas of middle-grade ecological risk area, while during the later periods (2008–2013 and 2013–2017) there was a substantial increase in the proportion of areas of high-grade ecological risk. Our results complement the official database of coastal landscape planning, and provide important information for assessing the potential effects of MMLL processes on coastal environments.

1. Introduction

Landscapes have been considered as a system composed of functionally interacting land units [1], which can be characterized based on the diversity and abundance of land cover and landscape pattern metrics [2]. Landscapes are also locations of human activities within which the sustainable use of resources is an important research topic [3], and thus, they are often regarded as a suitable scale for research on the effects of human activities on the environment [4]. Landscape ecology emphasizes the interaction between spatial patterns, ecological processes and scale. Landscape patterns often affect ecological processes (e.g., population dynamics, animal migration, biodiversity and ecological circadian condition) [5,6]. We can obtain a better understanding of ecological processes through the study of spatial patterns, due to the interaction between patterns and processes, and therefore ecological risk analysis based on landscape patterns can provide an integrated assessment of various types of potential ecological impacts and their cumulative consequences [7]. In addition, in landscape ecology, the scale effect refers to the phenomenon in which the spatial heterogeneity of landscape space changes when there are changes in extent or grain [8]. Studies have demonstrated that the spatial resolution (i.e., grain) of the analysis data and the spatial extent of the study area affect landscape analysis results [9]. Hence, to facilitate environmental management and decision making, it is necessary to consider the evolution of the landscape over a sufficiently long time interval.
With increasing influence of human activities on the evolution of the Earth’s surface environment system (ESES), and the overloading of its natural functioning, ESES has entered a stage of physical and chemical evolution within an artificially-driven framework [10]. In particular, since the onset of the Anthropocene [11], numerous man-made landform landscapes (MMLL) have been created via the direct or indirect combination of human activities and natural forces in the context of natural geography [11,12]. Moreover, in many cases, instantaneous information about the factors influencing MMLL processes within a specific region is difficult to obtain, or may not even exist momentarily. Fortunately, satellite remote sensing data may be an excellent substitute [13,14], and advances in satellite remote sensing can provide an unparalleled record of Earth’s status and dynamics [15]. The ever-increasing amount of updated and archived images offers an opportunity to determine the evolution and to assess the impacts of human activities on the eco-environment over large areas.
Ecological risk refers to the probability and possibility of the negative impact under external pressures (e.g., from natural changes or human activities) on ecosystem structure, functioning, and even ecosystem stability and sustainability [16,17], which are reflected in ecosystem health, productivity, and genetic structure, and economic and aesthetic values [18]. Ecological risk assessment is a process that assesses the adverse environmental impacts that may occur or are occurring due to one or more external factors [19]. Landscape ecology emphasizes the interaction between spatial patterns, ecological processes and scale. Landscape patterns often affect ecological processes (e.g., population dynamics, animal migration, biodiversity and ecological circadian conditions) [5,6]. We can obtain a better understanding of ecological processes through the study of spatial patterns, due to the interaction between patterns and processes. Therefore, ecological risk analysis based on landscape patterns can provide an integrated assessment of various types of potential ecological impacts and their cumulative consequences [7]. In landscape transition risk assessment studies, the landscape index used to measure the landscape pattern of an area was often included within the landscape transition risk assessment framework [20]. Several studies have already adopted a more meaningful ecological risk research framework in ecological assessment [21,22]. Specifically, comprehensive knowledge of landscape pattern processes is increasingly important for the quantitative assessment and analysis of the spatial distribution of ecological risk and to determine the evolutionary process of coastal ecosystems [23]. Most previous landscape ecological risk studies have focused on coastal patterns of ecological evolution [24,25,26,27] using land-use data or a range of landscape metrics [27,28,29,30]. However, the use of a landscape transition risk assessment framework is rare in studies of landscape pattern processes.
Human activities are exerting increasing pressure on coastal ecosystems worldwide, which is also true for the open coastal wetlands in Jiangsu, China. Coastal regions are exposed to a wide variety of hazards, such as landscape changes, habitat loss, seawater intrusion, coastal flooding and waterlogging, which are exacerbated by climate change [31,32,33]. Moreover, many marine and coastal projects have been developed to meet human needs [34,35], and the demand for such developments in coastal regions is often met by exploiting natural land (e.g., tidal flats, water, and open spaces), which may considerably reduce the area of ecological land [27]. Importantly, China’s growing economy is increasingly concentrated in coastal regions [36]; thus, the outline of the Jiangsu coastal reclamation development plan (2010–2020) was passed in 2009 [27], and landscapes such as tidal flats, and areas of cordgrass and open seawater are to be used for agriculture and urban development.
Landscape patterns that anthropogenic behaviors interact with represent the spatial consequences of underlying biophysical or socioeconomic processes [37]. At the same time, the evolution of coastal ecological risk patterns is a complicated process and coastal ecosystems are increasingly negatively impacted. Moreover, ecological losses mean that ecosystems in coastal areas are becoming increasingly sensitive and vulnerable under external pressures [38,39]. The recognition of these issues has resulted in increasing research emphasis on coastal dynamics, human development and coastal relationships, reflected by projects such as LOICZ (land-ocean interactions in the coastal zone, International Geosphere-Biosphere Programme core research program) and the FEC (Future Earth-Coasts, a project to support sustainability and adaptation to global change in the coastal zone; https://www.future earthcoasts.org/) [40]. Consequently, in the context of disaster risk reduction, the assessment of landscape features and the evolution of ecological risks are increasingly important for sustainable coastal planning and management, and it also contributes to the growing number of case studies of coastal changes worldwide.
In light of the above, the specific objectives of the present study of the open coastal wetlands in Jiangsu Province (OCWJ) were: (1) to construct landscape ecological risk indices using geographic information systems (GIS), landscape ecological metrics, Landsat-5 Thematic Mapper (TM) and Landsat-8 Operational Land Imagery (OLI) data; (2) to develop evolution models of landscape ecological risk pattern (LERP); and (3) to identify the spatiotemporal variation of landscape ecological risk patterns during the period of 2000–2017.

2. Materials and Methods

2.1. Study Area

The open coastal wetlands in Jiangsu (OCWJ) lie within the area bounded by 32°34′ and 34°28′ N, and 119°48′ and 120°56′ E (Figure 1). The wetlands were formed by sediment deposition within the old Yellow River Delta and the Yangtze River Delta under the impact of waves and tides in the Yellow Sea and the East China Sea. The study area is affected by the alternation of Tropical Ocean and polar continental air masses governed by the East Asian monsoon. The coastal zone has a maritime monsoon climate ranging from warm-temperate to northern subtropical [41]. The OCWJ are an important pathway and stopover of migratory birds in Northeast Asia–Australia. The terrestrial area of the province encompasses 101,800 km2 and it is one of the smallest provinces in China (only 1.06% of the total area). However, the province possesses an extensive coastal area and includes seven coastal cities (e.g., Xiangshui, Binhai, Sheyang, Dafeng, and Dongtai) (Figure 1). In this study, the scope of the Yancheng Wetland Biological Reserve (i.e., the core zone, the buffer zone, and the transition zone, China, http://sts.mep.gov.cn/zrbhq/) was taken as the study boundary.

2.2. Data Sets and Processing

Geographic information systems (GIS) and remote sensing (RS), as a mature Earth observation (EO) techniques, are powerful and cost-effective tools for studying the spatiotemporal dynamics of ecological and geographical information [43,44,45]. Landscape data were obtained from Thematic Mapper(TM), and Operational Land Imager(OLI) satellite images in the years 2000, 2004, 2008, 2013, and 2017 (Table 1) from the United States Geological Survey (USGS, http://earthexplorer.usgs.gov/). A total of 10 images were acquired in autumn and winter (October to February) to reflect the true wintering environments of the OCWJ. All the image dataset projections are the Universal Transverse Mercator (UTM) map projection system zone 50N and datum of the World Geodetic System 1984 (WGS84), ensuring consistency between datasets during analysis. Landsat TM/OLI image processing was done using ENVI (version 5.2, the Environment for Visualizing Images, Exiles Visual Information Solutions Company, Virginia, VA, USA.), including image geometric correction, atmospheric correction, image mosaicking, and masking-out the study area. The starting point for this study, especially evident since 2000, was the initial period of unprecedented socioeconomic transformation and intensified human activity which have dramatically modified ecosystems and landscapes across the nation [46]. Auxiliary data consisted of the administrative map of Jiangsu province, a land use map and land use history data, the topographic map of Jiangsu Province (1:50,000), and GPS field survey data and Google Earth samples.

2.3. Methodology

A framework was developed for assessing the ecological risks for open coastal wetlands in Jiangsu Province using evolution models of landscape ecological risk pattern (LERP). Figure 2 presents a structural illustration of the methodology applied to ecological risk assessment of the landscape. The structure was formed using a combination of risk assessment techniques—(version 4.2, Oregon State University, Corvallis, OR, USA, available online: http://www.umass.e du/landeco/research/fragstats/fragstats.html) and ArcGIS (version 10.3, ESRI, Redlands, CA, USA, available online: http://www.esri.com). The method involves the following steps:
(1) Step 1: Selecting landscape metrics. The landscape structure in our study area was analyzed using both class-level metrics (each patch type in the given mosaic) and landscape-level metrics (the landscape mosaic as a whole) to investigate and explore landscape scale variables. In addition, reference to previous publications [47,48,49] suggested key landscape metrics (LM) for the assessment of spatial patterns. Since no single LM can be considered [1] for its relevance to structural analysis, we computed 9 class-, 2 landscape-, and 2 class/landscape-level metrics for each map using FRAGSTATS (Table 2). A detailed technical description of each metric is provided for reference in FRAGSTATS [50].
(2) Step 2: Method of constructing optimal scale. Analysis of the grain effect is based on the grain size effect curve of the landscape metrics (LMs) to select the appropriate grain domain [51]. The curve is a horizontal axis with a different grain, and the corresponding LM value is determined on the vertical axis at this grain. Assessment of total area loss in the study area was based on regional land area change evaluation index models [8,52], as expressed in Equations (1) and (2). The conversion grain of 20 m was selected as the starting point, and 140 m as the end point, with an interval of 10 m, and the interpreted vector data were rasterized to obtain a grid map with different grains (Figure 3a–c). According to the area information conservation evaluation (AICE) method [52], using the grain as the abscissa and the regional landscape area change value as the ordinate, the corresponding landscape area change index was calculated separately (Figure 3d). It can be observed that when the grain is less than 50 m, the area loss accuracy value is smaller, and when the grain size is greater than 50 m, the index of area precision loss increases. Consequently, the best analysis grain (50 m) was chosen, which maximizes both the calculation reliability and efficiency.
L i = ( A i A b i ) / A b i × 100
S i = i = 1 n L i 2 n
Here, Ai is i-type landscape grid area; Abi is bi-type landscape vector area before scale conversion; Li denotes is the relative loss of area; Si is the area change index; and n is the number of landscape types.
First, based on the analysis of 50-m grain sizes, 26 moving windows of different sizes were obtained using a 50-m interval (i.e., 50, 150, 250 …, 2350, 2450, and 2550 m). Second, moving from the upper left corner of the study area, one grid at a time to the right, the value of the selected landscape metrics in the unit window was calculated, and then the landscape metrics value was assigned to the center grid of the pane, finally producing the relevant LM grid. Third, 5000 random points were generated using ArcGIS 10.3 for the open coastal wetlands in Jiangsu, and the value of the landscape metrics (Shannon’s evenness index, SHEI and the contagion index, CONTAG) representing the landscape was selected. The value of the landscape metrics of each point was imported into GS+ (version 7.0, Geostatistics for the Environmental Sciences, Gamma Design Software, Michigan, USA), and the block ratio of different moving windows was obtained (Figure 3e). It can be seen that when the radius of the moving window is less than 2150 m, the block ratio of the SHEI and CONTAG is erratic and irregular; and when the moving window radius is larger than or equal to 2150 m, the block ratio of the SHEI and CONTAG remains stable and tends to a certain value. Therefore, the selection of 2150 m was the characteristic scale used in this study.
(3) Step 3: Division into risk subareas. To realize the spatialization of landscape ecological risk index (ERI), based on previous research [53], the OCWJ were divided into subareas with a cell size of 7.29 × 7.29 km using an equidistance sampling method which resulted in 106 grid cells (risk subareas). Then, the LERI of each risk grid was calculated and used as the centroid value of each specific ecological risk subarea gird.
(4) Step 4: Constructing ecological risk indexes. This study introduces Disturbancei, Fragilityi, and Lossi to construct a landscape ecological risk index. Disturbancei is used to reflect the extent to which ecosystems represented by different landscapes are disturbed [20], as expressed in Equation 3. Fragilityi indicates the vulnerability of the ecosystem structure within different landscape types (i.e., the higher the index, the more vulnerable the ecosystem). Based on the characteristics of the OCWJ and previous studies [20,54], and previous studies [20,54] and the vulnerability of different landscape types, Fragilityi was first divided into ten levels, from low to high (building, reedbed, seepweed, cordgrass, saltwork, farmland, seawater, aquafarm, tidalflat, and dry pool), and Fragilityi was obtained by normalization processing (Table 3). Lossi reflects the degree of loss of natural attributes of different landscape types when subjected to external disturbances; it consists of Disturbancei and Fragilityi [55], as expressed in Equation (7). LERIi was constructed based on Lossi in a sample plot so that the landscape spatial pattern was transformed to a spatialized ecological risk variable using the sampling method, as expressed in Equation (8).
D i s t u r b a n c e i = a × F r a g m e n t a t i o n i + b × I s o l a t i o n i + c × D o min a n t i
Here, a, b and c are the weights of indices, set to 0.5, 0.3 and 0.2 according to their importance [20].
F r a g m e n t a t i o n i = N L i T A i
Here, NLi is the number for the ith landscape type, and TAi is the total area of the ith landscape type.
I s o l a t i o n i = T A 2 A L i N P i T A
Here, TA is the landscape total area; ALi is the area of ith land type; and NPi is the number of patches of the ith land type.
D o min a n t i = ( P D i + P F i ) 4 + P R i 2
Here, PDi is the patch density (i.e., the number of patches of the ith landscape type divided by the total number of patches); PFi is the patch frequency (i.e., the number of samples of the ith landscape type divided by the total number of samples); and PRi is the patch ratio (i.e., the area of samples of the ith landscape type be divided by total area of samples).
L o s s i = D i s t u r b a n c e i × F r a g i l i t y i
L E R I i = j = 1 N L A R k j A R k L o s s j
Here, Lossi is the ith landscape loss degree index; Disturbancei is the ith landscape disturbance index; Fragilityi is the landscape fragility index of the ith landscape type; LERIi is the ith subarea’s ecological risk index; Lossj is the jth landscape loss degree index; LARkj is the jth landscape’s area in the kth region; and ARk is an area of the kth region.
(5) Step 5: Method of spatial analysis. Geostatistical methods are statistical methods for detecting, simulating and estimating the correlation and distribution of variables within a particular area of study [5]. With the aid of the semi-variance function in geostatistical methods [20,52] it is used to fit the data of different models’ sampling points, and the spatial distribution map of landscape ecological risk index is plotted using the by kriging method. The natural breaks method was used to compare the changes in ecological risk for different periods in the OCWJ. The risks can be classified into categories such as “extremely low (LERIi < 130.10)”, “low (130.10 ≤ LERIi < 203.94)”, “relatively low (203.94 ≤ LERIi < 262.45)”, “medium (262.45 ≤ LERIi < 294.49)”, “relatively high (294.49 ≤ LERIi < 330.71)”, “high (330.71 ≤ LERIi < 372.50)”, and “extremely high (LERIi ≥ 372.50). On this basis, the distribution map of ecological risk in OCWJ can be used to assess the regional ecological risk.

3. Results

3.1. Landscape Classification and Landscape Pattern Mapping

Ten landscape classes were identified considering the specific local conditions, the national land resources classification system, and following the approach of Xu [27]. These landscape types were modified and verified using ENVI5.2 and ArcGIS10.3. The overall accuracy of each classified image is around 87%, which was verified using field surveys and Google Earth samples. The overall accuracies were within the acceptable range for further analysis [27,56]. The areal extent of ten landscape types (seawater, tidal flat, salt work, farmland, aquafarm, dry pool, building, reedbed, seepweed, and cordgrass) and their spatial distribution, for each time step, are presented in Figure 4. After classification, ten landscape types in different periods were treated as the OCWJ ecological risk receptors and the ecological risk indexes were constructed for them also. Finally, with these available spatial data sets and spatial analysis tools and techniques, the characteristics of a model of the evolution of landscape ecological risk patterns (LERP) of the OCWJ were investigated.

3.2. Landscape Pattern Analysis Using an Optinal Scale

With the intensification of human activities, the landscape pattern, as expressed by various landscape metrics of the individual landscape types of the coastal ecosystems they constitute, changed substantially during 2000–2017 (Figure 4). In 2000, tidal flat, farmland, and aquafarms were the main landscape types in the study area, accounting for 68.75% of the total area, followed by salt works, seawater, dry pool, reedbed, cordgrass, and seepweed. The area of tidalflat was reduced by 458.82 km2 during 2000–2017, with the most significant decrease during 2000–2004 (by 132.11 km2) and during 2008–2013 (by 152.32 km2). Most of the tidal flat area is occupied by aquafarms and farmland. From 2000 to 2017, the area of aquafarms increased by 479.88 km2 (with the largest increase of 14.10%). Evidently, therefore, besides natural factors, individual or artificial management activities were effective agents of landscape change in the study area.
Changes in different landscape types in the study area are also reflected in the landscape pattern metrics (Figure 5). Year mean size (YMS) and year percent (YP) values of the metrics used to define landscape pattern in the study area were calculated based on the standard method, using FRAGSTATS 4.2. Changes in the trend of YMS of the total area (CA), the percent of the landscape (PLAND), and the largest patch index (LPI), are consistent within 2000–2017, indicating the continued dominance of farmland and tidal flats. The patch density (PD) of the man-made landforms of the landscape remained low within the landscape types, indicating the encroachment upon and dramatic fragmentation of the non-man-made landforms landscape. The year mean size (YMS) of the PD, the edge density (ED), and the landscape shape index (LSI) exhibited similar trends, but there are differences in the amplitudes of the peaks. Patch density (PD) and landscape shape index (LSI) were highest in seepweed along the landscape transfer direction, and ED was highest in seawater. In the landscape class percent of the study area, the class year percent (YP) of saltworks in the CA decreased from 30.35% to 11.85% during 2000–2017, while the buildings area increased from zero to 42.06% in 2017. Moreover, the class year percent (YP) of ED and LSI of the landscape types was increasing for seawater, farmland, dry pool, and buildings landscapes, the most notable being buildings and dry pools. In the last 20 years, the class year percent (YP) of ED of the buildings area increased by 37.74%, and the landscape shape index (LSI) of the drypool increased by 40%. This shows that the continuous development of the study area during 2000–2017 and the continuous expansion of the built area have tended to render the landscape pattern of the area more fragmented; in addition, the boundary between buildings and other landscapes is more complex.
Heterogeneity, an inherent landscape attribute, is defined as the complexity and/or variability of a system property in space and/or time [5]. To further analyze the heterogeneity of the landscape, the central and southern regions, with a wide range of landscape categories, were selected for analysis (Figure 6). Using GIS technology, the raster data of 50 m was imported into FRAGSTATS 4.2, and the spatial range of the patch cohesion index (COHESION) and the Euclidean nearest neighbor distance mean index (ENN–MN) visualization were implemented with a moving window of 2150 m. The COHESION in the feature scale space is 88.06–100 (i.e., the degree of tessellation is from low to high) (Figure 6a–e). The main landscape types in the low-value areas of the open coastal wetlands in the Jiangsu (COWJ) landscape mosaic are farmlands, reedbeds, and aquafarms, where farmland and aquafarms are man-made landform landscapes (MMLLs). High-value areas are distributed on the landside of the study area and exhibit a north–south-aligned striped pattern.
From 2000–2004, the ENN–MN decreased within the entire study area (especially on the landward side of Sheyang and Xiangshui), indicating the diffusion of human activity (Figure 6f–j). This was mainly due to the addition of large man-made landform landscapes development units (i.e., salt works). However, from 2004–2008, the southeastern parts of the Sheyang and Dafeng areas exhibited coalescence characteristics in ENN–MN, while the western rural areas of Dafeng areas exhibited a diffusion mode (mainly farmlands and aquafarms). From 2008–2017, the entire center of the study area was again in coalescence mode, perhaps because of the renewed expansion of human activity beyond its original boundaries. Therefore, the diffusion–coalescence process in the OCWJ is unidirectional in time, but is an iterative process. It can be observed that the heterogeneity of the landscape in the OCWJ region was largely caused by human disturbance and the resulting impacts on the coastal wetland ecosystem after the landscape was modified.

3.3. Spatiotemporal Differences in Landscape Ecological Risk Patterns

The numbers of patches and areas of all landscape types in 2000, 2004, 2008, 2013, and 2017 were calculated using ArcGIS10.3 and FRAGSTATS 4.2. Using Equations (3)–(7), the landscape pattern indexes for all landscape types were calculated and summarized in the OCWJ for the five years (Table 3). Maps of landscape types are shown in Figure 7, and they demonstrate that the landscape ecological risk of the OCWJ changed substantially during 2000–2017. The specific statistics for the different risks grades are shown in Figure 8.
In terms of spatial characteristics, there are two main centers of extreme values in the areas of high-grade ecological risk, namely the offshore area of Dongtai in southern Jiangsu and the landward side of Dafeng. The level of ecological risk in the southern region of the open coastal wetlands in Jiangsu (OCWJ) is significantly higher than that of elsewhere. The intensity of human activity between Dafeng and Dongtai gradually increased, resulting in a more frequent succession of landscape types. The landscape types in the northern region of the OCWJ are relatively few, with saltwork as the main areas, accompanied by farmland. The process of development of the man-made landform landscapes (MMLL) was gradual. The disturbance of the landscape by human activities was less in the southeast, and thus the ecological risk level is relatively low. The centers of high values in the central and southern regions are near Dafeng and Dongtai, and the main landscape types are farmland, fish ponds, tidal flats, and construction land. Owing to the construction of MMLL and the intensive development of farmland and tidal flats, the ecological risk value rises. There is a trend of increasing ecological risk on the landward side of areas adjacent to the sea; the main reason for this is the continual movement of cofferdams towards the sea. Overall, diverse human activities have increased the MMLL and decreased the extent of natural landscapes, with a corresponding increase in ecological risk.
In general, in the areas of economic activity and human habitation, the areas with a high degree of landscape transformation also have high levels of ecological risk, while single landscape types have only a low level of ecological risk. The high-grade ecological risk zone remained dominant in the southern part and its area gradually increased. The low-grade (i.e., extremely low, low, and relatively low) ecological risk areas are distributed around the areas of medium- and high-grade ecological risk (i.e., relatively high, high, extremely high) and they generally expanded in a circular pattern. This shows that there was a diverse pattern of the formation of man-made landform landscapes (MMLL) following reclamation of the coastal wetland, landscape fragmentation was intensified, and the average patch area of natural landscapes such as seepweed was greatly reduced. At the same time, the reduction in the tidal flat area also impacted the coastal ecosystem. Therefore, it is important to pay attention to wetland protection during the development and utilization of open coastal wetlands in Jiangsu (OCWJ).
In the time series, the low-grade ecological risk area of the OCWJ decreased by 1058.55 km2 during 2000–2017. The area of extremely low ecological risk was reduced by 82.38 km2; the area of low ecological risk was reduced by 757.06 km2; and the area of relatively low ecological risk was reduced by 219.12 km2. Medium and relatively high are the dominant types of ecological risk in the study area, accounting for 20.51% and 28.99% of the total area in 2017, respectively. At the same time, the area of high ecological risk area increased by 168.88 km2.

3.4. Landscape Ecological Risk Transfer Analysis

To better analyze the transformation relationship between landscape ecological risk areas of open coastal wetlands in Jiangsu (OCWJ) in different periods, an ecological risk grade transfer matrix was used to study the changes in the areas of ecological risk area at each grade. ArcGIS10.3 was used to superimpose the ecological risk grade areas of different periods, and the results of grade area conversion of ecological risk areas of the OCWJ during 2000–2017 were analyzed (Table 4 and Figure 9).
Overall, the change in landscape ecological risk levels in the OCWJ from 2000–2017 is substantial, and the areas of low, relatively low, extremely low, and medium ecological risk have all undergone large changes. The total area of transformation from low-grade to high-grade ecological risk is approximately 1910.71 km2, accounting for 71.75% of the total area converted. The total level of ecological risk from high-grade to low-grade is about 752.23 km2, accounting for 28.25% of the total study area. The area converted from low-grade to high-grade ecological risk is 2.54 times that of the area converted from high-grade to low-grade, indicating that the ecological risk in the study area has increased. Moreover, the most significant areas of ecological risk grade area conversion are those converted from low to medium ecological risk and from low to relatively high ecological risk, with conversion areas of 331.74 km2 and 273.45 km2, respectively. These conversions mainly occur in areas where the landscape type conversion is relatively large (e.g., the southeast coastal area of Sheyang, Dafeng and Dongtai), where human activities are pronounced. The MMLL related to fisheries is constantly expanding and the ecological risk is increasing.
A comparison of the average annual conversion rate of ecological risk levels of the OCWJ for different time periods is illustrated in Figure 9. It is evident that in 2000–2004, 2004–2008, 2008–2013, and 2013–2017 the annual average conversion rate from high-grade to low-grade ecological risk (i.e., LOW-EXL, REL-EXL, REL-LOW, MED-EXL, MED-LOW, MED-REL, REH-EXL, REH-LOW, REH-REL, REH-MED, HIG-EXL, HIG-LOW, HIG-REL, HIG-MED, HIG-REH, EXH-EXL, EXH-LOW, EXH-REL, EXH-MED, EXH-REH, and HIG-REH) exhibited a decreasing trend. Extremely high to HIG changed the most, and the average annual transfer rate decreased by 48.37 km2/year. In contrast, the average annual conversion rate of ecological risk level from low-grade to high-grade (i.e., EXL-LOW, EXL-REL, EXL-MED, EXL-REH, EXL-HIG, EXL-EXH, LOW-REL, LOW-MED, LOW-REH, LOW-HIG, LOW-EXH, REL-MED, REL-REH, REL-HIG, REL-EXH, MED-REH, MED-HIG, MED-EXH, REH-HIG, REH-EXH, and HIG-EXH) exhibited a fluctuating but increasing trend. The average annual conversion rate from LOW-REL reached 103.17 km2/year during 2000–2004. The average annual conversions rates from LOW-MED and from REL-MED were 60.85 km2/year and 73.01 km2/year, respectively, during 2008–2013. However, the average annual conversion rate from 2013–2017 was generally low, indicating that impact of human activities and occupancy on the coastal zone will gradually decrease in the future. In general, the expansion of ecological risk levels from low-grade to high-grade areas increases the level of ecological risk in the OCWJ, and thus it is necessary to design and implement rational land use planning strategies and policies to protect the ecological environment.

4. Discussion

Remote sensing techniques (RS) and geographic information systems (GIS) are an effective means of characterizing the evolution of landscape ecological risk patterns. However, an important problem is how to develop evolution models of landscape ecological risk patterns (LERP) with an optimal scale, auxiliary landscape metrics, and geographic information technology. We used appropriate grain and extent factors to select a multi-metrics optimal scale method (MMOS) for decreasing scale effect, to determine the appropriate combinations of landscape ecological risk indexes for constructing optimal LERP models in the OCWJ region. We found that human activities have significantly impacted the coastal landscape (Figure 4) and the pattern of evolution of ecological risk (Figure 7). These findings extend those of previous assessments of ecological risk on a regional scale and they provide an additional reference study for studies of changes in coastal zones worldwide.
Landscape ecological risk assessment (LERA) is useful for representing landscape ecological security, which is a critical domain of landscape sustainability [57]. The construction of an ecological risk model is based on the following considerations. First, the determination of indicators in the ecological risk model of this study focuses on the characterization of landscape disturbances or risk exposure process, highlighting the geographical-ecological significance of the indicators. Second, more indicators are derived from remote sensing data or model-based calculations, which enhance the objectivity and accuracy of the landscape ecological functions reflected by the indicators. Third, the introduction of spatial analysis embodies the geographical research concept of time conversion space, enriches the landscape ecological risk assessment system, and embodies the dynamic concept of ecological sustainability. Therefore, the landscape ecological risk assessment model based on landscape pattern indexes minimizes the subjective cognitive uncertainty of the researcher and the sensitivity of the indicators to change. Previous work has documented the effectiveness of LERA in assessing the quality of the eco-environment of the same area [17,55]. However, these studies have either been of areas other than coastal zones, or they have not considered the evolution of the spatial characteristics of landscape and landscape ecological risk at optimal scales. We compared our LERP models results with conventional analyses of the eco-environment of the same area [24,27,58], and the results indicate that in terms of ecological significance, our LERP models have a similar or better performance for observing landscape ecological risk patterns. Consequently, a LERP model together with the MMOS method is a cost-effective means of assessing landscape transformation risk in open coastal wetlands, and our results demonstrate the effectiveness of this approach for studying the long-term effects of anthropogenic activities in the coastal zone.
In principle, LERA based on spatial patterns is suitable for all types of risk integrative assessment, but it may not apply to those assessments that emphasize a specific process. The source–sink landscape model that seeks a quantitative expression of the coupling of the landscape pattern and ecological processes is an effective way of integrating patterns and processes [59]. Source-sink landscape types and risk processes can avoid the drawbacks of the LERA method that only relies on landscape pattern indexes. Therefore, risk assessment based on the source-sink landscape model is an important trend in future use of the LERP model. However, for different types of risk stress, the contributions from the source-sink landscape can even be transformed into one other, which increases the quantitative difficulty of the evaluation. Therefore, the risk peeling of the dominant ecological processes on the landscape scale, and the correct identification of the source and sink landscape relationships, are the main difficulties in applying the source landscape model to solve practical problems [60].
Our study also has several limitations that should be addressed in future research. First, the method of weight assignment largely depends on our empirical statistical work and the analysis of macro risk factors. Theoretical explanations of specific scale effects cannot be fully determined and need to be further investigated. Second, only spatial patterns were addressed in our study, with the neglect of the possible effects of stress specific processes; in particular, a means of integrating landscape patterns and ecological factors is needed. In fact, source-sink landscape types, especially risky process variables, have great potential for improving landscape ecological risk assessment models [20]. We intend to conduct additional research on the more specific risks of stress and the role of the object in our landscape ecological risk assessment models. Third, due to the long-term development and utilization of human beings and the fragmentation of habitats, many landscapes have been severely damaged. The landscape types are both intermittent and connected in a spatial distribution, and species diversity becomes a major attribute of this connection [5]. Therefore, by studying the landscape patterns characteristics and the evolution of ecological risks, the differences between species diversity of landscape types, the relationship with environmental factors, and their spatial distribution patterns can be revealed to some extent. In addition, the limitations of the research data acquisition, a more detailed analysis of species biodiversity content cannot be conducted. In future work we will explore whether changes in landscape pattern diversity and species diversity can be detected with the combination of higher resolution images and vegetation information. Fourth, in future investigations the division of ecological risk units needs to pay more attention to the integrated geographical significance of the evaluation units to promote the application of the risk assessment LERP models in other coastal areas.

5. Conclusions

We have constructed an optimal scale method with an appropriate extent and grain using multi-temporal Landsat TM/OLI images acquired in the years 2000, 2004, 2008, 2013, and 2017. We then developed landscape ecological risk indices. Based on this model, the entire process of spatiotemporal evolution of ecological risk patterns of open coastal wetlands in Jiangsu (OCWJ) was mapped and analyzed. Our main conclusions are as follows: The main landscape types in the study area are tidal flats and farmland, and the main features of the landscape evolution are that the area of aquafarming has increased substantially at the expense of the tidal flat area, and the landscape heterogeneity has become more pronounced. The heterogeneity of the landscape in the OCWJ area is mainly due to human disturbance and the impact of landscape reconstruction on the coastal wetland ecosystem. In the past 20 years, the areas of low and relatively low ecological risk in the study area have been greatly reduced and the areas of medium, relatively high and high ecological risk have greatly increased. The areas of high-grade ecological risk areas are mainly distributed in Dongtai and in the area of Dafeng. The area transformed from low-grade to high-grade ecological risk accounted for 71.75% of the study area during period 2000–2017. In the previous periods (2000–2004 and 2004–2008), the area of low-grade ecological risk was transformed to middle-grade ecological risk, and in the later periods (2008–2013 and 2013–2017) there was an increase in the proportion of ecological risk areas which were transformed to areas of high-grade ecological risk.
In general, reducing the ecological risk protects the original natural landscape, and a reduction in the rate of increase of man-made landform landscapes (MMLL) is crucial. The results of this study potentially provide important insights and data for coastal landscape planners, researchers, and managers; for linking satellite data with relevant policy formulation; for enhancing the official database of coastal landscape planning; and for assessing the potential effects of the MMLL process on coastal environments. They also provide an important reference for studies of change in coastal zones worldwide.

Author Contributions

Y.L. (Yongchao Liu), Y.L. (Yongxue Liu) and J.L. conceived and designed the framework of this research work; Y.L. (Yongchao Liu) wrote the paper; W.L. and X.W. analyzed the data; C.S. helped to collect the data and edited the text.

Funding

This research was supported by the Jiangsu Provincial Natural Science Foundation (No. BK20160023), the National Natural Science Foundation of China (No. 41471068), NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization: No. U1609203, and the Project Supported by the Open Fund of Key Laboratory of Coastal Zone Exploitation and Protection, Ministry of Land and Resource (No: 2017CZEPK12).

Acknowledgments

We thank the two anonymous reviewers and the associate editor for their constructive suggestions for improving the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Inkoom, J.N.; Frank, S.; Greve, K.; Walz, U.; Fürst, C. Suitability of different landscape metrics for the assessments of patchy landscapes in West Africa. Ecol. Indic. 2018, 85, 117–127. [Google Scholar] [CrossRef]
  2. Martin, M.J.R.; De Pablo, C.L.; De Agar, P.M. Landscape changes over time: Comparison of land uses, boundaries and mosaics. Landsc. Ecol. 2006, 21, 1075–1088. [Google Scholar] [CrossRef]
  3. Wu, J. Landscape sustainability science: Ecosystem services and human well-being in changing landscapes. Landsc. Ecol. 2013, 28, 999–1023. [Google Scholar] [CrossRef]
  4. Peng, J.; Wang, A.; Liu, S.; Wu, J.; Li, W. Research on landscape ecology and sustainable land use. Acta Sci. Nat. Univ. Pekin. 2004, 40, 154–160. [Google Scholar] [CrossRef]
  5. Wu, J. Landscape Ecology Pattern, Process, Scale and Hierarchy; Higher Education Press: Beijing, China, 2007. [Google Scholar]
  6. Pickett, S.T.A.; Cadenasso, M.L. Landscape ecology—Spatial heterogeneity in ecological-systems. Science 1995, 269, 331–334. [Google Scholar] [CrossRef] [PubMed]
  7. Xie, H.; Wang, P.; Huang, H. Ecological risk assessment of land use change in the Poyang Lake Eco-economic Zone, China. Int. J. Environ. Res. Public Health 2013, 10, 328–346. [Google Scholar] [CrossRef] [PubMed]
  8. Zhang, L.; Zhao, Y.; Yin, S.; Fang, S.; Liu, X.; Pu, M. Gradient analysis of dry valley of minjiang river landscape pattern, based on moving window method. Acta Ecol. Sin. 2014, 34, 3276–3284. [Google Scholar] [CrossRef]
  9. Galpern, P.; Manseau, M. Finding the functional grain: Comparing methods for scaling resistance surfaces. Landsc. Ecol. 2013, 28, 1269–1281. [Google Scholar] [CrossRef]
  10. Crutzen, P.J. Geology of mankind. Nature 2002, 415, 23. [Google Scholar] [CrossRef] [PubMed]
  11. Crutzen, P.; Stoermer, E. The “anthropocene”. IGBP Newsl. 2000, 41, 17–18. [Google Scholar]
  12. Li, J.; Yang, L.; Pu, R.; Liu, Y. A review on anthropogenic geomorphology. J. Geogr. Sci. 2017, 27, 109–128. [Google Scholar] [CrossRef]
  13. Ju, C.; Jia, Y.; Shan, H.; Tang, C.; Ma, W. Gis-based coastal area suitability assessment of geo-environmental factors in laoshan district, Qingdao. Nat. Hazards Earth Syst. Sci. 2012, 12, 143–150. [Google Scholar] [CrossRef]
  14. Klemas, V. Remote sensing of coastal plumes and ocean fronts: Overview and case study. J. Coast. Res. 2012, 28, 1–7. [Google Scholar] [CrossRef]
  15. Cohen, W.B.; Goward, S.N. Landsat’s role in ecological applications of remote sensing. Bioscience 2004, 54, 535–545. [Google Scholar] [CrossRef]
  16. Munns, W.R., Jr. Assessing risks to wildlife populations from multiple stressors: Overview of the problem and research needs. Ecol. Soc. 2006, 11, 23. [Google Scholar] [CrossRef]
  17. Peng, J.; Zong, M.; Hu, Y.; Liu, Y.; Wu, J. Assessing landscape ecological risk in a mining city: A case study in Liaoyuan city, China. Sustainability 2015, 7, 8312–8334. [Google Scholar] [CrossRef]
  18. Zhou, Q.; Wang, R. Ecological risk assessment of urbanization processes—A case study. Acta Ecol. Sin. 1998, 18, 337–342. [Google Scholar]
  19. Heller, S. Managing industrial risk—Having a tested and proven system to prevent and assess risk. J. Hazard. Mater. 2006, 130, 58–63. [Google Scholar] [CrossRef] [PubMed]
  20. Li, J.; Pu, R.; Gong, H.; Luo, X.; Ye, M.; Feng, B. Evolution characteristics of landscape ecological risk patterns in coastal zones in Zhejiang province, China. Sustainability 2017, 9, 584. [Google Scholar] [CrossRef]
  21. Gaines, K.F.; Porter, D.E.; Dyer, S.A.; Wein, G.R.; Pinder, J.E.; Brisbin, I.L. Using wildlife as receptor species: A landscape approach to ecological risk assessment. Environ. Manag. 2004, 34, 528–545. [Google Scholar] [CrossRef]
  22. Tannenbaum, L.V. Detoxifying ecological risk assessment. Hum. Ecol. Risk Assess. 2005, 11, 469–472. [Google Scholar] [CrossRef]
  23. Xu, Y.; Cao, K.; Li, M.; Xu, Z. Coastal ecological risk assessment: Its research progress and prospect. Adv. Earth Sci. 2016, 31, 137–146. [Google Scholar]
  24. Chuai, X.; Huang, X.; Wu, C.; Li, J.; Lu, Q.; Qi, X.; Zhang, M.; Zuo, T.; Lu, J. Land use and ecosystems services value changes and ecological land management in coastal Jiangsu, China. Habitat Int. 2016, 57, 164–174. [Google Scholar] [CrossRef]
  25. Li, Y.; Zhu, X.; Sun, X.; Wang, F. Landscape effects of environmental impact on bay-area wetlands under rapid urban expansion and development policy: A case study of Lianyungang, China. Landsc. Urban Plan. 2010, 94, 218–227. [Google Scholar] [CrossRef]
  26. Li, Y.; Shi, Y.; Zhu, X.; Cao, H.; Yu, T. Coastal wetland loss and environmental change due to rapid urban expansion in Lianyungang, Jiangsu, China. Reg. Environ. Chang. 2014, 14, 1175–1188. [Google Scholar] [CrossRef]
  27. Xu, C.; Pu, L.; Zhu, M.; Li, J.; Chen, X.; Wang, X.; Xie, X. Ecological security and ecosystem services in response to land use change in the coastal area of Jiangsu, China. Sustainability 2016, 8, 816. [Google Scholar] [CrossRef]
  28. Boongaling, C.G.K.; Faustino-Eslava, D.V.; Lansigan, F.P. Modeling land use change impacts on hydrology and the use of landscape metrics as tools for watershed management: The case of an ungauged catchment in the philippines. Land Use Policy 2018, 72, 116–128. [Google Scholar] [CrossRef]
  29. Leonard, P.B.; Sutherland, R.W.; Baldwin, R.F.; Fedak, D.A.; Carnes, R.G.; Montgomery, A.P. Landscape connectivity losses due to sea level rise and land use change. Anim. Conserv. 2017, 20, 80–90. [Google Scholar] [CrossRef]
  30. Deng, J.; Wang, K.; Hong, Y.; Qi, J. Spatio-temporal dynamics and evolution of land use change and landscape pattern in response to rapid urbanization. Landsc. Urban Plan. 2009, 92, 187–198. [Google Scholar] [CrossRef]
  31. Torresan, S.; Critto, A.; Valle, M.D.; Harvey, N.; Marcomini, A. Assessing coastal vulnerability to climate change: Comparing segmentation at global and regional scales. Sustain. Sci. 2008, 3, 45–65. [Google Scholar] [CrossRef]
  32. Islam, M.A.; Mitra, D.; Dewan, A.; Akhter, S.H. Coastal multi-hazard vulnerability assessment along the ganges deltaic coast of bangladesh—A geospatial approach. Ocean Coast. Manag. 2016, 127, 1–15. [Google Scholar] [CrossRef]
  33. Addo, K.A. Assessing coastal vulnerability index to climate change: The case of Accra-Ghana. J. Coast. Res. 2013, 65, 1892–1897. [Google Scholar] [CrossRef]
  34. Wang, C.; Sun, G.; Dang, L. Identifying ecological red lines: A case study of the coast in Liaoning province. Sustainability 2015, 7, 9461–9477. [Google Scholar] [CrossRef]
  35. He, Q.; Bertness, M.D.; Bruno, J.F.; Li, B.; Chen, G.; Coverdale, T.C.; Altieri, A.H.; Bai, J.; Sun, T.; Pennings, S.C.; et al. Economic development and coastal ecosystem change in China. Sci. Rep. 2014, 4, 5995. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, A.; Groenewold, N. Reducing regional disparities in china: An evaluation of alternative policies. J. Comp. Econ. 2010, 38, 189–198. [Google Scholar] [CrossRef]
  37. Gong, J.; Yang, J.; Tang, W. Spatially explicit landscape-level ecological risks induced by land use and land cover change in a national ecologically representative region in China. Int. J. Environ. Res. Public Health 2015, 12, 14192–14215. [Google Scholar] [CrossRef] [PubMed]
  38. Ives, A.R.; Carpenter, S.R. Stability and diversity of ecosystems. Science 2007, 317, 58–62. [Google Scholar] [CrossRef] [PubMed]
  39. Carpenter, S.R.; Mooney, H.A.; Agard, J.; Capistrano, D.; De Fries, R.S.; Díaz, S.; Dietz, T.; Duraiappah, A.K.; Oteng-Yeboah, A.; Pereira, H.M. Science for managing ecosystem services: Beyond the millennium ecosystem assessment. Proc. Natl Acad. Sci. USA 2009, 106, 1305–1312. [Google Scholar] [CrossRef] [PubMed]
  40. Luo, Y. Sustainability associated coastal eco-environmental problems and coastal science development in China. Bull. Chin. Acad. Sci. 2016, 31, 1133–1142. [Google Scholar]
  41. Liu, Y.; Li, M.; Mao, L.; Cheng, L.; Chen, K. Seasonal pattern of tidal-flat topography along the Jiangsu middle coast, China, using HJ-1 optical images. Wetlands 2013, 33, 871–886. [Google Scholar] [CrossRef]
  42. Sun, C.; Liu, Y.; Zhao, S.; Li, H.; Sun, J. Saltmarshes response to human activities on a prograding coast revealed by a dual-scale time-series strategy. Estuar. Coasts 2017, 40, 522–539. [Google Scholar] [CrossRef]
  43. Herold, M.; Goldstein, N.C.; Clarke, K.C. The spatiotemporal form of urban growth: Measurement, analysis and modeling. Remote Sens. Environ. 2003, 86, 286–302. [Google Scholar] [CrossRef]
  44. Wulder, M.A.; Coops, N.C. Make earth observations open access. Nature 2014, 513, 30–31. [Google Scholar] [CrossRef] [PubMed]
  45. Kennedy, R.E.; Andrefouet, S.; Cohen, W.B.; Gomez, C.; Griffiths, P.; Hais, M.; Healey, S.P.; Helmer, E.H.; Hostert, P.; Lyons, M.B.; et al. Bringing an ecological view of change to landsat-based remote sensing. Front. Ecol. Environ. 2014, 12, 339–346. [Google Scholar] [CrossRef] [Green Version]
  46. Wu, J.; Xiang, W.; Zhao, J. Urban ecology in China: Historical developments and future directions. Landsc. Urban Plan. 2014, 125, 222–233. [Google Scholar] [CrossRef]
  47. Fan, C.; Myint, S. A comparison of spatial autocorrelation indices and landscape metrics in measuring urban landscape fragmentation. Landsc. Urban Plan. 2014, 121, 117–128. [Google Scholar] [CrossRef]
  48. Niesterowicz, J.; Stepinski, T.F. On using landscape metrics for landscape similarity search. Ecol. Indic. 2016, 64, 20–30. [Google Scholar] [CrossRef]
  49. Wang, Z.; Tang, L.; Qiu, Q.; Chen, H.; Wu, T.; Shao, G. Assessment of regional ecosystem health—A case study of the golden triangle of southern Fujian province, China. Int. J. Environ. Res. Public Health 2018, 15, 802. [Google Scholar] [CrossRef] [PubMed]
  50. McGarigal, K. Fragstats Help. Available online: http://www.umass.edu/landeco/research/fragstats/documents/fragstats.help.4.2.pdf (accessed on 12 January 2018).
  51. Zhao, W.; Fu, B.; Chen, L. The effects of grain change on landscape indices. Quat. Sci. 2003, 23, 326–333. [Google Scholar]
  52. Xu, Z.; Hu, Y.; Liu, Y.; Yan, Y. A review on the accuracy analysis of spatial scaling data. Prog. Geogr. 2012, 31, 1574–1582. [Google Scholar]
  53. Leuven, R.; Poudevigne, I. Riverine landscape dynamics and ecological risk assessment. Freshw. Biol. 2002, 47, 845–865. [Google Scholar] [CrossRef]
  54. Xu, Y.; Gao, J.; Zhao, J.; Chen, J. The research progress and prospect of watershed ecological risk assessment. Acta Ecol. Sin. 2012, 32, 284–292. [Google Scholar]
  55. Fan, J.; Wang, Y.; Zhou, Z.; You, N.; Meng, J. Dynamic ecological risk assessment and management of land use in the middle reaches of the Heihe river based on landscape patterns and spatial statistics. Sustainability 2016, 8, 536. [Google Scholar] [CrossRef]
  56. Kindu, M.; Schneider, T.; Teketay, D.; Knoke, T. Land use/land cover change analysis using object-based classification approach in munessa-shashemene landscape of the ethiopian highlands. Remote Sens. 2013, 5, 2411–2435. [Google Scholar] [CrossRef]
  57. Gao, Y.; Wu, Z.; Lou, Q.; Huang, H.; Cheng, J.; Chen, Z. Landscape ecological security assessment based on projection pursuit in Pearl River Delta. Environ. Monit. Assess. 2012, 184, 2307–2319. [Google Scholar] [CrossRef] [PubMed]
  58. Zang, Z.; Zou, X.; Zuo, P.; Song, Q.; Wang, C.; Wang, J. Impact of landscape patterns on ecological vulnerability and ecosystem service values: An empirical analysis of yancheng nature reserve in China. Ecol. Indic. 2017, 72, 142–152. [Google Scholar] [CrossRef]
  59. Stow, C.A.; Borsuk, M.E.; Stanley, D.W. Long-term changes in watershed nutrient inputs and riverine exports in the Neuse river, North Carolina. Water Res. 2001, 35, 1489–1499. [Google Scholar] [CrossRef]
  60. Jian, P.; Dang, W.; Liu, Y.; Zong, M.; Hu, X. Review on landscape ecological risk assessment. Acta Geogr. Sin. 2015, 70, 664–677. [Google Scholar]
Figure 1. (ac) Location of the study area. The main natural wetlands of the open coastal wetlands in Jiangsu (OCWJ) comprise tidalflat (d-1); cordgrass and seawater (d-2); seepweed and reedbed (d-3); and farmland (d-4) [42], an exotic introduction from North America which outcompeted the native wetlands and became a dominant species after 1990. The area bounded by the yellow line (c) is the spatial heterogeneity of landscape space.
Figure 1. (ac) Location of the study area. The main natural wetlands of the open coastal wetlands in Jiangsu (OCWJ) comprise tidalflat (d-1); cordgrass and seawater (d-2); seepweed and reedbed (d-3); and farmland (d-4) [42], an exotic introduction from North America which outcompeted the native wetlands and became a dominant species after 1990. The area bounded by the yellow line (c) is the spatial heterogeneity of landscape space.
Ijerph 15 01691 g001
Figure 2. Framework of the methodology used for the landscape ecological risk assessment.
Figure 2. Framework of the methodology used for the landscape ecological risk assessment.
Ijerph 15 01691 g002
Figure 3. Changes of landscape metric values with grain (ac) and changes in area landscape accuracy with different grain (d). With the increase in grain size, due to the boundary of the vector data and the changes of the properties of adjacent patches, scale turning points appear in the values of landscape metrics. Considering the distribution of inflection points and the calculation efficiency (ac), 20–100 m was selected as the appropriate scale domain to determine the analysis grain for the next step (d). The block ratio of different moving windows is obtained (e).
Figure 3. Changes of landscape metric values with grain (ac) and changes in area landscape accuracy with different grain (d). With the increase in grain size, due to the boundary of the vector data and the changes of the properties of adjacent patches, scale turning points appear in the values of landscape metrics. Considering the distribution of inflection points and the calculation efficiency (ac), 20–100 m was selected as the appropriate scale domain to determine the analysis grain for the next step (d). The block ratio of different moving windows is obtained (e).
Ijerph 15 01691 g003
Figure 4. Landscape classification during five periods from 2000–2017 in the open coastal wetlands in Jiangsu, China. (a) Map of landscape types in the year 2000; (b) map of landscape types in the year 2004; (c) map of landscape types in the year 2008; (d) map of landscape types in the year 2013; (e) map of landscape types in the year 2017.
Figure 4. Landscape classification during five periods from 2000–2017 in the open coastal wetlands in Jiangsu, China. (a) Map of landscape types in the year 2000; (b) map of landscape types in the year 2004; (c) map of landscape types in the year 2008; (d) map of landscape types in the year 2013; (e) map of landscape types in the year 2017.
Ijerph 15 01691 g004
Figure 5. Overall characteristics of landscape types. (a) CA; (b) PLAND; (c) PD; (d) LPI; (e) ED; and (f) LSI. To fully reflect the characteristics of the landscape pattern and to reduce the redundancy of the pattern information, the landscape metrics at landscape level were selected, and the characteristics of the landscape pattern were analyzed using the FRAGSTATS 4.2 selection method.
Figure 5. Overall characteristics of landscape types. (a) CA; (b) PLAND; (c) PD; (d) LPI; (e) ED; and (f) LSI. To fully reflect the characteristics of the landscape pattern and to reduce the redundancy of the pattern information, the landscape metrics at landscape level were selected, and the characteristics of the landscape pattern were analyzed using the FRAGSTATS 4.2 selection method.
Ijerph 15 01691 g005
Figure 6. (ae) Comparison of the landscape heterogeneity of COHESION during five intervals from 2000–2017. The degree of mosaicism of the landscape pattern is represented by the degree of COHESION. The greater the COHESION value, the higher the degree of patch connectivity of the landscape and the lower the degree of spatial mosaic; (fj) Comparison of the landscape heterogeneity of ENN–MN during five intervals from 2000–2017. Distance analysis of the landscape is expressed by ENN–MN. The larger the value of ENN–MN, the greater the distance between the landscapes type blocks and the more scattered the distribution.
Figure 6. (ae) Comparison of the landscape heterogeneity of COHESION during five intervals from 2000–2017. The degree of mosaicism of the landscape pattern is represented by the degree of COHESION. The greater the COHESION value, the higher the degree of patch connectivity of the landscape and the lower the degree of spatial mosaic; (fj) Comparison of the landscape heterogeneity of ENN–MN during five intervals from 2000–2017. Distance analysis of the landscape is expressed by ENN–MN. The larger the value of ENN–MN, the greater the distance between the landscapes type blocks and the more scattered the distribution.
Ijerph 15 01691 g006
Figure 7. Evolution of the spatial distribution of ecological risk index grades of the open coastal wetlands in Jiangsu (OCWJ) during the five periods from 2000–2017. (a) Map of the spatial distribution of ecological risk index grades in the year 2000; (b) map of the spatial distribution of ecological risk index grades in the year 2004; (c) map of the spatial distribution of ecological risk index grades in the year 2008; (d) map of the spatial distribution of ecological risk index grades in the year 2013; (e) map of the spatial distribution of ecological risk index grades in the year 2017.
Figure 7. Evolution of the spatial distribution of ecological risk index grades of the open coastal wetlands in Jiangsu (OCWJ) during the five periods from 2000–2017. (a) Map of the spatial distribution of ecological risk index grades in the year 2000; (b) map of the spatial distribution of ecological risk index grades in the year 2004; (c) map of the spatial distribution of ecological risk index grades in the year 2008; (d) map of the spatial distribution of ecological risk index grades in the year 2013; (e) map of the spatial distribution of ecological risk index grades in the year 2017.
Ijerph 15 01691 g007
Figure 8. Changes in the areas of different ecological risk grades from 2000–2017. Based on the analysis of the ecological risk level of open coastal wetlands in Jiangsu, the area occupied by the ecological risk grades of 106 risk communities in the study area were analyzed statistically.
Figure 8. Changes in the areas of different ecological risk grades from 2000–2017. Based on the analysis of the ecological risk level of open coastal wetlands in Jiangsu, the area occupied by the ecological risk grades of 106 risk communities in the study area were analyzed statistically.
Ijerph 15 01691 g008
Figure 9. Annual transfer rates of different ecological risk grades during 2000–2014, 2004–2008, 2008–2013, and 2013–2017 (km2/year). (a) Map of annual transfer rates of different ecological risk grades transposed from EXL, LOW, and REL; (b) map of annual transfer rates of different ecological risk grades transposed from MED, REH, HIG, and EXH. Notes: EXL: extremely low; LOW: low; REL: relatively low; MED: medium, REH: relatively high; HIG: high; EXH: extremely high.
Figure 9. Annual transfer rates of different ecological risk grades during 2000–2014, 2004–2008, 2008–2013, and 2013–2017 (km2/year). (a) Map of annual transfer rates of different ecological risk grades transposed from EXL, LOW, and REL; (b) map of annual transfer rates of different ecological risk grades transposed from MED, REH, HIG, and EXH. Notes: EXL: extremely low; LOW: low; REL: relatively low; MED: medium, REH: relatively high; HIG: high; EXH: extremely high.
Ijerph 15 01691 g009
Table 1. Details of the landsat remote sensing data used in this study (spatial resolution: 30 m).
Table 1. Details of the landsat remote sensing data used in this study (spatial resolution: 30 m).
SatelliteSensorsPath/RowDateSatelliteSensorsPath/RowDate
Landsat 5TM120/3613 December 2000Landsat 5TM119/3706 December 2000
Landsat 5TM120/3608 December 2004Landsat 5TM119/3715 November 2004
Landsat 5TM120/3619 December 2008Landsat 5TM119/3713 January 2009
Landsat 8OLI120/3601 December 2013Landsat 8OLI119/3710 December 2013
Landsat 8OLI120/3606 December 2017Landsat 8OLI119/3722 December 2017
Table 2. Commonly-used landscape metrics.
Table 2. Commonly-used landscape metrics.
LevelLandscape MetricSignificance/DescriptionThemes/Method
ClassNumber of patches (NP)Total number of patches in the landscapea Measure of configuration/standard method
Landscape shape index (LSI)Landscape boundary and total edge within the landscape divided by the total areaa/b Measure of configuration/standard method
Landscape division index (DIVISION)Degree of plaque fragmentation in the landscape, describing the complexity of the landscapea Diversity assessment/standard method
Shannon’s diversity index (SHDI)Heterogeneity of the landscape: the higher the landscape heterogeneity, the higher the diversity indexa Diversity assessment/standard method
Percentage of the landscape (PLAND)Percentage of the landscape comprising of the corresponding patch typeb Measure of configuration/standard method
Patch density (PD)Configuration metrics patch densityb Spatial differentiation/standard method
Total Area (CA)Sum of the areas of all patches belonging to a given classb Measure of configuration/standard method
Largest patch index (LPI)Percentage of the landscape comprised by the largest patch of the corresponding patch typeb Measure of configuration/standard method
Edge Density (ED)Amount of edge relative to the landscape areab Spatial differentiation/standard method
landscapePatch Cohesion Index (COHESION)Measures the physical connectedness of the corresponding patch typed Measure of configuration/moving window
Euclidean nearest neighbor distance mean index (ENN–MN)Distance to the nearest neighboring patch of the same type, based on shortest edge-to-edge distanced Measure of configuration /moving window
Class/landscapeShannon’s evenness index (SHEI)Expresses the uniformity of distribution at different landscape levelsa/c Measure of configuration/standard method and moving window
Contagion index (CONTAG)Measures the extent to which patch types are aggregated or clumpeda/c Measure of configuration/moving window
Note: a applicable to grain and extent; b applicable to overall landscape pattern characteristics; c applicable to block to base ratio; d applicable to landscape heterogeneity.
Table 3. Summary of landscape pattern indices for the study area.
Table 3. Summary of landscape pattern indices for the study area.
YearTypesSEATIDSALFARAQUDRYBUIREESEECOR
2000Fragmentationi0.040.010.010.000.010.010.000.050.130.01
Isolationi0.410.080.050.030.110.290.000.511.040.28
Dominanti2.138.373.289.403.500.750.001.621.121.39
Disturbancei0.571.700.671.890.730.240.000.500.600.37
Lossi0.070.250.060.210.120.040.000.020.030.03
2004Fragmentationi0.030.010.000.000.000.010.010.060.160.01
Isolationi0.370.100.040.040.070.381.410.741.350.29
Dominanti2.087.192.759.915.870.490.060.950.861.39
Disturbancei0.551.470.562.001.200.220.440.440.660.37
Lossi0.070.210.050.220.200.040.010.020.040.03
2008Fragmentationi0.030.010.000.000.000.010.010.020.100.02
Isolationi0.330.130.040.030.060.210.450.411.030.30
Dominanti2.126.282.129.716.871.010.221.010.881.33
Disturbancei0.541.300.441.951.390.270.180.340.530.36
Lossi0.070.190.040.210.230.050.000.010.030.03
2013Fragmentationi0.030.010.000.000.000.010.000.030.100.02
Isolationi0.300.130.090.040.050.180.170.441.340.37
Dominanti2.474.901.3910.347.611.030.671.270.521.35
Disturbancei0.601.020.312.081.540.260.190.400.550.39
Lossi0.080.150.030.230.250.050.000.020.030.03
2017Fragmentationi0.030.010.000.000.000.040.000.030.100.01
Isolationi0.300.120.110.040.050.450.220.451.420.27
Dominanti2.374.261.2810.237.811.740.671.310.521.51
Disturbancei0.580.890.292.061.580.500.200.410.580.39
Lossi0.070.130.030.220.260.090.000.020.030.03
Notes: Fragilityi is the landscape fragility index; the values are as follows: SEA: seawater (0.13), TID: tidalflat (0.15), SAL: saltwork (0.09), FAR: farmland (0.11), AQU: aquafarm (0.16), DRY: drypool (0.18), BUI: building (0.02), REE: reedbed (0.04), SEE: seepweed (0.06), COR: cordgrass (0.07).
Table 4. Transition matrix for changes in the area (km2) of landscape types in the OCWJ during 2000–2017.
Table 4. Transition matrix for changes in the area (km2) of landscape types in the OCWJ during 2000–2017.
Year 2017EXLLOWRELMEDREHHIGEXHTOT
Year 2000
EXL483.7166.1217.590000567.42
LOW1.636.98124.74331.74273.4588.150826.68
REL0070.72108.82262.80144.2023.22609.75
MED001.5636.99119.0199.8310.03267.42
REH0026.0518.4859.85121.2992.01317.68
HIG00112.1367.17164.5571.0227.72442.60
EXH0037.74132.85103.6486.442.78363.44
TOT485.3473.10390.52696.05983.29610.93155.763394.98
Notes: EXL: extremely low; LOW: low; REL: relatively low; MED: medium, REH: relatively high; HIG: high; EXH: extremely high; TOT: total.

Share and Cite

MDPI and ACS Style

Liu, Y.; Liu, Y.; Li, J.; Lu, W.; Wei, X.; Sun, C. Evolution of Landscape Ecological Risk at the Optimal Scale: A Case Study of the Open Coastal Wetlands in Jiangsu, China. Int. J. Environ. Res. Public Health 2018, 15, 1691. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081691

AMA Style

Liu Y, Liu Y, Li J, Lu W, Wei X, Sun C. Evolution of Landscape Ecological Risk at the Optimal Scale: A Case Study of the Open Coastal Wetlands in Jiangsu, China. International Journal of Environmental Research and Public Health. 2018; 15(8):1691. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081691

Chicago/Turabian Style

Liu, Yongchao, Yongxue Liu, Jialin Li, Wanyun Lu, Xianglin Wei, and Chao Sun. 2018. "Evolution of Landscape Ecological Risk at the Optimal Scale: A Case Study of the Open Coastal Wetlands in Jiangsu, China" International Journal of Environmental Research and Public Health 15, no. 8: 1691. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081691

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop