Long-Term Regional Environmental Risk Assessment and Future Scenario Projection at Ningbo , China Coupling the Impact of Sea Level Rise

Regional environmental risk (RER) denotes potential threats to the natural environment, human health and socioeconomic development caused by specific risks. It is valuable to assess long-term RER in coastal areas with the increasing effects of global change. We proposed a new approach to assess coastal RER considering spatial factors using principal component analysis (PCA) and used a future land use simulation (FLUS) model to project future RER scenarios considering the impact of sea level rise (SLR). In our study, the RER status was classified in five levels as highest, high, medium, low and lowest. We evaluated the 30 m × 30 m gridded spatial pattern of the long-term RER at Ningbo of China by assessing its 1975–2015 history and projecting this to 2020–2050. Our results show that RER at Ningbo has increased substantially over the past 40 years and will slowly increase over the next 35 years. Ningbo’s city center and district centers are exposed to medium-to-highest RER, while the suburban areas are exposed to lowest-to-medium lower RER. Storm surges will lead to strong RER increases along the Ningbo coast, with the low-lying northern coast being more affected than the mountainous southern coast. RER at Ningbo is affected principally by the combined effects of increased human activity, rapid population growth, rapid industrialization, and unprecedented urbanization. This study provides early warnings to support practical regulation for disaster mitigation and environmental protection.


Introduction
Regional environmental risk (RER) denotes potential threats to the natural environment, human health and socioeconomic development caused by specific risks [1,2].RER can be expressed by uncertainty in the consequence and degree of environmental damage [2].For example, industrial and domestic wastes create potential environmental risks that may threaten the sustainability of urban systems [3,4].Land-use change, accompanied by accelerating industrialization and urbanization, has led to the alteration of the natural environment and ecosystem degradation [5,6].Intensified human impacts have led to accumulated greenhouse effects and global change, ultimately resulting in global warming and sea level rise (SLR) [7].SLR increases the intensity and frequency of storm surges, intensifying RER in coastal areas that are significantly affected by inundation [8].Coastal  be exposed to increased RER because the ecosystems and environments are inherently vulnerable.It is therefore essential to consider the impact of SLR in long-term RER assessment and future scenario projection in coastal areas.RER assessment is to synthetically evaluate the risk sources, risk receptors and risk retardation from a macro perspective [9,10], offering possible countermeasures that can reduce risks.Efforts have been made to the spatial and quantitative assessment of RER as it relates to specific combinations of risks.Aiming to develop risk-mitigation measures, modelers have assessed RER in inland and coastal zones regarding hazardous substances [11] and pollutants [9,12].From a spatial perspective, a few publications have evaluated RER caused by multiple risk sources [13], metal in soil [14], vulnerable environmental factors [15] and toxic compounds [16].These studies have mapped RER patterns to identify high-risk areas and support the formulation of sustainable development plans.In RER studies, a Pressure-State-Response (PSR) framework has commonly been applied to estimate the risk level [10,17], with the support of the analytical hierarchy process (AHP) and geographical information systems (GIS).The entropy method [18], fuzzy method [19] and expert evaluation method (EEM) [20] have also been applied to the assessment of ecological security related to RER.While AHP and EEM have been applied to evaluate RER patterns with the absence of prior knowledge of risk levels [17,21], they strongly depend on subjective model parameterization that can influence the credibility of the results.In coastal areas, the impacts of sea level rise on environments have commonly been studied using three basic approaches, i.e., behavior-oriented models [22], process-based models [23] and hybrid models [24].McLeod et al. [25] summarized these coastal impact models and noted a few of them including inundation models (e.g., GIS), SLAMM, BTELSS, DIVA, and SimCLIM.These models can predict a wide range of coastal environments such as morphological evolution and estuarine sedimentation.
Despite considerable progress in the study of RER, improved assessment methods must consider other physical and socioeconomic driving forces, particularly because coastal environments are now substantially affected by global change and SLR.As a result, critical issues related to RER include: (1) objective parameterization for assessing spatial patterns of RER, (2) proper projection of future RER scenarios, and (3) comprehensive assessment of future RER patterns considering the influences of SLR and storm surges.In this study, storm surge refers to coastal flooding commonly associated with typhoons, heavy rainfall and extreme high tides [26], which may increase social vulnerability and environmental risks.To date, answers to these questions have not been adequately addressed, and further studies are needed for model development and integrated environmental management.
The objective of this study is to adequately address the above issues by developing new methods for assessing RER patterns and projecting their future scenarios.Our new assessment method is based on the principal component analysis (PCA).It is a dimensionality reduction method whose outputs are directly derived from the independent factors [27], leading to the objective definition of the effects of these factors.Land coverage indices and distance-based factors were applied as input variables for the modeling.These variables are representative of the land coverage status, human impacts, and socioeconomic aspect.Our new method also incorporates GIS to facilitate the spatial visualization [16,17,[28][29][30] and the quantification [31][32][33] of RER status, providing a clearer understanding.We used a bottom-up cellular automata (CA) model to project future scenarios by processing the RER maps into categorical figures.CA models have been widely applied to simulate and predict spatiotemporal, dynamic phenomena [34,35].The models can automatically update their transition rules during the simulation procedure, and have substantial advantages of categorical modeling in tessellated space but weak capabilities to model continuous value-based spatial phenomena [36].We applied the future land use simulation (FLUS) model, a state-of-art CA tool for modeling complex land use change and spatiotemporal phenomena [35,37], to predict the RER future scenarios.

Procedure
For modeling RER (Figure 1), we first used remote sensing imagery and vector maps to derive the factors impacting RER at Ningbo.Based on these factors, we then applied PCA to evaluate the present RER pattern in 2015 and to assess the past RER patterns for 1975, 1990 and 2006.Consequently, the FLUS model of spatial phenomena [35,37] was used to project before-inundation RER scenarios from 2020 to 2050.We then applied a GIS approach (inundation model) to evaluate after-inundation RER scenarios from 2020 to 2050, taking into account terrain (from a digital elevation model: DEM), SLR and storm surge.Our long-term RER assessment should help formulate more appropriate plans for regional development and environmental protection, and could help establish defensive measures to rehabilitate coastal environments by taking SLR into account.Consequently, the FLUS model of spatial phenomena [35,37] was used to project before-inundation RER scenarios from 2020 to 2050.We then applied a GIS approach (inundation model) to evaluate after-inundation RER scenarios from 2020 to 2050, taking into account terrain (from a digital elevation model: DEM), SLR and storm surge.Our long-term RER assessment should help formulate more appropriate plans for regional development and environmental protection, and could help establish defensive measures to rehabilitate coastal environments by taking SLR into account.

Study Area and Input Data
Ningbo is a port city situated in the middle of China's east coast (Figure 2a).The city covers 9816 km 2 and consists of two county-level cities (Cixi and Yuyao), two counties (Ninghai and Xiangshan), and six districts (Jinzhou, Haishu, Jiangbei, Zhenhai, Beilun and Fenghua) (Figure 2b).Haishu and Jinzhou constitute the Ningbo City center.Ningbo is an economic center of the southern wing of the Yangtze River Delta [38], and is the eastern departure port for the 'Marine Silk Road' on the Yongjiang River estuary.

Study Area and Input Data
Ningbo is a port city situated in the middle of China's east coast (Figure 2a).The city covers 9816 km 2 and consists of two county-level cities (Cixi and Yuyao), two counties (Ninghai and Xiangshan), and six districts (Jinzhou, Haishu, Jiangbei, Zhenhai, Beilun and Fenghua) (Figure 2b).Haishu and Jinzhou constitute the Ningbo City center.Ningbo is an economic center of the southern wing of the Yangtze River Delta [38], and is the eastern departure port for the 'Marine Silk Road' on the Yongjiang River estuary.
As of 2015, Ningbo's population was ~7.8 million [39], and the annual GDP was ~801 billion Chinese Yuan (~121 billion USD) [39].A vast population and a vibrant economy have resulted in high demand for commercial and residential space, leading to unprecedented land-use changes that pose serious challenges to ecosystems and the natural environment.Due to the major changes in land-use and environmental pressure, Ningbo has attracted attention from environmentalists and local authorities [40,41] and is of great interest to environmentalists, geographers and scientists [42,43].
We used graphic GIS datasets and remote sensing imagery (Table 1) to derive spatial factors influencing RER at Ningbo.The Ningbo City center and district/county centers were defined using administrative maps, which differ between 1975-1990 and 2006-2015.Four Landsat images were used to compute the normalized difference vegetation index (NDVI) and the normalized difference building index (NDBI).The terrain at Ningbo was stable over the past four decades, thus we used only the 2010 DEM to guide the impact of terrain on future RER scenarios.
Ningbo is a port city situated in the middle of China's east coast (Figure 2a).The city covers 9816 km 2 and consists of two county-level cities (Cixi and Yuyao), two counties (Ninghai and Xiangshan), and six districts (Jinzhou, Haishu, Jiangbei, Zhenhai, Beilun and Fenghua) (Figure 2b).Haishu and Jinzhou constitute the Ningbo City center.Ningbo is an economic center of the southern wing of the Yangtze River Delta [38], and is the eastern departure port for the 'Marine Silk Road' on the Yongjiang River estuary.

Input Variables
Increasing human demand for natural and land resources has put great pressure on the natural environment [44] from facility construction, landscape alteration, urban development and land reclamation [45,46].This has led to global land-use change, reduction in vegetation coverage, and growth in built-up areas [47].We selected four typical spatial influencing factors to evaluate RER patterns at Ningbo: proximity to Ningbo city center (PNCC), proximity to district centers (PDC), NDVI and NDBI, as defined in Table 2. Using the ArcGIS Euclidean Distance tool, we calculated PNCC and PDC to represent the distance from each pixel and its nearest city/district centers; using the ArcGIS Map Algebra tool, we computed NDVI and NDBI to describe the land coverage status.

Variable
Meaning Purpose PNCC Proximity to Ningbo City center Produce the distance to the Ningbo City center and evaluate its impact on the environments.

PDC
Proximity to district centers Produce the distances to the district centers and evaluate their impact on the environments.

NDBI
Normalized difference built-up index Produce the built-up density and evaluate its impact on the environments.

NDVI Normalized difference vegetation index
Produce the vegetation coverage and evaluate its impact on the environments.
Ningbo City center and the district centers have the highest density of impervious surfaces, suggesting the major human alteration of the natural environments.The use of proximity factors is therefore effective in delineating RER caused by human disturbance.Generally, a higher PNCC (or PDC) suggests a lower building density and higher vegetation coverage, hence a weaker human impact and less-altered natural environments [20,21].We produced one PNCC map and two PDC maps using the administrative maps in 1990 and 2006 (Figure 3a-c).These two factors are temporally stationary when including them to assess the past RER scenarios as well as projecting the before-inundation scenarios.The topo-bathymetry was considered for the after-inundation scenarios.
Sustainability 2019, 11, x FOR PEER REVIEW 5 of 18 intensity of human-disturbance and probably population density [49], and a higher built-up density suggests a greater human alteration of the natural environments [50].Figure 3h-k shows the NDBI maps that measure the built-up density as an important indicator of human impacts, reflecting the status of human-made surfaces.The maps in Figure 3 were produced following Feng et al. [21].
Table 2. Selected spatial influencing factors for assessing regional environmental risk (RER) patterns in Ningbo.

PNCC Proximity to Ningbo City center
Produce the distance to the Ningbo City center and evaluate its impact on the environments.

PDC Proximity to district centers
Produce the distances to the district centers and evaluate their impact on the environments.

NDBI Normalized difference built-up index
Produce the built-up density and evaluate its impact on the environments.

NDVI Normalized difference vegetation index
Produce the vegetation coverage and evaluate its impact on the environments.Ningbo City center and the district centers have the highest density of impervious surfaces, suggesting the major human alteration of the natural environments.The use of proximity factors is therefore effective in delineating RER caused by human disturbance.Generally, a higher PNCC (or PDC) suggests a lower building density and higher vegetation coverage, hence a weaker human NDVI reflects vegetation coverage and is an important indicator of the integrity of ecological environments, with a higher NDVI indicating a less-altered natural environment [48].Figure 3d-g shows a decline in vegetation coverage at Ningbo from 1975 to 2015.The built-up density reflects the intensity of human-disturbance and probably population density [49], and a higher built-up density suggests a greater human alteration of the natural environments [50].Figure 3h-k shows the NDBI maps that measure the built-up density as an important indicator of human impacts, reflecting the status of human-made surfaces.The maps in Figure 3 were produced following Feng et al. [21].
Ningbo City center and the district centers have the highest density of impervious surfaces, suggesting the major human alteration of the natural environments.The use of proximity factors is therefore effective in delineating RER caused by human disturbance.Generally, a higher PNCC (or PDC) suggests a lower building density and higher vegetation coverage, hence a weaker human impact and less-altered natural environments [20,21].We produced one PNCC map and two PDC maps using the administrative maps in 1990 and 2006 (Figure 3a-c).These two factors are temporally stationary when including them to assess the past RER scenarios as well as projecting the before-inundation scenarios.The topo-bathymetry was considered for the after-inundation scenarios.
NDVI reflects vegetation coverage and is an important indicator of the integrity of ecological environments, with a higher NDVI indicating less-altered natural environment [48].Figure 3d-g shows a decline in vegetation coverage at Ningbo from 1975 to 2015.The built-up density reflects the intensity of human-disturbance and probably population density [49], and a higher built-up density suggests a greater human alteration of the natural environments [50].Figure 3h-k shows the NDBI maps that measure the built-up density as an important indicator of human impacts, reflecting the status of human-made surfaces.The maps in Figure 3 were produced following Feng et al. [21].

RER Assessment and Projection Methods
We summarized the RER assessment and projection methods in Table 3, then described them in detail in the following sub-sections.The RER score was defined as a real number in the range of [0, 1], where a smaller score denotes better environmental status.A risk score of 0 means that the natural environment is at zero risk with very high vegetation coverage and no pollution sources; a risk score of 1 indicates that the natural environment is at high risk and seriously damaged by high-density buildings.Artificial impervious surface area causes great environmental risks by altering the natural environment, reducing infiltration capacity, changing urban runoff, and aggravating urban water pollution.We categorized risk in five equally-spaced levels (Table 3) following earlier publications [20,21].

Principal Component Analysis (PCA)
PCA identifies uncorrelated principal components from candidate factors using an orthogonal transformation that captures internal structural relationships among the factors [51].The method processes the dataset without considering the dependent variables, and then estimates a weighting factor for each.In this study, the RER score (the dependent variable) is unknown.As such, PCA is an objective parameterization method [52] that is particularly suitable for retrieving the weight of each RER factor, avoiding adverse negative effects that may be caused by subjective parameterization in methods like expert evaluation.Each principal component reflects one aspect of the independent variables, and all selected principal components should ideally explain more than 80% of these variables, with no information repeated [52].
In this study, the first three principal components of selected factors yielded the total cumulative variances of 82.2%, 82.7%, 90.0% and 87.7% for 1975,1990,2006 and 2015, successively.The weight of each factor was then defined by the first three principal components and represented as: where F year represents the comprehensive RER score each year, and the items in brackets represent the integrity and goodness of the natural environments.A smaller parameter indicates a higher impact in affecting RER at Ningbo.The coefficients in Equation ( 1) show that NDBI is the dominant RER factor followed by PDC, PNCC and NDVI.

Projection Method of RER Scenarios
We applied a FLUS model to project future RER scenarios at Ningbo from 2020 to 2050 at 5-year intervals.FLUS is an explicit modeling tool for simulating land-use change dynamics and spatial phenomena [35,37].This model integrates interactions between top-to-bottom system dynamics (SD) and a bottom-to-top CA model to simulate spatial dynamics.It first uses SD to project the future dimensions of different types considering the effects of various driving factors, then allocates spatial change according to the quantity of each type using CA modeling.Here, FLUS was calibrated using candidate factors to simulate the 2015 RER and project future RER scenarios.We applied six accuracy and error indices (Table 3) to assess the 2015 simulation and validate the FLUS model [53].

SLR Inundation Analysis
The IPCC fourth assessment report (AR4) shows that the average rate of global SLR was about 1.8 (1.3 to 2.3 range) mm per year from 1961 to 2003, and about 3.1 (2.4 to 3.8 range) mm per year from 1993 to 2003.The IPCC fifth assessment report (AR5) suggests that the global mean SLR rate was about 3 mm per year from 1993 to 2012.Both AR4 and AR5 demonstrate that global warming caused by human impacts is the main driver of SLR.The State Ocean Administration of China (SOAC) reported that the average sea level in China rose by 0.11 m from 1980 to 2005, and the rate was 3.7 mm per year from 1993 to 2010 [54].This is faster than the global average SLR rate reported by IPCC.We applied the SOAC SLR estimate to assume that the maximum inundation of 10 m in Ningbo, thus evaluating the inundation-induced RER.
SLR exacerbates the frequency and intensity of storm surges [55,56], which are likely causes of damage and disaster in coastal areas.Storm surges accompanied by rising tides and heavy rainfalls lead to coastal flooding, with longer inundation hours implying deeper bathymetric inundation, thus greater environmental damage [26].We assume that inundation time increases by 2 h if the surge-induced water level rises by 1-m bathymetry; conversely, it takes 2 h for the water level to drop by 1-m bathymetry.For flooded areas, more inundation hours imply higher RER.We also assume that the maximum water-level rise in Ningbo is 10 m, which is related to the inundation of 20 h.We categorized inundation risk into five levels according to the water-level rise and inundation hours with each 4-h (2-m bathymetry) inundation related to one inundation level.These risk levels are lowest, low, medium, high and highest, in ascending order.For instance, an inundation-time shorter than 4 h is the lowest risk.Meanwhile, inundation produces pressure on the coastal environments and intensifies RER.We produced the inundation maps and assessed the inundation-induced RER considering the elevation, producing comprehensive RER maps that take into account the combined effects of human activities and SLR.

RER Spatial Patterns and Disparity
The four panels in Figure 4 show good environmental status in most areas in Ningbo, indicating only weak environmental risk in these areas; however, the RER deterioration is evident in Ningbo City center and the eight district centers.(Figure 4d).Among these sub-areas, Ningbo City center has experienced the highest RER increase, while Ninghai and Xiangshan experienced only marginal increases.In 1975 and 1990, most areas had excellent environmental conditions, and only a small fraction (mainly in the Ningbo City center) had a higher (~0.8)RER. Figure 4a,b shows that the areas (in the district centers) with "medium" RER slightly expanded from 1975 to 1990.By 2006, there was further environmental deterioration, and areas with "lowest" RER expanded in the Ningbo City center (Figure 4c).Areas experiencing environmental degradation have spread in the district centers, including Cixi, Yuyao, Fenghua, and Zhenhai, all of which show evident increases in RER.By 2015, the environmental conditions in the Ningbo City center and the district centers of Cixi, Yuyao, Zhenhai and Beilun had deteriorated substantially, with the negative extent significantly expanded (Figure 4d).Among these sub-areas, Ningbo City center has experienced the highest RER increase, while Ninghai and Xiangshan experienced only marginal increases.

Changes in RER Patterns
During all three periods, RER increased substantially (by up to 0.4) in the Ningbo City center and district centers where artificial building surfaces have been greatly developed (Figure 6).RER in other areas was unchanged, or in some cases was reduced by up to −0.During all three periods, RER increased substantially (by up to 0.4) in the Ningbo City center and district centers where artificial building surfaces have been greatly developed (Figure 6).RER in other areas was unchanged, or in some cases was reduced by up to −0.

Future RER Scenarios
To forecast future RER scenarios for Ningbo, we calibrated the FLUS model using the five-level RER pattern in 2015 as the starting map and utilized four driving factors PDC, NDVI, NDBI, and PNCC as independent variables.A simulated 2015 RER map was then generated to validate the FLUS model.Visual inspection (Figure 7) shows the similarity between the assessed RER pattern and the simulated RER pattern.

Future RER Scenarios
To forecast future RER scenarios for Ningbo, we calibrated the FLUS model using the five-level RER pattern in 2015 as the starting map and utilized four driving factors PDC, NDVI, NDBI, and PNCC as independent variables.A simulated 2015 RER map was then generated to validate the FLUS model.Visual inspection (Figure 7) shows the similarity between the assessed RER pattern and the simulated RER pattern.

Future RER Scenarios
To forecast future RER scenarios for Ningbo, we calibrated the FLUS model using the five-level RER pattern in 2015 as the starting map and utilized four driving factors PDC, NDVI, NDBI, and PNCC as independent variables.A simulated 2015 RER map was then generated to validate the FLUS model.Visual inspection (Figure 7) shows the similarity between the assessed RER pattern and the simulated RER pattern.Table 4 shows that overall accuracy is 77.9% and Kappa coefficient is 77.3%, indicating good simulation accuracy of the FLUS model.The producer's and user's accuracies exceed 80% for the "highest," "low" and "lowest" RER levels, and are about 70% for the "high" and "medium" levels.The commission errors are negatively related to the producer's accuracies, while the omission errors are negatively related to the user's accuracies, confirming the effectiveness of the simulation.Overall, the calibrated FLUS model is relatively accurate in simulating the 2015 RER pattern.We therefore applied this model to project the RER scenarios (Figure 8) from 2020 to 2050 at five-year intervals.While the overall RER scenario patterns do not show significant changes with time, there is a noticeable trend of increasing RER (Figure 8).The enlarged Beilun area demonstrates the expansion of the "highest" and the "lowest" areas.From 2020 to 2050, there is at least one large patch of "highest" RER in the central area of Ningbo City center, Cixi, Yuyao, Zhenhai and Beilun, whereas the "highest" risk patches in Fenghua, Ninghai and Xiangshan are relatively smaller.The "lowest" RERs occur in southern Yuyao, western Fenghua, and the Ninghai and Xiangshan suburbs.
Our statistics show that the "highest" and "lowest" risk areas will increase, while the other three categories will decrease in the coming 35 years (Table 5).According to our projections, from 2015 to 2050 the "highest" area will increase by 6.8% and the "lowest" area will increase by 9.5%.In contrast, the "high," "medium" and "low" risk categories will decrease in acreage by 4.3%, 4.5%, and 7.6%, respectively.This indicates a possible transformation from the central three risk levels to the most secure (lowest risk) level and the most insecure (highest risk) level over the coming 35 years.

Impact of Storm Surge on RER
We used ArcGIS 10.3 to produce the map (Figure 9a) showing the areas affected by the storm surge in the 10 m inundation case, with different color denoting different flooding hours.The map shows that areas closer to the shoreline experience more flooding hours.In a storm surge, most areas of Ningbo are not inundated and are therefore not influenced by the RER.The areas in red are most  We overlaid the 10 m inundation-caused RER map on the FLUS-projected RER map for each year from 2020 to 2050, producing the overall RER patterns in Figure 10 that consider the combined influence of human activities, land development, SLR, and storm surges.Visual inspection shows no significant differences in RER over time.RER is the highest in the low-lying areas along the coastline in the north and northeast, as well as in the Ningbo City center and district centers.In contrast, there are relatively low RER in Fenghua, Xiangshan and Ninghai.The enlarged area in the Beilun District reveals more detail of the expansion of the "highest" RER.The summary statistics presented in Table 8 confirm this expansion.For example, the "highest" risk areas expand by about 2% from 2020 to 2045   6 shows the area corresponding to each flooding time and each inundation depth interval.The red areas most likely to be affected by storm surge are about 503 km 2 .In the case of a 10 m inundation, the total area affected by the storm surges reaches about 2104 km 2 , covering 21% of the study area, and the areas under the inundation of 8-14 h are all greater than 200 km 2 .Using the correlation between inundation and RER, we mapped the RER pattern caused by a maximum inundation of 10 m (Figure 9b) and calculated the acreage of each level (Table 7).This confirms the expectation that the region closest to the coastline has a higher RER, displaying a pattern similar to the inundation map.There is an apparent decrease in the RER from the coast to inland areas, indicating a decaying influence of storm surges.The "highest" RER areas (Figure 9b in red) are distributed along the coasts of Yuyao, Cixi and Beilun in the north, but also partially along the coasts of Ninghai and Xiangshan in the south.The area with the "highest" RER accounts for about 6% of the study area (Figure 9b).Inundation-caused RER spreads across about 58% in Zhenhai, with the "medium" to "lowest" levels; whereas RER spreads across about 42% area in Xiangshan, but only 17% in Ninghai.We overlaid the 10 m inundation-caused RER map on the FLUS-projected RER map for each year from 2020 to 2050, producing the overall RER patterns in Figure 10 that consider the combined influence of human activities, land development, SLR, and storm surges.Visual inspection shows no significant differences in RER over time.RER is the highest in the low-lying areas along the coastline in the north and northeast, as well as in the Ningbo City center and district centers.In contrast, there are relatively low RER in Fenghua, Xiangshan and Ninghai.The enlarged area in the Beilun District reveals more detail of the expansion of the "highest" RER.The summary statistics presented in Table 8 confirm this expansion.For example, the "highest" risk areas expand by about 2% from 2020 to 2045 and then remain at 17.7% in the last five years, while the "high" risk areas are reduced by 2% from 2020 to 2040 then remain at 11% from 2040 to 2050.The "medium" RER will stay at ~21% of the study area, while "low" RER areas will decrease by 5%, and the "lowest" RER areas will increase by ~5%.These suggest that human activities, land development, SLR and storm surges will substantially increase the RER at Ningbo.Significant differences exist between the RER scenarios before and after considering the effect of storm surge.As compared with the before-inundation RER (see Table 5), after-inundation RER in the coast of Yuyao, Cixi, Zhenhai and Beilun increases from "lowest" to "highest" level (Table 8).In each year, the after-inundation areas with the "highest" and "high" risks increase by about 4% and 3.4% respectively, and the after-inundation areas with the "lowest" and "low" risks are reduced by about 4% and 5.2% respectively, as compared with the before-inundation RER.This indicates the strong influence of storm surges leading to increased RER, particularly in the low-lying northern coast.storm surge.As compared with the before-inundation RER (see Table 5), after-inundation RER in the coast of Yuyao, Cixi, Zhenhai and Beilun increases from "lowest" to "highest" level (Table 8).In each year, the after-inundation areas with the "highest" and "high" risks increase by about 4% and 3.4% respectively, and the after-inundation areas with the "lowest" and "low" risks are reduced by about 4% and 5.2% respectively, as compared with the before-inundation RER.This indicates the strong influence of storm surges leading to increased RER, particularly in the low-lying northern coast.Compared to earlier publications on RER, our study is featured by the fine-scale assessment and objective parameterization, as well as long-term projection using a CA model considering the influences of SLR and storm surges.The PCA method provided an objective assessment without knowing RER levels prior to modeling, and avoided subjective evaluation that may be caused by AHP and EEM [17,21].In addition, our 30 m × 30 m grid-based assessment provides more spatial details of the RER as compared with those conducted under larger grids and the administrative boundaries [21,57].This may be attributed to the inclusion of the same spatial scale factors such as the land coverage indices and the distance to administrative centers.
Future scenario projection is a crucial issue in RER research.We applied FLUS to project future scenarios based on the classified RER pattern and its driving factors.FLUS is effective in simulating land-use and urban growth [35,37], and our study shows that it can be useful in modeling other spatial phenomena such as RER change.It should be noted that, however, the calibration of the CA-based FLUS model may affect RER simulation and projection.The influencing factors of calibration include the FLUS controlling parameterization, the sampling method, the number of samples, and the neighborhood configuration.While it is unclear how much the calibration may influence the 2015 simulation and future projection, there is a relatively high agreement between the 2015 simulation and assessment.The model shows lower accuracies (~70%) for the "high" and "medium" levels but higher accuracies (>85%) for the "highest" and "lowest" levels.This suggests FLUS may more accurately identify areas with the highest RER, and correctly simulate most of the "lowest" risk areas that should not be included in other categories with higher risks.The calibrated FLUS model is therefore useful in projecting future RER scenarios, and can accurately identify the "highest" RER areas for disaster prevention and the "lowest" RER areas for environmental protection.The future "highest" RER areas also represent an early warning of environmental status change.

Major Findings and Additional Explanations
Our results demonstrate that RER at Ningbo has increased substantially over the past 40 years, and will continue to increase slowly in the next 35 years.Ningbo City center and associated district centers experience high RER while the suburbs have relatively lower RER (Figure 4).The historical 40-year RER changes can be attributed to human impacts that are intensified by population growth, rapid urbanization, extensive industrialization, and pollutant emissions [46,58].RER assessment and ecological security evaluation share similarities, as both measure environmental conditions from a spatial perspective.Our results confirmed early studies that indicated deterioration in the ecosystem and natural environmental conditions at Ningbo [21].The relatively high RER in Ningbo City center and district centers is mainly caused by human activities and rapid urbanization, accompanied by land surface alteration [40,59].However, not all the areas studied have experienced increasing RER during the past 40 years.For example, the RER declined in parts of Jiangbei during 1990-2006 (Figure 6b), ascribed to environmental protection and urban park construction.
In addition to technical issues such as the FLUS calibration, the predicted future RER changes will probably be attributable to the aforementioned human impacts, intensifying natural disasters, and global warming.The before-inundation RER scenarios have more negative consequences in the Ningbo City center and district centers, and the after-inundation RER scenarios will significantly increase at the coast, and the RER will be lower in the suburban areas over time.The change may be caused by suburb-to-downtown population migration [60], urban sprawl, and industrialization that produces serious environmental pressures [46], and global change.

Limitations and Implications
Our research focused on RER scenario projection by coupling the impact of sea level rise loosely.We produced future scenarios affected by the decadal-scale inundation, considering the topography, bathymetry, and rising water level.This suggests our assessment of after-inundation scenarios is limited in extreme flood events at a specific time point.The process-based FLUS model is not the same category as those for modeling the impacts of the sea level rise.The former model focuses on modeling land-use change and spatial phenomena in terrestrial environments; in contrast, the latter models (e.g., SLAMM, BTELSS, DIVA, and SimCLIM) centralize the dynamic processes in coastal environments and are sophisticated in addressing the sea level rise.To address the dynamic processes in coastal environments, we instead applied a static GIS approach to spatially visualize the inundation by coupling the impacts of sea level rise.Therefore, future RER assessment work should include more unrelated risk factors and eco-environmental vulnerability on wider regional scales by considering the evolving topo-bathymetry.
This study identifies specific areas with different levels of potential environmental risks by model assessment and projection, providing valuable information for assessing the past urbanization and the environmental protection policies as well as for developing regulations for the management of urban development and coastal environments.For example, Ningbo authorities and urban planners can modify present regulations to achieve a more environmentally friendly city based on our research results.Landscape configuration [61] should be considered intently in residential development to reduce urban stormwater runoff and urban heat island effects, thus mitigating the environmental risk sources.The risks from sea level rise and storm-induced flooding should be closely considered in urban planning, land development, environmental protection and coastal management, thus reducing potential environmental risks and the loss of life and property.

Conclusions
We developed a PCA-based method for the spatially-explicit assessment of past RER patterns in coastal areas considering the impacts of spatial factors such as land coverage and spatial proximity, and then calibrated the FLUS model to project future RER scenarios by classifying RER patterns into categorical maps.Our evaluation and projection were at a 30 m resolution that provided appropriate spatial detail for effective environmental management.The PCA-based objective evaluation shows that historic RER at Ningbo has substantially increased over the past 40 years, while the FLUS-based projection shows that the future RER will slowly increase over the next 35 years.Inundation-impacted future scenarios can inform us about high RER areas that will probably be most affected by SLR and storm surges.Our research contributes to an improved understanding of the long-term (past and future) changes in RER patterns at Ningbo.
Our contribution towards RER research is the development of the PCA-based assessment method that features objective parameterization as well as the use of the FLUS model for the spatial prediction of future scenarios.The PCA used here is readily applicable to the assessment of RER elsewhere, and should be useful to evaluate the ecological security status.The current RER pattern could help authorities form environmental restoration measures, while the future scenarios are early warnings of impending changes in environmental status, supporting the development of practical policies and regulations for disaster mitigation and environmental protection.

Figure 1 .
Figure 1.Workflow used in this study.

Figure 1 .
Figure 1.Workflow used in this study.

Figure 2 .
Figure 2. Location and administrative boundaries of Ningbo City: (a) false-color Landsat composite image for 2006; and (b) a 2010 map for terrain aspect reflected in 10 azimuthal categories.

Figure 2 .
Figure 2. Location and administrative boundaries of Ningbo City: (a) false-color Landsat composite image for 2006; and (b) a 2010 map for terrain aspect reflected in 10 azimuthal categories.

Figure 5
Figure 5 presents histograms of the RER scores in each of the four years with a bin width of 0.05.All four years show left-skewed and platykurtic distributions, where the highest frequency in 1975, 1990, 2006 and 2015 occurs between 0.05-0.10,0.10-0.15,0.3-0.35,and 0.3-0.35,respectively.During 1975-2015, the highest RER score increased from 0.28 to 0.91, and the mean RER score substantially increased from 0.08 to 0.40.These indicate a movement of the highest frequency towards a higher risk score, suggesting an increasing RER in Ningbo.The standard deviation has slightly increased from 0.03 to 0.13, implying an increasing spatial heterogeneity in RER over the past 40 years.

Figure 5
presents histograms of the RER scores in each of the four years with a bin width of 0.05.All four years show left-skewed and platykurtic distributions, where the highest frequency in 1975, 1990, 2006 and 2015 occurs between 0.05-0.10,0.10-0.15,0.3-0.35,and 0.3-0.35,respectively.During 1975-2015, the highest RER score increased from 0.28 to 0.91, and the mean RER score substantially increased from 0.08 to 0.40.These indicate a movement of the highest frequency towards a higher risk score, suggesting an increasing RER in Ningbo.The standard deviation has slightly increased from 0.03 to 0.13, implying an increasing spatial heterogeneity in RER over the past 40 years.All four years show left-skewed and platykurtic distributions, where the highest frequency in1975,  1990, 2006  and 2015 occurs between 0.05-0.10,0.10-0.15,0.3-0.35,and 0.3-0.35,respectively.During 1975-2015, the highest RER score increased from 0.28 to 0.91, and the mean RER score substantially increased from 0.08 to 0.40.These indicate a movement of the highest frequency towards a higher risk score, suggesting an increasing RER in Ningbo.The standard deviation has slightly increased from 0.03 to 0.13, implying an increasing spatial heterogeneity in RER over the past 40 years.

4 .
There was a distinct risk-level change from 1975 to 1990, with higher risk levels found in the Ningbo City center and most district centers.The risk increases were spatially clustered and extended from 1990 to 2015, showing a similar tendency in risk change during all three time intervals.By comparison, those areas with increased risk during 2006-2015 were significantly larger than those during 1975-1990 and 1990-2006.Sustainability 2019, 11, x FOR PEER REVIEW 9 of 18

4 .
There was a distinct risklevel change from 1975 to 1990, with higher risk levels found in the Ningbo City center and most district centers.The risk increases were spatially clustered and extended from 1990 to 2015, showing a similar tendency in risk change during all three time intervals.By comparison, those areas with increased risk during 2006-2015 were significantly larger than those during 1975-1990 and 1990-2006.

Figure 6 .
Figure 6.RER change at Ningbo: red denotes increased risk, yellow denotes persistent (steady) risk, green denotes decreased risk.

Figure 6 .
Figure 6.RER change at Ningbo: red denotes increased risk, yellow denotes persistent (steady) risk, green denotes decreased risk.

Figure 6 .
Figure 6.RER change at Ningbo: red denotes increased risk, yellow denotes persistent (steady) risk, green denotes decreased risk.

Figure 7 .
Figure 7. Assessment and simulation of 2015 five-level-based RER patterns.Figure 7. Assessment and simulation of 2015 five-level-based RER patterns.

Figure 7 .
Figure 7. Assessment and simulation of 2015 five-level-based RER patterns.Figure 7. Assessment and simulation of 2015 five-level-based RER patterns.

Figure 8 .
Figure 8. RER future projections at Ningbo from 2020 to 2050.Figure 8. RER future projections at Ningbo from 2020 to 2050.

Figure 8 .
Figure 8. RER future projections at Ningbo from 2020 to 2050.Figure 8. RER future projections at Ningbo from 2020 to 2050.

Figure
Figure 9a also shows that more than half of Jiangbei and Beilun may endure a 2~10 h inundation, and the coast of Fenghua, Ninghai and Xiangshan will undergo lesser flooding.Table6shows the area corresponding to each flooding time and each inundation depth interval.The red areas most likely to be affected by storm surge are about 503 km 2 .In the case of a 10 m inundation, the total area affected by the storm surges reaches about 2104 km 2 , covering 21% of the study area, and the areas under the inundation of 8-14 h are all greater than 200 km 2 .

Figure 10 .
Figure 10.Future RER scenarios at Ningbo from 2020 to 2050 that include the influence of 10 m stormsurge inundation.

Figure 10 .
Figure 10.Future RER scenarios at Ningbo from 2020 to 2050 that include the influence of 10 m storm-surge inundation.

1 .
Our Contribution towards RER Assessment Spatially explicit assessment and projection of RER are challenging because RER is affected by many factors such as human impacts and global change.We developed a PCA-based approach to evaluate coastal RER patterns considering spatial factors and applied a bottom-up FLUS model to predict future RER scenarios.The historical evaluation (1975-2015) and the future scenario projection (2015-2050) have proved the effectiveness of both the proposed PCA assessment method and the FLUS prediction model.

Table 1 .
Spatial data and satellite imagery used in this study.

Table 2 .
Selected spatial influencing factors for assessing regional environmental risk (RER) patterns in Ningbo.

Table 3 .
A summary of the indicators of the RER assessment and projection methods.

Table 4 .
Accuracy and error of future land use simulation (FLUS) simulation of 2015 RER.

Table 5 .
Percentage area for each RER level through time.

Table 6 .
Flooding hour and area related to each inundation depth interval.

Table 7 .
RER levels and their area considering the influence of storm surges.

Table 6 .
Flooding hour and area related to each inundation depth interval.

Table 7 .
RER levels and their area considering the influence of storm surges.

Table 8 .
Percentage area at each RER level considering the influence of storm surge.