Effects of Urbanization Expansion on Landscape Pattern and Region Ecological Risk in Chinese Coastal City: A Case Study of Yantai City

Applied with remote sensing, GIS, and mathematical statistics, the spatial-temporal evolution characteristics of urbanization expansion of Yantai city from 1974 to 2009 was studied. Based on landscape pattern metrics and ecological risk index, the landscape ecological risk from the landscape pattern dynamics was evaluated. The results showed that the area of urban land increased by 189.77 km2 with average expansion area of 5.42 km2 y−1 from 1974 to 2009. The urbanization intensity index during 2004–2009 was 3.92 times of that during 1974–1990. The land use types of urban land and farmland changed greatly. The changes of landscape pattern metrics for land use patterns indicated that the intensity of human activities had strengthened gradually in study period. The landscape ecological risk pattern of Yantai city shaped half-round rings along the coastline. The ecological risk index decreased with increase of the distance to the coastline. The ratio of high ecological risk to subhigh ecological risk zones in 2009 was 2.23 times of that in 1990. The significant linear relationship of urbanization intensity index and regional ecological risk indicated that the anthropological economic activities were decisive factors for sustainable development of costal ecological environment.


Introduction
It is common that the urbanization expansion is rapid in the world, especially in the developing countries where there are about 70% of the world's largest cities [1][2][3]. The urbanization has boosted rapidly in China since the reform and opening-up in the 1980s. The urbanized area in Yangtze River Delta region increased exponentially during 1979-2008 [4]. According to the report of the National Bureau of Statistics in 2012, the urbanization ratio of China has increased by 1.35% per year from 2002 to 2011. By the end of 2011, the total urban population has been about 690.69 million people and the urbanization ratio was about 51.27%. The human development report 2013 of the United Nations Development Programme predicted that the urban population of China will be 1000 million people and the urbanization ratio will be 70% in 2030. Up to now, more than 40% of population and more than 70% of cities are widely distributed in China's coastal areas [5] where many coastal economic zones are set up and have become the preferred destination for millions of internal migrants and overseas investors [6]. Some studies indicated that the urbanization ratio will increase to 60% in 2020 in the eastern coastal provinces and cities of China [7].
The rapid development of urbanization expansion and heavy burden of human activities have impaired ecosystem functions in coastal areas [8][9][10][11][12]. From the view of landscape ecology, urbanization is the process in which land use/cover landscape changes from natural landscape which is mainly made of water, soil, and vegetation to manmade landscape which is mainly composed of cement, asphalt, chemical materials, and metal [13,14]. Meanwhile, the change of land use/cover landscape pattern could be interpreted as the 2 The Scientific World Journal changes of patch shape, area, quality, and spatial combination [15,16]. Rapid urbanization has been companied by more drastic changes of land use types, disorderly landscape layout, and fragile ecological environment [14,16,17]. Therefore, the problems of landscape pattern change and region ecological risk under urbanization expansion in China have drawn the increasing attention of the ecologists [4][5][6][7][8][9][10]. The technologies of RS and GIS have become important methods to identify the potential eco-environmental influences and guide urban planning management [3,15,[18][19][20][21][22]. The previous researches in China mainly focused on large cities distributed in the Yellow River [23], the Yangtze River [24,25], Pearl River Delta [9,26], and the Ring-Bohai Region [27][28][29]. Few studies have focused on landscape pattern and region ecological risk under urbanization expansion in the small and mediumsized coastal cities [30][31][32]. In the study, Yantai city, which was one of the 14 coastal open economy cities in China, was selected to study the impacts of urbanization expansion on the landscape pattern change and regional ecological risk for the small and medium-sized coastal cities. The objectives of the present study were to reveal (1) dynamic changes and spatial characteristics of landscape pattern and the regional ecological risk and (2) the relationship between urbanization intensity index and regional ecological risk based on the analysis of urbanization expansion in Yantai city.  Figure 1). It covers an area of 13746.50 km 2 with the coastline of 909.00 km. The mountain, mound, plain, and the marsh land are about 36.6%, 39.7%, 20.8%, and 2.9% of total area, respectively. Average annual precipitation is 651.90 mm and mean annual temperature is 11.8 ∘ C. The total population was 6.96 million in 2010 based on the 6th population census data of China. The gross domestic product (GDP) was about 490.68 billion Chinese YUAN according to the government work report of Yantai City in 2011. With the improvement of comprehensive competitiveness, Yantai city has been an important economic zone of the Ring-Bohai Region.

Materials and Methods
Five districts of Yantai city, that is, Zhifu District, Laishan District, Mouping District, Fushan District, and Economic and Technological Development Zone (Development Zone), were selected for the study (Figure 1). The study area (121 ∘ 15 -121 ∘ 56 E, 37 ∘ 04 -37 ∘ 38 N) covers an area of 2721.04 km 2 and the total population is about 2.22 million in 2010. With rapid economic development and enormous population growth, large area agricultural land has been occupied by nonagricultural lands and urban landscape pattern resulting from urbanization expansion in the past years.  The topographic maps of Yantai city, geometric correction of Landsat TM remote sensing images had been conducted by using ArcGIS 9.3. The error was controlled within 0.5 pixels. All images were classified by combing supervised classification with visual interpretation and the interpretation results were resampled in spatial resolution (30 m).

Data Sources and
According to the Chinese national standard of Status of Land Use Classification (GB/T 21020-2007), the land use types of the study area were divided into urban land, farmland, orchard, forestland, grassland, water body, other construction land (industrial land, mine land, and port land), and other object land (unused land, barren land). Based on land use interpretation results which were corrected by field survey, the classification error rate was lower than 5% of the land use actual change.

Urbanization Expansion.
The urbanization intensity index ( ue ) was used to analyze the urbanization expansion intensity from 1974 to 2009. The equation was as follows: The Scientific World Journal 3 To reflect the complexity of self-similarity of a patch Shannon's Diversity Index (SHDI) To what extent several principal landscape types control whole landscape To express the agglomeration degree among different landscape types : patch type; : number of patch; : the total area; : the area of patch ; : patch perimeter; : a constant, as to grid map, equals to 4; ℎ : contiguity number among patch ℎ and patch ; : the area ratio of class ; : quantity of region landscape types.
where ue is urbanization intensity index, Δ is expansion area of urban land during a certain period, TLA is the total area of the research area, and Δ is time span of a certain duration.

Landscape Pattern Metrics.
The six landscape metrics of Largest Patch Index (LPI), Fractal Dimension Index (FD), Shannon's Diversity Index (SHDI), Shannon's Evenness Index (SHEI), Dominance Index (Dominance), and Contagion Index (CONTAG) were selected for landscape pattern analysis. As a measure of landscape heterogeneity, SHDI is especially sensitive to the nonbalanced distribution of all patch types in the landscape. SHEI is applied to indicate the diversity of different landscapes or a certain landscape in different periods, which results in maximum evenness. As such, evenness could be the complement of dominance. Contagion Index is related to edge density and can measure the clumping trends of the patches. The higher the CONTAG is, the more homogeneous and contiguous the spatial pattern is. The Fragstats 3.3 was adopted to calculate and analyze landscape pattern metrics. Expressions and ecological significance of metrics in the present study are given in Table 1 [33][34][35][36][37][38].

Ecological Risk Assessment.
The ecological risk index (ERI) [39] (2) was adopted to assess the spatial difference of comprehensive ecological risk in the various land use types under urbanization expansion: where ERI is ecological risk index; is land use type; is area of land ; is the total area of research area; which is the comprehensive ecological risk intensity reflected by is determined by the method of the analytic hierarchy process (AHP) for each land use type. The values of for urban land, farmland, orchard, forestland, grassland, water body, other construction land, and other object land are 0.29, 0.13, 0.11, 0.03, 0.05, 0.04, 0.20, and 0.07, respectively.

Results and Discussion
The urbanization expansion and landscape pattern change were influenced by many factors, such as geographical location, population, economic policy, resources, and transportation [25,40]. The coastal location advantages played an important role in urbanization expansion [8,9]. The spatial distribution of urban land expanded along the coastline (Figure 2). Influenced by port economy, the distribution of urban land in the study area from 1974 to 1990 was formed with the spatial construction of "one center, " which was similar with the urbanization process of many coastal cities which depended on port economy in the world [26,41] (Table 2) and about 237 km 2 farmland lost during this period ( Table 3). The fast speed urbanization expansion started since 2000. From 2000 to 2009, the area of urban land added 92.20 km 2 and urbanization intensity index was 0.38. Meanwhile, the urban space had transformed to axial expansion that developed from compact nutty to dispersed groups (Figure 2), which was similar with the axial expansion of Shanghai city between 1979 and 2000 [25]. Furthermore, with the ideology changes of urban construction policy [42] and transportation construction [25], the urban land of Yantai city has expanded to inner regions with the T-shape along the coastline (Figure 2), then has gradually aggregated from dispersed. In recent years, the agglomeration growth pattern of urbanization expansion became more and more obvious.

Analysis of Land Use/Cover and Landscape Pattern
Changes. The dominant land use in study area in 1990 was farmland (48.45%), forestland (24.17%), and grassland (7.91%), the sum of which was more than 80% of the total area ( Table 3). The urban land covered only 82.24 km 2 (Figure 3), no more than 3% of the total area. In 2000, farmland and grassland had reduced to 1173.11 km 2 and 163.72 km 2 , respectively, while the areas of urban land, orchard, forestland, water body, and other construction land increased greatly. Compared to 1990, the orchard, urban land, and forestland increased by 92.44 km 2 , 44.76 km 2 , and 40.03 km 2 , respectively (Table 3)    of farmland and green land [14,31,40,43]. Farmland was the largest land use type in the study area (Table 3). During the period 1990-2009, the area of urban land increased greatly (about 137 km 2 ) and farmland decreased from 1318.33 km 2 in 1990 to 1081.05 km 2 in 2009 ( Figure 3 and Table 3). The result was similar with another Chinese coastal city, Longkou city [31]. Due to development of fruit economy, part of farmland was transformed into orchard ( Table 3). The area of forestland increased from 1990 to 2000 because of the implementation of "returning farmland to forestland and grassland" project, while decreasing from 2000 to 2009 due to human activities increase and ecology degradation.
In contrast, the forestland, water body, and other object land tented to decrease and most of lost areas had been transformed into maritime industrial land and tourism land. The calculation results of landscape metrics which indicate landscape pattern change showed that the LPI decreased from 23.32 in 1990 to 20.35 in 2009 (Table 4), indicating that the dominant landscape of farmland had declined. Due to the strengthen of human interferences, the study area showed more fragmented landscapes, fewer large patches, which caused the Dominance and CONTAG decrease and FD increase during 1990-2009 (Table 4). Meanwhile, the urban land rapidly expanded from the urban center to the surrounding areas and then aggregated to a new large patch. The increase of SHDI and SHEI from 1990 to 2009 indicated that the landscape shapes in the study area had become more and more complex because of human activities [8].
3.3. Analysis of Regional Ecological Risk. Using the Natural Breaks classification method [21], the ecological risk was grouped into five grades, that is, low ecological risk (I), The Scientific World Journal 7  sublow ecological risk (II), moderate ecological risk (III), subhigh ecological risk, and (IV) and high ecological risk (V), according to the ecological risk index. The results showed that the high and subhigh ecological risk areas are mainly distributed in the urban land, while low and sublow ecological risk areas appeared in forestland and grassland with high vegetation coverage (Figure 4). In addition, the ecological risk level decreased with the distance to the coastline increase. From 1990 to 2009, the moderate ecological risk level was major, and the proportions of subhigh and sublow ecological risk changed obviously in study area ( Figure 4 and Table 5). The areas of moderate and sublow ecological risk decreased to 145.12 km 2 and 266.47 km 2 from 1990 to 2009, respectively. The low ecological risk area reduced to 4.5 km 2 . On the contrary, the area of high ecological risk increased from 37.79 km 2 in 1990 to 61.47 km 2 in 2000 and expanded rapidly to 104.18 km 2 in 2009. The high ecological risk region mainly distributed along the coastline of Zhifu District and expanded westward to Fushan District and Development Zone (Figure 4). The area of subhigh ecological risk in 2009 was 2 times higher than that in 1990 and expanded to inner region from inshore region (Table 5 and Figure 4). The ecological security in the coastal zone was largely influenced by urbanization expansion and land use change [8,16,32,44]. As the sensitive area in land-ocean interactions zone, Yantai city is located in the integrated high ecological risk area of the Ring-Bohai Region [16,45]. In this study, the area of high and subhigh ecological risk region increased rapidly and extended eastward and westward along the coastline with layered distribution and the area of moderate and sublow ecological risk region reduced under the development of urbanization expansion (Figure 4). Such changes, especially the increasing trend of high ecological risk, were also observed in the coastal economic developing zone of Jinzhou Bay in Liaoning Province and Yancheng coastal area in Jiangsu Province [29,45]. Previous studies predicted that more and more farmland and coastal water of Yantai city will be lost in the next 40 years due to the acceleration of urbanization process; a series of coastal development activities and reclamation projects would enhance the regional ecological risk [32]. The significant positive linear relationship between the urbanization intensity index and regional ecological risk index both in 1990-2000 and 2000-2009 with the determination coefficient ( 2 ) of 0.831 and 0.951, respectively, ( Figure 5) was observed in the study, which obviously supported these prediction results. The results of this study suggested that the scientific urban development plan and rational coastal management measures should be implemented in order to reduce the ecological risk brought from urbanization expansion.

Conclusions
With economic development and population growth, the process of urbanization of Yantai city had accelerated obviously and the area of urban land increased quickly. The area of urban land increased about 137 km 2 and about 237 km 2 farmland lost during the period of 1990-2009. The distribution of urban land showed a T-shape along the coastline extent to inner regions. The area of high ecological risk increased rapidly and the moderate and sublow ecological 8 The Scientific World Journal risk decreased from 1990 to 2009. The high and subhigh ecological risk areas mainly distributed in the urban land and the ecological risk level decreased with the distance to the coastline increase. The significant positive relation between urbanization intensity and regional ecological risk was found in the study.