Transitional Behavior of a Flow Regime in Shoaling Tsunami Boundary Layers

The transitional flow regime of the bottom boundary layer under hypothetical shoaling tsunamis is investigated in the entire region from the tsunami source to the shallow sea area. In order to calculate the shoaling process of a tsunami, an analytical method based on Green’s law and the linear long wave theory are employed, and flow regime criteria for the wave boundary layer proposed by one of the authors are applied. It is found that the bottom boundary layer in a tsunami source area is located in the laminar regime. Subsequently, transition occurs to the smooth turbulence during the shoaling process, with a transition from the smooth to the rough turbulent region in the shallow area. For precise evaluation of bottom friction acting on the sea bed and the resulting energy dissipation beneath the tsunami, it is highly necessary to include such transitional behavior in sea bottom boundary layers.


Introduction
In general, the nature of the flow pattern can be classified into laminar and turbulent regimes depending on the relevant Reynolds number. In the case of a flow motion with a rigid boundary, the latter is further classified into smooth and rough turbulent flow depending on the roughness condition on the wall. Since the frictional resistance on the wall surface and resulting energy dissipation are significantly dependent on the mode of flow motion, identification and understanding of the flow regime is of practical importance for various flow phenomena.
In many practical applications in the field of hydraulic engineering, a Manning's roughness coefficient, which is valid only in a fully rough turbulent regime under a steady uniform condition, has been widely utilized in various flow situations. Even under tsunami waves, a Manning's coefficient or equivalent expression in terms of Darcy-Weisbach's f has commonly been utilized in order to evaluate bottom shear stress in tsunami numerical simulations for various hazard mitigation studies, simply assuming that the wave period of the tsunami is sufficiently long to satisfy the quasi-steady condition [1][2][3]. However, there have been quite few investigations on tsunami boundary layer development on the sea bottom, especially in regard to transitional behavior of flow regimes in tsunami boundary layers.
Tanaka et al. [4] carried out a detailed experimental study to explore the properties of a depth-limited wave boundary layer. A criterion has been proposed for the inception of the depth-limited state ( Figure 1a) on a rough bottom. According to the first reported field observation during the 2010 Chilean earthquake tsunami in U.S. by Lacy et al. [5], the depth-limited condition is not satisfied even at h = 10 m (h: the water depth). In addition, the acceleration effect is significant to the near-bed momentum, and bottom shear stress is not in phase with the velocity. This unsteady behavior cannot be simulated using a conventional tsunami simulation model with a friction term expressed in terms of quadratic friction law. William and Fuhrman [6] also concluded that over the range of conditions they considered, the tsunami boundary layer thickness grows to O (1 m), as demonstrated in Figure  1b and typically is not limited by the water depth as a steady flow in an open channel (Figure 1a). According to Sumer and Fuhrman [7], investigations for the 2004 Indian Ocean tsunami and the 2011 Tohoku tsunami suggest that the tsunami-induced boundary layers at the depths of respectively h = 14 m and h = 20 m would only span a fraction of the total water depth. In an experimental study of tsunami-induced scour around a monopole foundation carried out by Larsen et al. [8], it is shown that under the condition of an extremely long period, a current-like boundary layer is observed. In the recent investigations of authors, it was concluded that the bottom boundary layer behavior beneath the tsunami is similar to that induced by wind waves, rather than the steady flow boundary layer in the most computational domain from the tsunami source area to the shallow water region (Tinh and Tanaka [9]). However, Tinh and Tanaka [9] simply assumed in the conventional manner the tsunami boundary layer is under a rough turbulent condition, without detailed discussion of the inherent flow regime transition. The present study is a supplementary work to Tinh and Tanaka [9] in the sense that the transitional behavior of the flow regime in a bottom boundary layer beneath hypothetical shoaling tsunamis is focused on. The result of the present study enables us to predict precise bottom shear stress under a tsunami by identifying and understanding the flow regime for a given input tsunami condition.

Flow Regime in a Wave Boundary Layer
Jonsson [10,11] was the first to propose a flow regime diagram under wave motions as well as a friction factor diagram covering both laminar and turbulent flows. Later, Kamphuis [12] obtained similar diagrams for wave friction factors and flow regimes based on a laboratory experiment of direct shear stress measurement in an oscillating water tunnel. Myrhaug [13] proposed a theory for a wave friction coefficient and the phase lead of bottom shear stress for rough, smooth and transitional smooth-to-rough turbulent flow. Tanaka and Sana [14] applied the k-ε turbulence model to investigate the transitional regime in a wave boundary layer on a smooth bottom. An abrupt change in friction factor and phase difference is successfully predicted by the model, which is well in agreement with the available experimental data. For combined motion of the wave and current, In the recent investigations of authors, it was concluded that the bottom boundary layer behavior beneath the tsunami is similar to that induced by wind waves, rather than the steady flow boundary layer in the most computational domain from the tsunami source area to the shallow water region (Tinh and Tanaka [9]). However, Tinh and Tanaka [9] simply assumed in the conventional manner the tsunami boundary layer is under a rough turbulent condition, without detailed discussion of the inherent flow regime transition. The present study is a supplementary work to Tinh and Tanaka [9] in the sense that the transitional behavior of the flow regime in a bottom boundary layer beneath hypothetical shoaling tsunamis is focused on. The result of the present study enables us to predict precise bottom shear stress under a tsunami by identifying and understanding the flow regime for a given input tsunami condition.

Flow Regime in a Wave Boundary Layer
Jonsson [10,11] was the first to propose a flow regime diagram under wave motions as well as a friction factor diagram covering both laminar and turbulent flows. Later, Kamphuis [12] obtained similar diagrams for wave friction factors and flow regimes based on a laboratory experiment of direct shear stress measurement in an oscillating water tunnel. Myrhaug [13] proposed a theory for a wave friction coefficient and the phase lead of bottom shear stress for rough, smooth and transitional smooth-to-rough turbulent flow. Tanaka and Sana [14] applied the k-ε turbulence model to investigate the transitional regime in a wave boundary layer on a smooth bottom. An abrupt change in friction factor and phase difference is successfully predicted by the model, which is well in agreement with the available experimental data. For combined motion of the wave and current, similar investigations on flow regime transition have been carried out by Tanaka and Shuto [15], Christoffersen and Jonsson [16], Myrhaug and Slaattelid [17], and Tanaka and Thu [18].
Depending on the methodology for identifying and defining flow regime transition employed by each researcher, flow regime criteria are slightly different from one another. In this study, the flow regime criteria proposed by Tanaka and Thu [18] will be used, assuming that the steady current component is zero in their equations originally obtained for wave-current co-existent motion. The transition from laminar to smooth turbulent flow is given by the following equations.
Lower limit: R e = 2.5 × 10 5 Upper limit: R e = 6.0 × 10 5 Here, Reynolds number R e is defined by where U m : the amplitude of the wave-induced velocity outside the boundary layer, a m : the amplitude of water particle motion outside the boundary layer (= U m /σ), σ: the angular frequency (= 2π/T, T: the wave period), and ν: the molecular viscosity of the fluid. Williams and Fuhrman [6] carried out numerical analysis of tsunami-like wave boundary layers using a transitional variant of the k-ω turbulence model. According to their numerical result, the critical Reynold numbers similar to Equations (1) and (2) are valid for the transition to turbulence in a tsunami boundary layer.
With an increase in wave Reynolds number, another transition occurs in the wave boundary layer from a smooth to a rough turbulent regime. To recognize this type of transition under combined waves and current, both Myrhaug and Slaattelid [17] and Tanaka and Thu [18] applied Equations (4) and (5) proposed by Schlichting [19] for a steady current expressed in terms of roughness Reynolds number, k + s . Lower limit: Upper limit: where u * m : the maximum of friction velocity, and k s : the equivalent roughness length. The left side of Equations (4) and (5) represents the ratio of equivalent roughness to the viscous sublayer thickness. For practical convenience of engineering applications, Tanaka and Thu [18] further modified Equations (4) and (5) to express them in terms of Re and a m /k s as follows: Lower limit: Upper limit: Equations (4) and (5) or Equations (6) and (7) will be applied for recognition of transitional behavior in the present investigation on the tsunami boundary layer.

Tsunami Shoaling Process
Four cases of input tsunami characteristics in the source area are examined as listed in Table 1, having different wave periods T and seabed sediment diameters d. The equivalent roughness is obtained from the diameter of bed material using k s = 2d. In the case of sea bed with bed forms, vegetation or other forms of benthic aquatic life, the bed roughness will remarkably increase as compared with simple grain roughness considered in this study. The bathymetry of the sea area is assumed to have a 1/50 slope as depicted in Figure 2a, in which the x-coordinate is defined with the origin at the offshore boundary. Following Williams and Fuhrman [6] and Tinh and Tanaka [9], the wave shoaling process from the source to the shallow sea region is obtained by applying Green's law.
where : the wave height at an arbitrary point with the water depth ℎ, and the subscript "0" denotes quantity in the tsunami source area. Furthermore, the following linear long wave theory is applied Following Williams and Fuhrman [6] and Tinh and Tanaka [9], the wave shoaling process from the source to the shallow sea region is obtained by applying Green's law.
where H: the wave height at an arbitrary point with the water depth h, and the subscript "0" denotes quantity in the tsunami source area. Furthermore, the following linear long wave theory is applied to calculate the amplitude of tsunami-induced free-stream velocity U m .
where g: the gravitational acceleration.

Calculation of Friction Factor and Friction Velocity
In order to apply Equations (4) and (5), the maximum shear velocity under the tsunami is required to distinguish the flow regime from smooth and rough turbulent flow. In this study, the full-range equation proposed by Tanaka and Thu [18], originally proposed for wave-current combined motion, will be used to calculate wave friction coefficient under shoaling tsunami, f w , since it can provide the smoothly interpolated value for each flow regime even in the transitional regime. The definition of the friction factor and the relationship with shear velocity is where τ 0m is the maximum bottom shear stress under wave motion, and ρ is the fluid density. Calculation is made at selected x-coordinates depicted in Table 2.

Calculation of Boundary Layer Thickness
For computing wave boundary layer thickness under a shoaling tsunami, the full-range equation originally proposed by Sana and Tanaka [20,21] and recently modified by the authors [22] will be applied. The definition of wave boundary thickness proposed by Jensen et al. [23] is used in this study, which is the height of overshooting at σt = 0 under the wave crest as defined in Figure 1b. There are alternative definitions for wave boundary layer thickness, such as Jonsson [10] based on the minimum distance from the bottom to an elevation where the velocity equals the amplitude of the free-stream velocity, U m , and Sleath [24] and Yuan and Madsen [25] in terms of the distance where the defect velocity amplitude is 1% or 5% of the free-stream amplitude. Among these definitions, that of Jensen et al. [23] has successfully been applied for demarcating the friction factor under the tsunami [22].
Both of the full-range equations for wave friction factor and wave boundary layer thickness utilized in the present investigation are summarized in Table 3, in which Equation (12) is employed from Equation (6) for smooth interpolation in the transitional region between smooth and rough turbulent regime.  Table 3. Full-range equations for the wave friction coefficient and wave boundary layer thickness.

Transitional Flow Regime in Terms of Boundary Layer Thickness
The wave height and maximum velocity under the shoaling tsunami are computed using Equations (8) and (9), respectively, as shown Table 2, and are plotted in Figure 2b,c. Substituting the a m /k s and R e values in Table 2 into the full-range equation for wave boundary layer thickness shown in Table 3(b), Figure 3a,b are obtained for transitional variation of the boundary layer thickness under the shoaling tsunami. It is observed in Figure 3a (d = 0.3 mm) that transitional behavior occurs from laminar to smooth turbulence, and finally slight deviation from smooth turbulence to the smooth-rough transitional regime at h = 10 m is observed. Meanwhile, as shown in Figure 3b, with finer bed material (d = 0.1 mm), transition from smooth turbulence cannot be seen even at h = 10 m. The ratio of the boundary thickness to the water depth, /ℎ, in the entire region of tsunami shoaling is shown in Figure 2d. As seen in the figure, the boundary layer is very thin as compared with the water depth in all cases as explained in Figure 1b, and it can be confirmed that the steady flow friction coefficient cannot be applied, although almost all of the past studies on tsunami numerical simulation assumed so. As depicted in Figure 1, because of the steep velocity gradient in the wave boundary layer, the bottom shear stress from the wave friction factor ( ( ) ) is usually greater than that from steady flow theory ( ( ) ).     The ratio of the boundary thickness to the water depth, δ/h, in the entire region of tsunami shoaling is shown in Figure 2d. As seen in the figure, the boundary layer is very thin as compared with the water depth in all cases as explained in Figure 1b, and it can be confirmed that the steady flow friction coefficient cannot be applied, although almost all of the past studies on tsunami numerical simulation assumed so. As depicted in Figure 1, because of the steep velocity gradient in the wave boundary layer, the bottom shear stress from the wave friction factor τ 0(2) is usually greater than that from steady flow theory (τ 0(1) ).

Transitional Flow Regime in Terms of the R e vs. a m /k s Relationship
In Figure 4, the transitional flow regime under tsunami shoaling is plotted in a diagram in terms of R e and a m /k s (Tanaka and Thu [18]). The demarcation lines in this diagram are obtained from Equations (1), (2), (6) and (7). The relationship between R e and a m /k s is expresses as where δ S is the Stokes layer thickness defined as [19] δ S = 2ν σ  It is confirmed that in Figure 4, after departing the source area, transition occurs from laminar to transitional flow, from transitional to smooth turbulent flow, and finally from smooth turbulent to transitional area. In the condition considered in the preset study, however, a fully rough turbulent regime cannot be reached even at the water depth of 10 m.  Therefore, in Figure 4, the plotted data of each shoaling tsunami are located on the straight line with a slope of 1/2 in log-log plots.
The bottom boundary layer in the tsunami source area is located in laminar regime in Figure 4, as previously reported by Williams and Fuhrman [6]. According to the linear long wave theory, the maximum tsunami-induced velocity in the source area from Equation (9) is in the order of O (1 cm/s) as seen in Table 2, resulting in an extremely small Reynolds number-much smaller than the critical values given by Equation (1).
In Figure 4, field data obtained by Lacy et al. during the 2020 Chilean Tsunami are plotted, and are located slightly below Equation (7).
Under wind-generated wave boundary layers, there have been numerous laboratory experiments under full-scale sea wave conditions using an oscillating flow tunnel, such as those by Jonsson and Carlsen [26], Sleath [24], Jensen et al. [23], Tanaka et al. [4], and van der A et al. [27]. These experimental conditions for wind-generated waves are plotted in Figure 4 for comparison with tsunami boundary layers investigated in this study. It is confirmed that the conditions of wind-generated wave boundary layer studies with rough bottom are confined to the range in the left bottom corner of Figure 4 with a smaller value of R e and a m /k s , whereas the smooth bottom experiments that are plotted on the axis of a m /k s are infinity.
It is confirmed that in Figure 4, after departing the source area, transition occurs from laminar to transitional flow, from transitional to smooth turbulent flow, and finally from smooth turbulent to transitional area. In the condition considered in the preset study, however, a fully rough turbulent regime cannot be reached even at the water depth of 10 m.

Laminar-to-Turbulent Transition
Substituting Equation (9) into Equation (3), Reynolds number R e(S) at the tsunami source area is obtained.
Further modification of Equation (25) can be made by substituting constant values such as g, π and ν (= 1.0 × 10 −6 m 2 /s). R e(S) = 3.9 × 10 5 × H 2 0 T h 0 (SI units) (26) In order to discuss the Reynolds number at the tsunami source area, h 0 = 4000 m and H 0 = 1 m are substituted into Equation (26), and only the period T of tsunami will be treated as an independent variable. The relationship between Reynolds number and wave period is illustrated in Figure 5. When the period exceeds T = 42.7 min, transition occurs from laminar to turbulent flow according to Equation (1). In order to cross the upper limit of transition given by Equation (2), a much longer wave period is required at the tsunami source area.
Meanwhile, Reynolds number R e at an arbitrary location can be correlated with R e(S) defined by Equation (26) as Equation (27) indicates that reduction of water depth causes increase in the wave Reynolds number in the shallow area. The calculation result of the Reynolds number by Equation (27)  It is noteworthy that the Reynolds number relevant to a steady current, R, is defined using the water depth as a characteristic length.
Similar to Equation (27), R at an arbitrary place can be correlated with that at the tsunami source, R (s) , as where R (s) is defined as According to Equation (29), R decreases with decreasing water depth during the tsunami propagation, which is opposite to the behavior observed for Equation (27).

Meanwhile, Reynolds number
at an arbitrary location can be correlated with ( ) defined by Equation (26)  Equation (27) indicates that reduction of water depth causes increase in the wave Reynolds number in the shallow area. The calculation result of the Reynolds number by Equation (27) is shown in Figure 2e. Since the Reynolds number does not include the bottom roughness, the lines in Figure 2e are drawn separately for cases with different wave periods. It is observed that the Reynolds number for the T = 15 min tsunami (Case 1 and Case 3) exceeds the critical values given in Equations (1) and (2), at x = 101 km and x = 145 km, respectively, whereas for the T = 30 min tsunami (Case 2 and Case 4), this occurs at x = 42 km and x = 112 km, respectively.
It is noteworthy that the Reynolds number relevant to a steady current, , is defined using the water depth as a characteristic length.
Similar to Equation (27), at an arbitrary place can be correlated with that at the tsunami source, ( ) , as where ( ) is defined as According to Equation (29), decreases with decreasing water depth during the tsunami propagation, which is opposite to the behavior observed for Equation (27).
The steady current Reynolds number is obtained in Table 2 for the shoaling tsunami. According to Chow [28], the critical Reynolds number for transition to turbulence in a steady open channel flow is The steady current Reynolds number is obtained in Table 2 for the shoaling tsunami. According to Chow [28], the critical Reynolds number for transition to turbulence in a steady open channel flow is The R value in Table 2 is in the order of O (10 7 ), which remarkably exceeds the above critical value for a steady flow, as illustrated in Figure 2f. Hence, based on the steady current Reynolds number, it might be presumed that the boundary layer beneath the tsunami is fully turbulent. However, the steady current Reynolds number is not suitable under the tsunami for recognizing the flow regime, since the development of the sea bottom boundary layer is not limited by the water surface in most of the computational domain (see Figures 1 and 2d). In the wave Reynolds number defined by Equation (3), the representative length is a m , and its magnitude is distinctly smaller as compared with the water depth, especially in the deep sea region, resulting in a laminar regime in the tsunami source area as depicted in Figure 4.

Smooth-to-Rough Transition
In order to apply Equations (4) and (5) as a criterion for transition from smooth to rough turbulence, the maximum friction velocity, u * m , obtained from the full-range friction factor in Table 3 is used. According to the results shown in Figure 6, the roughness Reynolds numbers in four cases are all located in the region from the smooth turbulence to the transitional regime under shoaling tsunamis. The behavior has already been observed in Figure 4.

Smooth-to-Rough Transition
In order to apply Equations (4) and (5) as a criterion for transition from smooth to rough turbulence, the maximum friction velocity, * , obtained from the full-range friction factor in Table  3 is used. According to the results shown in Figure 6, the roughness Reynolds numbers in four cases are all located in the region from the smooth turbulence to the transitional regime under shoaling tsunamis. The behavior has already been observed in Figure 4.

Transitional Flow Regime in Terms of Friction Coefficient
In order to examine how the wave friction coefficient varies under the shoaling tsunami condition, values are plotted in the diagram proposed by Tanaka and Thu [18] for each computation point from the tsunami source to the shallow area (Figure 7(a1,b1)). Although transitional behavior is clearly observed from the laminar to the smooth turbulent regime in both diagrams, slightly different behavior is seen between these diagrams; at h = 10 m, the friction factor deviates from the smooth turbulent equation in Figure 7(a1), whereas the plots remain on the smooth turbulent equation up to h = 10 m in Figure 7(b1). This is due to the difference in grain size between the two test groups as summarized in Table 1. The same behavior has already been observed in diagrams for boundary layer thickness depicted in Figure 3a,b.

Comparison with Steady Friction Laws
For the purpose of comparison with a conventional shear stress evaluation method, Figure  7(a2,b2) is plotted using the following steady flow friction coefficient (derived using logarithmic velocity distribution) and (derived from Manning's formula) at each computational point.

Transitional Flow Regime in Terms of Friction Coefficient
In order to examine how the wave friction coefficient varies under the shoaling tsunami condition, f w values are plotted in the diagram proposed by Tanaka and Thu [18] for each computation point from the tsunami source to the shallow area (Figure 7(a1,b1)). Although transitional behavior is clearly observed from the laminar to the smooth turbulent regime in both diagrams, slightly different behavior is seen between these diagrams; at h = 10 m, the friction factor deviates from the smooth turbulent equation in Figure 7(a1), whereas the plots remain on the smooth turbulent equation up to h = 10 m in Figure 7(b1). This is due to the difference in grain size between the two test groups as summarized in Table 1. The same behavior has already been observed in diagrams for boundary layer thickness depicted in Figure 3a,b.

Comparison with Steady Friction Laws
For the purpose of comparison with a conventional shear stress evaluation method, Figure 7(a2,b2) is plotted using the following steady flow friction coefficient f c1 (derived using logarithmic velocity distribution) and f c2 (derived from Manning's formula) at each computational point.
where κ: the Karman constant (=0.4), z 0 : the roughness length (= k s /30), and n: Manning's roughness coefficient. Here Manning's n is assumed to be 0.025 s/m 1/3 as used in conventional computations. It is clear in Figure 7(a2,b2) that the increasing/decreasing tendency of f w from the tsunami source to the shallow area is totally opposite to that of f c1 and f c2 . In order to make a comparison among the three friction coefficients in Figure 7a,b, f c1 / f w and f c2 / f w are illustrated by taking h (the water depth) as a horizontal scale in Figure 8a,b for each computational case. By using the steady friction law, f c1 , the friction factor is underestimated in the entire computational domain, whereas f c2 gives overestimation in the shallow area. Similar analysis has been carried out by Larsen et al. [29] for their numerical analysis results of a full-scale tsunami, using a Reynolds-averaged Navier-Stokes model. They concluded that using n = 0.0125 s/m 1/3 , Equation (33) gives a tendency to under-predict the magnitude of bottom shear stress at the toe of a slope in their computational domain.
where : the Karman constant (=0.4), : the roughness length (= /30 ), and : Manning's roughness coefficient. Here Manning's is assumed to be 0.025 s/m 1/3 as used in conventional computations. It is clear in Figure 7(a2,b2) that the increasing/decreasing tendency of fw from the tsunami source to the shallow area is totally opposite to that of and .
(a) (b)  Re (2) (1) entire computational domain, whereas gives overestimation in the shallow area. Similar analysis has been carried out by Larsen et al. [29] for their numerical analysis results of a full-scale tsunami, using a Reynolds-averaged Navier-Stokes model. They concluded that using n = 0.0125 s/m 1/3 , Equation (33) gives a tendency to under-predict the magnitude of bottom shear stress at the toe of a slope in their computational domain.   Figure 9 indicates the transitional flow regime under the tsunami considered in the present study. As clarified by Tanaka et al. [4] and Tinh and Tanaka [9], depending on δ/h value, there is a transition with respect to the friction law, regardless of whether it is a wave friction zone or a steady friction zone. Moreover, according to the present study, the former can be classified into laminar, laminar-to-smooth turbulent transition, smooth turbulent regime, smooth-rough transition, and rough turbulent regime. For precise evaluation of the decay of tsunami waves, it is necessary to include such transitional behavior of bottom friction in tsunami numerical simulations. study. As clarified by Tanaka et al. [4] and Tinh and Tanaka [9], depending on /ℎ value, there is a transition with respect to the friction law, regardless of whether it is a wave friction zone or a steady friction zone. Moreover, according to the present study, the former can be classified into laminar, laminar-to-smooth turbulent transition, smooth turbulent regime, smooth-rough transition, and rough turbulent regime. For precise evaluation of the decay of tsunami waves, it is necessary to include such transitional behavior of bottom friction in tsunami numerical simulations.

Conclusions
In this study, development of a bottom boundary layer under shoaling tsunamis is investigated. In the previous study of Tinh and Tanaka [9], the flow condition was considered to be rough turbulent flow. However, this study revealed that transitional behavior prevails during the tsunami shoaling process. The most interesting outcome of the present study is that the boundary layer in the tsunami source region is located in the laminar regime, contrary to a common assumption of fully turbulent boundary layers existing in the entire range of depths. Then, as the wave propagates in the onshore direction, transition occurs from the laminar to the smooth turbulent regime. Under the condition of a hypothetical tsunami in this study, a fully rough bottom regime cannot be attained. When we evaluate wave height attenuation and sediment movement during tsunami propagation, precise identification of such behavior of the bottom boundary layer characteristic is highly required.

Conclusions
In this study, development of a bottom boundary layer under shoaling tsunamis is investigated. In the previous study of Tinh and Tanaka [9], the flow condition was considered to be rough turbulent flow. However, this study revealed that transitional behavior prevails during the tsunami shoaling process. The most interesting outcome of the present study is that the boundary layer in the tsunami source region is located in the laminar regime, contrary to a common assumption of fully turbulent boundary layers existing in the entire range of depths. Then, as the wave propagates in the onshore direction, transition occurs from the laminar to the smooth turbulent regime. Under the condition of a hypothetical tsunami in this study, a fully rough bottom regime cannot be attained. When we evaluate wave height attenuation and sediment movement during tsunami propagation, precise identification of such behavior of the bottom boundary layer characteristic is highly required.