Probabilistic characteristics of narrow-band long wave run-up onshore

8 The run-up of random long wave ensemble (swell, storm surge and tsunami) on the constant9 slope beach is studied in the framework of the nonlinear shallow-water theory in the 10 approximation of non-breaking waves. If the incident wave approaches the shore from deepest 11 water, runup characteristics can be found in two stages: at the first stage, linear equations are 12 solved and the wave characteristics at the fixed (undisturbed) shoreline are found, and, at the 13 second stage, the nonlinear dynamics of the moving shoreline is studied by means of the 14 Riemann (nonlinear) transformation of linear solutions. In the paper, detail results are obtained 15 for quasi-harmonic (narrow-band) waves with random amplitude and phase. It is shown that the 16 probabilistic characteristics of the runup extremes can be found from the linear theory, while the 17 same ones of the moving shoreline from the nonlinear theory. The role of wave breaking due to 18 large-amplitude outliers is discussed, so that it becomes necessary to consider wave ensembles 19 with non-Gaussian statistics within the framework of the analytical theory of non-breaking 20 waves. The basic formulas for calculating the probabilistic characteristics of the moving 21 shoreline and its velocity through the incident wave characteristics are given. They can be used 22 for estimates of the flooding zone characteristics in marine natural hazards. 23 24


Introduction
The flooded area size, the water flow depth, and its speed on the coast, and the coastal topography characteristics determine the consequences of marine natural disasters on the coast.The catastrophic events of recent years are well known, when tsunami waves and storm surges caused significant damage on the coast and many deaths.It is worth saying that only in 2018 two catastrophic tsunamis occurred in Indonesia, leading to the deaths of several thousand people (on Sulawesi in September and in the Sunda Strait in December).The calculations of the coastal flooding due to tsunamis and storm surges are mainly carried out within the framework of nonlinear shallow-water equations, taking into account the variable roughness coefficient for various areas of the coastal zone (Kaiser et al., 2011;Choi et al., 2012).The characteristics of the coastal destruction are determined either by using fragility curves (Macabuag et al., 2016;Park et al., 2017) or by using a direct calculation of the tsunami forces (Qi et al., 2014;Ozer et al., 2015a, b;Kian et al., 2016;Xiong et al., 2019).
The computational accuracy was tested on a series of benchmarks, including the idealized problem of the wave run-up onto the impenetrable slope of a constant gradient without friction (Synolakis et al., 2008).The nonlinear shallow-water equations for the bottom geometry of this kind are linearized by using the hodograph (Legendre) transformations.This step makes it possible to obtain a number of exact solutions describing the run-up on the coast.This approach, first suggested by Carrier and Greenspan (1958), was later on used to analyze the run-up of single and periodic waves of various shapes (Synolakis, 1987;Pelinovsky and Mazova, 1992;Carrier, 1995;Carrier et al., 2003;Tinti and Tonini, 2005;Madsen and Fuhrman, 2008;Madsen and Schaffer, 2010;Antuano andBrocchini, 2008, 2010;Didenkulova, 2009;Dobrokhotov et al., 2015;Aydin and Published by Copernicus Publications on behalf of the European Geosciences Union.S. Gurbatov and E. Pelinovsky: Probabilistic characteristics of narrow-band long-wave run-up onshore Kanoglu, 2017).Moreover, such an approach made it possible to determine the conditions for wave-breaking.The latter means the presence of steep fronts (the gradient catastrophe) within the hyperbolic shallow-water equation framework.The Carrier-Greenspan transformation was further generalized for the case of waves in an inclined channel of an arbitrary variable cross section (Rybkin et al., 2014;Pedersen, 2016;Shimozono, 2016;Anderson et al., 2017;Raz et al., 2018).In a number of practical cases, its use proves to be more efficient than the direct numerical computation within the 2-D shallow-water equation framework (Harris et al., 2015(Harris et al., , 2016)).
Due to the bathymetry variability and shoreline complexity, diffraction and scattering effects lead to an irregular shape of the waves approaching the coast.Moreover, very often the leading wave is not the maximum one.Such typical tsunami wave records on tide-gauges are well known and are not shown here.It is applied even more often to swell waves, which in some cases approach the coast without breaking (Huntley et al., 1977;Hughes et al., 2010).As a result, the statistical wave theory can be applied to such records and nonlinear shallow-water equations in the random function class can be solved.This approach was used to describe the statistical moments of the long-wave run-up characteristics in Didenkulova et al. (2008Didenkulova et al. ( , 2010Didenkulova et al. ( , 2011)).Special laboratory experiments were also conducted on irregular-wave run-up on a flat slope, the results of which are not very well described by theoretical dependencies (Denissenko et al., 2011(Denissenko et al., , 2013)).As for the field data, we are acquainted with two papers: Huntley et al. (1977) and Hughes et al. (2010), where the statistical characteristics of the moving shoreline were calculated on two Canadian beaches and one Australian beach.They confirmed the fact that the wave process on the coast is not Gaussian.In our opinion, the main problem in the theoretical model of describing the irregular-wave run-up on the shore is associated with the use of two hypotheses: (1) that the smallamplitude wave field (in the linear problem) is Gaussian and that (2) waves run-up on the shore without breaking.It is obvious, however, that in the nonlinear wave field some broken waves can always be present.They affect the distribution function tails and thus the statistical moments of the run-up characteristics.
The connection of the run-up parameters at the nonlinear stage with the linear field at a fixed point is described either in a parametric form or implicitly in the nonlinear equation (Didenkulova et al., 2010).This does not allow for using the standard methods of random processes.At the same time, it is known that this implicit equation is equivalent to a partial first-order differential equation (PDE), i.e., to the simple (the Riemann wave) equation (Rudenko and Soluyan, 1977).In statistical problems, this equation arises in nonlinear acoustics.This equation, or its generalization, the nonlinear diffusion equation called the Burgers' equation (Burgers, 1974), is the model equation in the hydrodynamic turbulence theory (Frisch, 1995).It should be noted that for the one-dimensional Burgers turbulence, as well as its three-dimensional version, used for the model description of the large-scale universe structure (Gurbatov et al., 2012), it is possible to give an almost comprehensive statistical description for certain initial conditions (Gurbatov et al., 1991(Gurbatov et al., , 1997(Gurbatov et al., , 2011;;Gurbatov and Saichev, 1993;Molchanov et al., 1995;Frisch, 1995;Woyczynski, 1998;Frisch and Bec, 2001;Bec and Khanin, 2007).In particular, the single-point and the two-point probability distributions of the velocity field and even the N-point probability distributions and, accordingly, the multipoint moment functions were found.This partially allows using a mathematical approach developed in statistical nonlinear acoustics.An experimental study of the nonlinear evolution of random quasi-monochromatic waves and their probability distributions and spectra analysis have been carried out in acoustics more than once.They confirmed the theoretical conclusions; see, for example, Gurbatov et al. (2018Gurbatov et al. ( , 2019)).
This paper is devoted to the analytical study of the probabilistic characteristics of the narrow-band long-wave run-up on the coast.Section 2 gives the basic equations of nonlinear shallow-water theory and the Carrier-Greenspan transformation, with the latter making it possible to linearize the nonlinear equations.Section 3 describes the moving shoreline dynamics when the deterministic sine wave approaches the slope.The probability characteristics of the deformed sine oscillations of the moving shoreline with a random phase are described in Sect. 4. Section 5 contains the probabilistic characteristics of the vertical displacement of the moving shoreline if the incident narrow-band wave has a random amplitude and phase.The discussion of the wave-breaking effects and their influence on the distribution of the run-up characteristics is given in Sect.6.The results obtained are summarized in Sect.7.

Basic equations and transformations
Here we will consider the classical formulation of the problem of a long-wave run-up on the constant-gradient slope in an ideal fluid (Fig. 1).The wave is one-dimensional and propagates along the x axis directed onshore.The basin depth is a linear depth function: h(x) = −αx, where α is the inclination angle tangent and point x = 0 corresponds to a fixed unperturbed water shoreline.L(t) and r(t) describe the horizontal and vertical displacement of the moving shoreline, and R(t) is the water level oscillations at x = 0.The bottom and the shore are assumed to be impenetrable.The long-wave dynamics are described by nonlinear shallow-water equations: Here, η(x, t) is the free surface elevation above the undisturbed water level, u(x, t) is the depth-averaged flow velocity (within the shallow-water theory, the flow velocity is the same on all horizons), and g is the gravity acceleration.Obviously, after introducing the total depth Eqs. ( 1) and ( 2) are a hyperbolic system with constant coefficients.This fact makes it possible to transform the system into a linear equation one by using a hodograph (Legendre) transformation, which was done in the pioneering work of Carrier and Greenspan in 1958.As a result, the wave field is described by a linear wave equation in the "cylindrical" coordinate system and all variables are expressed in terms of an auxiliary wave function (σ , λ) using explicit formulas It should be noted that the variable σ is proportional to the total water depth.
Therefore, the wave equation (Eq.4) is solved on the semiaxis σ ≥ 0, and this coordinate plays the radius role in the cylindrical coordinate system.We would like to emphasize that the point σ = 0 corresponds to a moving shoreline, and therefore the original problem, solved in the area with an unknown boundary, is reduced to the fixed area problem.
It is important to note that the hodograph transformation is valid if the Jacobian transformation is nonzero It is the case when a gradient catastrophe, identified in the framework of the shallow-water theory with the wavebreaking, does not occur.The necessary condition for the absence of wave-breaking is the boundedness and smoothness of all solutions; this question will be discussed below.We will assume that the wave approaches the coast from the area far from the shoreline (x → −∞), where the wave is linear.Then it is obvious that the function (σ , λ) can be completely found from the linear theory.The difficulty in finding the wave field in the near-shoreline area is due to the implicit transformation of the coordinates (x, t) to (σ , λ).However, for the most interesting point of the moving shoreline σ = 0 (its dynamics determines the size of the flooded area on the coast) all the formulas become explicit.In particular, from Eqs. ( 5) and ( 6) follows where r(t) and u(t) are the vertical displacement of the moving shoreline and its speed, and the functions R(t) and U (t) determine the field characteristics at the fixed point (x = 0) from the linear theory Then we add the obvious kinematic relations for the vertical displacement and velocity of the last sea point along the slope Let us note that Eq. ( 12) is identical to the so-called Riemann wave or a simple wave in a nonlinear, nondispersive medium (in particular, in nonlinear acoustics) if we consider the parameter 1/αg to be a "coordinate"; see, for example, Rudenko and Soluyan (1977) and Gurbatov et al. (1991Gurbatov et al. ( , 2011)).Moreover, Eq. ( 11) describes the integral over the Riemann wave.This analogy proves to be very useful when transferring the already known results in the wave nonlinear theory to the run-up characteristics described by Eqs. ( 11) and ( 12).Detailed calculations of the long-wave run-up on the coast were carried out repeatedly; see, for example, Carrier and Greenspan (1958), Synolakis (1987), Pelinovsky and Mazova (1992), Tinti and Tonini (2005), Madsen and Fuhrman (2008), Madsen and Schaffer (2010), Antuano andBrocchini (2008, 2010), Didenkulova (2009), Dobrokhotov et al. (2015), and Aydin and Kanoglu (2017).
It is worth mentioning that the nonlinear time transformation in Eqs. ( 11) and ( 12) leads to the shoreline oscillation distortion in comparison with the linear theory predictions.So, for large amplitudes the wave shape becomes multivalued (broken).The first moment of the wave breaking on the shoreline (the gradient catastrophe) is easily found from Eq. ( 12) by calculating the first derivative of the moving shoreline velocity and from this follows the wave-breaking condition where we have introduced the breaking parameter Br to designate the left-hand side in Eq. ( 16), which characterizes the nonlinear wave properties on the shoreline.The condition (Eq.16) can be given a physical meaning: that the breaking occurs when the last sea particle acceleration (α −1 d 2 R/dt 2 ) exceeds the component of the gravity acceleration along the shoreline (gα).As shown in Didenkulova (2009), the condition (Eq.16) coincides with Eq. ( 10) for Jacobian.It is important to emphasize that the breaking condition is unequivocally found through solving the linear problem of the wave run-up on the shore.It is determined only by the particle acceleration value on the shoreline, but it is not determined separately by the shoreline displacement or its velocity.The similar Carrier-Greenspan transformation is obtained for waves in narrow-inclined channels, fjords, and bays (Rybkin et al., 2014;Pedersen, 2016;Anderson et al., 2017;Raz et al., 2018); only the wave equation (Eq.4) and relations (Eqs.5-8) change.However, the moving shoreline dynamics are still described by Eqs. ( 11) and ( 12), valid for arbitrary cross section channels.

The moving shoreline dynamics at an initially monochromatic wave run-up
The monochromatic wave run-up on a flat slope by using the Carrier-Greenspan transformation has been studied in a number of papers cited above.Let us reproduce here the main features of the moving shoreline dynamics necessary for us to draw the statistical description below.Mathematically, the monochromatic wave run-up is described by an elementary solution of Eq. ( 4) where Q and l are arbitrary constants and J 0 is the zeroorder Bessel function.Far from the shoreline (σ → ∞) the Bessel function decreases, so the wave function becomes small.In this case, in Eqs. ( 5)-( 8) one can use approximate expressions (the "linear" Carrier-Greenspan transformation) and, using the asymptotic representation for the Bessel function, reduce Eq. ( 17) to the expression for the water surface displacement where The wave field away from the shoreline is a superposition of two waves of the same frequency and a variable amplitude a(x), which together form a standing wave.It immediately shows that the wave amplitude varies with depth according to Green's law (h −1/4 ), as it should be far from the coast.The same asymptotic result follows from the exact solution of linear shallow-water equations where R 0 is the wave amplitude at the fixed shoreline (x = 0), identified with the maximum run-up height in the linear theory.By connecting Eqs. ( 20) and ( 21), we obtain the formula for the run-up height obtained through the incident wave amplitude far from the coast Equation ( 22) allows for working further with the runup height R 0 instead of the wave amplitude far from the coast a(x).This run-up height will be considered the given value.Having determined Q and l through the incident wave parameters, we can calculate the run-up characteristics in the nonlinear theory, considering the limit of Eq. ( 17) with σ → 0 and using the Carrier-Greenspan transformation formulas (Eqs.5-8).The moving shoreline movement is determined by the parametric dependence It is convenient to introduce dimensionless variables and calculate the breaking parameter so Eqs.23 and 24 are finally rewritten in the form τ = ϕ − Br cos(ϕ), ( 27) which is another expression for Eqs. ( 11) and ( 12) if we take arising from Eq. ( 21) with x = 0. Let us note that the function z(τ , Br) is set in a parametric form, but after expressing ϕ from Eq. ( 28) and substituting it in Eq. ( 27), we can obtain the explicit expression for the function τ (z; Br).In this paper, we will use both explicit and implicit expressions of the functions describing the moving shoreline dynamics.
Figure 2 shows the moving shoreline dynamics at different wave height values in terms of the breaking parameter up to the limiting value (Br = 1).In the limit of small parameter values, the oscillations are close to sinusoidal (it is almost a linear problem).Then, with the increasing amplitude, the moving shoreline velocity gets a steep leading front, while at the moving shoreline vertical displacement a peculiar feature is formed at the wave run-down stage.As it is known, at the time of the Riemann wave breaking, the peculiarity like u ∼ t 1/3 is formed (Pelinovsky et al., 2013).Then, in the integral over the Riemann wave (at the moving shoreline displacement), this peculiar feature will have the form z ∼ t 4/3 .Thus, with the wave amplitude increase, the first breaking occurs at sea (at the run-down stage) and not on the coast.Then the breaking zone expands and moves on to the sea, but at this stage, analytical solutions based on the Carrier-Greenspan transformation become inapplicable.

Probabilistic characteristics of the initially sine wave run-up with a random phase
Let us now consider the probabilistic characteristics of the initially sine wave run-up with a random phase on the shore, assuming it to be uniformly distributed over the interval [0-2π ].These characteristics are found by using the geometric probability methods (Kendall and Stuart, 1969), so that for ergodic processes the probability density of the moving shoreline vertical displacement coincides with the relative location time of the function z(τ ) in the interval (z, z + dz) where the summation takes place at all intersection levels z(τ ).For harmonic disturbance, it is enough to restrict ourselves to considering the field over a half-period.So, for the moving shoreline vertical displacement in dimensionless variables, the derivative dτ/dz of the parametric curve (Eqs.27 and 28) can be calculated through the ratio of the derivatives dτ/dϕ and dz/dϕ We indicated here that the probability density depends on Br as a parameter.Finding cos ϕ from Eq. ( 28) for the vertical displacement, we obtain the final expression for the probability density which in the linear problem for a purely sinusoidal perturbation transforms into a well-known expression for the proba- bility distribution of a harmonic signal with a random phase (Kendall and Stuart, 1969) The probability distribution (Eq.32) for the three values of the parameter Br is shown in Fig. 3.As you can see, the probability density becomes an asymmetric function with a greater probability in the area of positive values corresponding to the wave run-up on the coast than at the run-down stage.At the ends of the interval, the probability density is unlimited throughout the entire range change of Br, since the shoreline oscillations near the maximum have a zero derivative (the moving shoreline velocity in it becomes zero).
The obtained probability density function can be used to calculate the statistical moments of the shoreline oscillations.Technically, however, it is easier to use the parametric equations (Eqs.27 and 28) and calculate all the moments.
Thus, the first moment, determines the average water level rise on the coast when the waves approach the shore (the setup phenomenon), which is commonly observed (Dean and Walton, 2009).
The second moment determines the dispersion characterizing the fluctuation range relative to the average value; it relatively weakly decreases with the growth of the parameter Br (less than 10 % for non-breaking waves).Finally, the total flooding time and its drainage time are easy to find from Eqs. ( 27) and ( 28) by finding from the Eq. ( 28) the value ϕ at which z = 0 and substituting the obtained values in Eq. ( 27) Both times change almost linearly with the increasing wave amplitude (parameter Br); see Fig. 4. It is worth noting that, in contrast to the vertical displacement, the moving shoreline velocity distribution [u = (ωR 0 /α)v], as it is easy to show, does not depend on the breaking parameter, and the probability density function is determined by the simple formula The distribution independence on the degree of nonlinearity is well known for the Riemann waves and is explained by the compensation of compression and rare fraction areas (Gurbatov et al., 1991(Gurbatov et al., , 2011)).
5 Probabilistic characteristics of a narrow-band wave run-up with a random amplitude and phase Let us consider the run-up of a quasi-harmonic wave with a random amplitude and phase on a flat slope.To do this, we will first rewrite Eqs. ( 32) and (37) for them to include the wave amplitude.It is convenient to enter the maximum height R max as the amplitude scales at which the breaking parameter turns into 1 and to use the dimensionless displacement (y = r/R max ).
Then the dimensionless amplitude is and Eq. ( 32) is converted to the form (−A < y < A) Assuming now that the wave amplitude A is a random variable, we average Eq.( 40) by using the amplitude distribution density W A (A) Equation ( 43) has an important practical meaning: by using the measured distribution of the wave amplitudes far from the coast (recomputed on run-up amplitudes in the linear theory), it is possible to obtain the distribution of the wave run-up characteristics on the coast.The only requirement imposed on the wave ensemble is that it should not contain breaking waves, which should be somehow removed from the record.It immediately follows that the Gaussian field containing large amplitude tails does not fit this requirement, and it should be modified.Therefore, we assume the amplitude distribution to be finite for A < A max = 1.The narrowband random wave field contains sine waves with almost constant frequency and random amplitude and phase.It means that if the wave amplitude is below the "breaking amplitude" A max = 1, the breaking will not be implemented in any way, and the random wave run-up will take place without any breaking.Further calculations depend on the specific type of the amplitude distribution.
Let us construct the finite amplitude distribution at which the linear field distribution is close to the Gaussian form and modify the Rayleigh distribution for wave heights in the area A < A max = 1 (Fig. 5) to make the density function distribution normalized.Here, A s is the so-called significant wave run-up height (an averaged value of one-third of the highest amplitudes).We would like to note here that it follows from Eqs. ( 11) and ( 12) that the extremal run-up characteristics in the nonlinear theory remain the same as in the linear theory.This means that the significant wave run-up height remains the same as in the nonlinear theory.When A s A max = 1, the distribution (Eq.44) transforms into the Rayleigh one, which is characteristic of the Gaussian initial distribution of a narrow-band random signal.With the help of Eq. ( 43), it becomes possible to calculate the distribution function of shoreline oscillations for the various wave energies.So, with the incident wave small amplitude (A s 1), the distribution (Eq.44) can be replaced by a simpler expression (Eq.33) and the answer is the run-up distribution characteristics in the linear theory: Besides, if A s A max = 1, the integral (Eq.45) is reduced to the Gaussian distribution where A s = 2σ y and σ 2 y is the moving shoreline oscillation dispersion.
Figure 6 shows the distribution of the run-up characteristics for different ratios of A s /A max values by Eqs. ( 43) and (44); they are shown in solid lines.Here the dashed lines show the calculation results according to the linear theory (Eq.46).As one can see, with A s /A max = 0.5 (Fig. 6a) and 0.7 (Fig. 6b), the linear distribution is close to the Gaussian one.Nonlinearity leads to the asymmetry of the distribution function density in the direction of positive values corresponding to the wave characteristics on the coast.If the undisturbed wave ensemble is made of relatively large waves (A s /A max = 1), their distribution is far from the Gaussian, both in the linear and in the nonlinear approximation.
The finite (A < A max ) power-law distribution concentrated mainly near the maximum amplitude A max can be considered to be another example of undisturbed large-amplitude waves.Figure 7 shows the graphs of the probabilistic density function of the moving shoreline displacement calculated by using Eq. ( 46) in the linear theory and Eq. ( 43) in the nonlinear theory.It is also seen in the figure that the nonlinear effects lead to a strong asymmetry towards the positive values, i.e., to the wave amplification at the run-up stage than at the rundown stage.
6 The wave-breaking effect on probabilistic run-up characteristics The theory described above is valid for non-breaking waves.
The mentioned wave ensemble, strictly speaking, cannot be the Gaussian one, as it always has unlimited tails in the probability density function.Let us briefly discuss what the formulas obtained for non-breaking waves lead to in the presence of broken waves.Figure 8 shows the parametric curve (Eqs.27 and 28) when Br = 2. Formally, the curve became multivalued in the range of negative values corresponding to the maximum water outflow from the coast.We have already indicated that the probability density function of the moving shoreline vertical displacement W (ξ ) coincides with the relative residence time ξ(t) of the function in the interval (ξ , ξ + dξ ), which is calculated by Eq. ( 17).In contrast to negative cut-off bias values, in the area of positive values there is no ambiguity, and therefore all the calculations can be carried out by using the formulas described above.The calculation example with Br = 2 and r > −0.5 (in the zone of one-value solution) is shown in Fig. 9.However, these results should be treated with caution.If Br > 1 the Jacobian breaks down seawards of the shoreline.This may affect the probabilistic distribution on the positive side.This important issue requires going beyond the theory discussed in this article.

Discussion and conclusion
In this paper, we study the run-up of irregular narrow-band waves with a random envelope (swell, storm surges, and tsunami) on a beach of a constant slope.The work was carried out in the framework of the nonlinear wave theory with one important assumption: there should be no breaking waves in the wave ensemble.This restriction is quite strict for field and laboratory conditions, but nevertheless, there are cases when it is performed.For instance, 75 % of historical tsunami waves climbed onto the coast with no breaking (Mazova et al., 1983).In the experiments performed in the Warwick University tank and in the Large Tank in Hanover (Denissenko et al., 2011(Denissenko et al., , 2013) ) this condition was fulfilled.
The wave nonlinearity at the run-up stage leads to increased deviations from Gaussianity, as might be expected from general considerations.Nevertheless, it is shown that the probability distribution of the moving shoreline velocity does not depend on the wave nonlinearity and can be calculated within the linear theory framework.The same conclusion can be drawn about the distribution of the extreme runup characteristics (the moving shoreline displacement and speed), which, in fact, has been discussed previously (Didenkulova et al., 2008).However, the probabilistic density function of the moving shoreline displacement differs from that predicted in the linear theory framework.It is described by Eq. ( 43) by using either the theoretical or the measured distribution of the incident wave amplitudes.The paper gives the calculation results of the probabilistic run-up characteris-tics with the modified Rayleigh distribution for wave amplitudes.
The wave-breaking leads to the inapplicability of the wave run-up theory based on the Carrier-Greenspan transformation.If, nevertheless, the share of large amplitude waves is small, the breaking occurs mainly at the run-down stage, having little effect on the long-wave coastal flooding characteristics (see Sect. 6).This question, however, requires a special study based on direct numerical solutions of the shallowwater equations or their nonlinear-dispersive generalizations.
Finally, it is worth noting that we considered the narrowband wave run-up with a random amplitude and phase.As far as the random waves with a wide spectrum are concerned, they may be the problem for further consideration.
The obtained probability density functions of the vertical displacement of the moving shoreline are useful to compute statistical characteristics of flooding time and force on coasts and constructions, which are necessity for the mitigation of natural marine hazards.Now, in practice, various generalizations of shallow-water equations are used to analyze tsunami run-up including the wave dispersion; see, for instance, Løvholt et al. (2012).The wave dissipation is a quadratic dissipative term that prevents us from getting analytical results, so its influence on statistical characteristics should be investigated in the future.

Figure 2 .
Figure 2. The moving shoreline dynamics (a) and its velocity (b) in the case of the incident monochromatic wave for different breaking parameter values Br (0: the dotted line; 0.5: the dashed line; 1: the solid line).

Figure 3 .
Figure 3.The probability density of the moving shoreline vertical displacement for the initially sine wave run-up at Br = 0 (the dotted line), 0.5 (the dashed line), and 1 (the solid line).

Figure 4 .
Figure 4.The total flooding time (the solid curve) and the drainage time (the dashed curve) depending on the parameter Br.

Figure 5 .
Figure 5.The modified Rayleigh distribution (Eq.43) for different distribution values A s /A max (0.5: the dotted curve; 0.7: the dashed line; 1: the solid line).

Figure 6 .
Figure 6.The probabilistic density function of the vertical shoreline displacement in the nonlinear theory (the solid lines) and in the linear theory (the dashed lines) for different A s /A max : (a) 0.5 values, (b) 0.7, and (c) 1 .

Figure 7 .
Figure 7.The probabilistic density function of the shoreline vertical displacement in the linear theory (the dashed line) and nonlinear theory (the solid line).

Figure 8 .
Figure 8.The parametric curve (Eqs.27 and 28) with Br = 2 (the solid curve) in comparison with the linear problem with Br = 0 (the dashed line).