The KPP boundary layer scheme for the ocean: revisiting its formulation and benchmarking one-dimensional simulations relative to LES
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- NOAA Geophysical Fluid Dynamics Lab., Princeton, NJ (United States); Princeton Univ., NJ (United States). Program in Atmospheric and Oceanic Sciences
- National Center for Atmospheric Research, Boulder, CO (United States). Climate and Global Dynamics Lab.
- Leibniz Inst. for Baltic Sea Research, Rostock (Germany)
We evaluate the Community ocean Vertical Mixing (CVMix) project version of the K-profile parameterization (KPP) for modeling upper ocean turbulent mixing. For this purpose, one-dimensional KPP simulations are compared across a suite of oceanographically relevant regimes against horizontally averaged large eddy simulations (LES). We find the standard configuration of KPP consistent with LES across many forcing regimes, supporting its physical basis. Our evaluation also motivates recommendations for KPP “best practices” within ocean circulation models, and identifies areas where further research is warranted. The original treatment of KPP recommends the matching of interior diffusivities and their gradients to the KPP predicted values computed in the ocean surface boundary layer (OSBL). However, we find that difficulties in representing derivatives of rapidly changing diffusivities near the base of the OSBL can lead to loss of simulation fidelity. To mitigate this difficulty, we propose and evaluate two computationally simpler approaches: (1) match to the internal predicted diffusivity alone, (2) set the KPP diffusivity to zero at the OSBL base. We find the KPP entrainment buoyancy flux to be sensitive to vertical grid resolution and details of how to diagnose the KPP boundary layer depth. We modify the KPP turbulent shear velocity parameterization to reduce resolution dependence. Additionally, an examination of LES vertical turbulent scalar flux budgets shows that the KPP parameterized non-local tracer flux is incomplete due to the assumption that it solely redistributes the surface tracer flux. This result motivates further studies of the non-local flux parameterization.
- Research Organization:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Los Alamos National Laboratory (LANL), Los Alamos, NM (United States); National Center for Atmospheric Research (NCAR), Boulder, CO (United States); Princeton Univ., NJ (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Biological and Environmental Research (BER); USDOE National Nuclear Security Administration (NNSA); National Science Foundation (NSF); Princeton Univ. (United States)
- Grant/Contract Number:
- AC52-06NA25396; AC05-00OR22725; SC0012605
- OSTI ID:
- 1480901
- Alternate ID(s):
- OSTI ID: 1475346; OSTI ID: 1480902
- Report Number(s):
- LA-UR-16-26992
- Journal Information:
- Journal of Advances in Modeling Earth Systems, Vol. 10, Issue 11; ISSN 1942-2466
- Publisher:
- American Geophysical Union (AGU)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
The GFDL Global Ocean and Sea Ice Model OM4.0: Model Description and Simulation Features
|
journal | October 2019 |
Influence of Mesoscale Features on Mixed Layer Dynamics in the Arabian Sea
|
journal | May 2019 |
Global Impacts of Subseasonal (<60 Day) Wind Variability on Ocean Surface Stress, Buoyancy Flux, and Mixed Layer Depth
|
journal | December 2019 |
Comparing Ocean Surface Boundary Vertical Mixing Schemes Including Langmuir Turbulence
|
journal | November 2019 |
Similar Records
Southern Ocean Uptake in the MPAS-Ocean Model
Transforming the representation of the boundary layer and low clouds for high-resolution regional climate modeling: Final report