Eulerian Description of Wave-Induced Stokes Drift Effect on Tracer Transport

: The wave-induced Stokes drift plays a signiﬁcant role on mass/tracer transport in the ocean and the evolution of coastal morphology. The tracer advection diffusion equation needs to be modiﬁed for Eulerian ocean models to properly account for the surface wave effects. The Eulerian description of Stokes drift effect on the tracer transport is derived in this study to show that this effect can be accounted for automatically in the wave-averaged advection-diffusion equation. The advection term in this equation is the wave-averaged concentration ﬂux produced by the interaction between ﬂuctuations of linear wave orbital velocity and tracer concentration, and the advection velocity is the same as the Stokes drift velocity. Thus, the effective dispersion of tracers by surface gravity waves is calculated due to the Stokes drift effect and the corresponding dispersion coefﬁcient in the depth-integrated equation is then derived. The Eulerian description of Stokes drift effect of tracer concentration is illustrated by the direct numerical simulation of the advection–diffusion equation under simple linear waves. The equivalence between both the Eulerian and Lagrangian descriptions is also veriﬁed by particle tracking method. The theoretical analysis is found to agree well with the wave-induced dye drift velocity observed outside the surf zone in a longshore current experiment.


Introduction
The Stokes drift velocity, which is the time mean velocity of a fluid particle in water waves [1,2], may play a significant role on the mass/tracer transport (e.g., larval, plankton, oil, plastic pollution, etc.) in the ocean and the evolution of coastal morphology. For instance, Monismith and Fong [3] showed through two field experiments that the onshore Stokes drift near the water surface could be stronger than measured Eulerian currents on the inner shelf, and it is necessary to include the Stokes drift to explain an observed shoreward migration of a dye plume trajectory. Although the Lagrangian Stokes drift cannot be detected directly by current measurements at a fixed position, it is suggested to be explicitly considered in the models that attempt to predict the transport of organisms and solutes in the nearshore environment as a surface wave effect. Rohrs et al. [4] presented that the Stokes drift exhibits a more constant direction in the long term than the Eulerian current, and the average distance to the coast of Northeast Arctic cod eggs and larvae decreases by 20% on average due to the wave effects. Reniers et al. [5] showed that when Stokes drift was considered in the model of the computed drifter trajectories, good agreement was found with the observed retention of surf zone material on a rip-channel beach. Without Stokes drift considered, however, the estimated results overpredicted the number of observed drifter exits by almost an order of magnitude during the drifter deployments. This emphasizes the importance of including Stokes drift when calculating the transport of floating materials on the ocean surface. Henderson et al. [6] showed that the Stokes drift could increase shoreward transport of suspended sediment and contribute to the shoreward bar migration. The Stokes drift is also a cross-shore exchange mechanism between the surf zone and the inner shore. Its importance relative to the transient rip currents was examined by using temperature and dye measurements at an approximately alongshore-uniform beach [7]. The dye release was simulated by the wave-resolving funwaveC model, but it does not include Stokes-drift driven exchange [8]. The tracer advection equation has been modified for Eulerian ocean models to account for Stokes drift effect by McWilliams et al. [9], but they did not consider the diffusion terms.
The variation of the Stokes drift velocity with water depth gives rise to the Taylor dispersive effect [10] in the longitudinal direction. Law [11] studied the dispersion of pollutants in the free surface waves. By assuming that the particles move passively with the water particles in the same phase and amplitude, the time-averaged velocity profile along the water depth was taken as the Stokes drift velocity profile. Taylor dispersion coefficient is obtained by using the formula derived by Van den Broeck [12]. A random walk method was also used to simulate the movement of pollutant particles in the surface waves. In the process of solving the particle position, the particle displacement at each time step is obtained by time integration of the instantaneous wave orbital velocities, so the effect of Stokes drift velocity is included. The results of the dispersion coefficients obtained by the above two methods are very close for the steady state, which indicates that the two methods are equivalent to describe the Stokes drift effect of pollutant diffusion under waves and the effects of wavy surface boundary can be ignored. For theoretical analysis, Shen and Zou [13] took the advection velocity in the vertical two-dimensional advective diffusion equation as the sum of the Stokes drift velocity and the linear wave orbital velocity, and derived the Taylor dispersion caused by the surface waves using a successive approximation method. The wave-averaged dispersion coefficient in their result is 35% larger than that of Law [11]. However, the sum of the Stokes drift velocity and the linear wave orbital velocity may be not properly used in the advection diffusion equation for the wave resolving concentration. Assuming the pollutant to be vertically well mixed, Winckler et al. [14] derived a two-dimensional depth-averaged advection diffusion equation for wave-averaged concentration in the shallow water of surf zone. Their theory shows that in addition to longshore current velocity induced by wave breakings, the advection velocities also include contributions from the depth-averaged Stokes drift velocities, which are due to the covariance of waves and fluctuations of the pollutant concentration. Within the vertically well-mixed surf zone, Stokes drift was found to be a negligible cross-shore dye dispersion mechanism relative to surf zone eddies [15]. However, outside the surf zone on the inner-shelf, vertically varying Stokes drift velocity could potentially be an important cross-shore dispersion mechanism.
In tracer transport models, the Stokes drift of gravity surface waves is usually considered using Lagrangian description [4,5]. However, the Eulerian concentration transport equation is also commonly used, so a Eulerian description for Stokes drift effect is necessary. This paper explores how to consider the effects of the Stokes drift in the Eulerian advective diffusion equation of tracer concentration and discusses the equivalence of the Eulerian description of the Stokes drift effect to the Lagrangian method. A multiple-scale analysis is used to derive the wave-averaged vertical two-dimensional and depth-averaged onedimensional advective diffusion equations of tracer transport in wave field. The Eulerian description of Stokes drift effect of tracer concentration is illustrated by numerically solving the advection diffusion equation in linear wave field directly. The theoretical results are also compared with the results from particle tracking method based on the Lagrangian description. A laboratory experiment is also presented to show that the dye drift induced by oblique incident waves was obvious outside the surf zone on a planar beach.

Eulerian Description of Stokes Drift Effect
The generation mechanism of Stokes drift effect on the tracer transport in the wave field by the Eulerian description is given in this section. A multiple-scale analysis [16] is used for advection diffusion equation of the tracer concentration in the surface gravity wave field to obtain the transport equation of the wave-averaged concentration.
We assume that no chemical reaction is presented, and the tracer density is the same as that of the surrounding fluid or the sea water. The concentration of passive and conservative tracer, c(x, z, t), in the gravity waves is assumed to be governed by the advection diffusion equation [17] ∂c ∂t where x and z represent the horizontal and vertical coordinates, u and w represent the horizontal and vertical velocity components, D x and D z are the horizontal and vertical diffusivities assumed to be constant. No-flux boundary conditions are applied to the impermeable bottom z = −h and water surface z = η: According to the characteristics of waves and tracer transportation outside the surf zone, the normalization of following variables is introduced: where A, k, T, h, c 0 , and D v are characteristic wave amplitude, wave number, wave period, water depth, tracer concentration, and vertical diffusivity, respectively. In terms of these dimensionless variables, Equation (1) can be rewritten non-dimensionally as: As the small-amplitude wave and the finite water depth outside the surf zone are considered, so ε = kA and kh = O(1) are assumed. Many field studies reported that the vertical diffusivity is small (e.g., D v = 10 −4~5 × 10 −3 m 2 /s [18]), so the ratio of wave period T to the mixing time h 2 /D v in a finite water depth(namely a time-scale for achieving a uniform condition between bottom and water surface) TD v /h 2~O (ε 2 ) is assumed. For example, considering the following characteristics of waves and vertical diffusivity: h = 3 m, H = 0.6 m, T = 5 s, D v = 0.005 m 2 /s, then kh = 0.75~O(1), TD v /h 2 = 0.002778, ε 2 = 0.005715, so we have TD v /h 2~O (ε 2 ).
Having identified the orders of these terms in Equation (5) and neglecting the higherorder terms, Equation (5) is then rewritten dimensionally, but in the simplified form of Equation (1), by introducing an ordering parameter ε, where the boundary conditions in Equations (2) and (3) are rescaled as: The order contrast among the terms in Equation (6) indicates that there are two distinct timescales in the advective diffusion process of passive tracer. One is wave period T, which characterizes the fast oscillations of tracer caused by surface gravity waves, and the other is mixing time h 2 /D v in the finite water depth, which characterizes the slow diffusion caused by the diffusivities. The ratio between the two timescales is TD v /h 2~O (ε 2 ). Because of the significant difference between the time scales of wave period and diffusion, we introduce multiple-scale coordinates for time: t for the wave motion and τ = ε 2 t for the effects of diffusion. Thus, the tracer concentration c(x, z, t, τ) depends on (t, τ) as independent variables. The original time derivative becomes, according to the chain rule, ∂c/∂t → ∂c/∂t + ε 2 ∂c/∂τ . The wave velocity, wave surface and tracer concentration are expanded as follows: where u (n) , w (n) and η (n) are nth harmonics of wave horizontal and vertical orbital velocities and surface wave profile, and c (n) = c (n) (x, z, t, τ). Stokes wave theory [19] is used, and we have and u (2) , w (2) and η (2) are second harmonics. Equating the coefficient of each power of ε to zero leads to a sequence of perturbation equations. At the leading order c (0) (x, z, τ) is expected to represent the period average and depends only on τ. The transport equation at O(1) is trivially satisfied by the period-averaged c (0) (x, z, τ) and ∂c (0) /∂z = 0 at z = 0 and z = −h. At O(ε), c (1) (x, z, t, τ) represents the perturbation due to the oscillating velocity field and satisfies ∂c (1) and the boundary conditions: and At O ε 2 , c (2) satisfies and the boundary conditions and ∂c (2) Substituting Equations (13), (14) into Equation (16), we get Equation (19) has a solution c (2) that is periodic in x and t if and only if the following solvability condition holds: where denoting the wave-average by overline, f = t+T t f dt/T . Substituting Equation (22) into Equation (19), we get This equation is a transport equation for wave-averaged concentration, c (0) , with the advection velocity Notice that the advection term, u c ∂c (0) /∂x, is the phase-averaged concentration flux produced by interaction between fluctuations of linear wave orbital velocities and tracer concentration and the advection velocity is identical in form to the Stokes drift velocity with the definition given by Longuet-Higgins [2] This means that the wave-averaged transport in wave velocity field can be expressed by the transport by Stokes drift. The former Equation (25) is a Eulerian description whereas the latter Equation (27) is a Lagrangian description, which means the two are equivalent for the case of considering the wave-averaged tracer transport problem.
The variation of the advection velocity with water depth gives rise to the Taylor dispersive effect [10] in the longitudinal direction, and a dispersion effect will also be generated in the depth-averaged equation where a is the depth-averaged tracer concentration, U is the depth-averaged velocity, K is the Taylor dispersion coefficient when the dispersing cloud becomes normally distributed approximately. Their expressions are: whereĉ (0) andû c are deviations of c (0) and u c from their depth-averaged values. The expression for K is the same as that was derived analytically by Law [11] using a different method of Van den Broeck [12]. Chatwin [20] has also shown that the dispersing cloud will become normally distributed approximately a time h 2 /D z after its injection. If the concentration at time is the initial condition for Equation (28) (M represents the total amount of the diffusing substance and σ 2 is the variance), the solution of Equation (28) can be expressed as At this stage, the deviation of concentration can be given byĉ (0) , then At this stage, a theoretical solution to Equation (1) can be derived as: The following two sections will give the results of the direct numerical solution of the advection diffusion Equation (1) and the Lagrangian description results of the particle tracking method, and further explain the equivalence of above Eulerian description of the Stokes drift effect to Lagrangian description.

Numerical Solution of the Advection Diffusion Equation
In the previous Section 2, the advection term (26) in the advection diffusion equation of the period-averaged concentration comes from the interaction of the first-order wave orbital velocities and corresponding concentration fluctuation. In order to visualize the physical process described by the Eulerian description, the concentration is directly simulated by numerically solving Equation (1) in the first-order wave field given by Equations (13)- (15).
In order to consider the variation of the calculation domain caused by the wave free surface, the calculation domain is transformed into a fixed calculation domain through the following σ coordinate transformation: By applying the chain rule to Equation (1), the advection diffusion equation in σcoordinates can be expressed as: The up and down boundary condition can be expressed as: and left and right boundary conditions are: The initial condition of a line source can be expressed by: where M represents the total amount of the diffusing substance, δ is a Dirac delta function, and the initial line source is at x = 20 m. The finite difference method is used to discrete the terms of Equation (37). ADI (alternating direction implicit) scheme is used to solve the equation, and the numerical scheme has been validated for a steady flow. In this calculation, water depth h = 3 m, wave height H = 0.6 m, wave period T = 5 s, and D x = D z = 0.005 m 2 /s are used. The length of model domain is 80 m. The grid spacing along the depth is 0.1 m, and the grid spacing along the direction of wave propagation is 0.1 m. The time step is 0.01 s. Figure 1 shows the concentration distribution in the flow field at time t = 800 s and the horizontal displacement of the concentration center caused by this Stokes drift effect. It can be seen that while the concentration diffuses under the action of waves, there is also a drift along the wave propagation direction. The drift velocity is the same as the Stokes drift velocity in Equation (27). This is illustrated here by estimating the depth-averaged velocity of this drift where C is the depth-averaged concentration, c(i, j) is the concentration at the grid point (i, j), and ∆ξ the horizontal grid spacing. The tracer concentration drift velocity obtained from Equation (41) is similar to the depth-averaged drift velocity of 0.03 m/s obtained by Equation (30) and that by the particle tracking method under the same conditions. The wave velocity field in the calculation is taken as the first-order wave orbital velocity components of Equations (13) and (14), but the second-order Stokes drift velocity is not considered. The numerical simulation has indicated that although the Stokes drift velocity is not added, the drift effect can be formed by the interaction of the first-order wave orbital velocity of the periodic wave and the concentration component caused by the periodic waves (i.e., caused by the convection term u∂c/∂x + w∂c/∂z in Equation (1)). where C is the depth-averaged concentration, ( , ) c i j is the concentration at the grid point ( , ) i j , and ξ Δ the horizontal grid spacing. The tracer concentration drift velocity obtained from Equation (41) is similar to the depth-averaged drift velocity of 0.03 m/s obtained by Equation (30) and that by the particle tracking method under the same conditions. The wave velocity field in the calculation is taken as the first-order wave orbital velocity components of Equations (13) and (14), but the second-order Stokes drift velocity is not considered. The numerical simulation has indicated that although the Stokes drift velocity is not added, the drift effect can be formed by the interaction of the first-order wave orbital velocity of the periodic wave and the concentration component caused by the periodic waves (i.e., caused by the convection term ∂ ∂ + ∂ ∂ / / u c x w c z in Equation (1)).

Modeling Tracer Transport by Particle Tracking Method
The principle of particle tracking method is to calculate the trajectories of particles in the flow field and can be used to determine the concentration distribution. In this Lagrangian method, the Stokes drift is reflected by the wave-averaged particle displacement. The particle tracking method [21] uses the governing equation: where  ( )=( ( ), ( )) Y t x t z t are the coordinate of particle positions,  ( , ) u Y t and  ( , ) w Y t are the first-order wave orbital velocities, A1 and A2 are independent random numbers from the standard normal distribution with zero mean and unit variance, and Dx and Dz are the horizontal and vertical diffusion coefficients.
An example of simulating tracer transport in a surface wave field by particle tracking method is given here to show how the Stokes drift phenomenon is realized in the Lagrangian description. The calculation conditions of water depth, wave conditions and diffusion coefficients are taken as the same as those in Section 3. There are 6000 particles released in the wave field. The initial horizontal coordinates of the particles are located at x = 0, and the random vertical coordinates are uniformly distributed along the water depth.

Modeling Tracer Transport by Particle Tracking Method
The principle of particle tracking method is to calculate the trajectories of particles in the flow field and can be used to determine the concentration distribution. In this Lagrangian method, the Stokes drift is reflected by the wave-averaged particle displacement. The particle tracking method [21] uses the governing equation: where Y (t) = (x(t), z(t)) are the coordinate of particle positions, u( Y, t) and w( Y, t) are the first-order wave orbital velocities, A 1 and A 2 are independent random numbers from the standard normal distribution with zero mean and unit variance, and D x and D z are the horizontal and vertical diffusion coefficients. An example of simulating tracer transport in a surface wave field by particle tracking method is given here to show how the Stokes drift phenomenon is realized in the Lagrangian description. The calculation conditions of water depth, wave conditions and diffusion coefficients are taken as the same as those in Section 3. There are 6000 particles released in the wave field. The initial horizontal coordinates of the particles are located at x = 0, and the random vertical coordinates are uniformly distributed along the water depth. The time step is 0.25 s, and the spatial positions of the particles at any time can be calculated according to Equation (42). The fourth order Runge-Kutta method is used to solve Equation (42), and the second term is processed using the random walk method. The numerical scheme is validated by theoretical results of Taylor Dispersion in steady channel flow. If the particles move beyond the free surface boundary and the bottom boundary due to random walk, they are mirror reflected into the water domain. Figure 2 shows the spatial distribution of particles and concentration distribution in the linear wave field at t = 800 s corresponding to Figure 1. The horizontal center position of the particle distribution x can be obtained from the spatial position distribution of the particles, and the average speed U of the overall particle drift can be obtained by dividing x by the corresponding time interval t. Figure 3 shows the value given by the Eulerian analytical solution of Equation (30) and that obtained by numerical simulation of U. The two values tend to be the same, which indicates that there is an overall drift with the velocity of Equation (30) as the advection velocity during the particle diffusion process in the wave field. erage speed U of the overall particle drift can be obtained by dividing x by the corresponding time interval t. Figure 3 shows the value given by the Eulerian analytical solution of Equation (30) and that obtained by numerical simulation of U. The two values tend to be the same, which indicates that there is an overall drift with the velocity of Equation (30) as the advection velocity during the particle diffusion process in the wave field.  The reason for the above consistency is that the overall drift velocity U of the particles is determined by the first term at the right end in Equation (42), while the second term does not affect the averaged motion in a statistical sense. Therefore, the Stokes drift velocity can be obtained only by the time integration of the first term (by tracking particles at different water depths). Figure 4 shows the distribution of the Stokes drift velocity along erage speed U of the overall particle drift can be obtained by dividing x by the corresponding time interval t. Figure 3 shows the value given by the Eulerian analytical solution of Equation (30) and that obtained by numerical simulation of U. The two values tend to be the same, which indicates that there is an overall drift with the velocity of Equation (30) as the advection velocity during the particle diffusion process in the wave field.  The reason for the above consistency is that the overall drift velocity U of the particles is determined by the first term at the right end in Equation (42), while the second term does not affect the averaged motion in a statistical sense. Therefore, the Stokes drift velocity can be obtained only by the time integration of the first term (by tracking particles at different water depths). Figure 4 shows the distribution of the Stokes drift velocity along The reason for the above consistency is that the overall drift velocity U of the particles is determined by the first term at the right end in Equation (42), while the second term does not affect the averaged motion in a statistical sense. Therefore, the Stokes drift velocity can be obtained only by the time integration of the first term (by tracking particles at different water depths). Figure   From the above results, the dispersion coefficient in the depth-averaged Equation (28) can also be determined. From the above results, the dispersion coefficient in the depth-averaged Equation (28) can also be determined.
where x is the horizontal center position of the particles, and N is the number of simulated particles. The condition for the above formula to be applicable with reasonable accuracy is t ≥ 0.4h 2 /D z [20], so the value calculated in the initial stage is not considered. Figure 5 shows the temporal variation of K calculated by the above formula, which approaches a stable constant (about K = 0.0108 m 2 /s) after the initial growth stage. The K value is found to be very close to K = 0.0115 m 2 /s, calculated from Equation (31). In this way, the dispersion effect expressed by Equation (31) caused by the linear surface wave field is illustrated well in Figures 2 and 5. From the above results, the dispersion coefficient in the depth-averaged Equation (28) can also be determined.
where x is the horizontal center position of the particles, and N is the number of simulated particles. The condition for the above formula to be applicable with reasonable ac- [20], so the value calculated in the initial stage is not considered.

Experimental Results and Discussion
A laboratory experiment was conducted to study the effects of wave-induced currents on the dye transport outside and inside the surf zone [22]. The experiment was carried out in the wave basin of the State Key Laboratory of Coastal and Offshore Engineering. Figure 7 shows the plain view of the experimental setup. The wave basin is 55 m long, 34 m wide, and 0.8 m deep. At one end of the basin was a piston-type wave maker, and the other end was the wave absorber. A concrete beach model was built with a 30° angle to the wave maker. The beach slope was 1:100, and the water depth at horizontal bottom was 0.18 m. The incident wave height used in this paper was 0.05 m, and the wave period

Experimental Results and Discussion
A laboratory experiment was conducted to study the effects of wave-induced currents on the dye transport outside and inside the surf zone [22]. The experiment was carried out in the wave basin of the State Key Laboratory of Coastal and Offshore Engineering. Figure 7 shows the plain view of the experimental setup. The wave basin is 55 m long, 34 m wide, and 0.8 m deep. At one end of the basin was a piston-type wave maker, and the other end was the wave absorber. A concrete beach model was built with a 30 • angle to the wave maker. The beach slope was 1:100, and the water depth at horizontal bottom was 0.18 m. The incident wave height used in this paper was 0.05 m, and the wave period was 1.5 s. The regular waves were propagating at an angle of 30 • to the shore. The three measuring carriages on which wave gauges (at the sections y = 4 m, 8 m, 12 m) and velocity meters (at the sections y = 8 m, 12 m) were mounted were perpendicular to the shoreline. tory tends to be parallel to the coast, which is shown in the last two photos in Figure 9e,f.
Long-term inner shelf observations [23] show that the offshore Eulerian flow can be roughly equal in magnitude but opposite in direction to the onshore Stokes drift at all depths, resulting in no net exchange. However, in the present experiment the effect of Stokes drift on dye transport is evident outside the surf zone when the incident waves are oblique to the coast, as the longshore current will take away the mass transport due to the Stokes drift. The measured mean Eulerian cross-shore velocity is smaller than the Stokes drift velocity outside the surf zone on inner shelf. This experiment is different from those experiments conducted in the wave flume [24,25] or basin [26] where incident waves are normal to the coast and the undertow or rip current velocity may decelerate the shoreward-moving tracers just outside the surf zone. Therefore, the above Stokes drift effect should be accounted for in the wave-averaged concentration flux of tracer along the wave direction besides the wave-averaged mean Eulerian velocity. This is also justified by the field observations given by Monismith [3] for the potential transport of solutes and organisms by surface waves.   Figure 8a). The dye was continuously released, and the dye-patch image was recorded by a chargecoupled device (CCD) at a height of about 8 m above the water surface. Figure 9 shows a series of CCD photos that was taken every 15 s. These photos show that the dye exhibits a shoreward drift under the unbroken waves. As the mean Eulerian velocity outside the surf zone is nearly zero (see Figure 8b), the drifting distance can be calculated by the Stokes drift velocity. The time interval between the fourth photo Figure 9d and the first photo Figure 9a is ∆t = 45 s, and the drift velocity at the corresponding water surface [u c ] z=0 ≈ 0.048 m/s. It is calculated according to Equation (25) by taking the wave height at the midpoint of the dye trajectory line as 0.055 m and water depth of h = 0.1 m. As the depth variation of Stokes effect is obvious outside the surf zone, the drift velocity at the corresponding water surface is used to calculate the horizontal length of the dye plume. The dye drifting distance from the two figures is ∆x = 45 0 u c dt ≈ 2.36 m. This distance is about 1.7 times the local wavelength of 1.44 m. This result can be also seen from Figure 9, because the diagonal bright lines in the figure are the wave crest lines, and the dye movement distance from the first photo to the fourth photo is two times the distance between two wave crest lines (wavelength). This demonstrates that the dye drift outside the surf zone is mainly controlled by the Stokes drift effect. However, the dye is not only driven by the Stokes drift but also by the longshore current (the direction is parallel to the coast, see Figure 8c) induced by the wave breaking after entering the surf zone. The dye trajectory tends to be parallel to the coast, which is shown in the last two photos in Figure 9e,f.

Conclusions
This study quantitatively investigates the Eulerian description of the Stokes drift effect in the advection diffusion of tracer concentration under surface water waves. The Stokes drift effect is produced by the interaction between the periodic wave velocities and the periodic concentration variation and can be naturally included when the wave field is taken as the first-order wave orbital velocity. This effect appears explicitly in the waveaveraged advection diffusion equation of concentration as the advection velocity is just the Stokes drift velocity. This velocity shows that the concentration drift in the direction of wave propagation in addition to the periodic fluctuations induced by the surface waves. This process of tracer advective diffusion is verified by numerical calculation of the Eulerian vertical two-dimensional advective diffusion equation and the Lagrangian particle tracking method in this research. The Taylor dispersion coefficient produced by the Stokes drift effect in the depth-averaged equation is also obtained from the above theoretical analysis and is also consistent with the results obtained in the particle tracking method and previous studies. The Stokes drift effect by surface water waves outside the surf zone is shown to be significant by the photos of shoreward drifting of dye in the

Conclusions
This study quantitatively investigates the Eulerian description of the Stokes drift effect in the advection diffusion of tracer concentration under surface water waves. The Stokes drift effect is produced by the interaction between the periodic wave velocities and the periodic concentration variation and can be naturally included when the wave field is taken as the first-order wave orbital velocity. This effect appears explicitly in the waveaveraged advection diffusion equation of concentration as the advection velocity is just the Stokes drift velocity. This velocity shows that the concentration drift in the direction of wave propagation in addition to the periodic fluctuations induced by the surface waves. This process of tracer advective diffusion is verified by numerical calculation of the Eulerian vertical two-dimensional advective diffusion equation and the Lagrangian particle tracking method in this research. The Taylor dispersion coefficient produced by the Stokes drift effect in the depth-averaged equation is also obtained from the above theoretical analysis and is also consistent with the results obtained in the particle tracking method and previous studies. The Stokes drift effect by surface water waves outside the surf zone is shown to be significant by the photos of shoreward drifting of dye in the laboratory longshore current experiment. The modified tracer advection diffusion equa- Long-term inner shelf observations [23] show that the offshore Eulerian flow can be roughly equal in magnitude but opposite in direction to the onshore Stokes drift at all depths, resulting in no net exchange. However, in the present experiment the effect of Stokes drift on dye transport is evident outside the surf zone when the incident waves are oblique to the coast, as the longshore current will take away the mass transport due to the Stokes drift. The measured mean Eulerian cross-shore velocity is smaller than the Stokes drift velocity outside the surf zone on inner shelf. This experiment is different from those experiments conducted in the wave flume [24,25] or basin [26] where incident waves are normal to the coast and the undertow or rip current velocity may decelerate the shoreward-moving tracers just outside the surf zone. Therefore, the above Stokes drift effect should be accounted for in the wave-averaged concentration flux of tracer along the wave direction besides the wave-averaged mean Eulerian velocity. This is also justified by the field observations given by Monismith [3] for the potential transport of solutes and organisms by surface waves.

Conclusions
This study quantitatively investigates the Eulerian description of the Stokes drift effect in the advection diffusion of tracer concentration under surface water waves. The Stokes drift effect is produced by the interaction between the periodic wave velocities and the periodic concentration variation and can be naturally included when the wave field is taken as the first-order wave orbital velocity. This effect appears explicitly in the wave-averaged advection diffusion equation of concentration as the advection velocity is just the Stokes drift velocity. This velocity shows that the concentration drift in the direction of wave propagation in addition to the periodic fluctuations induced by the surface waves. This process of tracer advective diffusion is verified by numerical calculation of the Eulerian vertical two-dimensional advective diffusion equation and the Lagrangian particle tracking method in this research. The Taylor dispersion coefficient produced by the Stokes drift effect in the depth-averaged equation is also obtained from the above theoretical analysis and is also consistent with the results obtained in the particle tracking method and previous studies. The Stokes drift effect by surface water waves outside the surf zone is shown to be significant by the photos of shoreward drifting of dye in the laboratory longshore current experiment. The modified tracer advection diffusion equation should be considered for Eulerian ocean models to account properly for the surface wave effects.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.