An Integrated Approach for the Simulation Modeling and Risk Assessment of Coastal Flooding

The evaluation of storm surge flood risk is vital to disaster management and planning at national, regional and local levels, particularly in coastal areas that are affected more severely by storm surges. The purpose of this paper is to propose a new method that includes two modules for the simulation modeling and risk assessment of coastal flooding. One is a hydrodynamic module for simulating the process of the flood inundation coastal inundation arising from storm surge, which is based on a cellular automata (CA) model. The other is a risk assessment module for quantitatively estimating the economic loss by using the inundation data and land use data. The coastal areas of Pearl River estuary in China were taken as a case study. Simulation results are compared to experimental results from MIKE 21 and depth data from a social-media-based dataset, which demonstrates the effectiveness of the CA model. By analyzing flood risk, the flood area and the direct economic losses predicted are close to the actual case incurred, further demonstrating the computational reliability of the proposed method. Additionally, an automatic risk assessment platform is designed by integrating the two modules in a Geographic Information System (GIS) framework, facilitating a more efficient and faster simulation of coastal flooding. The platform can provide the governments as well as citizens of coastal areas with user-friendly, real-time graphics for coastal flood disaster preparation, warning, response and recovery.


Introduction
A storm surge disaster is one of the most frequent and devastating natural disasters for human society, even surpassing earthquakes and tsunamis in terms of casualties and damage, and it has become one of the most serious constraints on the economic development of the coastal areas around the world [1]. As a result of global climate change and rapid urbanization, coastal areas are now more vulnerable to extreme events, such as typhoons, storm surges and floods [2,3]. Recent studies have indicated that as coastal zones suffer stronger storm surges, their risk of flooding increases significantly [4][5][6]. The quality of storm-tide inundation estimates has a tremendous influence on coastal cites' flood management [7]. Therefore, it is of vital significance to carry out coastal flooding simulation and risk analysis.
In the past decades, plenty of models have been widely used in the development and application of numerical simulation and Geographic Information System (GIS)-based analysis for simulating coastal flooding. Generally, current models are either static or dynamic in terms of modeling coastal flooding. Existing static models make simple assumptions that inundation occurs when the land area elevation is below the maximum water surface elevation and hydraulically connected to open water or submerged areas [8]. Due to its simplicity and computational efficiency, static methods have been widely used to

Case Background
Typhoon Hato (1713) storm surge in the coastal areas of Pearl River estuary in China was taken as a real event study. The Pearl River estuary lies in the Pearl River, the second largest river in China. The estuary is densely covered with river networks, and the Pearl River area is densely populated and economically developed compared to most other parts of China. The region is hit by typhoons frequently, so storm surges bring significant economic losses to the area. Typhoon Hato landed on the coast of Jinwan District, Zhuhai City, Guangdong Province around 12:50 on 23 August 2017. This has been the strongest typhoon landing in the Pearl River estuary since 1965 [35] in terms of the maximum tidal elevation of 2.79 m observed by the Zhuhai station. The typhoon track drawn according to the data provided by the National Meteorological Center is shown in Figure 1.
Water 2020, 12, x FOR PEER REVIEW 3 of 18

Case Background
Typhoon Hato (1713) storm surge in the coastal areas of Pearl River estuary in China was taken as a real event study. The Pearl River estuary lies in the Pearl River, the second largest river in China. The estuary is densely covered with river networks, and the Pearl River area is densely populated and economically developed compared to most other parts of China. The region is hit by typhoons frequently, so storm surges bring significant economic losses to the area. Typhoon Hato landed on the coast of Jinwan District, Zhuhai City, Guangdong Province around 12:50 on 23 August 2017. This has been the strongest typhoon landing in the Pearl River estuary since 1965 [35] in terms of the maximum tidal elevation of 2.79 m observed by the Zhuhai station. The typhoon track drawn according to the data provided by the National Meteorological Center is shown in Figure 1. The largest storm surge occurred in Xiangzhou District, Zhuhai City, in a way of spinning counterclockwise in the northern hemisphere, and the area on the right side of the typhoon track is most seriously affected. Therefore, this paper selected the coastal areas of Macao and Xiangzhou District, Zhuhai City (N 22.06°-N 22.30°, E 113.41°-E 113.60°) as the research area, which were more seriously affected. Topographic data were derived from the ASTGTM dataset provided by the Computer Network Information Center, Chinese Academy of Sciences (CNIC) (http://www.gscloud.cn). ArcGIS 10.5 software (maintained by Esri China Information Technology Co. Ltd. in Beijing, China) was used to preprocess the original data to obtain the range and topographic map of the area, see Figure 2.
The Landsat8 satellite remote sensing data LC81220452018139LGN00 provided by the geospatial data cloud was selected for georeferencing and geometric correction. Because the supervised classification has high precision and controllable classification type, it was selected as the extraction method of land use. The software Arcgis10.5 was imported for cutting and data conversion to obtain the land use distribution map, as shown in Figure 3. The largest storm surge occurred in Xiangzhou District, Zhuhai City, in a way of spinning counterclockwise in the northern hemisphere, and the area on the right side of the typhoon track is most seriously affected. Therefore, this paper selected the coastal areas of Macao and Xiangzhou District, Zhuhai City (N 22.06 • -N 22.30 • , E 113.41 • -E 113.60 • ) as the research area, which were more seriously affected. Topographic data were derived from the ASTGTM dataset provided by the Computer Network Information Center, Chinese Academy of Sciences (CNIC) (http://www.gscloud.cn). ArcGIS 10.5 software (maintained by Esri China Information Technology Co. Ltd. in Beijing, China) was used to preprocess the original data to obtain the range and topographic map of the area, see Figure 2.
The Landsat8 satellite remote sensing data LC81220452018139LGN00 provided by the geospatial data cloud was selected for georeferencing and geometric correction. Because the supervised classification has high precision and controllable classification type, it was selected as the extraction method of land use. The software Arcgis10.5 was imported for cutting and data conversion to obtain the land use distribution map, as shown in Figure 3.

Initial and Boundary Conditions for Model Testing
The initial and boundary conditions of model testing include three parts, which are the terrain roughness coefficient expressed as the Manning coefficient, the water-level boundary and the

Initial and Boundary Conditions for Model Testing
The initial and boundary conditions of model testing include three parts, which are the terrain roughness coefficient expressed as the Manning coefficient, the water-level boundary and the definitions of wet and dry. The Manning's roughness coefficient is another important parameter for calculation of coastal flooding. It reflects the influence of the rough surface of the terrain on the flow of water. Roughness values were associated with the type of surface coverings in different areas. Relative to the impact of terrain on the flood, the surface roughness has less effect on the water level. Only a rough estimate is made here. The Manning's roughness coefficient of each land use is estimated as follows: development area 0.017, forest 0.16, cropland 0.035, water 0.03, building area 0.015, and the distribution of land use is shown in Figure 3. Besides, the bottom friction coefficient adopts the uniform value (Nikuradse = 0.001). According to the information broadcast by the National Meteorological Center, the duration of typhoon influence in this area was about 8 h, so the simulated flooding process lasted for 8 h (10: 30-18: 30), with a peak hydrographic boundary condition of 2.79 m. The whole simulation was carried out in a closed area. In addition, in order to avoid that fluid retreating from a wet region produces negative water depth that will cause the model to crash. Dry grid cells are defined as cells where h ≤ 0.05 m.

Methodology
The proposed model is divided into two parts to conduct coastal flooding analysis: (1) a simplified two-dimensional shallow hydrodynamic module that rapidly and accurately simulates the process of the flood inundation; and (2) a risk assessment module that estimates the economic loss by combing the flood data and the land use categories in GIS. In the following subsections, the two parts are discussed in detail.

Simplified Two-Dimensional Shallow Hydrodynamic Module
CA is a grid-based dynamic system, which is discrete in time, space and state, and its spatial interaction and temporal causality are local [26]. It is generally composed of five parts: (1) cells, (2) cellular states, (3) cellular lattice, (4) cellular neighbors, and (5) transformational rules [36].
Traditional triangular mesh generation based on a finite element model can refine and process complex boundaries; however, it is composed of a large amount of data, complex structure and time-consuming calculation, which make it difficult for existing models to couple with other statistical models and carry out a risk analysis in ArcGIS. In this paper, rectangular grid division is used to establish cellular space, which can be easily converted and combined with raster data in spatial structure. The neighborhood structure is the Von Neumann type that consists of the central cell itself and its four adjacent cells, as is shown in Figure 4.
Water 2020, 12, x FOR PEER REVIEW 6 of 18 establish cellular space, which can be easily converted and combined with raster data in spatial structure. The neighborhood structure is the Von Neumann type that consists of the central cell itself and its four adjacent cells, as is shown in Figure 4. In the CA model, the state of the target cell at the next time steps is determined by both the state of itself and its adjacent cells at the current moment [26]. The core part of CA is to define transformation rules [37]. Therefore, to simulate flood routing more efficiently and accurately, this paper discretizes the two-dimensional shallow hydrodynamic equations and introduces it into the CA model as a transformation rule, in which the spatio-temporal variations in flood depth and velocity are governed by local rules in the neighborhood.
Saint-Venant two-dimensional shallow hydrodynamic Equations with the water depth and flow as main control variables are shown in Equations (1)-(3) [38].
Continuity Equations: In the CA model, the state of the target cell at the next time steps is determined by both the state of itself and its adjacent cells at the current moment [26]. The core part of CA is to define transformation rules [37]. Therefore, to simulate flood routing more efficiently and accurately, this paper discretizes the two-dimensional shallow hydrodynamic equations and introduces it into the CA model as a transformation rule, in which the spatio-temporal variations in flood depth and velocity are governed by local rules in the neighborhood.
Saint-Venant two-dimensional shallow hydrodynamic Equations with the water depth and flow as main control variables are shown in Equations (1)-(3) [38].
Continuity Equations: where h is the water depth (m), Z is the water stage (m), M and N are the vertical average unit width discharge in the direction of x and y, respectively (m 2 s −1 ), u and v are the velocity components of the vertical average velocity in the direction of x and y, respectively (ms −1 ),n is Manning's roughness coefficient, t is computing time (s) and g is the gravitational acceleration (m/s 2 ).
In momentum equations, the first term is the acceleration term, the second and third are convection terms, the fourth term is the gravity term, and the fifth is the resistance term. The second and third convection ones are nonlinear in differential equations, which often cause instability in calculations and accounts for a relatively small proportion; thus, they can be neglected [38]. In comparison, because of neglecting the convection term in the momentum equation, this greatly reduces the complexity of the equation. Therefore, Equations (1)-(3) are discretized using finite differences and an explicit scheme for the time dependence to obtain Equations (4)- (6). The cellular space difference diagram is shown in Figure 5.
In Equations (4)-(6), the computational grid of water depth and velocity is staggered in space and time. Solving Equations (4)-(6) simultaneously, we get Equations (7)-(9), where k is the number of cellular simulation changes, ∆t is the size of time step (s), ∆x and ∆y define the cell size (m), and i and j define the cell position. In the beginning, the water depth and the flow velocity should be initialized (w is the initial flow velocity caused by a typhoon (m s −1 )). To meet the convergence and stability conditions, the time step and distance step should be limited, as shown in Equation (10). Therefore, the time step will be adaptive and change during a simulation but is fixed in space at each time step. c is the wave speed (m/s).
Water 2020, 12, 2076 The cellular state at each moment can be obtained by repeatedly using Equations (7)- (9), assuming T is the total simulation time, the calculation process is shown in Figure 3, where the calculation conditions for assessing the water depth are as follows: according to Equations (7)-(9), when solving for M, there is a requirement that at least one water depth state in the x-direction neighborhood of the cell should be greater than zero. Similarly, the calculation condition should be satisfied in the y-direction of the cell.  In Equations (4)- (6), the computational grid of water depth and velocity is staggered in space and time. Solving Equations (4)-(6) simultaneously, we get Equations (7)-(9), where k is the number of cellular simulation changes, Δt is the size of time step (s), Δx and Δ y define the cell size (m), and i and j define the cell position. In the beginning, the water depth and the flow velocity should be initialized ( w is the initial flow velocity caused by a typhoon (ms −1 )). To meet the convergence and stability conditions, the time step and distance step should be limited, as shown in Equation (10). Therefore, the time step will be adaptive and change during a simulation but is fixed in space at each time step. c is the wave speed (m/s). Based on the above analysis, the algorithm mechanism of the CA model of coastal flooding is as follows: the study area is divided into several rectangular cells with the same size, which make up the cellular space. Each cell in the cellular space has a discrete state and neighborhoods. The state of cellular space itself and its neighborhoods at the current time step determine the cellular state at the next time step, while this process is controlled by the transformation rules. The main algorithm implemented is shown in Figure 6.
The whole simulation process is as follows: first, the original data are entered and the variables are initialized; second, the states of all cells in the cell space are calculated and updated according to the conversion rules, a new round of calculation is started after a one-time step is finished, and the process is repeated until the simulation is completed; finally, the information is the output for subsequent statistical risk analysis, such as the inundation extent, water depth, flow velocity and total duration. Based on the above analysis, the algorithm mechanism of the CA model of coastal flooding is as follows: the study area is divided into several rectangular cells with the same size, which make up the cellular space. Each cell in the cellular space has a discrete state and neighborhoods. The state of cellular space itself and its neighborhoods at the current time step determine the cellular state at the next time step, while this process is controlled by the transformation rules. The main algorithm implemented is shown in Figure 6. The whole simulation process is as follows: first, the original data are entered and the variables are initialized; second, the states of all cells in the cell space are calculated and updated according to the conversion rules, a new round of calculation is started after a one-time step is finished, and the process is repeated until the simulation is completed; finally, the information is the output for subsequent statistical risk analysis, such as the inundation extent, water depth, flow velocity and total duration.

Physically-Based and Real-Time Risk Assessment Module
The quantitative assessment of storm surge disaster loss is the basis and premise of storm surge disaster risk management. Based on the impact of disasters on human survival and development goals, losses can be divided into economic losses and non-economic losses [39]. Economic losses can be further divided into direct economic losses and indirect economic losses. Since the range covered by indirect economic losses is extremely wide, it is very difficult to calculate them thoroughly and

Physically-Based and Real-Time Risk Assessment Module
The quantitative assessment of storm surge disaster loss is the basis and premise of storm surge disaster risk management. Based on the impact of disasters on human survival and development goals, losses can be divided into economic losses and non-economic losses [39]. Economic losses can be further divided into direct economic losses and indirect economic losses. Since the range covered by indirect economic losses is extremely wide, it is very difficult to calculate them thoroughly and accurately. Therefore, in this paper, only direct economic losses are taken into account. According to the land use of the whole area under study, the calculation Equation of direct economic losses is [40]: where D d is the total direct economic loss, D ij is the loss value of the i-th disaster-bearing body within the j-th inundation, W ij is the value of the i-th disaster-bearing body within the j-th inundation, η ij is the loss rate of the i-th disaster-bearing body within the j-th inundation, n is the number of land use, and m is the number of inundation classifications. The determination of the storm surge disaster loss rate is the key to calculating the economic losses of storm surge disasters. The magnitude of the disaster loss rate is closely related to the disaster situation, the disaster-resistant ability of the disaster bearing body and the moisture-proof measures in the disaster area. In this paper, the disaster loss rate was determined through overlay analysis by studying the historical data of the study area and the loss rate of storm surge disasters in other areas. Land use mainly includes cropland area, forest area, water area, development area and building area. As stated in related reference work [40] and investigated data, the unit area value can be evaluated for the land use. The three types of cropland, forest and water were mainly considered as the economic losses caused by production. The building area was mainly considered the losses of family property, the destruction of houses and basic public facilities. The development area was mainly considered the destruction of houses and basic public facilities, and finally obtained the unit area value of land use. Table 1 shows the unit area value of each land-use type and the disaster loss rate corresponding to different water depths. In this paper, the unit area value of each land-use type is assumed to follow a normal distribution, with a mean and a standard deviation (SD) ( Table 1). In order to solve the long-standing problem of the integration of hydrodynamic models with large geographic systems. In this paper, by incorporating the flood simulation with geographic statistics and introducing the loss assessment method based on land use, a coupled risk analysis model of coastal flooding was established. To realize the automation of the whole process of storm surge flood disaster analysis, this study adopted a tightly coupled model to integrate the hydrodynamic module and the risk assessment module into the same GIS framework. Thereout, we develop a platform that integrates GIS with hydrological modeling, where the spatiotemporal simulation is integrated with loss analysis of coastal flooding. The platform can intuitively display the flooding process; meanwhile, we can obtain real-time water depth, flow velocity, inundation extent, loss distribution data at different times and so on.

Calibration and Testing
The hydraulic model used for comparison in this article is the HD module in MIKE 21, which is the most well-accepted software for assessing flood hazards in coastal areas [22][23][24]. It uses a two-dimensional hydrodynamic numerical model with an average vertical velocity assumption, for the changes of hydraulic elements such as water depth and flow velocity in the vertical direction are much smaller than the changes in the horizontal direction during the flooding process. It included complex models based on discretized, numerical solutions of general laws (such as the conservation of energy and momentum). The principle has been the idea that sufficient complexity in the model rules was necessary to mimic the complexity of actual flood dynamics. This article intended to discuss from a different perspective, that is spatial complexity can arise in a model with simple rules as long as the spatial interactions of these rules are allowed to be sufficiently complex.
To verify the accuracy and reliability, this paper compared the simulation results of the CA model and MIKE 21 model. Using the data from a real typhoon case, metrics were used to assess the degree of consistency and variation between model predictions, including the inundation extent, two statistics (Root Mean Square Error (RMSE) and Correlation Coefficient (R)) of water-depth data and a comparison of photos of typical sample points.
The RMSE is calculated using Equation (12), where W ca i,j represents the water level simulated by CA at the cell C i,j , and W mike i,j represents the water level simulated by MIKE 21 at the cell C i,j . The R statistic is calculated using Equation (13), where W ca i,j represents the mean value of water level simulated by CA at the cell C i,j , and W mike i,j represents the mean value of water level simulated by MIKE 21 at the cell C i,j .
According to Wang's research results [41], 300 valid photos were collected from three social media of Weibo, WeChat and Baidu. The photo shooting time, location and water-level information were determined according to the webpage description information and photo display, of which 150 were used for parameter adjustment (MIKE 21 and CA) and 150 for model verification (CA).

Comparison with MIKE 21 Simulation Results
The simulation was carried out by using the HD module of MIKE 21, and the model parameters were calibrated to achieve the best simulation effect with the collected 150 social media photos. We adopted the same computational conditions setting for comparing our modeling results and demonstrating the accuracy and efficiency of the proposed model.
The hydraulic elements can be expressed by an average value along the water depth only when the vertical scale of the grid is much smaller than the horizontal scale. The optimal grid size should be between 10 and 200 m. If it is too small, the above conditions cannot be met. If it is too large, the accuracy would be lost. So, the 50 m grid size was taken as an example to illustrate the simulation effect of the CA model. Figures 7-9 show the predicted water depths and inundation extent at three time points for both models and the absolute differences in water depth at these times. The inundation extent differs only marginally, and the Root Mean Square Error is only 0.03 m whilst the Correlation Coefficient is 0.95 (see Table 2). The maximum differences in water depth are ±0.1 m, but these are very localized and, over the majority of the domain, water depth differences are small. This demonstrates that the CA model is capable of reproducing results that are in agreement with the ones obtained from the MIKE 21. The reason for the discrepancy may be that MIKE 21 computed the full shallow water Equation while CA neglected the convective term, as errors may occur, especially in areas with sudden changes in altitude.      To test whether the model was sensitive to grid size, identical simulations were also run with the MIKE 21 and CA model at ∆x = 10, 25, 50, 100 and 200 m ( Table 2). It is known from Table 2 that the RMSE of the two models is below 0.05, which demonstrates that the inundation extent and water depth simulated by the CA model are spatially consistent with MIKE 21 results. Besides, their correlation coefficients are higher than 0.93, which illustrates that there is a high correlation and low anomaly between the two model simulation results, which is also reflected in Figures 7-9. In summary, the CA model performs well in the study area; in addition, the results are not greatly affected by the grid variation in the defined grid range as the extent of the study area is wide and the terrain slopes gently.
Runtime is an important factor for evaluating the efficiency of the CA model. In this regard, the CA model challenges MIKE 21 in simulating the research area at ∆x = 10, 25, 50, 100 and 200 m (see Table 3). In this comparative analysis, the simulation time of CA was shorter than MIKE 21 by more than 60% for all grid sizes. Furthermore, it needs to be recognized that MIKE 21 needs to carry out mesh processing and define a series of boundary conditions in the early stage, which increases the complexity and workload of the model simulation. The CA model can automatically generate cellular space and intelligently set partial parameters, which helps to realize simple and rapid evaluation of coastal flooding.

Comparison of CA Model with Social Media-Based Dataset
To evaluate the modeling with the actual water depth, 150 photos were collected from Weibo, WeChat and Baidu for further verification. The water depths were obtained from those photos for a comparison with the CA model at ∆x = 10, 25, 50, 100 and 200 m (see Table 4). The results show that errors are low for all grid sizes, but the CA model was marginally better at ∆x = 10, 25 and 50 m than ∆x = 100 and 200 m, where the solution quality is dominated by the effect of the grid size. To illustrate in detail, five typical location points were selected for a comparison of water depth whose locations are shown in Figures 10 and 11      According to the photos collected from the typical locations and related reports, it was found that when t = 105 min, the actual water depth of the Dama Road was about 1.1-1.4 m, and the water depth simulated by CA was 1.29 m; when t = 330 min, the actual water depth of the Honggang Harbor was around 0.5-0.8 m, and the water depth simulated by CA was 0.70 m (see Figure 12). Similarly, the water depth in other locations was consistent with the actual situation, illustrating the reliability and accuracy of the CA model. These findings confirm that the CA model can provide a real-time reference in terms of flooding and water depth data for the management of storm surge disasters. According to the photos collected from the typical locations and related reports, it was found that when t = 105 min, the actual water depth of the Dama Road was about 1.1-1.4 m, and the water depth simulated by CA was 1.29 m; when t = 330 min, the actual water depth of the Honggang Harbor was around 0.5-0.8 m, and the water depth simulated by CA was 0.70 m (see Figure 12). Similarly, the water depth in other locations was consistent with the actual situation, illustrating the reliability

Coastal Flood Risk Assessment
According to spatial geographic relationships, the maximum water depth data was coupled with land use to preliminarily evaluate the economic losses. The distribution of flood disaster losses in the study area was obtained, and the diagram in the platform interface is shown in Figure 13. The inundation extent and economic losses of land use corresponding to different water depths were counted and shown in Table 5. According to the photos collected from the typical locations and related reports, it was found that when t = 105 min, the actual water depth of the Dama Road was about 1.1-1.4 m, and the water depth simulated by CA was 1.29 m; when t = 330 min, the actual water depth of the Honggang Harbor was around 0.5-0.8 m, and the water depth simulated by CA was 0.70 m (see Figure 12). Similarly, the water depth in other locations was consistent with the actual situation, illustrating the reliability and accuracy of the CA model. These findings confirm that the CA model can provide a real-time reference in terms of flooding and water depth data for the management of storm surge disasters.

Coastal Flood Risk Assessment
According to spatial geographic relationships, the maximum water depth data was coupled with land use to preliminarily evaluate the economic losses. The distribution of flood disaster losses in the study area was obtained, and the diagram in the platform interface is shown in Figure 13. The inundation extent and economic losses of land use corresponding to different water depths were counted and shown in Table 5.  We added up the product of the loss rate and the inundation area for each land use type, and the five sums values were considered as a coefficient in normal distribution functions. Based on the assumption in Section 2.3.2, the probability density distribution of loss for each land use type was calculated. Loss for each land use type is an independent random variable that is normally distributed. Their sum is also normally distributed, which can be derived from Equations (14)-(16) (The sum of two independent normally distributed random variables is normal, with its mean being the sum of the two means, and its variance being the sum of the two variances.) Therefore, the total direct economic loss distribution function and the two relevant parameters (mean loss: CNY 3.17 billion, 95% confidence interval: CNY 1.03~5.32 billion) can be calculated. The actual direct economic loss was CNY 5.15 billion, as reported by the 2017 China Marine Disaster Bulletin [39], which indicates that the model could assess the loss of disasters, especially within an order of magnitude of the error. In the future, we will expand the database and optimize parameters to improve the model.
According to the summary of the data presented in Table 3, the total inundation extent of the study area was 66.63 km 2 . From Figure 13, we can see that the distribution of disaster losses shows a certain regularity: the disaster losses were higher in areas with higher value of per unit area, offshore areas and low marine topography, which reflects the main role of topography, location and unit area value in storm surge flood loss analysis. The order of the losses for each land-use type is as follows: building area > development area > forest > water > cropland. It is obvious that the largest proportion of disaster loss was in the building area. This is mainly because most of the building areas with high investment are located in the low coastal areas, which was severely affected by a storm surge. It can be inferred from Figure 13 and Table 5 that the real-time storm surge monitoring system and emergency plan system in building areas should be improved and strengthened in the process of disaster prevention and mitigation, and that the dike construction along the coastal areas should also be strengthened.

Conclusions
This study outlines a new model that integrates a hydrodynamic module and a risk assessment module for flood simulation and risk assessment in the coastal areas of the Pearl River estuary, China. The hydrodynamic module that based on the CA model, the effectiveness of which has been demonstrated by comparing the simulation results with MIKE 21 and a social-media-based dataset. Then, the simulated water depth and land-use data are utilized to estimate the storm surge risk that was visualized on the integrated platform.
The CA model in this study well simulates the historical storm events with good confidence and accuracy when comparing with MIKE 21 and a social-media-based dataset. In terms of the inundation data at three different times, the inundation and difference of the two models were drawn at a grid size of 50 m. Moreover, the RMSE (≤0.05 m) and R (≥0.93) were calculated at ∆x = 10, 25, 50, 100 and 200 m. The results show a little discrepancy of two models but that is reasonable. Similarly, the water depth obtained from the CA model was consistent with depth data from a social-media-based dataset at five grid sizes. The results clearly illustrate that the CA model performed well for the entire study area. At the same time, the simulation time of the CA model was over 60% shorter than that of MIKE 21, achieving a fast computational performance with high accuracy.
Analyzing flood risk, the results indicate that the total inundation extent of the study area was 66.63 km 2 , and the direct economic losses were 11.40 billion CNY, which is very close to the actual losses incurred and thus demonstrates computational reliability. They reflect the main role of topography, location and unit area value in storm surge flood risk analysis. Moreover, based on the GIS framework, this study integrated the hydrodynamic module with the statistical risk analysis model with a tight coupling, designing and developing a data-driven automatic assessment platform for coastal flood risk. It realized a refined risk assessment of coastal flood and will provide strong theoretical guidance and technical support for future storm surge disaster risk analysis and the development of disaster prevention measures.
The CA model established in this paper adopts the same cellular unit and neighborhood definition in the whole simulation process. Future work could improve the model to automatically adjust the grid division and select the neighborhood types according to the terrain complexity, and introduce a data assimilation algorithm to calibrate the model parameters. The spatio-temporal dynamic evolution of geographic phenomena and spatial heterogeneity of velocity can be introduced into the