The Three-Dimensional Wave-Induced Current Field: An Analytical Model

: Wave-induced currents play a critical role in coastal dynamics, influencing sediment transport and shaping bottom topography. Traditionally, long-and cross-shore currents in coastal zones were analyzed independently, often with two-dimensional models for longshore currents and undertow being used. The introduction of quasi-three-dimensional models marked a significant advancement toward a more holistic understanding. Despite recent proposals for fully three-dimensional models, none have achieved widespread acceptance, primarily due to challenges in accurately capturing depth-dependent radiation stress. This paper presents an innovative analytical model advocating for comprehensive three-dimensional approaches in coastal hydrodynamics. The model, based on novel simplification rules, refines relationships governing turbulent stress tensors and provides valuable insights into wave-induced stresses. It offers analytical solutions for both homogeneous and general coastal zones, laying the foundation for future advancements in numerical modeling techniques.


Introduction
Wave-induced currents play an important role in the dynamics of coastal zones dominated by wave activity, serving as a crucial mechanism for transferring energy from the atmosphere to the ocean and eventually to the seabed.Among these currents, the undertow and longshore currents stand out due to their significant contributions to the movement of bedload and, alongside wave-generated orbital velocities and turbulent movements near the seabed, to changes in the coastal zone bottom topography.Consequently, understanding these currents is essential for predicting erosional and accumulative processes in coastal zones, as well as for addressing key issues related to coastal protection and ensuring the safety of ships near ports.To gain a complete understanding of wave-induced currents, a full three-dimensional study is required, as this is the only method capable of accounting for all components influencing the velocity field within the coastal zone.Figure 1 illustrates the three-dimensional mean velocity distribution within the coastal zone.
Historically, the main currents within coastal zones, including both long-and cross-shore currents, have been analyzed independently.Since the early seventies, two-dimensional models for these currents have been established, with longshore currents modeled based on depth-integrated radiation stress and the undertow driven by the imbalance between the average depth-dependent momentum flux induced by surface waves and the slope of an average position of the water surface, described by the radiation stress [1][2][3].
By the late 1980s, the requirement for a three-dimensional representation of wave-induced current fields had been identified, prompting the creation of quasi-three-dimensional models [4].These models attempted to combine separate models of longshore currents and undertow to offer an approximate spatial distribution of the flow (see, e.g., [5,6]).Despite the popularity of these models in the 1990s, their reliance on two-dimensional solutions meant that in certain scenarios, such as within the coastal zone of a small island, the resulting velocity field could significantly deviate from reality.A comprehensive understanding of wave-induced currents in coastal zones necessitates a three-dimensional approach, ensuring the integration of all influencing elements into the flow dynamics.However, the primary challenge in developing such models lies in identifying depth-dependent radiation stress, which serves as the driving force behind these currents.Several attempts have been made to address this challenge.For instance, Nobuoka et al. [7] introduced a novel three-dimensional nearshore current model, incorporating key concepts such as the variability of onshore and offshore flow and a sole reliance on wave motion as the driving force.Mellor's ( [8][9][10]) research contributed to advancing comprehensive models integrating three-dimensional current dynamics and surface wave interactions, offering a holistic understanding of ocean dynamics.
Additionally, Kumar et al. [11] aimed to enhance three-dimensional radiation stress formulations to improve the accuracy of predicting wave-induced currents and rip currents in the surf zone.Chou et al. [12] focused on developing and applying a three-dimensional numerical model coupling wave dynamics with hydrodynamics in South San Francisco Bay, aiming to understand the complex interactions in the Bay Area.Ranasinghe [13] concentrated on developing advanced numerical models to simulate interactions between waves and submerged structures under regular waves, ultimately aiming to enhance the understanding of wave interactions in such scenarios.
It is worth noting that machine learning is rapidly advancing and being applied across various disciplines, including coastal zone hydrodynamics (see, e.g., [14]).However, it faces several limitations that currently hinder the successful description and resolution of wave-induced currents.
In terms of analytical model formulation, Xia et al. [15] focused on the vertical variation of radiation stress, providing insights into its depth-dependent changes and its influence on the generation and behavior of wave-induced currents.Additionally, an alternative formulation based on the Dirac delta function, presented by Kołodko and Gic-Grusza [16], offers a simplified and concise description of the vertical variability of excess horizontal momentum flux due to wave presence, presenting a novel approach to the analytical modeling of wave-induced currents.Chybicki [17] also contributed significantly to this discourse by extending the original concept of radiation stress in water waves to encompass its vertical variation, achieving a precise delineation of radiation stress density alongside straightforward computational formulas for assessing the vertical radiation stress profile.This work parallels the foundational principles posited by Dolata and Rosenthal [18], albeit with a more refined expression for the pressure function.The discourse on radiation stress theory was further elaborated by Mellor [19], who engaged in a comparative analysis between the concepts of "radiation stress" and "vortex force."In this work, Mellor advocates for the Longuet-Higgins ( [1,2]) and Mellor ([8,20]) theoretical framework, critiquing the "vortex force" theory as extended by McWilliams and Restrepo [21] and addresses the inaccuracies present in research based on the vortex force concept.The significance of radiation stresses within coastal zones was reaffirmed in Mellor's subsequent work [22], reinforcing their critical role in coastal dynamics.
The foundation of the research presented in this paper is anchored by numerous hypotheses, primarily focusing on the methodology for constructing the model's underlying equations.Central to contemporary hydrodynamics are the momentum transport equations (the Navier-Stokes equations), accompanied by the mass transport Equation (the continuity equation).However, it has become increasingly clear that transforming these equations into a structure that facilitates the creation of an exhaustive three-dimensional model for wave-induced currents in coastal areas poses substantial challenges and is yet to be fully accomplished.
Traditionally, the description of the flows in question utilized the wave periodaveraged mass transport equation (see, for example, Svendsen and Lorenz [23]; Nobuoka and Mimura [24]).However, some authors did not adequately consider the necessity of integrating this equation before its averaging (see, for example, Lee and Wang [25]; Smith [26]), resulting in a misinterpretation of the equations and, consequently, as observed by Cieślikiewicz and Gudmestad [27], the omission of the wave transport on the averaged free surface, i.e., the area between the wave crest and trough.The work presented here assigns crucial importance to this area, both in terms of mass and momentum transport.It seems that without accounting for the processes occurring on the averaged free surface, it is not feasible to meaningfully define the quasi-turbulent stress tensor (or, alternatively, quasi-Reynolds stress tensor), which acts as a local force-generating flow.
The main assumptions and hypotheses of this work are as follows: 1.
A consistent assumption that all terms of higher order than the square of the wave amplitude are negligibly small.This leads, among other things, to the exclusion from the model of wave-current and current-current interactions (unlike, for example, the work of Svendsen and Putrevu [28], in which certain terms of the order of the square of the wave amplitude, a, are omitted, while others-of the order of a 4 -are retained).

2.
The assumption is made that for the problem at hand, it is sufficient to apply the theory of small amplitude waves to describe wave motion.This is a widely accepted practice based on which, for instance, the components of the radiation stress tensor are defined.However, the applicability of this assumption is not suitable for the surfzone area.A further consequence of adopting linear wave theory in this work is the omission of the average value of products of orthogonal components of wave velocities (i.e., the product of horizontal velocity components with the vertical component), whose role is emphasized in the literature (see, for example, Rivero and Arcilla [29]).
For the same reason, the so-called "roller", defined as a mass of water moving between the crest and trough of a breaking wave at phase speed (see, e.g., [30]), which attempts to account for wave nonlinearity in the coastal zone, is also omitted.However, it is worth noting that Svendsen and Lorenz [23] estimated the contribution of the roller to the obtained flow velocity values at only 9%, somewhat justifying its exclusion in the prototype version of the constructed three-dimensional model.

3.
A new, simplified constitutive relationship defining turbulent stresses specifically for the coastal zone is proposed.To date, turbulence in this area has been poorly understood (see, for example, [31]).Therefore, the most straightforward approach to modeling the considered stresses using turbulent viscosity coefficients (according to the Boussinesq hypothesis) was chosen.However, taking into account the specifics of the coastal zone, certain additional simplifications were made, which consequently led to a significant "thinning" of the turbulent stress tensor, positively affecting the simplicity of the derived equations.

4.
A theory of distributions to describe momentum transport is adopted, allowing for a completely innovative formulation of the quasi-turbulent stress tensor, which, after integration into the vertical direction (with consideration to the part located on the averaged free surface of the distribution), yields the traditional (classical) radiation stress tensor.It should be noted that the proposed form of the tensor significantly deviates from the concept of vertical variability of radiation stresses published by Xia et al. [15] (as it seems incorrect).The authors of the discussed work omitted the stresses on the averaged free surface in their considerations and obtained vertical variability through the transformation of the traditional definition of radiation stresses proposed by Longuet-Higgins and Stewart [32] into a system of coordinates and through the calculation, in a somewhat convoluted manner, of the values of the integral expression.As for the radiation stresses on the averaged free surface, the results obtained in this work constitute a far-reaching generalization of the pioneering results of Stive and Wind ([33,34]) concerning the undertow.
It should be noted that without recourse to the theory of distributions, both in relation to the description of mass transport (the well-known Eulerian description of the Stokes drift) and momentum transport (the new, depth-distributed field of the so-called quasi-Reynolds radiation stresses and the partially known effective radiation friction on the averaged free surface), the presented model could not have been created.It should also be observed that only because purely distributive elements occur solely at the upper boundary (on the free surface) was it possible to use traditional partial differential equations, with only the conditions at the upper boundary being modified (distinguishing between conditions above and below the surface).Otherwise, integral equations would have had to be used to describe the considered problem.

Main Assumptions
The entire analysis is carried out in a rectangular, three-dimensional, right-handed coordinate system (x; y; z) (cf. Figure 1).All symbols used in this analysis are defined in Appendix A.
The model assumes that waves are described by the Airy theory and that no viscosity effects or capillary forces are taken into account.The water density is assumed to be constant (ρ = const).
Once turbulence is filtered out, there are two time scales: the short-term scale, associated with wave oscillations (and wave period), and the long-term one, in which wave characteristics may change.It seems thus reasonable to split the instantaneous velocity horizontally, as in U = (U; V), and in three dimensions, as in U = (U; V; W), with the free surface elevation (Z) and pressure (P) expressed as: where: Z = ⟨Z⟩, P = P + p, where: P = ⟨P⟩. ( where ⟨. ..⟩ denotes short-time averaging.This is a frequently used practice, although some authors express the instantaneous flow velocity as a sum of three components: the average velocity, wave fluctuations, and turbulent pulsations u ′ = (u ′ , v ′ , w ′ ) (see, e.g., [4,35]).Such an approach, although equally appropriate, seems hardly effective on account of a poor knowledge of the turbulence in the coastal zone.Some authors even point out that turbulence in the coastal zone is generated primarily by breaking waves, and thus the separation of certain turbulent components from the wave ones is virtually impossible ( [28,36]).In this work, turbulence is incorporated by using the Reynolds equations; turbulence parameterization was attained with coefficients of turbulent viscosity that are different in the horizontal and in the vertical direction.
In developing the equations, it was primarily assumed that all resultant quantities on the order higher than a 2 (where a is the wave amplitude) should be disregard; in this context is is worth noting the following: and that u ≫ U when a is small.

Mass Transport Equation
The mass transport equation for the current field is based on the continuity equation for constant density: The kinematic boundary conditions at the free surface (z = Z) and at the bottom (z = −D) are given as follows: Taking into account the decomposition of velocity (cf.Equation ( 1)), the continuity Equation (2) and the kinematic boundary conditions (3) can be expressed as follows: The averaging of the above equations on the time scale equal to the wave period leads to the continuity equation and the kinematic conditions for the current velocity in the area limited by the averaged free surface and the bottom: It should be noted that condition (7) at z = Z can be taken at z = 0 because Z is on the order a 2 .Similar equations can be obtained for wave oscillations by subtracting Equations ( 4)-( 7): At the next stage of model development, the volumetric transport Equation ( 2) is integrated with respect to the vertical coordinate z from the instantaneous free surface to the bottom and then averaged as follows: Assume the following: • The upper limit of the first integral can be replaced by 0, provided that integration is approaching zero from below (denoted by 0 − )' • The integral from Z to Z = Z + ζ practically equals the integral from 0 to ζ/ Once the terms on the order higher than a 2 are rejected, the equation can be written as follows: The volumetric transport equation can be written in another form using Leibniz's rule: where Q is the depth-integrated current-generated volumetric transport and q is the wave-generated volumetric transport In the Eulerian system, the depth-integrated wave-generated volumetric transport q is concentrated between the wave crest and trough, which, when averaged, corresponds to the mean surface level Z.As Z is on order a 2 , it may be assumed that the mass transport is located on z = 0 and that the horizontal orbital velocity at the mean surface level is equal to the horizontal orbital velocity on z = 0. Therefore, we obtain the following: where k is the wave vector, c is the phase velocity, and E is the wave energy.
The volumetric transport Equation ( 12) can be rewritten as follows: where Q is the depth-integrated total wave and current-generated volumetric transport: The physical quantity U is a generalized average velocity field: which is expressed by the Dirac delta function or distribution as follows: x ̸ = 0, and The use of the Dirac delta in the volumetric transport equation reveals that the average velocity field is a continuous function everywhere outside the averaged free surface where its value tends to infinity.This means that U is a form of distribution (a generalized function).
Given the degree of accuracy assumed (to the terms on the order a 2 ), the upper limit Z of the integrals ( 17) and ( 13) can be replaced by 0. A similar substitution can be applied to the Dirac delta function argument.However, as the Dirac delta coincides with the integration limit, the integration shall be approaching zero from above, which is denoted by 0 + as follows: Finally, assuming the constant density, the mass transport equation is the following: where M is the total averaged mass transport, M = ρQ.

Momentum Transport Equation
After filtration of turbulent pulsations denoted as ≪≫, the momentum transport equation takes the form of a Reynolds equation: where ν is the kinematic coefficient of molecular viscosity.The term ν∇ 2 U denotes the impact of viscous stress on the velocity, which is negligibly small in the situation studied.
The term ≪ u ′ u ′ ≫ multiplied by ρ is the Reynolds stress, which is of critical importance in the coastal zone because of a substantial contribution of the turbulence generated by breaking waves.
To simplify the model, the turbulent stress is calculated (according to the Boussinesq hypothesis) using an empirical relation with an average flow velocity (more precisely, with Cauchy's infinitesimal strain tensor): 1 /2[∇U + (∇U) T ] [37] where superscript T denotes the tensor transpose.The most general linear form of this constitutive relation can be written as follows: where N = N ijkl is the fourth-order tensor assigning an appropriate value of turbulent viscosity to each component, and •• is a double scalar product.It is worth noting the following: • As Equation ( 23) is linear and the wave field is known, the whole problem can be reduced to to entirely averaged flow velocities (after turbulent pulsations and wave oscillations have been filtered out) U = (U; V; W); • As parallel flow velocity changes are small, all the derivatives ∂ i U j (where i = j) can be disregarded such that the terms on the main diagonal of the Cauchy tensor vanish ; • As the coastal hydrodynamics is markedly diversified (the horizontal flow differs significantly from the vertical flow), Reynolds stress parameterization can be diversified as well by using the horizontal and vertical turbulent viscosities (ν h and ν v , respectively), where, in general, As the averaged vertical flow W is relatively low, it is omitted in the subsequent discussion.
The Reynolds stress tensor can be written as follows: The tensor divergence [a term in Equation ( 22)] is a vector.In compliance with the assumptions adopted, it can be written as follows: Because of the coastal flow diversification mentioned above, the horizontal momentum transport equation is written in vector notation, with scalar notation being used to write the vertical momentum transport equation.Consequently, an operator is introduced to the momentum transport Equation ( 22) in compliance with to the assumptions discussed above: In addition, boundary conditions for the instantaneous velocity vector ca be defined by the following: • Dynamic conditions: The model formulation requires defining the instantaneous pressure P, which is obtained from Equation (28).The left-hand side of the equation shows the second term to be negligibly small compared with the third term (horizontal changes are much smaller than those in the vertical plane); the first term can be disregarded as well, because when the equation is integrated with respect to z, a strong pressure fluctuation appears.
In its simplified form, Equation ( 28) is as follows: When averaged, it is as follows: Equation ( 31) is integrated with respect to z as the mass transport equation: The integration constant is determined based on the dynamic conditions on free surface (P = 0 at z = Z + ζ): Finally, decomposition into wave-and current-related parts and removal of all terms on the order higher than a 2 produce the following equation: The instantaneous pressure is therefore a sum of two values: • The average pressure (the accuracy is that of the terms on the order of a 2 , and only for z < Z): • Fluctuations of pressure related to the wave motion (the accuracy is that of the terms on the order of a): p = −ρgζ.
Just above the mean free surface level, z = Z + atmospheric pressure is, due to the low air density, equal to P = P + ≈ 0. However, at z = Z − (approaching the averaged level from below), the water pressure is lower than the atmospheric pressure.Therefore, a pressure difference occurs on the averaged free surface.
After decomposition, averaging and removal of nonlinear terms on the order higher than a 2 , and taking into account that ⟨uw⟩ = 0 (the linear wave theory shows the orbital velocities u and w to differ in phase by π /2), the form of the momentum transport equation is as follows: • Kinematic conditions: • Dynamic conditions: The term ∇ • ⟨uu⟩ in Equation (37) corresponds, in the three-dimensional space, to the product ∇ • ⟨uu⟩.It may be identified with the ∇• ≪ u ′ u ′ ≫ of the Reynolds Equation (22), i.e., with the divergence of kinematic turbulent stresses of the opposite signs.Consequently, quasi-turbulent wave stresses are identified and complemented with a wave-induced change in the average pressure.The stresses mentioned, along with the stresses occurring between the wave crest and trough, when integrated with respect to z, are known as the radiation stresses [38].
Substitution of Equation ( 35) in (37), the horizontal momentum transport equation results in the following form: Tensor R, which can be termed the quasi-turbulent stress tensor, is expressed as follows: where I is the unit tensor acting on the horizontal plane.Equation ( 41) is the averaged momentum transport equation which, together with the averaged mass transport Equation ( 12) and all the boundary conditions, is used to determine the three-dimensional velocity distributions of wave-induced currents in the coastal zone.
Tensor R corresponds to the local wave-induced force which generates the flow.For this reason, its components can be termed the depth-dependent radiation stresses.The components of R are calculated with the formulae defining individual wave motion velocities (cf.e.g., [39,40]): where Tensor F integrated with respect to z equals Therefore, the following can be obtained: S is the integrated form of the bottom to z = 0 − quasi-turbulent stress tensor or the vertical integrated internal radiation stress: where β = c g/c = 1 /2 1 + kD sinh −1 2kD .
To arrive at a completely balanced momentum, it is necessary to consider the momentum flow between the wave crest and trough.To this end, it is necessary to reiterate the following momentum transport equations before averaging: This produces the averaged form: Subtraction of the averaged Equations ( 50) and (51) from the non-averaged Equations ( 48) and (49), results in equations describing fluctuations (to the accuracy of a): With integrating (52) from the mean surface level Z (practically from 0) to the free surface Z + ζ (practically to ζ) and assuming that the pressure distribution in this region is hydrostatic p = ρg(ζ − z), the averaging with respect to the wave period results in the time derivative of the depth-dependent averaged volumetric transport as follows: where s = 1 /2EI.The transformation is indispensable at the next step of the model development.It involves integration of the basic horizontal momentum transport Equation ( 27) from the bottom to the instantaneous free surface, followed by the averaging of the wave motion scale.This is equivalent to the integration from the bottom to the average surface of the averaged horizontal momentum transport Equation ( 41), which obtains the following: Then, Equation ( 55) is added to the fluctuations Equation (54).The result is a wellknown (cf.Philips, [41]) momentum transport equation integrated vertically and averaged with respect to the wave period: where is the total wave-and current-generated mass transport divided by density, S is the radiation stress tensor, and G is the tensorial differential operator is the difference between the stress occurring just below the mean surface level τ − s and that on the bottom τ b τ + s is the stress acting just above the free surface which is practically equal to 0 because of the negligible impact of wind friction on the coastal zone hydrodynamics.All this leads to the final form of the depth-integrated momentum transport equations: or: In addition, τ − s is defined as follows: Stive and Wind [34] published a prototype of the relationship concerning the averaged stresses at the wave trough.It seems, however, that the prototype was not used on any larger scale in later research on wave-induced current field models, apart from some undertow analyses.In this study, the boundary condition (the value of ∂ z U) was defined based on Equation (62).
Similar to the mass transport Equation ( 20), the radiation stress in the left-hand part of Equation (61a,b) has a distribution that can be described as follows: where The sum of S [Equation (47)] and s results in a classical radiation stress: which shows that the modified quasi-turbulent stress tensor R in the form proposed, the stress at the mean surface level, the averaged stress at the mean surface level s, and the total radiation stress expressed using the Dirac delta function describe the vertical distribution of the current-inducing force correctly.This relationship in shown in Figure 2.

Special Solution
In a special case of the stationary state and the coastal zone where all the depth contours are parallel to the shore, it is possible to introduce additional simplifications which would lead to an analytical solution of the model.The simplifications involve the following:

•
Removal of all the time derivatives (∂ t = 0); • Removal of all derivatives with respect to coordinate y (∂ y = 0); • An assumption, based on the near-uniformity of alongshore velocity vertical distributions shown by measurements, that all the terms containing derivatives with respect to x only can be expressed as V ≈ Q y /D, resulting in the following: It is also assumed that With these assumptions, the momentum transport equation can be analytically integrated vertically.In preparation to the integration, the equation is written as two scalar equations: where α is the wave direction, and K = (2β − 1)D −1 E.

•
The first integration with respect to z, assuming a low bottom slope and ν h being constant vertically, results in the following: • The second integration (under the same assumptions) produces the following: • The third integration, within the boundaries from −D to 0 − , leads to the following: The final solution is obtained after determination of the constants C 1 , C 2 , C 3 , and C 4 .The solution assumes the bottom friction to be almost 0 and is developed with the knowledge that fields Z and Q = (Q x , Q y ).Substitution of z = −D after the first integration and removal of the friction terms results in the following: Other constants are determined from the last system of equations after inserting the known C 1 and C 2 : The knowledge of all the integration constants makes it possible to obtain the distributions of the velocities U and V.As no numerical solution of a three-dimensional problem is necessary, the situation is simplified.The final equations are as follows:

General Solution
Solving special case of three-dimensional model equations in the general case is very difficult.The major difficulty involves the necessity of the velocity field calculated to fulfill an additional integral condition.The solution was possible on account of a far-fetching reformulation of the problem.First of all, it was observed that the three-dimensional part of the model does not need to involve the continuity equation because the fields Q and Z from the stationary two-dimensional equations can be used and because the vertical velocity W can be, if necessary, calculated ex post (explicitly) from the equation.The following form of the stationary momentum transport Equation ( 41) was subject to transformation: where p subscript denotes the previous approximation.Then, the integral operator z = 0 − z . . .dz (in the case of multiple integrations: zz , zzz etc.) is used, which leads to the following equation: where τ − s = − 1 /2∇E.In addition, a new dependent variable Φ = z U is introduced to satisfy the following conditions: As a result, the equation discussed assumes a new form: All the terms on the right-hand side of the equation, except for the last one, are known.The last term is taken from the previous approximation in which it is assumed to be the following: The resultant equation can be integrated with respect to z in each point (x; y), and its subsequent integrals are as follows: Substitution of z = −D to all the three Equations ( 82), (84), and (86) produces the following: where Moreover, assuming a low vertical variation of U as confirmed by numerous coastal zone flow measurements, a depth-invariant term ρG • (Q/D) = X can be used instead of [ρG • U] p .It is then possible to calculate integrals of the expressions below: Consequently, the solution can be formulated as follows: These analytical equations allow for the solving of the general case based on numerical methods.

Results
The differential equations of the model for wave-generated currents in the coastal zone, along with the proposed solutions, have-due to containing quite significant elements of novelty-largely a prototype-exploratory character.Therefore, the computations performed in this research and the outcomes obtained are, to some extent, considered initial explorations.Their main purpose is to broadly confirm the validity of the model and its solutions.Given the limited scope of research concerning the calibration of the model and the challenges associated with the parameterization of turbulent viscosity-attributable to the notably incomplete characterization of turbulence within the coastal zone-the lack of relevant data in the existing literature, and the absence of adequate comparative measurement data, the results presented in this paper are to be interpreted as preliminary.
This study presents the calculation results for two conventional yet representative bathymetric configurations: (I) a coastal zone characterized by a uniformly sloping flat bottom and (II) a coastal zone with a multibar bathymetric profile (Figure 3).In both scenarios, it is assumed that all isobaths are parallel to the shoreline, which makes all derivatives with respect to the y-coordinate equal to zero.The model's input data were derived from measurements obtained during research expeditions at the Coastal Research Laboratory of the Institute of Hydro-Engineering of the Polish Academy of Sciences in Lubiatowo.For detailed insights into the environmental conditions of this area, see Dudkowska et al. [42].For the modeling of wave action in coastal zones I-II, the SWAN (Simulating WAves Nearshore) numerical wave model was used [43].The input data for the calculations in the presented examples consisted of actual deep-water wave conditions measured using a directional wave buoy WR, positioned at a distance of 1850 m from the shore at a of approximately 15 m during the Lubiatowo 2006 expedition (Table 1).Additionally, the calculations incorporated bathymetry and wind fields, treated as uniform (magnitude and direction) and measured during both expeditions.The calculations were performed on the outer grid, with a spatial step of 100 m × 100 m, which covered a strip with a width of 1850 m, characterized by a flat bottom with a constant slope.Additionally, calculations were carried out on an inner grid (smaller), encompassing an area with known bathymetry (I-II) with a spatial step of 10 m × 10 m.In a coastal zone characterized by a flat bottom with a constant slope and longitudinal homogeneity, the values of the averaged free surface elevation Z obtained from calculations should align with typical elevations associated with the presence of well-known phenomena such as set-up and set-down.Therefore, a primary significant test of the model was deemed to be the verification of calculation results conducted for bathymetric case I and various corresponding wave scenarios (C-F).Concurrently, a full range of bottom friction coefficient f values was employed in the calculations to assess its impact on the obtained values.The results of the calculations are illustrated in Figure 4.
For each wave scenario, based on the values of the averaged free surface elevation, it was possible to identify the location of the set-down zone (just before the wave breaking line) and the set-up zone (reaching its maximum near the shore).The value of set-down was close to zero with a bottom friction coefficient f = 0.01 and increases as the value of this coefficient decreases.This pattern allows us to conclude that the obtained results are generally consistent with the classic findings of Longuet-Higgins [32], obtained under the assumption of no bottom friction, and with the results of other studies on this topic (see, for example, Hsu et al. [44]).The Z values for scenarios C and D were similar, a consequence of the resemblance in wave height values for both cases.The strong dependency of the free surface elevation values on the function H(x) was confirmed by the results obtained for wave scenarios E and F. In these instances, a reduction in Z values was observed, as well as a shift of the minimum values (set-down) toward the shore.Based on the calculations performed, no distinct influence of the wave approach direction to the shore on the values of free surface elevation was detected.
To determine the impact of the cross-shore bathymetric variability on velocity values, the results obtained for coastal zones I-II and wave data A and B are summarized below.The presented results, for both wave data A and B, were achieved using f = 0.01.Given the certain similarities between wave fields A and B and the assignment of identical coefficient values (considering their significant impact on the outcomes), similar values for all presented magnitudes were obtained.The only observation worth noting is that the values obtained for case B are slightly higher.This is likely related to the less oblique wave approach to the shore.
The obtained values of free surface elevation (for conditions A and B) are consistent with the theoretical values and those measured under natural conditions (Figure 5).
Based on the results concerning flow velocity values for conditions A and B, it can be observed that they exhibit minimal vertical variability (see Figures 6 and 7).This can be attributed to the specified value of the vertical turbulent viscosity coefficient and is fairly consistent with the velocity distributions measured in the coastal zone.In coastal zones I and II, the velocities U are positive while those for V are negative, which for one, relates to the assumed longitudinal homogeneity and for another, to the wave approach direction to the shore.Furthermore, the values obtained for these cases (I and II) based on calculations A and B do not significantly differ.From this, it may be inferred that the analytical solution for specific conditions can be successfully applied in areas with isobaths parallel to the shore.

Summary and Conclusions
A three-dimensional model of wave-induced currents was developed from the following equations: They were also developed from the following boundary conditions: ) The results obtained in this work are expected to be of importance for studies on ocean dynamics.They are particularly applicable to the hydrodynamics of wave-dominated coastal zones.
The major novel aspects elucidated by the study can be summarized as follows: • Some new basic simplification rules appropriate for the problem considered were proposed; • A new simplified relationship (tailored to the nature of the coastal zone) that describes the turbulent stress tensor was formulated and applied to equations describing wave-current interactions; • Local (quasi-turbulent) wave-induced stresses were identified and treated as components of the three-dimensional radiation stress tensor; • The distribution of the wave-induced momentum transport (radiation stress) as a function of the depth-independent variable was described; • The presence of the virtual wave-induced stresses (the distributional part of the radiation stress tensor) at the averaged free surface was identified, and the stresses were quantified; • A fully three-dimensional model describing the coastal zone hydrodynamics was developed; • An approximate analytical solution for the special case of a coastal zone homogeneous in the longshore direction (with isobaths parallel to the shoreline) was proposed; • An approximate analytical solution for the general case was proposed.
Further research should focus on the development of a numerical model that solves the equations presented.

1.
A right-handed Cartesian coordinate system (x; y; z) was used; the x-axis is directed toward the open sea with 0 on the shoreline; the y-axis is parallel to the smoothed coastline; the z-axis is vertical, directed upward, with 0 at the mean sea level.

2.
Differential operators and a few integral operators are denoted in a compact format as and in the case of multiple uses as ∂ z ∂ z = ∂ zz , and as The case of multiple uses, they are denoted as z z = zz , z z z = zzz .3.
All the vectors and tensors are denoted by boldface characters, e.g., I, q, U.
Quantities averaged over the wave motion time scale are denoted by capital italicized characters (e.g., P, U); wave-generated oscillations are written with lower-case italicized characters (e.g., p, u); turbulent pulsations are written with lower-case italicized characters with apostrophes (e.g., u ′ ).Total quantities (the sum of averaged value and wave oscillations) are written with Roman characters (e.g., P, U).

7.
Distributions are denoted by a tilde over the Roman character (e.g., S), while their integrals are written with bold characters (e.g., S), as are the tensors.

Figure 1 .
Figure 1.Visual representation of the three-dimensional mean velocity distribution in the coastal zone.

Figure 2 .
Figure 2. Visual representation of the vertical distribution of the force generating the horizontal momentum transport.

Figure 3 .
Figure 3.The current model input data: water depth.I-flat-bottom, constant-slope profile; II-multibar profile with isobaths parallel to the shoreline.

Figure 4 .
Figure 4. Mean free surface elevation Z obtained from a general solution for wave parameters C-F and bathymetric data I. Calculations were performed for different values of bottom friction coefficient f .

Figure 5 .
Figure 5. Mean free surface elevation Z obtained from a general solution for wave parameters A-B and bathymetric data I-II.Calculations were performed for different values of bottom friction coefficient f .

Figure 6 .
Figure 6.Current velocity U (left) and V (right) for different values of σ (depth levels) obtained from calculations based on the general solution (upper figure) and the special solution (lower figure) corresponding to the wave parameters A and the bathymetric data I-II.

Figure 7 .
Figure 7. Current velocity U (left) and V (right) for different values of σ (depth levels) obtained from calculations based on the general solution (upper figure) and the special solution (lower figure) corresponding to the wave parameters B and the bathymetric data I-II.

Table 1 .
Deep water wave parameters measured during the Lubiatowo 2006 field campaign.H-wave height; T-wave period; α-angle of wave approach.