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

Abstract This 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.