Research on the Relationship between Urban Development Intensity and Eco-environmental Stresses in Bohai Rim Coastal Area, China

To realize sustainable urban development that minimizes environmental impacts, the relationship between urban development intensity and eco-environmental stresses should be clearly revealed. This paper focused on the Bohai Rim coastal area, where cities have experienced significant development in the last decade. An index system was developed to quantify the comprehensive urban development intensity (CDI) and comprehensive eco-environment stresses (CES). Remote sensing imagery and statistical data were used to provide indices for CDI and CES. Spatiotemporal analysis was carried out on the correlation between the two indices. The coupling between the CDI and CES was then investigated to explore the urban development characteristics of each city in the study area, its development level, and the trend of urban development. Results showed that human activities surrounding urban development were partly dependent on the use of ecological resources to a certain degree, and that the degree of dependence increased with year. To promote a sustainable level of urban development, the government should focus on not only the high development intensity, but also the high quality of the eco-environment. Dalian was a good model of how to achieve a balance between the two.


Introduction
The urbanization of China has accelerated greatly in recent decades, resulting in intense urban development that has negatively impacted ecological systems [1][2][3][4][5].The coordination of urban development activities, ecological and environmental needs, and influences on urban development level is a challenge for urban sustainable development and has been the focus of considerable research [6][7][8][9].Therefore, a great many scholars have studied the stability of urban eco-environmental systems from the perspective of the relationship between human development intensity and ecological environment [10][11][12][13].Liu et al. [14] used GIS spatial analysis and regression models to study the effect of development intensity factors and validated the results by ROC (Relative Operating Characteristics).Wang et al. [5] provided a new method to evaluate urban land development intensity from the perspective of land property by building a dynamic model influenced by development intensity and using the attribute and spatial characteristics of adjacent land.Gong et al. [15] studied the trend of land development intensity based on back-propagation neural network modeling, utilizing remote sensing image data of four phases and the related statistical data.Cropper and Griffiths [16] and Nyakaana et al. [17] studied the relationship between population and environment in the process of urbanization.After reading much literature about the interplay between transport, land use and environment, and showcasing new empirical approaches and investigations from different parts of the world, Yigitcanlar and Kamruzzaman [18] highlighted the importance of the relationship between transport, land use, and the environment.Coxhead [19] discussed the effects of economic growth, structural change, and industrialization on the environment, which identified the trends and phenomena likely to be associated with changes in the valuation and use of environment and natural resources assets.The above research mostly used land development intensity, urbanization development, or economic development to represent human development and construction activities.Moreover, Grossman and Krueger [20] proposed the hypothesis that economic growth and environmental damage present an inverted-U curve relationship, so-called the Environmental Kuznets Curve (EKC), and was restated by them in 1995 [21,22].From then on, many scholars have paid attention to it [21,23,24].However, the study of development intensity should not only be based on the perspective of urban development, but also on the perspective of harmony between society and nature to study the effects of human development on the ecological environment.In order to realize urban sustainable development, urban development should not destroy the eco-environment, and urban development and the protection of the ecological environment must be kept in a harmonious and orderly relationship, keeping the system in stable development [25,26].Further research is, therefore, needed to comprehensively evaluate the development level of a city or a region and its broader relationships and impacts.Coupling analysis between economic development and the environment were carried out in some cities in China, like Anhui [27] from 1990 to 2008, Xi'an [28] from 2001 to 2010, and Urumqi [29] for the past 20 years.These researchers mainly focused on the relationship between economic development and ecological environment, which is of important reference value in this research field.However, to comprehensively access the relationship between urban development intensity and ecological environment, more factors like population, transportation, and so on, should be taken into consideration.What is more, most of the research was based on time series, and more spatiotemporal comparative research between a group of cities in a region should be carried out.
In this context, the object of this content is to conduct a spatiotemporal contrastive analysis between the development intensity and eco-environmental stresses among the city of Bohai Rim coastal area and give suggestions for urban management and planning.

Study Area
The Bohai Sea (35 ˝42 1 N-42 ˝08 1 N and 115 ˝41 1 E-123 ˝32 1 E) is an inland sea surrounded by the Liaodong Peninsula, Shandong Peninsula and the North China Plain and is the major sea outlet of Northeast, North, and Northwest China (Figure 1).The region experiences four distinct seasons, but belongs to the humid and semi-humid monsoon climate zones, which differs from other coastal areas [30].The Bohai coastal area is about 125,588 km 2 with a shoreline of about 2647 km and a population of 55.69 million.In this paper, the study area of Bohai Rim coastal area (BRC) includes the new Binhai area of Tianjin (BH), Tangshan (TS), Qinhuangdao (QHD), Cangzhou (CZ) in the west, Dalian (DL), Jinzhou (JZ), Yingkou (YK), Panjin (PJ), Huludao (HLD) in the northern peninsula, and Dongying (DY), Yantai (YT), Weifang (WF), Binzhou (BZ) in the southern peninsula, for a total of 13 cities or districts.The BRC has become one of the key developing areas of the whole country since it was established as an economic region by the Chinese government in 1992.At the same time, the environment has faced increasingly serious problems of degradation.

Construction of Index System
Two indices were defined in this study, the Comprehensive Development Intensity (CDI) and the Comprehensive Eco-environmental Stresses (CES), in order to quantify the comprehensive information of development intensity and eco-environment stresses (Table 1).

Construction of Index System
Two indices were defined in this study, the Comprehensive Development Intensity (CDI) and the Comprehensive Eco-environmental Stresses (CES), in order to quantify the comprehensive information of development intensity and eco-environment stresses (Table 1).The CDI index shows the comprehensive development degree of different development and construction activities, which generally reflects the natural resources utility state, economic conditions, and social situation.Sub-indices are selected to support this index.Land resource is the basic natural resource upon which human development activities took place and different types of land use structure may show different land use degrees, so the land-use comprehensive degree (LCD) index [31] is selected to show the development intensity of land resources.From an economic point of view, economic activity intensity (EAI) is defined, which shows the average gross domestic product (GDP) produced on a per unit area.Urbanization shows another important aspect of human development activities.It can be mainly reflected on population, economy, urban expansion, and so on.Therefore, urbanization intensity (UI) is used to show the situation of urbanization of all these cities [32], including three sub-indices, (1) population urbanization ratio (PUR), showing the population growth or transformation; (2) economic urban ratio (EUR), showing the ratio of secondary industry and tertiary industry in total GDP; and (3) land urbanization ratio (LUR), the ratio of urban built up area to that of the whole area transportation construction is one of the most important work in human development activities, so transportation intensity (TI) is used to show the development situation of transportation.As sea reclamation is one of the specific activities in the coastal areas, sea reclamation intensity (SRI) is used to describe the intensity of sea reclamation in a time section [33].

CES and Its Sub-Indices
The CES index represents the comprehensive pressure on urban ecological environment caused by population density and pollutant emission [34].CES is calculated based on a series of sub-indices.First the population density (PD) is used to show the stress that population exerted on economic environment.Generally speaking, the higher the density of population, the heavier the stress is on the environment at the present situation of the study area.Pollutant emissions also produce great pressure on the ecological environment.There, air pollution (AP), water pollution (WP), and solid waste pollution (SWP) indices were selected as sub-indices for the CES [35,36].Here, AP is represented by two indices, SO 2 discharge per unit area (SODI), and dust discharge per unit area (DDI).WP is expressed by waste water discharge per unit area (WWDI), COD (Chemical Oxygen Demand) discharge per unit area (CODI) and AN discharge per unit area (ANDI).

Data Acquisition
In this paper, BRC was selected as the spatial scale with the years 2000, 2005, and 2010 chosen for investigation (reflecting the end of every National Five-year Plan of China).Remote sensing data from Landsat TM and ETM+ were downloaded from the data sharing platform Geospatial Data Cloud [37].Multi-spectral color images with a spatial resolution of 15 m were generated though band fusion.These images were interpreted and classified by the working group from the Institute of Remote Sensing and Digital Earth (RADI), Chinese Academy of Sciences (CAS) into several land-use classes, including urban build-up, road, and sea reclamation, to quantify their areas and get the coastline length according to visual interpretation under the referencing of Google Earth.Based on the interpretation results, LCA, LUR, TI, and SRI were calculated, respectively, where the sea reclamation area of 2000 is interpreted on the basis of the coastline of 1992.Economic and population data were gathered from the statistical yearbook of local governments and pollutant emission data of each city were acquired from China's ministry of environmental protection.These data were used to calculate each of the remaining indices.

Normalization Processing
Since not all the indices are in the same dimension, it can be difficult to compare and generalize them to produce and compare the comprehensive CDI and CES indices [38,39].Normalization of the indices was, therefore, conducted in order to keep all the indices for the three time nodes in the same dimensionless normalized state.The calculation formula was presented as follow ([12,28,40]): where Y i is the standardized value of index X i ; the value range is (0, 1); X i is the original data.X max is the maximum value of the same index of the three time nodes in the whole study area (13 cities).

Comprehensive Indices Computation by Principal Component Analysis
The indices selected above may have significant commonalities or are related to each other when they are used to describe CDI and CES.Therefore, simply stacking the indices randomly would make the information they carry overlapping and complex, obfuscating the results.A reasonable weight for each index should, therefore, be assigned prior the calculation of the CDI and CES.Principal component analysis (PCA) using SPSS 19.0 software (IBM, New York, NY, America) was used to determine the relative weight of each index.PCA is a statistical procedure that uses an orthogonal transformation to convert a set of observations of possibly-correlated variables into a set of values of linearly uncorrelated variables called principal components equal to or less than the number of original variables [41,42].
In order to keep the indices consistent and comparable in time series, the same weight group of CDI or CES would be acquired for each time node by taking all samples in the three time nodes (altogether 39 samples) as the original input array of each sub-indices.

Dynamic Degree
The dynamic degree method was introduced to calculate the yearly growth rate of CDI and CES.The formula is presented as follow: where DD is the comprehensive dynamic degree.DC i is the value of CDI or CES in the initial time node, while DC e is value at the end time node; n is the length of time series.The unit is the integral year.

Correlation Analysis
Correlation analysis with SPSS19.0 software was used to study the relationship between CDI and CES for all cities in each time node.The Pearson coefficient of product-moment correlation was set as the correlation coefficient.Results were tested using the t-statistic to inspect if there were significant relationships between the two variables and the relationship degree was divided into four classes when significant relationships existed (p value was less than 0.05) (Table 2).

Coupling Degree Model
Through correlation analysis, the relationship between the CDI and CES indices and its evolution can be determined.However, even having known the relationship between the CDI and CES, the interaction between development and environment was extremely complex.Once development was over-fasted, it would throw a huge pressure on environment.While the increased speed of CDI was proper, it would promote the environment in a virtuous cycle.There are all kinds of contradictions and interactions between the two subsystems; thus, the system which was composed of CDI and CES can be defined as a coupling system [43].Coupling degree is a physical quantity that describes the degree of influence of two objects on each other [35,44].In this paper, the coupling degree was used to quantitatively describe the interactional between the CDI and CES subsystems and to study the development characteristics and tendency of urban system at different coupling levels.Building on previous research [11,28,39,40,43,[45][46][47][48], the CDI-CES coupling model was constructed using the following equations: T " pCDI `CESq {2 where CD is the coupling degree of CDI and CES, and would be used to show the coupling relationship of the two subsystems in each time node for Bohai Rim coastal cities, C is the coordination degree of CDI and CES, and T is the comprehensive evaluation index of CDI and CES, indicating the overall benefit or level of development intensity and ecological environment.CD reaches 1 when CDI and CES both reach their highest value of 1.The model synthesizes the coordination status and the levels of CDI and CES, which is designed for quantitative evaluation and comparison of the coupling status of the two indices for different cities at the same time node or for a single city at different time nodes.

Urban Development Type and Classification System
The coupling degree not only shows the coupling level of CDI and CES, but also indicates the development type of each city.High CD values indicate that with rising CDI, the stress to the ecological environment rises, too.This means that the level of overall urban development is increasingly dependent on the consumption of the ecological environment, which may occur when both the CDI and CES are very high.In contrast, low CD values indicate the lower dependence, which usually appears at the beginning of urban development or at the mutual promotion stage of urban development intensity and ecological environment protection [39].To explore the priority order of urban development intensity and ecological environment protection in the promotion of the urban development level in terms of each city, values of CDI, CES, and CD are all reclassified and jointed to produce different urban development types.CDI and CES were classified into five grades by equal interval, from very low to very high for each index, as shown in Table 3. CD is reclassified into four classes, referring to the previous literatures [39,40,43,44] (Table 4).Based on Tables 3 and 4 different classes for the three indices were associated by permutations and combinations and generated various urban development types ("CD-CDI-CES").For example, "M-M-L" represents a city development style with a middle level coupling degree between CDI and CES, and CDI is at a medium level while CES is at a low level.Comparatively speaking, there is still space for the enhancing of urban development intensity without exerting too much stress on the ecological environment and trying to keep the coupling degree at a lower value.

Temporal and Spatial Distribution of CDI and CES
First, sub-indices of both CDI and CES were calculated.Then, all the results were turned into a normalized state to make further comparisons and computations.All of the sub-indices were drawn in Figure 2 to make a comparison among all the cities for the three time nodes.
"M-M-L" represents a city development style with a middle level coupling degree between CDI and CES, and CDI is at a medium level while CES is at a low level.Comparatively speaking, there is still space for the enhancing of urban development intensity without exerting too much stress on the ecological environment and trying to keep the coupling degree at a lower value.

Temporal and Spatial Distribution of CDI and CES
First, sub-indices of both CDI and CES were calculated.Then, all the results were turned into a normalized state to make further comparisons and computations.All of the sub-indices were drawn in Figure 2 to make a comparison among all the cities for the three time nodes.Principal component analysis shows that the first three principal components (PC) cumulate more than 85% of the variance in the sub-index correlation matrix when the initial Eigen value is larger than 1.The component matrix for the sub-index and the total variance is explained for every PC in SPSS, and the results were used to obtain the final weight group for CDI and CES, respectively.The comprehensive index CDI and CES were calculated by a weighted average for each city at every time node (year 2000, 2005, and 2010), and results were shown in Tables 5 and 6 and statistical graphs were made (Figures 3 and 4).In the time series, the CDI of the whole study area increased from low in 2000 to medium in 2005 and 2010, and the CDI of all of the cities increased monotonically.Among those cities, CDI of BH rose from medium in 2000 to high in 2005, and finally to very high in 2010 at rate of greater than all other cities for the same period.The development intensity at Tangshan was low in 2000, but then quickly rose to medium during the latter two time nodes, ranking second after Binhai in 2010.Huludao had the lowest CDI throughout the time series with the lowest dynamic degree.The CDI for Dalian, Panjin, and Dongying was medium for all the three time nodes.The remainder of the cities had low CDI values at the beginning of the time series, reaching medium at 2010.For the other cities, CES increased through the study period.

Correlation Analysis of CDI and CES in BRC
For each time node, correlation analysis was carried out between the CDI and CES values of all the cities in the study area, and the correlation coefficients (r) were listed in Table 7. Form Table 7, it can be seen that CDI and CES presented positive correlations for the entire time series.CES showed a rising tendency with the enhancement of CDI.In 2000, the correlation coefficient between CDI and CES was 0.389, with no obvious correlation, while it rose to 0.645 in 2005, when the two indices were significantly correlated (p = 0.05).In 2010, the CDI and CES were significantly correlated (p = 0.01) with r reaching 0.841.According to the correlation analysis, we can see that the correlation between CDI and CES in the whole study area became stronger with time, which has sufficiently demonstrated that there is a certain relationship between CDI and CES, though we cannot make a conclusion only by correlation analysis that the decrease of ecological environment quality is definitely caused by the enhancing of urban development intensity.Therefore, using the index of coupling degree, the interaction characteristic and coupling relationship of CDI and CES were studied next in this paper to reveal the urban development feature and type.This provides the reference advice for the urban development planning of the Bohai Rim region.

Coupling Degree of CDI and CES in BRC
The CDI and CES coupling degree was calculated for every city at all the time nodes, and results were shown in Figure 5.

Correlation Analysis of CDI and CES in BRC
For each time node, correlation analysis was carried out between the CDI and CES values of all the cities in the study area, and the correlation coefficients (r) were listed in Table 7. Form Table 7, it can be seen that CDI and CES presented positive correlations for the entire time series.CES showed a rising tendency with the enhancement of CDI.In 2000, the correlation coefficient between CDI and CES was 0.389, with no obvious correlation, while it rose to 0.645 in 2005, when the two indices were significantly correlated (p = 0.05).In 2010, the CDI and CES were significantly correlated (p = 0.01) with r reaching 0.841.According to the correlation analysis, we can see that the correlation between CDI and CES in the whole study area became stronger with time, which has sufficiently demonstrated that there is a certain relationship between CDI and CES, though we cannot make a conclusion only by correlation analysis that the decrease of ecological environment quality is definitely caused by the enhancing of urban development intensity.Therefore, using the index of coupling degree, the interaction characteristic and coupling relationship of CDI and CES were studied next in this paper to reveal the urban development feature and type.This provides the reference advice for the urban development planning of the Bohai Rim region.

Coupling Degree of CDI and CES in BRC
The CDI and CES coupling degree was calculated for every city at all the time nodes, and results were shown in Figure 5. From Figure 4, it can be seen that the coupling degree presents an increasing trend in the time series in the whole study area.CD values of Binhai and Tangshan rank first and second, From Figure 4, it can be seen that the coupling degree presents an increasing trend in the time series in the whole study area.CD values of Binhai and Tangshan rank first and second, respectively, in all the cities and rose fast during the time series.Tangshan had been at the medium level all the time, while Binhai reached a high level of coupling at 2010 (0.87), much higher than the other cities.Among the rest of the cities, CD values rose with time, but Dalian, Jinzhou, Panjin, and Dongying changed little, basically keeping stable, and Huludao and Dongying had been in the state of low level coupling, just reaching the lower limit value of medium level (0.5) in 2010.

Urban Development Type in BRC
According to Table 5, Table 6, and the coupling in each time node, we constructed the urban development type for all of the cities in the study area (Table 8).
From Table 8 we can see that for the whole study area, urban development type has a tendency of changing from M-L-L to M-M-L in the time series, with CDI rising to the medium level and CES remaining in the low level.
For the cities of Binhai, Tangshan, and Qinhuangdao, each experienced great changes in urban development types.Binhai's urban development type changed from M-M-L in 2000 to M-H-M in 2005 and rose to H-VH-M in 2010, where the coupling degree rose from the medium to the high level with CDI rising from medium to the very high and CES from low to medium.Tangshan changed from M-L-L in 2000 to M-M-M in 2005 and M-H-M in 2010, which shows that although the coupling degree was unchanged in the time series, CDI and CES both rose with time.With the quick growth of CDI, ecological environment in Binhai and Tangshan have suffered heavy stresses, which means that urban development was dependent on the destruction of the ecological environment to a certain degree.People in both cities must pay great attention to the protection of the ecological environment when CDI is rising.For Qinhuangdao, although it also experienced great changes in urban development types (from L-L-VL to M-L-VL and M-M-L), its CDI remained at low and medium levels and, at the same time, CES stayed at very low and low levels, which shows that urban development activities emphasized environmental protection.As such, there is still space for Qinhuangdao in the enhancing CDI.Dalian is the only city with CES dropping from low to very low and, at the same time, with CDI rising from low to medium, although the coupling degree remained at medium.The urban development types of Panjin, Huludao, and Dongying were unchanged during the study period (M-M-VL, L-L-VL, and L-M-VL, respectively).

Discussion
The results above show that during the last decade, all of the cities in the study area had experienced rapid growth of urban development intensity, and most of them had exerted stresses on the ecological environment.For the whole study area, CDI and CES, and their coupling degree, were all in a monotonically rising trend.And the relationships between CDI and CES for each time node were rising with time.It can be concluded that human activities of urban development were partly dependent on the consumption of the ecological environment, and the degree of dependence increased with year.This resulted in the low coordination degree between urban development and ecological environment, especially at the end of the study period.This phenomenon is more serious with a high coupling degree between CDI and CES.Binhai is a typical example of this finding.Since the beginning of this century, the development of Binhai has become a national strategy with a national economic policy supporting a rapid increase of the CDI that exceed its growth at other regional locations.The increase of CDI was bound to produce large amounts of pollutant emission, and if environmental governance measures were insufficient, CES would increase, leading to a high coupling degree between CDI and CES.Dalian is the only city in the study area that showed urban development, but with limited ecological impact.This indicates that Dalian was able to enhance urban development without seriously impacting their environment.On consideration of environment carrying capacity, Dalian strengthens the red line of ecological protection and environmental risk, rationally allocating economic and social environment resources, and gives full consideration of the collaborative relationship of development and environmental protection.Thus, a sound city environment management system has been established to guarantee the environmental protection while carrying on urban development.That is to say, Dalian had setup a typical model for sustainable development in the whole Bohai Rim coastal area.In order to enhance the whole urban development level, people should promote CDI, as well as decrease CES as much as possible, keeping the coupling degree at a lower level.
Low coupling degree does not definitely show a high urban development level.All three factors, including CDI, CES, and their coupling degree, should be taken into consideration in the study of urban development types.To ensure healthy urban development, the coupling degree should be kept down by maintaining higher CDI and lower CES.A low coupling degree often occurs at the initial stage of urban development or a high stage of urban development when urban development intensity and ecological environment protection keep in a mutual promotion stage [39].At the intermediate or middle stage of urban development, the coupling degree was usually higher.A high level of urban development should be based on not only the high development intensity, but also the high quality of the ecological environment.
Regarding the urban development type, most of the cities were of the M-L-L type (TS, CZ, DL, YK, YT, and WF) and L-L-VL type (QHD, JZ, HLD, and BZ) in 2000, which means most of the cities had a medium or low coupling degree with low CDI and low or very low CES at the beginning of the time series.With low urban development intensity, the urban ecological environment was less affected by human activities, keeping CDI and CES in a low or medium coupling degree.With the demand of economic development, the CDI increased in 2005 for most of those cities, and peaked in 2010 when most of the cities were concentrated in the M-M-L type (QHD, CZ, YK, YT, and WF) and the M-M-VL type (DL, PJ, and BZ).By 2010 most cities had discernable impacts on ecological environments.It reached the antagonism stage, when the environment stresses increased with the enhancement of urban development intensity.With further urban development, CES and CDI would run in a stage of mutually touching and adjusting.If the environment deterioration were left with no effective governance, it would develop towards a disorder state.On the other hand, if strong measures were taken to protect the ecological environment, CES would decrease and, in turn, the coupling degree between CES and CDI would decrease, too.Thus, urban development and environment protection would be running in a direction of high coordination.Therefore, in future urban development activities, more attention should be paid to the environment protection to keep CES at a lower level.The government should focus more on improving environmental quality than maximizing economic development [49].In order to promote urban development, urban development intensity could be enhanced on the premise of keeping eco-environment stresses at the low or very low level.

Conclusions
Taking Bohai Rim coastal area as the study area, this paper first made a spatiotemporal contrastive analysis on the correlation between CDI and CES.Then, by building a model of coupling degree, the authors comparatively analyzed the coupling degree of CDI and CES to show the urban development characteristic of each city and its development level, and further studied the orientation and trend of urban development in the study area.
For the whole study area, the CDI and CES increased monotonically during the study period.Within cities, the Binhai area of Tianjin and Tangshan ranked first and second with higher CDI and CES at the same time.The relationships between CDI and CES for each time node rose with time.It can be concluded that human activities of urban development were partly dependent on the consumption of the ecological environment, and the degree of dependence increased with year.All three factors, including CDI, CES, and their coupling degree, should be taken into consideration in the study of urban development levels.High levels of urban development should be based not only on the high development intensity, but also the high quality of the ecological environment.In order to promote sustainable urban development, urban development intensity could be enhanced on the premise of keeping eco-environment stresses at the low or very low level.Dalian is a good example of this.
Since many factors influence CES, they cannot be listed one by one.Based on limited impact factors, the stress level is relatively low in this paper, but according to the aim of in this paper, it can meet the requirements to make comparisons and reflect the change in the study area and in the time sequence based on the same method.

Figure 1 .
Figure 1.Location of the Bohai coastal area.

Figure 1 .
Figure 1.Location of the Bohai coastal area.

Figure 2 .
Figure 2. Statistical graph of sub-indices in the study area.Figure 2. Statistical graph of sub-indices in the study area.

Figure 2 .
Figure 2. Statistical graph of sub-indices in the study area.Figure 2. Statistical graph of sub-indices in the study area.

Figure 3 .
Figure 3. Statistical graph of CDI in the study area.

Figure 4 .
Figure 4. Statistical graph of CES in the study area.Over the study period, CES values for the whole study area increased monotonically.CES values for Binhai and Tangshan both increased from low in 2000 to medium in 2005 and 2010, and in the first two time nodes (2000 and 2005), Tangshan had the highest CES value while, in 2010, Binhai's CES value became the highest, reaching the lower limit of the high level (0.6).In Dalian and Dongying, CES rose in 2005 and dropped in 2010, indicating the stress of ecological environment was reduced in 2010.Jinzhou, Panjin, and Weifang remained steady for the study period with low dynamic degrees.For the other cities, CES increased through the study period.

Figure 3 .
Figure 3. Statistical graph of CDI in the study area.

Figure 3 .
Figure 3. Statistical graph of CDI in the study area.

Figure 4 .
Figure 4. Statistical graph of CES in the study area.Over the study period, CES values for the whole study area increased monotonically.CES values for Binhai and Tangshan both increased from low in 2000 to medium in 2005 and 2010, and in the first two time nodes (2000 and 2005), Tangshan had the highest CES value while, in 2010, Binhai's CES value became the highest, reaching the lower limit of the high level (0.6).In Dalian and Dongying, CES rose in 2005 and dropped in 2010, indicating the stress of ecological environment was reduced in 2010.Jinzhou, Panjin, and Weifang remained steady for the study period with low dynamic degrees.For the other cities, CES increased through the study period.

Figure 4 .
Figure 4. Statistical graph of CES in the study area.

Figure 5 .
Figure 5. Statistical graph of coupling degrees of CDI and CES.

Figure 5 .
Figure 5. Statistical graph of coupling degrees of CDI and CES.

Table 1 .
Assessment indices system of CDI and CES in the study area.
CDI LCD 100 Ai is classification value of No. i land-use degree; Ci is area percent of No. i land-use degree [31].EAI EAI = GDP/A Ais the total area; GDP is the total GDP of study area.UI PUR PUR = Pu/P Pu is the urban population; P is the total population of study area.EUR EUR = (GDP2 + GDP3)/GDP GDP2 and GDP3 represent secondary industry and tertiary industry; GDP is the total GDP of study area.LUR LUR = Au/A Au is the area of urban built-up; A is the total area.TI TI = At/A At is the total road area in study area; A is the total area.SRI SRI = Ar/L Ar is the total sea reclamation area of study area in every time span; L is the length of coastline in 1992.

Table 1 .
Assessment indices system of CDI and CES in the study area.
i is the ith study area.

Table 2 .
Reclassification of related coefficients.

Table 3 .
Grading table for CDI and CES.

Table 4 .
Coupling degree classification of CDI and CES.

Table 4 .
Coupling degree classification of CDI and CES.

Table 5 .
Statistical table of CDI and its dynamic degree.

Table 6 .
Statistical table of CES and its dynamic degree.

Table 7 .
Correlation statistical table of CDI and CES at each time node.

Table 7 .
Correlation statistical table of CDI and CES at each time node.

Table 8 .
Urban development type for each city in the study area.