On the Mellor–Yamada Turbulence Closure Scheme: The Surface Boundary Condition for q2

Michael W. Stacey Department of Physics, Royal Military College of Canada, Kingston, Ontario, Canada

Search for other papers by Michael W. Stacey in
Current site
Google Scholar
PubMed
Close
and
Stephen Pond Department of Oceanography, The University of British Columbia, Vancouver, British Columbia, Canada

Search for other papers by Stephen Pond in
Current site
Google Scholar
PubMed
Close
Restricted access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

A numerical model that uses a level-2½ turbulence closure scheme is used to compare two boundary conditions for the turbulent energy at the air–sea interface. One boundary condition, the most commonly used, sets the turbulent kinetic energy proportional to the friction velocity squared, while the other sets the vertical diffusive flux of turbulent kinetic energy proportional to the friction velocity cubed. The first boundary condition arises from consideration (simplification) of the turbulence closure scheme near boundaries, and the second arises from consideration of the influence of surface gravity waves on the transfer of turbulent kinetic energy from the wind to the water. Simulations using these two boundary conditions are compared to month-long observations of velocity, temperature, and salinity (as shallow as 2 m from the surface) from Knight Inlet, British Columbia, Canada. The circulation in the inlet is strongly influenced by the wind, tides, and freshwater runoff. The two boundary conditions produce simulations that are different down to a depth of at least 5 m. Somewhat more accurate simulations are produced by the second boundary condition. Also, simulations using the second boundary condition are more sensitive to variations in the roughness length. Based on the simulations, roughness lengths as large as 1 m (or greater) are possible.

Corresponding author address: Dr. Michael W. Stacey, Department of Physics, Royal Military College of Canada, Kingston, ON K7K 5LO, Canada.

Abstract

A numerical model that uses a level-2½ turbulence closure scheme is used to compare two boundary conditions for the turbulent energy at the air–sea interface. One boundary condition, the most commonly used, sets the turbulent kinetic energy proportional to the friction velocity squared, while the other sets the vertical diffusive flux of turbulent kinetic energy proportional to the friction velocity cubed. The first boundary condition arises from consideration (simplification) of the turbulence closure scheme near boundaries, and the second arises from consideration of the influence of surface gravity waves on the transfer of turbulent kinetic energy from the wind to the water. Simulations using these two boundary conditions are compared to month-long observations of velocity, temperature, and salinity (as shallow as 2 m from the surface) from Knight Inlet, British Columbia, Canada. The circulation in the inlet is strongly influenced by the wind, tides, and freshwater runoff. The two boundary conditions produce simulations that are different down to a depth of at least 5 m. Somewhat more accurate simulations are produced by the second boundary condition. Also, simulations using the second boundary condition are more sensitive to variations in the roughness length. Based on the simulations, roughness lengths as large as 1 m (or greater) are possible.

Corresponding author address: Dr. Michael W. Stacey, Department of Physics, Royal Military College of Canada, Kingston, ON K7K 5LO, Canada.

Save
  • Baker, P. D., 1992: Low frequency residual circulation in Knight Inlet, a fjord of coastal British Columbia. M.S. thesis, Dept. of Earth and Ocean Sciences, University of British Columbia, 184 pp.

  • ——, and S. Pond, 1995: The low-frequency residual circulation in Knight Inlet, British Columbia. J. Phys. Oceanogr.,25, 747–763.

  • Blumberg, A. F., and G. L. Mellor, 1987: A description of a three-dimensional coastal ocean circulation model. Three-Dimensional Coastal Ocean Models, N. S. Heaps, Ed., American Geophysical Union, 1–16.

  • Craig, P. D., and M. L. Banner, 1994: Modeling wave-enhanced turbulence in the ocean surface layer. J. Phys. Oceanogr.,24, 2546–2559.

  • Haney, R. L., 1991: On the pressure gradient force over steep topography in sigma coordinate ocean models. J. Phys. Oceanogr.,21, 610–619.

  • Lavelle, J. W., E. D. Cokelet, and G. A. Cannon, 1991: A model study of density intusions into and circulation within a deep, silled estuary: Puget Sound. J. Geophys. Res.,96, 16779–16800.

  • Mellor, G. L., and T. Yamada, 1982: Development of a turbulence closure model for geophysical fluid problems. Rev. Geophys. Space Phys.,20, 851–875.

  • Stacey, M. W., and S. Pond, 1992: A numerical model of the internal tide in Knight Inlet, British Columbia. Atmos.–Ocean,30, 383–418.

  • ——, ——, and Z. P. Nowak, 1995: A numerical model of the circulation in Knight Inlet, British Columbia, Canada. J. Phys. Oceanogr.,25, 1038–1062.