The Response of the Water Surface Layer to Internal Turbulence and Surface Forcing

: We have carried out an experimental study of the turbulence kinetic energy dissipation rate ( (cid:101) ), temperature dissipation rate ( χ ), and turbulent heat ﬂux (THF) within the water surface layer in the presence of non-breaking wave, surface convection, and horizontal heat and eddy ﬂuxes that play a prominent role in the front. We noted that the non-breaking wave dominates (cid:101) values within the surface layer. While analyzing the vertical (cid:101) variability, the presence of a wave-affected layer from the water surface to a depth of z ≈ 1.25 λ w is observed, where λ w is the wavelength. (cid:101) associated with non-breaking waves ranged to 4.9 × 10 − 6 –7 × 10 − 6 m 2 /s 3 for the wavelength range of 0.038 m < λ w < 0.098 m categorized as the gravity and gravity-capillary wave regimes. (cid:101) values increase for longer λ w and non-breaking wave (cid:101) values represent their signiﬁcant contribution to the ocean energy budget and dynamic of surface layer considering that the non-breaking wave covers the large fraction of ocean surface. We also found that the surface mean square slope (MSS) and wave generated (cid:101) have the same order of magnitude, i.e., MSS ∼ (cid:101) . Besides, we have documented that the small-scale temperature ﬂuctuation change (i.e., χ ) is consistent with the large-scale temperature gradient change (i.e., d < T > / dz ). The value of the THF is approximately constant within the surface layer. It represents that the measured THF near the water surface can be considered a surface water THF, challenging to measure directly.


Introduction
The ocean covers approximately 71 percent of the Earth's surface [1]. At the air-sea interface, the exchange of heat, moisture, momentum, and gas transfer is carried out by molecular transfer processes; however, at a depth greater than ∼1 mm, the turbulence dominates the mixing processes [2].
value has implications on heat flux across the ocean interface [3], the air-sea gas fluxes velocity, K g ∼ 1/4 [4] and it is essential for the oil industry [5] due to that the higher turbulence enhances the biodegradation of oil [6]. Understanding distribution and its role in near-surface water and ocean mixing are vital in studying ocean processes and the upper ocean boundary layer dynamic.
The value of can be then found from the relationship = 2υ < s ij s ij > [7], where s ij is the fluctuating rate of strain: s ij = 1 2 ( Here, υ is the kinematic molecular viscosity, i and j = 1, 2, and 3, and correspond to x, y and z coordinates, respectively. u i represents the fluctuating part of the velocity components, and < > is the ensemble average. The value of χ comes from the equation χ = 2K θ < ( ∂θ ∂x i ) 2 >= 2K θ < ( ∂θ ∂x ) 2 + ( ∂θ ∂y ) 2 + ( ∂θ ∂z ) 2 >, where K θ is the thermal molecular diffusivity, and θ is the fluctuating part of temperature. The vertical and χ variability have been investigated within the upper ocean boundary layer while the boundary layer is influenced by one of the surface forcing such as heat flux [4,8], internal sources of turbulence [9], surface gravity wave [8,[10][11][12][13], and capillary wave [14]. Fredriksson et al. [4] simulated a numerical model to calculate for a free surface flow driven by natural convection. They found that the oceanic free convection results in a sharp change of beneath the water surface. Wuest et al. [15] measured in the wind-forced stratified water and observed that ≈90% of turbulent kinetic energy was dissipated within the upper boundary layer. Terray et al. [16] investigated under breaking waves and observed a large uniform from the surface water to a depth of z = 0.6H s , where H s is the significant wave height.
Given that the percentage of the ocean surface covered by wave breaking under strong wind is less than 10% [17], it indicates the significant role of the non-breaking wave on the ocean budget. Babanin and Haus [18] conducted a laboratory experiment to measure beneath monochromatic non-breaking waves, which showed the presence of . Bogucki et al. [19] observed that associated with non-breaking solitary waves ranged to 3 × 10 −4 m 2 /s 3 for a wave amplitude of 50 cm.
The scarcity of field data hinders and χ universal parametrization within the upper ocean boundary layer, especially very near the surface. This paper presents laboratory experiments in the Air-Sea Interaction Saltwater Tank (ASIST) at the University of Miami. We try to simply simulate the ocean when the horizontal heat and eddy fluxes play a prominent role in the ocean, like the front in the Gulf of Mexico where the Mississippi River with massive horizontal eddy fluxes reach the ocean and generate the Front [13,20]. We investigate the THF, χ, and variability very close to the air-sea interface, 0.5 cm beneath the water surface, when subject to the internal turbulence, surface convection, and nonbreaking wave categorized as the gravity and gravity-capillary wave regimes [21][22][23]. The sources of turbulence and temperature flux in our experiment and the experimental setup are addressed in Section 2. Section 3 describes the experimental results of , χ, and THF. Finally, a conclusion is given in Section 4.

Experimental Setup, Data Acquisition, and Analysis
We simulated the oceanic-like forcing in our laboratory experiment by having the three turbulence sources, i.e., the internal sources of turbulence, surface convection, and nonbreaking surface waves.

Experimental Setup
The experiments were conducted in the ASIST tank at the University of Miami, equipped with the turbulence generating grid and a heated grid ( Figure 1). The tank walls are constructed of acrylic panels with a thickness of 0.024 m and have dimensions 15 m long, 1 m wide, and 1 m high. In our experiment, we analyzed data for the mean water flow velocity of < U >= 0.066 m/s, 0.125 m/s, 0.167 m/s, and 0.183 m/s. The heated grid was located at the tank entrance, and the turbulence-generating grid was mounted 0.5 m behind the heater. The freshwater depth during the experiment was kept constant, d t = 0.5 m. The thickness of the bar of the grid was d b = 0.02 m, and the distance between the centers of two cells on the grid horizontally and vertically was M = 0.1 m. The solidity of the grid was [24]. The grid Reynolds number in our experiment was given by Re M = (σM < U > /υ) [25] (Table 1). The minimum turbulence Reynolds number, Re λ = (< u i > rms λ)/υ [26], calculated based on the Taylor microscale, λ ≈ 2 15υ < u i > rms , was Re λ = 487 for < U >= 0.066 m/s, where < u i > rms is the root mean square of the velocity fluctuation.  The flow velocity was measured with a current meter (Infinity-EM, Model AEM-USB) with a sampling rate of 10 Hz. The Rockland Scientific vertical microstructure profiler (VMP200) [27][28][29][30] was used to measure the and χ with a sampling rate of 512 Hz. The VMP200 was equipped with two shear sensors that the sensors sampled the small-scale shear component. The shear probe, called the airfoil probe, was initially developed for the wind tunnel work [31], and Osborn [32] adopted it for oceanic measurements. The shear probe senses velocity fluctuations cross-stream to its travel direction. The VMP200 was also equipped with a fast thermistor FP07. The response time of the thermistor FP07 is 7 ms in water for the speed of 1 m/s [29]. The speed increase causes the decrease in thermistor response time. The temperature resolution of the FP07 is 0.0001 • C. The measured temperature temporal gradients ∂θ ∂t via Taylor's frozen hypothesis to ∂θ ∂x = 1 <U> ∂θ ∂t . The VMP was mounted horizontally in the tank (Figure 1), and it collected a time series of the velocity shear and temperature at a depth range of 0.5 cm < z < 25 cm. By considering the standard assumption that the oceanic flow can be approximated by idealized homogenous and isotropic turbulence, the values were calculated by the VMP200 measured shear spectrum Ψ(K) as [27,33] where K is the wavenumber. The shear probe's finite spatial size causes it to spatially average the smallest eddies for large K. The lost variance is corrected with a transfer function [28]. The shear spectrum is also fixed for the vibration-coherent portion by Goodman et al.'s technique [34].
Kolmogorov [35] derived the shear spectrum by assuming that the larges scales of turbulence are much larger than the Kolmogorov scale. Kolmogorov [35] presented that the shear spectrum is proportional to K 1/3 in the inertial subrange. Due to difficulties in resolving eddies scales smaller than the Kolmogorov length [36], the shear spectrum can be fitted with an empirical turbulence spectrum such as the Nasmyth spectrum [29,37] over the viscous and inertial subrange. The Nasmyth spectrum's integral over the wavenumber is considered as a value in this paper.
Temperature dissipation rate is estimated from temperature gradient spectra φ(K) [37,38] in one direction as χ = 6K θ ∞ 0 φ(K)dK [33]. For measuring χ values in this paper, we first find the fitting line for temperature gradient spectra along the dropping part of the spectrum. χ values were measured by multiplying 6K θ with the integral of the temperature gradient spectrum and fitting line along the wavenumber domain. A sophisticated thermistor signal processing, installed on the VMP200, minimizes the electronics noise. Therefore, the measured temperature spectra are only limited by the thermistor inertia [27]. When traveling through the water column is faster than 0.1 m/s, the thermistor (FP07) used in the VMP200 does not fully resolve the temperature variance of the temperature field (Lueck 1977), considering that we did not implement the thermistor spectral response function correction on the χ values [39]. The correction of χ values is difficult because the thermistor's response time has been found to depend on the VMP velocity [39] and thickness of the glass coating of the sensor tip [40] that varies for each individual thermistor. Nash et al. [40] observed that only 10% of the temperature gradient variance could be resolved at a profiler of 0.6 m/s for χ values larger than 1 × 10 −6 • C 2 /s. The VMP200measured χ values are the underestimates of the χ for the mean flow velocity of faster than 0.1 m/s. The imprecise-measured values present the χ behavior within the surface layer, and due to that, we present χ values in this paper.
The direction of the flow and the VMP probes generates an instantaneous angle of attack. The angle of attack was one of the most critical parameters that could affect the VMP results. The effects of this angle on the results of and χ were investigated in the ASIST tank. The shear and thermistor probe results were found to be consistent for an angle of attack <12 degree, with an error of about 10% and 20% for the shear and the thermistor probes, respectively. The angle of attack was kept less than 1 degree in experiment that gave an error of less than 2%.
The non-breaking waves were observed to be propagated and spread uniformly along the tank. We are unsure of the source of propagating waves; therefore, the generated-wave may not be entirely representative of the effects of actual ocean waves. We speculate that the surface stress created by the friction between the moving water and the stationary air contribute partially to generate a surface wave in addition to the grid. The sidewall also has effects on the wave's generation for larger flow velocity.
The laser wave slope instrument is used to measure the water surface slope, dη/dx (installed 1.35 m in front of the VMP). The η is the surface elevation that equals < η >= 0.5H, where H is wave height. The point height/slope gauge consisted of an Argon-Ion (488 nm-blue) laser transmitting 2 W of power, whose beam was directed upward through the water surface. Above the tank along the sidewall, a line-scan camera observed the surface spot and tracked the vertical movement [41]. The surface slope spectrum equals P(K) = K 2 S(K) , where S(K) is the elevation spectrum of η, and the mean square slope

Internal Source of Turbulence
The internal source of turbulence, or preexisting source of turbulence, in the ocean was simulated by the grid generated turbulence, which was created by passing water through a solid grid ( Figure 1). The values of , χ, and temperature variance, < θ 2 >, decay with distance from the grid proportional to ≈ (x) −(n+1) , χ ≈ (x) −(m+1) , and < θ 2 >≈ (x) −m , respectively [43] (see Appendix A). The power-law exponent of m and n are to be determined empirically. Antonia et al. [44] showed a value of n to be n = 1.28, and more recently Hearst and Lavoie [45] found a value of n = 1.37 and 1.39 behind a square-fractal-element grid. Warhaft and Lumley [46] found that the temperature decay rate varied over a wide range of 0.87 < m < 3.09.

Convection and Turbulent Heat Flux
In the experiment, the mean water temperature at a depth of z = 0.15 m was < T > = 26.82 • C (see Section 3.2) and the air temperature was < T > = 25.70 • C at 0.2 m above the water surface. The heat transfer from the freshwater to air resulted in thermal convection in our experiment as salinity was approximately zero within the water depth. Convection affects the vertical transport of heat, momentum, and other properties. Despite their importance to ocean circulation within the upper ocean boundary layer [47], the vertical fluxes of heat and momentum are only estimated indirectly. The turbulent heat flux (THF) was estimated as THF ≈ ρC P < w θ > [48], where ρ is water density and C p is the specific heat of water. We followed the Osborn and Cox [49] approach to determine the value of turbulent temperature flux < w θ >. For the steady and homogeneous turbulence, by retaining only vertical dependence and by neglecting the surface wave for simplicity, the turbulent temperature flux is given by [50] The values of χ and temperature gradient T z = ∂<T> ∂z were measured in our experiment; therefore, the value of the turbulent temperature flux can be estimated as < w θ >== − 1 2 χ/T z . Therefore, the THF equation is rewritten as

Turbulence Scaling
Normalizing of aids in understanding which processes have significant effects on the surface layer turbulence [4,51]. The upper ocean turbulence is predominantly generated due to atmosphere-ocean interaction by convection and surface wave [51]. The base of our understanding of the vertical turbulence variability within the surface boundary layers in the ocean is mainly based on the turbulent flow studies over the solid wall-layer turbulence with the corresponding "law of the wall" (LOW), = u * 3 (kz) , where k = 0.41 is von Karman's constant [52,53] and u * is water-side friction velocity defined as Here, τ is wind stress that assumes to be constant across air-sea interface so that τ = u 2 * ρ = u 2 * a ρ a [54], where u * a and ρ a are water-side friction velocity and air density, respectively. Terray et al. [16] suggested enhanced values for relative to the LOW in the upper ocean, and we used their scaling method for wave-generated turbulence, Wave (see Appendix B). In addition to wave, the transfer of heat between water and air increases within the upper ocean boundary layer [51] due to that the surface buoyancy flux, B 0 , is used to investigate the surface convection role on turbulence [4,55].
where g is gravity, α is the thermal expansion coefficient of water, β is the coefficient of salinity expansion, s is the surface salinity, and E is the evaporation rate.

Results and Discussion
Two sets of experiments were conducted. The first set was to establish the power-law exponent of m and n. The data were collected at a depth of z = 0.25 m downstream of the grid-generated turbulence, 25< x/M <38, and at an equal distance from the horizontal sidewalls. In the second set of experiments, the data were collected within the surface boundary layer at a x/M = 37 to clarify how the magnitudes of THF, χ, and change when approaching the water surface. The data were collected between depths of 0.005 m < z < 0.17 m vertically and at an equal distance from sidewalls. The heater was set to the high setting (P 0 = 8.3 kW) when χ and were measured.
As mentioned in Section 2.1, χ and values are measured by temperature spectrum and empirical Nasmyth spectrum. Comparing the shear spectrum with the Nasmyth spectrum for < U >= 0.167 m/s at x/M = 37 indicates an acceptable fitness between the graphs (Figure 2). The shear spectrum starts to rise at a wavenumber by a factor of 70 below its peak (Figure 2), which shows the spectrum contains more than 90 % of all shear variance as mentioned in Rockland scientific international note 28 [29].
This section uses four parts to represent the results of our experiment. Section 3.1 displays the power-law decay results and the modified grid-generated turbulence model. The changes in THF, χ, and within the surface boundary layer are presented in Sections 3.2 and 3.3. Finally, the effects of non-breaking waves on the near-surface are shown in Section 3.4.

Steady and Spatially Decaying Background Turbulence-Grid-Generated Turbulence
The VMP200 [27] measured , χ, and temperature variance < θ 2 > at a constant depth of z = 0.25 m and selected distances along the tank centerline (25 < x/M < 38) for < U >= 0.183 m/s. The data collected for < U >= 0.167 were measured between 28 < x/M < 38. We estimated a background heat flux of about 25 W/m 2 from the equation of Q =| (ρC p K θ ∂<T> ∂z )/A w | [56], where A w is the water surface area. It likely changes during the experiments (the measured background heat flux is the summation of heat transfer from all sides of the tank walls and the water surface through the whole length of the tank).
The linear fitting line is used to find the power-law exponent for χ and θ 2 that they were m = 1.13 and m = 1.25, respectively (Figure 3). We found that the power-law exponent for the decay of values was n = 1.04 and 1.16 for the flow velocities of 0.167 m/s and 0.183 m/s, respectively (Figure 3). Based on the observed and χ in our experiment ( Figure 3) and following Zhou et al. [43] and Bogucki et al. [57], we found the following set of equations for χ and as a function of mean flow velocity <U> and distance to turbulence-generating x: The set of equations are defined due to observed m and n in our experiment. The values are considered to be as m + 1 = 1.13 + 1 = 2.13 for χ and the average of n + 1 = (1.04 + 1. 16

Observations of Free Convective Flow and Associated Vertical Heat Flux
Data were collected at x/M = 37, where the air-water temperature difference was approximately uniform during the experiment. The average net surface vertical heat flux of Q ≈| 25 | W/m 2 in our experiment was not substantial in comparison to the mean range of net vertical heat flux in the ocean Q <| 150 | W/m 2 [58], and also the heated-grid generated a horizontal heat flux of 16.6 kW/m 2 , which is smaller than of the 5-year mean of the south China Sea front region ∼4.6 * 10 2 kW/m 2 [59] (for a mixed boundary layer of L m = 35 m).
The surface heat flux increased the temperature gradient near the water surface, z = 0.028 m up to the water surface (Figure 4a). Due to the water surface cooling, the positive temperature gradient results in the near-surface convection in our experiment. The temperature drives the density as, in our freshwater experiment, the salinity is close to zero in the entire water depth. The water density decreased about ≈0.132 kg m 3 from the water surface to a depth of 0.14 m, causing an unstratified boundary layer.
An appropriate normalization method is not identified for χ values; therefore, χ values are depicted as non-normalized values in Figure 4b similar to the Bogucki et al. [12,13] and Peterson and Fer [36] works. Comparing χ (Figure 4a) and temperature (Figure 4b) reveals that the gradient of χ varies when the temperature gradient changes (it is shown with red lines). χ value increases by approaching the water surface (Figure 4b) as the temperature gradient increases, which shows that the large-scale vertical temperature gradient change, d < T > /dz, is consistent with the value of small-scale temperature fluctuation, χ.
The average of the THF values (the red line shows it in Figure 4c) are approximately constant within the surface layer; however, the THF's variance increases up to 50% at 0 < z < 0.028 m, which is in the range of the wave-effected layer (see Figure 5c).

Vertical Profile Observations
The data were collected at x/M = 37 from a depth of z = 0.17 m up to z = 0.005 m to investigate the vertical profile. Scaling of aids in understanding and describing the boundary layer physics ( Figure 5). is normalized by the grid-generated turbulent kinetic energy dissipation rate called "background turbulence" 0 ( Figure 5a). The surface buoyancy flux B 0 (Equation (5)) and the wave-generated turbulent kinetic energy dissipation rate wave (Equation (3)) are also used to normalize (Figure 5b,c). 0 was measured at the tank center (z = 0.25 m). The background turbulence ranges between 0 ≈ 0.5 − 3 × 10 −5 m 2 /s 3 that are within the subset of ocean range, ≈ 10 −11 − 10 −3 m 2 /s 3 [60][61][62]. The normalized with background turbulence for different velocities converge together in deeper water; however, they diverge by approaching the water surface (z > 0.02 m). In Figure 5b, values are normalized with the summation of the surface buoyancy flux B 0 and background turbulence 0 , i.e., /(B 0 + 0 ) (Figure 5b). The Obukhov length scale L 0 = ku * 3 /B 0 [47] is used to normalize the water depth (Figure 5b), which characterizes the relative importance of the shear and buoyant convection in the boundary layer [47,63]. The change in B 0 for the mean flow velocity range of 0.066 m/s to 0.183 m/s was small ≈2.8% (Table 1) given that the air temperature, humidity, and water temperature did not change during the experiment. The similar rates for both scaling methods (Figure 5a,b) suggest that the background turbulence is larger than the turbulence generated with the surface convection, and the convergence of / 0 in depth indicates the dominance of the grid-generated turbulence.
The non-breaking surface wave is another source of turbulence in addition to the surface convection. Given that the source of propagating waves is not clear, we must note that our wave-generated measurements may not be truly representative of results that would be obtained in the field. The Doppler shifting correction [64,65] was performed on the surface elevation spectrum and graphed in Figure 6a for the mean flow velocity of < U >= 0.125 m/s, which represented the presence of the surface waves. In addition to the surface elevation spectrum, the visual observations show that the uniform surface waves are created along the whole tank when the water moves along the tank with the grid (Figure 6b) and without the grid (the figure does not show) inside the tank. The wave number of elevation spectrum peak, K max (Figure 6a), is considered to calculate the wavelength λ w of non-breaking waves. The wavelength ranges between 0.038 m < λ w < 0.098 m (Table 1) categorized as the gravity or gravity-capillary regimes [22]. The integration of elevation spectrum S(K) gives the < η 2 >= ∞ 0 S(K)dK [42], hereupon the wave height H = 2 < η > . The normalized with the summation of the wave-generated turbulence Wave , Equation (A9), and background turbulence, 0 , i.e., /( Wave + 0 ) are presented in Figure 5c. The normalized values converge together from water surface to depth of z = 1.25λ w (Figure 5c), where the average of normalized values is approximately one. This indicates that the wave-generated turbulence plays a dominant role on the near-surface value in comparison to the surface convection and internal turbulence, and that is why the layer from the water surface to the depth of z = 1.25λ w is called the "wave-affected layer" in our paper.

Wave-Generated
The mean flow velocity increase resulted in the increment of the mean square slope (MSS) and wavelength λ w (Table 1). has a larger value while the MSS increases (Figure 6c), and they have the same order of magnitude, i.e., MSS ∼ below the water surface z = 0.005 m. for the different mean flow velocities were compared at three different depths (Figure 6d) to investigate the effects of the non-breaking surface wave and the grid-generated turbulence on along the water column. Figure 6d indicates that at a depth of z = 0.145 m, change rate (red line) is higher than at z = 0.005 m (black line) while the surface heat flux was constant.
The turbulent kinetic energy dissipation rate comparison between the mean flow velocity range of 0.066 m/s to 0.183 m/s at the different depths shows that values increase ≈780%, 370%, and 95% at a depth of z = 0.145 m, z = 0.0175 m, and z = 0.005 m, respectively (Figure 6d). By considering constant heat flux during the experiment, the change rate should be the same when approaching the water surface if the grid was the only source of turbulence. The different change rates, like wave scaling results (Figure 5c), indicate that the non-breaking waves play a dominant role in on the surface layer. The subtraction of the averaged-within the wave-affected layer, 1.25λ w , from the background turbulent kinetic energy dissipation rate, 0 , is considered the wave-generated turbulent kinetic energy dissipation rate, wave−generated = 1.25λ w − 0 . The wave−generated for wavelength of 0.038 m <λ w < 0.098 m ranged to 4.9 × 10 −6 − 7 × 10 −6 m 2 /s 3 (Figure 6e). Given that the breaking wave is often related to the enhanced turbulent kinetic energy dissipation rate within the surface layer, wave−generated values represent the significant role of non-breaking waves on the upper-ocean mixing intensity. It indicates that welldocumented researches are necessary to shed light on the poor understanding of the non-breaking wave-generated turbulence. Furthermore, it depicts that the precise values of non-breaking wave-generated turbulent kinetic energy dissipation rate are required to quantify the air-sea interface processes.

Conclusions
Oceanic turbulence measurements are practically impossible when attempting to address processes within a few upper centimeters below the wave ocean surface. Our experiment aimed to investigate how the weak surface forcings and the horizontal heat and eddy fluxes affect the near-surface layer in a controlled laboratory setting like the front in the Gulf of Mexico where the Mississippi River with massive horizontal eddy fluxes reaches the ocean and generates the front. In our lab experiment, the internal and horizontal heat flux generated by the grids were the subsets of the ocean [60][61][62] and heat fluxes range [59].
While analyzing the vertical variability, we have observed that there is a "waveaffected layer" from the water surface to a depth of z ≈ 1.25λ w . Turbulence kinetic energy dissipation rate associated with non-breaking waves wave−generated ranged to 4.9 × 10 −6 -7 × 10 −6 m 2 /s 3 for the wavelength range of 0.038 m <λ w < 0.098 m categorized as the gravity and gravity-capillary regimes [21][22][23]. The increase in the MSS resulted in the larger , and they have the same order of magnitude MSS ∼ . Given that the non-breaking waves typically cover a larger fraction of the ocean surface, 90-100% [17], than breaking waves, the results indicate their significant contribution to the ocean energy budget. Therefore, the non-breaking wave's turbulence kinetic energy dissipation rate budget has to be considered to properly quantify the air-sea interface processes such as cool skin thickness [66], which is a fundamental parameter required for quantifying the physical process taking place at the air=-sea interface like the gas transfer [4] and heat transfer [67].
We also found that χ changes, which is the small-scale temperature fluctuation, are consistent with the large-scale temperature gradient, d < T > /dz, changes. The value of the THF is approximately constant within the surface layer. It represents that the measured THF near the water surface can be considered a surface water THF in the ocean, challenging to measure directly.
In addition, we observed that the power-law exponent of the tank is m = 1.13 and m = 1.25 for the decay of χ and temperature variance < θ 2 >, respectively, and the decay of equals n = 1.04 and n = 1.16 for the velocities 0.167 m/s and 0.183 m/s. The future work would be to connect the laboratory observation to the field observation. We will explore whether the THF value is constant within the upper ocean boundary layer as observed in this research. We would also like to study the importance of below the non-breaking gravity wave and investigate if the observed affected layer in this paper is common properties below the non-breaking and breaking waves in the ocean. Acknowledgments: We would like to thank RSMAS Ph.D. students Sanchit Mehta, Hanjing Dai, and Andrew W. Smith for help with the preparation of the experiment.

Conflicts of Interest:
The authors declare no conflicts of interest.

Appendix A. Transport Equation of and χ for Grid-Generated Turbulence
The internal turbulence was generated by passing water through the grids, see Section 4, and the internal turbulence rate was controlled by changing the mean flow velocity of passing water. Batchelor and Townsend [68], and Warhaft and Lumley [46] observed that the mean fluctuating turbulent kinetic energy, < q 2 >=< u 2 > + < v 2 > + < w 2 > and temperature variance, < θ 2 >, decay with distance from the grid. Tresso and Munoz [69] reported the existence of steady-state turbulence in the point behind the grid, where U is the flow velocity in the streamwise direction, x, and is represented by U =< U > +u (see Figure 2), where the mean flow velocity is represented by < U > and the fluctuation part of velocity is u . Similarly, the other velocity components are expressed by V =< V > +v in the y-direction, and W =< W > +w in the zdirection. The water temperature is defined by T =< T > +θ , where < T > is the mean temperature value and θ is the fluctuation part of the temperature. The governing equations for the evolution of mean turbulent kinetic energy and the temperature variance of the homogeneous and isotropic shear flow are [7] d < q 2 > 2dt ≈ P − (A1) and Here, P and P θ are the turbulent kinetic energy production and the temperature variance production rate, respectively. In the grid-generated turbulence, the production terms are equal to P = 0. Assuming a constant mean current, we get [43] − < U > 2 Numerous experiments [43,68] suggest that the downstream decay of the mean turbulent kinetic energy and temperature variance behind the grids are where n and m are the exponents of the mean turbulent kinetic energy and temperature variances, respectively. The constants A and B depend on the grid geometry and are typically determined empirically, and x is the horizontal distance from the turbulence generating grid. By substituting Equations (5) and (6) into (3) and (4), the power-law decay exponents for and χ are [43] < >≈