Development of a Tsunami Inundation Analysis Model for Urban Areas Using a Porous Body Model

To evaluate tsunami hazards with strong locality in urban areas, this study developed a novel tsunami inundation model based on nonlinear shallow water wave equations and a porous body model (PBM). By applying a kinematic boundary condition that includes both porosity and surface permeability of the porous medium, the proposed model could accurately incorporate geometric effects such as the flow anisotropy caused by the distributions of buildings. The proposed PBM demonstrated as good accuracy for the inundation heights around buildings near the coastline as with a conventional three-dimensional simulation with high resolution. In addition, the model showed its capability to reproduce a tsunami’s essential behaviors in urban areas. In particular, the amplification effect of flow velocity along straight roads surrounded by buildings was reasonably reproduced. It can be expected that the present model can become a useful tool to accurately evaluate the tsunami risks in urban areas.


Background
The damage induced in heavily populated areas by the tsunami generated during the 2011 Great East Japan Earthquake was enhanced because the tsunami struck an urban area [1,2].The flow of a tsunami is accelerated due to the contraction effect in urban areas where strong structures are lined up [2,3].The existence of buildings reduces the visibility of evacuees and causes complex tsunami flow that strikes from an unexpected direction (e.g., from behind).Therefore, precise evaluations of the tsunami risks in urban and populated areas caused by urban-type tsunamis are especially important for disaster risk reduction planning.
Tsunami predictions in urban area have been conducted by using several kinds of models [4].For example, Arikawa and Tomita [5,6] coupled a Storm surge and Tsunami simulator in Oceans and Coastal areas (STOC) [7] and a Super Roller Flume for Computer Aided Design of Maritime Structure in 3D (CADMAS-SURF/3D) [8] using a volume of fluid (VOF) method to efficiently simulate tsunamis from the source to run-up in both directions, and Suwa et al. [9] also developed a model integrating a Smoothed Particle Hydrodynamic (SPH) method and a two-dimensional nonlinear shallow water wave (NSWW) model (TUNAMI-N2 [10,11]).These three-dimensional tsunami inundation simulation models with high resolution (0.5 m-2 m) using the K computer, a Japanese supercomputer, are expected to provide detailed calculations not only of the reflections of tsunamis from buildings and the progression and confluence of a tsunami between buildings but also of the hydrodynamic forces acting on buildings because the individual shapes of buildings in urban areas are incorporated into the Geosciences 2018, 8, 12; doi:10.3390/geosciences8010012www.mdpi.com/journal/geosciencescalculations.However, to implement additional, more effective countermeasures against tsunamis, it is important to estimate the human and property damages caused by various assumed tsunami scenarios.To examine various scenarios and areas for pragmatic tsunami risk evaluations over wide areas, a more practical calculation method that can be executed without a supercomputer is required.Meanwhile, the vertically integrated two-dimensional calculation models for the investigation of near-field tsunamis (TUNAMI-N1, -N2 and -N3 [10,11]; MOST [12]; COMCOT [13,14]; NAMIDANCE [15]; TsunAWI [16]; ANUGA [17]), which are based on NSWW theory, have been widely applied for both research [18][19][20][21][22][23][24] and practical purposes.For example, the TUNAMI has been transferred to many countries in the Tsunami Inundation Modeling Exchange (TIME) project [11,25].Most of the models are often solved using a leap-frog method and a finite difference method (FDM) with a structured mesh because of higher scalability and experience stock.Recently, tsunami inundation calculations with a relatively high resolution (approximately 5 m) have been performed on general-purpose computers.However, it is difficult to precisely reproduce the shape of buildings in the model with limited resolution depending on the size of the buildings.It has been noted that representative length scale of a building needs to be divided into eight sections to adequately express the geometric effects of buildings [26].Namely, for a general house, a mesh size of approximately 1 m or less is required in models applying structured constant mesh system.It is obviously difficult to implement such high resolutions and large-scale calculations using a general-purpose computer.Indeed, such investigations are conducted only on a supercomputer; accordingly, Oishi et al. [22] concluded that a resolution higher than 1.67 m is necessary to have a well-converged solution for their simulation of Miyako City.
To apply the 2-D NSWW model with a practical resolution on a sub-grid scale to urban areas, various methods for incorporating the influences of groups of buildings have been proposed.Roughness-type models focusing on the mechanical effects of groups of buildings have been proposed:

•
(Model 1) the roughness map model, which uses Manning's roughness coefficient "n" according to the land-use type [27] (hereafter, Landuse-n model), • (Model 2) the equivalent roughness model (Equivalent-n model), which considers the drag forces acting on buildings and incorporates it in the model as an equivalent roughness [28,29].
The features of roughness-type models were compared by Muhari et al. [19].Since roughness-type models have high practicality, they are often applied to practical tasks.However, it is difficult to properly reproduce the geometric effects of buildings, such as rises in local water level around buildings [19].
The geometric effects of buildings have been modeled in some calculation models with a practical resolution.The models apply three-dimensional building shapes, referred to as the digital surface model (DSM), to topographic data:

•
(Model 4) the equivalent roughness and topographic model (Equivalent-n/Topography model) [29], which is a hybrid model combining the Equivalent-n model with a topographical model.
The differential equations of the coastal-forest model (Model 5) were derived using an assumption that trees exist in the center of the unit volume element, result of which is that the porous media is marked by vertically dependent volume space but not restriction of flow path area in the horizontal direction.The flood model (Model 6) introduced an isotropic unit porous media with a characteristic value of porosity in the fundamental differential equations, which correspond to Equations ( 1) and (2) in Section 2.1.The models by [34,35] use the value based on volume porosity of the cell mesh as the flow path area at edge of that in their process of discretization by FDM, enabling the anisotropy of flow to be numerically reproduce, wherein the building height effect is not considered.Sanders et al. [33] derived an integral form of a porous-type model that can take anisotropy into consideration by combining the Reynolds transport theorem with the conservation within the cell mesh in the finite volume method (FVM).In addition, the influence of shape of the unstructured mesh in the model has been investigated by Kim et al. [36,37].
The features of existing calculation methods employed for urban areas are summarized in Table 1.Not only the high-resolution simulation model but also the porous-type model using the FVM proposed by Sanders et al. [33] would be expected to calculate the geometric effect of buildings more accurately.Meanwhile, various 2-D NSWW numerical models have been applied to tsunami risk evaluations for practical purposes in terms of the available perspectives.However, differential equations based on NSWW theory, which includes the anisotropic geometric effects of building group as shown in Table 1 and which can be directly applied to the FDM, have not been proposed so far.Furthermore, in the current situation, wherein various numerical models exist, it has been increasingly important to understand the similarities and differences of the tsunami hazard results obtained via the various methods and models.

Objectives
In this study, a novel tsunami inundation analysis model using the FDM is developed to evaluate tsunami hazards with strong locality in urban areas, and characteristics of the tsunami hazards simulated by some conventional numerical models are also investigated.

Development of NSWW Theory and Numerical Model Based on a Porous Body Model
In order to consistently consider the anisotropy of building group in the process from theory to numerical discretization, NSWW differential equations based on the porous body model for anisotropic unit porous media proposed by Sha et al. [38] and Sakakiyama and Kajima [39] are derived.The voids in the porous media are continuous.Assuming tsunami inundation of an urban area with impermeable/non-destructive buildings, the equations are implemented to the conventional TUNAMI-N2 using the FDM.

Investigation of Characteristics of Tsunami Hazards Simulated by Numerical Models
Numerical simulations of the tsunami inundation for the high density of buildings in the town of Onagawa, Miyagi Prefecture, Japan, during the 2011 Great East Japan Earthquake are performed to compare to some conventional models.The Onagawa town was one of the affected areas that suffered tremendous damage during the tsunami [40][41][42][43][44][45].The tsunami hazards and damages in this area have been investigated numerically by some researchers using some methods to be available for inundation simulation in urban area [5,6,9,23,24], providing useful data for comparing various methods.However, complex and various phenomena in this area-e.g., destruction of houses/buildings and breakwaters [42][43][44][45], and morphological change [44]-are not included in these models nor in the proposed model.Tsunami hazards might be affected by various phenomena, e.g., destruction of breakwaters [23] and houses/buildings [46,47], floating debris [47][48][49] or morphological change due to tsunami-induced sediment transport [50][51][52][53], although understanding of the mechanism and the influence is limited.These phenomena, as well as the tsunami source, are uncertainties for the numerical simulations.Therefore, discussions of the accuracy of the models modeling specific physics (i.e., geometric effects of building groups) compared with field survey results (including any uncertainties) do not always result in general conclusions.Meanwhile, we focus on investigating the characteristics of various models using simulated tsunami hazards rather than validation of the models with the tsunami including uncertainties.: included; ( ): partly included; blank: not included.3-D NS: three-dimensional, high-resolution Navier-Stokes model (VOF [5,6] and SPH [9]); 2-D NSWW: two-dimensional, high-resolution nonlinear shallow water wave model [22]; Constant-n: constant roughness model; Landuse-n: roughness map model [27]; Equivalent-n: equivalent roughness model [28,29]; Landuse-n/Topography: roughness map and topographic model [18][19][20]23,30]; Equivalent-n/Topography: equivalent roughness and topographic model [29]; Coastal-Forest: coastal forest model [31]; Flood: flood model [FDM: 34,35,FVM: 32,33]; PBM: proposed model based on porous body model.In the coastal forest, flood models and PBM, the Landuse-n is assumed to be the bottom friction.The models with the circle ( ) are applied to tsunami inundation simulations in this study.

Governing Equations
The NSWW differential equations based on a porous body model are derived.For an incompressible fluid in a porous medium (as shown in Figure 1), where the voids are assumed to be continuous in the medium, the three-dimensional continuity and momentum equations [38,39] are expressed as: where γ V and γ i are the volume porosity (γ V = fluid volume space/dxdydz) and components in the i-th direction of the surface permeability (e.g., γ x = fluid surface area/dydz) in the unit element, respectively.The variables t, x i , u i , p, ρ and g i denote the time, components in the i-th direction of the spatial coordinates (x, y, z), components in the i-th direction of the fluid velocity coordinates (u, v, w), pressure, fluid density and gravity (0, 0, −g), respectively.The parameters M i and R i are the components in the i-th direction of the fluid inertia force and drag force, respectively, τ ij is the stress tensor acting on the surface of the unit volume element, and u i and γ V u i (or γ i u i ) are the pore velocity and the Darcy velocity in the porous medium, respectively.
incompressible fluid in a porous medium (as shown in Figure 1), where the voids are assumed to be continuous in the medium, the three-dimensional continuity and momentum equations [38,39] are expressed as: where γ V and γ i are the volume porosity (γ V = fluid volume space/dxdydz) and components in the i-th direction of the surface permeability (e.g., γ x = fluid surface area/dydz) in the unit element, respectively.The variables t, xi, ui, p, ρ and gi denote the time, components in the i-th direction of the spatial coordinates (x, y, z), components in the i-th direction of the fluid velocity coordinates (u, v, w), pressure, fluid density and gravity (0, 0, −g), respectively.The parameters Mi and Ri are the components in the i-th direction of the fluid inertia force and drag force, respectively, τ ij is the stress tensor acting on the surface of the unit volume element, and ui and γ V ui (or γ i ui) are the pore velocity and the Darcy velocity in the porous medium, respectively.Long waves, such as tsunamis propagating in a porous medium, are considered.Assuming excess pore water pressure in the medium to be negligible, the vertical distribution of pressure is expressed as the following hydrostatic approximation: ( ) where p0 and η are the atmospheric pressure and the water surface displacement in the porous medium, respectively.

Kinematic and Dynamic Boundary Conditions
To accurately evaluate the influences of flow path area (anisotropic permeability) and building height as shown in Tab.1, we must construct a theory based on adequate boundary conditions.Hence, we introduce the following kinematic boundary conditions at the water and bottom surfaces considering both the porosity and the surface permeability in the porous medium [54]: Long waves, such as tsunamis propagating in a porous medium, are considered.Assuming excess pore water pressure in the medium to be negligible, the vertical distribution of pressure is expressed as the following hydrostatic approximation: where p 0 and η are the atmospheric pressure and the water surface displacement in the porous medium, respectively.

Kinematic and Dynamic Boundary Conditions
To accurately evaluate the influences of flow path area (anisotropic permeability) and building height as shown in Table 1, we must construct a theory based on adequate boundary conditions.Hence, we introduce the following kinematic boundary conditions at the water and bottom surfaces considering both the porosity and the surface permeability in the porous medium [54]: where a fixed bed and the Darcy velocity are used.The depth-averaged values of γ V and γ i are defined using the integral values of γ V or i as follows: Meanwhile, the dynamic boundary condition at the water surface is p 0 = 0.The following derivation processes for the two-dimensional equations are explained using an x-z vertical plane for brevity.

Integration of the Continuity Equation
Equation ( 1) is integrated vertically from the bottom to the water surface using Equations ( 4)- (6).
Substituting the depth-averaged flow velocity in the x direction U x (x, t) into Equation ( 7), we can obtain the general form of the continuity equation for the shallow water wave equations based on the porous body model as follows: By using the flux in the , where D = η + h is flow depth), Equation ( 8) can be rewritten as: Next, assuming that the temporal variation in the porosity in the unit water column R V is sufficiently small, Equation ( 9) is rewritten as: where the two-dimensional horizontal form of Equation ( 10) is: The assumption imposed on Equations (10) and (11) holds when the vertical variation in the porosity is sufficiently small, i.e., when the temporal variation in the water level is similar to that of a long wave.

Integration of the Momentum Equation
We derive the momentum equation based on the porous body model.We substitute Equation (3) into Equation (2), which is then integrated vertically (see Matsutomi et al. [13] for descriptions of the terms M i , R i and τ ij ), leading to the following: Substituting U(x, t) into Equation (12), we can obtain the general form of the momentum equation of the shallow water wave equations for the porous body model as follows: Equation ( 13) is rewritten into the following form by using the flux Q(x, t): Using the assumption imposed on Equations ( 10) and (11), Equation ( 14) is rewritten as follows: where the two-dimensional horizontal forms of Equation ( 15) are: The proposed theory with anisotropy in porous medium would reproduce not only the local increase of the water level but also the flow acceleration around buildings due to a dam flow and a concentrative flow with exception a jet flow such as three-dimensional and non-hydrostatic flow.Moreover, the temporal variations, which are represented in their general form by Equations ( 8), ( 9), ( 13) and ( 14), in the properties of the porous medium are considered in the derived theory, and thus, an interaction between a solid and fluid due to the destruction of buildings and debris flow could be expected to be simulated by combining this model with a debris flow model [48,52,55], and the drag force coefficients, for example, in cases of debris dam, are also important to evaluate the decay of the energy in flow acceleration.
In this study, we focus on the geometric effects of buildings.Therefore, the drag force terms M i and R i are omitted, whereas the bottom friction due to τ ij is considered, although the terms represent important forces to simulate the hydraulics in urban areas where the effect of buildings suddenly changes spatially.Assuming that the temporal variation in the water level is sufficiently small similar to a long wave and that buildings are impermeable/non-destructive, Equations ( 11), ( 16) and ( 17) are applied as the basic equations in this study: by eliminating ∂R V /∂t, e.g., the equations would overestimate the water level and flow velocity if a tsunami rapidly overflowed the buildings, and tsunamis passing through and staying in the buildings are not considered.The equations are incorporated into TUNAMI-N2 [10,11] using the FDM, and the proposed model is hereafter labeled a PBM in the following sections.
In the proposed PBM, the effects of both the anisotropy of the surface permeability and the building height, which both constitute challenges in conventional models using the FDM, are properly considered by calculating the porosity and surface permeability in the water column based on information regarding the planar distribution of buildings, the building heights and the simulated flow depth.

Establishing the Porosity and Surface Permeability for a Group of Buildings
In this section, the procedure used to establish the porosity and surface permeability for a group of buildings by using the positional relationships between the buildings and the calculation mesh as shown in Figure 2a-d is explained.Since the 2011 disaster, 3-D building data have been compiled by the Geospatial Information Authority of Japan (GSI) to mitigate disaster.The building data employed in this study are polygonal data with median values consisting of the building height for each building, which apply the above ground level (AGL) obtained from DEM and DSM measured by airborne LiDAR surveys [30].The data of common and solid buildings was applied as building type.None of the buildings overlay.
in this study are polygonal data with median values consisting of the building height for each building, which apply the above ground level (AGL) obtained from DEM and DSM measured by airborne LiDAR surveys [30].The data of common and solid buildings was applied as building type.None of the buildings overlay.

Porosity of the Water Column
As shown in Figure 2a-c, we extract the occupancy area of each building An and the building height Hn from the calculation cell mesh, where n is the building ID within the cell mesh.In the case when the number of buildings that exist in the mesh is N, n = 0, 1,…, N, the integral values of the occupancy area and the characteristic building height on the cell mesh are defined as: where Ac (=ΔxΔy) is the area of the cell mesh and Δx and Δy are the lengths of the edge of the cell mesh in the x and y directions, respectively.
The porosity of the water column for every calculation time step V R as shown in Figure 2d is calculated by referencing the flow depth D during each time step:

Horizontal Surface Permeability of the Water Column
As shown in Figure 2a-c, we extract the occupancy length of each building Ln ' and the building height H'n ' within the edge of the cell mesh, where n' is the building ID in the edge of the cell mesh.In the case where the number of buildings that exist in the edge of the cell mesh is N', n' = 0, 1,…, N', the integral values of the occupancy length and the characteristic building height on the edge of the cell mesh are defined as:

Porosity of the Water Column
As shown in Figure 2a-c, we extract the occupancy area of each building A n and the building height H n from the calculation cell mesh, where n is the building ID within the cell mesh.In the case when the number of buildings that exist in the mesh is N, n = 0, 1, . . ., N, the integral values of the occupancy area and the characteristic building height on the cell mesh are defined as: where A c (=∆x∆y) is the area of the cell mesh and ∆x and ∆y are the lengths of the edge of the cell mesh in the x and y directions, respectively.The porosity of the water column for every calculation time step R V as shown in Figure 2d is calculated by referencing the flow depth D during each time step:

Horizontal Surface Permeability of the Water Column
As shown in Figure 2a-c, we extract the occupancy length of each building L n and the building height H n within the edge of the cell mesh, where n is the building ID in the edge of the cell mesh.In the case where the number of buildings that exist in the edge of the cell mesh is N , n = 0, 1, . . ., N , the integral values of the occupancy length and the characteristic building height on the edge of the cell mesh are defined as: The surface permeability values of the water column for every calculation time step R x and R y as shown in Figure 2d are calculated by referencing the flow depth D , which are the values in each edge of the cell mesh during each time step: In this simulation, we use digital elevation model (DEM) data for the topographic data.When A/A c > 0.95, we add the height H to the DEM to produce a DSM.Hence, the porous-type model is not applied to the cell meshes using DSM and is built by adding the porous medium to Landuse-n/Topography.

Differences between Conventional Porous-Type Models and the Proposed Model
Focusing on the geometric effects, the differences between the conventional porous-type models and the proposed PBM in terms of the R V and R x parameters defined in the previous section resulting from the kinematic boundary conditions of Equations ( 4) and ( 5) are clarified.
When R x = 1, the proposed model corresponds to the coastal-forest model introduced by Matsutomi et al. [31], which considers the height of the porous medium.Meanwhile, the PBM is equivalent to the flood model proposed by Miura et al. [34] and MLIT [35], wherein the porosity and the surface permeability are expressed using the following formula: where l and m are the mesh numbers in the xand ydirections, respectively.
In the derivation process of the flood model, the unit volume element is a homogeneous porous medium, and thus, in Equations ( 1), ( 2), ( 4) and ( 5), the porosity and the surface permeability are not distinguishable, namely, γ = γ V = γ i .The surface permeability in the flood model is defined by using the value of porosity as shown in Equations ( 23) or (24), and building height is not considered, as shown in Equation (22).However, since the proposed PBM is based on adequate boundary conditions considering both the porosity and the surface permeability [54], the PBM can more accurately reflect the geometric effects of urban areas on the tsunami inundation simulations.

Background of the Target Domain
The town of Onagawa, Miyagi Prefecture, Japan (shown in Figure 3), was one of the affected areas that suffered tremendous damage during the tsunami generated by the 2011 Great East Japan Earthquake, where the maximum inundation height and run-up height were 18.5 m and 34.7 m, respectively [40,41].Many houses and buildings in the inundation area were washed away.Six concrete structures were damaged (i.e., toppled and destroyed) [24,42,43,45].Furthermore, the breakwater in the mouth of Onagawa Bay was greatly damaged, and the seabed around the opening was eroded by as much as 9 m [44].This characteristic damage constituted enormous devastation.For these reasons, the tsunami hazards and damages in Onagawa have been investigated by many researchers, and thus, substantial amounts of data and knowledge have been accumulated to date.For example, time series of the flow depths and flow velocities in urban areas were obtained from an analysis of video images by Koshimura et al. [45], where the maximum flow depth and velocity in their analysis were approximately 15 m and 7.5 m/s, respectively.The tsunami hazards and disaster mechanisms in Onagawa have been investigated through three-dimensional, super-high-resolution simulations [5,6,9] and through two-dimensional simulations with the Landuse-n and Equivalent-n models [23,24].In this study, the characteristics of inundation in an urban area simulated using various two-dimensional models are examined, and points of concern are also extracted from each model.

Initial Setup of the Tsunami Simulation
The 2011 tsunami is simulated using a nesting grid system consisting of six layers with mesh sizes of 1215 m, 405 m, 135 m, 45 m, 15 m and 5 m, as shown in Figures 4 and 5.The influxes in the boundary of the smaller region are set to the linearly incorporated values from the larger region.The initial tide level is −0.42 m, and the tsunami for three hours after the earthquake is reproduced using the tsunami source model proposed by Satake et al. [56] as shown in Figure 3, including the spatial variability in the fault slip and time of rupture.The advantage to apply the tsunami source model is to be able to compare the numerical result by the proposed PBM with the three-dimensional, super-high-resolution simulations with same source [5,6,9].The area of focus is the town of Onagawa with ∆x = 5 m, as shown in Figure 5a.We compare two roughness-type models (Landuse-n and Landuse-n/Topography) and three porous-type models (the coastal-forest, the flood models using Equation ( 24) [34,35] in Section 4 and PBM).Manning's roughness coefficients corresponding to the land use conditions presented by Kotani et al. [27], which are water bodies, buildings/houses and paved surfaces, are similarly applied to each of the five models.In the four models other than Landuse-n, the DSM defined in Section 3.2 is used as the topographic data.The DSM employed in this study is shown in Figure 5a,b, where the black colors indicate the locations where building heights are incorporated into the DEM. Figure 6 shows the distribution of the occupancy area ratio A/A c , which is used in the porous-type models (i.e., the coastal-forest, the flood models and PBM).For these reasons, the tsunami hazards and damages in Onagawa have been investigated by many researchers, and thus, substantial amounts of data and knowledge have been accumulated to date.For example, time series of the flow depths and flow velocities in urban areas were obtained from an analysis of video images by Koshimura et al. [45], where the maximum flow depth and velocity in their analysis were approximately 15 m and 7.5 m/s, respectively.The tsunami hazards and disaster mechanisms in Onagawa have been investigated through three-dimensional, super-high-resolution simulations [5,6,9] and through two-dimensional simulations with the Landuse-n and Equivalent-n models [23,24].In this study, the characteristics of inundation in an urban area simulated using various two-dimensional models are examined, and points of concern are also extracted from each model.

Initial Setup of the Tsunami Simulation
The 2011 tsunami is simulated using a nesting grid system consisting of six layers with mesh sizes of 1215 m, 405 m, 135 m, 45 m, 15 m and 5 m, as shown in Figure 4 and 5.The influxes in the boundary of the smaller region are set to the linearly incorporated values from the larger region.The initial tide level is −0.42 m, and the tsunami for three hours after the earthquake is reproduced using the tsunami source model proposed by Satake et al. [56] as shown in Figure 3, including the spatial variability in the fault slip and time of rupture.The advantage to apply the tsunami source model is to be able to compare the numerical result by the proposed PBM with the three-dimensional, superhigh-resolution simulations with same source [5,6,9].The area of focus is the town of Onagawa with Δx = 5 m, as shown in Figure 5a.We compare two roughness-type models (Landuse-n and Landusen/Topography) and three porous-type models (the coastal-forest, the flood models using Equation ( 24) [34,35] in section 4 and PBM).Manning's roughness coefficients corresponding to the land use conditions presented by Kotani et al. [27], which are water bodies, buildings/houses and paved surfaces, are similarly applied to each of the five models.In the four models other than Landuse-n, the DSM defined in section 3.2 is used as the topographic data.The DSM employed in this study is shown in Figure 5a,b, where the black colors indicate the locations where building heights are incorporated into the DEM. Figure 6 shows the distribution of the occupancy area ratio A/Ac, which is used in the porous-type models (i.e., the coastal-forest, the flood models and PBM).[24], where the circle indicates the location of Onagawa, Miyagi Prefecture, Japan, and the star is the epicenter of the earthquake.The units of distance are km.
The buildings/houses and the breakwater in the bay are assumed to be intact; thus, the roughness coefficients and the parameters of the unit volume element in the porous medium, namely, A, H, L and H , are temporally constant values.As previously mentioned, the drag force terms are not included within the five models to focus on the geometric effect.Furthermore, we perform the simulation with a fixed bed, although it has been noted that morphological changes due to tsunami-induced sediment transport can affect tsunami hazards [50][51][52][53].
Geosciences 2018, 8, 12 11 of 23 The buildings/houses and the breakwater in the bay are assumed to be intact; thus, the roughness coefficients and the parameters of the unit volume element in the porous medium, namely, A, H, L and H', are temporally constant values.As previously mentioned, the drag force terms are not included within the five models to focus on the geometric effect.Furthermore, we perform the simulation with a fixed bed, although it has been noted that morphological changes due to tsunamiinduced sediment transport can affect tsunami hazards [50][51][52][53].

Results and Discussions
Figure 7a,b show the maximum water level obtained by Landuse-n and the proposed PBM, respectively, where the green lines are the field survey result for the inundation area by Haraguchi and Iwamatsu [57].During the tsunami, the building destruction could lead to change the flow path area and inundation volume, while the generation of rubble might affect the drag force and the apparent viscosity in the multiphase flow.Furthermore, destruction of the breakwater and morphological change also occurred.However, despite the fact that the proposed PBM does not consider the destruction of buildings etc., the inundation area within the present model is relatively similar to the results of the field survey.These results imply that the inundation area was strongly limited by the steep slope behind the urban area and the large flow depth which is much greater than the building height rather than these uncertainties.Therefore, Figure 7 supports that notion that the numerical method and conditions including the tsunami source model are available to the tsunami inundation simulation in this area.

Results and Discussions
Figure 7a,b show the maximum water level obtained by Landuse-n and the proposed PBM, respectively, where the green lines are the field survey result for the inundation area by Haraguchi and Iwamatsu [57].During the tsunami, the building destruction could lead to change the flow path area and inundation volume, while the generation of rubble might affect the drag force and the apparent viscosity in the multiphase flow.Furthermore, destruction of the breakwater and morphological change also occurred.However, despite the fact that the proposed PBM does not consider the destruction of buildings etc., the inundation area within the present model is relatively similar to the results of the field survey.These results imply that the inundation area was strongly limited by the steep slope behind the urban area and the large flow depth which is much greater than the building height rather than these uncertainties.Therefore, Figure 7 supports that notion that the numerical method and conditions including the tsunami source model are available to the tsunami inundation simulation in this area.

Results and Discussions
Figure 7a,b show the maximum water level obtained by Landuse-n and the proposed PBM, respectively, where the green lines are the field survey result for the inundation area by Haraguchi and Iwamatsu [57].During the tsunami, the building destruction could lead to change the flow path area and inundation volume, while the generation of rubble might affect the drag force and the apparent viscosity in the multiphase flow.Furthermore, destruction of the breakwater and morphological change also occurred.However, despite the fact that the proposed PBM does not consider the destruction of buildings etc., the inundation area within the present model is relatively similar to the results of the field survey.These results imply that the inundation area was strongly limited by the steep slope behind the urban area and the large flow depth which is much greater than the building height rather than these uncertainties.Therefore, Figure 7 supports that notion that the numerical method and conditions including the tsunami source model are available to the tsunami inundation simulation in this area.Table 2 shows the inundation areas and the ratios to that using Landuse-n, where the target region is in the northwest range of the breakwater in the mouth of Onagawa Bay and is indicated by the dashed white lines as shown in Figure 7.The area using Landuse-n, which is one of the most practical methods, is applied as the reference value to compare between the models.The ratios for the field survey, Landuse-n, Landuse-n/Topography, the coastal-forest model, the flood model and PBM are 0.98, 1.00, 0.93, 0.95, 0.88 and 0.93, respectively.The numerical results except Landuse-n underestimate the field survey result within less than 10%, which probably is a consequence of the uncertainties mentioned above.For example, according to Adriano et al. [23], it is implied that their numerical results of inundation height and area might become closer to the field survey results by taking account of the destruction of the breakwater, which is one of uncertainties.Therefore, the differences between the ratios do not directly explain the accuracy of the modeling methods for the building effect.The improvements by considering uncertainties are the challenges in the future work.The inundation area acquired by Landuse-n is the largest among the five numerical models, and the model reveals the most dangerous tsunami hazards in terms of the inundation area.In addition, although the features of the porous medium are added to Landuse-n/Topography in the coastal-forest model, the flood model and PBM, the inundation area represented by the coastal-forest model is larger than that obtained using Landuse-n/Topography.The porosity can increase the water level and the potential energy.Meanwhile, the surface permeability limits the flow and decreases the inundation volume.For these reasons, the coastal-forest model including the porosity but excluding Table 2 shows the inundation areas and the ratios to that using Landuse-n, where the target region is in the northwest range of the breakwater in the mouth of Onagawa Bay and is indicated by the dashed white lines as shown in Figure 7.The area using Landuse-n, which is one of the most practical methods, is applied as the reference value to compare between the models.The ratios for the field survey, Landuse-n, Landuse-n/Topography, the coastal-forest model, the flood model and PBM are 0.98, 1.00, 0.93, 0.95, 0.88 and 0.93, respectively.The numerical results except Landuse-n underestimate the field survey result within less than 10%, which probably is a consequence of the uncertainties mentioned above.For example, according to Adriano et al. [23], it is implied that their numerical results of inundation height and area might become closer to the field survey results by taking account of the destruction of the breakwater, which is one of uncertainties.Therefore, the differences between the ratios do not directly explain the accuracy of the modeling methods for the building effect.The improvements by considering uncertainties are the challenges in the future work.The inundation area acquired by Landuse-n is the largest among the five numerical models, and the model reveals the most dangerous tsunami hazards in terms of the inundation area.In addition, although the features of the porous medium are added to Landuse-n/Topography in the coastal-forest model, the flood model and PBM, the inundation area represented by the coastal-forest model is larger than that obtained using Landuse-n/Topography.The porosity can increase the water level and the potential energy.Meanwhile, the surface permeability limits the flow and decreases the inundation volume.For these reasons, the coastal-forest model including the porosity but excluding the surface permeability might reproduce larger inundation areas than both Landuse-n/Topography and the other porous-type models.
Figure 8 shows the differences in the maximum water level between the porous-type models and Landuse-n/Topography.In urban areas, the flow depths calculated using the coastal-forest model are larger because of the effects of the porosity.The results using the PBM are similar to those in Landuse-n/Topography since the effects of the porosity and the surface permeability are moderately balanced.When using the flood model, since the porosity and surface permeability are vertically constant, the inundation volume within the urban area is greatly restricted, resulting in a in the inundation height and an increase in the water level over the sea due to the reflection of the tsunami from the urban area.
Geosciences 2018, 8, 12 14 of 23 the surface permeability might reproduce larger inundation areas than both Landuse-n/Topography and the other porous-type models.
Figure 8 shows the differences in the maximum water level between the porous-type models and Landuse-n/Topography.In urban areas, the flow depths calculated using the coastal-forest model are larger because of the effects of the porosity.The results using the PBM are similar to those in Landuse-n/Topography since the effects of the porosity and the surface permeability are moderately balanced.When using the flood model, since the porosity and surface permeability are vertically constant, the inundation volume within the urban area is greatly restricted, resulting in a decrease in the inundation height and an increase in the water level over the sea due to the reflection of the tsunami from the urban area.Figure 9 shows the distributions of the maximum flow velocity obtained by the five models; those simulated using the porous-type models (i.e., the coastal-forest model, the flood model and PBM) describe the Darcy velocity.In the porous-type models, the flow velocities between the buildings become large, thereby forming a reticulated high-flow velocity field.The differences in the maximum flow velocity between the porous-type models and Landuse-n/Topography are shown in Figure 10.The flow velocities calculated by the coastal-forest model in urban areas are smaller than those calculated using Landuse-n/Topography because of the increased flow depth with larger flow paths.Meanwhile, the flow velocities obtained using the flood model are relatively small due to the Figure 9 shows the distributions of the maximum flow velocity obtained by the five models; those simulated using the porous-type models (i.e., the coastal-forest model, the flood model and PBM) describe the Darcy velocity.In the porous-type models, the flow velocities between the buildings become large, thereby forming a reticulated high-flow velocity field.The differences in the maximum flow velocity between the porous-type models and Landuse-n/Topography are shown in Figure 10.The flow velocities calculated by the coastal-forest model in urban areas are smaller than those calculated using Landuse-n/Topography because of the increased flow depth with larger flow paths.Meanwhile, the flow velocities obtained using the flood model are relatively small due to the decreased inundation volume.Furthermore, the flow velocities from the PBM are locally larger than those from the other models.As shown in Table 2, the mean value of 100 % maximum velocity from Landuse-n is the largest among the five models.However, the values for the PBM in the 10% and 30% conditions are larger than those for the Landuse-n model as well as the other model, and it shows that the locally change in the velocity using the PBM is large.In the flood model, the inundation is restricted even if the flow depth is larger than building height as shown in Equations ( 22) and (24).For this reason, inundation volume is decreased, resulting in the decay of flow velocity across a wide area.Meanwhile, since both the flow path and the building height are considered the present model, the local increase in flow velocity is simulated.
Geosciences 2018, 8, 12 15 of 23 decreased inundation volume.Furthermore, the flow velocities from the PBM are locally larger than those from the other models.As shown in Table 2, the mean value of 100 % maximum velocity from Landuse-n is the largest among the five models.However, the values for the PBM in the 10% and 30% conditions are larger than those for the Landuse-n model as well as the other model, and it shows that the locally change in the velocity using the PBM is large.In the flood model, the inundation is restricted even if the flow depth is larger than building height as shown in Equations ( 22) and (24).For this reason, inundation volume is decreased, resulting in the decay of flow velocity across a wide area.Meanwhile, since both the flow path and the building height are considered in the present model, the local increase in flow velocity is simulated.Next, we examine the characteristics of the inundation processes calculated by the various models.Figure 11 shows the time series of the flow depth at the point near Marine Pal Onagawa shown in Figure 5b.The video analysis data were obtained by Koshimura et al. [45].The three-dimensional calculation results using a volume of fluid (VOF) method with a 1-m resolution obtained by Arikawa and Tomita [5,6] are labeled "VOF (1 m, TT)" and "VOF (1 m, FS8)", respectively.The tsunami sources are applied the models proposed by Takagawa & Tomita [58] and Satake et al. [56], which are hereafter called "TT" and "FS8" (Fujii and Satake model version 8.0), respectively.The "SPH (2 m, FS8)" calculated results were obtained by Suwa et al. [9] using SPH method with a 2-m resolution from the FS8 tsunami source model.The "Landuse-n/Topography (5 m, COGJ)" and "Equivalent-n (5 m, COGJ)" calculated results were acquired by Adriano et al. [23] using TUNAMI-N2 by employing the Landuse-n/Topography and Equivalent-n models for 5-m resolution, respectively, with the tsunami source model proposed by the Cabinet Office, Government of Japan [59], which is hereafter called the "COGJ" model.Comparing the results simulated by this study with those from three-dimensional simulations (the VOF and SPH) by using the FS8 tsunami source model, the results from the PBM agree well with the 3-D results obtained by Arikawa and Tomita [5,6].It implies that the two-dimensional simulations demonstrate good accuracies for the inundation heights around buildings near the coastline.On the other hand, the maximum water level calculated by Suwa et al. [9] is larger than in the other models (i.e., the VOF, Landuse-n and PBM).As noted in their study, the reflections along boundaries, which connects the SPH method with the two-dimensional calculations, and the artificial viscosity term might result in an excessive increase in the water level.
However, the results from the PBM are different from the observation in term with the phase of the tsunami.It appears to be result of the applied tsunami source model.The accuracies of the numerical models shown in Figure 11 cannot be compared unconditionally with one another because the applied tsunami source models are different.The phases of the COGJ and TT tsunami source models are relatively similar to the video analysis results, whereas those of the FS8 model are relatively early.Meanwhile, the collapsed breakwater and the erosion of the seabed affect the phase of the tsunami in this area, and these damages could have accelerated the arrival time of the tsunami.For the tsunami in Onagawa town, the COGJ and TT models are more adequate than the FS8 model in terms of the phase of the tsunami, which also depends on these damage processes.
Comparing the two results obtained by this study or those reported by Adriano et al. [23], the differences in the modeled building effects do not significantly affect the maximum water level near the coastline.In fact, as far as we can compare the five sets of calculated results in this study (two from Landuse-n, Landuse-n/Topography, the coastal-forest model, the flood model and PBM are shown in Figure 11), the differences in the maximum water level at the point near Marine Pal Onagawa are less than 0.5 m.These results imply that differences in the modeling of buildings has a limited effect on the maximum water level near the coastline, which may be related to the topography, i.e., the aforementioned steep slope behind the urban area, and the vicinity of the observation point to the coastline.Namely, since the tsunami is blocked by the steep slope, which results in an increase in the water level, the flow at the maximum water level is relatively weak.Thus, an increase in the water level due to run-up with strong flow around a building, which results from the geometric effect, is unlikely at this location during the period with the maximum water level.In addition, in the case of a similar tsunami scenario wherein the city area is completely submerged, the water level around the coastline during the stagnation period of the flow strongly depends on the exterior scale of the tsunami, and thus, the building conditions within the city area hardly influence the maximum water level at the observation point.However, the differences in the modeled building effects might become substantial when, in contrast to a region with a steep slope, a progressive tsunami over a plain area is simulated, where the flow speed of a tsunami could become large when the water level is high.In fact, comparing the results using the PBM with those from Landuse-n in this study for the period of run-up before the water level reaches its maximum value (35 min-45 min), the water level in the PBM is higher than that in Landuse-n, which does not consider the geometric Comparing the results simulated by this study with those from three-dimensional simulations (the VOF and SPH) by using the FS8 tsunami source model, the results from the PBM agree well with the 3-D results obtained by Arikawa and Tomita [5,6].It implies that the two-dimensional simulations demonstrate good accuracies for the inundation heights around buildings near the coastline.On the other hand, the maximum water level calculated by Suwa et al. [9] is larger than in the other models (i.e., the VOF, Landuse-n and PBM).As noted in their study, the reflections along boundaries, which connects the SPH method with the two-dimensional calculations, and the artificial viscosity term might result in an excessive increase in the water level.
However, the results from the PBM are different from the observation in term with the phase of the tsunami.It appears to be result of the applied tsunami source model.The accuracies of the numerical models shown in Figure 11 cannot be compared unconditionally with one another because the applied tsunami source models are different.The phases of the COGJ and TT tsunami source models are relatively similar to the video analysis results, whereas those of the FS8 model are relatively early.Meanwhile, the collapsed breakwater and the erosion of the seabed affect the phase of the tsunami in this area, and these damages could have accelerated the arrival time of the tsunami.For the tsunami in Onagawa town, the COGJ and TT models are more adequate than the FS8 model in terms of the phase of the tsunami, which also depends on these damage processes.
Comparing the two results obtained by this study or those reported by Adriano et al. [23], the differences in the modeled building effects do not significantly affect the maximum water level near the coastline.In fact, as far as we can compare the five sets of calculated results in this study (two from Landuse-n, Landuse-n/Topography, the coastal-forest model, the flood model and PBM are shown in Figure 11), the differences in the maximum water level at the point near Marine Pal Onagawa are less than 0.5 m.These results imply that differences in the modeling of buildings has a limited effect on the maximum water level near the coastline, which may be related to the topography, i.e., the aforementioned steep slope behind the urban area, and the vicinity of the observation point to the coastline.Namely, since the tsunami is blocked by the steep slope, which results in an increase in the water level, the flow at the maximum water level is relatively weak.Thus, an increase in the water level due to run-up with strong flow around a building, which results from the geometric effect, is unlikely at this location during the period with the maximum water level.In addition, in the case of a similar tsunami scenario wherein the city area is completely submerged, the water level around the coastline during the stagnation period of the flow strongly depends on the exterior scale of the tsunami, and thus, the building conditions within the city area hardly influence the maximum water level at the observation point.However, the differences in the modeled building effects might become substantial when, in contrast to a region with a steep slope, a progressive tsunami over a plain area is simulated, where the flow speed of a tsunami could become large when the water level is high.In fact, comparing the results using the PBM with those from Landuse-n in this study for the period of run-up before the water level reaches its maximum value (35 min-45 min), the water level in the PBM is higher than that in Landuse-n, which does not consider the geometric effect; therefore, this result is reasonable.Meanwhile, comparing the results from Landuse-n/Topography with those from Equivalent-n acquired by Adriano et al. [23] for the corresponding period, the mechanical effects modeled by Equivalent-n serve as geometric effects similar to found in Landuse-n/Topography.However, Muhari et al. [19] clarified that the flow depth and velocity calculated by Equivalent-n is locally smaller than that calculated by Landuse-n/Topography because these localities are gradually strengthened inside an urban area where buildings are lined up.Namely, the locality of tsunami hazards could become stronger (like an internal area of an urban setting).
Focusing on the inside area of urban setting, the characteristics of locally tsunami hazards from each model are examined.Figure 12a,b show snapshots of the water level calculated by the Landuse-n model and the PBM, respectively.The gradient of the water surface simulated using the PBM in the direction normal to the shoreline is larger than that obtained by Landuse-n because of the geometric effect on the interception of the tsunami.This indicates that larger pressure gradients will lead to an amplification of the flow velocity.As a result, a high-speed flow between the buildings will occur, especially to the west of Marine Pal Onagawa, where a large flow velocity is generated along the straight road, as shown in Figure 13.
Figure 14a shows the time series of the inundation height and the absolute value of the flow velocity at Point A shown in Figure 5b.Due to the geometric effects of the buildings in the Landuse-n/Topography model, the coastal-forest model, the flood model and PBM, the arrival times of the tsunami in those models are delayed coincident with a gradual rise in the water level (approximately 2300 s), after which the water levels are elevated suddenly (after 2400 s), leading to high flow velocities, especially in the flood model and the PBM (Figure 14b).However, the run-up speed calculated by the coastal-forest model is larger than that calculated by the Landuse-n/Topography model because the increase in the water level due to the porosity makes it easier for the tsunami to spread out over the land.
In the flood model, the surface permeability is expressed as the mean value of two porosities, as shown in Equation (24).Thus, in the case when D < H, the surface permeability in the flood model tends to be larger than in the PBM over the straight road, where the porosity spatially changes discontinuously.A large surface permeability value means that a tsunami over the road will tend to disperse.As a consequence, the flood model also exhibits smaller flow velocities than the PBM over the road.Namely, since the PBM considers the geometric effects of buildings in detail, the PBM will restrict the progress of the tsunami according the building distribution and reproduce an adequate flow field concentrating along the straight road.
According to the above results, the calculated tsunami hazards show that the results are reasonable and commensurate with the degree of the geometric effects introduced in the models.Therefore, consideration of the geometric effect demands evaluating the tsunami hazards with respect to the locality.The porosity is important for avoiding an underestimation of not only the inundation height but also the inundation area (based on the comparison of the Landuse-n model to the coastal-forest model, which adds porosity to the Landuse-n model), and the surface permeability plays an important role in representing the local high flow velocity.We must investigate the validity of the modeling and the applicability to fields including various uncertainties through comparisons to an ideal experiment or theory in the future.

Conclusions and Recommendations
To evaluate tsunami hazards with strong locality in urban areas, NSWW equations based on the porous body model were derived, and a two-dimensional horizontal tsunami inundation analysis model was developed.The proposed PBM introduces an adequate kinematic boundary condition that includes both porosity and surface permeability of the porous body model within the derivation process.Meanwhile, geometric effects such as the flow anisotropy and building heights of groups of buildings, which present a challenge in conventional models, were incorporated into the numerical model more accurately.The proposed numerical model was applied to a simulation of the tsunami inundation in the town of Onagawa, Miyagi Prefecture, Japan, during the 2011 Great East Japan Earthquake, and the characteristics of the tsunami hazards calculated by various numerical models were investigated.Our results suggest the following conclusions.Earthquake, and the characteristics of the tsunami hazards calculated by various numerical models were investigated.Our results suggest the following conclusions.Earthquake, and the characteristics of the tsunami hazards calculated by various numerical models were investigated.Our results suggest the following conclusions.Earthquake, and the characteristics of the tsunami hazards calculated by various numerical models were investigated.Our results suggest the following conclusions.• The proposed PBM exhibited as good accuracy for the inundation heights around building near the coastline as with conventional three-dimensional simulation with high resolution.

•
In the case where the inundation area is restricted by the topography with a steep slope behind the urban area and the flow depth is much greater than the building height, the differences in the modeled building effects and computational resolution do not significantly affect the maximum water level near the coastline.Therefore, two-dimensional simulations with a practical resolution demonstrate good accuracy while estimating the inundation height near the coastline, although we must examine the applicability of these models to regions where the locality can become stronger in future work.

•
A model that includes porosity can yield an increase in the water level, which makes it easier for the tsunami to spread out over the land, thereby decreasing the flow velocity.

•
A model that includes the surface permeability restricts the progress of the tsunami according to the distribution of the buildings and reproduces a flow field concentrated along the straight road.

•
By properly incorporating the porosity and the surface permeability into the theory and the numerical model, the model can adequately reproduce high flow velocities due to the increase in the gradient of the water surface and the concentration of the flow during the inundation process within the urban area.In addition, the water level at the time will increase due to the flow and geometric effects.Therefore, a numerical model that does not consider geometric effects could underestimate the local hydrodynamic force.
The proposed PBM showed its capability to reproduce the tsunami's essential behaviors in urban areas.It can be expected that the present model can become a useful tool to accurately evaluate the tsunami risks in urban areas.In future works, we will investigate the accuracy of the proposed model using the drag force term and its applicability in the field.

Figure 1 .
Figure 1.Unit volume element in the porous body model, where the black and white are solid and fluid, respectively.

Figure 1 .
Figure 1.Unit volume element in the porous body model, where the black and white are solid and fluid, respectively.

Figure 2 .
Figure 2. Modeling process for buildings using the porous body model.(a) Buildings on computational mesh; (b) Buildings on cell mesh; (c) Occupancy area An on cell mesh and length Ln' at the edge of cell mesh; and (d) porosity of the water column on cell mesh and horizontal surface permeability at the edge of the cell mesh.

Figure 2 .
Figure 2. Modeling process for buildings using the porous body model.(a) Buildings on computational mesh; (b) Buildings on cell mesh; (c) Occupancy area A n on cell mesh and length L n at the edge of cell mesh; and (d) porosity of the water column on cell mesh and horizontal surface permeability at the edge of the cell mesh.

Figure 3 .
Figure 3. Accumulation of crustal deformation in the tsunami source model proposed by Satake et al.[24], where the circle indicates the location of Onagawa, Miyagi Prefecture, Japan, and the star is the epicenter of the earthquake.The units of distance are km.

Figure 3 .
Figure 3. Accumulation of crustal deformation in the tsunami source model proposed by Satake et al.[24], where the circle indicates the location of Onagawa, Miyagi Prefecture, Japan, and the star is the epicenter of the earthquake.The units of distance are km.

Figure 5 .
Figure 5. Topography in Onagawa, Miyagi Prefecture, Japan (Region 6: Δx6 = 5 m).(a) The blue line is the inundation area measured by Haraguchi and Iwamatsu [57], and the black colors are the locations where the building height data are incorporated into the topographic data, i.e., it indicates A/Ac > 0.95 in Equation (18).(b) The ground elevation around Marine Pal Onagawa, which is the region in the red box shown in Figure 4a.

Figure 6 .
Figure 6.The occupancy area ratio of buildings around Marine Pal Onagawa in Figure 4b.The black colors indicate the locations where building height data are incorporated into the topographic data, i.e., it indicates A/Ac > 0.95 in Equation (18).

Figure 5 .
Figure 5. Topography in Onagawa, Miyagi Prefecture, Japan (Region 6: ∆x 6 = 5 m).(a) The blue line is the inundation area measured by Haraguchi and Iwamatsu [57], and the black colors are the locations where the building height data are incorporated into the topographic data, i.e., it indicates A/A c > 0.95 in Equation (18).(b) The ground elevation around Marine Pal Onagawa, which is the region in the red box shown in Figure 5a.

Figure 5 .
Figure 5. Topography in Onagawa, Miyagi Prefecture, Japan (Region 6: Δx6 = 5 m).(a) The blue line is the inundation area measured by Haraguchi and Iwamatsu [57], and the black colors are the locations where the building height data are incorporated into the topographic data, i.e., it indicates A/Ac > 0.95 in Equation (18).(b) The ground elevation around Marine Pal Onagawa, which is the region in the red box shown in Figure 4a.

Figure 6 .
Figure 6.The occupancy area ratio of buildings around Marine Pal Onagawa in Figure 4b.The black colors indicate the locations where building height data are incorporated into the topographic data, i.e., it indicates A/Ac > 0.95 in Equation (18).

Figure 6 .
Figure 6.The occupancy area ratio of buildings around Marine Pal Onagawa in Figure 5b.The black colors indicate the locations where building height data are incorporated into the topographic data, i.e., it indicates A/A c > 0.95 in Equation (18).

Figure 7 .
Figure 7. Maximum water level obtained by the proposed PBM, where the green line is the field survey result for the inundation area [57].

Figure 7 .
Figure 7. Maximum water level obtained by the proposed PBM, where the green line is the field survey result for the inundation area [57].
(a) Difference between Coastal-forest and Landuse-n/Topography (b) Difference between Flood and Landuse-n/Topography (c) Difference between PBM and Landuse-n/Topography

Figure 8 .
Figure 8. Difference in the maximum water level between the porous-type models (Coastal-forest: coastal-forest model, Flood: flood model, PBM: proposed model based on the porous body model) and Landuse-n/Topography, where the green lines are the field survey result for the inundation area [57].

Figure 8 .
Figure 8. Difference in the maximum water level between the porous-type models (Coastal-forest: coastal-forest model, Flood: flood model, PBM: proposed model based on the porous body model) and Landuse-n/Topography, where the green lines are the field survey result for the inundation area [57].

Figure 9 .
Figure 9. Maximum flow velocities obtained using the (a) Landuse-n: roughness map model, (b) Landuse-n/Topography: roughness map and topographic model, (c) Coastal-forest: coastal forest model, (d) Flood: flood model, and (e) PBM: proposed model based on the porous body model, where the porous-type models (i.e., the coastal-forest model, the flood model and PBM) describe the Darcy velocity and the pink lines are the field survey result for the inundation area [57].

Figure 9 .
Figure 9. Maximum flow velocities obtained using the (a) Landuse-n: roughness map model, (b) Landuse-n/Topography: roughness map and topographic model, (c) Coastal-forest: coastal forest model, (d) Flood: flood model, and (e) PBM: proposed model based on the porous body model, where the porous-type models (i.e., the coastal-forest model, the flood model and PBM) describe the Darcy velocity and the pink lines are the field survey result for the inundation area [57].

Figure 10 .
Figure 10.Difference in the maximum flow velocity between the porous-type models (i.e., the coastal-forest model, the flood model and PBM) and Landuse-n/Topography, where the green lines are the field survey result for the inundation area [57].

Figure 11 .
Figure 11.Time series of the flow depth at the point near Marine Pal Onagawa shown in Figure 5b.

Figure 11 .
Figure 11.Time series of the flow depth at the point near Marine Pal Onagawa shown in Figure 5b.

Figure 12 .Figure 13 .Figure 14 .
Figure 12.Water level at 2400 s after the earthquake: (a) Landuse-n (roughness map model); (b) the proposed PBM based on the porous body model.

Figure 12 .
Figure 12.Water level at 2400 s after the earthquake: (a) Landuse-n (roughness map model); (b) the proposed PBM based on the porous body model.

Figure 12 .Figure 13 .Figure 14 .
Figure 12.Water level at 2400 s after the earthquake: (a) Landuse-n (roughness map model); (b) the proposed PBM based on the porous body model.

Figure 13 .
Figure 13.Flow velocity at 2400 s after the earthquake: (a) Landuse-n (roughness map model); (b) the proposed PBM based on the porous body model, where the pore velocity is described.

Figure 12 .Figure 13 .Figure 14 .
Figure 12.Water level at 2400 s after the earthquake: (a) Landuse-n (roughness map model); (b) the proposed PBM based on the porous body model.

Figure 14 .
Figure 14.Time series of the inundation height and flow velocity at Point A in Figure 5b: (a) inundation height; (b) flow velocity, where the pore velocity is described (Landuse-n: roughness map model; Landuse-n/Topography: roughness map and topographic model; Coastal-forest: coastal forest model; Flood: flood model; PBM: proposed model based on porous body model).

Table 1 .
Features of existing calculation methods for urban areas.

Table 2 .
Inundation area and maximum flow velocities in Onagawa town as shown in Figure.5a(field survey [57]; Landuse-n: roughness map model; Landuse-n/Topography: roughness map and topographic model; Coastal-forest: coastal forest model; Flood: flood model; PBM: proposed model based on porous body model).

Table 2 .
[57]dation area and maximum flow velocities in Onagawa town as shown in Figure5a(field survey[57]; Landuse-n: roughness map model; Landuse-n/Topography: roughness map and topographic model; Coastal-forest: coastal forest model; Flood: flood model; PBM: proposed model based on porous body model).