Wetland Transitions and Protection under Rapid Urban Expansion: a Case Study of Pearl River Estuary, China

Understanding wetland changes under urbanization is important for wetland management. In this study, net transition intensity (NTI) and total transition intensity (TTI) are presented to characterize wetland transitions based on spatial data obtained from Landsat satellite images of Pearl River estuary in South China. NTI is commonly used to represent changes in absolute amounts for each class of wetland, while TTI reflects the internal transition activities and amounts. The third index, the urbanization intensity index (UII), is used to investigate the intensity and velocity of urban land expansion at the same time periods. The results show that one-third of the total wetland area was lost from 1979 to 2009 in the study area and seven types of estuarine wetlands were degraded. The basic pattern of wetland transition is from natural wetlands to constructed wetlands and then to urban lands. Intertidal mud and sand and paddy fields were the major natural and constructed wetlands, respectively, transferred to urban lands. The TTI value was generally greater than the NTI value for these wetlands. TTI >> NTI is an important indicator for wetland transitions under rapid urban expansion in the Pearl River estuary. Based on the integration of the two indices (NTI and TTI), a protection and management plan framework for the Pearl River estuary wetlands is proposed. This plan emphasizes the key and important zones and their different features, and includes actions that can be implemented in and around natural and constructed wetlands.


Introduction
Global urbanization significantly affects estuary wetlands by altering, destroying or fragmenting natural habitats, which alters hydrological and geomorphological processes and degrades wetland functions [1][2][3][4].Generally, the location, size, and shape of wetlands significantly influence ecosystem functions.An increase in wetland area enhances ecosystem functions such as flood control, pollution filtration, carbon entrapment, and fishing, and provides migratory pathways for native species by improving hydrological and biological connectivity, which determines the total habitable area [5].However, a decrease in wetland area reduces the available habitat for native fauna and flora [6] and ecosystem functions [7].During the past few decades, a significant area of estuarine wetlands has disappeared in some large river deltas due to rapid urbanization [8].The remaining wetlands are suffering from population pressure, water shortage, and pollution, as well as biological invasions, and are becoming increasingly vulnerable [9][10][11].
The changes to estuarine wetlands include not only direct occupation and loss, but also construction or increase of new wetlands.Those wetlands occupied or lost are completely replaced by agricultural use, aquaculture lands, infrastructure, and tidal control structures; the remaining wetlands are bounded by development to varying extents [12][13][14].These conditions reduce the hydrological and biological connectivity among wetlands, such that the wetlands lose their resilience when facing different stressors [15].Newly built wetlands, such as salt fields, shrimp ponds, reservoirs, and restored sites, can partially offset the negative effects arising from the decrease in wetland area.
Different methods, such as rate of change of total wetland area [16][17][18], transition matrices for the process of wetland loss [9,19,20], and landscape pattern indices for wetland development [21][22][23], have been used to measure and evaluate wetland changes with absolute amounts by comparing their state at two time points.However, because the differences between newly increased wetlands and well-preserved ones, as well as between the same wetland types, are significant in structure and function due to their development period and location [13], understanding the internal transitions together with analyses of wetland changes plays an important role in conservation and restoration of wetlands [24].Nevertheless, fewer indices have been studied to understand wetland interchange processes that occur at various stages, especially in the urban expansion area, during wetlands transferring to other land uses or vice versa.
The main objectives of this study were: (1) to develop a method to detect the characteristics of wetland transitions under rapid urban expansion in the Pearl River estuary; and (2) to quantify the spatial relationship between wetlands and urban lands for future restoration and conservation of estuarine wetlands.

Study Area
The Pearl River estuary (21 ˝34 1 N to 23 ˝56 1 N, 111 ˝58 1 E to 114 ˝38 1 E) with an area of 31,831 km 2 (Figure 1) is located on the southern coast of China.The Pearl River enters into the southern China Sea through the estuary by means of eight river branches, which results in different types of large wetlands.The climate of this region is warm and moist with a mean annual temperature of 22.1 ˝C.The weather is dominated by the northeast monsoon and the mean rainfall is low (30-40 mm month ´1) from October to April.From May to September, the southwest monsoon brings up large amounts of rainfall (300-400 mm month ´1) from lower latitudes [16,25].Over the past three decades, the Pearl River estuary wetlands have suffered from rapid industrialization and urbanization.Seven major large cities, including 35 districts, have been formed [26,27].
The changes to estuarine wetlands include not only direct occupation and loss, but also construction or increase of new wetlands.Those wetlands occupied or lost are completely replaced by agricultural use, aquaculture lands, infrastructure, and tidal control structures; the remaining wetlands are bounded by development to varying extents [12][13][14].These conditions reduce the hydrological and biological connectivity among wetlands, such that the wetlands lose their resilience when facing different stressors [15].Newly built wetlands, such as salt fields, shrimp ponds, reservoirs, and restored sites, can partially offset the negative effects arising from the decrease in wetland area.
Different methods, such as rate of change of total wetland area [16][17][18], transition matrices for the process of wetland loss [9,19,20], and landscape pattern indices for wetland development [21][22][23], have been used to measure and evaluate wetland changes with absolute amounts by comparing their state at two time points.However, because the differences between newly increased wetlands and well-preserved ones, as well as between the same wetland types, are significant in structure and function due to their development period and location [13], understanding the internal transitions together with analyses of wetland changes plays an important role in conservation and restoration of wetlands [24].Nevertheless, fewer indices have been studied to understand wetland interchange processes that occur at various stages, especially in the urban expansion area, during wetlands transferring to other land uses or vice versa.
The main objectives of this study were: (1) to develop a method to detect the characteristics of wetland transitions under rapid urban expansion in the Pearl River estuary; and (2) to quantify the spatial relationship between wetlands and urban lands for future restoration and conservation of estuarine wetlands.

Study Area
The Pearl River estuary (21°34′ N to 23°56′ N, 111°58′ E to 114°38′ E) with an area of 31,831 km 2 (Figure 1) is located on the southern coast of China.The Pearl River enters into the southern China Sea through the estuary by means of eight river branches, which results in different types of large wetlands.The climate of this region is warm and moist with a mean annual temperature of 22.1 °C.The weather is dominated by the northeast monsoon and the mean rainfall is low (30-40 mm month −1 ) from October to April.From May to September, the southwest monsoon brings up large amounts of rainfall (300-400 mm month −1 ) from lower latitudes [16,25].Over the past three decades, the Pearl River estuary wetlands have suffered from rapid industrialization and urbanization.Seven major large cities, including 35 districts, have been formed [26,27].

Data Sources
In this study, wetland information was extracted from the satellite data of three Landsat instruments (Multispectral Scanner [MSS], Thematic Mapper [TM], and Enhanced Thematic Mapper plus [ETM + ]).Forty-nine scenes of Landsat satellite data were obtained from the Center for Earth Observation and Digital Earth, Chinese Academy of Science and the Global Land Cover Facility at the University of Maryland [28].Autumn and winter images, with less cloud and rain, were chosen to better differentiate the spectral information.The characteristics of the image data are presented in Table 1.A contour of ´6 m obtained from Chinese digital elevation model (1985) was used to mark off the wetland class of shallow marine water (NW1) based on the Ramsar Convention.The other used data mainly include: (1) topographic maps (1/100,000 scaled); (2) detailed vegetation maps (2005, Forestry and Water Affairs Minister in Guangzhou); (3) ground data (land survey and fieldwork); and (4) socioeconomic data (statistical year books).

Satellite Images Pre-Processing
All satellite images were geometrically corrected based on topographic map by linear polynomial method with 8 control points [29].The total root mean square (RMS) error for this study was 4.5 m approximately, which is less than half the size of the pixel (15 m) [30,31].The images were georeferenced using the Universal Transverse Mercator projection (UTM-49N, WGS-84) and resampled to 30 m spatial resolution [32].The processing of remote-sensing satellite was carried out using the software package of Erdas Image V9.3 [33].
Related literature and image data were gained for the sake of better knowledge of the study region.Meanwhile, a fieldwork was developed with interviews, questionnaire and photography [32,34].From August to October in 2009, field transects were carried out in order to verify the classification results of classification, and the hand-held GPS was used to record all the locations of interesting points and field transects.

First Class Level Second Class Level Meaning
Natural wetlands (NW) Shallow marine water (NW1) Ocean area, intertidal and shallow (<6 m at low tide) sub-tidal soft sediment areas, including sea bays and straits, including water of estuaries and estuarine systems of deltas.
Intertidal mud and sand (NW3) Sand, shingle or pebble shores; sand bars, spits and sandy islets, including dune systems and humid dune slacks; salt marshes, salt meadows, saltings, raised salt marshes; tidal brackish and freshwater marshes.
Aquaculture pond (CW2) Land used to feed fish and shrimp, including fishponds and aquafarms.
Paddy fields (CW3) Irrigation channels and rice fields.

Urban land (UL)
Urban land (UL1) Land covered by buildings, residential areas, and other man-made structures.
Construction land (UL2) Land, outside the town, is independent of mines, large industrial zones, oil fields, salt fields, quarry lands, traffic roads, airports, etc.
Other land (OL) Forest (OL1) Land dominated by trees, with a percent canopy cover of >10% Grassland (OL2) Land dominated by herbaceous types, or tree and shrub cover is <10%.

Dry farmland (OL3)
Land covered with temporary crops followed by harvest and a bare soil period (e.g., single and multiple cropping systems).Note that perennial woody crops will be classified as the appropriate forest or shrub land cover type.

Classification of Satellite Images
The classification criteria for the main types of wetlands and other land uses in the Pearl River estuary according to the Ramsar Convention on wetland classification are shown in Table 2.A multi-step classification procedure using the masked images was employed to characterize the different wetlands.At the beginning, the coastline was directly digitized from the satellite data, and NW1 wetlands were determined using digitized coastline and the ´6 m contour.A supervised maximum likelihood classification was then performed [32,33].
After classification, the overall accuracy and a Kappa analysis could be used to perform a classification accuracy assessment based on an error matrix analysis [35].Independently classified images were compared with each other to determine the changes of landuse types.Accuracy levels of more than 80% were considered enough for reliable classification of landuse [32,36].In this study, the Kappa parameter was computed as an accuracy assessment for the different years, and were given as 0.85 (1979), 0.84 (1986), 0.85 (1990), 0.89 (1995), 0.85 (2000), 0.89 (2005), and 0.85 (2009).

Monitoring Indices of Wetland Transition and Urban Expansion
Detection and monitoring of wetland transitions included the use of multi-date images to evaluate the differences in wetland patterns [17,20].Overlay analysis was implemented for continuous years to examine the inherent dynamics of wetlands during different periods.Two images on different dates were labeled and classified independently.Then, the transferred areas were extracted by the comparison of classification results.To understand the dynamics of estuarine wetland transitions, the mutual transfer extent among land uses was obtained from the overlay analysis and transition matrices using the spatial analysis module in ArcGIS 10.0.
Two indices, NTI i,tÑt+n (Net Transition Intensity) and TTI i,tÑt+n (Total Transition Intensity), were developed to evaluate wetland transitions.NTI i,tÑt+n is the index of net transition intensity of wetlands during the time span t, t + n; TTI i,tÑt+n is the total transition intensity of wetlands in the transfer process during the time span t, t + n.They are defined as: where wa in_i,t is the wetland area of the ith wetland class transferred from other land use during the time span t, t + n; and wa out_i,tÑt+n is the area of other land use that is transferred from the ith wetland class during the time span t, t + n.WA i,t+n and WA i,t are the area of wetlands in the year t + n and t, respectively.N is the total number of wetland classes.Generally speaking, TTI can reflect the internal transition activities and total transition amounts for wetlands, but NTI is just used to represent changes in absolute amounts.
To examine the spatial relationship between wetland transitions and urban expansion in the Pearl River estuary, the urbanization intensity index (UII) [37] was used to characterize the intensity and velocity of urban land expansion at different time periods, which is given as: where UII i,tÑt+n is the index of urban expansion intensity for spatial unit i in the period of t; t + n, ULA i,t+n and ULA i,t are the urban land use area in the year t + n and t, respectively; and TLA i is the total area of spatial unit i.

Patterns of Wetland Transitions
Seven wetland classes and other land use (Figure 2) were obtained from the interpretation of Landsat images.Mutual transfer quantities among land uses were evaluated using the overlay analysis and transition matrices (Figure 3).One-third of the total wetlands were lost in the past 30 years in the Pearl River estuary (decreased from 16,986 km 2     It was found that NW1 and CW3 were the main natural and constructed wetland classes, respectively, in the study area.About 2% of NW1 was transferred to NW3, CW3, and CW2 from 1979 to 1995; 0.5% of NW1 was transferred to NW3; and 0.1% of NW1 was transferred to CW3 and UL.NW2 gradually decreased and was transferred into CW3 and UL1 over the past 30 years.Both NW3 and NW4 were closely related to NW1, and their area was transferred from NW1. NW3 was mainly transferred to NW4, CW2, and CW3, and NW4 was transferred to CW2 and CW3.The most extensive and intensive transitions of NW1, NW3, and NW4 occurred along the shoreline of Huangmao Ocean and the western shoreline of Lingding Ocean. Figure 4a,b shows two specific locations from these shorelines.The three classes of constructed wetlands were mostly transferred into other land and urban land (Figure 4d).The transfer intensity of CW3 was the largest among the seven classes of wetlands.Three main modes of transition patterns were determined:  It was found that NW1 and CW3 were the main natural and constructed wetland classes, respectively, in the study area.About 2% of NW1 was transferred to NW3, CW3, and CW2 from 1979 to 1995; 0.5% of NW1 was transferred to NW3; and 0.1% of NW1 was transferred to CW3 and UL.NW2 gradually decreased and was transferred into CW3 and UL1 over the past 30 years.Both NW3 and NW4 were closely related to NW1, and their area was transferred from NW1. NW3 was mainly transferred to NW4, CW2, and CW3, and NW4 was transferred to CW2 and CW3.The most extensive and intensive transitions of NW1, NW3, and NW4 occurred along the shoreline of Huangmao Ocean and the western shoreline of Lingding Ocean. Figure 4a,b shows two specific locations from these shorelines.The three classes of constructed wetlands were mostly transferred into other land and urban land (Figure 4d).The transfer intensity of CW3 was the largest among the seven classes of wetlands.

Amounts of Wetland Transitions (NTI and TTI)
Total and absolute transition amounts for wetlands were calculated (Figure 5).If NTI < 0, the absolute amount for wetland was decreased.Otherwise, the absolute amount increased.Two net transition patterns among the seven classes of wetlands based on NTI values were determined.One was NTI < 0 in each stage within a five-year interval, such as NW1, NW2, CW2, and CW3, where their area decreased continually during the 30 years.The other was NTI > 0 at certain stages, such as NW3, NW4, and CW1, which showed a fluctuating trend over the six stages.The negative NTI values of NW, CW, and total wetland from 1990 to 1995 were the largest out of the six stages, indicating the maximum net area loss in the past 30 years.TTI values for most wetlands were much bigger than their NTI absolute values.For example, TTI of NW3 was four times its NTI at the same stage.This indicates that the internal transition, which includes both decreases and increases in areas, for most wetlands was frequent and complex.The TTI values of NW1 and NW2 were equal to their NTI absolute values at the same stage, indicating only a decrease in area occurred in NW1 and NW2.

Amounts of Wetland Transitions (NTI and TTI)
Total and absolute transition amounts for wetlands were calculated (Figure 5).If NTI < 0, the absolute amount for wetland was decreased.Otherwise, the absolute amount increased.Two net transition patterns among the seven classes of wetlands based on NTI values were determined.One was NTI < 0 in each stage within a five-year interval, such as NW1, NW2, CW2, and CW3, where their area decreased continually during the 30 years.The other was NTI > 0 at certain stages, such as NW3, NW4, and CW1, which showed a fluctuating trend over the six stages.The negative NTI values of NW, CW, and total wetland from 1990 to 1995 were the largest out of the six stages, indicating the maximum net area loss in the past 30 years.TTI values for most wetlands were much bigger than their NTI absolute values.For example, TTI of NW3 was four times its NTI at the same stage.This indicates that the internal transition, which includes both decreases and increases in areas, for most wetlands was frequent and complex.The TTI values of NW1 and NW2 were equal to their NTI absolute values at the same stage, indicating only a decrease in area occurred in NW1 and NW2.
In Figure 6a,b, the relationships between wetland transitions and urban land expansions are obvious in the Pearl River estuary based on the integrated TTI and NTI values and the UII values.The highest intensities of urban expansion were observed in the center of Guangzhou City, the whole of Shenzhen City, and the northeast region of Zhuhai City, where no high rapid wetland transition zones were found.The highest intensity of wetland transitions only appeared in rapid urban expansion zones, such as Nanhai District, Shancheng District, and Nanshan District, where a large number of CW3 wetlands exist.According to the hierarchy maps of wetland transitions and urban expansions, a conservation plan to restore and maintain a variety of wetlands was devised.This plan divides the whole study area into five zones: (1) wetland restoration zones; (2) important wetland protection zones; (3) general wetland maintenance zones; (4) optimized urban development zones; and (5) natural wetland restoration zones.The scope and conservation measures for each zone are given in Table 3 and Figure 6c.The integration of transition indices and a wetland conservation plan will facilitate the development of a long-term assessment and management system for estuarine wetlands in the Pearl River delta.

Characteristics of Wetland Transitions under Rapid Urban Expansion
Transition intensity is an important indicator for assessing land-use changes [19,20].In this study, two transition indices were presented to represent the wetland transition intensity between various wetlands and urban lands.NTI is commonly used to represent changes in absolute amounts, while TTI reflects the internal transition activities for each class of wetland.When wetlands are in their original state or at the beginning of a transition, the values of NTI and TTI are the same due to a single change (e.g., NW1 to NW2).When wetlands are in the middle stages of the transition process, bi-directional conversions might occur [17,22].Generally, the TTI value is much greater than the NTI value.For example, the NTI value of CW3 is <5% from 1990 to 1995, but TTI values of CW3 exceed 10% due to the change in transfer patterns.This phenomenon is an important hint of wetland degradation in the rapid urbanized area.
Different transition intensity indices can be explained as the degree of 'reclamation activities' in a specific geographical region and different social communities.A large number of NW1 wetlands were occupied by newly built intertidal mud due to reclamation activities.The intertidal mud was then quickly occupied by CW3.During this process, the loss and newly built intertidal mud existed at the same stage, and only 10% of old intertidal mud was left at its original location in the Pearl River estuary.In the case of wetlands located on the reclaimed landscape, specific criteria have not yet been developed.The combined analysis of the two indices could be used to explain the characteristics of the detailed transitions and reclamation activities of wetlands.

Implications for Wetland Management
The management of natural and constructed wetlands usually focuses on one or several objectives, such as wetland functions and service values.In fact, it is difficult to operate wetlands with fragmented habitats induced by urbanization.Coastal zone management is considered a better way to strengthen the conservation of wetlands [38].This study proposes a conservation plan (Table 3 and Figure 6c) that integrates the indices of wetland transitions and urban expansion intensities to facilitate the management of the degraded and constructed wetlands.This plan involves activities that can be conducted in and around both natural and constructed wetlands.The most important measure is controlling the expansion of urban land or limiting the occupation of the natural wetlands.This plan depends on the government actions and the participation of local people.The Ramsar Convention has been playing a leading role, and guidelines on the participatory management of wetlands have been adopted [39,40].
In the study area, the wetland transition activities were very intense in the first 15 years (1979-1995) due purely to the pursuit of economic development.In the second 15 years (1995-2009), wetland conservation was considered seriously by the local government, and many activities affecting natural wetlands were forbidden by the government regulations.Some wetland reserves were built along Lingding Ocean, e.g., Futian Mangrove Nature Reserve (2000) and Nanshan Wetland Nature Reserve (2007).The aforementioned two transition periods indicate the importance of rational management policies and measures for wetlands.Some individual and governmental activities, such as digging ditches and impoundments, stream channelization, and placement of impervious surfaces, which lead to the loss and degradation of wetlands, have not been completely controlled so far.Further measures are required for conservation of wetlands in the estuary wetlands.

Conclusions
In this study, the spatial characteristics of estuarine wetland transitions under rapid urban expansion from 1979 to 2009 in the Pearl River estuary were investigated.Two indices, NTI and TTI, were calculated to characterize the wetland transitions at the regional scale.The combined analysis of NTI and TTI provided useful information on wetland transitions and the patterns of each wetland class.According to the hierarchy maps of wetland transitions and urban expansion intensities, a conservation plan was devised to restore and maintain the various wetlands.More wetland transition indices are required to forecast future trends and facilitate decision-making for wetland conservation.

Figure 1 .
Figure 1.Layout of the study region (Seven large cities and two major river branches in the Pearl River estuary, i.e., Guangzhou City, Foshan City, Zhongshan City, Dongguan City, Jiangmen City, Shenzhen City, Zhuhai City and Huangmao Ocean, Lingding Ocean).

Figure 1 .
Figure 1.Layout of the study region (Seven large cities and two major river branches in the Pearl River estuary, i.e., Guangzhou City, Foshan City, Zhongshan City, Dongguan City, Jiangmen City, Shenzhen City, Zhuhai City and Huangmao Ocean, Lingding Ocean).
in 1979 to 10,465 km 2 in 2009).The decreased area for wetlands from 1979 to 1995 and from 1995 to 2009 is 4665 km 2 and 1851 km 2 , respectively.Sustainability 2016, 8, 471 7 of 15 decreased area for wetlands from 1979 to 1995 and from 1995 to 2009 is 4665 km 2 and 1851 km 2 , respectively.

Figure 3 .
Figure 3.The transition process of various wetland classes from 1979 to 2009.

Figure 3 .
Figure 3.The transition process of various wetland classes from 1979 to 2009.

Sustainability 2016, 8 , 471 9 of 15 →
CW3 → UL (Figure4d).The transition processes of all wetland classes showed that the dominant transition pattern was NW → CW → UL or OL in the six stages from 1979 to 2009 (Figure3).

Figure 5 .
Figure 5. NTI (net transition intensity) and TTI (total transition intensity) of the main wetland classes in six periods (TW: total wetlands).

Table 1 .
List of the satellite data used.

Table 2 .
The classification of main wetlands and other land uses.

Table 3 .
Protection plan for wetlands of the Pearl River estuary.