Simulating the Impact of Economic and Environmental Strategies on Future Urban Growth Scenarios in Ningbo , China

Coastal cities in China are challenged by multiple growth paths and strategies related to demands in the housing market, economic growth and eco-system protection. This paper examines the effects of conflicting strategies between economic growth and environmental protection on future urban scenarios in Ningbo, China, through logistic-regression-based cellular automata (termed LogCA) modeling. The LogCA model is calibrated based on the observed urban patterns in 1990 and 2015, and applied to simulate four future scenarios in 2040, including (a) the Norm-scenario, a baseline scenario that maintains the 1990–2015 growth rate; (b) the GDP-scenario, a GDP-oriented growth scenario emphasizing the development in city centers and along economic corridors; (c) the Slow-scenario, a slow-growth scenario considering the potential downward trend of the housing market in China; and (d) the Eco-scenario, a slow-growth scenario emphasizing natural conservation and ecosystem protections. The CA parameters of the Normand Slow-scenarios are the same as the calibrated parameters, while the parameters of proximities to economic corridors and natural scenery sites were increased by a factor of 3 for the GDPand Eco-scenarios, respectively. The Normand GDP-scenarios predicted 1950 km2 of new growth for the next 25 years, the Slow-scenario predicted 650 km2, and the Eco-scenario predicted less growth than the Slow-scenario. The locations where the newly built-up area will emerge are significantly different under the four scenarios and the Slowand Eco-scenarios are preferable to achieve long-term sustainability. The scenarios are not only helpful for exploring sustainable urban development options in China, but also serve as a reference for adjusting the urban planning and land policies.


Introduction
The coastal cities are the most densely populated urban settlements and are also growing very rapidly.Particularly in China, rapid growth and globalization of the domestic economy have dramatically accelerated urban growth [1,2].In 2014, there were 88 municipalities and prefecture-level cities each with more than 5 million people, nearly 40% of which were located in coastal areas [3].Many cities are increasing uncontrollably in size in terms of both their non-agricultural population and the size of the built-up areas [4].This is largely driven by the growing demand for higher living standards, higher incomes and better-living conditions [5].The increasing effects of human activities on the physical environment of coastal cities have resulted in many unplanned or underplanned land use and land cover changes [6], undermining long-term sustainability.For instance, rapid urbanization and industrialization have led to serious environmental problems such as air and water pollutions [7][8][9][10].Inevitably, there is growing public awareness of and demand for environmental protection.Different strategic directions for coastal city development have been proposed under the somehow conflicting demands for economic growth and environmental protection.Consequently, there is an urgent need to predict and assess the potential impact of different development strategies and the urban growth scenarios that these strategies would produce.
Among the above models, the logistic CA model (termed as LogCA hereafter) was the earliest approach proposed by Wu [13] to capture and calibrate the CA transition rules to simulate the rural-to-urban land transition in Guangzhou of China.This type of CA model has been widely adopted to simulate dynamic urban growth due to its effectiveness in addressing the urban dynamics [16,18].Liu and Feng [17] extended the LogCA model by applying continuous cell states to represent the progressive rural-to-urban land transition, with a case study of Gold Coast, Australia.Reference [28] incorporated a Markov chain into the LogCA model to build a hybrid CA model for modeling dynamic land change in the Poyang Lake region of China.These studies show that the structure of the LogCA model allows for flexibility in incorporating changes into development regulations; it also allows for detailed spatial interaction settings to model urban growth drivers and processes.While geographical information system (GIS) based CA models have been applied in a range of urban contexts to simulate urban growth, most of the existing research focuses on model calibration and testing.It is more challenging to predict future growth, particularly in the rapidly changing contemporary Chinese economy and its urban landscape [36].It is necessary for policymakers to use science-and evidence-based scenario modeling to guide future urban development and be able to adequately assess the impact of different strategies, whether the focus is on economic growth or environmental protection [37].Although CA-based urban and land use modeling in Chinese cities have been extensively documented [19,22,28,38,39], few are concerned with predicting possible future scenarios, especially at the current turning point in time concerning urbanization, population growth, economic prosperity, and environmental protection [40].Ningbo is an important harbor city situated in the east coastal area of Zhejiang province, China.During the past 25 years, from 1990 to 2015, the city has witnessed rapid urban growth that has been accelerated by the prosperous housing market.Urbanization encroachment on green space, agricultural land, wetlands and forests has led to the degradation of environments and ecosystems [41].This stimulated growing public awareness of and demand for environmental protection.Different strategic directions for coastal development have been proposed under somehow conflicting demands for economic growth and environmental protection.As a result, a crucial question has been aroused as to what the future scenarios would be like in Ningbo.
This paper calibrates the LogCA model to simulate future urban scenarios of Ningbo, a rapidly evolving coastal city in eastern China.The LogCA model was built using land use patterns classified from remotely sensed imagery during 1990-2015.Considering the impacts of economic development and environmental protection, four scenarios of urban growth to the year 2040 in Ningbo were constructed: (a) the Norm-scenario, a baseline scenario that maintains the urban growth rate from 1990 to 2015 and rate and bases on the normal transitions in economy; (b) the GDP-scenario, a GDP growth oriented scenario considering the economic corridors; (c) the Slow-scenario, the growth rate is much lower than in the past considering the potential downward trend of the housing market in China; and (d) the Eco-scenario, an ecosystem protection centered scenario considering the needs of environmental protection.By modeling the outcomes of different urban growth strategies, the impact of these strategies on Ningbo's urban growth will be assessed and compared.The scenarios are not only helpful for exploring future urban development options in Ningbo, but also useful as a reference for decision makers to adjust relevant planning and land management policies in Chinese coastal areas.

Study Area
Ningbo is a sub-provincial and port city located in Zhejiang Province, about 220 km south of Shanghai in eastern China.It has a total area of 9671 km 2 .The area under study is 8653 km 2 which includes inland waters and coastal wetlands but exclude sea water.Ningbo comprises the Ningbo City, three satellite cities (i.e., Cixi, Yuyao, and Fenghua), and two counties (Ninghai and Xiangshan) (Figure 1).With a total population of 7.6 million, there were 3.5 million urban inhabitants in Ningbo in 2014 [42].A beautiful port city with natural landscapes and many cultural attractions, Ningbo is also an important economic hub of Zhejiang Province and a major exporter of food, textiles, electrical products, and industrial tools [43].In 2014, the total gross domestic product (GDP) of Ningbo was 760.2 billion RMB Yuan (approximately $117 billion USD), 80 percent of which was from the private sector located in Ningbo City and its three satellite cities [42].
environmental protection.By modeling the outcomes of different urban growth strategies, the impact of these strategies on Ningbo's urban growth will be assessed and compared.The scenarios are not only helpful for exploring future urban development options in Ningbo, but also useful as a reference for decision makers to adjust relevant planning and land management policies in Chinese coastal areas.

Study Area
Ningbo is a sub-provincial and port city located in Zhejiang Province, about 220 kilometers south of Shanghai in eastern China.It has a total area of 9671 km 2 .The area under study is 8653 km 2 which includes inland waters and coastal wetlands but exclude sea water.Ningbo comprises the Ningbo City, three satellite cities (i.e., Cixi, Yuyao, and Fenghua), and two counties (Ninghai and Xiangshan) (Figure 1).With a total population of 7.6 million, there were 3.5 million urban inhabitants in Ningbo in 2014 [42].A beautiful port city with natural landscapes and many cultural attractions, Ningbo is also an important economic hub of Zhejiang Province and a major exporter of food, textiles, electrical products, and industrial tools [43].In 2014, the total gross domestic product (GDP) of Ningbo was 760.2 billion RMB Yuan (approximately $117 billion USD), 80 percent of which was from the private sector located in Ningbo City and its three satellite cities [42].Like many other Chinese cities, Ningbo has experienced rapid urban growth over the past 25 years (1990-2015).This growth has resulted in dramatic land use changes with a tenfold expansion of its urban extent, from 213 km 2 in 1990 to 2163 km 2 in 2015.This remarkable change has been accelerated since 1998 when property commercialization commenced.Rapid urban expansion and the boom in the housing market have threatened the natural landscape and cultural attractions in Ningbo [44].Subsequently, there are conflicting demands for land resources by various community groups which could result in different development options and affect its future urban growth scenarios.Planners and government authorities must explore and address the question of how the city should grow in the future under different strategies.Like many other Chinese cities, Ningbo has experienced rapid urban growth over the past 25 years (1990-2015).This growth has resulted in dramatic land use changes with a tenfold expansion of its urban extent, from 213 km 2 in 1990 to 2163 km 2 in 2015.This remarkable change has been accelerated since 1998 when property commercialization commenced.Rapid urban expansion and the boom in the housing market have threatened the natural landscape and cultural attractions in Ningbo [44].Subsequently, there are conflicting demands for land resources by various community groups which could result in different development options and affect its future urban growth scenarios.Planners and government authorities must explore and address the question of how the city should grow in the future under different strategies.

Data
Three types of data were collected: vector data (maps), remotely sensed images and socio-economic datasets (Table 1).The vector data include the administrative boundary map of the city, the transportation map in 2000 and 2015 Ningbo Municipal Tourism Administration [45].Although many local roads have been constructed in Ningbo after the year 2000, its major road network was constructed prior to the year 2000.Consequently, the 2000 transportation map was used to identify the main road network and quantify the spatial proximity to main roads as one of the key factors that impact urban growth.The tourism map was used to identify natural scenery sites for modeling the ecosystem protection oriented scenario.Two Landsat images, one in 1990 and one in 2015, were collected and classified using the supervised Mahalanobis distance classifier in ENVI version 5.5 to identify the past and current land use patterns.Using the region of interest (ROI) identified for training the classifier, the overall accuracies were 96.3% and 95.7% for 1990 and 2015, respectively.Digital elevation model (DEM) was assembled from the Shuttle Radar Topography Mission (SRTM) to extract land slope and measure land development suitability [46].The socio-economic data collected from Ningbo Statistical Yearbook [42] were used to estimate the total quantity of developable land and identify the economic corridors for modeling the GDP growth oriented scenario.Three urban centers (including Ningbo City, Cixi, Yuyao), three town centers (including Zhenhai, Guangcheng and Zhouxiang), and the Port of Ningbo were identified as key areas that have contributed to GDP growth (Figure 2, left box).The five-star natural scenery sites whose environment and ecosystem status have been well maintained in the past two decades were identified and used to model the ecosystem protection oriented scenario (Figure 2, right box).All data were resampled to 60 m spatial resolution except for the SRTM DEM to construct the LogCA model.

Variables Used in the LogCA Model
Ten driving factors, including accessibility and proximity-based variables, neighborhood, slope, global and local constraints, and a stochastic perturbation factor (Table 2), were identified to simulate urban growth in Ningbo from 1990 to 2015.These factors were selected considering the geographical and socio-economic conditions of Ningbo as well as by referencing factors affecting the urban growth dynamics of similar coastal cities in China, such as Shanghai [19,22,25,26].Of the six distance-based variables, Duc, Dtc, Dmr, and Dw have been commonly used in CA models [13,19,47], while Dec and Dns were selected to reflect the economic condition and tourism attraction features in Ningbo.Specifically, Dec was used to model the GDP growth oriented scenario while Dns was used for simulating the ecosystem protection oriented scenario.A stochastic factor was included to represent the factors that might not be clear and the departure from a well-defined process of urban [13,19].Water bodies and wetlands were considered as local constraints while the total quantity available for development was served as global constraints.Among such factors, the spatial variables and neighborhood factor represent the agglomerative effects of existing built-up areas and supports

Stochastic perturbation due to unexpected errors and policy change
Returned randomly using R = 1 + (−ln r) a , where r is a random number in the range from 0 to 1 and α = 4 Constraint Constraints for urban development Computed in ENVI from remote sensing imagery Of the six distance-based variables, D uc , D tc , D mr , and D w have been commonly used in CA models [13,19,47], while D ec and D ns were selected to reflect the economic condition and tourism attraction features in Ningbo.Specifically, D ec was used to model the GDP growth oriented scenario while D ns was used for simulating the ecosystem protection oriented scenario.A stochastic factor was included to represent the factors that might not be clear and the departure from a well-defined process of urban [13,19].Water bodies and wetlands were considered as local constraints while the total quantity available for development was served as global constraints.Among such factors, the spatial variables and neighborhood factor represent the agglomerative effects of existing built-up areas and supports of infrastructure in contrast to the stochastic and other constraint factors [48][49][50][51].All spatial and slope variables were normalized using the following equation [26,40]: where D norm is the normalized value of a variable in the range of [0, 1], D orig is the original value, D max is the maximum value, and D min is the minimum value.To prepare the spatial variables as input layers, they were spatially visualized and mapped in ArcGIS (Figure 3). of infrastructure in contrast to the stochastic and other constraint factors [48][49][50][51].All spatial and slope variables were normalized using the following equation [26,40]: where is the normalized value of a variable in the range of [0, 1], is the original value, is the maximum value, and is the minimum value.To prepare the spatial variables as input layers, they were spatially visualized and mapped in ArcGIS (Figure 3).

CA Model Calibrated by Logistic Regression
A LogCA model based on logistic regression was built for modeling the spatiotemporal process of urban growth in Ningbo.The parameters of transition rules of the LogCA model were calculated

CA Model Calibrated by Logistic Regression
A LogCA model based on logistic regression was built for modeling the spatiotemporal process of urban growth in Ningbo.The parameters of transition rules of the LogCA model were calculated using a logistic regression method [13,17,18,28,52].In a tessellated space, CA determines the state of a non-urban cell i at time t + 1 as the integrated effects of itself and its neighboring cells at time t according to land transition rules.One-directional rural-to-urban conversion across space and over time was modeled in this study.The probability of the cell being converted into an urban land use state is defined as [22]: where Pg i,t is the overall probability of possible conversion of a cell at location i at time t; this is determined by the accumulated effects of its spatial distances, the neighborhood effects, local and global constraints, and the stochastic factors.Pl i,t denotes the local probability of land conversion for cell i at time t as a result of a set of local spatial factors, which can be given by: where x j (j = 1, 2, . . ., k) are driving factors affecting urban growth, a 0 is a constant or intercept, and a j (j = 1, 2, . . ., k) are CA parameters indicating the extent to which each spatial factor affects urban growth.Spatial factors commonly include the proximity of a cell to the urban center, district centers, main roads and other urban facilities.These parameters were derived using logistic regression given that this modeling approach has been proven as reliable in capturing CA parameters [17,28,53].Pn i,t in Equation ( 2) represents the conversion probability being affected by the state of other cells within its neighborhood.A square neighborhood within m × m cells is commonly used, and the probability of the conversion is given as [54]: where ∑ m×m j (S i,t=Urban ) indicates the sum of urbanized cells within the neighborhood, excluding the central cell i.
The conversion probability of a non-urban cell is also affected by other constraints, for example, large-scale water bodies and steep areas.These areas would be less likely to be developed due to the high development costs.Primary farmlands are another kinds of areas being constrained from development due to land use regulations.All constraint factors in Equation ( 2) are represented as Con and can be written as: where Con is either 0 or 1, with 0 meaning the cell i is unavailable, while 1 is available to be developed.R in Equation ( 2) is a stochastic factor indicating unforeseen policy changes and other unknown errors [19,25,55].The factor is defined as: where r is a randomly generated real number between 0 and 1, and α is a control parameter between 0 and 10 to adjust the impact of R on land conversion probability.The overall conversion probability Pg i,t is compared with a pre-defined threshold value P thd to determine whether the cell i at time t can be developed [13,56].If the overall conversion probability is larger than P thd , the non-urbanized cell will be developed into an urban use in the next time step; otherwise, the cell i will not change its state.The state of cell i at time t + 1, denoted as S i,t+1 , can be written as: Urban, Pg i,t ≥ P thd NonUrban, Pg i,t < P thd (7) 3. Results

Transition Rules Defined in LogCA
Using the stratified sampling method, 8664 sample points were selected from each of the spatial distance data layers to construct the CA transition rules.The CA parameters were retrieved at a p = 0.01 significance level using R-Gui (Table 3).Six parameters were involved in Equation (3) to compute the conversion probability of each cell in the simulation process, therefore to model land use change in Ningbo.
A negative a j (j = 0, 1, . . ., 6) results in a larger Pg i,t , i.e., a higher potential for cell i to transit from a non-urban into an urban state.In contrast, a positive a i leads to in a lower Pg i,t and hence, a lower potential for the cell to develop further in a subsequent time step.While D uc , D tc , D mr , and D ec are all negative, the absolute value of D uc is the smallest, indicating that the impact of Ningbo's city center is weaker than the other three spatial variables.This is because Ningbo is not a mono-center city like other coastal cities such as Shanghai [26].Rather, the satellite cities of Cixi and Yuyao have experienced remarkable urban expansion from 1990 to 2015.Among the spatial variables, the distances to the economic corridors (D ec ), as well as to main roads (D mr ), have the largest impact on development.Conversely, D w and D ns are positive, meaning that a cell closer to the water body or natural scenery has a lower probability of being developed into an urban state.D w has the smallest absolute value amongst all spatial variables, indicating that this factor is least influential on urban growth.When predicting the future urban expansion scenarios to the year 2040, the values of D ec and D ns were increased by a factor of 3 to represent the GDP and Eco-scenarios, respectively.The parameter is increased by a factor of 3 as the calibrated value to reflect different scenarios.

Model Calibration
For setting the threshold, we selected 16 candidate values between 0.50 and 0.80 with an interval of 0.02.The LogCA model was implemented based on each candidate threshold with 25 iterations, with each iteration representing one year from 1990 to 2015.The LogCA model therefore generated 16 replicates for Ningbo in 2015.Among all replicates, a threshold of 0.62 produced the result with the highest overall accuracy.We therefore simulated the land use pattern in 2015 with a threshold of 0.62 and 25 iterations, which was compared with the observed urban pattern in 2015.The observed urban growth from 1990 to 2015 was the overlaying map of the two classified patterns in 1990 and 2015 (Figure 4a), and the simulated urban growth was the overlaying map of the 1990 classified pattern and the 2015 simulated pattern (Figure 4b).
To assess simulation accuracy, we applied the measures of quantity and allocation agreement and disagreement developed by Pontius, et al. [57].The results show that the model has achieved an overall agreement of 92.4% with 7.6% allocation disagreement (Figure 5).No quantity disagreement was observed, as the total land available for development was used to restrict the quantity of urban growth in the model.Of the 92.4% accuracy achieved by the LogCA model, 71.1% is due to the correct rejection of non-urban cells while 21.3% is due to the correct simulation of growth.Figure 5 also shows that both the commission and omission errors for urban and non-urban areas are 3.8%, resulting in a total of 7.6% errors.Such simulation accuracies are considered high and acceptable; hence, the model can be used for projecting future urban growth scenarios for the next 25 years under different development strategies.

Simulated Scenarios of Ningbo for 2040
Four different development strategies for Ningbo to the year 2040 were simulated based on the calibrated CA model from 1990 to 2015.These include (a) the Norm-scenario, a baseline scenario that maintains the urban growth rate at the 1990-2015 period and assumes a continuous economic transition as from the past; (b) the GDP-scenario, a GDP growth oriented scenario considering the accelerated growth along its economic centers and corridors; (c) the Slow-scenario, where urban growth will slow down and the growth rate is much lower than in the past, considering the potential downward trend of the housing market in China; and (d) the Eco-scenario, an ecosystem and environmental protection oriented growth scenario.
The Norm-and GDP-scenarios assume that urban growth in the next 25 years (2015-2040) will continue as in the past 25 years.According to the two classified Landsat images, urban land use has

Simulated Scenarios of Ningbo for 2040
Four different development strategies for Ningbo to the year 2040 were simulated based on the calibrated CA model from 1990 to 2015.These include (a) the Norm-scenario, a baseline scenario that maintains the urban growth rate at the 1990-2015 period and assumes a continuous economic transition as from the past; (b) the GDP-scenario, a GDP growth oriented scenario considering the accelerated growth along its economic centers and corridors; (c) the Slow-scenario, where urban growth will slow down and the growth rate is much lower than in the past, considering the potential downward trend of the housing market in China; and (d) the Eco-scenario, an ecosystem and environmental protection oriented growth scenario.
The Norm-and GDP-scenarios assume that urban growth in the next 25 years (2015-2040) will continue as in the past 25 years.According to the two classified Landsat images, urban land use has

Simulated Scenarios of Ningbo for 2040
Four different development strategies for Ningbo to the year 2040 were simulated based on the calibrated CA model from 1990 to 2015.These include (a) the Norm-scenario, a baseline scenario that maintains the urban growth rate at the 1990-2015 period and assumes a continuous economic transition as from the past; (b) the GDP-scenario, a GDP growth oriented scenario considering the accelerated growth along its economic centers and corridors; (c) the Slow-scenario, where urban growth will slow down and the growth rate is much lower than in the past, considering the potential downward trend of the housing market in China; and (d) the Eco-scenario, an ecosystem and environmental protection oriented growth scenario.
The Norm-and GDP-scenarios assume that urban growth in the next 25 years (2015-2040) will continue as in the past 25 years.According to the two classified Landsat images, urban land use has increased by 1950 km 2 from 1990 to 2015, hence it is assumed that another 1950 km 2 of urban area will be developed in the next 25 years under the Norm-and GDP-scenarios, resulting in a total urban area of 4113 km 2 in 2040.For the Norm-scenario, the same set of model parameters as defined during the model calibration process was used to run the model until the desired quantity of urban growth is achieved.However, for the GDP-scenario, the value of D ec was increased by a factor of three to reflect the assumption that strong growth will occur along the economic corridors to achieve higher GDP growth (see Table 3 under the GDP-scenario column).
The Slow-and Eco-scenarios were defined under the assumption that urban growth in the region will slow down in the next 25 years, given the booming building construction in China and some clear signals indicating that the rate of urban growth will slow down [4].Considering the increasing tight regulation on the property market as well as the strong call for ecosystem protection, we assumed that the maximum allowable area of urban growth for the next 25 years will be limited to one-third of that during the 1990-2015 interval, i.e., about 650 km 2 .In the Eco-scenario, urban expansion will be tightly constrained at locations near important natural scenery sites.This scenario is reflected by increasing the D ns values by a factor of three, as compared to the calibrated value from past development (see Table 3, under Eco-scenario).The urban growth scenarios under the four development strategies are illustrated in Figure 6.
increased by 1950 km 2 from 1990 to 2015, hence it is assumed that another 1950 km 2 of urban area will be developed in the next 25 years under the Norm-and GDP-scenarios, resulting in a total urban area of 4113 km 2 in 2040.For the Norm-scenario, the same set of model parameters as defined during the model calibration process was used to run the model until the desired quantity of urban growth is achieved.However, for the GDP-scenario, the value of Dec was increased by a factor of three to reflect the assumption that strong growth will occur along the economic corridors to achieve higher GDP growth (see Table 3 under the GDP-scenario column).
The Slow-and Eco-scenarios were defined under the assumption that urban growth in the region will slow down in the next 25 years, given the booming building construction in China and some clear signals indicating that the rate of urban growth will slow down [4].Considering the increasing tight regulation on the property market as well as the strong call for ecosystem protection, we assumed that the maximum allowable area of urban growth for the next 25 years will be limited to one-third of that during the 1990-2015 interval, i.e., about 650 km 2 .In the Eco-scenario, urban expansion will be tightly constrained at locations near important natural scenery sites.This scenario is reflected by increasing the Dns values by a factor of three, as compared to the calibrated value from past development (see Table 3, under Eco-scenario).The urban growth scenarios under the four development strategies are illustrated in Figure 6.

Norm-Scenario: A Baseline Scenario with a Growth Rate Observed during 1990-2015
The Norm-scenario assumes steady urban growth and a normal transition in the economy, resulting in the same scale of new development for the coming 25 years.Although such an assumption may not be supported by China's current economic situation and housing market, it provides an insight to examine how the city would be if it continues at the same rate of urban growth.Under the Norm-scenario, Ningbo will gain an urban area of 1950 km 2 , resulting in a total urban area of 4113 km 2 in 2040.Most urban development will occur near the existing Ningbo urban center and its three satellite cities: Cixi, Yuyao and Fenghua (Figure 6a).The Fenghua city center will be merged into the Ningbo city as a fringe suburb, while both Cixi and Yuyao will grow to be connected to the Ningbo city center, except for the junctions of the three regions, where the natural scenery is characterized by a steep topographic slope.

GDP-Scenario with Strong Growth along Its Economic Corridors
Urban centers typically contribute a much larger portion to the overall economy in most cities in China [58].In addition, the private sector in Ningbo has been a significant player in its economic growth and is becoming increasingly important [59].The economic corridors (city center, town centers in Zhenhai, Cixi and Yuyao) also play significant roles to its overall economic growth.Therefore, the city center, town centers and economic corridors would have a higher probability of being developed to boost higher GDP growth.The CA model was run by applying the distance to the economic corridors parameter (D ec ) three times larger than the Norm-scenario to generate the same quantity of urban land in 2040, hence the model reallocates new urban growth to unconstrained areas.As in the Norm-scenario, Figure 6b shows that most urban development occurs near the existing Ningbo city center and its three satellite cities of Cixi, Yuyao and Fenghua.However, the GDP-scenario generated more urban area in the economic corridors that also contain some of the natural sceneries at the junction of Ningbo city center, Cixi and Yuyao.Gains are smaller at Fenghua as compared to the Norm-scenario.This suggests that GDP-oriented urban growth will occur at the expense of the natural environment.

Slow-Scenario Accommodating Lower Housing and Property Demand
China's housing market has reached a turning point with falling housing prices and demand as well as an oversupply of commercial properties [60].Considering the potential slowdown of the booming housing market, property development companies have lowered expectations of new residential housing demand [61].Consequently, the Slow-scenario reduces the allowable land for urban growth to 650 km 2 , resulting in a total urban area of 2813 km 2 in 2040.Under this scenario, urban growth will occur mainly around the existing urban area of Ningbo City and Fenghua, and a small portion of growth will occur in Cixi (Figure 6c).We anticipate that this scenario is most likely to occur, given the tightening policy on investment properties, as well as the increasing concern about environmental issues and demand for ecosystem protections.

Eco-Scenario: A Slow-Growth Scenario that Protects the Environment and Its Ecosystems
A total of 24 natural sceneries located in Ningbo were identified along the coast and in the mountain areas with relatively high land slope, for example, the junction of Ningbo city center, Cixi and Yuyao.These areas are of high ecological value and should be protected from urban development under the Eco-scenario.This scenario was realized by increasing the value of D ns (distance to the natural scenery sites) by a factor of three to relocate new urban growth to other unconstrained areas.As a result, this scenario only generated 522 km 2 of new urban area in 2040 compared to the base year in 2015-about 128 km 2 less than in the Slow-scenario. Figure 6d shows a similar pattern to the Slow-scenario, where urban growth mainly occurs around the existing urban area of Ningbo city center and Fenghua, but gains less in Cixi.

Discussion
Significant increases in the urban area for Ningbo by the year 2040 were predicted under each of the four urban growth scenarios.Even with the conservative strategies which led to the Eco-and Slow-scenarios, the urban areas in Ningbo will increase by 520 to 650 km 2 .The increase in urban land use results from increases in population in the Ningbo city core, rapid development in the private sector in Cixi, Yuyao and Fenghua, and development of the logistics service industry at the port of Ningbo.
On the other hand, the four simulations produced few changes in Ninghai and Xiangshan, where the natural ecological environment remains in good condition and the tourism industry is operating at a high level.To visualize the differences amongst the four scenarios, two comparison maps were produced based on the differences between the Norm-vs.GDP-scenarios and between the Slow-vs.Eco-scenarios (Figure 7).Percentage of urban land change is presented in Figure 8 for the four scenarios using the map comparison method proposed by Pontius and colleagues [57,62,63] and is used to assess and compare the results of CA models [24,26].Figure 7 shows that, compared to the GDP-scenario, the Norm-scenario simulated more urban growth in Fenghua and to the north of Cixi, while there was less growth in Yuyao and at the junction of Ningbo city center, Cixi and Yuyao, where most of the private industry is located.Compared with the Eco-scenario, the Slow-scenario generated more urban growth to the northwest of Ningbo city, and the administrative centers of Cixi and Yuyao where the areas were considered as protected areas under the Eco-scenario.

Discussion
Significant increases in the urban area for Ningbo by the year 2040 were predicted under each of the four urban growth scenarios.Even with the conservative strategies which led to the Eco-and Slow-scenarios, the urban areas in Ningbo will increase by 520 to 650 km 2 .The increase in urban land use results from increases in population in the Ningbo city core, rapid development in the private sector in Cixi, Yuyao and Fenghua, and development of the logistics service industry at the port of Ningbo.
On the other hand, the four simulations produced few changes in Ninghai and Xiangshan, where the natural ecological environment remains in good condition and the tourism industry is operating at a high level.To visualize the differences amongst the four scenarios, two comparison maps were produced based on the differences between the Norm-vs.GDP-scenarios and between the Slow-vs.Eco-scenarios (Figure 7).Percentage of urban land change is presented in Figure 8 for the four scenarios using the map comparison method proposed by Pontius and colleagues [57,62,63] and is used to assess and compare the results of CA models [24,26].Figure 7 shows that, compared to the GDP-scenario, the Norm-scenario simulated more urban growth in Fenghua and to the north of Cixi, while there was less growth in Yuyao and at the junction of Ningbo city center, Cixi and Yuyao, where most of the private industry is located.Compared with the Eco-scenario, the Slow-scenario generated more urban growth to the northwest of Ningbo city, and the administrative centers of Cixi and Yuyao where the areas were considered as protected areas under the Eco-scenario.Both the Norm-and GDP-scenarios predicted that nearly half (i.e., 47.5%, or 4113 km 2 out of its total land areas of 8653 km 2 ) will become urban areas by the year 2040.This includes 25.0% of existing urban areas as in 2015 (i.e., persistence) and 22.5% of new urban gains (i.e., 1950 km 2 ).Under these two scenarios, Ningbo will become a highly urbanized area except for Ninghai and Xiangshan counties.However, this will be at the cost of deteriorating its beautiful natural scenery, the environment and ecosystem.On the other hand, urban areas will increase by 7.5% and 6.0%, or 650 km 2 and 522 km 2 , by the year 2040 under the Slow-and Eco-scenarios, respectively.Considering the Both the Norm-and GDP-scenarios predicted that nearly half (i.e., 47.5%, or 4113 km 2 out of its total land areas of 8653 km 2 ) will become urban areas by the year 2040.This includes 25.0% of existing urban areas as in 2015 (i.e., persistence) and 22.5% of new urban gains (i.e., 1950 km 2 ).Under these two scenarios, Ningbo will become a highly urbanized area except for Ninghai and Xiangshan counties.However, this will be at the cost of deteriorating its beautiful natural scenery, the environment and ecosystem.On the other hand, urban areas will increase by 7.5% and 6.0%, or 650 km 2 and 522 km 2 , by the year 2040 under the Slow-and Eco-scenarios, respectively.Considering the size of non-urbanized land being converted into urban use and the long-term sustainability of the region, the Slow-and Eco-scenarios are preferable given that these scenarios preserve more land from urban use and put more emphasis on environmental protection and the harmony between humans and nature.The locations where new urban growth will occur differ in the four scenarios, in response to different emphases such as the housing market, GDP growth and protection of the environment and the ecosystem.The scenarios generated by the model are useful for decision makers to assess different urban development options before any planning decision is made.size of non-urbanized land being converted into urban use and the long-term sustainability of the region, the Slow-and Eco-scenarios are preferable given that these scenarios preserve more land from urban use and put more emphasis on environmental protection and the harmony between humans and nature.The locations where new urban growth will occur differ in the four scenarios, in response to different emphases such as the housing market, GDP growth and protection of the environment and the ecosystem.The scenarios generated by the model are useful for decision makers to assess different urban development options before any planning decision is made.

Conclusions
Coastal city development is being increasingly challenged in an era of rapid urbanization and intensive human activity.It is important to understand the physical processes in urban systems and the human-environment interactions to plan for long-term sustainable development of cities.China's population increase and urban development in the coastal cities have been phenomenal and are continuously growing very rapidly, however, there are signals showing that a turning point is emerging in the housing market in China.The well-known air pollution and worsening environments in China-mainly caused by the large-scale urbanization and industrialization-have drawn public attention, resulting in a growing awareness of the protection of the environment and ecosystem.To simulate the future scenarios, CA modeling serves as a diagnosis for the urban development in China as well as assesses the impacts of economic development and environmental protection strategies on urban growth.
This paper applies a LogCA model to predict the urban scenarios of a coastal city, Ningbo, in southeast China.The LogCA model was validated based on the urban patterns of 1990 and 2015 and used to predict the scenarios in 2040.Four scenarios were considered including a normal scenario with an urban growth rate as observed 1990-2015, a GDP growth oriented scenario, a scenario with a lower expectation of urban growth rate, and an ecosystem protection oriented scenario.The Normand GDP-scenarios predicted 1950 km 2 of urban growth for the coming 25 years, the Slow-scenario predicted 650 km 2 of urban growth, and the Eco-scenario predicted less than the Slow-scenario.The locations where the new built-up cells occur are quite different for the four scenarios in response to different demands such as the housing market, the growth of the GDP and protection of the ecosystem.The Slow-and Eco-scenarios are preferable in terms of sustainability and urban planning policy could pay more attention to environmental protection and restoration and put more emphasis on the harmony between man and nature.These scenarios are helpful for decision makers who are exploring various sustainable urban development options before making any planning adjustment.

Conclusions
Coastal city development is being increasingly challenged in an era of rapid urbanization and intensive human activity.It is important to understand the physical processes in urban systems and the human-environment interactions to plan for long-term sustainable development of cities.China's population increase and urban development in the coastal cities have been phenomenal and are continuously growing very rapidly, however, there are signals showing that a turning point is emerging in the housing market in China.The well-known air pollution and worsening environments in China-mainly caused by the large-scale urbanization and industrialization-have drawn public attention, resulting in a growing awareness of the protection of the environment and ecosystem.To simulate the future scenarios, CA modeling serves as a diagnosis for the urban development in China as well as assesses the impacts of economic development and environmental protection strategies on urban growth.
This paper applies a LogCA model to predict the urban scenarios of a coastal city, Ningbo, in southeast China.The LogCA model was validated based on the urban patterns of 1990 and 2015 and used to predict the scenarios in 2040.Four scenarios were considered including a normal scenario with an urban growth rate as observed 1990-2015, a GDP growth oriented scenario, a scenario with a lower expectation of urban growth rate, and an ecosystem protection oriented scenario.The Normand GDP-scenarios predicted 1950 km 2 of urban growth for the coming 25 years, the Slow-scenario predicted 650 km 2 of urban growth, and the Eco-scenario predicted less than the Slow-scenario.The locations where the new built-up cells occur are quite different for the four scenarios in response to different demands such as the housing market, the growth of the GDP and protection of the ecosystem.The Slow-and Eco-scenarios are preferable in terms of sustainability and urban planning policy could pay more attention to environmental protection and restoration and put more emphasis on the harmony between man and nature.These scenarios are helpful for decision makers who are exploring various sustainable urban development options before making any planning adjustment.

Figure 1 .
Figure 1.The study area: Ningbo city of Zhejiang Province, China.(a) Map of China and (b) Map of Ningbo City.

Figure 1 .
Figure 1.The study area: Ningbo city of Zhejiang Province, China.(a) Map of China and (b) Map of Ningbo City.

Figure 2 .
Figure 2. Economic corridors and natural scenery sites in Ningbo.Left: administrative centers and economic corridor; right: sites of natural scenery.

Table 2 ., 1 R
Dependent and independent variables extracted to retrieve CA transition rules.variable y is assigned to 1 if a land cell is converted from non-urban to urban; otherwise, y is assigned to 0. Duc Distance to urban centers Measured in ArcGIS from administrative map Dtc Distance to town centers Dmr Distance to main roads Measured in ArcGIS from road networks map Dw Distance to water body Measured in ArcGIS and ENVI from remote sensing imagery Dec Distance to economic corridors Measured in ArcGIS from the map of economic corridors Dns Distance to natural scenery sites Measured in ArcGIS from the map of natural scenery sites Slope Developing suitability based on DEM Computed based on SRTM 90 m DEM V4 Stochastic perturbation due to unexpected errors and policy change Returned randomly using = 1 + (−ln ) , where r is a random number in the range from 0 to 1 and α = 4 Constraint Constraints for urban development Computed in ENVI from remote sensing imagery

Figure 2 .
Figure 2. Economic corridors and natural scenery sites in Ningbo.Left: administrative centers and economic corridor; right: sites of natural scenery.
model can be used for projecting future urban growth scenarios for the next 25 years under different development strategies.

Figure 4 .
Figure 4. Comparison of urban growth observed from remote sensing and simulated using LogCA, from 1990 to 2015.(a) The observed urban growth derived from the two classified Landsat images; and (b) The simulated urban growth using the LogCA model.

Figure 5 .
Figure 5. Accuracies and errors in the LogCA model.

Figure 4 .
Figure 4. Comparison of urban growth observed from remote sensing and simulated using LogCA, from 1990 to 2015.(a) The observed urban growth derived from the two classified Landsat images; and (b) The simulated urban growth using the LogCA model.

Figure 4 .
Figure 4. Comparison of urban growth observed from remote sensing and simulated using LogCA, from 1990 to 2015.(a) The observed urban growth derived from the two classified Landsat images; and (b) The simulated urban growth using the LogCA model.

Figure 5 .
Figure 5. Accuracies and errors in the LogCA model.

Figure 5 .
Figure 5. Accuracies and errors in the LogCA model.

Figure 8 .
Figure 8. Percentage of land use change under four urban growth strategies from 2015 to 2040.

Figure 8 .
Figure 8. Percentage of land use change under four urban growth strategies from 2015 to 2040.

Table 1 .
Spatial and socio-economic datasets used to produce land use patterns and variables.

Table 2 .
Dependent and independent variables extracted to retrieve CA transition rules.
y Conversion label Dependent variable y is assigned to 1 if a land cell is converted from non-urban to urban; otherwise, y is assigned to 0. w Distance to water body Measured in ArcGIS and ENVI from remote sensing imagery D ec Distance to economic corridors Measured in ArcGIS from the map of economic corridors D ns Distance to natural scenery sites Measured in ArcGIS from the map of natural scenery sites

Table 3 .
Parameters generated using the logistic regression method for Ningbo.