1. Introduction
a. Initial information
Surface drift currents play a significant role in marine shipping, the spread of surface waste, and marine industry safety. Thus, studying and understanding this phenomenon, which was comprehensively described in the recent review by Bremer and Breivik (2017), is of physical and practical importance. From the physical perspective, the drift currents are interesting due to their observability, availability for measurements, and possibility of theoretical description (e.g., Stokes 1847; Longuet-Higgins 1953; Craik 1982; Bremer and Breivik 2017, and references therein). For this reason, a large number of surface drift current measurements have been obtained over long periods, started by Shemdin (1972), Wu (1975, 1983), Churchill and Csanady (1983), and Bye (1988), and continued at present (e.g., Longo 2012; Longo et al. 2012; Zavadsky and Shemer 2017). The earlier work was well described by Tsahalis (1979), and an extended list of recent work was reported by Bremer and Breivik (2017). In addition, extensive theoretical work has addressed numerous aspects of the dynamics of currents induced by waves in water, namely, the generation of current vorticity due to viscosity, particle trajectories in water, turbulence induced by waves, the oscillation of boundary layers, the occurrence of return currents in a tank, the origin of Langmuir circulations, and so on. Most of these topics have been well described in certain papers (e.g., Longuet-Higgins 1953; Craik 1982; Teixeira and Belcher 2002; Babanin and Haus 2009; Teixeira 2018), surveys (e.g., Sullivan and McWilliams 2010; Bremer and Breivik 2017), and monographs (e.g., Phillips 1977; Janssen 2004; Bühler 2014). However, in this paper, we consider measurements of time-averaged surface currents induced by waves in water. The focus of the study is on the experimental aspects of the surface current problem rather than particle trajectories near the surface, turbulence generation by waves, or the depth profile of induced currents.
Despite the evident simplicity of surface drift current measurements, accurate measurements of a wavy water surface can have many sources of error. Field measurements (e.g., Churchill and Csanady 1983; Babanin 1988; Malinovsky et al. 2007; Kudryavtsev et al. 2008) are influenced by the nonstationarity of the wind field, uncontrolled background currents Ub, the stratification of air and water in situ, numerous technical difficulties in field measurements, and other factors. In laboratory (tank) measurements (Toba 1973; Wu 1975; Tsahalis 1979; Longo 2012; Longo et al. 2012; Paprota et al. 2016; Zavadsky and Shemer 2017), the sources of errors include the influence of the upper and lateral boundaries of the tank, presence of return currents, short wind and wave fetches, and difficulty of placing equipment in a tank. Some specific problems related to surface current measurements will be discussed in this paper.
All these limitations influence the wind and current profiles and establish dependencies of the surface drift on the wind and wave parameters of the air–water interface. Besides, small tank dimensions limit the range of wave state parameters, thereby decreasing the ability to identify the abovementioned dependencies. However, tank measurements have the evident advantage that all parameters in the air–water system can be completely controlled during an experiment. Therefore, the most exact and reliable results for surface drift currents have been obtained in laboratory measurements (Wu 1975, 1983; Longo 2012; Longo et al. 2012; Zavadsky and Shemer 2017). Nevertheless, many features of the drift current phenomenon are still not completely understood. In this paper, some of these issues will be addressed based on our own tank measurements.
b. Specification of the tasks
For completeness, consider two principally different cases: (i) mechanical waves and (ii) wind waves in water.
Regarding mechanical waves, one may expect that the surface drift current Ud should be of the same order as the Stokes mass-transport current USt if the ideal fluid assumptions are applicable (Stokes 1847). In the potential approximation and the case of a monochromatic nonlinear gravity wave with amplitude a, frequency ω, and wavenumber k = ω2/g (g is the gravity acceleration), the theory predicts the ratio (Stokes 1847; Phillips 1977)
where the first equality corresponds to the scaling of USt by the phase velocity, ω/k, and the second equality corresponds to the scaling of USt by the orbital velocity, ωa. The factor ka is the steepness of the wave (typically ~0.1), which determines the scale of the Stokes drift with respect to the mentioned velocities. For the wave spectrum S(ω), Eq. (1) for USt can be generalized. In the case of quasi-monochromatic (narrowband) waves, Eq. (1) can be replaced by the following formula:
where ωp and kp are the spectral peak frequency and the peak wavenumber, respectively. In the case of mechanical waves with wideband spectra, such as Pierson–Moscowitz (PM) or Joint North Sea Wave Observation Project (JONSWAP; JW) spectra (Komen et al. 1994), the Stokes drift can be estimated by the following formula (Churchill and Csanady 1983):
Here, [ωmin, ωmax] is the available frequency band of the wave spectrum S(ω), and z is the vertical coordinate directed positively upward. Equation (3) allows the Stokes drift for a given spectrum S(ω) to be estimated at both the surface and at any depth z; besides, it shows that the Stokes drift decays exponentially with depth z. Note that ratios (2) and (3) are also used for the estimation of USt(0) in the presence of wind waves (Wu 1983; Toba 1973; Clarke and Gorder 2018).
Unfortunately, direct measurements of real surface drift currents Ud, induced by mechanical waves, are rather limited (according to Bremer and Breivik 2017). For example, Craik (1982) described the depth profile formation of currents induced by different kinds of waves in viscous and contaminated water. Herterich and Hasselmann (1982) provided results for the diffusion of tracers by surface waves. Monismith et al. (2007) presented results for return current formation in a water body of finite depth. Paprota et al. (2016) and Breivik et al. (2014) measured Stokes drift profiles. In all these works, the authors noted a discrepancy between the measured currents USt(z) and the theoretical results (e.g., Longuet-Higgins 1953; Craik 1982; Constantin 2006). However, none of the researchers did quantify the measured surface drift current Ud in comparison with the theoretical value USt(0). This relation needs further specification, and this issue will be addressed in this paper based on tank measurements.
In the case of wind waves, the total surface drift (in the absence of background currents) includes two terms (the case of one-dimensional tank measurements is assumed):
where Udw is the wind-induced part of the surface drift current Ud. To estimate Udw, the Stokes drift needs to be excluded from Ud in Eq. (4); to this aim, the values of USt(0) are often estimated using Eq. (2). Based on such estimates, previous studies (Wu 1975, 1983; Churchill and Csanady 1983; Longo 2012; Longo et al. 2012; Zavadsky and Shemer 2017) showed that the Stokes drift at the surface USt(0) has a typical empirical value of approximately 10%–15% of the wind-induced drift Udw. This indicates that USt should be taken into account in measurements and practical tasks, especially in the presence of intensive long waves that have high phase velocities (e.g., swell). Moreover, in the field measurements by Clarke and Gorder (2018), the Stokes drift was considered as the total surface drift current, with no attention given to the wind-induced part Udw.
In practice, the wind-induced drift currents Udw have values on the order of the friction velocity in air
where coefficient cd varies in the range from 0.3 in the laboratory measurements (Zavadsky and Shemer 2017) to 0.53 (Wu 1975), and from 0.2 (Babanin 1988) to 1.5 (Tsahalis 1979) in the field measurements. In all cases, no authors considered the dependence of Udw on specific wave parameters. In the present paper, we address this relation in detail.
Until recently, it was unclear why the surface drift current Udw has values on the order of the friction velocity in air
In Polnikov (2019), a simple semi-phenomenological model was constructed based on the theoretical concept of the three-layer structure of the air–water interface (Polnikov 2011). Processing the numerical simulation data obtained by Chalikov and Rainchik (2011) and Polnikov (2011) showed that the profile of the total drift current, Ud(z), corresponding to Eq. (4), is linearly related to depth z in the so-called “wave zone,” where the air and the water alternatively exist. This finding was empirically confirmed by the measurements of Longo (2012) and Longo et al. (2012) (see details in Polnikov 2019). This means that the wave zone plays the role of a friction layer between the air and the water layers.
Polnikov (2019) proposed that in the wave zone the following ratio is valid
where
where Cp is a constant of the order of air-to-water density ratio ρa/ρw. Thus, the empirical ratio (5) describing the relation between Udw and
One may expect that the wind-induced drift velocity Udw is dependent on the following dimensionless parameters of the wave state: the steepness σ = akp, the wave age
Based on our own direct drift current measurements obtained with surface floats in a tank, the above dependencies were successfully revealed. In parallel, we quantified the ratio between the measured surface drift currents Ud, induced by mechanical waves in water, and the expected Stokes drift USt(0), provided by Eq. (3). Obtaining a solution for the first task is the main objective of our experiment, while the second task is an auxiliary consideration following from the limitations of the studies listed above. On this way, we encountered some methodological difficulties related to measuring the surface drift currents with floats, and these issues are worthy of discussion.
The layout of the paper is as follows. Section 2 describes the equipment and measurement methods used. The results and analysis for mechanical waves and wind waves are presented in sections 3 and 4, respectively. Section 5 contains a general discussion, and the conclusions are given in section 6.
2. Equipment and measurement description
In our experiment, we used the wind-wave tank of the First Institute of Oceanography of the Ministry of Natural Resources (FIO) located in Qingdao, China. A sketch of the tank is shown in Fig. 1.

A sketch of the tank. WGs are the wave gauges, PT are the Pitot tubes, ADV is the set of the acoustic Doppler velocimeters, and IP is the set for injecting in water float particles or a color dye. The width of the tank is 1 m.
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1

A sketch of the tank. WGs are the wave gauges, PT are the Pitot tubes, ADV is the set of the acoustic Doppler velocimeters, and IP is the set for injecting in water float particles or a color dye. The width of the tank is 1 m.
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
A sketch of the tank. WGs are the wave gauges, PT are the Pitot tubes, ADV is the set of the acoustic Doppler velocimeters, and IP is the set for injecting in water float particles or a color dye. The width of the tank is 1 m.
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
For the wave measurements, capacity wave gauges (WGs) were used. Pitot tubes (PTs) and three acoustic Doppler velocimeters (ADV) were used to measure the wind profile and current fluctuations, respectively. The dimensions of the tank, water depth, and air space, as well as the equipment locations, are given in Fig. 1. The locations of WG1–WG4 are denoted below as points P1–P4.
The wind is provided by the fan located at the outlet of the tank. Closed air circulation is designed (not shown). The PT set has five channels located at horizons of 7, 10, 14, 25, and 40 cm. For winds greater than 10 m s−1, all the PTs were located 1.2 cm above these levels.
Two kinds of mechanical waves were studied: (i) regular quasi-monochromatic waves and (ii) stochastic waves with a wide spectral band similar to those of PM or JW wave spectra (Komen et al. 1994). These waves were generated by a wavemaker with three initial dominant frequencies f0 of 1.5, 1.0, and 0.7 Hz and 3–5 amplitudes at each initial frequency, ranging from HS = 1 cm to HS = 20 cm (here, HS is the programmed significant wave height for generated waves). All these parameters were established in the wavemaker computer.
As another type of wave, wind waves were generated by the fan. The fan wind W is approximately equal to the mean wind speed at the middle of the air part of the tank (i.e., at z = 40 cm). Five values of wind W were applied: 4, 6, 8, 10, and 12 m s−1.
The wave records were 10 min long with a sampling frequency of 50 Hz, and each run had different wave or wind parameters. The same duration was used for the wind measurements with PT, for which the record averaging over 1-min intervals was performed. The measurements of the surface currents were obtained during these 10 min for each run. The 10-min measurements of the current fluctuations in the water body (at three horizons: −10, −20, and −30 cm) were made by an ADV with a sampling of 100 Hz. The description of these fluctuations is beyond the scope of this paper and will be presented in a separate work.
All the measurements were made 5 min after switching on the wavemaker or the fan to allow for system stabilization. In total, 42 runs were executed.
The surface drift current was measured based on the visual timing of the surface floats passing through the fixed length L. The surface floats were spherical plastic balls with bright colors (red, orange, and yellow, to distinguish them on the water surface) and two different diameters of 3 and 5 mm. The specific gravity of the floats was approximately 500 kg m−3. The fixed length L varied from 1 to 5 m based on the expected drift velocity. For each run, the floats’ timing was recorded 3–7 times until getting a minimal scattering of the results. An acceptable scattering range is approximately 5%–10% of the mean timing, which defines the empirical accuracy of the drift-current measurements.
For mechanical waves, the drift measurements were made at two drift-starting points: P1 = 9 m and P2 = 12 m. It was assumed that the drift current features for mechanical waves depend weakly on location. The most intensive measurements were made at point P2; the measurements at point P1 were mainly made to check the homogeneity of the current field. For wind waves, intensive measurements were made at three points: P1, P2, and P4 = 22 m. Small shifts in the locations of P1, P2, and P4 from the actual locations of WGs occurred due to technical reasons.
The floats were injected to the tank space with special equipment, IP, which included a 3-mm aluminum tube that was vertically fixed at the roof of the tank and continued to the level of 10 cm above the still water surface. Inside this tube, there was a 1 mm (capillary) flexible tube that was used to inject dye into the water at a fixed depth to study the wave-induced turbulent diffusion (this topic is beyond the scope of this paper). For each set of drift measurements, 2–3 plastic floats were simultaneously added to the tank at a point located 5–10 cm in front of the drift-starting point for each run. Several technical problems related to using spherical floats in this task will be discussed in section 5. The processing of the drift current data was performed without any software, as no records were collected with the equipment used.
The wave and wind data were processed with the MATLAB software. For the wave spectra evaluation, the autoregression (AR) and Welch methods were used. Because the data series had 30 000 points (more than 1000 dominant periods), the 95% confidence intervals were very small (approximately 1500 degrees of freedom). In the bilogarithmic coordinates, they were [+10%, −12%] and [+15%, −20%] for the AR method and the Welch method, respectively. The frequency spectra S(f), shown in this paper, were estimated with the AR method and are presented based on the circular frequency f given in hertz (instead of the theoretical giving spectrum in the radian frequency, ω = 2πf) (shown in Fig. 2).

Spectra at point P2 for waves with Hs = 3 cm, fp =1.5 Hz: (a) regular waves and (b) irregular waves. The 95% confidence intervals, shown by small bars, are on the order of [+10%, −12%].
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1

Spectra at point P2 for waves with Hs = 3 cm, fp =1.5 Hz: (a) regular waves and (b) irregular waves. The 95% confidence intervals, shown by small bars, are on the order of [+10%, −12%].
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
Spectra at point P2 for waves with Hs = 3 cm, fp =1.5 Hz: (a) regular waves and (b) irregular waves. The 95% confidence intervals, shown by small bars, are on the order of [+10%, −12%].
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
The intensity of wave breaking (in percent) was visually fixed based on the relative numbering of the broken wave crests with respect to the total number of dominant crests. This numbering process was executed by two persons in parallel watching the dominant crests moving between points P2 = 12 m and P4 = 22 m during one minute of observation just before the main measurements were made. The intensity of the microbreaking of small-size crests was not fixed.
3. Mechanical waves: Results and analysis
Regular (narrowband) and irregular (wideband) mechanical waves were considered separately. The results for the first kind of wave are described in detail, and these details are not repeated for the second kind of wave. The analysis of all the results is combined.
a. Results for regular waves
The results are summarized in Table 1. First, we should note a reasonable (up to 30%) difference between the estimations of the surface Stokes drift USt found with Eqs. (2) and (3) (hereafter, the notation of the surface level, z = 0, is omitted for simplicity). An example of the relation between these estimates of USt for regular waves is shown in Fig. 3a, and an example of the spectrum for narrowband waves, used to get the above estimations, is shown in Fig. 2a. Importance of using Eq. (3) is evident. Notably, in the integral of Eq. (3), we restricted the upper boundary at fmax = 15 Hz due to the accuracy restrictions of the WGs with a 1-mm wave height. The difference between the mentioned estimates suggests that Eq. (3) is the more accurate equation. Hereafter, only such estimations of USt are presented.
Parameters of waves and surface currents for regular waves. Note that the two values for the Stokes drift, separated by a slash, correspond to the estimations with Eqs. (2) and (3), respectively; “n-dim” indicates nondimensional.



Relation between the Stokes drift estimations with Eq. (2), Ust2, and with Eq. (3), Ust3. (a) Regular mechanical waves (Table 1): the solid trend line (y1 = 1.30x1) corresponds to fp = 1.5 Hz; the dashed trend line (y2 = 1.27x2) corresponds to fp = 1.0 Hz; the dotted trend line (y3 = 1.18x3) corresponds to fp = 0.7 Hz. (b) Irregular waves with PM spectra (Table 2): the solid trend line (y1 = 1.98x1) corresponds to fp = 1.5 Hz; the dashed trend line (y2 = 2.73x2) corresponds to fp = 1.0 Hz.
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1

Relation between the Stokes drift estimations with Eq. (2), Ust2, and with Eq. (3), Ust3. (a) Regular mechanical waves (Table 1): the solid trend line (y1 = 1.30x1) corresponds to fp = 1.5 Hz; the dashed trend line (y2 = 1.27x2) corresponds to fp = 1.0 Hz; the dotted trend line (y3 = 1.18x3) corresponds to fp = 0.7 Hz. (b) Irregular waves with PM spectra (Table 2): the solid trend line (y1 = 1.98x1) corresponds to fp = 1.5 Hz; the dashed trend line (y2 = 2.73x2) corresponds to fp = 1.0 Hz.
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
Relation between the Stokes drift estimations with Eq. (2), Ust2, and with Eq. (3), Ust3. (a) Regular mechanical waves (Table 1): the solid trend line (y1 = 1.30x1) corresponds to fp = 1.5 Hz; the dashed trend line (y2 = 1.27x2) corresponds to fp = 1.0 Hz; the dotted trend line (y3 = 1.18x3) corresponds to fp = 0.7 Hz. (b) Irregular waves with PM spectra (Table 2): the solid trend line (y1 = 1.98x1) corresponds to fp = 1.5 Hz; the dashed trend line (y2 = 2.73x2) corresponds to fp = 1.0 Hz.
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
The second peculiarity of the measurements is that the ratio Ud/USt varies in the range of 0.5–0.8 (the second-to-last column in Table 1), reflecting a permanent underestimation of the expected surface Stokes drift. Specifically, there is no evident dependence of this ratio on the intensity of wave breaking (presented in the rightmost column in Table 1). Following the review by Bremer and Breivik (2017), the estimates of Ud/USt for the tank measurements are pioneering, despite numerous empirical reports of the abovementioned discrepancies between measurements and theoretical values. A more detailed analysis will be made in section 3c.
b. Results for irregular waves
These results are summarized in Table 2. In this case, the peculiarities are similar to those in the previous case. The role of spectral shape is well shown in Figs. 2a,b and 3a,b. From these figures, one can see that for a wideband spectrum, the estimates of Ust based on Eqs. (2) and (3) can be more than 2 times different.
Parameters of waves and surface currents for irregular waves; “n-dim” indicates nondimensional.


The additional difference in the surface drift for irregular waves is a slightly greater ratio Ud/USt varying from 0.6 to 0.9. In this case, the integral in Eq. (3) should be greater due to the consideration of a wider spectral band than in the previous case (Fig. 2b); therefore, the only explanation for the higher values of Ud/USt is accounting for the small steepness Ϭ of irregular waves (see Table 2). Theoretically, it is reasonable to assume that for smaller wave steepness values, the observations of Ud should be closer to the theoretical values (Stokes 1847).
c. Combined analysis of mechanical waves
The presented results are novel, despite numerous empirical reports on the discrepancies between measurements and theoretical values of Stokes drift [e.g., Craik (1982), Breivik et al. (2014), Paprota et al. (2016), and so on, to the review by Bremer and Breivik (2017)]. Notably, for the first time, we obtained empirically quantified ratios Ud/USt for regular and irregular mechanical waves in a tank. These ratios are interesting from three points of view: (i) understanding the reasons for the empirical underestimation of the surface drift currents, (ii) methodical points of measurements, and (iii) identifying possible applications of these results.
Underestimating Ud is not new; nevertheless, this issue requires separate consideration. The following reasons for the empirical underestimation of surface drift currents in experiments are well known: (i) the viscosity of the water (which disrupts the potentiality of wave motion), (ii) return currents, (iii) water surface contamination, and (iv) the sailing of floats, including dynamic-pressure braking (McNally and White 1985). The first three reasons were discussed in previous studies (e.g., Longuet-Higgins 1953; Craik 1982; Herterich and Hasselmann 1982; see also Breivik et al. 2014; Paprota et al. 2016). We cannot add more to this discussion, noting that the first two reasons are inevitable and cannot be regulated; however, the third influence can be minimized by cleaning the water surface after each change in frequency, thus the third item should be excluded from the list of reasons discussed. Therefore, we do not dwell on the first three reasons, but the phenomenon of the float sailing is worthy of discussion.
Besides the return current influence (Monismith et al. 2007), which cannot be quantified in our case, in our mind, the sailing of floats on the water surface is the main reason for the abovementioned discrepancy between the experimental and theoretical values of the Stokes drift. To validate this assumption, we tried to regulate the effect of sailing by choosing different sizes of floats. Indeed, when we used plastic balls with a diameter of 5 mm as the floats, the drift current values were fixed at values in 20%–30% smaller than those for floats that were 3 mm in diameter. Thereby, the assumption stated above is justified.
The nature of the sailing process is relatively evident. In the case of mechanical waves and spherical floats, the sailing has two negative constituents, namely, the air and the water components. The air component is associated with the air-induced resistance to float propulsion due to the viscosity and negative dynamic pressure of the ambient motionless air. On the other hand, the water sailing is influenced by the same manner due to the rapid decrease in the total drift current with depth, described by profile Ud(z). For mechanical waves, this decrease follows from Eq. (3). In the case of wind waves, the rapid decrease in the wind-induced current Udw(z) has been empirically shown by Longo et al. (2012) and described in Polnikov (2011, 2019).
These considerations account for mismatch between theoretical and experimental values, and one may go on to methodical points of the surface-drift current measurements.
From a methodological perspective, a number of previous experiments (e.g., Wu 1975, 1983; Zavadsky and Shemer 2017), in which the volumetric surface floats were used, have suffered from the abovementioned problems, resulting in the underestimation of surface drift in these experiments. Therefore, it is important to make certain methodological changes when obtaining future measurements. To avoid the sailing of floats, it is reasonable to use floats that are fully flat. For example, floats could be pieces of thin flexible tissue (similar to gauze) with a net density less than that of water. Such experiments have been performed in the past (see references in Tsahalis 1979), and they should be planned in the future to further improvement the results presented above.
Another alternative approach involves using PIV methods (e.g., Longo et al. 2012; Paprota et al. 2016) or the laser-warming method (thermal-marking velocimetry) (Savelyev et al. 2012) to measure surface drift currents. However, due to technical reasons related to the optical equipment location, such methods are mainly appropriate for small tanks. But in such cases, the tank boundaries could result in additional sources of errors. Another problem associated with using the optical methods to obtain current measurements in tanks is the restriction to equipment location and movement along a tank. Ideally, technical progress will help overcome the mentioned methodological problems.
In practical applications, the obtained results can be used as the first guess for expected drift current estimates in the engineering of harbors and other marine infrastructure, taking in mind the influence by swell. On other hand, from the physical point of view, the presented results are useful for accurately estimating the observed wind-induced drift currents Udw, based on Eq. (4). It is this idea (rather than simple estimating ratio Ud/USt) was stimulating our measurements of surface drift currents Ud for mechanical waves, the results of which will be actively used below in the data analysis for wind waves.
4. Wind waves: Results and analysis
The results for wind waves are summarized in Table 3, where the following notations are used: Ϭ = akp is the mean wave steepness,
Parameters of wind, waves and surface currents for wind waves. Note that the ↔ sign separates the empirical and parametrical values


a. Analysis of the accuracy of the results
First, let us estimate the accuracy of the results starting with the wind parameter
The wave parameters ωp = 2πfp, Ϭ, and A were estimated with very high accuracy (errors less than 1%) from the wave spectra calculated within the MATLAB software (see section 2). For the completeness of the wind-wave spectra description, their spatial evolution is presented in Fig. 4. Due to a wide spectral band for wind waves, using Eq. (3) is quite important for estimation of theoretical value for the Stokes drift. In our case, these values have errors of approximately 5%–7%, what is defined mainly by the value of the upper limit of the integral in Eq. (3). It was chosen to be 15 Hz due to the technical restrictions of the WGs used. As shown in Fig. 4, the wave spectra are radically changing along the tank (this feature was ignored for mechanical waves), therefore the theoretical estimates of USt are also varying, which plays a crucial role in the correct estimating the wind-induced drift currents Udw based on using Eq. (3).

Wind-wave spectra at points P1, P2, and P4 for the run with W = 8 m s−1. The straight line with digits shows the frequency-decay law of spectra in the tail range (the error is ~ 3%).
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1

Wind-wave spectra at points P1, P2, and P4 for the run with W = 8 m s−1. The straight line with digits shows the frequency-decay law of spectra in the tail range (the error is ~ 3%).
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
Wind-wave spectra at points P1, P2, and P4 for the run with W = 8 m s−1. The straight line with digits shows the frequency-decay law of spectra in the tail range (the error is ~ 3%).
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
The visually obtained float movement timing resulted in the error on the order of 10% for the measured total mean drift Ud. In turn, the intensity of wave breaking was estimated (by two persons in parallel) with the same error. These values influence the accuracy of the empirical results for the wind-induced drift currents Udw.
b. Analysis of the main results
In the case of wind waves, we estimated the wind-induced drift currents as the difference between the measured value Ud and the theoretically estimated value USt by assuming that the positive floats sailing in the air (due to the wind) is compensated for by the negative floats sailing in the water. This assumption was verified by the measurements with plastic balls of different diameters (3 and 5 mm) used as the floats. This means that the values of Udw presented in Table 3 are the smallest possible values, taking in mind the results for mechanical waves discussed above in section 3a.
Based on Table 3, the determined empirical ratio
According to the theoretical model of the wind-induced drift currents proposed in a previous study (Polnikov 2019), one may expect
In a graphic form, the relation between the empirical and parametrical values of

Correspondence between empirical and parametrical values of ratio
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1

Correspondence between empirical and parametrical values of ratio
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
Correspondence between empirical and parametrical values of ratio
Citation: Journal of Physical Oceanography 50, 10; 10.1175/JPO-D-20-0009.1
Equation (8) provides the dependence relation between the wind-induced surface drift currents and the wind and wave parameters, which is the main aim of this paper. This dependence is established for the first time during a long history of empirical investigations (section 1a). Below, it will be shown that ratio (8) corresponds relatively well to the results published earlier for the tank measurements (Wu 1975; Longo 2012; Longo et al. 2012; Zavadsky and Shemer 2017). In the future, parameterization (8) could be checked for the field data based on accurate measurements of wind and wave parameters.
5. Discussion
In this section, we discuss the following points: (i) the impact of tank dimensions and wind variations along a tank on the accuracy of estimates of
First, we should mention that the small dimensions of tanks used in some studies (Wu 1975; Longo et al. 2012; Zavadsky and Shemer 2017) could reasonably restrict the accuracy of the wind and wave measurements. Evidently, in such cases, the upper and lateral boundaries of tanks are the sources of errors in the wind profile measurements and, consequently, the estimations of
Regarding the Stokes drift estimations, in the presence of mechanical waves, the current measurements inevitably underestimate the wave-induced drift at the water surface with respect to the theoretical value USt. In turn, as shown in section 3, the theoretical estimates of USt with Eq. (2) that are typically used in practice (e.g., Toba 1973; Wu 1975; Longo 2012; Longo et al. 2012) are incorrect. Therefore, the reported ratio between USt and Udw was underestimated, the value of which was fixed at approximately 15% (Toba 1973; Wu 1975; Longo et al. 2012). The more accurate estimates, especially in the presence of wind waves with wideband spectra, can be obtained with using Eq. (3). In such a case, as shown in Table 3, the ratio USt/Udw could reach values of 35%–40%. This circumstance is very important from a practical perspective. Of course, all these results should be checked again in future experiments conducted with highly accurate measurements.
Regarding the wind-induced drift parameterization, the general assumption is of the form
where factors Cb, Cw, and power parameters k, m, and n are the fitting dimensionless constants. Based on the trial-and-error method, we found that the simplest version of Eq. (9) has the form of Eq. (8), when all the parameters are equal to unity. According to Table 3, Eq. (8) provides a relatively good approximation, with a mean error for the value
At meanwhile, it is interesting to apply parameterization (8) to the other known results. For example, for results of the classical work by Wu (1975), we obtained the following findings. From Figs. 2 and 3 in Wu (1975), typical values for the wave parameters are a = 0.02–0.03 m, ωp = 12.5 s−1, while
Addressing the data from Longo (2012) and Longo et al. (2012), we obtained the following findings. From Tables 1 and 3 of Longo (2012), the typical wave parameters are a ≈ 0.002–0.003 m and ωp = 30–40 s−1, which resulting in Ϭ ≈ 0.2, for a tunnel wind of W = 10.9 m s−1. From Table 1 of Longo et al. (2012), in this case,
In Zavadsky and Shemer (2017), the analogous data for wind W = 10.5 m s−1 are as follows: a ≈ 0.01 m, ωp ≈ 15 s−1, and
As one can see, parameterization (8) works well in all these cases. Similar examples can be explored for the other experimental works mentioned above if the proper data are supplied.
Finally, the methodology should be discussed. Naturally, using hard surface floats to measure surface drift currents seems obsolete. The more advanced methods are PIV (e.g., Longo 2012; Longo et al. 2012; Paprota et al. 2016) or laser-warming methods (Savelyev et al. 2012), although they cannot be used in large tanks due to the need to apply many optical devices through transparent tank boundaries and restrictions related to moving optical equipment throughout the tank. Most laboratories have no such equipment at their disposal. Nevertheless, traditional floats can still provide reasonable results.
The main problem, in this case, is the sailing of floats, as discussed in detail in section 3c. This problem could be overcome if one uses appropriate floats. The solution is to use flat floats in the form of pieces of thin flexible tissue (like gauze) with a net density less than that of water. If the width of the tissue is 0.1 mm, the sailing of the floats will be negligibly small; thus, they would not be disturbed by the wind or the vertical shear of drift currents. In such a case, the results of the drift measurements would be more accurate than in most other cases, what suggests that one could plan similar experiments in the near future. Thus, the best choice is to use tanks with lengths greater than 30 m, widths greater than 1 m and air part of approximately 1 m or more to ensure correct estimations of the friction velocity
6. Conclusions
Finally, we state the following conclusions.
First, the surface drift currents Ud induced by mechanical waves in a tank were measured with plastic surface floats and compared with the Stokes drift at the surface USt estimated by Eq. (3). The ratio Ud/USt varies in the range of 0.5–0.93 and slightly increases with decreasing wave steepness, with no visible dependence on the breaking intensity. The reduction in empirical values of Ud with respect to the theoretical value USt is mainly due to the following reasons: return currents, fluid viscosity, and the sailing of floats in the air and in the water. The significant effect of the latter reason could be minimized by using flat and flexible tissue floats or by employing modern optical methods to collect surface current measurements.
Second, in the case of wind waves, the wind-induced part of the surface drift Udw was estimated by subtracting the theoretical value for the Stokes drift USt found with Eq. (3) from the measured total drift Ud at three locations along the tank for a series of wind speeds. These values of Udw were compared with the friction velocity
Third, the obtained ratio
Acknowledgments
The authors are grateful to PhD student Shumin Jiang and graduate students Xue Wang and Chao Li for their help with measurements and data processing. The work was supported by the National Key Research and Development Program of China under Grant 2017YFC1404005, the Russian Foundation for Basic Research Grant 19-04-00161, and the Natural Science Foundation of China under Grants 41821004 and U160640.
Data availability statement
The raw measurement data are available from the corresponding author (V. Polnikov).
REFERENCES
Anis, A., and J. N. Moum, 1995: Surface wave–turbulence interactions: Scaling ε(z) near the sea surface. J. Phys. Oceanogr., 25, 2025–2045, https://doi.org/10.1175/1520-0485(1995)025<2025:SWISNT>2.0.CO;2.
Babanin, A. V., 1988: Relationship between the parameters of the surface drift current and the determining wind in the northwestern part of the Black Sea (in Russian). Mar. Hydrophys. J., 4, 55–58.
Babanin, A. V., and B. K. Haus, 2009: On the existence of water turbulence induced by non-breaking surface waves. J. Phys. Oceanogr., 39, 2675–2679, https://doi.org/10.1175/2009JPO4202.1.
Babanin, A. V., and D. Chalikov, 2012: Numerical investigation of turbulence generation in non-breaking potential waves. J. Geophys. Res., 117, C00J17, https://doi.org/10.1029/2012JC007929.
Benilov, A. Y., 2012: On the turbulence generated by the potential surface waves. J. Geophys. Res., 117, C00J30, https://doi.org/10.1029/2012JC007948.
Breivik, O., J. Bidlot, and P. A. E. M. Janssen, 2014: Approximate Stokes drift profiles in deep water. J. Phys. Oceanogr., 44, 2433–2445, https://doi.org/10.1175/JPO-D-14-0020.1.
Bremer, T. S., and Ø. Breivik, 2017: Stokes drift. Philos. Trans. Roy. Soc., 376A, 20170104, http://doi.org/10.1098/rsta.2017.0104.
Bühler, O., 2014: Waves and Mean Flows. 2nd ed. Cambridge University Press, 360 pp.
Bye, J., 1988: The coupling of wave drift and wind velocity profiles. J. Mar. Res., 46, 457–472, https://doi.org/10.1357/002224088785113559.
Chalikov, D., and S. Rainchik, 2011: Coupled numerical modelling of wind and waves and theory of the wave boundary layer. Bound.-Layer Meteor., 138, 1–41, https://doi.org/10.1007/s10546-010-9543-7.
Churchill, J. H., and G. T. Csanady, 1983: Near-surface measurements of quasi-Lagrangian velocities in open water. J. Phys. Oceanogr., 13, 1669–1680, https://doi.org/10.1175/1520-0485(1983)013<1669:NSMOQL>2.0.CO;2.
Clarke, A. J., and S. V. Gorder, 2018: The relationship of near-surface flow, Stokes drift and the wind stress. J. Geophys. Res. Oceans, 123, 4680–4692, https://doi.org/10.1029/2018JC014102.
Constantin, A., 2006: The trajectories of particles in Stokes waves. Invent. Math., 166, 523–535, https://doi.org/10.1007/s00222-006-0002-5.
Craik, A. D. D., 1982: The drift velocity of water waves. J. Fluid Mech., 116, 187–205, https://doi.org/10.1017/S0022112082000421.
Herterich, K., and K. Hasselmann, 1982: The horizontal diffusion of tracers by surface waves. J. Phys. Oceanogr., 12, 704–711, https://doi.org/10.1175/1520-0485(1982)012<0704:THDOTB>2.0.CO;2.
Huang, N. E., 1979: On surface drift currents in the ocean. J. Fluid Mech., 91, 191–208, https://doi.org/10.1017/S0022112079000112.
Janssen, P. A. E. M., 2004: The Interaction of Ocean Waves and Wind. Cambridge University Press, 308 pp.
Kitaigorodskii, S., M. A. Donelan, J. L. Lumley, and E. A. Terray, 1983: Wave–turbulence interactions in the upper ocean. Part II: Statistical characteristics of wave and turbulent components of the random velocity field in the marine surface layer. J. Phys. Oceanogr., 13, 1988–1999, https://doi.org/10.1175/1520-0485(1983)013<1988:WTIITU>2.0.CO;2.
Komen, G. I., L. Cavaleri, M. Donelan, K. Hasselmann, S. Hasselmann, and P. A. E. M. Janssen, 1994: Dynamics and Modelling of Ocean Waves. Cambridge University Press, 554 pp.
Kudryavtsev, V. N., V. Shrira, V. Dulov, and V. V. Malinovsky, 2008: On the vertical structure of wind-driven sea currents. J. Phys. Oceanogr., 38, 2121–2144, https://doi.org/10.1175/2008JPO3883.1.
Longo, S., 2012: Wind-generated water waves in a wind tunnel: Free surface statistics, wind friction and mean air flow properties. Coastal Eng., 61, 27–41, https://doi.org/10.1016/j.coastaleng.2011.11.008.
Longo, S., L. Chiapponi, M. Clavero, T. Mäkel, and D. Liang, 2012: Study of the turbulence in the air-side and the water-side boundary layers in experimental laboratory wind induced surface waves. Coastal Eng., 69, 67–81, https://doi.org/10.1016/j.coastaleng.2012.05.012.
Longuet-Higgins, M., 1953: Mass transport in water waves. Philos. Trans. Roy. Soc. London, 245A, 535–581, https://doi.org/10.1098/rsta.1953.0006.
Lumley, J. L., and E. A. Terray, 1983: Kinematics of turbulence convected by a random wave field. J. Phys. Oceanogr., 13, 2000–2007, https://doi.org/10.1175/1520-0485(1983)013<2000:KOTCBA>2.0.CO;2.
Malinovsky, V. V., V. A. Dulov, A. E. Korinenko, A. N. Bol’shakov, and V. E. Smolov, 2007: Field investigations of the drift of artificial thin films on the sea surface. Izv. Atmos. Oceanic Phys., 43, 103–111, https://doi.org/10.1134/S0001433807010124.
McNally, G. J., and W. B. White, 1985: Wind driven flow in the mixed layer observed by drifting buoys during autumn-winter in the midlatitude North Pacific. J. Phys. Oceanogr., 15, 684–694, https://doi.org/10.1175/1520-0485(1985)015<0684:WDFITM>2.0.CO;2.
Mellor, G. L., 2001: One-dimensional, ocean surface layer modelling: A problem and a solution. J. Phys. Oceanogr., 31, 790–809, https://doi.org/10.1175/1520-0485(2001)031<0790:ODOSLM>2.0.CO;2.
Monismith, S. G., E. A. Cowen, H. M. Nepf, J. Magnaudet, and L. Thais, 2007: Laboratory observations of mean flows under surface gravity waves. J. Fluid Mech., 573, 131–147, https://doi.org/10.1017/S0022112006003594.
Paprota, M., W. Sulisz, and A. Reda, 2016: Experimental study of wave-induced mass transport. J. Hydraul. Res., 54, 423–434, https://doi.org/10.1080/00221686.2016.1168490.
Phillips, O. M., 1977: The Dynamics of the Upper Ocean. 2nd ed. Cambridge University Press, 336 pp.
Polnikov, V. G., 2011: Integrated model for a wave boundary layer. Mar. Sci., 1, 10–21, https://doi.org/10.5923/j.ms.20110101.02.
Polnikov, V. G., 2019: A semi-phenomenological model for wind-induced drift currents. Bound.-Layer Meteor., 172, 417–433, https://doi.org/10.1007/s10546-019-00456-1.
Rascle, N., and F. Ardhuin, 2009: Drift and mixing under the ocean surface revisited: Stratified conditions and model-data comparisons. J. Geophys. Res., 114, C02016, https://doi.org/10.1029/2007JC004466.
Savelyev, I. B., E. Maxeiner, and D. Chalikov, 2012: Turbulence production by nonbreaking waves: Laboratory and numerical simulations. J. Geophys. Res., 117, C00J13, https://doi.org/10.1029/2012JC007928.
Shemdin, O. H., 1972: Wind-generated current and phase speed of wind waves. J. Phys. Oceanogr., 2, 411–419, https://doi.org/10.1175/1520-0485(1972)002<0411:WGCAPS>2.0.CO;2.
Stokes, G. G., 1847: On the theory of oscillatory waves. Trans. Cambridge Philos. Soc., 8, 441–455.
Sullivan, P. P., and J. C. McWilliams, 2010: Dynamics of winds and currents coupled to surface waves. Annu. Rev. Fluid Mech., 42, 19–42, https://doi.org/10.1146/annurev-fluid-121108-145541.
Teixeira, M. A. C., 2018: A model for the wind-driven current in the wavy oceanic surface layer: Apparent friction velocity reduction and roughness length enhancement. J. Phys. Oceanogr., 48, 2721–2736, https://doi.org/10.1175/JPO-D-18-0086.1.
Teixeira, M. A. C., and S. E. Belcher, 2002: On the distortion of turbulence by a progressive surface wave. J. Fluid Mech., 458, 229–267, https://doi.org/10.1017/S0022112002007838.
Toba, Y., 1973: Local balance in the air–sea boundary processes. III: On the spectrum of wind waves. J. Oceanogr. Soc. Japan, 29, 209–220, https://doi.org/10.1007/BF02108528.
Troitskaya, Y. I., D. A. Sergeev, A. A. Kandaurov, G. A. Baidakov, M. A. Vdovin, and V. I. Kazakov, 2012: Laboratory and theoretical modeling of air-sea momentum transfer under severe wind conditions. J. Geophys. Res., 117, C00J21, https://doi.org/10.1029/2011JC007778.
Tsahalis, D. T., 1979: Theoretical and experimental study of wind- and wave-induced drift. J. Phys. Oceanogr., 9, 1243–1257, https://doi.org/10.1175/1520-0485(1979)009<1243:TAESOW>2.0.CO;2.
Wu, J., 1975: Wind-induced drift currents. J. Fluid Mech., 68, 49–70, https://doi.org/10.1017/S0022112075000687.
Wu, J., 1983: Sea-surface drift currents induced by wind and waves. J. Phys. Oceanogr., 13, 1441–1451, https://doi.org/10.1175/1520-0485(1983)013<1441:SSDCIB>2.0.CO;2.
Zavadsky, A., and L. Shemer, 2017: Water waves excited by near-impulsive wind forcing. J. Fluid Mech., 828, 459–495, https://doi.org/10.1017/jfm.2017.521.