Turbulent motions in the thin ocean surface boundary layer control exchanges of momentum, heat and trace gases between the atmosphere and ocean. However, present parametric equations of turbulent motions that are applied to global climate models result in systematic or substantial errors in the ocean surface boundary layer. Significant mixing caused by surface wave processes is missed in most parametric equations. A Large Eddy Simulation model is applied to investigate the wave-induced mixed layer structure. In the wave-averaged equations, wave effects are calculated as Stokes forces and breaking waves. To examine the effects of wave parameters on mixing, a series of wave conditions with varying wavelengths and heights are used to drive the model, resulting in a variety of Langmuir turbulence and wave breaking outcomes. These experiments suggest that wave-induced mixing is more sensitive to wave heights than to the wavelength. A series of numerical experiments with different wind intensities-induced Stokes drifts are also conducted to investigate wave-induced mixing. As the wind speed increases, the influence depth of Langmuir circulation deepens. Additionally, it is observed that breaking waves could destroy Langmuir cells mainly at the sea surface, rather than at deeper layers.Atmosphere 2020, 11, 207 2 of 13 is also affected by Langmuir turbulence and is primarily generated by surface waves. The depth at which LC the upper ocean is largely determined by the Stokes drift e-folding depth and the mixed layer thickness [10,11]. The injection of turbulent kinetic energy from breaking surface gravity waves provides another significant source of upper-ocean turbulence [12]. Increased turbulence due to breaking waves (BW) decays following reductions in significant wave heights [13,14].The importance of wave-induced mixing is demonstrated in upper-ocean mixed layer studies. However, most turbulence schemes such as the Mellor-Yamada [15] or K-profile [16] parameterizations do not account for wave processes. These turbulence models are handicapped by the failure to include mixed layer physical processes such as BW and LC. This results in ocean mixed layer models producing too shallow mixed layer depths in high-latitude oceans (especially the Southern Ocean) [17][18][19][20][21][22] and highly diffused thermoclines in tropical oceans [8,[23][24][25][26], leading to weaker El Niño in climate models [27]. Furthermore, the principal way in which waves influence the oceanic circulation is through wave-induced vertical mixing, which is currently limited by model resolution, thus necessitating parameterization [28][29][30]. To more realistically reproduce the upper ocean in ocean general circulation models, a number of ocean mixed layer models have been proposed [28][29][30][31][32][33][34][35][36][37][38][39]. For example, a parameterization scheme including wave-induced mixing is proposed to modify current Mellor-Yamada and K-profile parameterization turbulence models [36,40]. The comparison of numerical results with observat...