Enhancing an Analysis Method of Compound Flooding in Coastal Areas by Linking Flow Simulation Models of Coasts and Watershed

Rainfall, tide, and wave overtopping can together cause inundation in some coastal areas. However, there have been few cases of inundation analysis considering compound flooding due to the difficulty of the interdisciplinary technique. The purpose of this study is to suggest an analysis technique of compound flooding in coastal areas that links flow simulation tools of coasts and watershed. XP-SWMM was used to analyze rainfall–runoff simulation and 2D surface inundation analysis. A coupled model of ADCIRC and SWAN (ADCSWAN) was used for analyzing expected flow, including wave spectrum, tide, and storm surge from the distant ocean, and the FLOW-3D model was used to estimate the flow and wave overtopping in the coastal area. The analysis technique was applied to Marine City located in Busan, the Republic of Korea, where rainfall and wave overtopping occurred simultaneously during typhoons, and the results were reasonable. The study results are meaningful as they can contribute to the improvement of reproducibility of real phenomena of compound flooding in coastal areas through linked application of flow simulation tools of coasts and watershed.


Introduction
Many urban areas with various infrastructures are located in coastal areas, as the land of the Republic of Korea is situated in the Korean Peninsula and is mountainous [1]. Climate change due to recent global warming raises the risk of disasters occurring in those coastal regions. These disasters may be a result of a rise in the average seawater level and increases in rainfall intensity and typhoon strength [2]. Under these circumstances, the development of coastal areas through reclamation has increased the risk of flooding, which has caused loss of life and damage to property [3].
Coastal areas can be flooded by compound factors, such as rainfall, rises in seawater level at high tide and wave overtopping by storm surge. The phenomenon was defined as compound flooding (CF), which is an extreme event taking place in low-lying coastal areas as a result of the interaction between co-occurring high sea levels and large amounts of rainfall [4]. When the two hazards interact, the resulting impact can be worse than when they occur individually [5]. In fact, many areas were flooded as the tide of Masan Bay increased more than 2 m above the forecast during Typhoon Maemi in Korea in 2003. As a result, the disaster caused significant loss of life and property [6]. The heavy rainfall combined with high tide in Busan on July 7, 2009 caused poor drainage, and the flood damage was aggravated [7]. Typhoon Chaba in 2016 also caused massive flooding in Marine City in Haeundae-gu, Busan, due to the rainfall and wave overtopping which was caused by storm surge [1]. Therefore, these compound external forces must be included for reasonable simulations of flooding in the above coastal areas. However, there is a critical knowledge gap that exists in the attempt to address the effects of compound flooding due to the combined effects of different flood sources. Accordingly, these processes have been studied separately using different types of models [8].
The previous studies in hydrologic field related to inundation analysis in coastal areas were conducted considering only rainfall and tides as boundary conditions. For example, Pokharel et al. [9] analyzed an inundation situation of Nagoya Bay in Japan for the heavy rainfall and seawater level rise due to global warming by a synthetic analysis model that includes a sea model, an inundation model, a river model and a sewer system model. However, wave overtopping was not considered in this study. Similarly, Ray et al. [10] analyzed the combined effect of storm surge and inland rainfall on the floodplain of a coastal bayou in the Houston area by using Hydrologic Engineering Center's (HEC) hydrologic modeling system (HEC-HMS) and river analysis system (HEC-RAS), excluding wave overtopping. Karamouz et al. [11] analyzed floodplain delineation in the coastal area of the Bronx River in New York City by using a storm water management model (SWMM) to simulate the runoff from rainfall and the HEC-RAS. Choi and Jun [12] conducted flooding analysis considering the rainfall and seawater level rise due to storm surge by using XP-SWMM.
Meanwhile, most of the coastal engineering studies were carried out, analyzing coastal inundation using physical phenomena from the sea. Leroy et al. [13] proposed a methodology to simulate coastal flooding by wave overtopping with spatio-temporal variability in an urban area with a very high-resolution digital elevation model (DEM). They simulated a flood event induced by overtopping during the Johanna storm in the village of Gâvres in France using the SURF-WB (two-dimensional SURFace water flow with Well-Balanced scheme) model. Gallien [14] estimated maximum wave run-up elevation and overtopping volumes using numerical (XBeach) and empirical (EurOtop) models to validate the coastal flood model. Suh and Kim [15] studied inundation analysis, calculating discharge of wave overtopping by incorporating the empirical formulas of the EurOtop module into the storm surge models of ADCIRC (ADvanced CIRCulation model) and SWAN (Simulating WAves Nearshore) during Typhoon Chaba in Marine City, Busan, Korea. O'Grady et al. [16] demonstrated the combined effect of overtopping with a rising surge-tide inundation and highlighted the importance of wave overtopping to inundation by using the numerical model, that is, the shallow water integrated flood tool (SWIFT).
Recently, studies have been conducted considering both rainfall-runoff and wave overtopping. Song et al. [1] compared an inundation area with or without application of the wave overtopping in Suyoung Bay, Busan, the Republic of Korea by XP-SWMM. However, there was a limit in the analysis method of this study in that wave overtopping discharge was calculated by using the dimensionless discharge formula [17]. As a case similar to the current study, Chen et al. [18] developed a coupled model for predicting and analyzing inundation in coastal nuclear power plants (CNPP) in China under the combined action of extreme rainfall, wave overtopping and seawater level based on TELEMAC-2D and SWMM. The focus of this study was to present an imaginary extreme scenario causing the most dangerous situation to the CNPP. In other words, they did not consider the actual physical phenomena for each of the external forces, and the results of the inundation analysis were not real but possible.
Most of the previous studies of coastal inundation analysis focusing on coastal areas took account of a single external force condition corresponding to their own research field. In the cases of studies considering various external force conditions, some external force conditions were incorporated by some simple empirical formulas rather than the physical mechanism. In other words, few studies realized the actual phenomena for the inundation analysis of the coastal areas by physically configuring all external force conditions. This is because most studies of inundation analysis focusing on the coast were conducted independently by hydrological engineers or coastal engineers. The current study is different from previous studies in that it was conducted by hydrological engineers in cooperation with coastal engineers to analyze real phenomena by using physical models for hydrological and coastal engineering.
XP-SWMM, ADCIRC, SWAN and FLOW-3D models were used to simulate inundation analysis due to complex factors in the coastal area. The adequacy of analysis was validated by applying the models to the area where the inundation damage occurred.

Research Area
The research area is the Marine City basin located in Busan, the Republic of Korea (Figure 1). The area is a residential-oriented town area that was built on a part of the land reclamation area of Suyoung Bay. hydrological engineers in cooperation with coastal engineers to analyze real phenomena by using physical models for hydrological and coastal engineering. XP-SWMM, ADCIRC, SWAN and FLOW-3D models were used to simulate inundation analysis due to complex factors in the coastal area. The adequacy of analysis was validated by applying the models to the area where the inundation damage occurred.

Research Area
The research area is the Marine City basin located in Busan, the Republic of Korea ( Figure 1). The area is a residential-oriented town area that was built on a part of the land reclamation area of Suyoung Bay. Although the major buildings are adjacent to the coastline, active installation of prevention facilities against wave overtopping has been restricted due to the right of the local residents to view the sea. Because of these regional characteristics, many roads and shopping districts in Marine City were inundated due to overflow of seawater caused by wave overtopping during Typhoon Chaba in 2016 ( Figure 2) and Typhoon Kong-Rai in 2018. That is, Marine City frequently flooded when typhoons accompanied by heavy rainfall and storm surge hit the area.

Characteristics of Watershed and Storm Sewer Networks
The watersheds and storm sewer networks of the study area were analyzed by collecting the digital topographic map, planned floor plan for stormwater facilities, etc. The analyzed results were also checked through a field investigation. The total watershed area of Marine City is 0.53 ㎢ and the  Although the major buildings are adjacent to the coastline, active installation of prevention facilities against wave overtopping has been restricted due to the right of the local residents to view the sea. Because of these regional characteristics, many roads and shopping districts in Marine City were inundated due to overflow of seawater caused by wave overtopping during Typhoon Chaba in 2016 ( Figure 2) and Typhoon Kong-Rai in 2018. That is, Marine City frequently flooded when typhoons accompanied by heavy rainfall and storm surge hit the area.

Characteristics of Watershed and Storm Sewer Networks
The watersheds and storm sewer networks of the study area were analyzed by collecting the digital topographic map, planned floor plan for stormwater facilities, etc. The analyzed results were also checked through a field investigation. The total watershed area of Marine City is 0.53 km 2 ; and the region is a typical urban area consisting mostly of residential areas, commercial areas and roads, with the exclusion of some bare lands.
The total watershed was divided into 25 subcatchments, and they were connected by 25 elements of storm sewers. The watershed runoff discharges to eight outfalls. The main drainage system of Marine City watershed is the system with the outfall element of Out 1 as the outlet. All the eight drainage systems directly discharge runoffs from rainfall to the sea. Figure 3 shows the watershed including the subcatchments connected by storm sewers and a schematic diagram of the drainage systems.

ADCSWAN (ADCIRC+SWAN) model
To predict the impact of storm surge caused by typhoons, the prediction model should reproduce the following comprehensive phenomena: sea-level variation induced by typhoons incurring pressure drop and sea surface winds; topographic effects; tides; and waves [19,20]. The coupled model of ADCIRC and SWAN, named as ADCSWAN, was used in this study. The ADCSWAN had various application and verification cases simulating the possible phenomena caused by typhoons [21][22][23][24][25]. The ADCSWAN calculation data became the boundary condition data for the FLOW-3D model.

ADCSWAN (ADCIRC+SWAN) Model
To predict the impact of storm surge caused by typhoons, the prediction model should reproduce the following comprehensive phenomena: sea-level variation induced by typhoons incurring pressure drop and sea surface winds; topographic effects; tides; and waves [19,20]. The coupled model of ADCIRC and SWAN, named as ADCSWAN, was used in this study. The ADCSWAN had various application and verification cases simulating the possible phenomena caused by typhoons [21][22][23][24][25]. The ADCSWAN calculation data became the boundary condition data for the FLOW-3D model.
We used the 2D depth-integrated version of ADCIRC [26] that solves the vertically integrated generalized wave continuity equation (GWCE) and the momentum equation (Equations (1) and (2), respectively) ∂h ∂t where h is the bathymetric depth, t is the time, ζ is the surge elevation above MSL (mean sea level), → U is the depth-averaged horizontal velocity vector, p is the barometric pressure, f is the Coriolis force coefficient,k is a vertical unit vector, τs is the free surface shear stress, τb is the bottom shear stress, ρ is the water density, and g is the gravitational acceleration. The wave model of SWAN follows a fully implicit finite difference method solving the wave action balance equation [27]. SWAN computes the wave action density spectrum N(x, t, σ, θ) in geographical space x and time t, with σ the relative frequency and θ the wave direction, as governed by the action balance equation [28].
where N is the wave action density spectrum, → C g is the wave group velocity, → U is the ambient current vector, C σ is the propagation velocity and C θ is the wave propagation velocity. The source term, S tot represents all physical processes that generate, dissipate or redistribute wave energy. Detailed information regarding governing equations related to ADCIRC and SWAN can be found in [26,28].
ADCIRC uses the SWAN results to extrapolate forward the wave forcing in time. After completion of the coupled time interval, ADCIRC passes on the wind velocity, water levels and currents to SWAN to further update the radiation stress. In this way, the radiation stress gradients used by ADCIRC and the wind speeds, water levels and currents used by SWAN are mutually exchanged between both models [25]. Further details about the coupled ADCSWAN model can be found in [22].

FLOW-3D Model
The FLOW-3D model solves three-dimensional Navier-Stokes equations and is a non-hydrostatic model to simulate high-speed fluids with excellence considering the advection to the vertical direction. The model is a full three-dimensional model; it may simulate complicated free surface water motion such as scattering and wave breaking by adopting the numerical scheme of the volume of fluid [29]. Here, in order to express the free surface, the surface shape of the liquid in contact with the gas is defined as the VOF (volume of fluid) function F(x, y, z, t) as in Equation (4), and F is per unit volume when it contains the free surface. It represents the volume ratio occupied by the fluid.
where, F SOR is the volume change rate of the fluid due to the generation and disappearance of mass, and the turbulence diffusion term F DIF is as follows (Equation (5)): where, V F is the turbulent diffusion coefficient. Turbulent diffusion coefficients can be selectively simulated vortices near structures by employing the Prandtl mixing length model, the one-equation, the two-equation k-ε, the renormalization group (RNG) k-ε model, and large eddy simulation (LES). The RNG k-ε model was applied in this study [30]. The RNG k-ε model is a method of adjusting the coefficients of the existing k-ε model through statistical techniques, removing the motion of the fluid with small-scale turbulent motion, and influencing the effects of large-scale fluid motion and viscosity. It is expressed using coefficient adjustment and a modified viscous force.
Notably, the model incorporates flows, waves, and irregular waves considering the two-dimensional wave spectrum [31], which can yield reasonable simulation results of storm surge heights and storm wave generations during a typhoon event and dynamic exercises of sea water [32,33].

XP-SWMM
Simulations of inundation in an urban area require an urban runoff model coupled with a hydrodynamic model. The urban runoff model must be able to simulate water flows from subcatchments and through storm sewers. The necessary functions of the urban runoff model are overflow computation out of storm sewer nodes and upstream propagation of water through storm sewers. The hydrodynamic model must be able to simulate inundation flow induced by water wave overtopping and overflow from storm sewer nodes.
XP-SWMM is a fully, dynamic hydraulic and hydrologic modeling software that combines 1D calculations for upstream to downstream flow with 2D overland flow calculations [34]. XP-SWMM, developed by XP Solutions Inc. [35][36][37], is highly utilized for urban flooding simulations around the world [38][39][40][41]. Therefore, XP-SWMM was used to analyze runoff of watersheds and inundation in coastal areas in this study. Notably, XP-SWMM can input external flows to nodes and thereby consider the wave overtopping calculated by the FLOW-3D model.

Connections Among the Models
The hydrostatic approximation of the ADCSWAN model needs the experimental equation suggested by a manual [15,42] to simulate wave overtopping with low computational burdens. The FLOW-3D model, however, directly simulates wave overtopping over shore protection structures or breakwaters without an experimental equation, allowing high computational burdens. Hence, the current study applied the ADCSWAN model to the areas of distant sea and the FLOW-3D model to calculate wave overtopping rates around the coastal areas.
Meanwhile, XP-SWMM is also used to analyze rainfall-runoff of watersheds and ground surface inundation in coastal areas. Figure 4 represents the connection relations of the models used in compound flooding analysis in coastal areas and the results calculated by each model. coastal areas in this study. Notably, XP-SWMM can input external flows to nodes and thereby consider the wave overtopping calculated by the FLOW-3D model.

Connections among the Models
The hydrostatic approximation of the ADCSWAN model needs the experimental equation suggested by a manual [15,42] to simulate wave overtopping with low computational burdens. The FLOW-3D model, however, directly simulates wave overtopping over shore protection structures or breakwaters without an experimental equation, allowing high computational burdens. Hence, the current study applied the ADCSWAN model to the areas of distant sea and the FLOW-3D model to calculate wave overtopping rates around the coastal areas.
Meanwhile, XP-SWMM is also used to analyze rainfall-runoff of watersheds and ground surface inundation in coastal areas. Figure 4 represents the connection relations of the models used in compound flooding analysis in coastal areas and the results calculated by each model.

Estimation of the Wave Overtopping Rate
The storm surge and wave models were built using the ADCSWAN model. Table 1 provides important information on how to simulate the ADCSWAN model. The depth of the model used

Estimation of the Wave Overtopping Rate
The storm surge and wave models were built using the ADCSWAN model. Table 1 provides important information on how to simulate the ADCSWAN model. The depth of the model used KHOA (Korea Hydrographic and Oceanographic Agency, http://www.khoa.go.kr/eng/) and GEBCO (General Bathymetric Chart of the Oceans) data, and air pressure and wind used the typhoon best track information of KMA (Korea Meteorological Administration, http://www.kma.go.kr/eng/index.jsp). For the sea surface wind, the Holland gradient wind vortex algorithm was applied, and the calculation area was set wide enough considering the effect of typhoons, and the grid resolution was composed of a triangular grid of at least 30 m with a maximum of 30 km (See Table 1). Figure 5 shows the grid and water depth map of the storm surge model constructed by the ADCSWAN model. Figures 6 and 7 shows the analysis results of the tidal level and storm surge height and significant wave height of Typhoon Chaba. The results were compared with the observed data of the tide station in Busan and Haeundae Buoy for verification of the ADCSWAN model.      The wave overtopping rates, which were input data to XP-SWMM, were calculated by the FLOW-3D model. The area of the FLOW-3D model for the target area is 4.8 km north-south and 4.0 km east-west, and the grid resolution is composed of a Cartesian coordinate system with a minimum of 3.0 m and a maximum of 12.0 m. The duration of the experiment was simulated for a total of 5 h based on when the maximum significant wave height occurred. The input boundary of the FLOW-  The wave overtopping rates, which were input data to XP-SWMM, were calculated by the FLOW-3D model. The area of the FLOW-3D model for the target area is 4.8 km north-south and 4.0 km east-west, and the grid resolution is composed of a Cartesian coordinate system with a minimum of 3.0 m and a maximum of 12.0 m. The duration of the experiment was simulated for a total of 5 h based on when the maximum significant wave height occurred. The input boundary of the FLOW-3D model for wave and storm surge was divided into three zones, and the results of the ADCSWAN  The wave overtopping rates, which were input data to XP-SWMM, were calculated by the FLOW-3D model. The area of the FLOW-3D model for the target area is 4.8 km north-south and 4.0 km east-west, and the grid resolution is composed of a Cartesian coordinate system with a minimum of 3.0 m and a maximum of 12.0 m. The duration of the experiment was simulated for a total of 5 h based on when the maximum significant wave height occurred. The input boundary of the FLOW-3D model for wave and storm surge was divided into three zones, and the results of the ADCSWAN model were assigned as boundary conditions to the target area. The wave and storm surge boundary data should be considered in both the south and east directions, but as a result of analyzing the observations before and after the overtopping occurs in the target area, the wave is proceeding in the south. Therefore, for convenience of calculations, boundary values were applied only to the south direction. The number of grids created in the target area where overtopping occurs was 32, and the interval was 30 m along the breakwater in front of Marine City. Each flow rate passing through a grid point was calculated and determined as a wave overtopping rate. Figure 8a represents the area enclosed by the boundary conditioned by the input of the wave and storm surge of the FLOW-3D model; Figure 8b represents the boundary of the output-producing wave overtopping rate. Figure 9 shows the time series of the total wave overtopping rate per hour for

Rainfall-Runoff and the 2D Simulation Model By XP-SWMM
The drainage system of Marine City watershed shown in Figure 3 was modeled by using XP-SWMM ( Figure 10). Table 2 includes major information and options used in the rainfall-runoff simulation by XP-SWMM. The observed rainfall data were collected from the Haeundae automatic weather system (AWS), from which the Marine City watershed is located within 3 ㎞ (Figure 11). The nonlinear storage equation and the dynamic wave model were used for watershed and channel

Rainfall-Runoff and the 2D Simulation Model By XP-SWMM
The drainage system of Marine City watershed shown in Figure 3 was modeled by using XP-SWMM ( Figure 10). Table 2 includes major information and options used in the rainfall-runoff simulation by XP-SWMM. The observed rainfall data were collected from the Haeundae automatic weather system (AWS), from which the Marine City watershed is located within 3 ㎞ (Figure 11).

Rainfall-Runoff and the 2D Simulation Model By XP-SWMM
The drainage system of Marine City watershed shown in Figure 3 was modeled by using XP-SWMM ( Figure 10). Table 2 includes major information and options used in the rainfall-runoff simulation by XP-SWMM. The observed rainfall data were collected from the Haeundae automatic weather system (AWS), from which the Marine City watershed is located within 3 km (Figure 11). The nonlinear storage equation and the dynamic wave model were used for watershed and channel routings, respectively. The Horton equation was applied to calculate infiltration of the ground surface. The time interval of rainfall and calculation time steps of watershed and channel routings were set by considering the sizes of the subcatchments and the lengths of the pipes.  XP-SWMM requires topographical data to analyze two-dimensional surface inundation. A digital terrain model (DTM) was constructed by using a digital map with the scale of 1:5,000 and survey results by a drone, as shown in Figure 10a. The DTM was used to create two-dimensional simulation grids, as shown in Figure 10b. The size of the grid is 3 × 3 m and the number of generated grids is 72,728. The time step of numerical simulation of inundation was set to 1 s for stable calculation.  XP-SWMM requires topographical data to analyze two-dimensional surface inundation. A digital terrain model (DTM) was constructed by using a digital map with the scale of 1:5,000 and survey results by a drone, as shown in Figure 10a. The DTM was used to create two-dimensional simulation grids, as shown in Figure 10b. The size of the grid is 3 × 3 m and the number of generated grids is 72,728. The time step of numerical simulation of inundation was set to 1 s for stable calculation.
On the other hand, a number of junctions were located along the shoreline, as shown in Figure 10b. The 32 junctions were set to accept the water from wave overtopping calculated by the FLOW-3D model.

Results of Inundation Simulation
The inundation simulation model constructed for Marine City was verified during Typhoon Chaba in 2016. The inundation simulation period was from 8 a.m. to 12 a.m. on October 5, 2016. The followings were the rainfall records of the Haeundae AWS for the simulation period: the total amount of rainfall was 59.0 mm; the maximum rainfall intensity was 40.5 mm/hr. Figure 12 shows the simulation results of flooding in Marine City over time during Typhoon Chaba. Wave overtopping started at 09:20, and it was shown that the inundation area increased over time. After the time of maximum wave overtopping, the wave overtopping rate decreased, and the floodwater on the road drained through the storm sewers into the sea, which resulted in a reduction in the flooding area.

Results of Inundation Simulation
The inundation simulation model constructed for Marine City was verified during Typhoon Chaba in 2016. The inundation simulation period was from 8 a.m. to 12 a.m. on October 5, 2016. The followings were the rainfall records of the Haeundae AWS for the simulation period: the total amount of rainfall was 59.0 mm; the maximum rainfall intensity was 40.5 mm/h. Figure 12 shows the simulation results of flooding in Marine City over time during Typhoon Chaba. Wave overtopping started at 09:20, and it was shown that the inundation area increased over time. After the time of maximum wave overtopping, the wave overtopping rate decreased, and the floodwater on the road drained through the storm sewers into the sea, which resulted in a reduction in the flooding area. Figure 13 shows the five collected photographs for the flooding event and the corresponding spots of the maximum flooding area obtained from the simulation results. An area with inundation depth less than 10 cm was not considered as a flooding area. The matching of flooding between the simulation results and photographs reveals the adequacy of the simulation results for compound flooding. A field survey was also conducted for three locations to analyze the inundation results quantitatively. Figure 14 shows the images of the field survey for the sites to estimate maximum inundation depths, and the spots correspond to A, F and G in Figure 13. The analysis accuracy ranged from 88% to 92%, and the errors were less than 10 cm for the maximum inundation depths at the sites (Table 3) Figure 13 shows the five collected photographs for the flooding event and the corresponding spots of the maximum flooding area obtained from the simulation results. An area with inundation depth less than 10 cm was not considered as a flooding area. The matching of flooding between the simulation results and photographs reveals the adequacy of the simulation results for compound flooding. A field survey was also conducted for three locations to analyze the inundation results quantitatively. Figure 14 shows the images of the field survey for the sites to estimate maximum inundation depths, and the spots correspond to A, F and G in Figure 13. The analysis accuracy ranged from 88% to 92%, and the errors were less than 10 cm for the maximum inundation depths at the sites (Table 3).
Meanwhile, Figure 15 is the archived inundation map of the Marine City area provided by the public data portal (https://www.data.go.kr/). Compared with the analysis results, it can be confirmed Meanwhile, Figure 15 is the archived inundation map of the Marine City area provided by the public data portal (https://www.data.go.kr/). Compared with the analysis results, it can be confirmed that the flooding area is quite similar to the inundation map. Specifically, the areas of inundation estimated by the archived inundation map and simulation results were 35,500 m 2 and 34,776 m 2 , respectively; the simulated inundation area was 98.0% coincident with the archived inundation map. That is, the analysis procedures of compound flooding in coastal areas by linking flow simulation models of coasts and watershed were reasonably validated. that the flooding area is quite similar to the inundation map. Specifically, the areas of inundation estimated by the archived inundation map and simulation results were 35,500 m 2 and 34,776 m 2 , respectively; the simulated inundation area was 98.0% coincident with the archived inundation map. That is, the analysis procedures of compound flooding in coastal areas by linking flow simulation models of coasts and watershed were reasonably validated.

Discussion
The compound flooding in coastal areas was analyzed by using the XP-SWMM, ADCIRC, SWAN and FLOW-3D models. As the wave overtopping during Typhoon Chaba was the main cause of the inundation, we were able to accurately simulate such inundation by using only the ADCIRC, SWAN, and FLOW-3D models. (Figure 16).

Discussion
The compound flooding in coastal areas was analyzed by using the XP-SWMM, ADCIRC, SWAN and FLOW-3D models. As the wave overtopping during Typhoon Chaba was the main cause of the inundation, we were able to accurately simulate such inundation by using only the ADCIRC, SWAN, and FLOW-3D models. (Figure 16). However, under conditions where rainfall is relatively large, the flooding phenomenon can be clearly reproduced by reflecting the rainfall-runoff only. The following analyses were performed in cases of increased rainfall under a certain storm surge condition to examine the inundation effect. The return periods of rainfall were 30-year, 50-year and 100-year under the 50-year frequency storm surge yielding wave overtopping and tide. The rainfall duration was 1 h and the time distribution followed a Huff curve for third-quartile storms. The center of rainfall hyetograph was set to coincide with that of the wave overtopping series. Figure 17 represents the analysis results. The higher the return period of rainfall, the wider the inundation area. Hence, one may obtain an enhanced result of flooding of coastal areas by linking flow simulation models of coasts with a watershed runoff model. However, under conditions where rainfall is relatively large, the flooding phenomenon can be clearly reproduced by reflecting the rainfall-runoff only. The following analyses were performed in cases of increased rainfall under a certain storm surge condition to examine the inundation effect. The return periods of rainfall were 30-year, 50-year and 100-year under the 50-year frequency storm surge yielding wave overtopping and tide. The rainfall duration was 1 h and the time distribution followed a Huff curve for third-quartile storms. The center of rainfall hyetograph was set to coincide with that of the wave overtopping series. Figure 17 represents the analysis results. The higher the return period of rainfall, the wider the inundation area. Hence, one may obtain an enhanced result of flooding of coastal areas by linking flow simulation models of coasts with a watershed runoff model.

Conclusions
Rainfall, tide, and wave overtopping can together cause inundations in coastal areas. While analyzing compound flooding has been restricted due to difficulties of interdisciplinary analyses between hydrological and coastal engineering sectors, our study analyzed the connected processes of storm surge, waves, wave overtopping, and rainfall-runoff.
A coupled model of ADCIRC and SWAN (ADCSWAN) was used for analyzing the expected flow, including the elements of wave spectrum, tide, and storm surge from the distant ocean. Then, the FLOW-3D model was used to estimate the flow and wave overtopping rate in the coastal area. Under the boundary conditions of the sea level and wave overtopping rate, the XP-SWMM model yielded the results of the rainfall-runoff and 2D surface inundation simulations. The analysis technique was applied to Marine City located in Busan, the Republic of Korea, where the rainfall and wave overtopping occurred simultaneously during the period of Typhoon Chaba.
We visually presented the state of variation in the flooding area and depth over time from the inundation analysis. Besides, the analysis results were validated by comparing them with the archived inundation map and image materials, namely, newspapers, photos, and videos. As a result, the simulated inundation area was approximately 98% coincident with the archived inundation map, and the maximum inundation depths were about 90% accurate at the three surveyed sites. These results supported the notion that the analysis technique was reasonable.
The procedures in this study are meaningful as they can contribute to the improvement of reproducibility of real phenomena of compound flooding in coastal areas through linked application of flow simulation tools of coasts and watershed.