Ocean Boundary Layer Research Articles

Submesoscales are a significant turbulence source in global ocean surface boundary layer

The turbulent ocean surface boundary layer is a key part of the climate system affecting both the energy and carbon cycles. Accurately simulating the boundary layer is critical in improving climate model performance, which deeply relies on our understanding of the turbulence in the boundary layer. Turbulent energy sources in the boundary layer are traditionally believed to be dominated by waves, winds and convection. Recently, submesoscale phenomena with spatial scales of 0.1~10 km at ocean fronts have been shown to also make a contribution. Here, by applying a non-dimensional turbulent kinetic energy budget equation, we show that the submesoscale geostrophic shear production at fronts is a significant turbulent energy source within the ocean boundary layer away from the sea surface. The contribution reaches 34% of the total dissipation in winter and 17% in summer at the mid-depth of the boundary layer, despite its intermittency in space and time. This work indicates fundamental deficiencies in previous conceptions of ocean boundary layer turbulence, and invites a reappraisal of the sampling scale in observations, model resolution and parameterizations, and other consequences of the global energy budget.

Impact of ocean vertical-mixing parameterization on Arctic sea ice and upper-ocean properties using the NEMO-SI3 model

Abstract. We evaluate the vertical turbulent-kinetic-energy (TKE) mixing scheme of the NEMO-SI3 ocean–sea-ice model in sea-ice-covered regions of the Arctic Ocean. Specifically, we assess the parameters involved in TKE mixed-layer-penetration (MLP) parameterization. This ad hoc parameterization aims to capture processes that impact the ocean surface boundary layer, such as near-inertial oscillations, ocean swells, and waves, which are often not well represented in the default TKE scheme. We evaluate this parameterization for the first time in three regions of the Arctic Ocean: the Makarov, Eurasian, and Canada basins. We demonstrate the strong effect of the scaling parameter that accounts for the presence of sea ice. Our results confirm that TKE MLP must be scaled down below sea ice to avoid unrealistically deep mixed layers. The other parameters evaluated are the percentage of energy penetrating below the mixed layer and the length scale of its decay with depth. All these parameters affect mixed-layer depth and its seasonal cycle, surface temperature, and salinity, as well as underlying stratification. Shallow mixed layers are associated with stronger stratification and fresh surface anomalies, and deeper mixed layers correspond to weaker stratification and salty surface anomalies. Notably, we observe significant impacts on sea-ice thickness across the Arctic Ocean in two scenarios: when the scaling parameter due to sea ice is absent and when the TKE mixed-layer-penetration process vanishes. In the former case, we observe an increase of several meters in mixed-layer depth, along with a reduction in sea-ice thickness ranging from 30 to 40 cm, reflecting the impact of stronger mixing. Conversely, in the latter case, we notice that a shallower mixed layer is accompanied by a moderate increase in sea-ice thickness, ranging from 10 to 20 cm, as expected from weaker mixing. Additionally, interannual variability suggests that experiments incorporating a scaling parameter based on sea-ice concentration display an increased mixed-layer depth during periods of reduced sea ice, which is consistent with observed trends. These findings underscore the influence of enhanced ocean mixing, through specific parameterizations, on the physical properties of the upper ocean and sea ice.

Rapid Downwelling of Tracer Particles Across the Boundary Layer and Into the Pycnocline at Submesoscale Ocean Fronts

AbstractA neutrally buoyant float deployed in an atmospherically driven turbulent ocean boundary layer on the dense side of a submesoscale front was repeatedly carried across the boundary layer by the turbulence and then trapped beneath the slumping front. Lagrangian particles in a large‐eddy simulation of a similar baroclinically unstable front forced by surface cooling move along convergent surface filaments toward filament junctions. They are also caught by convective plumes that downwell them at speeds similar to those of the float. Subsequently, some are trapped in the pycnocline by frontal slumping due to ageostrophic secondary frontal circulations. In both observations and simulations, boundary layer turbulence and frontal circulations work together to trap and subduct particles from the mixed layer. The small‐scale boundary layer motions move them vertically within the boundary layer and larger, submesoscale frontal circulations move them laterally out of the boundary layer and under the slumping fronts.

The K‐Profile Parameterization Augmented by Deep Neural Networks (KPP_DNN) in the General Ocean Turbulence Model (GOTM)

AbstractThis study utilizes Deep Neural Networks (DNN) to improve the K‐Profile Parameterization (KPP) for the vertical mixing effects in the ocean's surface boundary layer turbulence. The deep neural networks were trained using 11‐year turbulence‐resolving solutions, obtained by running a large eddy simulation model for Ocean Station Papa, to predict the turbulence velocity scale coefficient and unresolved shear coefficient in the KPP. The DNN‐augmented KPP schemes (KPP_DNN) have been implemented in the General Ocean Turbulence Model (GOTM). The KPP_DNN is stable for long‐term integration and more efficient than existing variants of KPP schemes with wave effects. Three different KPP_DNN schemes, each differing in their input and output variables, have been developed and trained. The performance of models utilizing the KPP_DNN schemes is compared to those employing traditional deterministic first‐order and second‐moment closure turbulent mixing parameterizations. Solution comparisons indicate that the simulated mixed layer becomes cooler and deeper when wave effects are included in parameterizations, aligning closer with observations. In the KPP framework, the velocity scale of unresolved shear, which is used to calculate ocean surface boundary layer depth, has a greater impact on the simulated mixed layer than the magnitude of diffusivity does. In the KPP_DNN, unresolved shear depends not only on wave forcing, but also on the mixed layer depth and buoyancy forcing.

A preliminary observational study on the characteristics of surface turbulent fluxes over the South China Sea Islands

AbstractIn recent years, there has been a rise in human activities in oceanic areas, making the land–atmosphere interactions over islands a major scientific concern on a global scale. Examining the observation data from offshore areas enables a more comprehensive understanding of the turbulent fluxes in offshore atmospheric environments, patterns of momentum, energy and material exchange between the atmosphere and underlying surface in an oceanic boundary layer, and development of a heterogeneous atmospheric boundary layer. The related findings will assist in developing theoretical models and parameterization schemes to simulate the influence of heterogeneous surfaces on land–atmosphere interactions on the South China Sea Islands. Existing studies on the turbulent fluxes over the South China Sea Islands were mainly conducted on the Nansha Islands, whereas studies on the waters of the South China Sea are scarce. In this study, we used 10 Hz high‐frequency turbulence measurements to calculate the latent and sensible heat fluxes over the South China Sea Islands using the eddy correlation method. These findings were then compared with data from the Dunhuang Gobi, Ordos desert, and Xilingol grassland regions in inland China, along with the observed net radiation and surface heat fluxes. The findings indicate that the energy fluxes over the South China Sea in summer exhibit prominent diurnal variations. The magnitude of either latent or soil heat flux is low, and the net radiation is predominantly transformed into sensible heat flux, which warms the atmosphere. Furthermore, the daily variation curves of sensible and latent heat fluxes are influenced by intermittent turbulence on the islands and reefs, resulting in a less smooth pattern compared with soil heat flux. Although the South China Sea Islands have small land areas and are surrounded by the sea, the land–atmosphere interactions over the underlying surface of this region are similar to those over the underlying surface of grasslands in inland China during summer. The daily mean sensible heat flux on the islands is higher than that in an inland area, and the time lag in its response to sunrise is longer than that in inland areas by approximately 1 h. The overall energy balance ratio is approximately 0.75, c which is in line with the average level, but an energy balance residual of approximately 25% still exists. Furthermore, extreme weather conditions, such as typhoons, can disrupt the diurnal variations of sensible and latent heat fluxes, and the cyclical patterns are subsequently restored.

The Wave-Coherent Stress and Turbulent Structure over Swell Waves

Abstract The generation of ocean surface waves by wind has been studied for a century, giving rise to wave forecasting and other crucial applications. However, the reacting force of swell waves on the turbulence in the marine atmospheric boundary layer (ABL) remains unknown partly due to the unclear magnitude and profile of wave-coherent (WC) stress. In this study, the intersection frequency between the energy-containing range and inertial subrange range in the turbulent spectra is identified based on the attached eddy model (AEM), as the intersection modulated by swell wave could help to comprehend the physical process between the ocean surface wave and the marine ABL. Using observations from a fixed platform located in the South China Sea, this study shows that the intersection when the WC stress accounts for a lower proportion of the total wind stress (<10%) follows U/(2πz) given by AEM, where U is the wind speed and z is the height. While the intersection depends on the drag coefficient of WC stress for the case, WC stress accounts for a large part of the total wind stress (>10%). Considering the unclear magnitude and profile of WC stress, this study derives a new function to depict the WC stress. Significance Statement Turbulence located in the energy-containing range of the spectra is the primary source of momentum, heat, and water vapor fluxes between the ocean and marine atmospheric boundary layer (ABL). Thus, a better understanding of the turbulent structure within the wave boundary layer can help us accurately parameterize those fluxes. In addition to absorbing energy from the marine ABL, waves, especially swell waves, influence the turbulent structure. However, until now, the turbulent structure is unclear due to the unclear wave-coherent stress size and profile. In this study, a function that can depict the wave-coherent stress is proposed. We found that the wave-coherent stress can modulate the intersection between the energy-containing range and the inertial subrange range in the turbulent spectra.

Turbulent Ice–Ocean Boundary Layers in the Well-Mixed Regime: Insights from Direct Numerical Simulations

Abstract The meltwater mixing line (MML) model provides a theoretical prediction of near-ice water mass properties that is useful to compare with observations. If oceanographic measurements reported in a temperature–salinity diagram overlap with the MML prediction, then it is usually concluded that the local dynamics are dominated by the turbulent mixing of an ambient water mass with near-ice meltwater. While the MML model is consistent with numerous observations, it is built on an assumption that is difficult to test with field measurements, especially near the ice boundary, namely, that the effective (turbulent and molecular) salt and temperature diffusivities are equal. In this paper, this assumption is tested via direct numerical simulations of a canonical model for externally forced ice–ocean boundary layers in a uniform ambient. We focus on the well-mixed regime by considering an ambient temperature close to freezing and run the simulations until a statistical steady state is reached. The results validate the assumption of equal effective diffusivities across most of the boundary layer. Importantly, the validity of the MML model implies a linear correlation between the mean salinity and temperature profiles normal to the interface that can be leveraged to construct a reduced ice–ocean boundary layer model based on a single scalar variable called thermal driving. We demonstrate that the bulk dynamics predicted by the reduced thermal driving model are in good agreement with the bulk dynamics predicted by the full temperature–salinity model. Then, we show how the results from the thermal driving model can be used to estimate the interfacial heat and salt fluxes and the melt rate. Significance Statement We investigate the turbulent dynamics and thermodynamical properties of water masses below ice shelves using new data from high-resolution simulations. This is important because there are currently too few observations of ice–ocean boundary layer dynamics to construct reliable models of ice-shelf melting as a function of ocean conditions. Our results demonstrate that the turbulent diffusivities of salt and temperature are approximately equal in the well-mixed regime. This implies a linear correlation between the mean temperature and salinity profiles, consistent with the meltwater mixing line prediction that is often used to interpret polar observations. We take advantage of this correlation to propose a reduced model of ice–ocean boundary layers that can predict ice-shelf melt rates at relatively low computational cost.

Global variability in atmospheric new particle formation mechanisms

A key challenge in aerosol pollution studies and climate change assessment is to understand how atmospheric aerosol particles are initially formed1,2. Although new particle formation (NPF) mechanisms have been described at specific sites3–6, in most regions, such mechanisms remain uncertain to a large extent because of the limited ability of atmospheric models to simulate critical NPF processes1,7. Here we synthesize molecular-level experiments to develop comprehensive representations of 11 NPF mechanisms and the complex chemical transformation of precursor gases in a fully coupled global climate model. Combined simulations and observations show that the dominant NPF mechanisms are distinct worldwide and vary with region and altitude. Previously neglected or underrepresented mechanisms involving organics, amines, iodine oxoacids and HNO3 probably dominate NPF in most regions with high concentrations of aerosols or large aerosol radiative forcing; such regions include oceanic and human-polluted continental boundary layers, as well as the upper troposphere over rainforests and Asian monsoon regions. These underrepresented mechanisms also play notable roles in other areas, such as the upper troposphere of the Pacific and Atlantic oceans. Accordingly, NPF accounts for different fractions (10–80%) of the nuclei on which cloud forms at 0.5% supersaturation over various regions in the lower troposphere. The comprehensive simulation of global NPF mechanisms can help improve estimation and source attribution of the climate effects of aerosols.

Using NIW Observations to Assess Mixed Layer Parameterizations: A Case Study in the Tropical Atlantic

AbstractTropical sea surface temperature (SST) biases can cause atmospheric biases on global scales, hence SST needs to be represented well in climate models. A major source of uncertainties is the representation of turbulent mixing in the oceanic boundary layer, or mixed layer (ML). In the present study we focus on near‐inertial wave (NIW) induced mixing. The performance of two mixing schemes, Turbulent Kinetic Energy and K‐profile parameterization (KPP), is assessed at two sites (11.5°N, 23°W and 15°N, 38°W) in the tropical Atlantic. At 11.5°N, turbulence observations (eddy diffusivities, shear and stratification) are available for comparison. We find that the schemes differ in their representation of NIWs, but both under‐represent the observed enhanced diffusivities below the observed ML. However, we find that the models do mix below the ML at 15°N when a storm passes nearby. The near‐inertial oscillations remain below the ML for the following 10 days. Near‐inertial kinetic energy (NIKE) biases in the models are not directly correlated with the wind speed, the MLD biases, or the stratification at the ML base. Instead, NIKE biases are sensitive to the vertical mixing scheme parameterization. NIKE biases are lowest when the KPP scheme is used.

Sensitivity of boundary layer features to depth-dependent baroclinic pressure gradient and turbulent mixing in an ocean of finite depth

The present numerical study builds on Ekman (1905)’s work in surface boundary layer and extends the boundary value problem to overcome some of its limitations. Previous studies addressed model’s limitations by assuming that deviations from observations are usually ascribed to different eddy viscosity shapes, but seldom to the presence of baroclinic pressure gradients and shallow seas, which are the mainstays of this work. Improved solutions in the ocean boundary layer are obtained considering both depth-dependent wind-induced eddy viscosity and horizontal density gradients, ranging from well-mixed to highly-stratified conditions in a finite-depth ocean. High-order numerical solutions extend those in previous analytical and numerical works in the literature and widens the parameter space analyzed. Remarkably, the current profiles are obtained without ambiguity as a truly superposition of a geostrophic and a ageostrophic terms. Results indicate that, for a vertically-uniform eddy viscosity without density gradients and in shallow waters, currents are practically aligned with wind. As depth increases, misalignment between currents and wind increases and the complexity of the vertical structure increases. At large depths, Ekman’s values are attained, i.e., deflection angles relative to wind direction, θW, are θS−θW=−45∘ at the surface, where the current is maximum, and θT−θW=−90° for the depth-integrated transport (negative for deflections to the right in the Northern Hemisphere). These features remain regardless of the magnitude of the eddy-viscosity. For non-uniform eddy viscosity, θS−θW decreases from −45∘ up to −90∘ from low to high stratification level, respectively, whereas θT−θW is rather insensitive (θT−θW≈−90∘). Contrary to wind effects, the presence of the only vertically-uniform density gradient forcing, with constant eddy viscosity, deflects the surface angle θS−θD=+45∘ relative to the density-gradient direction, θD, in deep waters. Maximum currents no longer occur at the surface in this case. For non-uniform density-gradient profiles, current magnitude decreases overall while θS−θD≈+45∘ as long as the gradient affects the entire surface boundary layer. The deflection angle θT−θD∼+95∘ remains less sensitive to changes in density-gradient profile.

How Currents Trigger Extreme Sea Waves. The Roles of Stokes Drift, Eulerian Return Flow, and a Background Flow in the Open Ocean

AbstractA deterministic system of ocean surface waves and flow in the oceanic boundary layer is key to understanding the dynamics of the upper ocean. For the description of such complex systems, a higher‐order shear‐current modified nonlinear Schrödinger equation is newly derived and then used to physically interpret the interplay between Stokes drift, Eulerian return flow due to a passing wave group, and an open‐ocean vertically sheared flow in the extreme sea wave generation. The conditions for the suppression or enhancement of the modulation instability in the rogue wave dynamics in the presence of a background flow are reported, whose relevance and influence to the Craik‐Leibovich type 2 instability in triggering a Langmuir‐type circulation is discussed. The findings highlight the need for future studies to establish and assess the energy transfer from waves to currents or in the reversing order, asserting a plausible physical mechanism for the dissipation of the surface wave energy through wave‐current interactions in the open ocean.

Submesoscale Variability and Basal Melting in Ice Shelf Cavities of the Amundsen Sea

AbstractMelting of ice shelves can energize a wide range of ocean currents, from three‐dimensional turbulence to relatively large‐scale boundary currents. Here, we conduct high‐resolution simulations of the western Amundsen Sea to show that submesoscale eddies are prevalent inside ice shelf cavities. The simulations indicate energetic submesoscale eddies at the top and bottom ocean boundary layers, regions with sharp topographic slopes and strong lateral buoyancy gradients. These eddies play a substantial role in the vertical and lateral (along‐isopycnal) heat advection toward the ice shelf base, enhancing the basal melting in all simulated cavities. In turn, the meltwater provides strong buoyancy gradients that energize the submesoscale variability, forming a positive loop that could affect the overall efficiency of heat exchange between the ocean and the ice shelf cavity. Our study implies that submesoscale‐induced enhancement of basal melting may be a ubiquitous process that needs to be parameterized in coarse‐resolution climate models.

Observations of the Upper Ocean from Autonomous Platforms During the Passage of Extratropical Cyclone Epsilon (2020)

Hurricane Epsilon (2020) was a late-season, category-3 tropical cyclone that underwent extratropical transition and became Extratropical Cyclone Epsilon on October 26. The upper ocean response to the passage of the storm was observed by three types of autonomous platforms: an eXpendable Spar buoy, an Air-Launched Autonomous Micro-Observer profiling float, and two Seagliders. Taken together, this array enabled the rare collection of contemporaneous observations of the upper ocean, air-sea interface, and atmospheric boundary layer before, during, and after the passage of the storm. The evidence presented suggests that Extratropical Cyclone Epsilon contributed to breaking down the residual North Atlantic summer stratification regime and accelerated the shift to the prolonged ocean cooling associated with winter. The synergistic capabilities of the observational array are significant for two reasons: (1) by enabling the comparison of complementary atmosphere and ocean observations, taken from different platforms, they permit a comprehensive approach to better understand how storm-induced momentum, heat, and moisture fluxes alter upper ocean structure, and (2) they demonstrate the ability of future, targeted deployments of similar observational arrays to assess the fidelity of coupled ocean-atmosphere-wave numerical prediction models.

A New Hybrid Mass‐Flux/High‐Order Turbulence Closure for Ocean Vertical Mixing

AbstractWhile various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first‐order closure schemes such as the K‐Profile Parameterization (KPP) lack energetic constraints; two‐equation models like directly solve an equation for the turbulent kinetic energy but do not account for non‐diffusive fluxes, and high‐order closures that include the high‐order transport terms are computationally expensive. To address these, we extend the Assumed‐Distribution Higher‐Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second‐order and higher‐order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non‐diffusive fluxes without the computational cost of a full higher‐order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and for surface buoyancy‐driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

A novel (2+1)-dimensional nonlinear evolution equation for weakly stratified free-surface boundary layers

To get an insight into the dynamics of the oceanic surface boundary layer we develop an asymptotic model of the nonlinear dynamics of linearly decaying three-dimensional long-wave perturbations in weakly stratified boundary-layer flows. Although in nature the free-surface boundary layers in the ocean are often weakly stratified due to solar radiation and air entrainment caused by wave breaking, weak stratification has been invariably ignored. Here, we consider an idealized hydrodynamic model, where finite-amplitude three-dimensional perturbations propagate in a horizontally uniform unidirectional weakly stratified shear flow confined to a boundary layer adjacent to the water surface. Perturbations satisfy the no-stress boundary condition at the surface. They are assumed to be long compared with the boundary-layer thickness. Such perturbations have not been studied even in a linear setting. By exploiting the assumed smallness of nonlinearity, wavenumber, viscosity and the Richardson number, on applying triple-deck asymptotic scheme and multiple-scale expansion, we derive in the distinguished limit a novel essentially two-dimensional nonlinear evolution equation, which is the main result of the work. The equation represents a generalization of the two-dimensional Benjamin–Ono equation modified by the explicit account of viscous effects and new dispersion due to weak stratification. It describes perturbation dependence on horizontal coordinates and time, while its vertical structure, to leading order, is given by an explicit analytical solution of the linear boundary value problem. It shows the principal importance of weak stratification for three-dimensional perturbations.

Parameterizing Vertical Mixing Coefficients in the Ocean Surface Boundary Layer Using Neural Networks

AbstractVertical mixing parameterizations in ocean models are formulated on the basis of the physical principles that govern turbulent mixing. However, many parameterizations include ad hoc components that are not well constrained by theory or data. One such component is the eddy diffusivity model, where vertical turbulent fluxes of a quantity are parameterized from a variable eddy diffusion coefficient and the mean vertical gradient of the quantity. In this work, we improve a parameterization of vertical mixing in the ocean surface boundary layer by enhancing its eddy diffusivity model using data‐driven methods, specifically neural networks. The neural networks are designed to take extrinsic and intrinsic forcing parameters as input to predict the eddy diffusivity profile and are trained using output data from a second moment closure turbulent mixing scheme. The modified vertical mixing scheme predicts the eddy diffusivity profile through online inference of neural networks and maintains the conservation principles of the standard ocean model equations, which is particularly important for its targeted use in climate simulations. We describe the development and stable implementation of neural networks in an ocean general circulation model and demonstrate that the enhanced scheme outperforms its predecessor by reducing biases in the mixed‐layer depth and upper ocean stratification. Our results demonstrate the potential for data‐driven physics‐aware parameterizations to improve global climate models.

A Finite-Time Ensemble Method for Mixed Layer Model Comparison

Abstract This work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

Cool skin effect and warm skin phenomenon observed by shipboard radiometer in the Northwest Pacific

Sea surface temperature (SST) is an important variable in the study of ocean boundary layers and heat exchange. The accurate simulation and measurement of skin effects are vital to air–sea model processing and satellite SST retrieval. Shipboard measurements from eleven cruises in the Northwest Pacific between August 2015 and October 2018 were used to estimate the cool skin effect and compare model results. The temperature difference ΔT between the sea surface skin temperature (SSTskin), as measured by an infrared radiometer, and the sea surface depth temperature (SSTdepth) at around 4 meters showed a mean difference and a standard deviation of the same 0.2 K, with a total of 5-min 39909 measurements. Both daytime and nighttime ΔT values were compared to physical model simulations and were found to have relatively larger mean values. A set of new coefficients for an exponential parameterization of the cool skin effect was derived in the research area, which performed well in comparison to previous empirical models. In nighttime observations from two summer cruises, the reverse process of heat flux transfer from the air to the sea in the form of a warm skin was distinguished. There were 667 positive ΔT values out of the 1917 nighttime observations, with magnitudes ranging from around 0 to 0.3 K. A high proportion of the cases of the warm skin phenomenon occurred when the air was very humid and much warmer than the sea surface.

Observation and parameterization of turbulent Reynolds stress in the ocean surface boundary layer under swell-dominated conditions

The turbulent Reynolds stress in the ocean surface boundary layer (OSBL) is investigated using two sets of turbulence measurements conducted on an offshore tower-based platform in the coastal region of the northern South China Sea. Observations reveal that the turbulent Reynolds stress in the OSBL (τw) is mainly determined by surface waves rather than by wind and ocean currents under swell-dominated conditions. This stress is closely related to the Stokes drift velocity (us); a parameterization scheme is proposed for τw as a function of us. This scheme indicates that τw increases with wave amplitude, and decays exponentially with water depth. Further analysis shows that the turbulent stress estimated using this scheme has a high correlation with the observed stress, evidenced by a correlation coefficient of 0.92 during our two measurement periods.

Langmuir Turbulence Controls on Observed Diurnal Warm Layer Depths

AbstractThe turbulent ocean surface boundary layer (OSBL) shoals during daytime solar surface heating, developing a diurnal warm layer (DWL). The DWL significantly influences OSBL dynamics by trapping momentum and heat in a shallow near‐surface layer. Therefore, DWL depth is critical for understanding OSBL transport and ocean‐atmosphere coupling. A great challenge for determining DWL depth is considering wave‐driven Langmuir turbulence (LT), which increases vertical transport. This study investigates observations with moderate wind speeds (4–7 m/s at 10 m height) and swell waves for which breaking wave effects are less pronounced. By employing turbulence‐resolving large eddy simulation experiments that cover observed wind, wave, and heating conditions based on the wave‐averaged Craik‐Lebovich equation, we develop a DWL depth scaling unifying previous approaches. This scaling closely agrees with observed DWL depths from a year‐long mooring deployment in the subtropical North Atlantic, demonstrating the critical role of LT in determining DWL depth and OSBL dynamics.