Next Article in Journal
Plankton Resting Stages Distribution in Bottom Sediments along the Confinement Gradient of the Taranto Sea System (Ionian Sea, Southern Italy)
Previous Article in Journal
Resting Cysts of the Toxic Dinoflagellate Gymnodinium catenatum (Dinophyceae) Ubiquitously Distribute along the Entire Coast of China, with Higher Abundance in Bloom-Prone Areas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Bed Variations on Linear Wave Propagation beyond the Mild Slope Condition

Instituto de Ciencias del Mar (CSIC), 08003 Barcelona, Spain
J. Mar. Sci. Eng. 2024, 12(9), 1652; https://doi.org/10.3390/jmse12091652
Submission received: 30 July 2024 / Revised: 27 August 2024 / Accepted: 11 September 2024 / Published: 14 September 2024
(This article belongs to the Section Ocean Engineering)

Abstract

:
As water waves travel from deep to shallow waters, they experience increased nonlinearity and decreased dispersion due to the reduced water depth. While the impact of bed slope on wave propagation celerity is documented, it is often overlooked in commonly used depth-integrated wave models. This study uses the WKB approximation and solves higher-order slope-related terms to analyze the influence of varying depth, including the gradient and laplacian of water depth. One result is an extended longwave for linear wave reflection and transmission on a ramp to deeper waters. The main outcome and focus of this work, however, is a new, simple analytical expression for linear dispersion that includes bed variations. The results are applied to two cases: wind-generated water wave propagation in the nearshore, emphasizing corrections to bathymetry inversion methods, and tsunami propagation over the continental slope, highlighting the limitations of neglecting slope on dispersion and the significant role of the Laplacian of water depth.

1. Introduction

As wind-generated water waves travel from the depths of the ocean toward shallower coastal regions, the diminishing depth has significant implications for their propagation. In deep waters, nonlinear effects are usually minimal, and the waves remain unaffected by the seabed, meaning that their propagation is independent of depth. In this regime, dispersive effects dominate, leading to different celerities for different wave frequencies. However, as waves enter intermediate to shallow waters, where they begin to feel the bottom—i.e., where wave propagation is influenced by water depth—nonlinear effects can become significant, and dispersive effects are reduced. In fact, in shallow waters, celerity depends solely on water depth and becomes independent of wave frequency (i.e., there is no dispersion).
For both practical and scientific purposes, water wave propagation models must be computationally efficient. This necessity leads to a key simplification: depth-integrating the governing equations to derive a system of two-dimensional equations [1,2,3,4]. Additionally, it is common practice to assume that the flow is inviscid and irrotational, which is a good approximation except within the thin bottom boundary layer and during wave breaking. The effects of the boundary layer and wave breaking are typically incorporated into the energy equation [5,6,7]. Further simplifications involve time-averaging the governing equations to bypass resolving the wave phase, allowing a focus on wave statistics over large spatial domains [8].
Regarding phase-resolving equations, which is the focus here, although depth-integrated and irrotational, they must still effectively capture nonlinear, dispersive, and bed slope effects. These models can be broadly categorized into two largely complementary families: those derived from Airy theory and those based on the Shallow-Water Equations (SWEs).
Airy theory is fully dispersive, i.e., valid for arbitrary depths, but it is applicable only to small-amplitude waves [2,9]. Further, it is limited to flat beds. Despite its simplicity, Airy theory provides linear equations that yield useful approximate expressions for velocity profiles, pressure fields, and other wave characteristics. One fundamental result from Airy theory is the dispersion relationship, which relates the wave angular frequency, ω = 2 π / T ( T = wave period), the wavenumber, k = 2 π / L ( L = wavelength), and the water depth, h. This relationship, valid in principle for flat beds, is given by [10] (Burnside)
ω 2 = g k tanh k h ,
where g is the gravitational acceleration. Because of the asymptotic behavior of the hyperbolic tangent, it is commonly considered that k h π indicates deep water conditions (where ω 2 = g k ), while k h π / 10 indicates shallow waters (where ω 2 = g k 2 h ) [2].
There are two main extensions of the Airy linear theory. The first is Stokes theory [1,11], which introduces a perturbative approach based on wave amplitude to account for nonlinear effects to some extent. While Stokes theory has limitations in shallow waters, it is valuable for applications such as the harmonic decomposition of waves [1]. The second extension involves the Mild-Slope Equations (MSEs) [12,13] and related variants, including the “Modified” and “Complementary” MSEs, e.g., [14,15,16,17]. This family of equations is particularly useful for addressing bed-related phenomena such as wave shoaling, refraction, and diffraction. Given that the MSEs implicitly account for two characteristic lengths—namely, the wavelength and the length scale over which the bed bathymetry changes—the WKB (Wentzel–Kramers–Brillouin) perturbative method is particularly well suited for their derivation [1,18]. Notably, although the MSEs account for bed slope effects, the dispersion Equation (1) remains unchanged. This is because, at the perturbative order considered, the bed can be considered locally flat with respect to dispersion.
On the other hand, the SWEs can handle arbitrary nonlinearity (large-amplitude waves) but only in shallow waters, where the water depth is small compared to the wavelength [19]. Starting from the SWEs and using perturbative techniques, Boussinesq-type equations (BTEs) are derived to incorporate dispersive effects, extending the range of applicability to deeper waters [20,21,22,23,24]. The ability of the BTEs to manage dispersive conditions is typically assessed by comparing the performance of the linearized BTEs against the linear Airy theory, specifically in terms of wave celerity and shoaling behavior [24,25].
The above models do not fully capture the influence of bed slope on wave dispersion, particularly in fully dispersive conditions. This work addresses this gap by introducing a new formulation that incorporates bed slope effects into the dispersion Equation (1). In order to achieve this, we apply the WKB approximation to a sufficient order, enabling us to accurately capture the effects of bed slope on the dispersion. Actually, the limitations of the Burnside equation are already well documented. Broadly speaking, if σ represents the slope, its influence can be neglected only when [18]
σ k h ,
is small. Importantly, it is not just σ that must be small, but rather σ / k must be small compared to h, where σ / k = 2 π σ L is a measure of the depth change over a wavelength. Thus, the same slope σ can be considered mild or steep depending on k h , with slopes appearing milder in deep waters (large k h ) and steeper in shallow waters (small k h ). For values of k h ranging from approximately 0.46 to 0.86 , Booij [26] found that the mild slope hypothesis holds when σ 1 / 3 (expressed in terms of ω 2 h / g instead of k h in the original work). Therefore, the mild slope condition can be more precisely stated as
σ k h 0.38 = 1 / 3 0.86 .
For very shallow waters and a sloped plane bed, the problem has already been addressed by Ehrenmark [27] (see also [28]). The expression by Ehrenmark [27] reads
ω 2 = g k tanh arctan σ σ k h ,
through an asymptotic analysis. Since the above holds only for very shallow waters, the hyperbolic tangent tends to the identity, and
k = ω g h σ arctan σ = k 0 1 + σ 2 6 + O ( σ 4 ) ,
where k 0 = ω / g h is the Burnside solution for shallow waters. According to Equation (3), k k 0 , so the celerity c = ω / k is reduced by the slope.
The previous approach is valid only for the asymptotic case of very shallow waters. For the fully dispersive regime, Ge et al. [29] recently presented a dispersion relationship that incorporates bed slope effects, but it is limited to the specific case of sinusoidal bathymetry. In contrast, we derive a more general expression that is fully dispersive and applicable to arbitrary bathymetry. We employ a WKB approach following [1,18] but extend the analysis to higher-order terms. The resulting equation bears a striking resemblance to that proposed by Ge et al. [29], with the critical difference of an opposite sign—a distinction that is highly significant. We then apply the proposed dispersion equation to several important cases to demonstrate the impact of accounting for (or neglecting) bed slope variations in the dispersion relationship.

2. Governing Equations

Throughout the work, x = ( x , y ) will refer to the horizontal coordinates, and z to the vertical one, with z = 0 at the mean sea water level, z = η at the free surface, and z = h at the bottom (Figure 1 for the one-dimensional case). We use the notation usually employed in water-wave textbooks as, e.g., those by Dean and Dalrymple [2] or Svendsen [3]. Also, as usual in depth-integrated models and many other wave propagation models, the flow is considered inviscid and irrotational [1,2,3,18], and the potential of the velocity is denoted with Φ . In assuming, hereinafter, a rigid bed, the linearized equations for Φ are [18]
2 Φ + 2 Φ z 2 = 0 , h z 0 ,
Φ · h + Φ z = 0 ,     z = h ,
2 Φ t 2 + g Φ z = 0 z = 0 ,
where ∇ stands for the horizontal gradient ( / x , / y ) .
Equations (4)–(6) are, respectively, the mass conservation, the kinematic boundary condition at the bottom, and a combination of the kinematic and the dynamic boundary conditions at the free surface. The Bernoulli equation, as well as the kinematic and dynamic boundary conditions at the free surface, is avoided so that the equations are expressed exclusively using the potential Φ [3,18]. The linearized kinematic and dynamic boundary conditions at the free surface are, for future use,
η t Φ z = 0 z = 0 ,
Φ t + g η = 0   z = 0 .
The depth-integrated continuity and momentum equations are, e.g., [18],
· h 0 Φ + η t = 0 ,
g η + Φ z = 0 t = 0 ,
and the depth-integrated and time-averaged energy equation is
· h 0 Φ t Φ d z = 0 ,
where · stands for the time average (i.e., over one wave period anticipating the time-harmonic flow conditions).

2.1. Equations for Time-Harmonic Free-Surface Flows

Assuming a time-harmonic solution, we consider
Φ = Φ 1 exp ( i ω t ) + complex   conjugate ,
η = η 1 exp ( i ω t ) + complex   conjugate ,
with ω as the wave angular frequency; Φ 1 and η 1 , two complex valued functions of x and z; and i = 1 , the imaginary unit hereinafter. From free-surface boundary condition (8), it is
η 1 = i ω g Φ 1 , z = 0 ,
so that η 1 is known if Φ 1 , which is the focus herein, is resolved.
As any complex variable, the time-independent Φ 1 ( x , z ) can be expressed in its polar form as
Φ 1 = A ( x , z ) exp ( i θ ( x , z ) ) ,
with A and θ as two real functions (the amplitude and the phase, respectively). Introducing the above Equations (12) and (15) into Equations (4)–(6), we obtain a set of equations with the following real part:
2 A A k 2 + 2 A z 2 A θ z θ z = 0 , h z 0 ,
A · h + A z = 0 ,   z = h ,
ω 2 A + g A z = 0 ,     z = 0 ,
with k 2 = θ · θ . The imaginary part is
2 A · θ + A 2 θ + 2 A z θ z + A 2 θ z 2 = 0 , h z 0 ,
A θ · h + A θ z = 0 , z = h ,
g A θ z = 0 ,         z = 0 .
Similarly, in recalling Equation (14), the depth-integrated continuity and energy Equations (9) and (11) read
· h 0 ( A + i A θ ) exp ( i θ ) d z + ω 2 g A z = 0 exp ( i θ z = 0 ) = 0 ,
· h 0 A 2 θ d z = 0 ,
The momentum Equation (10) is automatically satisfied in this time-harmonic flow case. In recalling Equations (12) and (13), Equation (10) becomes
g η 1 i ω Φ 1 , z = 0 = 0 ,
which is automatically satisfied if, as stated in Equation (14), η 1 = i g 1 ω Φ 1 , z = 0 .

2.2. The Imaginary Part: Vertical Profile of the Phase

In multiplying by A, Equation (19) is
· ( A 2 θ ) + z A 2 θ z = 0 ,
so that, by integrating in z from h to z and using bottom boundary condition (20), we obtain
θ z = 1 A 2 · h z A 2 θ d z .
Taking into account the free surface boundary condition (21), the above reduces to (energy) Equation (23) for z = 0 . Equation (25) provides insight into the vertical variation in the phase θ . In the derivation of the Mild-Slope Equation and relatives, the phase is typically assumed to be uniform with respect to the vertical coordinate, i.e., θ / z = 0 . Solving the (non-uniform) vertical profile to a higher-order accuracy will, however, be necessary to obtain the modified dispersion relationship.

3. The WKB Approximation

Equations in Section 2.1, which do not require any condition on the bottom slope being small, do not have analytical solutions other than trivial solutions. To obtain analytical non-trivial solutions that allow us to understand some relevant aspects of wave propagation (wave celerity, wave shoaling, ...), the bed is usually considered to vary slowly to obtain MSEs and relatives. In this Section 3, we introduce the WKB approximation, the same method used to obtain MSEs, which allows us to rigorously apply this assumption.
Following the usual WKB approximation in water wave propagation problems [18,30], with χ as a characteristic small value of the bed slope, the slow spatial coordinates x = χ x are introduced. The space derivatives relative to this variable are = χ 1 , and consequently,
h = χ 1 h = χ 1 𝒪 ( χ ) = 𝒪 ( 1 ) ,
i.e., x allows for the depth variations for that coordinate to be order one. Further, A and θ are expanded in WKB as (see Dingemans [30], pp. 60–63, for a detailed rationale behind these equations)
A = A 0 + χ 2 A 1 + 𝒪 ( χ 4 ) , θ = χ 1 ( θ 0 + χ 2 θ 1 + 𝒪 ( χ 4 ) ) ,
assuming that A i and θ i are slow variables, so that their derivatives and 2 are order one. From the above expansions, we obtain, e.g.,
A = χ A 0 + χ 3 2 A 1 + 𝒪 ( χ 5 ) , 2 A = χ 2 2 A 0 + 𝒪 ( χ 4 ) ,
and
θ = θ 0 + χ 2 θ 1 + 𝒪 ( χ 4 ) .
From the last equation, k = θ and k 0 = θ 0 are the wavenumber vector and its leading order approximation, respectively, with k = k 0 + 𝒪 ( χ 2 ) , and further,
k = θ · θ = θ 0 · θ 0 + 2 χ 2 θ 0 · θ 1 + 𝒪 ( χ 4 ) = θ 0 · θ 0 + χ 2 θ 0 · θ 1 θ 0 · θ 0 + 𝒪 ( χ 4 ) ,
i.e.,
k = k 0 + χ 2 k 1 + 𝒪 ( χ 4 ) ,
where k 0 2 = θ 0 · θ 0 and k 1 = k 0 1 θ 0 · θ 1 . The main goal of this work is, actually, to obtain the expressions for k 0 and k 1 , and therefore the (linear) dispersion relationship including the influence of the bed variations.
For future use, it is also of interest that, from Equation (25),
θ z = χ A 0 2 · h z A 0 2 θ 0 d z + 𝒪 ( χ 3 ) = 𝒪 ( χ 1 ) ,
so that, in recalling the expansions in Equation (26),
θ 0 z = ( χ + 𝒪 ( χ 2 ) ) θ z = 𝒪 ( χ 2 ) ,
and therefore,
k 0 z = z θ 0 = θ 0 z = 𝒪 ( χ 2 ) .
In order to obtain Equation (28), it has been considered, e.g., that
1 A 2 = 1 A 0 2 + 2 χ 2 A 0 A 1 + 𝒪 ( χ 4 ) = 1 A 0 2 + 𝒪 ( χ 2 ) .
These kinds of manipulations will often be carried out without being explicitly shown.

4. The Leading-Order Solution

By introducing Expansion (26) into Equations (16)–(18) and noting θ / z = 𝒪 ( χ 1 ) from Equation (28), we obtain, to the leading order,
A 0 k 0 2 + 2 A 0 z 2 = 0 , h z 0 ,
A 0 z = 0 ,   z = h ,
ω 2 A 0 + g A 0 z = 0 , z = 0 .

4.1. Solution for the Equations

The solution of the above equations, and all the results in this Section 4.1, are well known and given here without detailed explanations [3]. The amplitude A 0 is
A 0 = K c 0 cosh q 0 , [ q 0 = k 0 ( z + h ) ] .
From the free surface condition (32), we obtain
ω 2 = g k 0 tanh μ 0 ,
which is the leading-order Burnside [10] dispersion relationship. Equations (33) and (34), with constant K c 0 , constitute the Airy solution for water-wave propagation over flat beds [2]. Many interesting phenomena can already be understood from these two equations (deep and shallow water conditions, good approximations of the velocity profile and the pressure field, standing waves, ...).
Even within this leading-order solution, the variations in K c 0 , which is a representation of the wave height, can be estimated due to the water depth changes. From Equation (23), the leading order is
· ( K c 0 2 ( sinh 2 μ 0 + 2 μ 0 ) s 0 ) = 0 , [ μ 0 = k 0 h ] ,
where s 0 = k 0 / k 0 is the unit vector in the direction of the wave propagation. Expression (35) relates K c 0 with the water depth, embedded in μ 0 . It is a fundamental equation in the MSE family.
From Equation (15), Φ 1 = K c 0 cosh q 0 exp ( i θ ) , and from Equation (14),
η 1 = i ω g Φ 1 , z = 0 = i ω g K c 0 cosh μ 0 exp ( i θ z = 0 ) ,
or
η 1 = i H 0 4 exp ( i θ z = 0 ) , with H 0 = 4 ω g K c 0 cosh μ 0 ,
which is the wave height. In using the wave height, energy Equation (35) is
· ( H 0 2 c g 0 ) = 0 , with c g 0 = ω k 0 s 0 ,
which is the group celerity vector. The latter equation has been widely used to understand the shoaling in terms of the wave height as it approaches the shore, and Green’s law is a particular case for shallow waters [2]. Introducing a drag term on the right-hand side has been used to understand the wave attenuation associated with bottom friction [7].

4.2. Spatial Derivatives

To obtain higher-order approximations of the solution in the following Section 5 and Section 6, the gradient and the Laplacian of the main functions obtained in the leading-order solution will be required. For any function f of h, it is
f = d f d h h , 2 f = d f d h 2 h + d 2 f d h 2 | h | 2 ,
and therefore, d f / d h and d 2 f / d h 2 suffice to obtain f and 2 f as functions of h and 2 h . Following Schäffer and Madsen [31], we define
α f = h f d f d h , β f = h 2 f d 2 f d h 2 ,
which are dimensionless, by construction, and it turns out that for the functions we are interested in, they depend solely on μ 0 . The gradient and Laplacian can be expressed as
f = α f f h h , 2 f = α f f h 2 h β f f h 2 | h | 2 .
As shown in Appendix A,
α k 0 = G 0 1 + G 0 , G 0 = 2 μ 0 sinh 2 μ 0 ,
α K c 0 = G 0 M 0 2 ( 1 + G 0 ) 2 , M 0 = cosh 2 μ 0 + 1 .
Other functions, such as α μ 0 , β k 0 , and β K c 0 , are also included in Appendix A. In addition, as q 0 = k 0 ( z + h ) , we have
d q 0 d h = d k 0 d h ( z + h ) + k 0 = d k 0 d h q 0 k 0 + k 0 = α k 0 q 0 + μ 0 h ,
and after manipulation,
A 0 h = K c 0 h ( δ c 0 cosh q 0 + ( δ s 1 q 0 + δ s 0 ) sinh q 0 ) ,
2 A 0 h 2 = K c 0 h 2 ( ( ϵ c 2 q 0 2 + ϵ c 1 q 0 + ϵ c 0 ) cosh q 0 + ( ϵ s 1 q 0 + ϵ s 0 ) sinh q 0 ) ,
where the functions δ and ϵ are shown in Table 1.

5. Solution to 𝒪 ( χ 1 )

Section 5 presents results for θ / z based on the leading-order solution, which, as shown in Equation (28), is 𝒪 ( χ 1 ) . The implications of this solution on the classical problem of the wave reflection over a ramp are hinted at.

5.1. The Vertical Profile of the Phase

In recalling the expansions and introducing the solution into Equation (25),
θ 0 z = χ 2 K c 0 2 cosh 2 q 0 · h z K c 0 2 cosh 2 q 0 k 0 d z + 𝒪 ( χ 4 ) .
In taking into account that k 0 varies slowly in z—see Equation (29)—the above can be integrated to obtain θ 0 / z and, further, to obtain
θ 0 = θ 0 b + χ 2 q 0 2 μ 0 ( 2 α K c 0 tanh q 0 + α k 0 q 0 2 μ 0 ) h · s 0 + 𝒪 ( χ 4 ) ,
with s 0 being the above-mentioned unit vector in the direction of the wave propagation, and θ 0 b as a function of x only—not of z. In recalling that θ = χ 1 θ 0 +
θ = θ b + q 0 2 μ 0 ( 2 α K c 0 tanh q 0 + α k 0 q 0 2 μ 0 ) h · s 0 = θ b + 𝒪 ( χ ) ,
with θ b = χ 1 θ 0 b .

5.2. Transmission and Reflections on a Ramp

The above results, particularly the vertical profile obtained in Equation (43) for the phase, which has not been described in the literature to the author’s knowledge, allow us here to extend the longwave theory of wave transmission and reflection on a ramp, e.g., ref. [1] to include the dispersion effects. This section is not required for the rest of the work and can be skipped if the reader’s interest is solely in the dispersion relationship including the bed slope influence. Moreover, it is not the author’s intention to thoroughly analyze all the implications of the previous result beyond a brief illustration of the reflection problem over a ramp, as the primary objective here is to establish a dispersion relation that incorporates the effects of the bed slope. Analyzing the impact of the above results on the Bragg resonance analysis, e.g., falls out of the scope of this work.
In introducing the results in Section 4.1 and Section 5.1 into Expansion (26), the solution is
Φ 1 = g H 0 4 ω cosh q 0 cosh μ 0 exp ( i θ ) , η 1 = i H 0 4 exp ( i θ z = 0 ) ,
where θ is given in Equation (43). This is the solution to order 𝒪 ( χ 1 ) .
For the analysis of wave transmission and reflection on a ramp for the one-dimensional case, we now follow the procedure as in Mei [1], but considering the vertical variations in the phase θ . A key issue is the fact that the water flux,
Q 1 = h 0 Φ 1 x d z ,
introduces the function θ , which is now a function of z through q 0 , into Φ 1 . In using the results above, the continuity and momentum depth-integrated equations, instead of the Saint-Venant equations used by Mei [1] for longwaves, are as follows:
Q 1 x i ω η 1 = 0 ,
i ω Q 1 μ 0 tanh μ 0 + i f s d h d x + g h d η 1 d x = 0 .
where s = + 1 if the wave propagates in the direction of x, and s = 1 otherwise, and
f = ( α k 0 G 0 α K c 0 ) tanh μ 0 + μ 0 α K c 0 ( tanh 2 μ 0 1 ) tanh 2 μ 0 .
Note that in using the definitions of the functions involved in f, in shallow waters, we have, for the first term in Equation (46),
lim μ 0 0 μ 0 tanh μ 0 + i f s d h d x = 1 i μ 0 3 s d h d x 1 ,
so that the bed slope effect becomes negligible in shallow waters, with Equations (45) and (46) reducing to the Saint-Venant equations. Departing now from Equations (45) and (46) and following the same reasoning as that by Mei [1], the solution is (here “+” stands for the incident component traveling in the direction of x, and “−” stands for the reflected one)
η 1 = E + exp ( i ( + θ b + S + ) ) + E exp ( i ( θ b + S ) ) c g 0 ,
where S + and S are constants of integration used to arbitrarily set lags between the incident and the reflected components (often omitted), and c g 0 = ω / k 0 is the group celerity. Also, E + and E , the energies, satisfy the equations
d E + d x = r h d h d x exp ( i ( 2 θ b + S S + ) ) E ,
d E d x = r h d h d x exp ( i ( + 2 θ b + S + S ) ) E + .
with
r = μ 0 α K c 0 ( 1 tanh 2 μ 0 ) tanh μ 0 .
In Equations (47)–(49), the phase θ b can be replaced with θ , because the differences between them are 𝒪 ( χ 1 ) , therefore corresponding to the following order of accuracy in the WKB perturbative framework. Figure 2 shows the new function r ( μ 0 ) ; note that in shallow waters, r tends to 1 / 4 , and Equations (48) and (49) reduce to Equation (4.5.13) in Mei [1]—after some rewriting of Mei’s equations.

6. Solution to 𝒪 ( χ 2 )

The main result of this work, to find a dispersion relationship to express k = k 0 + χ 2 k 1 as a function of the wave angular velocity and the water depth while taking into account the varying bed, is presented in this section by solving the order 𝒪 ( χ 2 ) of Equations (16)–(18).

6.1. Equations for A 1

In introducing Expansion (26) into the partial differential Equation (16) and retaining all the terms up to order 𝒪 ( χ 2 ) , and after some rearrangements,
A 0 k 0 2 + 2 A 0 z 2 I + χ 2 2 A 0 2 A 0 k 0 k 1 A 1 k 0 2 + 2 A 1 z 2   χ 2 A 0 θ 0 z θ 0 z = 𝒪 ( χ 4 ) .
The last term on the left-hand side includes θ 0 / z = 𝒪 ( χ 2 ) , so it is 𝒪 ( χ 2 ) and is to be retained. Furthermore, according to the leading-order solution in Equation (30), the first two terms, I = A 0 k 0 2 + 2 A 0 / z 2 , cancel out if the are evaluated using k 0 . However, we are to find the solutions using (and evaluating A 0 at) k = k 0 + χ 2 k 1 . A Taylor expansion gives
I = A 0 ( k = k 0 + χ 2 k 1 ) k 0 2 + 2 A 0 z 2 ( k = k 0 + χ 2 k 1 ) =   A 0 k 0 2 + 2 A 0 z 2 χ 2 k 1 A 0 k k 0 2 + χ 2 k 1 k 2 A 0 z 2 ,
where all the evaluations on the right-hand side are at k 0 , so that A 0 k 0 2 + 2 A 0 / z 2 = 0 on the right-hand side of Equation (51), and therefore,
I = χ 2 A 0 k k 0 2 + k 2 A 0 z 2 k 1 ,
with the evaluations of A 0 and its derivatives performed at k 0 . Note, however, that the evaluations of these terms could also be carried out at k = k 0 + χ 2 k 1 , since they add terms that are, overall, 𝒪 ( χ 4 ) and negligible.
Since A 0 = K c 0 cosh q 0 , the above expression (52) is also, after manipulation,
I = 2 χ 2 A 0 k 0 k 1 ,
so that Equation (50) is also (recall that θ 0 / z = 𝒪 ( χ 2 ) , and thus, all the terms below are the same order)
A 1 k 0 2 + 2 A 1 z 2 = 2 A 0 + χ 4 A 0 θ 0 z θ 0 z .
The above corresponds to the terms 𝒪 ( χ 2 ) of continuity Equation (16). In proceeding similarly with boundary conditions (17) and (18), the equations for A 1 are
A 1 k 0 2 + 2 A 1 z 2 = 2 A 0 χ 4 A 0 θ 0 z θ 0 z , h z 0 ,
A 1 z = A 0 · h ,   z = h ,
ω 2 A 1 + g A 1 z = g K c 0 m k 1 , z = 0 ,
with
m = μ 0 tanh μ 0 sinh μ 0 sinh μ 0 μ 0 cosh μ 0 .
Equations (53) and (54) will allow us to solve for A 1 , and then k 1 will be obtained from Equation (55).

6.2. Analytical solutions for θ 0 / z = 0

The forcing term including the squared θ 0 / z in Equation (53) prevents analytical solutions to exist. For illustration and also practical purposes, it is useful to first introduce the solution of the problem ignoring this forcing term.
The gradient and the Laplacian of A 0 in Equations (53) and (54) are
A 0 = A 0 h h , 2 A 0 = A 0 h 2 h + 2 A 0 h 2 | h | 2 ,
where A 0 / h and 2 A 0 / h 2 are known from Equations (39) and (40). Differential Equation (53) happens to have an analytical solution. Already imposing boundary condition (54), the solution is
A 1 K c 0 = δ s 1 4 q 0 2 δ s 0 2 q 0 cosh q 0 h 2 h μ 0 2 + 2 δ c 0 + δ s 1 4 q 0 + μ 0 2 sinh q 0 h 2 h μ 0 2 + ϵ s 1 + ϵ c 2 4 q 0 2 + 2 ϵ s 0 + ϵ c 1 4 q 0 cosh q 0 | h | 2 μ 0 2 + ϵ c 2 6 q 0 3 ϵ c 1 4 q 0 2 + 2 ϵ c 0 + ϵ s 1 ϵ c 2 4 q 0 α k 0 μ 0 2 sinh q 0 | h | 2 μ 0 2 + K c 1 cosh q 0 .
This new contribution of 𝒪 ( χ 2 ) to the velocity profile is not considered in the work by Ge et al. [29], and it is also different than the depth functions considered in Modified and Complementary MSEs [14,15,16]. Solving the horizontal variations om K c 1 (constant in z) in Equation (57) requires, just as in the leading order for K c 0 , the use of energy Equation (23), even though, again as with the leading-order case, the dispersion relationship does not require this constant K c 1 to be solved. The wavenumber k 1 is obtained from Equation (55), i.e.,
k 1 = 1 g K c 0 m ω 2 A 1 + g A 1 z , at   z = 0 ,
or, using the leading-order dispersion relationship ω 2 = g k 0 tanh μ 0 —allowed since the errors introduced in k 1 are 𝒪 ( χ 2 ) , i.e., 𝒪 ( χ 4 ) in k = k 0 + χ 2 k 1 —and using the above solution (57) for A 1 , after manipulation, we have
k 1 = k 0 ( F 1 h 2 h + F 2 | h | 2 ) ,
where F i denotes the expressions of μ 0 , which are shown in Appendix B and plotted in Figure 3. As anticipated, K c 1 is not involved in Equation (59) for k 1 . By introducing the above into k = k 0 + χ 2 k 1 , we finally obtain
k = k 0 ( 1 + F ) ,
where F, which can be seen as a correction factor due to the bed slope, is
F = F 1 h 2 h + F 2 | h | 2 .
The proposed expression (60), with F in Equation (61), trivially reduces to the Burnside solution when 2 h = 0 and | h | = 0 , for flat beds. Functions F 1 and F 2 , in Figure 3 ( F 1 shown, since F 1 < 0 and logarithmic scales are used), are so that | F i | < 10 2 for μ 0 2.3 , decreasing exponentially as μ 0 increases (Figure 3B). In other words, in deep waters, the dispersion relationship is not affected by the variations in the bed, as it is not by the depth itself. For practical purposes, we will consider that the influence of the slope is negligible for μ 0 > 2.3 .
On the other hand, for μ 0 0.1 , the slope in the log–log plot (Figure 3A) is 2 , i.e., both functions are μ 0 2 . Specifically, they are
F 1 1 8 μ 0 2 , F 2 + 1 32 μ 0 2 , [ μ 0 0.1 ] ,
so that in shallow waters,
F h 2 h 8 μ 0 2 + | h | 2 32 μ 0 2 , [ μ 0 0.1 ] .
Given that the perturbative approximations are valid as long as the perturbation created, here F, is relatively small, say | F | < F c (a critic value), the above expression (62) also allows us to theoretically find the limit for the application of the present analysis in shallow waters. For the particular case of the bed being a sloped plane, when 2 h = 0 , from the above expression,
| h | μ 0 32 F c | h | μ 0 2.5   if   F c 0.2 .
This condition has the same structure as the one for MSEs in Equation (2) and represents a significant increase in the validity domain. Also, for sloped planes, Figure 4 includes two relevant isolines for F = F 2 | h | 2 , showing the zone of present interest in the plane { | h | , μ 0 } , which corresponds to 0.001 F F c with F c 0.2 ; for F > F c , the influence is too high, while for F < 0.001 , it can be considered negligible. The domain of the asymptotic result obtained by Ehrenmark [27] is illustrated in Figure 4 as μ 0 10 3 , only as an indication to show that the domains are complementary. While both solutions recover the Burnside expression for | h | = 0 , no quantitative comparison is possible for small slopes since the solution of Ehrenmark is valid for infinitesimally small values of μ 0 , where the present theory has a singularity (except for flat beds).
For practical purposes, and although the exact expressions for F 1 and F 2 are provided in Appendix B, the simpler expressions
F 1 = exp ( 0.315 μ 0 2 0.039 μ 0 3 + 0.013 μ 0 4 ) 8 μ 0 2
F 2 = + exp ( + 2.430 μ 0 2 2.173 μ 0 3 + 0.641 μ 0 4 0.069 μ 0 5 ) 32 μ 0 2 ,
provide approximations with errors below 0.5 % for the range of interest, i.e., 0 μ 0 2.3 .

6.3. Numerical Solution for the General Case ( θ 0 / z 0 )

In the general case, when θ 0 / z is not null, i.e., recalling Equation (43), when h · s 0 0 , the system of Equations (53)–(55) has no analytical solutions for A 1 . Now, the Equation (53) with bottom boundary condition (54) is to be solved numerically to obtain A 1 , and k 1 is recovered from (55), also numerically. However, given the linear nature of the equations, it can be foreseen that the solution is still k = k 0 ( 1 + F ) , now with
F = F 1 h 2 h + F 2 | h | 2 + F 3 ( h · s 0 ) 2 ,
where functions F 1 and F 2 are the ones presented above, and F 3 , obtained numerically, is shown in Figure 3. In fact, in solving the system of equations numerically, F 1 and F 2 can also be estimated numerically, and they have been used to check the quality of the numerical integration and the overall validity of Expression (66). For the range of interest, 0 μ 0 2.3 , the relative errors of the numerical estimations of F 1 and F 2 , compared to the exact expressions, are below 10 6 , which also guarantees the quality of the numerical estimation of F 3 . The function F 3 , included in Figure 3, can be again approximated in the range of interest through an exponential expression as
F 3 = exp ( 0.396 μ 0 2 0.143 μ 0 3 0.047 μ 0 4 + 0.034 μ 0 5 ) 6 ,
with errors below 0.5 % relative to the numerical solution. Function F 3 is non-negligible, when compared to F 1 and F 2 , only for μ 0 0.3 (see Figure 3A), i.e., where all F i are small and the influence of the bottom variations is expected to be nearly negligible.

6.4. Comparison to the Expression by Ge et al. [29]

The expression by Ge et al. [29], for sinusoidal bathymetries, is obtained using Fredholm’s alternative theorem (FAT) and by coupling the governing equation with the wave number and the varying seabed effects, without explicitly using perturbative approaches (see the original paper for details). Their equation is
ω 2 g k = ( tanh k h + G 1 h 2 h + G 2 | h | 2 ) ,
where
G 1 ( μ ) = ξ 4 μ , G 2 ( μ ) = ξ 2 ( 9 σ 2 μ 1 μ + μ 2 μ 2 + μ 3 μ 3 + μ 4 μ 4 ) 12 ( σ + μ ξ ) 3 ,
with σ = tanh μ , ξ = 1 σ 2 , and
μ 1 = 12 σ + 9 σ 3 , μ 2 = 9 + 45 σ 2 18 σ 4 , μ 3 = 28 σ 78 σ 3 + 30 σ 5 , μ 4 = 10 40 σ 2 + 42 σ 4 12 σ 6 .
Equation (68) can be rewritten in the same form as our Equation (60). First, note that G 1 h 2 h and G 2 | h | 2 are the small perturbations due to the bed shape, 𝒪 ( χ 2 ) . By following the perturbation theory procedures, we can assume that k is k 0 + χ 2 k 1 + 𝒪 ( χ 4 ) . In doing so, the above Equation (68) yields, with μ 0 = k 0 h ,
k 0 = ω 2 g tanh k 0 h ,
i.e., the Burnside equation, and
χ 2 k 1 = k 0 G 1 ( μ 0 ) h 2 h + G 2 ( μ 0 ) | h | 2 μ 0 tanh 2 μ 0 tanh μ 0 μ 0 ,
so that we can write (compare to Equations (60) and (61))
k = k 0 ( 1 + F 1 , G h 2 h + F 2 , G | h | 2 ) + 𝒪 ( χ 4 ) ,
where
F i , G = G i ( μ 0 ) μ 0 tanh 2 μ 0 tanh μ 0 μ 0 , i = 1 , 2 .
Note that, from Figure 5, and for μ 0.3 , in spite of the different approaches considered, the results collapse for both i = 1 , 2 . However, the signs are opposite, i.e., F i , G F i .

6.5. Numerical Checking of the Solution

The results for the second-order solution involve a considerable amount of algebra that deserves some kind of fully independent checking, especially after the results in Section 6.4. Here, we will numerically check the case of θ 0 / z = 0 , i.e., when an analytical solution is available (Section 6.2). We consider this case because the analytical solution simplifies the analysis, and at the same time, it includes the core of the work ( F 1 and F 2 , which dominate F, are obtained in this case). Additionally, we will consider a one-dimensional (1D) case for convenience and without loss of generality.
In recalling Equation (60), the second-order solution is
k = k 0 + χ 2 k 1 = k 0 ( 1 + F ) = k 0 ( 1 + F 1 h h x x + F 2 h x 2 ) ,
where F 1 and F 2 have analytical expressions in Appendix B. From Equations (33) and (57) for A 0 and A 1 , the second-order solution for A is the following; in recalling Section 6.1, note that the evaluations are performed particularly for A 0 and K c 0 , using k and q = k ( z + h ) ,
A = A 0 + χ 2 A 1 = K c 0 cosh q + K c 0 δ s 1 4 q 2 δ s 0 2 q cosh q h h x x μ 2 + K c 0 2 δ c 0 + δ s 1 4 q + μ 2 sinh q h h x x μ 2 + K c 0 ϵ s 1 + ϵ c 2 4 q 2 + 2 ϵ s 0 + ϵ c 1 4 q cosh q h x 2 μ 2 + K c 0 ϵ c 2 6 q 3 ϵ c 1 4 q 2 + 2 ϵ c 0 + ϵ s 1 ϵ c 2 4 q α k 0 μ 2 sinh q h x 2 μ 2 + K c 0 K c 1 cosh q ,
where K c 0 is, from Equation (35),
K c 0 = a sinh 2 μ + 2 μ ,
where a is a constant representative of the deep-water wave amplitude.
The numerical checking is carried using the profile by Yu and Slinn [32], i.e.,
h = 3.5 tanh 6 x l 1.2 x l 2.5 exp 100 x l + 0.25 2 ,
where x < 0 . Figure 6 shows the profile for two different values of l; it is straightforward from the expression of the profile that changing l has the effect of stretching or shrinking in the horizontal axis direction, and it can be readily shown that, for constant x 0 / l (e.g., dashed lines in Figure 6), it is h x l 1 and h x x l 2 , so that l 1 plays the role of χ here.
The original conditions, prior to the use of WKB, the governing Equations (16)–(18) read as follows in 1D and ignoring θ / z :
( PDE ) A x x A k 2 + A z z = 0 , h z 0 ,
( BBC ) A x h x + A z = 0 ,       z = h ,
( FSBC ) ω 2 A + g A z = 0 , z = 0 ,
where, as in Figure 6 and Figure 7, PDE = partial differential equation, BBC = bottom boundary condition, and FSBC = free-surface boundary condition. The correctness of the results is checked by introducing the above solutions (71) and (72) into a finite difference discretization of the governing Equations (73)–(75). To ensure that errors due to the finite difference discretization do not spoil the analysis, Δ x = Δ z = 0.01 and very-high-order finite difference approximations are used, with truncation errors 𝒪 ( Δ x 8 , Δ z 8 ) .
According to the perturbative approach followed in this work, when introducing the obtained solutions in the original (raw) governing equations, there will be errors 𝒪 ( χ 4 ) = 𝒪 ( l 4 ) in all three equations (PDE, BBC, and FSBC), which correspond to unresolved orders—recall that θ / z = 0 . Figure 7A shows how the errors do actually behave as expected—the evaluations are always performed at x / l = 0.2 , as highlighted in Figure 6.
To illustrate how any change in the proposed solution can spoil it, making the convergence not be 𝒪 ( χ 4 ) , we consider two (wrong-on-purpose) cases. First, the Burnside solution k 0 is considered instead of k = k 0 ( 1 + F ) in all the equations involved; the result, in Figure 7B, shows that the convergence rate for both the PDE and the BBC is correct, but it is 𝒪 ( χ 2 ) for the FSBC, precisely the equation from which k is derived. This emphasizes the well-known relevance of the free-surface boundary in obtaining the dispersion relationship. Finally, as shown in Figure 7C, if the result is considered to be δ c 0 = 1.01 α K c 0 instead of δ c 0 = α K c 0 , as it should be according to Table 1, in this case, the convergence of both the PDE and FSBC is spoiled.

7. Two Cases of Application

Water wave propagation numerical models rely on depth-integrated equations such as the linear and fully dispersive MSEs (modified or not) or the nonlinear and weakly dispersive Boussinesq-type equations. These equations often neglect the effects of bottom variations or utilize velocity profile assumptions that differ from the solution presented in this study, and do not provide an expression for the corresponding dispersion equation. In using the new results presented here, particularly in Section 6 for 𝒪 ( χ 2 ) , the consequences of including (or neglecting) the bottom variations on the wave celerity are shown in this Section 7 with two one-dimensional (1D) examples—although the results provided here are valid for 2D cases. Both examples correspond to shallow-water conditions since, as already seen, the influence of the bed shape is negligible in deep waters.

7.1. Nearshore Propagation and Depth Inversion Problem

A first example deals with the propagation of wind-generated water waves as they approach the shore over a barred beach profile. We consider the propagation of linear waves with a period T = 7 s over the barred beach profile shown in Figure 8D, given by
h = 3.5 tanh 6 x l 1.2 x l 2.5 exp 100 x l + 0.25 2 ,
where x < 0 and l = 100 m . Note that this profile has relatively strong slopes to enhance the influence of the bottom variations; a more realistic case would consider l > 200 m .
The parameter μ 0 = k 0 h , which is of paramount importance in the present analysis, is obtained from the leading-order dispersion equation and included in Figure 8A; it ranges from ∼0.6 in the offshore deepest zone to ∼0.3 over the bar and ∼0.2 in the shallower part. The factor F, which is required to obtain k = k 0 ( 1 + F ) , reads in this 1D case
F = F 1 h h x x + F 2 h x 2 + F 3 h x 2 ,
with h x = d h / d x and h x x = d 2 h / d x 2 . Figure 8B shows the evolution of the three terms of the sum: the term F 1 h h x x dominates over the bar, while F 2 h x 2 dominates very near the shore. As expected, the term F 3 h x 2 is generally small. The wave celerities c 0 = ω / k 0 and c = ω / k are shown in Figure 8C, with the differences reaching values up to ± 7 % .
The above differences, due to F, have consequences on a relatively recent branch of algorithms for bathymetry estimation from video images. These algorithms, e.g., refs. [33,34,35] obtain the bathymetry from videos of the wave propagation in the shallow-water region, i.e., where the waves feel the bathymetry. For this purpose, these algorithms first obtain estimates of the wave period and the wavenumber using different techniques such as Fourier analysis or mode decomposition, and once ω and k have been estimated, all of them rely on the leading-order dispersion relationship to retrieve the water depth h.
Figure 9 is meant to illustrate the consequences of using the leading-order dispersion equation in bathymetric inversion methods, by simulating the type of errors that can be incurred when disregarding the slope h x and curvature h x x of the bottom. Figure 9 shows a bathymetry h, the same as in Figure 8D, the bathymetry retrieved using the leading-order dispersion equation (Burnside), and the wavenumber values obtained through the modified equation (i.e., the ones that would actually be observed). Linear waves of T = 7 s have again been considered. The maximum errors, i.e., the differences between the proposed bathymetry and the retried one are ∼20 cm in this case. If the characteristic length of the bathymetry, l, is increased to 200 m in Equation (76), i.e., in milder conditions, these errors are reduced to ∼5 cm. It should be noted that the bathymetric inversion algorithms have other difficulties related to the nonlinearity of the waves or the obliquity of the images that can be more relevant.

7.2. Tsunami Propagation over the Continental Slope

A second example is used to specifically show the influence of the slope, h x , and also to show the limitations of the present analysis in very shallow waters. We consider the propagation of a tsunami over a continental slope as it travels to the shore. The continental slope is the abrupt rise from the ocean floor to the edge of the continental shelf. It marks the boundary between the continental shelf and the abyssal plain, and the slope angles typically range from 3 to 6 degrees (i.e., 0.05 | h x | 0.10 ), reaching higher values in subduction zones. The idealized continental slope considered in this example has | h x | = 0.075 , taking ∼25 km to go from a depth of h = 2000 m to 100 m (see Figure 10D). The junctions between the slope and flat zones (the abyssal plain and the continental shelf, respectively) have been smoothed over a length l = 4 km using polynomials that ensure the continuity of h, h x , and h x x . We initially consider the propagation of a tsunami with a wavelength L = 200 km in the deep zone; this would correspond, using the dispersion relationship, to a period of T 24 min and a celerity, in the deep zone, of ∼500 km/h.
Figure 10 shows the most relevant results as in Figure 8. From Figure 10A, μ 0 0.07 in all the domain, with μ 0 0.01 in the shallower zone (continental shelf). Here, the values of μ 0 are an order of magnitude smaller than in the previous example, which implies a more limiting condition on the maximum slope that can be handled. Figure 10B shows again the influence of each of the three terms summing up factor F in Equation (77). Owing to the way in which the bathymetry has been defined, h x x is null over the whole slope (except in the junction areas), and the term F 1 h h x x is zero. Further, F 3 F 2 in the slope, as expected due to the small values of μ 0 (Figure 3). In summary, the term F 2 h x 2 dominates the slope, with increasing values as the depth decreases and reaching values of ∼0.4. The impact on the wave celerity is shown in Figure 10C; the traveling time of the wave over the part of the slope where h x x = 0 is 231 s using c 0 and 262 s using the celerity c (i.e., + 13.1 % ).
Table 2 shows results obtained for the traveling time over the zone of the slope where h x x = 0 for different combinations of slope | h x | and wavelength L in the deep zone. The results show, quantitatively, that the greater the slope steepness and the longer the wavelength, the larger the differences in the traveling time become. At the same time, as the slope or the wavelength L increases, the maximum value of | F | at that zone, | F | max , increases, and the condition | F | max < F c is violated. This is indicated with symbols in Table 2 (the results marked with † and ‡ should be regarded with caution). Tsunamis do usually have values of μ 0 so small that even relatively small slopes of the continental shelf can violate the slope condition, both for the original condition in Equation (2), and also the new one in Equation (63). This point highlights the potential errors in tsunami propagation over the continental slope using depth-integrated models.
A second relevant aspect of Figure 10 is that the term F 1 h h x x is dominant in the junction zones and, moreover, forces the violation of the condition | F | < F c . Given that the condition is not satisfied, the results of the current theory cannot be applied to the junctions, and the behavior of the celerity remains unknown in that zone.

8. Concluding Remarks

This study explored the impact of the spatial variability in the seabed on key wave properties relevant to linear wave propagation employing a WKB approximation. Through solving the equations up to second-order 𝒪 ( χ 2 ) , where χ represents a characteristic value of the slope, two significant outcomes emerged.
Firstly, the well-established longwave theory governing linear transmission and reflection over a ramp was expanded to intermediate and deep waters. Secondly, and more prominently, a new dispersion relationship that incorporates the gradient and Laplacian of the mean water depth was derived. This resulting expression is readily applicable, and its implications for the problems of depth inversion and tsunami propagation are illustrated with two examples.
For the specific case of sloped planes, the analysis broadens the applicability range of the dispersion relationship over uneven bottoms from the classical mild-slope condition ( | h | 0.38 k h ) to an extended range ( | h | 2.5 k h ). Importantly, our findings underscore potential pitfalls associated with employing depth-integrated models under specific circumstances. The challenge posed by the tsunami propagation over the continental slope is a case that may surpass the validity of this study, clearly showing the influence of the Laplacian term.

Funding

This research was funded by the Spanish Ministry of Science, Innovation and Universities—National Research Agency and EU “NextGenerationEU/PRTR” grant numbers PID2021-124272OB-C21/C22 (MOLLY) and TED2021-130321B-I00 (SOLDEMOR).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The author expresses gratitude for the valuable insights gained from discussions with Daniel Calvete, Asunción Baquerizo, Juan Simarro, and Jorge Guillén. This work is dedicated to Irene, Elena, and Nara.

Conflicts of Interest

The author declares no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BTEsBoussinesq-type equations;
MSEsMild-slope equation;
SWEsShallow-water equations;
WKBWentzel, Kramers, and Brillouin (approach);
PDEPartial differential equation;
BBCBottom boundary condition;
FSBCFree-surface boundary condition.

Appendix A. Expressions for α and β

In this appendix, we avoid, for ease, the subindex “0” to denote that we are dealing with the leading-order solution. From dispersion relationship (34), deriving relative to h, we have
0 = g d k d h tanh μ + g k ( 1 tanh 2 μ ) d k d h h + k ,
so that after manipulation,
α k = h k d k d h = 2 μ sinh 2 μ + 2 μ = G 1 + G , G = 2 μ sinh 2 μ .
Since c = ω / k and μ = k h , it can be easily shown that α c = α k and α μ = α k 1 . It can also be shown that α G = ( 1 G cosh 2 μ ) α μ . From energy Equation (35), we have the following manipulations:
α K c = G M 2 ( 1 + G ) 2 , [ M = cosh 2 μ + 1 ] ,
Regarding β , related to second-order derivatives, it is, e.g.,
β k = G 2 ( M + 2 ( 1 + G ) ) ( 1 + G ) 3 ,
β K c = 2 μ 2 ( 1 + G ) 3 + G 2 M ( 2 G M 4 ( 1 + G ) 3 M ) 4 ( 1 + G ) 4 .

Appendix B. Functions F1 and F2 for k1

The expressions for both F 1 and F 2 are ( j = 1 , 2 )
F j = 1 12 μ 0 2 F j a μ 0 tanh μ 0 sinh μ 0 F j b sinh μ 0 F j c μ 0 cosh μ 0 μ 0 tanh μ 0 sinh μ 0 sinh μ 0 μ 0 cosh μ 0 ,
where the coefficients F j a , F j b , and F j c are
F 1 a = 3 ( α + 2 ) , F 1 b = 3 α , F 1 c = 3 α + .
and
F 2 a = + 2 ( α k 0 2 3 α k 0 + 3 ) μ 2 3 ( β + α α k 0 2 α k 0 ) , F 2 b = + 6 ( α k 0 2 2 α k 0 + 1 ) μ 2 3 ( β + α α k 0 ) , F 2 c = + 2 ( α k 0 2 3 α k 0 + 3 ) μ 2 3 ( β + α α k 0 2 α k 0 + 2 α + ) ,
where α ± = 2 α K c 0 ± α k 0 and β ± = 2 β K c 0 ± β k 0 .

References

  1. Mei, C.C. The Applied Dynamics of Ocean Surface Waves; World Scientific: Singapore, 1989. [Google Scholar]
  2. Dean, R.G.; Dalrymple, R.A. Water Wave Mechanics for Engineers and Scientists; World Scientific: Singapore, 1991. [Google Scholar]
  3. Svendsen, I.A. Introduction to Nearshore Hydrodynamics; Advanced Series on Ocean Engineering; World Scientific: Singapore, 2005. [Google Scholar]
  4. Kim, D.; Powers, T. Vertical integration of the Navier-Stokes equations for wave propagation in coastal waters. J. Fluid Mech. 2006, 566, 425–455. [Google Scholar]
  5. Chen, Q.; Kirby, J.T.; Dalrymple, R.A.; Kennedy, A.B.; Chawla, A. Boussinesq modeling of wave transformation, breaking, and runup. II: 2D. J. Waterw. Port Coast. Ocean Eng. 2000, 126, 48–56. [Google Scholar] [CrossRef]
  6. Mendez, F.J.; Losada, I.J. An empirical model to estimate the propagation of random breaking and nonbreaking waves over vegetation fields. Coast. Eng. 2004, 51, 103–118. [Google Scholar] [CrossRef]
  7. Infantes, E.; Orfila, A.; Simarro, G.; Terrados, J.; Alomar, C.; Tintoré, J. Effect of a seagrass (Posidonia oceanica) meadow on wave propagation. Mar. Ecol. Prog. Ser. 2012, 456, 63–72. [Google Scholar] [CrossRef]
  8. Mentaschi, L.; Besio, G.; Cassola, F.; Mazzino, A. Performance evaluation of Wavewatch III in the Mediterranean Sea. Ocean Model. 2015, 90, 82–94. [Google Scholar] [CrossRef]
  9. Airy, G.B. Tides and waves. Encycl. Metrop. 1845, 5, 241–396. [Google Scholar]
  10. Burnside, W. On the modification of a train of waves as it advances into shallow water. Proc. Lond. Math. Soc. 1915, s2_14, 131–133. [Google Scholar] [CrossRef]
  11. Stokes, G.G. On the theory of oscillatory waves. Trans. Camb. Philos. Soc. 1847, 8, 441–455. [Google Scholar]
  12. Berkhoff, J.C.W. Computation of combined refraction-diffraction. In Proceedings of the 13th Conference on Coastal Engineering, Vancouver, BC, Canada, 10–14 June 1972; Volume 2, pp. 471–490. [Google Scholar]
  13. Berkhoff, J.C.W. Mathematical Models for Simple Harmonic Linear Waves: Wave Diffraction and Refraction. Ph.D. Thesis, TU Delft, Delft, The Netherlands, 1976. [Google Scholar]
  14. Chamberlain, P.G.; Porter, D. The modified mild-slope equation. J. Fluid Mech. 1995, 291, 393–407. [Google Scholar] [CrossRef]
  15. Hsu, T.W.; Lin, T.Y.; Wen, C.C.; Ou, S.H. A complementary mild-slope equation derived using higher-order depth function for waves obliquely propagating on sloping bottom. Phys. Fluids 2006, 18, 087106. [Google Scholar] [CrossRef]
  16. Porter, R. An extended linear shallow-water equation. J. Fluid Mech. 2019, 876, 413–427. [Google Scholar] [CrossRef]
  17. Porter, D. The mild-slope equations: A unified theory. J. Fluid Mech. 2020, 887, A29. [Google Scholar] [CrossRef]
  18. Dingemans, M.W. Water Wave Propagation over Uneven Bottoms; World Scientific: Singapore, 1997. [Google Scholar] [CrossRef]
  19. Griffiths, D.J. Introduction to Fluid Mechanics, 1st ed.; Prentice Hall: Englewood Cliffs, NJ, USA, 1976. [Google Scholar]
  20. Boussinesq, J. Théorie des vagues propagées le long d’un canal ou d’une mer peu profonde. Comptes Rendus L’AcadéMie Des Sci. 1872, 74, 891–897. [Google Scholar]
  21. Peregrine, D.H. Long waves on a beach. J. Fluid Mech. 1967, 27, 815–827. [Google Scholar] [CrossRef]
  22. Liu, P.L.F. Model equations for wave propagation from deep water to shallow water. In Advances in Coastal and Ocean Engineering; World Scientific: Singapore, 1994; Chapter 3; pp. 125–158. [Google Scholar]
  23. Wei, G.; Kirby, J.T.; Grilli, S.T.; Subramanya, R. A fully nonlinear Boussinesq model for surface waves. Part 1: Highly nonlinear unsteady waves. J. Fluid Mech. 1995, 294, 71–92. [Google Scholar] [CrossRef]
  24. Simarro, G.; Orfila, A.; Galan, A. Linear shoaling in Boussinesq-type wave propagation models. Coast. Eng. J. 2013, 80, 100–106. [Google Scholar] [CrossRef]
  25. Nwogu, O. Alternative form of Boussinesq equations for nearshore wave propagation. J. Waterw. Port Coastal Ocean Eng. 1993, 119, 618–638. [Google Scholar] [CrossRef]
  26. Booij, N. A note on the accuracy of the mild-slope equation. Coast. Eng. 1983, 7, 191–203. [Google Scholar] [CrossRef]
  27. Ehrenmark, U. An alternative dispersion equation for water waves over an inclined bed. J. Fluid Mech. 2005, 543. [Google Scholar] [CrossRef]
  28. Williams, P.; Ehrenmark, U. A note on the use of a new dispersion formula for wave transformation over Roseau’s curved beach profile. Wave Motion 2010, 47, 641–647. [Google Scholar] [CrossRef]
  29. Ge, H.; Liu, H.; Zhang, L. Accurate depth inversion method for coastal bathymetry: Introduction of water wave high-order dispersion relation. J. Mar. Sci. Eng. 2020, 8, 153. [Google Scholar] [CrossRef]
  30. Dingemans, M.W. Water Wave Propagation over Uneven Bottoms. Ph.D. Thesis, Delft University of Technology, Delft, The Netherlands, 1994. [Google Scholar]
  31. Schäffer, H.A.; Madsen, P.A. Further enhancements of Boussinesq-type equations. Coast. Eng. 1995, 26, 1–14. [Google Scholar] [CrossRef]
  32. Yu, J.; Slinn, D. Effects of wave-current interaction on rip currents. J. Geophys. Res 2003, 108, 87. [Google Scholar] [CrossRef]
  33. Holman, R.; Bergsma, E.W.J. Updates to and performance of the cBathy algorithm for estimating nearshore bathymetry from remote sensing imagery. Remote Sens. 2021, 13, 3996. [Google Scholar] [CrossRef]
  34. Gawehn, M.; de Vries, S.; Aarninkhof, S. A self-adaptive method for mapping coastal bathymetry on-the-fly from wave field video. Remote Sens. 2021, 13, 4742. [Google Scholar] [CrossRef]
  35. Simarro, G.; Calvete, D. UBathy (v2.0): A software to obtain the bathymetry from video imagery. Remote. Sens. 2022, 14, 6139. [Google Scholar] [CrossRef]
Figure 1. Coordinates x and z (one-dimensional case), and elevations of the mean sea level ( z = 0 ), free surface ( z = η ), and bed ( z = h ).
Figure 1. Coordinates x and z (one-dimensional case), and elevations of the mean sea level ( z = 0 ), free surface ( z = η ), and bed ( z = h ).
Jmse 12 01652 g001
Figure 2. Function r = r ( μ 0 ) for the transmission and reflection on a ramp.
Figure 2. Function r = r ( μ 0 ) for the transmission and reflection on a ramp.
Jmse 12 01652 g002
Figure 3. Functions F i = F i ( μ 0 ) for i = 1 , 2 , 3 using log-log (A) and semi-log (B) scales. Dashed line is at μ 0 = 2.3 .
Figure 3. Functions F i = F i ( μ 0 ) for i = 1 , 2 , 3 using log-log (A) and semi-log (B) scales. Dashed line is at μ 0 = 2.3 .
Jmse 12 01652 g003
Figure 4. Sloped plane: isolines 0.2 and 0.001 for F = F 2 | h | 2 (sloped plane).
Figure 4. Sloped plane: isolines 0.2 and 0.001 for F = F 2 | h | 2 (sloped plane).
Jmse 12 01652 g004
Figure 5. Functions F i = F i ( μ 0 ) and F i , G = F i , G ( μ 0 ) for i = 1 , 2 using log-log (A) and semi-log (B) scales.
Figure 5. Functions F i = F i ( μ 0 ) and F i , G = F i , G ( μ 0 ) for i = 1 , 2 using log-log (A) and semi-log (B) scales.
Jmse 12 01652 g005
Figure 6. Yu and Slinn [32] profile for numerical checkings. FSBS: point to check the free-surface boundary condition; PDE: point to check the partial differential equation; BBC: point to check the bottom boundary condition.
Figure 6. Yu and Slinn [32] profile for numerical checkings. FSBS: point to check the free-surface boundary condition; PDE: point to check the partial differential equation; BBC: point to check the bottom boundary condition.
Jmse 12 01652 g006
Figure 7. Evolution of the errors for equations using the proposed solutions. (A) Proposed solution (thin lines have slope 4 , corresponding to 𝒪 ( l 4 ) ); (B) incorrectly using Burnside equation for k (thin lines have slopes 2 and 4 ); (C) using the proper k but introducing an error on purpose— δ c 0 = 1.01 α K c 0 instead of δ c 0 = α K c 0 as it should according to Table 1 (thin lines have slopes 2 and 4 ).
Figure 7. Evolution of the errors for equations using the proposed solutions. (A) Proposed solution (thin lines have slope 4 , corresponding to 𝒪 ( l 4 ) ); (B) incorrectly using Burnside equation for k (thin lines have slopes 2 and 4 ); (C) using the proper k but introducing an error on purpose— δ c 0 = 1.01 α K c 0 instead of δ c 0 = α K c 0 as it should according to Table 1 (thin lines have slopes 2 and 4 ).
Jmse 12 01652 g007
Figure 8. Wind-generated water waves of period T = 7 s propagation over a barred beach profile: evolution of μ 0 (A), components of F (B), celerities c 0 and c (C) and water depth h (D).
Figure 8. Wind-generated water waves of period T = 7 s propagation over a barred beach profile: evolution of μ 0 (A), components of F (B), celerities c 0 and c (C) and water depth h (D).
Jmse 12 01652 g008
Figure 9. The inversion problem in a barred beach profile (h retrieved with the Burnside equation).
Figure 9. The inversion problem in a barred beach profile (h retrieved with the Burnside equation).
Jmse 12 01652 g009
Figure 10. Tsunami propagation over a continental shelf: evolution of μ 0 (A), components of F (B), celerities c 0 and c (C) and water depth h (D).
Figure 10. Tsunami propagation over a continental shelf: evolution of μ 0 (A), components of F (B), celerities c 0 and c (C) and water depth h (D).
Jmse 12 01652 g010
Table 1. Expressions for δ and ϵ —see Appendix A for α and β .
Table 1. Expressions for δ and ϵ —see Appendix A for α and β .
δ s 1 = α k 0
δ c 0 = α K c 0 δ s 0 = μ 0
ϵ c 2 = α k 0 2
ϵ c 1 = 2 μ α k 0 ϵ s 1 = 2 α K c 0 α k 0 β k 0
ϵ c 0 = μ 0 2 β K c 0 ϵ s 0 = 2 μ 0 ( α K c 0 + α k 0 )
Table 2. Increment in traveling time over the slope (the zone where h x x = 0 ) using the modified dispersion relationship instead of using the leading-order relationship. The cases with 0.2 < | F | max 0.4 are marked with † ( | F | max = maximum value of | F | in the ramp), those with 0.4 < | F | max 0.6 with ‡, and those with | F | max > 0.6 are ignored (—).
Table 2. Increment in traveling time over the slope (the zone where h x x = 0 ) using the modified dispersion relationship instead of using the leading-order relationship. The cases with 0.2 < | F | max 0.4 are marked with † ( | F | max = maximum value of | F | in the ramp), those with 0.4 < | F | max 0.6 with ‡, and those with | F | max > 0.6 are ignored (—).
L in the Deep Zone [km]
Slope | h x | [-] 1050100200500
0.025 + 0.0 % + 0.1 % + 0.5 % + 1.8 % + 11.4 %
0.050 + 0.0 % + 0.4 % + 1.6 % + 6.4 %
0.075 + 0.0 % + 0.8 % + 3.2 % + 13.1 %
0.100 + 0.0 % + 1.3 % + 5.3 % + 21.5 %
0.150 + 0.1 % + 2.5 % + 10.6 %
0.200 + 0.0 % + 4.1 % + 17.4 %
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Simarro, G. Influence of Bed Variations on Linear Wave Propagation beyond the Mild Slope Condition. J. Mar. Sci. Eng. 2024, 12, 1652. https://doi.org/10.3390/jmse12091652

AMA Style

Simarro G. Influence of Bed Variations on Linear Wave Propagation beyond the Mild Slope Condition. Journal of Marine Science and Engineering. 2024; 12(9):1652. https://doi.org/10.3390/jmse12091652

Chicago/Turabian Style

Simarro, Gonzalo. 2024. "Influence of Bed Variations on Linear Wave Propagation beyond the Mild Slope Condition" Journal of Marine Science and Engineering 12, no. 9: 1652. https://doi.org/10.3390/jmse12091652

APA Style

Simarro, G. (2024). Influence of Bed Variations on Linear Wave Propagation beyond the Mild Slope Condition. Journal of Marine Science and Engineering, 12(9), 1652. https://doi.org/10.3390/jmse12091652

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop